In order to study the distribution and evolution of the transient particle and heat fluxes during edge-localized mode (ELM) bursts on the Experimental Advanced Superconducting Tokamak (EAST), the BOUT++ six-field two-fluid model is used to simulate the pedestal collapse. The profiles from the EAST H-mode discharge #56129 are used as the initial conditions. Linear analysis shows that the resistive ballooning mode and drift-Alfven wave are two dominant instabilities for the equilibrium, and play important roles in driving ELMs. The evolution of the density profile and the growing process of the heat flux at divertor targets during the burst of ELMs are reproduced. The time evolution of the poloidal structures of Te is well simulated, and the dominant mode in each stage of the ELM crash process is found. The studies show that during the nonlinear phase, the dominant mode is 5, and it changes to 0 when the nonlinear phase goes to saturation after the ELM crash. The time evolution of the radial electron heat flux, ion heat flux, and particle density flux at the outer midplane (OMP) are obtained, and the corresponding transport coefficients Dr, χir, and χer reach maximum around 0.3 ∼ 0.5 m2 s−1 at ΨN = 0.9. The heat fluxes at outer target plates are several times larger than that at inner target plates, which is consistent with the experimental observations. The simulated profiles of ion saturation current density (js) at the lower outboard (LO) divertor target are compared to those of experiments by Langmuir probes. The profiles near the strike point are similar, and the peak values of js from simulation are very close to the measurements.
Six-field two-fluid simulations of peeling–ballooning modes using BOUT++