An All-Speed Roe-type scheme and its asymptotic analysis of low Mach number behaviour

A new scheme, All-Speed-Roe scheme, was proposed for all speed flows. Compared with traditional preconditioned Roe scheme, All-Speed-Roe scheme changes non-linear eigenvalues in the numerical dissipation terms of Roe-type schemes. With an asymptotic analysis, the low Mach number behaviour of the scheme is studied theoretically in two ways. In one way, All-Speed-Roe scheme is regarded as finite magnification of Low-Speed-Roe scheme in the low Mach number limit. In the other way, a general form of All-Speed-Roe scheme is analyzed. Both ways demonstrate that All-Speed-Roe scheme has the same low Mach number behaviour as the original governing equation in the continuous case, which includes three important features: pressure variation scales with the square of the Mach number, the zero order velocity field is subject to a divergence constraint, and the second order pressure satisfies a Poisson-type equation in the case of constant-entropy. The analysis also leads to an unexpected conclusion that the velocity filed computed by traditional preconditioned Roe scheme does not satisfy the divergence constraint as the Mach number vanishes. Moreover, the analysis explains the reason of checkerboard decoupling and shows that momentum interpolation method provides a similar mechanism as traditional preconditioned Roe scheme inherently possesses to suppress checkerboard decoupling. In the end, general rulers for modifying non-linear eigenvalues are obtained. Finally, several numerical experiments are provided to support the theoretical analysis. All-Speed-Roe scheme has a sound foundation and is expected to be widely studied and applied to all speed flow calculations.

[1]  H. Guillard,et al.  On the behavior of upwind schemes in the low Mach number limit: II. Godunov type schemes , 2004 .

[2]  E. Turkel,et al.  PRECONDITIONING TECHNIQUES IN COMPUTATIONAL FLUID DYNAMICS , 1999 .

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

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

[5]  P. Sagaut,et al.  Large Eddy Simulation of Flow Around an Airfoil Near Stall , 2002 .

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

[7]  R. Klein Semi-implicit extension of a Godunov-type scheme based on low Mach number asymptotics , 1995 .

[8]  A. Meister Asymptotic based preconditioning technique for low Mach number flows , 2003 .

[9]  XueSong Li,et al.  Preconditioning method and engineering application of large eddy simulation , 2008 .

[10]  Pierre Sagaut,et al.  An algorithm for unsteady viscous flows at all speeds , 2000 .

[11]  Sang-Hyeon Lee,et al.  Cancellation problem of preconditioning method at low Mach numbers , 2007, J. Comput. Phys..

[12]  David L. Darmofal,et al.  The Importance of Eigenvectors for Local Preconditioners of the Euler Equations , 1996 .

[13]  J. Edwards,et al.  Low-Diffusion Flux-Splitting Methods for Flows at All Speeds , 1997 .

[14]  Philipp Birken,et al.  Stability of Preconditioned Finite Volume Schemes at Low Mach Numbers , 2005 .

[15]  Chun-wei Gu,et al.  Research of a Preconditioned Roe Scheme for All Speed Flows in Turbomachinery , 2007 .

[16]  A. Majda,et al.  Compressible and incompressible fluids , 1982 .