Development of Anisotropic Magnetized Viscosity for Magnetized Liner Inertial Fusion Simulations in FLASH

This paper presents the first implementation of the full Braginskii magnetized viscosity tensor in the FLASH code for Magnetized Liner Inertial Fusion (MagLIF) simulations, demonstrating through verification and application that this anisotropic transport mechanism significantly damps vortical structures, mitigates Rayleigh-Taylor instabilities, and preserves fusion yield.

The Gist

Imagine you are trying to build a miniature star in a bottle to generate clean, infinite energy. This is the goal of Magnetized Liner Inertial Fusion (MagLIF).

Here's how it works in simple terms: You take a tiny cylinder of metal (the "liner"), fill it with fuel (hydrogen gas), and blast it with a massive magnetic field. Then, you squeeze the cylinder inward with incredible force, hoping to crush the fuel so hard that it fuses, releasing energy.

The problem? Nature hates being squeezed. As you crush the fuel, it gets messy. It swirls, it mixes, and it leaks heat. It's like trying to squeeze a water balloon; if you squeeze too fast or unevenly, the water splashes out, and the balloon doesn't pop the way you want.

The Problem: The "Slippery" Fuel
In the past, scientists simulated this process using computer models that treated the fuel like a simple, smooth fluid. They assumed that if the fuel started swirling or spinning (vortices), it would just keep spinning forever, wasting energy that could have been used for fusion. They also assumed that when the fuel got crushed, it would mix with the cold metal liner, cooling down the star before it could ignite.

But in reality, the fuel in these experiments is super-hot plasma (a gas so hot the electrons are stripped from atoms) and it is super-magnetic.

The Discovery: The "Magnetic Honey"
This paper introduces a new, crucial ingredient to the computer models: Magnetized Viscosity.

Think of viscosity as "thickness" or "stickiness." Honey has high viscosity; water has low viscosity.

• Without magnets: The fuel is like water. If you stir it, it swirls easily and mixes with the container walls.
• With strong magnets: The fuel acts like honey that only flows in one direction. The magnetic field acts like a set of invisible train tracks. The fuel can slide easily along the tracks, but it is incredibly sticky and resistant if it tries to move across them.

The authors of this paper built a new "engine" for their computer simulation (called FLASH) that understands this "magnetic honey" physics. Before this, the computer models were missing this rule, so they were predicting the fuel would be much messier and cooler than it actually is.

What Happened When They Turned It On?
The researchers ran two sets of simulations: one with the old "water-like" physics and one with the new "magnetic honey" physics. The results were dramatic:

1. The Swirls Stopped: In the old models, the fuel developed wild, chaotic swirls (vortices) that wasted energy. In the new models, the "magnetic honey" acted like a brake, smoothing out these swirls.
2. The Heat Stayed In: Those swirling motions didn't just disappear; the "stickiness" of the magnetic honey converted that wasted spinning energy into heat. It's like rubbing your hands together to warm them up; the friction (viscosity) turned motion into warmth. This made the fuel hotter, which is exactly what you need for fusion.
3. The "Hot Spot" Survived: When the fuel was crushed, the old models showed the cold metal liner mixing into the hot center, killing the reaction. The new models showed the "magnetic honey" acting as a shield, keeping the cold liner out and preserving a clean, hot core.
4. More Energy: Because the fuel stayed hotter and cleaner, the simulated fusion reaction produced significantly more energy. In one test, the energy output jumped by 134% compared to the old models!

The Big Picture
Think of this like upgrading a video game.

• Old Version: The physics engine was too simple. The characters (fuel particles) slid around too easily, crashed into walls, and the game ended in a mess.
• New Version: The developers added a new "friction" setting that accounts for the magnetic field. Suddenly, the characters move more realistically. They don't slide into the walls as easily, they generate heat when they rub against each other, and the game (the fusion reaction) actually works.

Why Does This Matter?
This paper proves that magnetized viscosity is not just a tiny detail; it is a major player in making fusion work. By adding this "magnetic honey" physics to their simulations, scientists can now predict how fusion reactors will behave much more accurately.

This is a huge step forward for companies like Pacific Fusion and others trying to build real fusion power plants. It means they can design better targets, squeeze the fuel more effectively, and get closer to the day when we have a clean, limitless energy source.

In a Nutshell:
They found that the fuel in fusion experiments acts like sticky, magnetic honey rather than watery fluid. By teaching their computers to understand this stickiness, they discovered that the fuel stays hotter, cleaner, and produces twice as much energy as previously thought. This brings us one step closer to building a star on Earth.

Technical Summary

1. Problem Statement

Magnetized Liner Inertial Fusion (MagLIF) is a promising pathway to controlled thermonuclear fusion, operating in a regime where strong magnetic fields (10–20 T) significantly alter plasma transport properties. While thermal conduction is well-studied in MagLIF, viscosity has historically been neglected or treated as an isotropic scalar.

However, in strongly magnetized plasmas, the viscous stress tensor becomes highly anisotropic due to ion gyromotion. The magnetic field suppresses cross-field momentum transport while leaving parallel transport largely unaffected. The absence of this physics in current simulations leads to:

• Overestimation of hydrodynamic instabilities (specifically Rayleigh-Taylor and Kelvin-Helmholtz).
• Inaccurate predictions of stagnation temperatures (discrepancies of 10–15% compared to experiments).
• Failure to account for the conversion of kinetic energy in vortical structures into thermal energy (viscous heating).

There was no prior implementation of the full Braginskii magnetized viscosity tensor for arbitrary magnetic field orientations within MagLIF simulation frameworks.

2. Methodology

The authors developed and implemented the full Braginskii viscosity tensor into the Pacific Fusion branch of the FLASH code, a multi-physics radiation-magnetohydrodynamics (MHD) solver.

• Mathematical Formulation:

  • The implementation uses the full Braginskii constitutive relation, which expresses the viscous stress tensor ($\Pi$) as a linear combination of five basis tensors ($W_0$ to $W_4$) weighted by five viscosity coefficients ($\eta_0$ to $\eta_4$).
  • These coefficients depend on local plasma state ($n_i, T_i$) and the ion Hall parameter ($\omega_i \tau_i$).
  • The model accounts for field-parallel compression, perpendicular shear, and gyroviscous effects (non-dissipative momentum transfer).
  • Only the dissipative coefficients ($\eta_0, \eta_1, \eta_2$) contribute to viscous heating ($Q_{visc} = -\Pi : \nabla v$); gyroviscous terms ($\eta_3, \eta_4$) are non-dissipative.

• Numerical Implementation:

  • Implicit Solver: To handle the extreme range of viscosity values (spanning orders of magnitude between cold liner material and hot fuel) without restrictive time-step constraints, the authors employed an implicit backward Euler scheme for velocity diffusion.
  • HYPRE Integration: The momentum equation was reformulated into a canonical form compatible with the HYPRE library's structured grid solvers. This required mapping the Braginskii tensor divergence into a sparse linear system involving rank-4 coefficient tensors.
  • Geometry: The implementation supports arbitrary magnetic field orientations within a 2D cylindrical (axisymmetric) framework.

3. Verification and Validation

The implementation was rigorously verified through four distinct test cases:

1. Aligned Field Diffusion: Compared against an analytical solution for velocity diffusion where the magnetic field and velocity gradients are parallel. Results showed excellent agreement in damping rates.
2. Out-of-Plane Field Limit: Verified that the general implementation correctly reduces to Braginskii's original analytic expressions for a purely azimuthal magnetic field.
3. Method of Manufactured Solutions (MMS): Applied to a fully 3D magnetic field topology to verify spatial and temporal convergence. The solver demonstrated second-order spatial accuracy and first-order temporal accuracy, confirming the discretization is correct.
4. Magnetized Viscous Shock Profiles: Compared against semi-analytic solutions for collision-dominated shocks. The simulation correctly reproduced the shock thickness narrowing caused by magnetization and the transition between weakly and strongly magnetized regimes.

4. Key Results

The authors applied the new module to two MagLIF-relevant configurations: a "Pool Heated" target and a traditional seeded-perturbation implosion.

• Damping of Vortical Structures:

  • In viscous simulations, magnetized viscosity effectively damps vortical flow structures (vorticity) that persist in inviscid runs.
  • This damping converts kinetic energy from rotational motions into thermal energy, raising the ion temperature ($T_i$) by 100–200 eV (relative increases of 50–100% in cooler regions).

• Stabilization of Instabilities:

  • In simulations with seeded Rayleigh-Taylor perturbations, the viscous model significantly suppressed the growth of deceleration-phase instabilities.
  • Hot Spot Integrity: In inviscid runs, the hot spot collapsed due to hydrodynamic mixing. In viscous runs, a coherent, low-density hot spot was preserved.
  • Temperature Retention: At stagnation, the central ion temperature in the viscous case was more than an order of magnitude higher than in the inviscid case.

• Fusion Yield Preservation:

  • Viscous simulations consistently produced higher fusion yields across all perturbation amplitudes.
  • Yield Enhancement: At the largest perturbation amplitude (40 $\mu$m), the viscous simulation preserved 134% of the yield compared to the inviscid case (which suffered catastrophic yield loss).
  • Alpha Heating Feedback: The yield improvement was amplified when alpha-particle deposition was enabled, creating a positive feedback loop where better confinement leads to higher temperatures and further viscosity effects.
  • Neutron-Averaged Metrics: Including viscosity increased the burn-region pressure by 35.7% and density by 25.1% at high perturbation amplitudes.

5. Significance and Conclusion

This work establishes magnetized viscosity as a non-negligible physical mechanism essential for predictive modeling of MagLIF and related magnetized ICF concepts.

• Scientific Impact: It resolves discrepancies between simulations and experiments regarding stagnation temperatures and instability growth. It demonstrates that anisotropic transport physics can fundamentally alter implosion dynamics by stabilizing instabilities and enhancing thermal energy via viscous dissipation.
• Engineering Impact: The validated capability allows for more accurate performance predictions for the Pacific Fusion Demonstration System (targeting 60 MA).
• Future Outlook: The authors plan to extend this capability to 3D simulations and conduct systematic parameter sweeps to optimize MagLIF target designs, leveraging the new understanding of how viscosity stabilizes high-energy-density plasmas.

In summary, the successful integration of the full Braginskii viscosity tensor into FLASH provides a critical tool for understanding and optimizing magnetized fusion implosions, revealing that viscosity acts as a stabilizing and heating mechanism that significantly improves fusion yield.