Semi-conservative reduced speed of sound technique for low Mach number flows with large density variations

Abstract

The reduced speed of sound technique (RSST) has been used for efficient simulation of low Mach number flows in solar and stellar convection zones. The basic RSST equations are hyperbolic, and are suitable for parallel computation by domain decomposition. The application of RSST is limited to cases where density perturbations are much smaller than the background density. In addition, non-conservative variables are required to be evolved using this method, which is not suitable in cases where discontinuities like shock waves co-exist in a single numerical domain. In this study, we suggest a new semi-conservative formulation of the RSST that can be applied to low Mach number flows with large density variations. We derive the wave speed of the original and newly suggested methods to clarify that these methods can reduce the speed of sound without affecting the entropy wave. The equations are implemented using the finite volume method. Several numerical tests are carried out to verify the suggested methods. The analysis and numerical results show that the original RSST is not applicable when mass density variations are large. In contrast, the newly suggested methods are found to be efficient in such cases. We also suggest variants of the RSST that conserve momentum in the machine precision. The newly suggested variants are formulated as semi-conservative equations, which reduce to the conservative form of the Euler equations when the speed of sound is not reduced. This property is advantageous when both high and low Mach number regions are included in the numerical domain. The newly suggested forms of RSST can be applied to a wider range of low Mach number flows.

arXiv comment12 pages, 10 figures, accepted for publication in Astronomy & Astrophysics

References in corpus (13)

• An Unsplit Godunov Method for Ideal MHD via Constrained Transport
• Radiative MHD simulation of sunspot structure
• Solar differential rotation and meridional flow: The role of a subadiabatic tachocline for the Taylor-Proudman balance
• Structure and Evolution of Giant Cells in Global Models of Solar Convection
• A Well-Posed Kelvin-Helmholtz Instability Test and Comparison
• Three-dimensional Numerical Simulations of Rayleigh-Taylor Unstable Flames in Type Ia Supernovae
• High-resolution calculation of the solar global convection with the reduced speed of sound technique: II. Near surface shear layer with the rotation
• MAESTRO: An Adaptive Low Mach Number Hydrodynamics Algorithm for Stellar Flows
• Three-dimensional Magnetohydrodynamic Simulation of the Formation of Solar Chromospheric Jets with Twisted Magnetic Field Lines
• High-resolution calculations of the solar global convection with the reduced speed of sound technique: I. The structure of the convection and the magnetic field without the rotation
• Numerical Calculation of Convection with Reduced Speed of Sound Technique
• Effect of coronal temperature on the scale of solar chromospheric jets
• Magnetic flux concentrations from turbulent stratified convection

Cited by in corpus (6)

• On rising magnetic flux tube and formation of sunspots in a deep domain
• Spontaneous Generation of delta-sunspots in Convective Magnetohydrodynamic Simulation of Magnetic Flux Emergence
• MAESTROeX: A Massively Parallel Low Mach Number Astrophysical Solver
• A Multistate Low-dissipation Advection Upstream Splitting Method for Ideal Magnetohydrodynamics
• Impact of subsurface convective flows on the formation of sunspot magnetic field and energy build-up
• Low-Frequency Internal Gravity Waves are Pseudo-incompressible