🔬 physics

Parameter Scan of Multi-Fluid Equilibria in Rotating p-11B Plasmas: Effects on Fusion Power and Bremsstrahlung Losses

This paper introduces VEQ-MF, a fast spectral framework for computing two-dimensional multi-fluid equilibria with species-dependent rotation, and applies it to rotating p-11B spherical tokamaks to reveal how competing centrifugal and polarization effects influence the fusion-to-bremsstrahlung power ratio.

Original authors: Xingyu Li, Huasheng Xie, Lai Wei, Zhengxiong Wang

Published 2026-07-17

📖 4 min read☕ Coffee break read

CC BY 4.0

Original authors: Xingyu Li, Huasheng Xie, Lai Wei, Zhengxiong Wang

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine trying to cook a meal so hot that the ingredients themselves start to glow with energy, but the heat keeps escaping faster than you can cook. This is the daily struggle of scientists trying to build a fusion reactor, a machine that mimics the sun to create clean, limitless energy. The specific recipe they are testing here involves smashing protons (hydrogen nuclei) into boron atoms. It's a "neutron-free" dream because, unlike other fusion recipes, it doesn't spit out dangerous radioactive particles. However, there's a catch: to make these atoms fuse, they need to be incredibly hot, and at those temperatures, they start screaming out energy in the form of light (radiation) faster than they can fuse. It's like trying to fill a bucket with a hose while someone is drilling a hole in the bottom; the water (energy) leaks out as fast as you pour it in.

The key to winning this game isn't just making things hotter; it's about how the ingredients move. In a normal pot, everything jiggles randomly. But in a fusion reactor, scientists can spin the plasma (the super-hot soup of atoms) like a giant centrifuge. When you spin something heavy, it gets pushed to the outside. In this specific recipe, the boron atoms are much heavier than the protons. If you spin them all at the same speed, the heavy boron gets flung to the outer edge, bunching up and making the plasma glow even brighter, which wastes more energy. The big question this paper asks is: What happens if we spin the heavy boron and the light protons at different speeds? Could we trick the plasma into keeping its energy better?

This paper, written by a team of physicists, acts like a high-speed video game simulator to answer that question. They built a digital model called VEQ-MF to test a specific type of fusion reactor design (spherical tokamaks) filled with protons and boron. Instead of building a real machine, they ran thousands of virtual experiments, changing the rotation speeds of the protons and boron to see how the energy balance shifted. They found that if you spin the heavy boron slower than the light protons, two cool things happen. First, the heavy boron doesn't pile up as much on the outside, which stops the plasma from glowing as brightly and wasting energy. Second, the difference in speed between the two types of atoms creates a "drift" that actually helps them smash together more effectively, boosting the fusion power.

The results are exciting but come with a big "simulated only" warning. In their computer models, specifically for a larger reactor design called EHL-3B, they found a "sweet spot" where the fusion power could actually beat the energy lost to radiation alone. In the standard spinning scenario, the fusion power was only about 18% of the energy lost to radiation. But by spinning the protons and boron at different speeds, they simulated a scenario where fusion power rose to about 130% of the radiation loss. It's like finding a way to patch the hole in the bucket while simultaneously turning up the hose.

However, the authors are very careful not to call this a solved problem. They point out that their main calculation of this ratio ignores the massive amount of energy required to keep the protons and boron spinning at different speeds in the first place. In fact, they performed a separate, rough "order-of-magnitude" estimate (detailed in the paper's Appendix C) suggesting that the friction between the spinning layers might require about 5.2 gigawatts of power to sustain—far more than the 120 megawatts of fusion power they generated. This friction cost wasn't included in the main ratio because the model focused on the equilibrium state itself, but the estimate serves as a reality check. So, while the paper proves that this "differential rotation" trick works beautifully in a computer simulation to improve the balance between fusion and radiation, it doesn't mean we can build a working reactor tomorrow. The energy cost to maintain the spin currently dwarfs the energy gained. It's a promising clue for future scientists, showing that if we can figure out how to spin the ingredients differently without spending a fortune on energy, we might finally get fusion to work.

Technical Summary: Parameter Scan of Multi-Fluid Equilibria in Rotating p-11B Plasmas

Problem Statement
The proton-boron-11 (p-11B) fusion reaction is attractive for its neutron-lean nature but faces significant challenges due to the high ion temperatures (hundreds of keV) required to overcome the Coulomb barrier. At these temperatures, relativistic bremsstrahlung radiation becomes a dominant energy loss channel, often exceeding fusion power in conventional thermal scenarios. Furthermore, maintaining non-equilibrium plasma states to mitigate these losses incurs momentum-exchange power costs. While single-fluid models are insufficient for p-11B plasmas due to the large mass disparity between protons and boron ions ( $m_B \approx 11m_p$ ), existing multi-fluid equilibrium models often lack the capability to efficiently scan complex parameter spaces involving species-dependent toroidal rotation. This paper addresses the need for a robust, fast spectral framework to analyze how prescribed differential rotation affects density redistribution, bremsstrahlung losses, and fusion power in rotating spherical tokamak (ST) configurations.

Methodology
The authors developed VEQ-MF, a fast spectral parameter-scan framework for two-dimensional axisymmetric multi-fluid equilibria with prescribed species-dependent toroidal rotation. The methodology integrates three core components:

Governing Equations: The model couples generalized Boltzmann density responses (accounting for centrifugal forces), quasineutral electrostatic polarization, and a generalized Grad–Shafranov (GS) equation. It assumes dominant toroidal rotation, flux-surface isothermality, and a massless electron limit, neglecting poloidal flow inertia to avoid transonic singularities.
Numerical Solver: VEQ-MF extends reduced-parameter GS representations and the VEQ spectral solver. It employs a two-layer nonlinear structure: an inner loop solving local quasi-neutrality via a vectorized Newton–Raphson method, and an outer loop solving the global GS operator using shifted Chebyshev polynomials and a hybrid Broyden method. Homotopy continuation is used to incrementally increase rotation frequencies.
Power Evaluation: Post-processing calculates fusion power ( $P_{fus}$ ) using a drift-Maxwellian reaction-rate coefficient ( $\langle \sigma v \rangle_{DM}$ ) to account for relative toroidal flow, and bremsstrahlung power ( $P_{brem}$ ) using a relativistic analytical radiation model. The primary metric is the bremsstrahlung-only power ratio $R_{fb} = P_{fus}/P_{brem}$ .
Validation: The solver was benchmarked against a high-resolution finite-difference (FDM) solver and subjected to convergence checks regarding spectral parameters and grid refinement.

Key Contributions

VEQ-MF Framework: The introduction of a spectral solver capable of handling coupled multi-species rotating equilibria with high computational efficiency (0.032–0.050 s per point for nonzero cases in MATLAB).
Diagnostic Decomposition: A novel method to separate the contributions of density redistribution (geometric overlap) and the drift-induced reaction-rate coefficient to the total fusion power.
Parameter Scans: Comprehensive scans of the $(\Omega_p, \Omega_B)$ plane for two ST geometries, EHL-2 ( $R_0=1.1$ m) and EHL-3B ( $R_0=3.2$ m), identifying regimes where differential rotation alters the fusion-to-bremsstrahlung balance.

Results

Iso-Rotation ( $\Omega_p = \Omega_B$ ): Under iso-rotation, the heavy boron ions accumulate strongly on the low-field side (LFS) due to centrifugal forces. This induces an electrostatic polarization that shifts electrons outward, increasing the volume-integrated $n_e^2$ . Consequently, bremsstrahlung losses rise faster than fusion power, causing $R_{fb}$ to decrease significantly with increasing rotation (e.g., dropping from ~0.40 to <0.10 in EHL-2).
Differential Rotation ( $\Omega_p \neq \Omega_B$ ):
- Bremsstrahlung Reduction: Reducing the boron rotation frequency ( $\Omega_B$ ) relative to protons ( $\Omega_p$ ) weakens the centrifugal polarization. This flattens the electron density profile, reducing the $n_e^2$ integral and lowering $P_{brem}$ .
- Fusion Power Recovery: In the larger EHL-3B geometry, the relative toroidal flow generates significant drift kinetic energy ( $E_d$ ). This shifts the relative velocity distribution, increasing the drift-Maxwellian reaction-rate coefficient $\langle \sigma v \rangle_{DM}$ .
- Competing Effects: A diagnostic decomposition reveals two competing trends as differential rotation increases: (1) the proton-boron density overlap integral decreases monotonically (reducing fusion potential), and (2) the density-overlap-weighted effective reaction-rate coefficient rises steadily. In EHL-3B, the rise in the reaction-rate coefficient eventually offsets the density penalty, allowing $P_{fus}$ to recover from a minimum.
Power Ratio Crossing: In EHL-3B, the combined effect of reduced bremsstrahlung and enhanced fusion reactivity drives the bremsstrahlung-only power ratio $R_{fb}$ from 0.18 (at iso-rotation) to approximately 1.3 at the largest imposed relative toroidal flow ( $\gamma_{Bp} \approx 0$ ). In EHL-2, the effect is limited to bremsstrahlung reduction due to insufficient drift energy to significantly alter $\langle \sigma v \rangle_{DM}$ .
Numerical Stability: Benchmarks show global differences of 1.7–3.4% between VEQ-MF and FDM for stored energy, bremsstrahlung, and fusion power. Convergence checks indicate sub-percent sensitivity to spectral and grid refinements.

Significance and Claims
The paper explicitly states that its purpose is not to assess reactor-level feasibility for p-11B fusion. Instead, it serves as a demanding model case to test a parameter-scan multi-fluid equilibrium calculation. The significance lies in providing a benchmarked equilibrium-scan workflow that successfully separates the density-overlap and reaction-rate-coefficient contributions to fusion power in rotating multi-fluid plasmas.

The authors emphasize that the reported $R_{fb} = 1$ crossing is a result of a "reduced bremsstrahlung-only power ratio" under prescribed profiles. It does not account for auxiliary power, transport losses, or the substantial power cost (estimated at ~5.2 GW for momentum exchange in the largest case) required to sustain the differential rotation. Therefore, the results characterize the equilibrium consequences of imposed differential rotation rather than the transport-limited sustainment of such states. The work establishes a foundation for future studies that would couple these equilibrium calculations with transport-level analyses of profile evolution and energetic costs.

Technical Summary: Parameter Scan of Multi-Fluid Equilibria in Rotating p-11B Plasmas

More like this