Gradient-based optimization of 3D MHD equilibria

Using recently developed adjoint methods for computing the shape derivatives of functions that depend on magnetohydrodynamic (MHD) equilibria (Antonsen et al., J. Plasma Phys., vol. 85, issue 2, 2019; Paul et al., J. Plasma Phys., vol. 86, issue 1, 2020), we present the first example of analytic gradient-based optimization of fixed-boundary stellarator equilibria. We take advantage of gradient information to optimize figures of merit of relevance for stellarator design, including the rotational transform, magnetic well and quasi-symmetry near the axis. With the application of the adjoint method, we reduce the number of equilibrium evaluations by the dimension of the optimization space ( ${\sim }50\text {--}500$ ) in comparison with a finite-difference gradient-based method. We discuss regularization objectives of relevance for fixed-boundary optimization, including a novel method that prevents self-intersection of the plasma boundary. We present several optimized equilibria, including a vacuum field with very low magnetic shear throughout the volume.

Global coronal and heliospheric magnetic field modelling for Solar Orbiter

<p>Computing the solar coronal magnetic field and plasma<br>environment is an important research topic on it's own right<br>and also important for space missions like Solar Orbiter to<br>guide the analysis of remote sensing and in-situ instruments.<br>In the inner solar corona plasma forces can be neglected and<br>the field is modelled under the assumption of a vanishing<br>Lorentz-force. Further outwards (above about two solar radii)<br>plasma forces and the solar wind flow has to be considered.<br>Finally in the heliosphere one has to consider that the Sun<br>is rotating and the well known Parker-spiral forms.<br>We have developed codes based on optimization principles<br>to solve nonlinear force-free, magneto-hydro-static and<br>stationary MHD-equilibria. In the present work we want to<br>extend these methods by taking the solar rotation into account.</p>

An Optimization Principle for Computing Stationary MHD Equilibria with Solar Wind Flow

In this work we describe a numerical optimization method for computing stationary MHD equilibria. The newly developed code is based on a nonlinear force-free optimization principle. We apply our code to model the solar corona using synoptic vector magnetograms as boundary condition. Below about two solar radii the plasma $\beta $ β and Alfvén Mach number $M_{A}$ M A are small and the magnetic field configuration of stationary MHD is basically identical to a nonlinear force-free field, whereas higher up in the corona (where $\beta $ β and $M_{A}$ M A are above unity) plasma and flow effects become important and stationary MHD and force-free configuration deviate significantly. The new method allows for the reconstruction of the coronal magnetic field further outwards than with potential field, nonlinear force-free or magnetostatic models. This way the model might help to provide the magnetic connectivity for joint observations of remote sensing and in-situ instruments on Solar Orbiter and Parker Solar Probe.

Magnetohydrodynamic Equilibria

Magnetohydrodynamic equilibria are time-independent solutions of the full magnetohydrodynamic (MHD) equations. An important class are static equilibria without plasma flow. They are described by the magnetohydrostatic equations j×B=∇p+ρ∇Ψ,∇×B=μ0j,∇·B=0. B is the magnetic field, j the electric current density, p the plasma pressure, ρ the mass density, Ψ the gravitational potential, and µ0 the permeability of free space. Under equilibrium conditions, the Lorentz force j×B is compensated by the plasma pressure gradient force and the gravity force. Despite the apparent simplicity of these equations, it is extremely difficult to find exact solutions due to their intrinsic nonlinearity. The problem is greatly simplified for effectively two-dimensional configurations with a translational or axial symmetry. The magnetohydrostatic (MHS) equations can then be transformed into a single nonlinear partial differential equation, the Grad–Shafranov equation. This approach is popular as a first approximation to model, for example, planetary magnetospheres, solar and stellar coronae, and astrophysical and fusion plasmas. For systems without symmetry, one has to solve the full equations in three dimensions, which requires numerically expensive computer programs. Boundary conditions for these systems can often be deduced from measurements. In several astrophysical plasmas (e.g., the solar corona), the magnetic pressure is orders of magnitudes higher than the plasma pressure, which allows a neglect of the plasma pressure in lowest order. If gravity is also negligible, Equation 1 then implies a force-free equilibrium in which the Lorentz force vanishes. Generalizations of MHS equilibria are stationary equilibria including a stationary plasma flow (e.g., stellar winds in astrophysics). It is also possible to compute MHD equilibria in rotating systems (e.g., rotating magnetospheres, rotating stellar coronae) by incorporating the centrifugal force. MHD equilibrium theory is useful for studying physical systems that slowly evolve in time. In this case, while one has an equilibrium at each time step, the configuration changes, often in response to temporal changes of the measured boundary conditions (e.g., the magnetic field of the Sun for modeling the corona) or of external sources (e.g., mass loading in planetary magnetospheres). Finally, MHD equilibria can be used as initial conditions for time-dependent MHD simulations. This article reviews the various analytical solutions and numerical techniques to compute MHD equilibria, as well as applications to the Sun, planetary magnetospheres, space, and laboratory plasmas.