Treatment of solid objects in the Pencil Code using an immersed boundary method and overset grids

ABSTRACT Two methods for solid body representation in flow simulations available in the Pencil Code are the immersed boundary method and overset grids. These methods are quite different in terms of computational cost, flexibility and numerical accuracy. We present here an investigation of the use of the different methods with the purpose of assessing their strengths and weaknesses. At present, the overset grid method in the Pencil Code can only be used for representing cylinders in the flow. For this task, it surpasses the immersed boundary method in yielding highly accurate solutions at moderate computational costs. This is partly due to local grid stretching and a body-conformal grid, and partly due to the possibility of working with local time step restrictions on different grids. The immersed boundary method makes up the lack of computational efficiency with flexibility in regard to application to complex geometries, due to a recent extension of the method that allows our implementation of it to represent arbitrarily shaped objects in the flow.

[1]  Tomasz G. Drozda,et al.  Computational Fluid Mechanics and Heat Transfer , 2013 .

[2]  J. Quirk.  An alternative to unstructured grids for computing gas dynamic flows around arbitrarily complex two-dimensional bodies , 1994 .

[3]  Carlos Pantano,et al.  A Cartesian-based embedded geometry technique with adaptive high-order finite differences for compressible flow around complex geometries , 2014, J. Comput. Phys..

[4]  Sanjay Mittal,et al.  Excitation of shear layer instability in flow past a cylinder at low Reynolds number , 2005 .

[5]  E. Balaras,et al.  A general reconstruction algorithm for simulating flows with complex 3D immersed boundaries on Cartesian grids , 2003 .

[6]  Dan S. Henningson,et al.  High Order Accurate Solution of Flow Past a Circular Cylinder , 2006, J. Sci. Comput..

[7]  X. Gloerfelt,et al.  Study of interpolation methods for high-accuracy computations on overlapping grids , 2012 .

[8]  Derek M. Causon,et al.  Developments in Cartesian cut cell methods , 2003, Math. Comput. Simul..

[9]  Robert L. Meakin.  The Chimera Method of Simulation for Unsteady Three-Dimensional Viscous Flow , 1996 .

[10]  Gianluca Iaccarino,et al.  IMMERSED BOUNDARY METHODS , 2005 .

[11]  J. Benek,et al.  A 3-D Chimera Grid Embedding Technique , 1985 .

[12]  Jungwoo Kim,et al.  An immersed-boundary finite-volume method for simulations of flow in complex geometries , 2001 .

[13]  Ralph Noack,et al.  SUGGAR: A General Capability for Moving Body Overset Grid Assembly , 2005 .

[14]  Derek M. Causon,et al.  Calculation of shallow water flows using a Cartesian cut cell approach , 2000 .

[15]  C. Peskin.  The immersed boundary method , 2002, Acta Numerica.

[16]  Rajat Mittal,et al.  A versatile sharp interface immersed boundary method for incompressible flows with complex boundaries , 2008, J. Comput. Phys..

[17]  R. Mittal,et al.  A Computational Study of the Effect of False Vocal Folds on Glottal Flow and Vocal Fold Vibration During Phonation , 2009, Annals of Biomedical Engineering.

[18]  Vanessa Hertzog,et al.  Computational Fluid Mechanics And Heat Transfer , 2016 .

[19]  W. Henshaw,et al.  Composite overlapping meshes for the solution of partial differential equations , 1990 .

[20]  Fu-jun Wang,et al.  Quantitative numerical analysis of flow past a circular cylinder at Reynolds number between 50 and 200 , 2013 .

[21]  Miguel R. Visbal,et al.  On the use of higher-order finite-difference schemes on curvilinear and deforming meshes , 2002 .

[22]  Jianren Fan,et al.  A ghost-cell immersed boundary method for simulations of heat transfer in compressible flows under different boundary conditions , 2016 .

[23]  N. E. L. Haugen,et al.  High-order overset grid method for detecting particle impaction on a cylinder in a cross flow , 2018, International Journal of Computational Fluid Dynamics.

[24]  J. Williamson.  Low-storage Runge-Kutta schemes , 1980 .

[25]  C. Peskin.  Flow patterns around heart valves: A numerical method , 1972 .

[26]  F. Durst,et al.  Heating effect on steady and unsteady horizontal laminar flow of air past a circular cylinder , 2004 .

[27]  J. L. Steger,et al.  On the use of composite grid schemes in computational aerodynamics , 1987 .

[28]  Cheng Liu,et al.  An efficient immersed boundary treatment for complex moving object , 2014, J. Comput. Phys..

[29]  Steven J. Owen,et al.  A Survey of Unstructured Mesh Generation Technology , 1998, IMR.

[30]  T. Poinsot.  Boundary conditions for direct simulations of compressible viscous flows , 1992 .

[31]  G. Lauder,et al.  Computational modelling and analysis of the hydrodynamics of a highly deformable fish pectoral fin , 2010, Journal of Fluid Mechanics.

[32]  J. Ferziger,et al.  A ghost-cell immersed boundary method for flow in complex geometry , 2002 .

[33]  Richard Shock,et al.  Prediction of vortex shedding from a circular cylinder using a volumetric Lattice-Boltzmann boundary approach , 2009 .

[34]  Thomas Rung,et al.  Analysis of non-conservative interpolation techniques in overset grid finite-volume methods , 2017 .

[35]  Haecheon Choi,et al.  Numerical solutions of flow past a circular cylinder at Reynolds numbers up to 160 , 1998 .

[36]  N. Haugen,et al.  Numerical investigation of free-stream turbulence effects on the transition-in-wake state of flow past a circular cylinder , 2018 .

[37]  Jianren Fan,et al.  A ghost-cell immersed boundary method for the simulations of heat transfer in compressible flows under different boundary conditions Part-II: Complex geometries , 2017 .

[38]  Nils Erland L. Haugena.  Hydrodynamic and hydromagnetic energy spectra from large eddy simulations , 2006 .

[39]  J. L. Steger,et al.  A chimera grid scheme , 2011 .

[40]  A. Kageyama,et al.  ``Yin-Yang grid'': An overset grid in spherical geometry , 2004, physics/0403123.

[41]  Jung Hee Seo,et al.  A high-order immersed boundary method for acoustic wave scattering and low-Mach number flow-induced sound in complex geometries , 2011, J. Comput. Phys..

[42]  O. Posdziech,et al.  A systematic approach to the numerical calculation of fundamental quantities of the two-dimensional flow over a circular cylinder , 2007 .

[43]  Jianren Fan,et al.  A ghost-cell based high-order immersed boundary method for inter-phase heat transfer simulation , 2014 .

[44]  Stuart E. Rogers,et al.  Pegasus 5: An Automated Pre-Processor for Overset-Grid Cfd , 2013 .

[45]  Matthias Meinke,et al.  A Cartesian cut-cell method for sharp moving boundaries , 2011 .

[46]  W. Dobler,et al.  Magnetic Field Generation in Fully Convective Rotating Spheres , 2006 .

[47]  Abdellah Hadjadj,et al.  On the use of immersed boundary methods for shock/obstacle interactions , 2011, J. Comput. Phys..

[48]  Danesh K. Tafti,et al.  A new approach for conjugate heat transfer problems using immersed boundary method for curvilinear grid based solvers , 2014, J. Comput. Phys..

[49]  Odd M. Faltinsen,et al.  A local directional ghost cell approach for incompressible viscous flow problems with irregular boundaries , 2008, J. Comput. Phys..

[50]  B. Strand.  Summation by parts for finite difference approximations for d/dx , 1994 .

[51]  Wolfgang Schröder,et al.  An accurate moving boundary formulation in cut-cell methods , 2013, J. Comput. Phys..

[52]  Miguel R. Visbal,et al.  High-Order-Accurate Methods for Complex Unsteady Subsonic Flows , 1999 .

[53]  R. Pletcher,et al.  Computational Fluid Mechanics and Heat Transfer. By D. A ANDERSON, J. C. TANNEHILL and R. H. PLETCHER. Hemisphere, 1984. 599 pp. $39.95. , 1986, Journal of Fluid Mechanics.

[54]  Lars Schiøtt Sørensen,et al.  An introduction to Computational Fluid Dynamics: The Finite Volume Method , 1999 .

[55]  An Immersed Boundary Method on Unstructured Cartesian Meshes for Incompressible Flows with Heat Transfer , 2006 .

[56]  D. Mavriplis.  UNSTRUCTURED GRID TECHNIQUES , 1997 .

[57]  Robert L. Street,et al.  A composite multigrid method for calculating unsteady incompressible flows in geometrically complex domains , 1995 .

[58]  J. Nordström,et al.  Summation by Parts Operators for Finite Difference Approximations of Second-Derivatives with Variable Coefficients , 2004, Journal of Scientific Computing.

[59]  Stuart E. Rogers,et al.  PEGASUS 5: An Automated Preprocessor for Overset-Grid Computational Fluid Dynamics , 2003 .

[60]  D. Gaitonde,et al.  Pade-Type Higher-Order Boundary Filters for the Navier-Stokes Equations , 2000 .

[61]  Steinar Kragset,et al.  Particle impaction on a cylinder in a crossflow as function of Stokes and Reynolds numbers , 2010, Journal of Fluid Mechanics.

[62]  Björn Sjögreen,et al.  Conservative and Non-Conservative Interpolation between Overlapping Grids for Finite Volume Solutions of Hyperbolic Problems , 1994 .

[63]  Scott E. Sherer,et al.  High-order compact finite-difference methods on general overset grids , 2005 .