Volume 10, Issue 1
All Speed Scheme for the Low Mach Number Limit of the Isentropic Euler Equations

Pierre Degond & Min Tang

Commun. Comput. Phys., 10 (2011), pp. 1-31.

Published online: 2011-10

Preview Full PDF 519 1919
Export citation
  • Abstract

An all speed scheme for the Isentropic Euler equations is presented in this paper. When the Mach number tends to zero, the compressible Euler equations converge to their incompressible counterpart, in which the density becomes a constant. Increasing approximation errors and severe stability constraints are the main difficulty in the low Mach regime. The key idea of our all speed scheme is the special semi-implicit time discretization, in which the low Mach number stiff term is divided into two parts, one being treated explicitly and the other one implicitly. Moreover, the flux of the density equation is also treated implicitly and an elliptic type equation is derived to obtain the density. In this way, the correct limit can be captured without requesting the mesh size and time step to be smaller than the Mach number. Compared with previous semi-implicit methods [11,13,29], firstly, nonphysical oscillations can be suppressed by choosing proper parameter, besides, only a linear elliptic equation needs to be solved implicitly which reduces much computational cost. We develop this semi-implicit time discretization in the framework of a first order Local Lax-Friedrichs (or Rusanov) scheme and numerical tests are displayed to demonstrate its performances.

  • Keywords

  • AMS Subject Headings

  • Copyright

COPYRIGHT: © Global Science Press

  • Email address
  • BibTex
  • RIS
  • TXT
@Article{CiCP-10-1, author = {}, title = {All Speed Scheme for the Low Mach Number Limit of the Isentropic Euler Equations}, journal = {Communications in Computational Physics}, year = {2011}, volume = {10}, number = {1}, pages = {1--31}, abstract = {

An all speed scheme for the Isentropic Euler equations is presented in this paper. When the Mach number tends to zero, the compressible Euler equations converge to their incompressible counterpart, in which the density becomes a constant. Increasing approximation errors and severe stability constraints are the main difficulty in the low Mach regime. The key idea of our all speed scheme is the special semi-implicit time discretization, in which the low Mach number stiff term is divided into two parts, one being treated explicitly and the other one implicitly. Moreover, the flux of the density equation is also treated implicitly and an elliptic type equation is derived to obtain the density. In this way, the correct limit can be captured without requesting the mesh size and time step to be smaller than the Mach number. Compared with previous semi-implicit methods [11,13,29], firstly, nonphysical oscillations can be suppressed by choosing proper parameter, besides, only a linear elliptic equation needs to be solved implicitly which reduces much computational cost. We develop this semi-implicit time discretization in the framework of a first order Local Lax-Friedrichs (or Rusanov) scheme and numerical tests are displayed to demonstrate its performances.

}, issn = {1991-7120}, doi = {https://doi.org/10.4208/cicp.210709.210610a}, url = {http://global-sci.org/intro/article_detail/cicp/7432.html} }
TY - JOUR T1 - All Speed Scheme for the Low Mach Number Limit of the Isentropic Euler Equations JO - Communications in Computational Physics VL - 1 SP - 1 EP - 31 PY - 2011 DA - 2011/10 SN - 10 DO - http://doi.org/10.4208/cicp.210709.210610a UR - https://global-sci.org/intro/article_detail/cicp/7432.html KW - AB -

An all speed scheme for the Isentropic Euler equations is presented in this paper. When the Mach number tends to zero, the compressible Euler equations converge to their incompressible counterpart, in which the density becomes a constant. Increasing approximation errors and severe stability constraints are the main difficulty in the low Mach regime. The key idea of our all speed scheme is the special semi-implicit time discretization, in which the low Mach number stiff term is divided into two parts, one being treated explicitly and the other one implicitly. Moreover, the flux of the density equation is also treated implicitly and an elliptic type equation is derived to obtain the density. In this way, the correct limit can be captured without requesting the mesh size and time step to be smaller than the Mach number. Compared with previous semi-implicit methods [11,13,29], firstly, nonphysical oscillations can be suppressed by choosing proper parameter, besides, only a linear elliptic equation needs to be solved implicitly which reduces much computational cost. We develop this semi-implicit time discretization in the framework of a first order Local Lax-Friedrichs (or Rusanov) scheme and numerical tests are displayed to demonstrate its performances.

Pierre Degond & Min Tang. (2020). All Speed Scheme for the Low Mach Number Limit of the Isentropic Euler Equations. Communications in Computational Physics. 10 (1). 1-31. doi:10.4208/cicp.210709.210610a
Copy to clipboard
The citation has been copied to your clipboard