Effect of convective transport in edge/SOL plasmas of ADITYA-U tokamak

By integrating a custom limiter geometry into the UEDGE fluid transport code, this study demonstrates that modeling the electron density profiles of the ADITYA-U tokamak requires both a constant inward convective velocity and a perpendicular diffusion coefficient that lies between neoclassical and Bohm values.

The Gist

The Great Plasma Barrier: A Story of Heat, Wind, and Walls

Imagine you are trying to run a high-tech, ultra-hot furnace inside a building. This furnace is so powerful that if the heat touches the walls of the building, it will melt them instantly. To prevent this, you build a "buffer zone"—a specialized layer of air and shielding that sits between the intense core of the fire and the actual walls of the building.

In the world of fusion energy, scientists are trying to build a "star in a bottle" using a machine called a tokamak (in this case, one called ADITYA-U). The "star" is a ball of superheated gas called plasma. The problem? Plasma is incredibly hot and wants to escape. If it touches the machine's metal parts, it destroys them.

This paper is about how scientists used a sophisticated computer simulation to understand that "buffer zone" (called the Edge/SOL region) to make sure the "walls" stay safe.

---

1. The Problem: The Leaky Fence

Think of the plasma core as a crowded, high-energy dance floor. In the center, everyone is dancing in tight, organized circles (this is the core). But as you move toward the edge of the room, the crowd gets disorganized. Some people start stumbling toward the walls (the limiter).

Scientists want to know exactly how fast these "people" (particles) are moving toward the walls and how much "heat" they are carrying with them. If they can predict this, they can design better "fences" to protect the machine.

2. The Tool: The Digital Twin

To study this without melting a real machine, the researchers used a computer program called UEDGE.

Think of UEDGE as a highly advanced weather forecasting simulator. Just as meteorologists use math to predict how wind and heat move through a city, these scientists used UEDGE to predict how plasma particles and heat move through the edge of the tokamak.

However, UEDGE was originally designed for a different "room shape" (a divertor). The researchers had to act like digital architects, rewriting the code to create a custom "3D map" (a computational mesh) that accurately represented the specific shape of the ADITYA-U machine.

3. The Discovery: It’s Not Just a Leak; There’s a Wind!

This is the most important part of the paper.

Before this study, scientists thought the particles moved toward the walls mostly through diffusion.

• Analogy for Diffusion: Imagine a drop of ink spreading slowly in a glass of water. It moves from where there is a lot of ink to where there is very little. It’s a slow, passive process.

The researchers tried to model the plasma using only this "ink-spreading" logic, but the computer's prediction didn't match what they saw in the real experiments. The real plasma was moving toward the walls much more aggressively than the "ink" model suggested.

They realized they were missing something: Convection.

• Analogy for Convection: Instead of just ink spreading in still water, imagine a strong gust of wind blowing the ink across the room.

The researchers discovered that there is an inward "wind" (convective velocity) of about 1.5 meters per second. This "wind" pushes the particles toward the edge much faster than simple spreading would. By adding this "wind" into their computer simulation, the digital model finally matched the real-world measurements perfectly.

4. The Heat: The Hot Spot

Finally, the researchers looked at the "heat flux"—essentially how much "thermal energy" is being slammed into the machine's parts.

They found that the heat doesn't hit the walls evenly. It’s like a flashlight beam: the heat is most intense right at the "limiter tip" (the sharp edge of the physical barrier). They calculated exactly how much heat was being carried by "conduction" (the slow soak of heat) versus "convection" (the fast blast of heat).

Why does this matter?

If we want to use fusion to power our homes, we need to build massive machines like ITER. To do that, we can't just guess how the plasma behaves. We need to know exactly how much "wind" and "heat" the edge of the plasma will throw at our walls.

This paper proves that by using smart computer models and accounting for both "spreading" and "blowing," we can accurately predict the behavior of the plasma, helping us build safer, more powerful fusion reactors for the future.

Technical Summary

Technical Summary: Effect of Convective Transport in Edge/SOL Plasmas of ADITYA-U Tokamak

1. Problem Statement

Understanding the transport of particles and energy in the edge and Scrape-Off Layer (SOL) of a tokamak is critical for protecting the machine's material boundaries and maintaining core plasma stability. In the ADITYA-U tokamak, which utilizes a limiter configuration (rather than a divertor), the dynamics of the edge region are complex due to sharp gradients in density and temperature, open magnetic field lines, and interactions with material surfaces.

A specific challenge addressed in this study is that standard diffusive transport models (using only a perpendicular diffusion coefficient, $D_\perp$) often fail to accurately reproduce experimental electron density ($n_e$) profiles in the SOL. There is a need to determine whether non-diffusive transport mechanisms, specifically convective velocity ($v_{conv}$), are required to achieve a high-fidelity match between simulation and experimental Langmuir probe measurements.

2. Methodology

The researchers employed a multi-step computational and comparative approach:

• Code Modification: The study utilized UEDGE, a 2-D edge plasma fluid transport code. Since UEDGE is natively designed for divertor geometries (with X-points), the authors developed a custom MATLAB/PYTHON-based mesh-generator. This routine uses plasma equilibrium flux surfaces (from the IPREQ code) to create a non-uniform computational mesh tailored to the circular limiter geometry of ADITYA-U.
• Mesh Design: To capture sharp gradients near the limiter tip, the mesh was designed with high resolution: a minimum poloidal cell arm of $\sim 5\text{ mm}$ and a radial arm of $\sim 1\text{ mm}$ near the limiter surfaces.
• Simulation Parameters: The UEDGE code solved fluid equations for particle continuity, ion/electron temperatures, and momentum. The researchers tested various combinations of:
  • Perpendicular Diffusion ($D_\perp$): Tested as a radially constant value.
  • Convective Velocity ($v_{conv}$): Tested in various forms, including constant, step functions, and radially increasing profiles.
• Validation: The simulated radial profiles of $n_e$ and $T_e$ were compared against experimental data obtained from Langmuir probe arrays during a typical ADITYA-U discharge ($I_p \approx 120\text{ kA}$).

3. Key Contributions

• Geometric Adaptation of UEDGE: Successfully integrated a custom mesh-generation routine to allow a divertor-centric code (UEDGE) to simulate a limiter-based tokamak configuration.
• Identification of Convective Necessity: Demonstrated that diffusive transport alone is insufficient to model the $n_e$ profile in the SOL of ADITYA-U.
• Transport Coefficient Characterization: Provided a quantitative comparison between the simulated $D_\perp$ and various theoretical transport regimes (Neoclassical, Bohm, ITG, and Resistive Ballooning modes).

4. Results

• Optimal Transport Parameters: The best match with experimental electron density profiles was achieved using a radially constant inward convective velocity ($v_{conv}$) of $1.5\text{ m/s}$ and a radially constant perpendicular diffusion coefficient ($D_\perp$) of $\sim 0.2\text{ m}^2\text{/s}$.
• Comparison with Theory:
  • The obtained $D_\perp$ ($\sim 0.2\text{ m}^2\text{/s}$) is higher than the neoclassical estimate ($\sim 0.02\text{ m}^2\text{/s}$).
  • However, it is significantly lower than the Bohm ($\sim 1\text{ m}^2\text{/s}$), ITG ($\sim 1.5\text{ m}^2\text{/s}$), and Resistive Ballooning ($\sim 320\text{ m}^2\text{/s}$) estimates. This suggests that fluctuation-induced transport is not the dominant driver of particle diffusivity in this specific edge region.
• Heat Flux Dynamics:
  • The maximum radial electron conductive heat flux was $\sim 2.1 \times 10^{21}\text{ eV m}^{-2}\text{ s}^{-1}$.
  • The maximum radial electron convective heat flux was $\sim 1.3 \times 10^{21}\text{ eV m}^{-2}\text{ s}^{-1}$.
  • While conduction dominates in the core-edge region, the total radial electron heat flux reaches its maximum near the limiter tip location, where convective processes become highly significant.
• Convection Nature: The study found that the $n_e$ profile is highly sensitive to the convective velocity within the core-edge region ($r < r_{limiter-tip}$), but less sensitive to the specific magnitude of $v_{conv}$ once it exceeds $1.5\text{ m/s}$ in the SOL.

5. Significance

This work provides a robust framework for modeling the edge plasma of limiter-configured tokamaks, which is relevant for the early operational phases of ITER (which will use temporary limiters). By proving that convective transport must be included to accurately predict density profiles, the study improves the predictive capability of fluid codes for plasma-wall interactions. Furthermore, the successful integration of the mesh generator paves the way for coupling UEDGE with neutral transport codes (like DEGAS2) to study complex neutral-plasma dynamics.