Development of a Reduced Multi-Fluid Equilibrium Model and Its Application to Proton-Boron Spherical Tokamaks

This paper develops a computationally efficient reduced multi-fluid model to capture rotation-induced centrifugal separation and electrostatic polarization in proton-boron spherical tokamaks, demonstrating that these effects significantly alter plasma equilibrium and must be accounted for in reactor design.

The Gist

The "Spinning Salad" Problem: Making Clean Energy with Proton-Boron Fusion

Imagine you are trying to build the world’s cleanest, most powerful engine. Instead of burning gasoline, this engine uses a "nuclear soup" made of protons and boron. This is called Proton-Boron (p-B11) fusion.

Unlike current nuclear technology, this "soup" doesn't produce messy, radioactive neutrons. It’s clean, safe, and potentially endless. But there is a massive catch: to get this soup to fuse, you have to make it incredibly hot—way hotter than anything we’ve ever achieved—and you have to hold it in place using massive, spinning magnetic fields.

This paper describes a new mathematical "map" (a model) designed to help scientists design the perfect container for this high-speed, high-heat soup.

---

The Problem: The "Centrifuge Effect"

To understand why we need this new model, imagine a salad spinner.

When you spin a salad spinner, the heavy bits—like chunks of cucumber or olives—get flung to the very edges of the bowl, while the light bits—like tiny droplets of dressing—stay closer to the center.

In a fusion reactor (specifically a "Spherical Tokamak"), we use magnetic fields to hold the plasma. To keep the plasma stable, we often make it spin incredibly fast.

• The Protons are like the light salad dressing.
• The Boron ions are like the heavy olives.

Because the boron is much heavier than the protons, the spinning motion flings the boron to the outer edges of the magnetic "bowl."

Why is this a disaster? For fusion to work, the protons and the boron need to be bumping into each other in the center to create energy. If the "olives" (boron) are all stuck at the edges and the "dressing" (protons) is in the middle, the engine won't start. It’s like trying to make a sandwich where the meat is in the kitchen and the bread is in the living room.

---

The Old Map vs. The New Map

Until now, scientists have mostly used a "Single-Fluid Model." Think of this like looking at a salad spinner through a blurry lens where you just see a "green blob." It assumes everything in the bowl moves together as one single liquid. This is easy to calculate, but it’s wrong for Proton-Boron fusion because it ignores the fact that the heavy and light parts are separating.

On the other hand, there are "Full Multi-Fluid Models." These are like looking at every single individual leaf of lettuce and every drop of oil. They are incredibly accurate, but they are so mathematically complex that even the world's fastest supercomputers struggle to solve them quickly. They are like trying to predict the movement of every single molecule in a hurricane—it takes too long to be useful for designing a machine.

The authors of this paper created a "Reduced Multi-Fluid Model."
Think of this as a High-Definition Map. It’s not quite as detailed as looking at every molecule, but it is much smarter than the "blurry blob" model. It specifically focuses on the two things that matter most:

1. The Spin (Centrifugal Force): How the heavy boron gets flung to the edges.
2. The Electric Charge (Electrostatic Potential): When the heavy boron moves to the edge, it leaves a "charge imbalance" behind, creating an internal electric field (like a giant battery inside the plasma).

---

What did they find?

The researchers tested their new "map" on two designs: a small experimental machine (EHL-2) and a massive future power plant (EHL-3B).

They discovered a "tipping point":

• Slow Spin: If the plasma spins slowly, the "salad" stays mixed. The old, simple models work fine.
• Fast Spin: Once the spin hits a certain speed, the separation becomes violent. In the large reactor design, the boron gets pushed so far to the edge that it forms a "crescent shape," leaving the center empty. They also found that this creates a massive electrical "wall" (about 10,000 volts!) inside the plasma.

Why does this matter?

If we want to build a real Proton-Boron power plant, we can't use the old, blurry maps—they would lead us to build a machine that doesn't work.

This new model gives engineers a "Goldilocks" tool: it is fast enough to use for designing reactors, but smart enough to realize that the "olives" are flying to the edges. It provides the blueprint needed to ensure the protons and boron stay close enough to shake hands and create clean, limitless energy.

Technical Summary

Technical Summary: Development of a Reduced Multi-Fluid Equilibrium Model for Proton-Boron Spherical Tokamaks

1. Problem Statement

Proton-Boron ($p\text{-}^{11}\text{B}$) fusion is a highly desirable aneutronic energy pathway, but it requires extreme ion temperatures ($>100\text{ keV}$) and robust magnetic confinement, typically provided by Spherical Tokamaks (STs). In these devices, high-power Neutral Beam Injection (NBI) drives significant toroidal rotation.

The research identifies a critical gap in existing modeling:

• Single-Fluid MHD Models: Fail to account for the extreme mass disparity between protons ($1m_p$) and boron ions ($\sim 11m_p$). Under rotation, centrifugal forces cause heavy ions to migrate toward the low-field side (LFS), creating species separation and electrostatic polarization that single-fluid models cannot capture.
• Comprehensive Multi-Fluid Models: While physically accurate, they include poloidal flows and pressure anisotropy, which introduce mathematical "stiffness" (transonic singularities) and excessive complexity, making them impractical for routine engineering and reactor design.

2. Methodology

The authors propose a "reduced" multi-fluid equilibrium model designed to balance physical fidelity with computational robustness.

• Core Assumptions:
  • Dominant Toroidal Rotation: The model prioritizes toroidal inertia while neglecting poloidal flow inertia to ensure the governing Grad-Shafranov operator remains strictly elliptic (avoiding numerical singularities).
  • Isothermal Surfaces: Each species' temperature is assumed to be a function of the magnetic flux ($\psi$).
  • Massless Electrons: Electron inertia is neglected, simplifying the centrifugal potential.
• Mathematical Framework:
  • The model utilizes a generalized Grad-Shafranov (GS) equation coupled with species-specific Bernoulli relations to describe density distributions.
  • Quasi-neutrality constraint: A self-consistent electrostatic potential ($\Phi$) is calculated at every spatial point to balance the charge separation caused by centrifugal forces.
• Numerical Implementation: The authors developed the NFEQ (N-Fluid EQuilibrium) code, which uses a nested iterative scheme: an outer Picard iteration loop for the magnetic flux ($\psi$) and an inner Newton-Raphson loop to solve for the electrostatic potential ($\Phi$). A feedback mechanism is also implemented to maintain a fixed total plasma current ($I_p$).

3. Key Contributions

• Development of a Tractable Model: Created a theoretical baseline that captures the primary physics of species separation and electrostatic coupling without the numerical instability of full multi-fluid models.
• New Governing Equations: Formulated a closed system of equations (Eq. 2–9 in the paper) that links magnetic geometry, species density, and electrostatic potential through a generalized pressure gradient.
• Validation via Engineering Scales: Applied the model to two specific, real-world ENN Group designs: the experimental EHL-2 and the reactor-scale EHL-3B.

4. Results

The simulations reveal that the impact of multi-fluid effects is governed by the ion Mach number ($M$):

• Low-Rotation Regime ($M < 0.5$): Multi-fluid effects are negligible; the solution converges to the standard single-fluid MHD limit.
• High-Rotation Regime ($M > 1$):
  • Centrifugal Separation: Heavy $^{11}\text{B}$ ions are driven toward the LFS, creating a "crescent-shaped" density distribution. This causes a spatial mismatch between protons and boron, potentially reducing local fusion power density in the core.
  • Electrostatic Polarization: In the reactor-scale EHL-3B, a massive internal electrostatic potential of approximately 10–15 kV is generated to maintain quasi-neutrality.
  • Profile Modification: The redistribution of pressure significantly modifies the safety factor ($q$) profile and the Shafranov shift.

5. Significance

• For $p\text{-}^{11}\text{B}$ Reactor Design: The study proves that multi-fluid modeling is a necessity, not an option, for designing $p\text{-}^{11}\text{B}$ reactors. Ignoring these effects would lead to inaccurate predictions of fusion power, stability, and transport.
• Broader Fusion Applications: The model is generic and can be applied to:
  • High-Z Impurity Transport: Predicting how heavy impurities (like Tungsten) accumulate in D-T tokamaks.
  • Isotopic Separation: Assessing fuel distribution in D-T or D-$^3\text{He}$ plasmas.
  • Advanced Tokamak (AT) Modes: Analyzing internal transport barriers (ITBs) in strongly rotating plasmas.