TY - JOUR
T1 - An approximate riemann solver for ideal magnetohydrodynamics
AU - Dai, Wenlong
AU - Woodward, Paul R.
PY - 1994/4
Y1 - 1994/4
N2 - To construct numerical schemes of the Godunov type for solving magnetohydrodynamical (MHD) problems, an approximate method of solving the MHD Riemann problem is required in order to calculate the time-averaged fluxes at the interfaces of numerical zones. Such an MHD Riemann solver is presented here which treats all waves emanating from the initial discontinuity as themselves discontinuous. Thus shock jump conditions are used for rarefactions, which limits the applicability of this work to weak rarefactions, the case most important for computation. The solutions from our approximate MHD Riamnn solver consist of two fast waves (either shock or rarefaction) two rotational discontinuities, two rarefaction waves (either shock or rarefaction), and one contact discontinuity for a general MHD Riemann problem. In order to display rotational discontinuities, a three-component model is necessary. Only under very limited circumstances is there no rotational discontinuity involved and thus the two component approximation may be used in the MHD Riemann problem. The solutions of the MHD Riemann problem in the shock tube problem which generates the compound wave in the earlier work contain two fast rarefaction waves, two slow shocks, one contact discontinuity, and one rotational discontinuity in our formalism.
AB - To construct numerical schemes of the Godunov type for solving magnetohydrodynamical (MHD) problems, an approximate method of solving the MHD Riemann problem is required in order to calculate the time-averaged fluxes at the interfaces of numerical zones. Such an MHD Riemann solver is presented here which treats all waves emanating from the initial discontinuity as themselves discontinuous. Thus shock jump conditions are used for rarefactions, which limits the applicability of this work to weak rarefactions, the case most important for computation. The solutions from our approximate MHD Riamnn solver consist of two fast waves (either shock or rarefaction) two rotational discontinuities, two rarefaction waves (either shock or rarefaction), and one contact discontinuity for a general MHD Riemann problem. In order to display rotational discontinuities, a three-component model is necessary. Only under very limited circumstances is there no rotational discontinuity involved and thus the two component approximation may be used in the MHD Riemann problem. The solutions of the MHD Riemann problem in the shock tube problem which generates the compound wave in the earlier work contain two fast rarefaction waves, two slow shocks, one contact discontinuity, and one rotational discontinuity in our formalism.
UR - http://www.scopus.com/inward/record.url?scp=28144445934&partnerID=8YFLogxK
UR - http://www.scopus.com/inward/citedby.url?scp=28144445934&partnerID=8YFLogxK
U2 - 10.1006/jcph.1994.1069
DO - 10.1006/jcph.1994.1069
M3 - Article
AN - SCOPUS:28144445934
VL - 111
SP - 354
EP - 372
JO - Journal of Computational Physics
JF - Journal of Computational Physics
SN - 0021-9991
IS - 2
ER -