Gyrokinetic turbulent transport simulations on steady burning condition in D-T-He plasmas

Ion temperature gradient(ITG) and trapped electron modes(TEM) driven turbulent transport in an ITER-like plasma is investigated by means of multi-species gyrokinetic Vlasov simulations with D, T, He, and real-mass kinetic electrons including their inter-species collisions. Beyond the conventional zero-dimensional power balance analysis presuming the global energy and particle confinement times, gyrokinetic-simulation-based evaluation of a steady burning condition with He-ash exhaust and D-T fuel inward pinch is demonstrated. It is clarified that a significant imbalance appears in the turbulent particle flux for the fuel ions of D and T, depending on the D-T density ratio and the He-ash accumulation. Then several profile regimes to satisfy Reiter’s steady burning condition are, for the first time, identified by the gyrokinetic simulation. Also, the impacts of zonal flows and nonthermal He-ash on the optimal profile regimes are examined.

💡 Research Summary

This paper presents a first‑principles, multi‑species gyrokinetic Vlasov study of turbulent transport in an ITER‑like deuterium‑tritium‑helium (D‑T‑He) burning plasma. The authors employ the electromagnetic gyrokinetic code GKV, solving the full five‑dimensional gyrokinetic equation for deuterium, tritium, thermalized helium ash, and real‑mass kinetic electrons, including inter‑species collisions. Linear analysis shows that the ITG‑TEM hybrid instability persists in the presence of up to 10 % helium ash, with only minor changes in growth rates and frequencies when normalized by the effective charge‑density fraction.

Non‑linear simulations, performed with a (169 × 49) mode Fourier grid and (64 × 64 × 24) phase‑space resolution, reveal that energy fluxes for D, T, and He approach gyro‑Bohm scaling (Q_D ≈ Q_T ≈ 0.5 Q_gB, Q_He < 0.5 Q_gB, Q_e ≈ 1.1 Q_gB). Crucially, particle fluxes exhibit a pronounced species‑dependent asymmetry: Γ_D and Γ_T have opposite signs in different poloidal locations, and the net helium‑ash flux is outward (Γ_He > 0) while the fuel ions experience an inward pinch (Γ_D, Γ_T < 0). This asymmetry is absent in single‑effective‑mass (A_eff) models, underscoring the importance of treating each ion species explicitly.

Parameter scans varying the D‑T ratio (n_T/(n_T + n_D)) and helium‑ash concentration (0, 5, 10 %) demonstrate that while the ion energy transport remains roughly symmetric around a 50‑50 D‑T mix, the particle transport becomes increasingly unbalanced as helium ash accumulates. The imbalance grows with deviation from the 50‑50 mix, indicating that the D‑T ratio can be used as a control knob for fuel‑pinch strength.

The authors then reinterpret Reiter’s steady‑burning condition (τ_He < α τ_E, with α≈7–15) in terms of turbulent fluxes. By expressing the helium‑ash confinement time τ_He ≈ a²/D_eff,He and the energy confinement time τ_E ≈ a²/χ_eff,i, the condition becomes η_i T_i (n_i/n_He) Γ_He > Q_i/α together with Γ_He > 0 and Γ_D, Γ_T < 0. Gyrokinetic results identify a family of profile regimes satisfying these inequalities: D‑T ratios between 0.4 and 0.6, helium‑ash fractions ≤5 %, and temperature‑gradient‑to‑density‑gradient ratios η_i ≈ 0.8–1.0. In these regimes the effective α lies in the range 10–12, comfortably within Reiter’s requirement.

Additional simulations explore the role of zonal flows and non‑thermal helium ash (high‑energy alphas). Strong zonal flows suppress the turbulence amplitude, reducing both energy and particle fluxes and weakening the inward fuel pinch. Non‑thermal helium ash flattens the electron temperature gradient, damping the trapped‑electron mode and further lowering turbulent transport. These findings suggest that active control of zonal flow activity and helium‑ash energy distribution can be leveraged to maintain the steady‑burning condition.

Overall, the study demonstrates that multi‑species gyrokinetic turbulence simulations can quantitatively capture the delicate balance between fuel‑pinch and helium‑ash exhaust that is invisible to conventional zero‑dimensional power‑balance models. The identified optimal profile windows provide concrete guidance for fueling strategies, helium‑ash removal systems, and turbulence‑control techniques in ITER and future burning‑plasma devices.