Comparative Study of Global MHD Simulations of the Terrestrial Magnetosphere With Different Numerical Schemes

• Published in: IEEE transactions on plasma science Vol. 38; no. 9; pp. 2229 - 2235
• Main Authors: Matsumoto, Y, Terada, N, Miyoshi, T, Fukazawa, K, Umeda, T, Ogino, T, Seki, K
• Format: Journal Article
• Language: English
• Published: New York IEEE 01-09-2010; The Institute of Electrical and Electronics Engineers, Inc. (IEEE)
• Subjects: Algorithms; Comparative analysis; Computational modeling; Computer simulation; Electric shock; Equations; Finite difference method; Finite difference methods; Finite volume methods; Interplanetary magnetic fields; Interpolation; Lagrange multiplier; Magnetic fields; Magnetohydrodynamics; Magnetosphere; Mathematical analysis; Mathematical models; MHD; Numerical analysis; numerical model; Numerical models; Plasma simulation; Simulation; Studies
• ISSN: 0093-3813; 1939-9375
• DOI: 10.1109/TPS.2010.2056704

We compare recent global MHD simulation models of the terrestrial magnetosphere based on different numerical schemes. The schemes include the finite-difference method based on the modified leapfrog (MLF) scheme, and the semi-Lagrangian scheme based on the constrained interpolation profile (CIP) algorithm. With the two models, we examined the simulation under a northward interplanetary magnetic field (IMF) condition. As a result, we found out that the two simulation models give consistent results on the magnetopause locations at the subsolar point and the terminator, and the overall structures of the cusp in the meridian plane. However, discrepancies are also found in the location and jump conditions of the bow shock. The MLF model showed higher thermal pressure value and weaker magnetic field strength in the downstream than those in the CIP model. The difference in the jump condition across the shock is also reflected in the difference in the length of the magnetotail in the two models. The magnetotail is shorter in the CIP model than that in the MLF model. We conclude that further comparative studies with finite-volume methods are necessary to verify the solution of the bow shock formation and the location of the last closed field line under northward IMF conditions.