Axisymmetric reduction of the adjoint operator in ideal MHD equilibrium

This paper presents adjoint formulations for ideal magnetohydrodynamic equilibrium to enable sensitivity analysis and optimization, demonstrating that the axisymmetric reduction of the three-dimensional adjoint operator coincides with the perturbed Grad-Shafranov operator, thereby unifying the adjoint frameworks for both axisymmetric and fully three-dimensional equilibria.

The Gist

Imagine you are trying to balance a complex, floating sculpture made of magnetic fields and hot gas (plasma) inside a fusion reactor. To keep it stable, you need to solve a massive set of math equations. But sometimes, you need to know: "If I tweak this one tiny knob on the outside, how does the whole sculpture wobble?"

Doing this calculation for every single knob is like trying to count every grain of sand on a beach one by one. It takes forever.

This paper introduces a "shortcut" method called the Adjoint Method. Think of the Adjoint Method as a special mirror. Instead of pushing the sculpture from every angle to see how it reacts, you look in the mirror. The mirror instantly tells you how the sculpture would react to any push, all at once. This is incredibly useful for designing better reactors or figuring out what's happening inside a plasma just by looking at data.

Here is the core story of the paper, broken down into simple concepts:

1. The Two Different Worlds

For decades, scientists have treated plasma in two different ways:

• The Simple World (Axisymmetric): Imagine a perfect donut (a tokamak). It looks the same no matter how you spin it around its center. In this world, the math is simpler and well-understood.
• The Complex World (3D): Imagine a twisted pretzel or a star-shaped object (a stellarator). It looks different from every angle. The math here is much harder and requires powerful computers.

The problem was that scientists had a great "mirror" (the Adjoint Method) for the Simple World, and they were building a new "mirror" for the Complex World. But nobody had proven that the Complex Mirror was actually just the Simple Mirror stretched out and twisted. They were worried the two mirrors might show different reflections.

2. The Big Discovery: The Mirrors Match

The authors of this paper did the heavy lifting to prove that the Complex Mirror is exactly the Simple Mirror in disguise.

They started with the most complex, 3D version of the equations (the "parent" equations) and built the Adjoint Method from scratch. Then, they forced the math to behave like a perfect donut (imposing axisymmetry).

The Result: When they did this, the complex 3D math collapsed perfectly into the simple, known math for the donut shape. It wasn't just similar; it was identical.

3. The "Weighted" Scale Analogy

There was a small snag in how the math was written. In the simple donut world, the math uses a special "scale" or "weight" to measure things correctly. This weight is based on the distance from the center of the donut (mathematically, a factor of $1/R^2$).

• The Old Way: Some previous studies ignored this weight to make the math look cleaner. It worked for practical calculations, but it hid the deep connection between the simple and complex worlds. It was like weighing an apple in pounds and an orange in kilograms and saying, "They are the same weight," without converting the units.
• The New Way: This paper insists on keeping that "weight" (the $1/R^2$ factor). By doing so, they showed that the 3D mirror and the 2D mirror are structurally identical. The "weight" isn't an artificial rule; it naturally comes from how magnetic energy is stored in a donut shape.

4. Why This Matters (According to the Paper)

The paper doesn't claim to have built a new reactor or cured a disease. Instead, it provides a theoretical bridge.

• Confidence: It proves that the advanced 3D tools scientists are building are consistent with the trusted 2D tools they've used for years.
• The "Weak" Twist: The authors suggest this bridge is especially helpful for "quasi-axisymmetric" designs. These are reactors that are mostly perfect donuts but have tiny, weak twists to them. Because the math is now proven to connect, scientists can use the simple 3D mirror to understand how those tiny twists affect the plasma, without having to start from scratch.

Summary

Think of the paper as a translator. It took the complex, 3D language of plasma physics and showed that when you speak it in a "donut accent" (axisymmetry), it sounds exactly like the simple, 2D language we already know. This confirms that the shortcuts (Adjoint methods) used for simple shapes are mathematically valid foundations for the complex shapes of the future.

Technical Summary

Technical Summary: Axisymmetric Reduction of the Adjoint Operator in Ideal MHD Equilibrium

Problem Statement
Ideal magnetohydrodynamic (MHD) equilibrium calculations are fundamental to the design of magnetic confinement devices and the reconstruction of plasma states from diagnostic data. While axisymmetric equilibria (e.g., tokamaks) are governed by the Grad–Shafranov equation, three-dimensional (3D) equilibria (e.g., stellarators) require numerical solvers. Both regimes increasingly utilize optimization-based approaches where sensitivity analysis with respect to numerous parameters is required. The adjoint method offers an efficient framework for such sensitivity analysis. However, a precise mathematical relationship between the adjoint formulation of the full 3D equilibrium problem and the adjoint of the axisymmetric Grad–Shafranov equation has not been explicitly demonstrated. Specifically, it remains unclear how the 3D adjoint operator reduces under the assumption of axisymmetry and whether the resulting structure aligns with established axisymmetric adjoint formulations, particularly regarding the geometric weighting factors inherent to axisymmetric magnetic energy.

Methodology
The authors employ an operator-theoretic approach to derive the adjoint of the ideal MHD equilibrium equations without initially assuming symmetry.

1. General Adjoint Formulation: The MHD equilibrium is treated as a constrained optimization problem where the state variables are the current density ($J$), magnetic field ($B$), and pressure ($p$). A Lagrangian is constructed using the equilibrium constraints and an objective functional.
2. 3D Adjoint Derivation: By taking variations of the Lagrangian with respect to the state variables and applying the divergence theorem under fixed-boundary conditions ($B \cdot n = 0$), the authors derive the 3D adjoint equations. The adjoint variables ($\alpha, \beta, \gamma$) are introduced as Lagrange multipliers dual to the static equilibrium constraints, distinct from physical Lagrangian displacements.
3. Axisymmetric Reduction: The authors impose axisymmetry on the derived 3D adjoint system. They utilize standard representations for divergence-free vector fields in toroidal geometry, expressing the magnetic field and adjoint field $\alpha$ in terms of poloidal flux functions ($\psi$) and toroidal angles ($\phi$).
4. Projection and Simplification: The reduced system is projected along and across the toroidal direction to isolate scalar components. The authors specifically analyze the scalar quantity $\eta = -\alpha \cdot \nabla \psi$, interpreting it as a flux-like perturbation.
5. Operator Comparison: The resulting scalar equation is compared against the perturbed Grad–Shafranov operator. The authors also examine the inner product structure, specifically the role of the geometric weight $1/R^2$ arising from the variation of poloidal magnetic energy.

Key Contributions and Results

• Explicit Reduction: The paper demonstrates that the axisymmetric reduction of the full 3D homogeneous adjoint system yields a scalar operator that coincides exactly with the perturbed Grad–Shafranov operator (Zakharov & Pletzer 1999). This operator governs infinitesimal variations of the equilibrium flux function.
• Variable Correspondence: The study establishes a direct correspondence between the 3D adjoint variables and the reduced scalar variable. The reduced scalar adjoint variable is identified as $\eta = -\alpha \cdot \nabla \psi$, where $\alpha$ is the incompressible adjoint field associated with the force-balance equation.
• Weighted Inner Product Structure: The analysis clarifies that the reduced scalar operator is formally self-adjoint only with respect to a weighted inner product $(u, v)_w = \int u v R^{-2} dV$. This weight factor $1/R^2$ is shown to arise naturally from the axisymmetric expression of poloidal magnetic energy, rather than being an arbitrary imposition.
• Structural Consistency: The work resolves discrepancies with previous formulations (e.g., Sanpei et al. 2021, 2023) that omitted the geometric weight in the inner product. The authors show that while unweighted formulations are practically viable for numerical optimization, the weighted formulation is necessary to maintain structural compatibility and self-adjointness between the 3D adjoint and its axisymmetric reduction.

Significance
The paper provides a rigorous operator-level consistency statement linking the full 3D adjoint formulation to the axisymmetric Grad–Shafranov adjoint problem. By clarifying how the axisymmetric adjoint structure is embedded within the 3D framework, the work offers a unified theoretical basis for sensitivity analysis. This is particularly relevant for:

• Weakly Non-Axisymmetric Configurations: The 3D adjoint formulation serves as a "parent" framework for analyzing equilibria that are small perturbations of an axisymmetric state, allowing for the systematic organization of axisymmetric responses and 3D corrections.
• Adjoint-Based Reconstruction and Optimization: The explicit correspondence between 3D and axisymmetric variables aids in interpreting sensitivity results and may assist in bifurcation analysis or characterizing singular responses in linearized equilibrium systems.

The authors conclude that this reduction clarifies the mathematical relationship between 3D and axisymmetric ideal MHD adjoint problems, ensuring that the structural properties (such as self-adjointness) are preserved across dimensions when the appropriate geometric weights are applied.