### A Preconditioning Mass Matrix to Accelerate the Convergence to the Steady Euler Solutions Using Explicit Schemes

#### Abstract

When explicit time muching algorithms are used to reach the steady state of problems governed by the Euler eqns, the rate of convergence is strongly impaired both in the zones with low Mach number and in the zones with transonic flow, let say Mach ≤α and |Mach-1|≤α, with α≤2. The rate oi convergence becomes slower as α diminishes.

We show in this paper, with analytical and numerical results, how the use of a preconditioning ma.ss 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-dimensional cases are presented.

In section 1, using the one-dimensional Euler eqns, we first explain why there exists & slaw rate of convergence when the plain lumping of mass is used. Then the convergence rate to steady solutions is analyzed 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 shaw the effectiveness oi the propolled mass matrix when compared with the plain scheme.

In section 2 we give the multidimensional version of the preconditioning mass matrix. We make a stability analysis and compare the group velocities and damping

with aDd without the new mass matrix. To finish, we show the velocity of convergence for a common test problem.

We show in this paper, with analytical and numerical results, how the use of a preconditioning ma.ss 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-dimensional cases are presented.

In section 1, using the one-dimensional Euler eqns, we first explain why there exists & slaw rate of convergence when the plain lumping of mass is used. Then the convergence rate to steady solutions is analyzed 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 shaw the effectiveness oi the propolled mass matrix when compared with the plain scheme.

In section 2 we give the multidimensional version of the preconditioning mass matrix. We make a stability analysis and compare the group velocities and damping

with aDd without the new mass matrix. To finish, we show the velocity of convergence for a common test problem.

#### Full Text:

PDF

Güemes 3450

S3000GLN Santa Fe, Argentina

Phone: 54-342-4511594 / 4511595 Int. 1006

Fax: 54-342-4511169

E-mail: amca(at)santafe-conicet.gov.ar

**Asociación Argentina de Mecánica Computacional**Güemes 3450

S3000GLN Santa Fe, Argentina

Phone: 54-342-4511594 / 4511595 Int. 1006

Fax: 54-342-4511169

E-mail: amca(at)santafe-conicet.gov.ar

**ISSN 2591-3522**