An implicit relaxation scheme is derived for the simulation of multidimensional flows at all Mach numbers, ranging from very small to order unity. An analytical proof of the asymptotic-preserving property is proposed and the divergence-free condition on the velocity in the incompressible regime is respected. The scheme possesses a general structure, which is independent of the considered state law and thus can be adopted to solve gas and fluid flows, but also deformations of elastic solids. This is achieved by adopting the Jin-Xin relaxation technique in order to get a linear transport operator. The spatial derivatives are thus independent of the equation of state and an easy implementation of fully implicit time discretizations is possible. Several validations on multidimensional tests are presented, showing that the correct numerical viscosity is recovered in both the fully compressible and the low Mach regimes. An algorithm to perform grid adaptivity is also proposed, via the computation of the entropy residual of the scheme.
An asymptotic-preserving all-speed scheme for fluid dynamics and nonlinear elasticity / Abbate, E.; Iollo, A.; Puppo, G.. - In: SIAM JOURNAL ON SCIENTIFIC COMPUTING. - ISSN 1064-8275. - 41:5(2019), pp. A2850-A2879. [10.1137/18M1232954]
An asymptotic-preserving all-speed scheme for fluid dynamics and nonlinear elasticity
Puppo G.
2019
Abstract
An implicit relaxation scheme is derived for the simulation of multidimensional flows at all Mach numbers, ranging from very small to order unity. An analytical proof of the asymptotic-preserving property is proposed and the divergence-free condition on the velocity in the incompressible regime is respected. The scheme possesses a general structure, which is independent of the considered state law and thus can be adopted to solve gas and fluid flows, but also deformations of elastic solids. This is achieved by adopting the Jin-Xin relaxation technique in order to get a linear transport operator. The spatial derivatives are thus independent of the equation of state and an easy implementation of fully implicit time discretizations is possible. Several validations on multidimensional tests are presented, showing that the correct numerical viscosity is recovered in both the fully compressible and the low Mach regimes. An algorithm to perform grid adaptivity is also proposed, via the computation of the entropy residual of the scheme.File | Dimensione | Formato | |
---|---|---|---|
Abbate_An-asymptotic-preserving_2019.pdf
accesso aperto
Tipologia:
Versione editoriale (versione pubblicata con il layout dell'editore)
Licenza:
Tutti i diritti riservati (All rights reserved)
Dimensione
2.86 MB
Formato
Adobe PDF
|
2.86 MB | Adobe PDF |
I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.