Solving 3D magnetohydrostatics with RBF-FD: Applications to the solar corona

We present a novel magnetohydrostatic numerical model that solves directly for the force-balanced magnetic field in the solar corona. This model is constructed with Radial Basis Function Finite Differences (RBF-FD), specifically 3D polyharmonic splines plus polynomials, as the core discretization. This set of PDEs is particularly difficult to solve since in the limit of the forcing going to zero it becomes ill-posed with a multitude of solutions. For the forcing equal to zero there are no numerically tractable solutions. For finite forcing, the ability to converge onto a physically viable solution is delicate as will be demonstrated. The static force-balance equations are of a hyperbolic nature, in that information of the magnetic field travels along characteristic surfaces, yet they require an elliptic type solver approach for a sparse overdetermined ill-conditioned system. As an example, we reconstruct a highly nonlinear analytic model designed to represent long-lived magnetic structures observed in the solar corona.