Kramer, G. J.; Bortolon, A.; Ferraro, N. M.; Spong, D. A.; Crocker, N. A.; Darrow, D. S.; Fredrickson, E. D.; Kubota, S.; Park, J.-K.; Podesta, M.; Heidbrink, W. W.
In this paper we present initial simulations of pedestal control by Lithium Granule Injection (LGI) in NSTX. A model for small granule ablation has been implemented in the M3D-C1 code [1], allowing the simulation of realistic Lithium granule injections. 2D simulations in NSTX L-mode and H-mode plasmas are done and the effect of granule size, injection angle and velocity on the pedestal gradient increase are studied. For H-mode cases, the amplitude of the local pressure perturbation caused by the granules is highly dependent on the solid granule size. In our simulations, reducing the granule injection velocity allows one to inject more particles at the pedestal top.
We discuss a novel diagnostic allowing direct measurements of the local electric field in the edge region in NSTX/NSTX-U. This laser based diagnostic's principle consists of depleting the naturally populated $n=3$ level to a Rydberg state --sensitive to electric fields-- that will result in a suppression of part of the $D_{\alpha}$ emission. We refer to this approach as Laser-Induced Rydberg Spectroscopy (LIRyS). It is shown that the local electric field can be measured through the Stark induced resonances observed as dips in the $D_\alpha$ emission. Using forward-modeling of simulated absorption spectra, we show precisions reaching \SI{\pm 2}{\kilo\volt\per\meter} in regions with a local electric field of \SI{15}{\kilo\volt\per\meter}.
Z. R. Wang; A. H. Glasser; D. Brennan; Y. Q. Liu; J-K. Park
Abstract:
The method of solving linear resistive plasma response, based on the asymptotic matching approach, is developed for full toroidal tokamaks by upgrading the Resistive DCON code [A.H. Glasser, Z.R. Wang and J.-K. Park, Physics of Plasmas, \textbf{23}, 112506 (2016)]. The derived matching matrix, asymptotically matching the outer and inner regions, indicates that the applied three dimension (3-D) magnetic perturbations contribute additional small solutions at each resonant surface due to the toroidal coupling of poloidal modes. In contrast, the resonant harmonic only affects the corresponding resonant surface in the cylindrical plasma. Since the solution of ideal outer region is critical to the asymptotic matching and is challenging to be solved in the toroidal geometry due to the singular power series solution at the resonant surfaces, systematic verification of the outer region $\Delta^\prime$ matrix is made by reproducing the well known analytical $\Delta^{\prime}$ result in [H.P. Furth, P.H. Rutherford and H. Selberg, The Physics of Fluids, \textbf{16}, 1054-1063 (1073)] as well as by making a quantitative benchmark with the PEST3 code [A. Pletzer and R.L. Dewar, J. Plasma Physics, \textbf{45}, 427-451 (1991)]. Finally, the reconstructed numerical solution of resistive plasma response from the toroidal matching matrix is presented. Comparing with the ideal plasma response, the global structure of the response can be affected by the small finite island at the resonant surfaces.
The Electromagnetic Particle Injector (EPI) concept is advanced through the simulation of ablatant deposition into ITER H-mode discharges with calculations showing penetration past the H-mode pedestal for a range of injection velocities and granule sizes concurrent with the requirements of disruption mitigation. As discharge stored energy increases in future fusion devices such as ITER, control and handling of disruption events becomes a critical issue. An unmitigated disruption could lead to failure of the plasma facing components resulting in financially and politically costly repairs. Methods to facilitate the quench of an unstable high current discharge are required. With the onset warning time for some ITER disruption events estimated to be less than 10 ms, a disruption mitigation system needs to be considered which operates at injection speeds greater than gaseous sound speeds. Such an actuator could then serve as a means to augment presently planned pneumatic injection systems. The EPI uses a rail gun concept whereby a radiative payload is delivered into the discharge by means of the JxB forces generated by an external current pulse, allowing for injection velocities in excess of 1 km/s. The present status of the EPI project is outlined, including the addition of boost magnetic coils. These coils augment the self-generated rail gun magnetic field and thus provide a more efficient acceleration of the payload. The coils and the holder designed to constrain them have been modelled with the ANSYS code to ensure structural integrity through the range of operational coil cu
These GROMACS trajectories show the existence of a critical point in deeply supercooled WAIL water. Also included is the code necessary to reproduce the figures in the corresponding paper from these trajectories. From this data the critical temperature, pressure, and density of the model can be found, and critical fluctuations in the deeply supercooled liquid can be directly observed (in a computer-simulation sense).
Movies of relativistic reconnection and particle acceleration in relativistic reconnection accompanying the article "Relativistic Reconnection: an Efficient Source of Nonthermal Particles" by Lorenzo Sironi and Anatoly Spitkovsky.
Magnetic reconnection in partially ionized plasmas is a ubiquitous and important phenomena in both laboratory and astrophysical systems. Here, simulations of partially ionized magnetic reconnection with well-matched initial conditions are performed using both multi-fluid and fully-kinetic approaches. Despite similar initial conditions, the time-dependent evolution differs between the two models. In multi-fluid models, the reconnection rate locally obeys either a decoupled Sweet-Parker scaling, where neutrals are unimportant, or a fully coupled Sweet-Parker scaling, where neutrals and ions are strongly coupled, depending on the resistivity. In contrast, kinetic models show a faster reconnection rate that is proportional to the fully-coupled, bulk Alfv\'en speed, $v_A^\star$. These differences are interpreted as the result of operating in different collisional regimes. Multi-fluid simulations are found to maintain $\nu_{ni}L/v_A^\star \gtrsim 1$, where $\nu_{ni}$ is the neutral-ion collision frequency and $L$ is the time-dependent current sheet half-length. This strongly couples neutrals to the reconnection outflow, while kinetic simulations evolve to allow $\nu_{ni}L/v_A^\star < 1$, decoupling neutrals from the reconnection outflow. Differences in the way reconnection is triggered may explain these discrepancies.
Lunsford, R.; Bortolon, A.; Roquemore, A.L.; Mansfield, D.K.; Jaworski, M.A.; Kaita, R.; Maingi, R.; Nagy, A.
Abstract:
By employing a neutral gas shielding (NGS) model to characterize impurity granule
injection the pedestal atomic deposition for three different species of granule:
lithium, boron, and carbon are determined. Utilizing the duration of ablation
events recorded on experiments performed at DIII-D to calibrate the NGS model we
are able to quantify the ablation rate and mass deposition location with respect
to the plasma density profile. The species specific granule shielding constant
is then used to model granule ablation within NSTX-U discharges. Simulations of
300, 500 and 700 micron diameter granules injected at 50 m/sec are presented for
NSTX-U L-mode type plasmas as well as H-mode discharges with low natural ELM
frequencies. Additionally, ablation calculations of 500 micron granules of each
species are presented at velocities ranging from 50 � 150 m/sec. In H-mode type
discharges these simulations show that the majority of the injected granule is
ablated within or just past the steep gradient region of the discharge. At this
radial position, the perturbation to the background plasma generated by the ablating
granule can lead to conditions advantageous for the rapid triggering of an ELM crash
event.