Gyrokinetic simulation of the effect of transient fueling on plasma turbulence in ADITYA-U tokamak

Global electrostatic gyrokinetic simulations of the ADITYA-U tokamak reveal that transient gas puffing flattens the radial density profile, thereby suppressing trapped electron mode turbulence and enhancing core temperature and energy confinement time as an active control mechanism.

The Gist

The Big Picture: Taming the "Boiling Pot"

Imagine a tokamak (a device used to create fusion energy) as a giant, super-hot pot of soup. To get the soup to cook (fusion), you need to keep it extremely hot and dense in the center.

However, there's a problem: the soup is constantly "boiling" in a chaotic way. In physics terms, this is called plasma turbulence. This turbulence acts like a leaky lid, letting all the heat and energy escape from the center to the edges. If the heat escapes too fast, the fusion reaction dies out.

The scientists in this paper wanted to figure out how to stop this "leaking" and keep the heat trapped inside.

The Experiment: The "Quick Spritz"

The researchers used a machine called ADITYA-U (a tokamak in India). They tried a clever trick: instead of just adding more fuel slowly, they gave the plasma a short, sharp burst of gas (like a quick spritz from a spray bottle) at the edge of the machine.

What happened?

1. The Density Changed: The gas puff didn't just make the whole pot denser. It specifically made the middle part of the plasma "flatter." Imagine the density profile was like a steep mountain peak; the gas puff turned it into a gentle hill.
2. The Temperature Rose: Surprisingly, right after this gas spritz, the temperature in the very center of the plasma shot up.
3. The Result: The machine held onto its energy much better.

The Secret Sauce: Stopping the "Trapped Electron"

Why did this happen? The scientists used a super-computer simulation (a digital twin of the experiment) to look under the hood. They found the culprit was a specific type of chaos called the Trapped Electron Mode (TEM).

The Analogy: The Bouncing Ball
Think of the electrons in the plasma as tiny, hyperactive balls bouncing around inside a hallway.

• Before the gas puff: The hallway has steep walls (a steep density gradient). The balls get trapped in the corners, bouncing back and forth wildly. This bouncing creates a lot of noise and chaos (turbulence), which pushes heat out of the center.
• After the gas puff: The gas puff smoothed out the walls of the hallway. Now, the balls don't get trapped in the corners as easily. They stop bouncing chaotically.
• The Result: With less chaotic bouncing, the "noise" stops. The heat stays in the center, and the temperature rises.

The Simulation: The Digital Detective Work

The researchers used a powerful code called GTC (Gyrokinetic Toroidal Code) to simulate this. It's like running a video game where they can freeze time and watch how the invisible particles move.

• What they saw: Before the gas, the "chaos waves" (turbulence) were strong and spread all the way to the center, stealing the heat.
• After the gas: The simulation showed that the "chaos waves" were pushed out of the center. The core became a calm, quiet zone where heat could stay trapped. The turbulence was essentially "expelled" from the middle of the machine.

The Takeaway: A New Way to Control Fusion

The most exciting part of this paper is the conclusion: Gas puffing isn't just about adding fuel; it's a control knob.

Usually, scientists think of adding gas just to get more particles. But this study shows that by timing a short burst of gas perfectly, you can actually smooth out the plasma's shape to stop the turbulence.

In simple terms:
Imagine you are trying to keep a campfire hot.

• Old way: You just keep throwing more wood on it (adding fuel), but the wind (turbulence) keeps blowing the heat away.
• New way (this paper): You use a quick puff of air to arrange the logs so they sit tighter together. This stops the wind from blowing through the gaps. The fire gets hotter, not because you added more wood, but because you stopped the heat from escaping.

This discovery gives scientists a new tool to build better fusion reactors in the future, helping us get closer to clean, limitless energy.

Technical Summary

1. Problem Statement

Achieving enhanced energy confinement in tokamaks is limited by gradient-driven microturbulence, which causes anomalous particle and heat transport. Traditional methods to mitigate this often involve impurity injection (altering collisionality and dilution) or edge heating/current drive. However, these methods can compromise core plasma purity or lack precise temporal control.

The specific challenge addressed in this study is understanding the causal link between transient neutral gas puffing and the suppression of micro-instabilities in the ADITYA-U tokamak. While experiments showed that short gas puffs increase core temperature and confinement time, the underlying physical mechanism—specifically how a peripheral gas injection modifies core turbulence—remained to be fully quantified and theoretically explained.

2. Methodology

The study employs a combined approach of experimental data analysis and advanced numerical simulation:

• Experimental Data:

  • Device: ADITYA-U tokamak (Shot #36136).
  • Technique: Periodic hydrogen gas puffing (short pulses, ~1 ms duration, injecting $\sim 10^{17}$ molecules).
  • Diagnostics: Microwave interferometry for radial density profiles and Soft X-ray (SXR) chords combined with Langmuir probes for electron temperature ($T_e$) reconstruction.
  • Observation Window: Analysis focused on the current flat-top phase (133–155 ms), specifically comparing conditions before and after gas injection.

• Numerical Simulation:

  • Code: GTC (Gyrokinetic Toroidal Code), a global electrostatic gyrokinetic simulation code.
  • Inputs: Realistic equilibrium profiles (density and temperature) derived from experimental data for both "Before" and "After" gas puff scenarios.
  • Physics Model:
    • Kinetic treatment of both passing and trapped electrons.
    • Global simulation domain spanning flux coordinates $\psi/\psi_X \in [0.1, 0.95]$.
    • Resolution: 120 flux surfaces, 4000 poloidal grid points, 32 parallel grid points, and 50 marker particles per cell.
  • Analysis: Linear regime analysis to identify dominant instabilities and nonlinear turbulent regime analysis to calculate transport coefficients (diffusivity and heat flux).

3. Key Contributions

• Mechanism Identification: The paper establishes that transient gas puffing acts as an active control mechanism for microturbulence, not merely a particle source. It demonstrates that flattening the core density profile directly suppresses the drive for Trapped Electron Modes (TEM).
• Profile-Turbulence Coupling: It quantifies how a localized density modification in the mid-radius ($0.3 < \rho < 0.8$) propagates to suppress turbulence in the core, shifting the turbulence spectrum from electron-scale (TEM) to more benign ion-scale regimes.
• Validation of GTC: The study successfully validates the GTC code's ability to reproduce experimental observations (core temperature rise and confinement improvement) by modeling the specific profile changes induced by gas puffing in ADITYA-U.

4. Key Results

Experimental Observations

• Temporal Evolution: Following a gas puff, the central line-averaged density rises rapidly (~1 ms). Subsequently, the core electron temperature ($T_e$) increases, peaking 3–4 ms after injection.
• Profile Modification: The gas puff causes a significant flattening of the radial density profile in the mid-radius region. The core density ($\rho < 0.3$) remains relatively constant, while the gradient ($\kappa_n$) in the mid-radius decreases significantly.
• Confinement Improvement: This profile modification correlates with a rise in core temperature and an increase in the energy confinement time, which exceeds the timescale of the temperature rise, indicating a reduction in transport losses.

Simulation Findings

• Dominant Instability: Both before and after gas puffing, the Trapped Electron Mode (TEM) is identified as the dominant instability (over Ion Temperature Gradient or ITG modes), confirmed by the sign of the diamagnetic frequency.
• Suppression Mechanism:
  • Before Puff: Strong TEM-driven turbulence exists with high growth rates ($\gamma \approx 0.82 C_s/R_0$). The mode structure is localized closer to the core ($\psi/\psi_X \approx 0.5$).
  • After Puff: The flattened density profile reduces the density gradient drive. This leads to:
    1. A slight reduction in linear growth rates ($\gamma \approx 0.79 C_s/R_0$).
    2. A radial shift of the mode structure outward (peak moves to $\psi/\psi_X \approx 0.75$), effectively expelling the TEM from the core.
    3. A reduction in the poloidal mode number (from $m \approx 160$ to $120$).
• Transport Reduction:
  • The nonlinear simulations show a drastic reduction in turbulent transport.
  • Electron Diffusivity ($D_e$): Reduced by ~84%.
  • Ion Diffusivity ($D_i$): Reduced by ~94%.
  • Heat Transport: The electron heat diffusivity ($\chi_e$) is reduced by an order of magnitude, directly explaining the observed core temperature rise.
• Propagation: The inward propagation of turbulence fluctuations from the edge to the core is curtailed after gas puffing.

5. Significance

This work provides a critical link between experimental control actuators (gas puffing) and fundamental plasma physics (gyrokinetic turbulence).

• Active Control Strategy: It proves that neutral fueling can be used as a sophisticated tool to actively regulate turbulence and optimize confinement, rather than just refueling the plasma.
• Predictive Capability: The self-consistent picture provided by the simulations offers a predictive model for how profile shaping influences transport, which is essential for the design of future fusion reactors (like ITER or DEMO) where active control of turbulence will be vital for sustained high-performance operation.
• Scientific Insight: It clarifies that the suppression of TEM turbulence via density profile flattening is a viable pathway to enhance core confinement, offering a new "recipe" for transport mitigation in magnetic confinement fusion.