Conservative finite-volume framework and pressure-based algorithm for flows of incompressible, ideal-gas and real-gas fluids at all speeds

A conservative finite-volume framework, based on a collocated variable arrangement, for the simulation of flows at all speeds, applicable to incompressible, ideal-gas and real-gas fluids is proposed in conjunction with a fully-coupled pressure-based algorithm. The applied conservative discretisation and implementation of the governing conservation laws as well as the definition of the fluxes using a momentum-weighted interpolation are identical for incompressible and compressible fluids, and are suitable for complex geometries represented by unstructured meshes. Incompressible fluids are described by predefined constant fluid properties, while the properties of compressible fluids are described by the Noble-Abel-stiffened-gas model, with the definitions of density and specific static enthalpy of both incompressible and compressible fluids combined in a unified thermodynamic closure model. The discretised governing conservation laws are solved in a single linear system of equations for pressure, velocity and temperature. Together, the conservative finite-volume discretisation, the unified thermodynamic closure model and the pressure-based algorithm yield a conceptually simple, but versatile, numerical framework. The proposed numerical framework is validated thoroughly using a broad variety of test-cases, with Mach numbers ranging from 0 to 239, including viscous flows of incompressible fluids as well as the propagation of acoustic waves and transiently evolving supersonic flows with shock waves in ideal-gas and real-gas fluids. These results demonstrate the accuracy, robustness and the convergence, as well as the conservation of mass and energy, of the numerical framework for flows of incompressible and compressible fluids at all speeds, on structured and unstructured meshes.

[1]  Eli Turkel,et al.  Preconditioning and the Limit to the Incompressible Flow Equations , 1993 .

[2]  A. A. Amsden,et al.  A numerical fluid dynamics calculation method for all flow speeds , 1971 .

[3]  M. Darwish,et al.  A high-resolution pressure-based algorithm for fluid flow at all speeds , 2001 .

[4]  Feng Xiao,et al.  Unified formulation for compressible and incompressible flows by using multi-integrated moments II: Multi-dimensional version for compressible and incompressible flows , 2006, J. Comput. Phys..

[5]  S. I. Rubinow,et al.  The transverse force on a spinning sphere moving in a viscous fluid , 1961, Journal of Fluid Mechanics.

[6]  P. Woodward,et al.  The numerical simulation of two-dimensional fluid flow with strong shocks , 1984 .

[7]  Sukumar Chakravarthy,et al.  High Resolution Schemes and the Entropy Condition , 1984 .

[8]  T. Hou,et al.  Why nonconservative schemes converge to wrong solutions: error analysis , 1994 .

[9]  Sophie Papst,et al.  Computational Methods For Fluid Dynamics , 2016 .

[10]  I. Demirdzic,et al.  A finite volume method for computation of fluid flow in complex geometries , 1982 .

[11]  E. Dick,et al.  Mach-uniformity through the coupled pressure and temperature correction algorithm , 2005 .

[12]  Linear stability, transient energy growth, and the role of viscosity stratification in compressible plane Couette flow. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

[13]  Hrvoje Jasak,et al.  Error analysis and estimation for the finite volume method with applications to fluid flows , 1996 .

[14]  C. Munz,et al.  Multiple pressure variables methods for fluid flow at all Mach numbers , 2005 .

[15]  M. Darwish,et al.  A Fully Coupled Navier-Stokes Solver for Fluid Flow at All Speeds , 2014 .

[16]  R. W. MacCormack,et al.  A Numerical Method for Solving the Equations of Compressible Viscous Flow , 1981 .

[17]  Eli Turkel,et al.  Assessment of Preconditioning Methods for Multidimensional Aerodynamics , 1997 .

[18]  J. Ferziger.  Interfacial transfer in Tryggvason's method , 2003 .

[19]  P. Wesseling,et al.  A conservative pressure-correction method for flow at all speeds , 2003 .

[20]  F. Moukalled,et al.  A coupled finite volume solver for the solution of incompressible flows on unstructured grids , 2009, J. Comput. Phys..

[21]  J. Anderson,et al.  Modern Compressible Flow: With Historical Perspective , 1982 .

[22]  S. Patankar,et al.  Pressure based calculation procedure for viscous flows at all speeds in arbitrary configurations , 1988 .

[23]  William Gropp,et al.  Efficient Management of Parallelism in Object-Oriented Numerical Software Libraries , 1997, SciTools.

[24]  Yeng-Yung Tsui,et al.  A Pressure-Based Unstructured-Grid Algorithm Using High-Resolution Schemes for All-Speed Flows , 2008 .

[25]  A. Harten.  High Resolution Schemes for Hyperbolic Conservation Laws , 2017 .

[26]  S. A. Mohsen Karimian,et al.  Discretization and parallel performance of an unstructured finite volume Navier–Stokes solver , 2006 .

[27]  Erik Dick,et al.  Solving low Mach number Riemann problems by a momentum interpolation method , 2015, J. Comput. Phys..

[28]  M. Taeibi-Rahni,et al.  Three-dimensional compressible–incompressible turbulent flow simulation using a pressure-based algorithm , 2008 .

[29]  P. Colella,et al.  A second-order projection method for the incompressible navier-stokes equations , 1989 .

[30]  Stanley Osher,et al.  A second order primitive preconditioner for solving all speed multi-phase flows , 2005 .

[31]  A. Toutant.  General and exact pressure evolution equation , 2017, 2005.05148.

[32]  J. Marsden,et al.  A mathematical introduction to fluid mechanics , 1979 .

[33]  G. Soave.  Equilibrium constants from a modified Redlich-Kwong equation of state , 1972 .

[34]  Hester Bijl,et al.  A Unified Method for Computing Incompressible and Compressible Flows in Boundary-Fitted Coordinates , 1998 .

[35]  Fabian Denner,et al.  Fully-Coupled Balanced-Force VOF Framework for Arbitrary Meshes with Least-Squares Curvature Evaluation from Volume Fractions , 2014 .

[36]  A. Majda,et al.  Singular limits of quasilinear hyperbolic systems with large parameters and the incompressible limit , 1981 .

[37]  Peter Maurer,et al.  Introduction to Statistical Thermodynamics , 1960 .

[38]  P. Wesseling.  Principles of Computational Fluid Dynamics , 2000 .

[39]  Yeng-Yung Tsui,et al.  A Pressure-Correction Method for Incompressible Flows Using Unstructured Meshes , 2006 .

[40]  G. Sod.  A survey of several finite difference methods for systems of nonlinear hyperbolic conservation laws , 1978 .

[41]  D. Hoff.  The Zero-Mach Limit of Compressible Flows , 1998 .

[42]  Z. Lilek,et al.  A fourth-order finite volume method with colocated variable arrangement , 1995 .

[43]  Stéphane Popinet,et al.  An all-Mach method for the simulation of bubble dynamics problems in the presence of surface tension , 2018, J. Comput. Phys..

[44]  Richard Saurel,et al.  The Noble-Abel Stiffened-Gas equation of state , 2016 .

[45]  A. Harten,et al.  The artificial compression method for computation of shocks and contact discontinuities. I - Single conservation laws , 1977 .

[46]  Fabian Denner,et al.  Fully-coupled pressure-based finite-volume framework for the simulation of fluid flows at all speeds in complex geometries , 2017, J. Comput. Phys..

[47]  K. H. Chen,et al.  Primitive variable, strongly implicit calculation procedure for viscous flows at all speeds , 1991 .

[48]  Eli Turkel,et al.  Numerical Methods and Nature , 2006, J. Sci. Comput..

[49]  R. F. Warming,et al.  An Implicit Factored Scheme for the Compressible Navier-Stokes Equations , 1977 .

[50]  D. Peng,et al.  A New Two-Constant Equation of State , 1976 .

[51]  Pierre Degond,et al.  An Asymptotic-Preserving all-speed scheme for the Euler and Navier-Stokes equations , 2011, J. Comput. Phys..

[52]  J. Murthy,et al.  A PRESSURE-BASED METHOD FOR UNSTRUCTURED MESHES , 1997 .

[53]  B. Wachem,et al.  A coupled solver approach for multiphase flow calculations on collocated grids , 2006 .

[54]  E. Turkel,et al.  Preconditioned methods for solving the incompressible low speed compressible equations , 1987 .

[55]  M. Darwish,et al.  The Finite Volume Method in Computational Fluid Dynamics: An Advanced Introduction with OpenFOAM® and Matlab , 2015 .

[56]  Jayathi Y. Murthy,et al.  Pressure-Based Finite-Volume Methods in Computational Fluid Dynamics , 2007 .

[57]  N. N. Yanenko,et al.  Methods for the Localization of Singularities in Numerical Solutions of Gas Dynamics Problems , 1989 .

[58]  Fabian Denner.  Fully-coupled pressure-based algorithm for compressible flows: Linearisation and iterative solution strategies , 2018, Computers & Fluids.

[59]  Berend van Wachem,et al.  A coupled solver approach for multiphase flow problems , 2007 .

[60]  S. G. Rubin,et al.  A diagonally dominant second-order accurate implicit scheme , 1974 .

[61]  Ian A Johnston,et al.  The Noble-Abel Equation of State: Thermodynamic Derivations for Ballistics Modelling , 2005 .

[62]  Fabian Denner,et al.  Pressure-based algorithm for compressible interfacial flows with acoustically-conservative interface discretisation , 2018, J. Comput. Phys..

[63]  With Invariant Submanifolds,et al.  Systems of Conservation Laws , 2009 .

[64]  Fabian Denner,et al.  Unified formulation of the momentum-weighted interpolation for collocated variable arrangements , 2018, J. Comput. Phys..

[65]  R. Kunz,et al.  Development of an Implicit Method for Multi-fluid Flow Simulations , 1999 .

[66]  C. Rhie,et al.  Numerical Study of the Turbulent Flow Past an Airfoil with Trailing Edge Separation , 1983 .

[67]  Tarik Kousksou,et al.  Pressure-velocity coupling allowing acoustic calculation in low Mach number flow , 2012, J. Comput. Phys..

[68]  A. Chorin.  Numerical solution of the Navier-Stokes equations , 1968 .

[69]  Katja Bachmeier,et al.  Numerical Heat Transfer And Fluid Flow , 2016 .

[70]  S. M. H. Karimian,et al.  Pressure-based control-volume finite element method for flow at all speeds , 1995 .

[71]  O. Metayer,et al.  Élaboration des lois d'état d'un liquide et de sa vapeur pour les modèles d'écoulements diphasiques Elaborating equations of state of a liquid and its vapor for two-phase flow models , 2004 .

[72]  A. D. Gosman,et al.  The computation of compressible and incompressible recirculating flows by a non-iterative implicit scheme , 1986 .

[73]  William J. Rider,et al.  On sub-linear convergence for linearly degenerate waves in capturing schemes , 2008, J. Comput. Phys..

[74]  Eli Turkel,et al.  Review of preconditioning methods for fluid dynamics , 1993 .

[75]  Erik Dick,et al.  A combined momentum-interpolation and advection upstream splitting pressure-correction algorithm for simulation of convective and acoustic transport at all levels of Mach number , 2019, J. Comput. Phys..

[76]  The Lagrange multiplier and the stationary Stokes equations , 2017, 1703.08216.

[77]  A. A. Amsden,et al.  Numerical calculation of almost incompressible flow , 1968 .

[78]  Z. J. Chen,et al.  A coupled pressure-based computational method for incompressible/compressible flows , 2010, J. Comput. Phys..

[79]  Thomas J. R. Hughes,et al.  A comparative study of different sets of variables for solving compressible and incompressible flows , 1998 .

[80]  H. Swann.  The convergence with vanishing viscosity of nonstationary Navier-Stokes flow to ideal flow in ₃ , 1971 .

[81]  Feng Xiao,et al.  Unified formulation for compressible and incompressible flows by using multi-integrated moments I: one-dimensional inviscid compressible flow , 2004 .

[82]  C. Munz,et al.  The extension of incompressible flow solvers to the weakly compressible regime , 2003 .

[83]  M. Darwish,et al.  The Finite Volume Method , 2016 .

[84]  D. Spalding,et al.  A calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows , 1972 .

[85]  Fabian Denner,et al.  TVD differencing on three-dimensional unstructured meshes with monotonicity-preserving correction of mesh skewness , 2015, J. Comput. Phys..

[86]  J. Pain.  Fluid Dynamics , 1967, Nature.

[87]  S. M. H. Karimian,et al.  Pressure-Based Computational Method for Compressible and Incompressible Flows , 1994 .

[88]  A. Gosman,et al.  Solution of the implicitly discretised reacting flow equations by operator-splitting , 1986 .

[89]  I. Demirdzic,et al.  Numerical method for coupled fluid flow, heat transfer and stress analysis using unstructured moving meshes with cells of arbitrary topology , 1995 .

[90]  Tim Colonius,et al.  Finite-volume WENO scheme for viscous compressible multicomponent flows , 2014, J. Comput. Phys..

[91]  G. D. Raithby,et al.  The segregated approach to predicting viscous compressible fluid flows , 1986 .

[92]  C. Rhie.  Pressure-based Navier-Stokes solver using the multigrid method , 1986 .

[93]  M. Banholzer,et al.  A hybrid pressure‐based solver for nonideal single‐phase fluid flows at all speeds , 2018 .

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

[95]  Richard Saurel,et al.  Shock jump relations for multiphase mixtures with stiff mechanical relaxation , 2007 .

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