Abstract:The edge plasma turbulence and transport dynamics, as well as the divertor power loads during the thermal quench phase of tokamak disruptions are numerically investigated with BOUT++'s flux-driven, six-field electromagnetic turbulence model. Here a transient yet intense particle and energy sources are applied at the pedestal top to mimic the plasma power drive at the edge induced by a core thermal collapse, which flattens core temperature profile. Interesting features such as surging of divertor heat load (up to 50 times), and broadening of heat flux width (up to 4 times) on the outer divertor target plate, are observed in the simulation, in qualitative agreement with experimental observations. The dramatic changes of divertor heat load and width are due to the enhanced plasma turbulence activities inside the separatrix. Two cross-field transport mechanisms, namely the $E\times B$ turbulent convection and the stochastic parallel advection/conduction, are identified to play important roles in this process. Firstly, elevated edge pressure gradient drives instabilities and subsequent turbulence in the entire pedestal region. The enhanced turbulence not only transports particles and energy radially across the separatrix via $E\times B$ convection which causes the initial divertor heat load burst, but also induces an amplified magnetic fluctuation $\tilde{B}$. Once the magnetic fluctuation is large enough to break the magnetic flux surface, magnetic flutter effect provides an additional radial transport channel. In the late stage of our simulation, $|\tilde{B}_r/B_0|$ reaches to $10^{-4}$ level that completely breaks magnetic flux surfaces such that stochastic field-lines are directly connecting pedestal top plasma to the divertor target plates or first wall, further contributing to the divertor heat flux width broadening.

What problem does this paper attempt to address?

The problem that this paper attempts to solve is the changes in edge plasma turbulence dynamics and divertor heat load during the thermal quench phase in tokamak plasma discharge. Specifically, the author uses the six - field electromagnetic turbulence model under the BOUT++ framework to study, through numerical simulation, the turbulence and transport characteristics of tokamak edge plasma during thermal quench, and how these changes affect the heat load on the divertor. ### Main research questions: 1. **Dynamics of edge plasma turbulence and transport**: During the thermal quench process, how does the plasma energy release due to the rapid temperature drop in the core region cause an increase in the pressure gradient of the edge plasma, thereby exciting turbulence activities? 2. **Changes in divertor heat load**: During the thermal quench process, how does the heat load on the divertor increase significantly (up to 50 times), and how does the heat flux width expand (up to 4 times)? 3. **Role of turbulence mechanisms**: Identify and analyze the roles of two cross - field transport mechanisms - E×B turbulence convection and stochastic parallel convection/ conduction - in this process. ### Specific research contents: - **Turbulence enhancement**: Through numerical simulation, it is observed that during the thermal quench process, the turbulence activity of the edge plasma is significantly enhanced, especially the pressure - gradient - driven instability, which leads to a significant increase in the turbulence intensity and range in the entire pedestal region. - **Impact of magnetic perturbations**: Turbulence activities not only enhance radial transport but also cause the amplification of magnetic field perturbations, especially the increase in the radial magnetic field component $\tilde{B}_r$. When $\left|\frac{\tilde{B}_r}{B_0}\right|$ reaches the level of $10^{-4}$, the magnetic surface structure is destroyed, causing more plasma to be directly connected to the divertor or the first wall, further increasing the heat load on the divertor. - **Heat load response**: It is found that the heat load response on the divertor lags behind the turbulence activity in the mid - plane, but overall there is a good correspondence between the two. Electron parallel conduction transfers energy from the upstream to the downstream in a short time, resulting in a rapid rise in the divertor heat load. ### Mathematical formulas: - **Normalized radial magnetic field component**: \[ \left|\frac{\tilde{B}_r}{B_0}\right| \] - **Parallel gradient operator**: \[ \nabla_\parallel f=\left(\hat{b}_0+\tilde{b}\right)\cdot\nabla f = \hat{b}_0\cdot\nabla f-\frac{\hat{b}_0}{B}\times\nabla A_\parallel\cdot\nabla f \] - **Normalized pressure perturbation**: \[ \frac{\tilde{p}}{p_0} \] Through these studies, the author aims to provide important references for the design and operation of future tokamak reactors (such as ITER), especially in dealing with thermal quench events, how to effectively manage and reduce the thermal impact on the divertor and the first wall.

Electromagnetic turbulence simulation of tokamak edge plasma dynamics and divertor heat load during thermal quench

Advanced Simulation of Electron Heat Transport in Fusion Plasmas

Gyrokinetic Simulation of Dissipative Trapped Electron Mode in Tokamak Edge

The effect of divertor particle sources on scrape-off-layer turbulence

Studies of edge poloidal rotation and turbulence momentum transport during divertor detachment on HL-2A tokamak

Tokamak edge-SOL turbulence in H-mode conditions simulated with a global, electromagnetic, transcollisional drift-fluid model

Role of turbulent separatrix tangle in the improvement of the integrated pedestal and heat exhaust issue for stationary-operation tokamak fusion reactors

Broadening of the Divertor Heat Flux Profile in High Confinement Tokamak Fusion Plasmas with Edge Pedestals Limited by Turbulence in DIII-D

Fluid turbulence simulations of divertor heat load for ITER hybrid scenario using BOUT++

Simulations of tokamak edge plasma turbulent fluctuations based on a minimal 3D model

Prediction of Divertor Heat Flux Width for ITER Using BOUT++ Transport and Turbulence Module

Corrigendum: Prediction of Divertor Heat Flux Width for ITER Using BOUT++ Transport and Turbulence Module (2019 Nucl. Fusion 59 046014)

How turbulence spreading improves power handling in quiescent high confinement fusion plasmas

Effects of runaway electrons on plasma facing components

Turbulent features near the X point of a DTT-like tokamak plasma

Dependence of divertor turbulence on plasma density and current in TCV

Simulating the Impact of Baffling on Divertor Performance Using SOLPS-ITER

Divertor turbulent transport in the single null and snowflake in the TCV tokamak

Self-consistent multi-component simulation of plasma turbulence and neutrals in detached conditions