Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userCoupling JOREK and STARWALL for Non-linear Resistive-wall Simulations
127
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The implementation of a resistive-wall extension to the non-linear MHD-code JOREK via a coupling to the vacuum-field code STARWALL is presented along with first applications and benchmark results. Also, non-linear saturation in the presence of a resistive wall is demonstrated. After completion of the ongoing verification process, this code extension will allow to perform non-linear simulations of MHD instabilities in the presence of three-dimensional resistive walls with holes for limited and X-point plasmas.
rate research
Read More
A collisional-radiative non-equilibrium impurity treatment for JOREK 3D nonlinear magneto-hydrodynamic (MHD) simulations has been developed. The impurities are represented by super-particles flowing along the fluid velocity field lines, while ionizing and recombining independently according to ADAS data and local fluid density and temperature. The non-equilibrium impurity contributions are then projected back to the fluid field for self-consistent time evolution. A 2D test case is used to compare the new non-equilibrium impurity model against previous Coronal Equilibrium (CE) impurity treatment, as well as to compare the non-equilibrium impurity behavior between the single and the two temperature model. Further, we conduct benchmark with previously published coronal non-equilibrium results by other 3D nonlinear MHD codes such as M3D-C1 and NIMROD. The new non-equilibrium treatment is shown to successfully capture the early phase cooling by weakly ionized impurities which the CE model missed. The benchmarks with M3D-C1 and NIMROD show general agreement in both the integrated quantities and the 2D profile evolution, despite the difference in the atomic model used. The above comparison and benchmark cases demonstrate the capability of the non-equilibrium impurity model for JOREK, paving the way for more sophisticated 3D non-linear Massive Material Injection (MMI) simulations which have important applications in disruption mitigation studies.
JOREK is a massively parallel fully implicit non-linear extended MHD code for realistic tokamak X-point plasmas. It has become a widely used versatile code for studying large-scale plasma instabilities and their control developed in an international community. This article gives a comprehensive overview of the physics models implemented, numerical methods applied for solving the equations and physics studies performed with the code. A dedicated section highlights some of the verification work done for the code. A hierarchy of different physics models is available including a free boundary and resistive wall extension and hybrid kinetic-fluid models. The code allows for flux-surface aligned iso-parametric finite element grids in single and double X-point plasmas which can be extended to the true physical walls and uses a robust fully implicit time stepping. Particular focus is laid on plasma edge and scrape-off layer (SOL) physics as well as disruption related phenomena. Among the key results obtained with JOREK regarding plasma edge and SOL, are deep insights into the dynamics of edge localized modes (ELMs), ELM cycles, and ELM control by resonant magnetic perturbations, pellet injection, as well as by vertical magnetic kicks. Also ELM free regimes, detachment physics, the generation and transport of impurities during an ELM, and electrostatic turbulence in the pedestal region are investigated. Regarding disruptions, the focus is on the dynamics of the thermal quench and current quench triggered by massive gas injection (MGI) and shattered pellet injection (SPI), runaway electron (RE) dynamics as well as the RE interaction with MHD modes, and vertical displacement events (VDEs). Also the seeding and suppression of tearing modes (TMs), the dynamics of naturally occurring thermal quenches triggered by locked modes, and radiative collapses are being studied.
Injecting frozen deuterium pellets into an ELMy H-mode plasma is a well established scheme for triggering edge localized modes (ELMs) before they naturally occur. Based on an ASDEX Upgrade H-mode plasma, this article presents a comparison of extended MHD simulations of spontaneous type-I ELMs and pellet-triggered ELMs allowing to study their non-linear dynamics in detail. In particular, pellet-triggered ELMs are simulated by injecting deuterium pellets into different time points during the pedestal build-up described in [A. Cathey et al. Nuclear Fusion 60, 124007 (2020)]. Realistic ExB and diamagnetic background plasma flows as well as the time dependent bootstrap current evolution are included during the build-up to capture the balance between stabilising and destabilising terms for the edge instabilities accurately. Dependencies on the pellet size and injection times are studied. The spatio-temporal structures of the modes and the resulting divertor heat fluxes are compared in detail between spontaneous and triggered ELMs. We observe that the premature excitation of ELMs by means of pellet injection is caused by a helical perturbation described by a toroidal mode number of n = 1. In accordance with experimental observations, the pellet-triggered ELMs show reduced thermal energy losses and narrower divertor wetted area with respect to spontaneous ELMs. The peak divertor energy fluency is seen to decrease when ELMs are triggered by pellets injected earlier during the pedestal build-up.
The dynamics of large scale plasma instabilities can strongly be influenced by the mutual interaction with currents flowing in conducting vessel structures. Especially eddy currents caused by time-varying magnetic perturbations and halo currents flowing directly from the plasma into the walls are important. The relevance of a resistive wall model is directly evident for Resistive Wall Modes (RWMs) or Vertical Displacement Events (VDEs). However, also the linear and non-linear properties of most other large-scale instabilities may be influenced significantly by the interaction with currents in conducting structures near the plasma. The understanding of halo currents arising during disruptions and VDEs, which are a serious concern for ITER as they may lead to strong asymmetric forces on vessel structures, could also benefit strongly from these non-linear modeling capabilities. Modeling the plasma dynamics and its interaction with wall currents requires solving the magneto-hydrodynamic (MHD) equations in realistic toroidal X-point geometry consistently coupled with a model for the vacuum region and the resistive conducting structures. With this in mind, the non-linear finite element MHD code JOREK has been coupled with the resistive wall code STARWALL, which allows to include the effects of eddy currents in 3D conducting structures in non-linear MHD simulations. This article summarizes the capabilities of the coupled JOREK-STARWALL system and presents benchmark results as well as first applications to non-linear simulations of RWMs, VDEs, disruptions triggered by massive gas injection, and Quiescent H-Mode. As an outlook, the perspectives for extending the model to halo currents are described.
The low-frequency rotating plasma instability (spoke) in the ISCT200 thruster operating in the wall-less configuration was simulated with a 3 dimensional PIC MCC code. In the simulations an m = 1 spoke rotating with a velocity of 6.5 km/s in the ExB direction was observed. The rotating electron density structure in the spoke is accompanied by a strongly depleted region of the neutral gas, which clearly shows that the spoke instability is of an ionization nature, similar to the axial breathing mode oscillations. In the simulation the electron cross-field transport through the spoke core was caused by diffusion in the high-frequency (4-10 MHz), short-scale (3 mm) electric field oscillations. These short-scale oscillations play a crucial role in the thruster discharge as over 70% of the electron current to the anode originates from the spoke core. The rest of the current originates from the spoke front where the electron cross-field transport toward the anode is due to the ExB drift in the spoke macroscopic azimuthal electric field.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
