2013/08/28 by S. A. Gilchrist, M. S. Wheatland · 1 citation
Biochemistry, Genetics and Molecular Biology · Engineering · Physics and Astronomy · #Advanced Numerical Methods in Computational Mathematics #Electromagnetic Simulation and Numerical Methods #Geomagnetism and Paleomagnetism Studies #Magnetic field #Nonlinear system #Polar #Polar coordinate system #Spherical coordinate system #Spherical geometry #Spherical harmonics #Spherical shell #Spin-weighted spherical harmonics #Vector spherical harmonics #astro-ph.SR
paper · pdf · doi:10.1007/s11207-013-0406-5
published as Sol Phys, 289, 1153 (2014) · Accepted for publication in Solar Physics
arxiv created 2013/08/28 · openalex publication_date 2013/09/20 · openalex created_date 2016/06/24 · arxiv updated 2016/10/28 · openalex updated_date 2026/08/05
We present a code for solving the nonlinear force-free equations in spherical polar geometry, with the motivation of modeling the magnetic field in the corona. The code is an implementation of the Grad-Rubin method. Our method is applicable to a spherical domain of arbitrary angular size. The implementation is based on a global spectral representation for the magnetic field which makes no explicit assumptions about the form of the magnetic field at the transverse boundaries of the domain. We apply the code to a bipolar test case with analytic boundary conditions, and we demonstrate the convergence of the Grad-Rubin method, and the self-consistency of the resulting numerical solution.