Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “NSTX”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

Modeling of lithium granule injection in NSTX using M3D-C1

In this paper, we present 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 (Jardin et al 2012 Comput. Sci. Discovery 5 014002), allowing the simulation of realistic lithium granule injections. 2D and 3D simulations of Li injections in NSTX H-mode plasmas are performed and the effect of granule size, injection angle and velocity on the pedestal gradient increase is studied. The amplitude of the local pressure perturbation caused by the granules is found to be highly dependent on the solid granule size. Adjusting the granule injection velocity allows one to inject more particles at the pedestal top. 3D simulations show the destabilization of high order MHD modes whose amplitude is directly linked to the localized pressure perturbation, which is found to depend on the toroidal localization of the granule density source.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Energy exchange dynamics across L–H transitions in NSTX

Here, we studied the energy exchange dynamics across the low-to-high-confinement (L–H) transition in NSTX discharges using the gas-puff imaging (GPI) diagnostic. The investigation focused on the energy exchange between flows and turbulence to help clarify the mechanism of the L–H transition. We applied this study to three types of heating schemes, including a total of 17 shots from the NSTX 2010 campaign run. Results show that the edge fluctuation characteristics (fluctuation levels, radial and poloidal correlation lengths) measured using GPI do not vary just prior to the H-mode transition, but change after the transition. Using a velocimetry approach (orthogonal-dynamics programming), velocity fields of a $24\times 30$ cm GPI view during the L–H transition were obtained with good spatial (~1 cm) and temporal (~2.5 μs) resolutions. Analysis using these velocity fields shows that the production term is systematically negative just prior to the L–H transition, indicating a transfer from mean flows to turbulence, which is inconsistent with the predator–prey paradigm. Moreover, the inferred absolute value of the production term is two orders of magnitude too small to explain the observed rapid L–H transition. These discrepancies are further reinforced by consideration of the ratio between the kinetic energy in the mean flow to the thermal free energy, which is estimated to be much less than 1, suggesting again that the turbulence depletion mechanism may not play an important role in the transition to the H-mode. Although the Reynolds work therefore appears to be too small to directly deplete the turbulent free energy reservoir, order-of-magnitude analysis shows that the Reynolds stress may still make a non-negligible contribution to the observed poloidal flows.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Multi-species impurity granule injection and mass deposition projections in NSTX-U discharges

Here, by employing a neutral gas shielding (NGS) model to characterize impurity granule injection, the ablation rates 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 quantify the ablation rate 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 s –1 are presented for NSTX-U L-mode type plasmas, as well as H-mode discharges with low natural ELM frequency. Additionally, ablation calculations of 500 micron granules of each species are presented at velocities ranging from 50–150 m s –1 . In H-mode discharges these simulations show that the majority of the injected granule is ablated within or just past the edge steep gradient region. 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 ELM crashes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Study of the likelihood of Alfvénic mode bifurcation in NSTX and predictions for ITER baseline scenarios

Rare Alfvénic wave transitions between fixed-frequency and chirping phases are identified in NSTX, where Alfvénic waves are normally observed to exhibit either chirping or avalanching responses. For those transitions, we apply in this paper a criterion (Duarte et al 2017 Nucl. Fusion 57 054001) to predict the nature of fast ion redistribution in tokamaks to be in the convective or diffusive nonlinear regimes. For NSTX discharges in which the transition is not accompanied by changes in the beam deposited power or modifications in the injected radiofrequency power, it has been found that the anomalous fast ion transport is a likely mediator of the bifurcation between the fixed-frequency mode behavior and rapid chirping. For a quantitative assessment, global gyrokinetic simulations of the effects of electrostatic ion temperature gradient turbulence and trapped electron mode turbulence on chirping were pursued using the GTS code. The investigation is extended by means of predictive studies of the probable spectral behavior of Alfvénic eigenmodes for baseline ITER cases consisting of elmy, advanced and hybrid scenarios. Finally, it has been observed that most modes are found to be borderline between the steady and the chirping phases.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Effect of sawtooth crashes on fast ion distribution in NSTX-U

During the 2016 experimental campaign of NSTX-U (Menard et al 2012 Nucl. Fusion 52 083015), long L-mode and reproducible sawtoothing plasmas have been achieved that were previously not accessible before the upgrade. This provides a good opportunity to study the effect of sawtooth crashes on fast ion confinement and redistribution in spherical tokamaks. The solid-state neutral particle analyzer (SSNPA) and fast ion D-alpha (FIDA) diagnostics on NSTX-U each have a tangentially-viewing instrument and a radially- or vertically-viewing instrument, which are mainly sensitive to passing and trapped fast ions, respectively. It has been observed on both diagnostics that passing particles are strongly expelled from the plasma core to the edge during sawtooth crashes, while trapped fast ions are weakly affected. The tangentially-viewing SSNPA observes large signal spikes at the sawtooth crashes because fast ions move to the edge and charge exchange with edge neutrals. The radially-viewing SSNPA data suggest that there is a small drop of trapped particles in the core. The tangentially-viewing FIDA (t-FIDA) system observes a depletion as large as 25% in the region inside the inversion radius, while an increase at the outer region. There is almost no change in the signals of the vertically-viewing FIDA system. The neutron emission can drop as much as 15% at the sawtooth crashes, accompanied by an increase of edge D α light. Simulations with the Kadomstev and Porcelli sawtooth models have also been performed and compared with the measurements. The full reconnection Kadomtsev model overestimates the neutron rate drop at each sawtooth crash. When tuning the sawtooth input parameters, the partial reconnection Porcelli sawtooth model can qualitatively reproduce the neutron rate drop and t-FIDA signal drop in the core, but it fails to predict the t-FIDA signal increase at the edge.

43 PARTICLE ACCELERATORS↗

Divertor leg filaments in NSTX-U

This study presents the first observation of divertor localized turbulence in the NSTX-Upgrade spherical tokamak. Previous work in NSTX discharges described fluctuations on the divertor target due to blobs generated on the low field side midplane. In the region that shows disconnection from upstream turbulence (private flux region and proximity of the strike point), divertor-localized fluctuations are observed via imaging of C III and D-α emission in diverted L-mode discharges and are characterized in this paper. Field-aligned filaments connected to the divertor target plate are radially localized at the separatrix on the outer divertor leg and in the private flux region on the inner divertor leg. These are limited to the region below the X-point with fluctuation levels up to 10%–20%. The filaments have comparable poloidal and radial correlation lengths (10–100 ion gyroradii ρ i ) and parallel correlation lengths of several meters. Toroidal mode numbers are in the range of 10–30 and 2–10 for outer and inner leg filaments. The poloidal motion of outer leg filaments has the same direction and comparable magnitude to advection due to E x B drift calculated by the multi-fluid edge transport code UEDGE with inclusion of cross field drifts. Disconnection between inner and outer leg filaments as well as the absence of correlation with upstream turbulence support the hypothesis of X-point disconnection for divertor leg filaments. Simulations with the ArbiTER linear eigenvalue solver and a resistive-balooning model were performed with a simulation grid limited to the divertor legs. Instabilities localized to the bad curvature side of the divertor legs were observed, in qualitative agreement with the experiment.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Initial transport and turbulence analysis and gyrokinetic simulation validation in NSTX-U L-mode plasmas

Transport analysis, ion-scale turbulence measurements, and initial linear and nonlinear gyrokinetic simulations are reported for a transport validation study based on low aspect ratio NSTX-U L-mode discharges. The relatively long, stationary L-modes enabled by the upgraded centerstack provide a more ideal target for transport validation studies that were not available during NSTX operation. Transport analysis shows that anomalous electron transport dominates energy loss while ion thermal transport is well described by neoclassical theory. Linear gyrokinetic GYRO analysis predicts that ion temperature gradient (ITG) modes are unstable around normalized radii ρ=0.6-0.8, although E×B shearing rates are larger than the linear growth rates over much of that region. Deeper in the core (ρ=0.4-0.6), electromagnetic microtearing modes (MTM) are unstable as a consequence of the relatively high beta and collisionality in these particular discharges. Consistent with the linear analysis, local, nonlinear ion-scale GYRO simulations predict strong ITG transport at ρ=0.76, whereas electromagnetic MTM transport is important at ρ=0.47. The prediction of ion-scale turbulence is consistent with 2D beam emission spectroscopy (BES) that measures the presence of broadband ion-scale fluctuations. Interestingly, the BES measurements also indicate the presence of bi-modal poloidal phase velocity propagation that could be indicative of two different turbulence types. Furthermore, in the region between (ρ=0.56, 0.66), ion-scale simulations are strongly suppressed by the locally large E×B shear. Instead, electron temperature gradient (ETG) turbulence simulations predict substantial transport, illustrating electron-scale contributions can be important in low aspect ratio L-modes, similar to recent analysis at conventional aspect ratio. However, agreement within experimental uncertainties has not been demonstrated, which requires additional simulations to test parametric sensitivities. The potential need to include profile-variation effects (due to the relatively large value of ρ∗=ρ i /a at low aspect ratio), including electromagnetic and possibly multi-scale effects, is also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling of carbon pellets disruption mitigation in an NSTX-U plasma

Single carbon pellet disruption mitigation simulations using M3D-C1 were conducted for 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.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Real-time capable modeling of ICRF heating on NSTX and WEST via machine learning approaches

Abstract A real-time capable core Ion Cyclotron Range of Frequencies (ICRF) heating model on NSTX and WEST is developed. The model is based on two nonlinear regression algorithms, the random forest ensemble of decision trees and the multilayer perceptron neural network. The algorithms are trained on TORIC ICRF spectrum solver simulations of the expected flat-top operation scenarios in NSTX and WEST assuming Maxwellian plasmas. The surrogate models are shown to successfully capture the multi-species core ICRF power absorption predicted by the original model for the high harmonic fast wave and the ion cyclotron minority heating schemes while reducing the computational time by six orders of magnitude. Although these models can be expanded, the achieved regression scoring, computational efficiency and increased model robustness suggest these strategies can be implemented into integrated modeling frameworks for real-time control applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

A pulse-burst laser system for Thomson scattering on NSTX-U

A pulse-burst laser system has been built for Thomson scattering on NSTX-U, and is currently being integrated into the NSTX-U Thomson scattering diagnostic system. The laser will be operated in three distinct modes. The base mode is continuous 30 Hz rep rate, and is the standard operating mode of the laser. The base mode will be interrupted to produce a "slow burst" (specified 1 kHz rep rate for 50 ms) or a "fast burst" (specified 10 kHz rep rate for 5 ms). The combination of base mode→ interruption→ burst mode is new and has not been implemented on any previous pulse-burst laser system. Laser pulsing is halted for a set period (~ 1 minute) following a burst to allow the YAG rods to cool; this type of operation is called a heat-capacity laser. The laser is Nd:YAG operated at 1064 nm, q-switched to produce ≥ 1.5 J pulses with ~ 20 ns FWHM. It is flashlamp pumped, with dual-rod oscillator (9 mm) and dual-rod amplifier (12 mm). Variable pulsewidth drive of the flashlamps is accomplished by IGBT (insulated gate bipolar transistor) switching of electrolytic capacitor banks. Direct control of the laser Pockels cell drive enables optimal pulse energy extraction. The laser system has demonstrated compliance with all specifications, and is capable of exceeding design specifications by significant margins, e.g., higher rep rates for longer burst periods. Burst operation of this laser system will be used to capture fast time evolution of the electron temperature and density profiles during events such as ELMs, the L-H transition, and various MHD modes.

47 OTHER INSTRUMENTATION↗

A poloidal high- k scattering system for NSTX-U

A previous 5-channel tangential high-k scattering system is being replaced by an 8-channel, poloidal high-k scattering system on the National Spherical Torus eXperiment Upgrade (NSTX-U) device located in Princeton, NJ, USA. The 693 GHz poloidal scattering system replaces a 280 GHz tangential scattering system to study high-k electron density fluctuations on NSTX-U, thereby considerably enhancing planned turbulence physics studies by providing a measurement of the k θ -spectrum of both electron temperature gradient (ETG) and ion temperature gradient (ITG) modes. Two approaches to generating the 693 GHz probe beam are under development: an optically-pumped far-infrared (FIR) laser that generates ~50 mW, and a compact gyrotron that can potentially generate in excess of 5 W. Large aperture optics collect radiation scattered from density fluctuations in the plasma core at 8 simultaneous scattering angles ranging from 2 to 15° corresponding to poloidal wavenumbers that extend to >40 cm -1 . Finally, steerable launch optics coupled with receiver optics mounted on a 5-axis receiver carriage allow the scattering volume to be placed radially from r/a = 0.3 out to the pedestal region (r/a ~ 0.99) and translated horizontally as needed to satisfy wavenumber matching.

47 OTHER INSTRUMENTATION↗

NSTX-U Advances in Real-Time C++11 on Linux

Programming languages like C and Ada combined with proprietary embedded operating systems have dominated the real-time application space for decades. The new C++11standard includes native, language-level support for concurrency, a required feature for any nontrivial event-oriented real-time software. Threads, Locks, and Atomics now exist to provide the necessary tools to build the structures that make up the foundation of a complex real-time system. The National Spherical Torus Experiment Upgrade (NSTX-U) at the Princeton Plasma Physics Laboratory (PPPL) is breaking new ground with the language as applied to the needs of fusion devices. A new Digital Coil Protection System (DCPS) will serve as the main protection mechanism for the magnetic coils, and it is written entirely in C++11 running on Concurrent Computer Corporation's real-time operating system, RedHawk Linux. It runs over 600 algorithms in a 5 kHz control loop that determine whether or not to shut down operations before physical damage occurs. To accomplish this, NSTX-U engineers developed software tools that do not currently exist elsewhere, including real-time atomic synchronization, real-time containers, and a real-time logging framework. Together with a recent (and carefully configured) version of the GCC compiler, these tools enable data acquisition, processing, and output using a conventional operating system to meet a hard real-time deadline (that is, missing one periodic is a failure) of 200 microseconds.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Virtual Tokamak for Test and Development of Plasma Control Applied to NSTX-U

Safe and efficient operation of tokamak experiments depends on shot preparation prior to the experiment to maximize performance and avoid issues that are predictable consequences of known physics. The General Atomics TokSys toolbox is used to test and develop the plasma control system (PCS) on numerous tokamaks around the world. Tokamaks that use a version of the DIII-D PCS can connect it to a TokSys simulation and control a virtual version of the tokamak with the real PCS. Recent upgrades to the simulation include more detailed profile and transport modeling to support design and implementation of optimized scenarios. The simulation combines a module called “Profiles” that simulates the 1D profiles of density, pressure, and current and a module called “GSevolve'' that evolves the 2D Grad-Shafranov equilibrium. The “Profiles'' module primarily uses simplified models for fueling, heating, and current. The TRANSP code provides precomputed diffusion coefficients and heating/current drive profiles from radio-frequency and neutral beam sources for the scenarios being simulated. The fidelity of a simulation can be enhanced by analyzing it with TRANSP and making corrections to the precomputed values and iterating until convergence. Simulations of the National Spherical Torus Upgrade (NSTX-U) are presented. It is shown that the simulation reproduces the experimental profile evolution. Here, the simulation has also been upgraded to reproduce the experimental outcome at points of bifurcation. This is demonstrated with a simulation of vertical displacement in NSTX-U.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quasi-linear gyrokinetic predictions of the Coriolis momentum pinch in NSTX

This paper presents quasi-linear gyrokinetic predictions of the Coriolis momentum pinch for low aspect-ratio NSTX H-modes where previous experimental measurements were focused. Local, linear calculations predict that in the region of interest (just outside the mid-radius) of these relatively high-beta plasmas, profiles are most unstable to microtearing modes that are only effective in transporting electron energy. However, sub-dominant electromagnetic and electrostatic ballooning modes are also unstable, which are effective at transporting energy, particles and momentum. The quasi-linear prediction of transport from these weaker ballooning modes, assuming they contribute transport in addition to that from microtearing modes in a nonlinear turbulent state, leads to a very small or outward convection of momentum, inconsistent with the experimentally measured inward pinch, and opposite to predictions in conventional aspect ratio tokamaks. Additional predictions of a low beta L-mode plasma, unstable to more traditional electrostatic ion temperature gradient-trapped electron mode instability, show that the Coriolis pinch is inward but remains relatively weak and insensitive to many parameter variations. The weak or outward pinch predicted in NSTX plasmas appears to be at least partially correlated to changes in the parallel mode structure that occur at finite beta and low aspect ratio, as discussed in previous theories. The only conditions identified where a stronger inward pinch is predicted occur either in the purely electrostatic limit or if the aspect ratio is increased. As the Coriolis pinch cannot explain the measured momentum pinch, additional theoretical momentum transport mechanisms are discussed that may be potentially important.

Tokamak Turbulence Transport Momentum pinch↗

Hydrogen Retention in Lithium on Metallic Walls from “In Vacuo” Analysis in LTX and Implications for High-Z Plasma-Facing Components in NSTX-U

The application of lithium to plasma-facing components (PFCs) has long been used as a technique for wall conditioning in magnetic confinement devices to improve plasma performance. Determining the characteristics of PFCs at the time of exposure to the plasma, however, is difficult because they can only be analyzed after venting the vacuum vessel and removing them at the end of an operational period. The Materials Analysis and Particle Probe (MAPP) addresses this problem by enabling PFC samples to be exposed to plasmas, and then withdrawn into an analysis chamber without breaking vacuum. The MAPP system was used to introduce samples that matched the metallic PFCs of the Lithium Tokamak Experiment (LTX). Lithium that was subsequently evaporated onto the walls also covered the MAPP samples, which were then subject to LTX discharges. In vacuo extraction and analysis of the samples indicated that lithium oxide formed on the PFCs, but improved plasma performance persisted in LTX. The reduced recycling this suggests is consistent with separate surface science experiments that demonstrated deuterium retention in the presence of lithium oxide films. Since oxygen decreases the thermal stability of the deuterium in the film, the release of deuterium was observed below the lithium deuteride dissociation temperature. This may explain what occurred when lithium was applied to the surface of the NSTX Liquid Lithium Divertor (LLD). The LLD had segments with individual heaters, and the deuterium-alpha emission was clearly lower in the cooler regions. The plan for NSTX-U is to replace the graphite tiles with high-Z PFCs, and apply lithium to their surfaces with lithium evaporation. Experiments with lithium coatings on such PFCs suggest that deuterium could still be retained if lithium compounds form, but limiting their surface temperatures may be necessary.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Modeling of Lithium Granule Injection in NSTX with M3D-C1

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.

Fil, A↗

Elemental and topographical imaging of microscopic variations in deposition on NSTX-U and DIII-D samples

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.

Surface analysis↗

Suppression of Alfvén modes on NSTX-U with outboard beam injection

In this paper we present data from experiments on NSTX-U where it is shown for the first time that small amounts of high pitch-angle beam ions can strongly suppress the counter-propagating Global Alfvén Eigenmodes (GAE). GAE have been implicated in the redistribution of fast ions and modification of the electron power balance in previous experiments on NSTX. The ability to predict the stability of Alfvén modes, and developing methods to control them, is important for fusion reactor like the International Tokamak Experimental Reactor (ITER) which are heated by a large population of non-thermal, super-Alfvénic ions consisting of fusion generated alphas and beam ions injected for current profile control. We present a qualitative interpretation of these observations using an analytic model of the Doppler-shifted ion-cyclotron resonance drive responsible for GAE instability which has an important dependence on k⊥ρL. A quantitative analysis of this data with the HYM stability code predicts both the frequencies and instability of the GAE prior to, and suppression of the GAE after the injection of high pitch-angle beam ions.

Energetic particles↗