A new model of heating, current drive, torque and other effects of neutral beam injection on NSTX-U that uses neural networks has been developed. The model has been trained and tested on the results of the Monte Carlo code NUBEAM for the database of experimental discharges taken during the first operational campaign of NSTX-U. By projecting flux surface quantities onto empirically derived basis functions, the model is able to efficiently and accurately reproduce the behavior of both scalars, like the total neutron rate and shine through, and profiles, like beam current drive and heating. The model has been tested on the NSTX-U real-time computer, demonstrating a rapid execution time orders of magnitude faster than the Monte Carlo code that is well suited for the iterative calculations needed to interpret experimental results, optimization during scenario development activities, and real-time plasma control applications. Simulation results of a proposed design for a nonlinear observer that embeds the neural network calculations to estimate the poloidal flux profile evolution, as well as effective charge and fast ion diffusivity, are presented.
Active control of the toroidal current density profile is critical for the upgraded National Spherical Torus eXperiment device (NSTX-U) to maintain operation at the desired high-performance, MHD-stable, plasma regime. Initial efforts towards current density profile control have led to the development of a control-oriented, physics-based, plasma-response model, which combines the magnetic diffusion equation with empirical correlations for the kinetic profiles and the non-inductive current sources. The developed control-oriented model has been successfully tailored to the NSTX-U geometry and actuators. Moreover, a series of efforts have been made towards the design of model-based controllers, including a linear-quadratic-integral optimal control strategy that can regulate the current density profile around a prescribed target profile while rejecting disturbances. In this work, the tracking performance of the proposed current-profile optimal controller is tested in numerical simulations based on the physics-oriented code TRANSP. These high-fidelity closed-loop simulations, which are a critical step before experimental implementation and testing, are enabled by a flexible framework recently
developed to perform feedback control design and simulation in TRANSP.
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
Abstract: Tokamak plasma facing components have surface roughness that can cause microscopic spatial variations in erosion and deposition and hence influence material migration, erosion lifetime, dust and tritium accumulation, and plasma contamination. However high spatial resolution measurements of deposition on the scale of the surface roughness have been lacking to date. We will present elemental images of graphite samples from NSTX-U and DIII-D DiMES experiments performed with a Scanning Auger Microprobe at sub-micron resolution that show strong microscopic variations in deposition and correlate this with 3D topographical maps of surface irregularities. The NSTX-U samples were boronized and exposed to deuterium plasmas and the DiMES samples had localized Al and W films and were exposed to dedicated helium plasmas. Topographical maps of the samples were performed with a 3D confocal optical microscope and compared to the elemental deposition pattern. The results revealed localized deposition concentrated in areas shadowed from the ion flux, incident in a direction calculated (for the DiMES case) by taking account of the magnetic pre-sheath.
Cole M; Hager R; Moritaka T; Dominski J; Kleiber R; Ku S; Lazerson S; Riemann J; Chang C
XGC (X-point Gyrokinetic Code) is a whole-volume, total-f gyrokinetic particle-in-cell code developed for modelling tokamaks.In recent work, XGC has been extended to model more general 3D toroidal magnetic configurations, such as stellarators.These improvements have resulted in the XGC-S version.In this paper, XGC-S is benchmarked in the reduced delta-f limit for linear electrostatic ion temperature gradient-driven microinstabilities, which can underlie turbulent transport in stellarators.An initial benchmark of XGC-S in tokamak geometry shows good agreement with the XGC1, ORB5, and global GENE codes.A benchmark between XGC-S and the EUTERPE global gyrokinetic code for stellarators has also been performed, this time in geometry of the optimised stellarator Wendelstein 7-X.Good agreement has been found for the mode number spectrum, mode structure, and growth rate.
\f0\fs38 \cf2 The dynamics of the radial envelope of a weak coherent drift wave is approximately governed by a nonlinear Schr\'f6dinger equation, which emerges as a limit of the modified Hasegawa\'97Mima equation. The nonlinear Schr\'f6dinger equation has well-known soliton solutions, and its modulational instability can naturally generate solitary structures. In this paper, we demonstrate that this simple model can adequately describe the formation of solitary zonal structures in the modified Hasegawa\'97Mima equation, but only when the amplitude of the coherent drift wave is relatively small. At larger amplitudes, the modulational instability produces stationary zonal structures instead. Furthermore, we find that incoherent drift waves with beam-like spectra can also be modulationally unstable to the formation of solitary or stationary zonal structures, depending on the beam intensity. Notably, we show that these drift waves can be modeled as quantumlike particles (\'93driftons\'94) within a recently developed phase-space (Wigner\'97Moyal) formulation, which intuitively depicts the solitary zonal structures as quasi-monochromatic drifton condensates. Quantumlike effects, such as diffraction, are essential to these condensates; hence, the latter cannot be described by wave-kinetic models that are based on the ray approximation.\
In homogeneous drift-wave (DW) turbulence, zonal flows (ZFs) can be generated via a modulational instability (MI) that either saturates monotonically or leads to oscillations of the ZF energy at the nonlinear stage. This dynamics is often attributed as the predator-prey oscillations induced by ZF collisional damping; however, similar dynamics is also observed in collisionless ZFs, in which case a different mechanism must be involved. Here, we propose a semi-analytic theory that explains the transition between the oscillations and saturation of collisionless ZFs within the quasilinear Hasegawa-Mima model. By analyzing phase-space trajectories of DW quanta (driftons) within the geometrical-optics (GO) approximation, we argue that the parameter that controls this transition is N ~ \gamma_MI/\omega_DW, where \gamma_MI is the MI growth rate and \omega_DW is the linear DW frequency. We argue that at N << 1, ZFs oscillate due to the presence of so-called passing drifton trajectories, and we derive an approximate formula for the ZF amplitude as a function of time in this regime. We also show that at N >~ 1, the passing trajectories vanish and ZFs saturate monotonically, which can be attributed to phase mixing of higher-order sidebands. A modification of N that accounts for effects beyond the GO limit is also proposed. These analytic results are tested against both quasilinear and fully-nonlinear simulations. They also explain the earlier numerical results by Connaughton et al. [J. Fluid Mech. 654, 207 (2010)] and Gallagher et al. [Phys. Plasmas 19, 122115 (2012)] and offer a revised perspective on what the control parameter is that determines the transition from the oscillations to saturation of collisionless ZFs.
Schwartz, Jacob; Emdee, Eric; Goldston, Robert; Jaworski, Michael
The lithium vapor box divertor is a potential solution for power exhaust in toroidal confinement devices. The divertor plasma interacts with a localized, dense cloud of lithium vapor, leading to volumetric radiation, cooling, recombination, and detachment. To minimize contamination of the core plasma, lithium vapor is condensed on cool (300°C to 400°C) baffles upstream of the detachment point. Before implementing this in a toroidal plasma device with a slot divertor geometry, we consider an experiment with a scaled baffled-pipe geometry in the high-power linear plasma device Magnum-PSI. Three 15 cm-scale open cylinders joined by 6 cm diameter ‘nozzles’ are positioned on the plasma beam axis upstream of a target. The central box may be loaded with several tens of grams of lithium, which can be evaporated at 650°C to produce a vapor predicted, using a simple plasma-neutral interaction model, to be dense enough to cause volumetric detachment in the plasma. The power delivered to the target and box walls as measured by increases in their temperatures after a 10 s plasma pulse can be compared to determine the effectiveness of the vapor in detaching the plasma. Direct Simulation Monte Carlo simulations are performed to estimate the flow rates of lithium vapor between the boxes and to estimate the trapping of H2 delivered by the plasma in the boxes, which could inadvertently lead to detachment. Details of the geometry, simulations, and possible diagnostic techniques are presented.
F. M. Laggner, A. Diallo, B. P. LeBlanc, R. Rozenblat, G. Tchilinguirian, E.Kolemen, the NSTX-U team
A detailed description of a prototype setup for real-time (rt) Thomson scattering (TS) analysis is presented and implemented in the multi-point Thomson scattering (MPTS) diagnostic system at the National Spherical Torus Experiment Upgrade(NSTX-U). The data acquisition hardware was upgraded with rt capable electronics (rt-analog digital converters (ADCs) and a rt server) that allow for fast digitization of the laser pulse signal of eight radial MPTS channels. In addition, a new TS spectrum analysis software for a rapid calculation of electron temperature (Te) and electron density (ne) was developed. Testing of the rt hardware and data analysis soft-ware was successfully completed and benchmarked against the standard, post-shot evaluation. Timing tests were performed showing that the end-to-end processing time was reproducibly below 17 ms for the duration of at least 5 s, meeting a 60 Hz deadline by the laser pulse repetition rate over the length of a NSTX-U discharge. The presented rt framework is designed to be scalable in system size, i.e. incorporation of additional radial channels by solely adding additional rt capable hardware. Furthermore, it is scalable in its operation duration and was continuously run for up to 30 min, making it an attractive solution for machines with long discharge duration such as advanced, non-inductive tokamaks or stellarators.