Time-dependent runaway electron simulations: Ampere–Faraday equations implemented in CQL3D

Nuclear Fushion

article

R.W. Harvey¹  Yu.V. Petrov¹  Charlson C. Kim²  C.B. Forest³  L.L. Lao⁴  P.B. Parks⁴ 

[1] CompX;SLS2 Consulting;University of Wisconsin;General Atomics

关键词: runaway electrons;    Fokker–Planck;    simulation;    Ampere–Faraday equations;   

DOI  :  10.1088/1741-4326/ab38cb

来源: Institute of Physics Publishing Ltd.

【 摘 要 】

The runaway electron (RE) distributions driven by a large toroidal electric field induced by the drop in the temperature profile due to disruption or pellets are comprehensively simulated by the 3D Fokker–Planck (FP) solver CQL3D (Harvey and McCoy 1992Proc. of IAEA TCM ), recently coupled to the Ampere–Faraday (AF) equations. The evolution of the toroidal current in a plasma occurs on a resistive time scale,τ _res  =  4 πa ²/( c ² η ), which is typically of the order of seconds in present tokamaks. Here,aandηare respectively plasma radius or radial extent of a current density perturbation, and Ohmic resistivity. From the Faraday EM equation, the toroidal electric field is proportional to the time derivative of the poloidal magnetic field, which, from the Ampere equation, is proportional to the toroidal current. Thus, the toroidal electric field rapidly increases due to an abrupt temperature drop decrease in conductivity, to prevent change in the toroidal current faster thanτ _res. This is a example of Lenz's law. For example, in simulations with KPRAD (Whyteet al2003J. Nucl. Mater . 313–6 1239) of neon pellet injection into a DIII-D shot,T _e drops from 2 keV to 10 eV in 0.1 ms andZ _eff increases 1–4, giving that the electric field increases 3500 ×   to 0.8 V cm⁻¹. As described in Harveyet al(2000 PoP 7 4590), this places much of the tail electron distribution beyond the Dreicer runaway velocity, giving so-called 'hot-tail runaways' which for a time are the dominant source of runaways, more so than the knockon source. In this prior calculation, performed for a single flux surface, the toroidal current density is held constant, on the basis thatτ _res is large. Most of the initial current can be converted to runaway current, which is then dangerous, particularly for ITER. A more comprehensive A–F model recently implemented in CQL3D, taking into account the time-development of the full-plasma-width toroidal electric field on time-scales of orderτ _res applies an iterative technique for the toroidal field previously developed for a different application (Kupferet al1996 PoP 3 3644), maintaining the implicit-in-time evolution of CQL3D. The degree of runaway current formation is reduced in AF augmented CQL3D, but the basic mechanism of 'hot-tail runaways' remains a dominant contribution to the REs at early times after theT _e drop in these simulations. On the other hand, a NIMROD (Sovinecet al2004J. Comput. Phys. 195 355 ) simulation of shattered-pellet shutdown of DIII-D plasma (Kim 2018 APS/DPP Meeting), gives a slower thermal quench; when the plasma profiles and electric field are coupled one-way to CQL3D, the 'hot-tail' REs are much less, and growth of RE is dominated by the knockon process.

【 授权许可】

Unknown