The impact of plasma turbulence on atomic reaction rates in detached divertors

Basic Information

• Paper ID: 2510.09579
• Title: The impact of plasma turbulence on atomic reaction rates in detached divertors
• Authors: K. Eder, W. Zholobenko, A. Stegmeir, K. Zhang, F. Jenko (Max-Planck-Institut für Plasmaphysik)
• Classification: physics.plasm-ph (Plasma Physics)
• Publication Date: October 13, 2025
• Paper Link: https://arxiv.org/abs/2510.09579

Abstract

This study analyzes the impact of plasma turbulence fluctuations on atomic reaction rates in detached divertors. Through GRILLIX edge turbulence simulations, differences in ionization, recombination, and radiation rates between turbulent and mean-field conditions were investigated. Under attached conditions with low fluctuation amplitudes (<50%), the impact is negligible; however, under detached conditions with high fluctuation amplitudes (>300%), significant differences emerge. Ionization and radiation rates obtained from turbulent inputs are approximately 2 times lower than those from mean-field inputs, attributed to the anti-correlation between electron density and temperature under detached conditions. When this correlation is eliminated, turbulent rates recover to mean-field levels; when arranged as positive correlation, turbulent ionization rates increase by approximately 3 times.

Research Background and Motivation

Problem Definition

Numerical plasma models of the edge and scrape-off layer (SOL) rely on reaction rates to describe critical atomic processes such as ionization, recombination, charge exchange, and line radiation. These rates exhibit nonlinear dependence on local density and temperature, making them sensitive to turbulent fluctuations.

Significance

1. Reactor Relevance: Detached conditions represent the preferred operational state for future reactors, where turbulence effects on reaction rates are more pronounced
2. Systematic Errors: Transport codes (such as SOLPS-ITER, EMC3) employ mean-field approximations that may introduce systematic errors in particle and energy balance
3. Knowledge Gap: Comprehensive global simulations in realistic divertor X-point geometry are lacking

Limitations of Existing Methods

• Transport codes adopt mean-field approximations for computational efficiency, neglecting turbulent fluctuation effects
• Previous studies have primarily focused on attached conditions, with limited investigation of detached conditions
• Comprehensive analysis in complex geometries is absent

Core Contributions

1. First systematic analysis of turbulence effects on atomic reaction rates under detached divertor conditions
2. Discovery that density-temperature correlation is the key factor determining the direction of reaction rate deviations
3. Quantification of differences between turbulent and mean-field reaction rates (up to 2-50 times)
4. Provision of realistic geometry simulation results based on ASDEX Upgrade experiments
5. Revelation of potential systematic errors in conventional mean-field methods under detached conditions

Methodology Details

Task Definition

Comparison of atomic reaction rates calculated from turbulent fluctuation inputs and mean-field inputs, quantifying differences and analyzing influencing factors.

Inputs: Plasma density n, electron temperature Te, neutral particle density nn Outputs: Ionization source siz, recombination sink src, impurity radiation prad Constraints: Quasi-neutrality condition, relevant reaction rate coefficients

Model Architecture

GRILLIX Code Framework

• Plasma Model: Global electromagnetic full-f drift-fluid model
• Neutral Particle Model: Three-moment fluid model
• Geometry Processing: Complex stellarator and divertor tokamak geometries

Reaction Rate Model

1. Ionization Source:

   siz = nn·n·⟨σv⟩iz
   

2. Recombination Sink:

   src = -n²·⟨σv⟩rc
   

3. Impurity Radiation:

   prad = n²·cimp·LZ(Te)
   

where reaction rate coefficients employ polynomial fits from OPEN-ADAS data:

⟨σv⟩ = exp(∑∑ αij(ln n)^m(ln Te)^n)

Technical Innovations

Post-processing Analysis Method

• Turbulent Average Rate: ⟨f⟩ = ⟨f(n, Te, nn)⟩t,φ (compute then average)
• Mean-field Rate: f⟨◦⟩ = f(⟨n⟩t,φ, ⟨Te⟩t,φ, ⟨nn⟩t,φ) (average then compute)

Correlation Analysis Method

1. Realistic Input: Using original turbulent data
2. Decorrelated Input: Random shuffling in time and toroidal angle
3. Ordered Input: Monotonically increasing arrangement producing positive correlation

Quantification Metrics

• Maximum Norm: max(⟨f⟩)/max(f⟨◦⟩)
• Volume Integral Norm: ∫⟨f⟩dV/∫f⟨◦⟩dV

Experimental Setup

Simulation Cases

1. Detached Case: Detached L-mode plasma from ASDEX Upgrade discharge #40333 with large fluctuation amplitudes (up to 500%)
2. Attached Case: Attached L-mode plasma from ASDEX Upgrade discharge #38839 with low fluctuation amplitudes (<50%)

Control Volume

Toroidal continuous region near X-point:

• Poloidal angle range: -2.042 rad < θpol < -1.728 rad
• Magnetic flux surface range: 1.0 < ρpol < 0.98

Impurity Model

Nitrogen impurity with concentration cimp = 5%, employing coronal equilibrium approximation

Data Sampling

• Temporal averaging: 50 snapshots spanning 0.1 ms
• Toroidal averaging: 16 poloidal planes

Experimental Results

Main Results

Attached Conditions

• Turbulent-to-mean-field ratio for all reaction rates ≈ 1
• Relative differences approximately 1%, negligible impact

Detached Conditions

Key Findings:

• Turbulent ionization and radiation rates are approximately 2 times lower than mean-field values
• Recombination rates increase significantly in turbulence (locally up to 50 times)
• Density-temperature correlation coefficient: Rn,Te = -0.46 (negative correlation)

Correlation Impact Experiments

Decorrelated Input (Rn,Te = 0.0)

• Ionization rate ratio recovers to ≈ 1
• Turbulent effects essentially disappear upon correlation elimination

Positive Correlation Input (Rn,Te = +0.93)

• Ionization rate ratio increases to 3.37
• Consistent with turbulence enhancement effects observed in previous studies

Physical Mechanism Analysis

Two-dimensional histogram analysis in n-Te phase space reveals:

1. Attached Conditions: Narrow fluctuation range, limited rate variation
2. Detached Conditions: Wide fluctuation range, covering 10^20-10^26 m^-3s^-1 rate range
3. Negative Correlation Effect: High temperature corresponds to low density, low temperature to high density, resulting in reduced overall reaction rates

Related Work

Historical Research

• Marandet et al. (2011): Turbulence effects on neutral particle transport
• Guzman et al. (2015): Turbulent scrape-off layer impurity ionization balance
• Thrysøe et al. (2016): Blob effects on neutral particles
• Fan et al. (2019): Turbulence effects in TOKAM3X-EIRENE code

Contributions of This Work

• First comprehensive analysis in realistic X-point geometry
• Systematic investigation of density-temperature correlation effects
• Quantitative results under detached conditions

Conclusions and Discussion

Main Conclusions

1. Condition Dependence: Turbulence effects on reaction rates are negligible under attached conditions but significant under detached conditions
2. Correlation Determinism: Density-temperature correlation determines the direction and magnitude of rate deviations
3. Systematic Errors: Mean-field methods may introduce systematic errors up to 2 times under detached conditions
4. Reactor Relevance: Given that future reactors will operate in detached mode, this effect is of significant importance

Limitations

1. Post-processing Analysis: Self-consistent turbulence-transport coupling simulations were not performed
2. Geometric Constraints: Analysis concentrated on specific control volumes
3. Model Assumptions: Simplified impurity model and quasi-neutrality assumptions employed

Future Directions

1. Development of coupled iterative schemes between turbulence and transport codes
2. Extension to additional impurity species and reaction types
3. Validation of results in different tokamak geometries
4. Development of corrected transport models accounting for turbulence effects

In-Depth Evaluation

Strengths

1. High Innovation: First systematic investigation of turbulence effects on atomic reaction rates under realistic detached conditions
2. Rigorous Methodology: Physical mechanisms revealed through controlled correlation analysis with scientific approach
3. High Practical Value: Provides important reference for improving transport codes
4. Reliable Results: Based on mature GRILLIX code and ASDEX Upgrade experimental data

Weaknesses

1. Limited Analysis Scope: Primarily concentrated on specific control volumes; global effects require further investigation
2. Model Simplifications: Simplifications in impurity model and boundary conditions may affect result generalizability
3. Insufficient Verification: Lack of direct comparison with experimental measurements

Impact

1. Academic Value: Provides important insights for plasma physics and fusion research
2. Engineering Significance: Offers guidance for ITER and future reactor divertor design
3. Methodological Contribution: Proposed analysis framework applicable to other turbulence-reaction coupling problems

Applicable Scenarios

• Tokamak edge plasma modeling
• Divertor physics research
• Fusion reactor design optimization
• Turbulent transport code development

References

This paper cites 27 important references covering key fields including plasma boundary physics, numerical simulation codes, and atomic databases, providing a solid theoretical foundation for the research.