Nonlinear ELM simulations based on a nonideal peeling–ballooning model using the BOUT++ code

Journal Article · Fri Sep 23 00:00:00 EDT 2011 · Nuclear Fusion

DOI:https://doi.org/10.1088/0029-5515/51/10/103040· OSTI ID:1313556

Xu, X. Q. ^[1]; Dudson, B. D. ^[2]; Snyder, P. B. ^[3]; Umansky, M. V. ^[1]; Wilson, H. R. ^[2]; Casper, T. ^[4]

---

1. Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
2. Univ. of York, York (United Kingdom)
3. General Atomics, San Diego, CA (United States)
4. ITER Organization, St. Paul lez Durance (France)

A minimum set of equations based on the peeling–ballooning (P–B) model with nonideal physics effects (diamagnetic drift, E × B drift, resistivity and anomalous electron viscosity) is found to simulate pedestal collapse when using the BOUT++ simulation code, developed in part from the original fluid edge code BOUT. Linear simulations of P–B modes find good agreement in growth rate and mode structure with ELITE calculations. The influence of the E × B drift, diamagnetic drift, resistivity, anomalous electron viscosity, ion viscosity and parallel thermal diffusivity on P–B modes is being studied; we find that (1) the diamagnetic drift and E × B drift stabilize the P–B mode in a manner consistent with theoretical expectations; (2) resistivity destabilizes the P–B mode, leading to resistive P–B mode; (3) anomalous electron and parallel ion viscosities destabilize the P–B mode, leading to a viscous P–B mode; (4) perpendicular ion viscosity and parallel thermal diffusivity stabilize the P–B mode. With addition of the anomalous electron viscosity under the assumption that the anomalous kinematic electron viscosity is comparable to the anomalous electron perpendicular thermal diffusivity, or the Prandtl number is close to unity, it is found from nonlinear simulations using a realistic high Lundquist number that the pedestal collapse is limited to the edge region and the ELM size is about 5–10% of the pedestal stored energy. Furthermore, this is consistent with many observations of large ELMs. The estimated island size is consistent with the size of fast pedestal pressure collapse. In the stable α-zones of ideal P–B modes, nonlinear simulations of viscous ballooning modes or current-diffusive ballooning mode (CDBM) for ITER H-mode scenarios are presented.

• Research Organization: Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
• Sponsoring Organization: USDOE
• Grant/Contract Number: AC52-07NA27344
• OSTI ID: 1313556
• Report Number(s): LLNL-JRNL-464898
• Journal Information: Nuclear Fusion, Vol. 51, Issue 10; ISSN 0029-5515
• Publisher: IOP Science
• Country of Publication: United States
• Language: English

Related Subjects

70 PLASMA PHYSICS AND FUSION