On physical consistency of HLL-type Riemann solvers for ideal MHD
On physical consistency of HLL-type Riemann solvers for ideal MHD
Approximate Riemann solvers are widely used for solving hyperbolic conservation laws, including those of magnetohydrodynamics (MHD). However, due to the nonlinearity and complexity of MHD, obtaining accurate and robust numerical solutions to MHD equations is non-trivial, and it may be challenging for an approximate MHD Riemann solver to preserve the positivity of scalar variables, particularly when the plasma beta is low. As we have identified that the inconsistency between the numerically calculated magnetic field and magnetic energy may be responsible for the loss of positivity of scalar variables, we propose a physical consistency condition and implement it in HLL-type MHD Riemann solvers, thus alleviating the erroneous magnetic field solutions breaking scalar positivity. In addition, (I) for the HLLC-type scheme, we have designed a revised two-state approximation, specifically reducing numerical error in magnetic field solutions, and (II) for the HLLD-type scheme, we replace the constant total pressure assumption by a three-state assumption for the intermediate thermal energy, which is more consistent with our other assumptions. The proposed schemes perform better in numerical examples with low plasma \b{eta}. Moreover, we explained the energy error introduced during time integration.
Fan Zhang、Andrea Lani、Stefaan Poedts
电工基础理论高电压技术
Fan Zhang,Andrea Lani,Stefaan Poedts.On physical consistency of HLL-type Riemann solvers for ideal MHD[EB/OL].(2025-07-14)[2025-07-25].https://arxiv.org/abs/2507.10420.点此复制
评论