Convergence acceleration for computing steady-state compressible flow at low Mach numbers

A novel technique is introduced to accelerate the convergence of compressible low Mach number flow computations to the steady state. The stiffness due to the large disparity of flow velocity and acoustic wave speeds is bypassed by artificially reducing the speed of sound and thereby increasing the Mach number. This Mach number transformation is achieved by subtracting a constant value from the pressure in the entire flow field. Only the inviscid terms of the energy equation are influenced by that pressure decrease. The steady-state error is corrected by solving a scalar equation after each time step such that the steady-state solutions of the modified and non-modified schemes coincide. Thus, for each low Mach number simulation, one can obtain a convergence performance comparable to the corresponding simulation with a Mach number of about 0.4. This convergence acceleration is demonstrated for premixed laminar flames. If the present technique is implemented without time splitting, it corresponds to a novel low Mach number preconditioning.

[1]  Bertil Gustafsson,et al.  Navier-Stokes equations for almost incompressible flow , 1991 .

[2]  Philip L. Roe,et al.  Characteristic time-stepping or local preconditioning of the Euler equations , 1991 .

[3]  Bertil Gustafsson,et al.  A numerical method for incompressible and compressible flow problems with smooth solutions , 1986 .

[4]  K. Kuo.  Principles of combustion , 1986 .

[5]  M. Hafez,et al.  A unified approach for numerical simulation of viscous compressible and incompressible flows over adiabatic and isothermal walls , 1992 .

[6]  S. Venkateswaran,et al.  DEVELOPMENT OF LINEARIZED PRECONDITIONING METHODS FOR ENHANCING ROBUSTNESS AND EFFICIENCY OF EULER AND NAVIER-STOKES COMPUTATIONS , 1997 .

[7]  L. Eriksson.  A preconditioned Navier-Stokes solver for low Mach number flows , 1996 .

[8]  Jörn Sesterhenn,et al.  Flux-vector splitting for compressible low Mach number flow , 1993 .

[9]  Patrick Jenny,et al.  Rankine-Hugoniot-Riemann Solver Considering Source Terms and Multidimensional Effects , 1998 .

[10]  P. J O'Rourke,et al.  Two scaling transformations for the numerical computation of multidimensional unsteady laminar flames , 1979 .

[11]  Jörn Sesterhenn,et al.  On the Cancellation Problem in Calculating Compressible Low Mach Number Flows , 1999 .

[12]  V. Mousseau,et al.  Damped artificial compressibility method for steady-state low-speed flow calculations , 1991 .

[13]  Robert J. Kee,et al.  PREMIX :A F ORTRAN Program for Modeling Steady Laminar One-Dimensional Premixed Flames , 1998 .

[14]  C. L. Merkle,et al.  The application of preconditioning in viscous flows , 1993 .

[15]  R. J. Kee,et al.  Chemkin-II : A Fortran Chemical Kinetics Package for the Analysis of Gas Phase Chemical Kinetics , 1991 .

[16]  D. Gottlieb,et al.  Splitting methods for low Mach number Euler and Navier-Stokes equations , 1989 .

[17]  F. Wubs.  Notes on numerical fluid mechanics , 1985 .

[18]  Jürgen Warnatz,et al.  Numerical Methods in Laminar Flame Propagation , 1982 .

[19]  P. Jenny.  On the numerical solution of the compressible navier-stokes equations for reacting and non-reacting gas mixtures , 1997 .

[20]  G. Volpe.  Performance of compressible flow codes at low Mach numbers , 1993 .

[21]  Time-marching computations in the complete Mach/Reynolds number domain , 1993 .

[22]  Charles L. Merkle,et al.  Computation of low-speed compressible flows with time-marching procedures , 1988 .

[23]  Eli Turkel,et al.  PRECONDITIONING-SQUARED METHODS FOR MULTIDIMENSIONAL AERODYNAMICS , 1997 .

[24]  Bram van Leer,et al.  A local Navier-Stokes preconditioner for all Mach and cell Reynolds numbers , 1997 .