Unified formulation for compressible and incompressible flows by using multi-integrated moments I: one-dimensional inviscid compressible flow

This paper presents a novel numerical formulation for computing flows of any Mach number. The following features make the present scheme different from the existing works: (1) The spatial discretizations are constructed by using two kinds of moments which are defined as the volume integrated average (VIA) and the surface integrated average (SIA), both of which are computed as the prognostic variables; (2) the CIP-CSL3 (constrained interpolation profile-conservative semi-lagrangian with third-order polynomial function) scheme is used for the advection transport; (3) a simple artificial compression can be devised by modifying the slope parameter in the CIP-CSL3 scheme. Moreover, the conservation is exactly assured for the VIA quantities. As the first of the series, this paper mainly focuses on the computations of inviscid compressible flows. Our numerical results show that the present method gives correct shock speed, well-resolved shock front, contact discontinuity and rarefaction waves for a wide spectrum of test problems.

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

[2]  P. Moin,et al.  Application of a Fractional-Step Method to Incompressible Navier-Stokes Equations , 1984 .

[3]  Huanan Yang,et al.  An artificial compression method for ENO schemes - The slope modification method. [essentially nonoscillatory , 1990 .

[4]  Takashi Yabe,et al.  A universal solver for hyperbolic equations by cubic-polynomial interpolation I. One-dimensional solver , 1991 .

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

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

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

[8]  F. Harlow,et al.  Numerical Calculation of Time‐Dependent Viscous Incompressible Flow of Fluid with Free Surface , 1965 .

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

[10]  Philip L. Roe,et al.  A Well-Behaved TVD Limiter for High-Resolution Calculations of Unsteady Flow , 1997 .

[11]  T. Yabe,et al.  The constrained interpolation profile method for multiphase analysis , 2001 .

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

[13]  Wei Shyy,et al.  Adaptive grid computation for inviscid compressible flows using a pressure correction method , 1988 .

[14]  Wayne A. Smith,et al.  Preconditioning Applied to Variable and Constant Density Flows , 1995 .

[15]  Mark L. Wilkins,et al.  Use of artificial viscosity in multidimensional fluid dynamic calculations , 1980 .

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

[17]  Takashi Yabe,et al.  A universal solver for hyperbolic equations by cubic-polynomial interpolation. II, Two- and three-dimensional solvers , 1991 .

[18]  T. Yabe,et al.  Completely conservative and oscillationless semi-Lagrangian schemes for advection transportation , 2001 .

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

[20]  Takashi Yabe,et al.  Shock capturing with improved numerical viscosity in primitive Euler representation , 1999 .

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

[22]  Chi-Wang Shu,et al.  Efficient Implementation of Weighted ENO Schemes , 1995 .

[23]  Barry Koren,et al.  Analysis of preconditioning and multigrid for Euler flows with low-subsonic regions , 1995, Adv. Comput. Math..

[24]  Feng Xiao,et al.  Profile-modifiable Conservative Transport Schemes and a Simple Multi-integrated Moment Formulation for Hydrodynamics , 2003 .

[25]  J. McGregor,et al.  Economical Determination of Departure Points for Semi-Lagrangian Models , 1993 .

[26]  Karkenahalli Srinivas,et al.  Computational Fluid Dynamics 2002 , 2003 .

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

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

[29]  H. Guillard,et al.  On the behaviour of upwind schemes in the low Mach number limit , 1999 .

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

[31]  T. Yabe,et al.  An Exactly Conservative Semi-Lagrangian Scheme (CIP–CSL) in One Dimension , 2001 .

[32]  Feng Xiao,et al.  An efficient method for capturing free boundaries in multi‐fluid simulations , 2003 .

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

[34]  G. D. van Albada,et al.  A comparative study of computational methods in cosmic gas dynamics , 1982 .

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

[36]  Takashi Yabe,et al.  Unified Numerical Procedure for Compressible and Incompressible Fluid , 1991 .

[37]  T. Yabe,et al.  Conservative and oscillation-less atmospheric transport schemes based on rational functions , 2002 .

[38]  P. Roe.  Approximate Riemann Solvers, Parameter Vectors, and Difference Schemes , 1997 .

[39]  C. W. Hirt.  Heuristic stability theory for finite-difference equations☆ , 1968 .

[40]  James J. McGuirk,et al.  Shock capturing using a pressure-correction method , 1989 .

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