2016/09/30 by S. A. Gilchrist, D. C. Braun, G. Barnes · 1 citation
Earth and Planetary Sciences · Physics and Astronomy · #Cosmology and Gravitation Theories #Geophysics and Gravity Measurements #Solar and Space Plasma Dynamics #astro-ph.SR
paper · pdf · doi:10.1007/s11207-016-0992-0
Accepted for publication in Solar Physics
arxiv created 2016/10/27 · arxiv updated 2016/10/28
Magnetohydrostatic models of the solar atmosphere are often based on idealized analytic solutions because the underlying equations are too difficult to solve in full generality. Numerical approaches, too, are often limited in scope and have tended to focus on the two-dimensional problem. In this article we develop a numerical method for solving the nonlinear magnetohydrostatic equations in three dimensions. Our method is a fixed-point iteration scheme that extends the method of Grad and Rubin (Proc. 2nd Int. Conf. on Peaceful Uses of Atomic Energy 31, 190, 1958) to include a finite gravity force. We apply the method to a test case to demonstrate the method in general and our implementation in code in particular.