A weakly-compressible Cartesian grid approach for hydrodynamic flows

Abstract The present article aims at proposing an original strategy to solve hydrodynamic flows. In introduction, the motivations for this strategy are developed. It aims at modeling viscous and turbulent flows including complex moving geometries, while avoiding meshing constraints. The proposed approach relies on a weakly-compressible formulation of the Navier–Stokes equations. Unlike most hydrodynamic CFD (Computational Fluid Dynamics) solvers usually based on implicit incompressible formulations, a fully-explicit temporal scheme is used. A purely Cartesian grid is adopted for numerical accuracy and algorithmic simplicity purposes. This characteristic allows an easy use of Adaptive Mesh Refinement (AMR) methods embedded within a massively parallel framework. Geometries are automatically immersed within the Cartesian grid with an AMR compatible treatment. The method proposed uses an Immersed Boundary Method (IBM) adapted to the weakly-compressible formalism and imposed smoothly through a regularization function, which stands as another originality of this work. All these features have been implemented within an in-house solver based on this WCCH (Weakly-Compressible Cartesian Hydrodynamic) method which meets the above requirements whilst allowing the use of high-order ( > 3 ) spatial schemes rarely used in existing hydrodynamic solvers. The details of this WCCH method are presented and validated in this article.

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

[2]  S. Armfield,et al.  A representation of curved boundaries for the solution of the Navier-Stokes equations on a staggered three-dimensional Cartesian grid , 2003 .

[3]  A. S. Grove,et al.  An experimental investigation of the steady separated flow past a circular cylinder , 1964, Journal of Fluid Mechanics.

[4]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme V. A second-order sequel to Godunov's method , 1979 .

[5]  Salvatore Marrone,et al.  Numerical diffusive terms in weakly-compressible SPH schemes , 2012, Comput. Phys. Commun..

[6]  Madeleine Coutanceau,et al.  The early stage of development of the wake behind an impulsively started cylinder for 40 < Re < 104 , 1980, Journal of Fluid Mechanics.

[7]  Stephen M. Ruffin,et al.  An Embedded Boundary Cartesian Grid Scheme for Viscous Flows using a New Viscous Wall Boundary Condition Treatment , 2004 .

[8]  David Le Touzé,et al.  On distributed memory MPI-based parallelization of SPH codes in massive HPC context , 2016, Comput. Phys. Commun..

[9]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .

[10]  J. Morris,et al.  Modeling Low Reynolds Number Incompressible Flows Using SPH , 1997 .

[11]  Peter MacNeice,et al.  Paramesh: A Parallel Adaptive Mesh Refinement Community Toolkit , 2013 .

[12]  Eric Lamballais,et al.  A reference solution of the flow over a circular cylinder at Re = 40 , 2013, 1310.6641.

[13]  M. Uhlmann An immersed boundary method with direct forcing for the simulation of particulate flows , 2005, 1809.08170.

[14]  L. Sirovich,et al.  Modeling a no-slip flow boundary with an external force field , 1993 .

[15]  E. Toro Riemann Solvers and Numerical Methods for Fluid Dynamics , 1997 .

[16]  Nikolaus A. Adams,et al.  A conservative interface method for compressible flows , 2006, J. Comput. Phys..

[17]  P. Colella,et al.  Local adaptive mesh refinement for shock hydrodynamics , 1989 .

[18]  A. Colagrossi,et al.  Prediction of energy losses in water impacts using incompressible and weakly compressible models , 2015 .

[19]  P. Roe CHARACTERISTIC-BASED SCHEMES FOR THE EULER EQUATIONS , 1986 .

[20]  R. Verzicco,et al.  Large Eddy Simulation in Complex Geometric Configurations Using Boundary Body Forces , 2000 .

[21]  E. Balaras Modeling complex boundaries using an external force field on fixed Cartesian grids in large-eddy simulations , 2004 .

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

[23]  W. Shyy,et al.  Regular Article: An Accurate Cartesian Grid Method for Viscous Incompressible Flows with Complex Immersed Boundaries , 1999 .

[24]  D. Tritton Experiments on the flow past a circular cylinder at low Reynolds numbers , 1959, Journal of Fluid Mechanics.

[25]  B. V. Leer,et al.  Towards the ultimate conservative difference scheme. II. Monotonicity and conservation combined in a second-order scheme , 1974 .

[26]  S. Dennis,et al.  Numerical solutions for steady flow past a circular cylinder at Reynolds numbers up to 100 , 1970, Journal of Fluid Mechanics.

[27]  R. Verzicco,et al.  Combined Immersed-Boundary Finite-Difference Methods for Three-Dimensional Complex Flow Simulations , 2000 .

[28]  H. S. Udaykumar,et al.  A Sharp Interface Cartesian Grid Methodfor Simulating Flows with ComplexMoving Boundaries , 2001 .

[29]  Salvatore Marrone,et al.  SPH accuracy improvement through the combination of a quasi-Lagrangian shifting transport velocity and consistent ALE formalisms , 2016, J. Comput. Phys..

[30]  J. Monaghan Smoothed particle hydrodynamics , 2005 .

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

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

[33]  U. Ghia,et al.  High-Re solutions for incompressible flow using the Navier-Stokes equations and a multigrid method , 1982 .

[34]  Eleuterio F. Toro,et al.  Finite-volume WENO schemes for three-dimensional conservation laws , 2004 .

[35]  M. Berger,et al.  Progress Towards a Cartesian Cut-Cell Method for Viscous Compressible Flow , 2012 .

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

[37]  O. Botella,et al.  BENCHMARK SPECTRAL RESULTS ON THE LID-DRIVEN CAVITY FLOW , 1998 .