No Arabic abstract
Single carbon pellet disruption mitigation simulations using M3D-C1 were conducted in an NSTX-U-like plasma to support the electromagnetic pellet injection concept (EPI). A carbon ablation model has been implemented in M3D-C1 and tested with available data. 2D simulations were conducted in order to estimate the amount of carbon needed to quench the plasma, finding that the content in a $1,$mm radius vitreous carbon pellet (~ 3.2x10E20 atoms) would be enough if it is entirely ablated. 3D simulations were performed, scanning over pellet velocity and parallel thermal conductivity, as well as different injection directions and pellet concepts (solid pellets and shell pellets). The sensitivity of the thermal quench and other related quantities to these parameters has been evaluated. A 1 mm radius solid pellet only partially ablates at velocities of 300 m/s or higher, thus being unable to fully quench the plasma. To further enhance the ablation, approximations to an array of pellets and the shell pellet concept were also explored. 3D field line stochastization plays an important role in both quenching the center of the plasma and in heat flux losses, thus lowering the amount of carbon needed to mitigate the plasma when compared to the 2D case. This study constitutes an important step forward in `predict-first simulations for disruption mitigation in NSTX-U and other devices, such as ITER.
A dynamic mitigation mechanism of the two-stream instability is discussed based on a phase control for plasma and fluid instabilities. The basic idea for the dynamic mitigation mechanism by the phase control was proposed in the paper [Phys. Plasmas 19, 024503(2012)]. The mitigation method is applied to the two-stream instability in this paper. In general, instabilities appear from the perturbations, and normally the perturbation phase is unknown. Therefore, the instability growth rate is discussed in fluids and plasmas. However, if the perturbation phase is known, the instability growth can be controlled by a superimposition of perturbations imposed actively. For instance, a perturbed driver induces a perturbation to fluids or plasmas; if the perturbation induced by the perturbed driver is oscillated actively by the driver oscillation, the perturbation phase is known and the perturbation amplitude can be controlled, like a feedforward control. The application result shown in this paper demonstrates that the dynamic mitigation mechanism works well to smooth the non-uniformities and mitigate the instabilities in plasmas.
Current models predict the hose instability to crucially limit the applicability of plasma-wakefield accelerators. By developing an analytical model which incorporates the evolution of the hose instability over long propagation distances, this work demonstrates that the inherent drive-beam energy loss, along with an initial beam energy spread detune the betatron oscillations of beam electrons, and thereby mitigate the instability. It is also shown that tapered plasma profiles can strongly reduce initial hosing seeds. Hence, we demonstrate that the propagation of a drive beam can be stabilized over long propagation distances, paving the way for the acceleration of high-quality electron beams in plasma-wakefield accelerators. We find excellent agreement between our models and particle-in-cell simulations.
An internal disruption is simulated in a large tokamak plasma with monotonic safety factor profile close to unity. The domain and the time scale of the event are set to match observations. The simulation follows passive alpha particles with energies 35 keV-3.5 MeV, whose initial density peak is localized in the disrupting domain. While the 35 keV profile flattens, a synergy of multiple physical factors allows the 3.5 MeV profile to remain peaked, motivating the use of moderate internal disruptions in a fusion reactor to expel helium ash while preserving good confinement of fast alphas.
Numerical simulation of plasma turbulence in the Large Plasma Device (LAPD) [Gekelman et al, Rev. Sci. Inst., 62, 2875, 1991] is presented. The model, implemented in the BOUndary Turbulence (BOUT) code [M. Umansky et al, Contrib. Plasma Phys. 180, 887 (2009)], includes 3-D collisional fluid equations for plasma density, electron parallel momentum, and current continuity, and also includes the effects of ion-neutral collisions. In nonlinear simulations using measured LAPD density profiles but assuming constant temperature profile for simplicity, self-consistent evolution of instabilities and nonlinearly-generated zonal flows results in a saturated turbulent state. Comparisons of these simulations with measurements in LAPD plasmas reveal good qualitative and reasonable quantitative agreement, in particular in frequency spectrum, spatial correlation and amplitude probability distribution function of density fluctuations. For comparison with LAPD measurements, the plasma density profile in simulations is maintained either by direct azimuthal averaging on each time step, or by adding particle source/sink function. The inferred source/sink values are consistent with the estimated ionization source and parallel losses in LAPD. These simulations lay the groundwork for more a comprehensive effort to test fluid turbulence simulation against LAPD data.
A fusion boundary-plasma domain is defined by axisymmetric magnetic surfaces where the geometry is often complicated by the presence of one or more X-points; and modeling boundary plasmas usually relies on computational grids that account for the magnetic field geometry. The new grid generator INGRID (Interactive Grid Generator) presented here is a Python-based code for calculating grids for fusion boundary plasma modeling, for a variety of configurations with one or two X-points in the domain. Based on a given geometry of the magnetic field, INGRID first calculates a skeleton grid which consists of a small number of quadrilateral patches; then it puts a subgrid on each of the patches, and joins them in a global grid. This domain partitioning strategy makes possible a uniform treatment of various configurations with one or two X-points in the domain. This includes single-null, double-null, and other configurations with two X-points in the domain. The INGRID design allows generating grids either interactively, via a parameter-file driven GUI, or using a non-interactive script-controlled workflow. Results of testing demonstrate that INGRID is a flexible, robust, and user-friendly grid-generation tool for fusion boundary-plasma modeling.