A method for incorporating Gauss' lasw into electromagnetic pic codes
Journal of Computational Physics
Numerical computation of internal & external flows: fundamentals of numerical discretization
Numerical computation of internal & external flows: fundamentals of numerical discretization
Notes on the eigensystem of magnetohydrodynamics
SIAM Journal on Applied Mathematics
Maintaining pressure positivity in magnetohydrodynamic simulations
Journal of Computational Physics
A solution-adaptive upwind scheme for ideal magnetohydrodynamics
Journal of Computational Physics
A positive conservative method for magnetohydrodynamics based on HLL and Roe methods
Journal of Computational Physics
The &Dgr; • = 0 constraint in shock-capturing magnetohydrodynamics codes
Journal of Computational Physics
A note on magnetic monopoles and the one-dimensional MHD Riemann problem
Journal of Computational Physics
AN APPROXIMATE RIEMANN SOLVER FOR MAGNETOHYDRODYNAMICS (That Works in More than One Dimension)
AN APPROXIMATE RIEMANN SOLVER FOR MAGNETOHYDRODYNAMICS (That Works in More than One Dimension)
Solution-Adaptive Magnetohydrodynamics for Space Plasmas: Sun-to-Earth Simulations
Computing in Science and Engineering
Adaptive numerical algorithms in space weather modeling
Journal of Computational Physics
Classical and semirelativistic magnetohydrodynamics with anisotropic ion pressure
Journal of Computational Physics
Hi-index | 31.46 |
We derive a system of equations for semirelativistic magnetohydrodynamics (MHD) in which the bulk speed and the sound speed of the plasma are nonrelativistic, but the Alfvén speed can be relativistic. The characteristic wave speeds of the modified equation set are determined and compared to the wave speeds in "classical" (MHD). The stability conditions of the semirelativistic MHD equations are also investigated in detail.This form of the MHD equations has a use beyond modeling flows with high Alfvén speeds. Even in cases with moderate Alfvén speeds, the semirelativistic form or certain approximations of it can be used to achieve accelerated numerical convergence to steady-state solutions by artificially reducing the speed of light, provided that the steady-state solutions of these equations are fully independent of the speed of light. Numerical tests are presented that demonstrate the behavior of solutions at high Alfvén speeds and the convergence acceleration that can be achieved when a steady-state solution is desired.