Nonlinear Two-Fluid Simulations of Type-I Edge Localized Mode Crashes in MAST Upgrade

The Type-I edge localized mode (ELM) crash dynamics is investigated for the MAST Upgrade (MAST-U) spherical tokamak plasma through nonlinear two-fluid BOUT++ simulations. Contrary to conventional peeling-ballooning paradigms, in this work the ELM crashes in MAST-U are found to be triggered by resistive ballooning modes with growth rates scaling as γ∝η1/3, accompanied by formation of n=20–25 (n is the toroidal mode number) filamentary structures during nonlinear evolution. Parametric scans of the plasma resistivity (η) reveal an inverse proportionality between the ELM size and resistivity, aligning with prior MAST experimental observations but contrary to the classical resistive ballooning theory. Linear two-fluid simulations corroborate this anomalous scaling, while analysis of the numerical data identifies a critical stabilizing mechanism: electromagnetic-to-thermal energy exchange reduces the ELM size at high resistivity. Furthermore, nonlinear simulations demonstrate that the bootstrap current—though inactive in linear instability drives—amplifies pedestal collapse in the MAST-U scenario.