Premium
A preconditioning mass matrix to accelerate the convergence to the steady Euler solutions using explicit schemes
Author(s) -
Storti M. A.,
Baumann C. E.,
Idelsohn S. R.
Publication year - 1992
Publication title -
international journal for numerical methods in engineering
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.421
H-Index - 168
eISSN - 1097-0207
pISSN - 0029-5981
DOI - 10.1002/nme.1620340210
Subject(s) - mach number , rate of convergence , steady state (chemistry) , convergence (economics) , transonic , mathematics , matrix (chemical analysis) , euler's formula , mass matrix , euler equations , backward euler method , mathematical analysis , mechanics , physics , computer science , aerodynamics , materials science , computer network , channel (broadcasting) , chemistry , nuclear physics , neutrino , economics , composite material , economic growth
When explicit time marching algorithms are used to reach the steady state of problems governed by the Euler equations, the rate of convergence is strongly impaired both in the zones with low Mach number and in the zones with transonic flow, e.g. Mach ⩽ α and | Mach − 1| ⩽ α, with α ⩽ 0·2. The rate of convergence becomes slower as α diminishes. We show in this paper, with analytical and numerical results, how the use of a preconditioning mass matrix accelerates the convergence in the aforementioned ranges of Mach numbers. The preconditioning mass matrix (PMM) we advocate in this paper can be applied to any FEM/FVM that uses an explicit time marching scheme to find the steady state. The method's rate of convergence to the steady state is studied, and results for the one‐ and two‐dimiensional cases are presented. In Sections 1‐3, using the one‐dimensional Euler equations, we first explain why there exists a slow rate of convergence when the plain lumping of mass is used. Then the convergence rate to steady solutions is analysed from its two constituents, that is, convergence by absorption at the boundaries and by damping in the domain. Next we give the natural solution to this problem, and with several examples we show the effectiveness of the proposed mass matrix when compared with the plain scheme. In Sections 4‐8 we give the multidimensional version of the preconditioning mass matrix. We make a stability analysis and compare the group velocities and damping with and without the new mass matrix. To finish, we show the velocity of convergence for a common test problem.