Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userSimulations of COMPASS Vertical Displacement Events with a self-consistent model for halo currents including neutrals and sheath boundary conditions
94
0
0.0
(
0
)
Added by
Francisco Javier Artola F.J.
Publication date
2021
fields
Physics
and research's language is
English
Ask ChatGPT about the research
No Arabic abstract
The understanding of the halo current properties during disruptions is key to design and operate large scale tokamaks in view of the large thermal and electromagnetic loads that they entail. For the first time, we present a fully self-consistent model for halo current simulations including neutral particles and sheath boundary conditions. The model is used to simulate Vertical Displacement Events (VDEs) occurring in the COMPASS tokamak. Recent COMPASS experiments have shown that the parallel halo current density at the plasma-wall interface is limited by the ion saturation current during VDE-induced disruptions. We show that usual MHD boundary conditions can lead to the violation of this physical limit and we implement this current density limitation through a boundary condition for the electrostatic potential. Sheath boundary conditions for the density, the heat flux, the parallel velocity and a realistic parameter choice (e.g. Spitzer $eta$ and Spitzer-Harm $chi_parallel$ values) extend present VDE simulations beyond the state of the art. Experimental measurements of the current density, temperature and heat flux profiles at the COMPASS divertor are compared with the results obtained from axisymmetric simulations. Since the ion saturation current density ($J_{sat}$) is shown to be essential to determine the halo current profile, parametric scans are performed to study its dependence on different quantities such as the plasma resistivity and the particle and heat diffusion coefficients. In this respect, the plasma resistivity in the halo region broadens significantly the $J_{sat}$ profile, increasing the halo width at a similar total halo current.
rate research
Read More
We describe a hierarchy of stochastic boundary conditions (SBCs) that can be used to systematically eliminate finite size effects in Monte Carlo simulations of Ising lattices. For an Ising model on a $100 times 100$ square lattice, we measured the specific heat, the magnetic susceptibility, and the spin-spin correlation using SBCs of the two lowest orders, to show that they compare favourably against periodic boundary conditions (PBC) simulations and analytical results. To demonstrate how versatile the SBCs are, we then simulated an Ising lattice with a magnetized boundary, and another with an open boundary, measuring the magnetization, magnetic susceptibility, and longitudinal and transverse spin-spin correlations as a function of distance from the boundary.
The Particle-in-Cell (PIC) method was used to study two different ion thruster concepts - Stationary Plasma Thrusters (SPT) and High Efficiency Multistage Plasma Thrusters (HEMP-T), in particular the plasma properties in the discharge chamber due to the different magnetic field configurations. Special attention was paid to the simulation of plasma particle fluxes on the thrusters channel surfaces. In both cases, PIC proved itself as a powerful tool, delivering important insight into the basic physics of the different thruster concepts. The simulations demonstrated that the new HEMP thruster concept allows for a high thermal efficiency due to both minimal energy dissipation and high acceleration efficiency. In the HEMP thruster the plasma contact to the wall is limited only to very small areas of the magnetic field cusps, which results in much smaller ion energy flux to the thruster channel surface as compared to SPT. The erosion yields for dielectric discharge channel walls of SPT and HEMP thrusters were calculated with the binary collision code SDTrimSP. For SPT, an erosion rate on the level of 1 mm of sputtered material per hour was observed. For HEMP, thruster simulations have shown that there is no erosion inside the dielectric discharge channel.
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.
Multifluid simulations of plasma sheaths are increasingly used to model a wide variety of problems in plasma physics ranging from global magnetospheric flows around celestial bodies to plasma-wall interactions in thrusters and fusion devices. For multifluid problems, accurate boundary conditions to model an absorbing wall that resolves a classical sheath remains an open research area. This work justifies the use of vacuum boundary conditions for absorbing walls to show comparable accuracy between a multifluid sheath and lower moments of a continuum-kinetic sheath.
This paper discusses the capture of an ion beam in a magnetized plasma of an Electron Cyclotron Resonance Ion Source based Charge Breeder, as modelled by numerical simulations. As a relevant step forward with respect to previous works, here the capture is modeled by considering a plasma structure determined in a self-consisent way. The plasmoid-halo structure of the ECR plasma - that is consisting of a dense core (the plasmoid) surrounded by a rarefied halo - is further confirmed by the self-consistent calculations, having also some fine structures affected by the electromagnetic field distribution and by the magnetostatic field profile. The capture of Rb1+ ions has been investigated in details, vs. various plasma parameters, and then compared to experimental results.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
