Search

Start Over Domain Natural Sciences Community Princeton Plasma Physics Laboratory

1. A Reduced Resistive Wall Mode Kinetic Stability Model for Disruption Forecasting

Author(s):
Berkery, J.W.; Sabbagh, S.A.; Bell, R.E.; Gerhardt, S.P.; LeBlanc, B.P.
Abstract:
Kinetic modification of ideal stability theory from stabilizing resonances of mode-particle interaction has had success in explaining resistive wall mode (RWM) stability limits in tokamaks. With the goal of real-time stability forecasting, a reduced kinetic stability model has been implemented in the new Disruption Event Characterization and Forecasting (DECAF) code, which has been written to analyze disruptions in tokamaks. The reduced model incorporates parameterized models for ideal limits on beta, a ratio of plasma pressure to magnetic pressure, which are shown to be in good agreement with DCON code calculations. Increased beta between these ideal limits causes a shift in the unstable region of delta W_K space, where delta W_K is the change in potential energy due to kinetic effects that is solved for by the reduced model, such that it is possible for plasmas to be unstable at intermediate beta but stable at higher beta. Gaussian functions for delta W_K are defined as functions of E cross B frequency and collisionality, with parameters reflecting the experience of the National Spherical Torus Experiment (NSTX). The reduced model was tested on a database of discharges from NSTX and experimentally stable and unstable discharges were separated noticeably on a stability map in E cross B frequency, collisionality space. The reduced model only failed to predict an unstable RWM in 15.6% of cases with an experimentally unstable RWM and performed well on predicting stability for experimentally stable discharges as well.
Type:
Dataset
Issue Date:
May 2017

2. A Riccati Solution for the Ideal MHD Plasma Response with Applications to Real-time Stability Control

Author(s):
Glasser, A.; Kolemen, E.; Glasser, A.H.
Abstract:
To effectuate near real-time feedback control of ideal MHD instabilities in a tokamak geometry, a rapid solution for stability analysis is a prerequisite. Toward this end, we reformulate the δW stability method with a Hamilton-Jacobi theory, elucidating analytical and numerical features of the generic tokamak ideal MHD stability problem. The plasma response matrix is demonstrated to be the solution of an ideal MHD matrix Riccati differential equation (MRDE). Since Riccati equations are prevalent in the control theory literature, such a shift in perspective brings to bear a range of numerical methods that are well-suited to the robust, fast solution of control problems. We discuss the usefulness of Riccati techniques in solving the stiff ODEs often encountered in ideal MHD stability analyses-—for example, in tokamak edge and stellarator physics. We then demonstrate the applicability of such methods to an existing 2D ideal MHD stability code—DCON—enabling its parallel operation in near real-time. Output is shown to match with high accuracy, and with wall-clock time ≪ 1s. Such speed may help enable active feedback ideal MHD stability control, especially in tokamak plasmas whose ideal MHD equilibria evolve with inductive timescale τ > 1s-—as in ITER.
Type:
Dataset
Issue Date:
March 2018

3. A dynamic magnetic tension force as the cause of failed solar eruptions

Abstract:
Coronal mass ejections are solar eruptions driven by a sudden release of magnetic energy stored in the Sun’s corona. In many cases, this magnetic energy is stored in long-lived, arched structures called magnetic flux ropes. When a flux rope destabilizes, it can either erupt and produce a coronal mass ejection or fail and collapse back towards the Sun. The prevailing belief is that the outcome of a given event is determined by a magnetohydrodynamic force imbalance called the torus instability. This belief is challenged, however, by observations indicating that torus-unstable flux ropes sometimes fail to erupt. This contradiction has not yet been resolved because of a lack of coronal magnetic field measurements and the limitations of idealized numerical modelling. Here we report the results of a laboratory experiment that reveal a previously unknown eruption criterion below which torus-unstable flux ropes fail to erupt. We find that such ‘failed torus’ events occur when the guide magnetic field (that is, the ambient field that runs toroidally along the flux rope) is strong enough to prevent the flux rope from kinking. Under these conditions, the guide field interacts with electric currents in the flux rope to produce a dynamic toroidal field tension force that halts the eruption. This magnetic tension force is missing from existing eruption models, which is why such models cannot explain or predict failed torus events.
Type:
Dataset
Issue Date:
11 December 2015

4. A multi-machine scaling of halo current rotation

Author(s):
Myers, C.E.; Eidietis, N.W.; Gerasimov, S.N.; Gerhardt, S.Pl.; Granetz, R.S.; Hender, T.C.; Pautasso, G.
Abstract:
Halo currents generated during unmitigated tokamak disruptions are known to develop rotating asymmetric features that are of great concern to ITER because they can dynamically amplify the mechanical stresses on the machine. This paper presents a multi-machine analysis of these phenomena. More specifically, data from C-Mod, NSTX, ASDEX Upgrade, DIII-D, and JET are used to develop empirical scalings of three key quantities: (1) the machine-specific minimum current quench time, tauCQ; (2) the halo current rotation duration, trot; and (3) the average halo current rotation frequency, <fh>. These data reveal that the normalized rotation duration, trot/tauCQ, and the average rotation velocity, <vh>, are surprisingly consistent from machine to machine. Furthermore, comparisons between carbon and metal wall machines show that metal walls have minimal impact on the behavior of rotating halo currents. Finally, upon projecting to ITER, the empirical scalings indicate that substantial halo current rotation above <fh> = 20 Hz is to be expected. More importantly, depending on the projected value of tauCQ in ITER, substantial rotation could also occur in the resonant frequency range of 6-20 Hz. As such, the possibility of damaging halo current rotation during unmitigated disruptions in ITER cannot be ruled out.
Type:
Dataset
Issue Date:
October 2017

5. A novel scheme for error field correction in permanent magnet stellarators

Author(s):
Rutkowski, Adam; Hammond, Kenneth; Zhu, Caoxiang; Gates, David; Chambliss, Amelia
Abstract:
Stellarators offer a promising path towards fusion reactors, but their design and construction are complicated by stringent tolerance requirements on highly complex 3D coils. A potential way to simplify the engineering requirements for stellarators is to use simple planar toroidal field coils along with permanent magnet arrays to generate shaping fields. In order to ensure sufficient field accuracy while minimizing engineering complexity and system cost, new techniques are required to correct the field produced by the permanent magnet arrays to within requirements set by plasma physics. This work describes a novel correction method developed for this purpose. This analysis is applied to the design of a quasi-axisymmetric stellarator that employs a combination of permanent magnets and planar toroidal field coils to generate its magnetic field. Analysis techniques and initial results using the method for error correction on a proposed permanent magnet stellarator are shown, and it is demonstrated that the method successfully meets the design requirements of the project.
Type:
Dataset
Issue Date:
7 December 2022

6. A scalable real-time framework for Thomson scattering analysis: Application to NSTX-U

Author(s):
F. M. Laggner, A. Diallo, B. P. LeBlanc, R. Rozenblat, G. Tchilinguirian, E.Kolemen, the NSTX-U team
Abstract:
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.
Type:
Dataset
Issue Date:
March 2019

7. A thermodynamic phase transition in magnetic reconnection

Author(s):
Jara-Almonte, Jonathan; Hantao, Ji
Abstract:
Data supporting the manuscript "A thermodynamic phase transition in magnetic reconnection" published in Physical Review Letters.
Type:
Dataset
Issue Date:
7 July 2021

8. Advances in boronization on NSTX-Upgrade

Author(s):
Skinner, C.H.; Bedoya, F.; Scotti, F.; Allain, J.P.; Blanchard, W.; Cai, D.; Jaworski, M.; Koel, B.E.
Abstract:
Boronization has been effective in reducing plasma impurities and enabling access to higher density, higher confinement plasmas in many magnetic fusion devices. The National Spherical Torus eXperiment, NSTX, has recently undergone a major upgrade to NSTX-U in order to develop the physics basis for a ST-based Fusion Nuclear Science Facility (FNSF) with capability for double the toroidal field, plasma current, and NBI heating power and increased pulse duration from 1–1.5 s to 5–8 s. A new deuterated tri-methyl boron conditioning system was implemented together with a novel surface analysis diagnostic. We report on the spatial distribution of the boron deposition versus discharge pressure, gas injection and electrode location. The oxygen concentration of the plasma facing surface was measured by in-vacuo XPS and increased both with plasma exposure and with exposure to trace residual gases. This increase correlated with the rise of oxygen emission from the plasma.
Type:
Dataset
Issue Date:
January 2017

9. Analytic stability boundaries for compressional and global Alfven eigenmodes driven by fast ions. I. Interaction via ordinary and anomalous cyclotron resonances.

Author(s):
Lestz J.B., Gorelenkov N.N., Belova E.V., Tang S.X., Crocker N.A.
Abstract:
Conditions for net fast ion drive are derived for beam-driven, sub-cyclotron compressional (CAE) and global (GAE) Alfven eigenmodes, such as those routinely observed in spherical tokamaks such as NSTX(-U) and MAST. Both co- and counter-propagating CAEs and GAEs are investigated, driven by the ordinary and anomalous Doppler-shifted cyclotron resonance with fast ions. Whereas prior results were restricted to vanishingly narrow distributions in velocity space, broad parameter regimes are identified in this work which enable an analytic treatment for realistic fast ion distributions generated by neutral beam injection. The simple, approximate conditions derived in these regimes for beam distributions of realistic width compare well to the numerical evaluation of the full analytic expressions for fast ion drive. Moreover, previous results in the very narrow beam case are corrected and generalized to retain all terms in omega/omega_{ci} and k_{||}/kperp, which are often assumed to be small parameters but can significantly modify the conditions of drive and damping when they are non-negligible. Favorable agreement is demonstrated between the approximate stability criterion, simulation results, and a large database of NSTX observations of cntr-GAEs.
Type:
Dataset
Issue Date:
September 2019

10. Application of Benchmarked Kinetic Resistive Wall Mode Stability Codes to ITER, Including Additional Physics

Author(s):
Berkery, J.W.; Wang, Z.R.; Sabbagh, S.A.; Liu, Y.Q.; Betti, R.; Guazotto, L.
Abstract:
Leading resistive wall mode (RWM) stability codes MARS-K [Y. Liu, et al., Phys. Plasmas 15, 112503 (2008)] and MISK [B. Hu, et al., Phys. Plasmas 12, 057301 (2005)] have been previously benchmarked. The benchmarking has now been extended to include additional physics, and used to project the stability of ITER in a realistic operating space. Due to ITER's relatively low plasma rotation and collisionality, collisions and non-resonance rotational effects were both found to have little impact on stability, and these non-resonance rotational effects also will not self-consistently affect the ITER RWM eigenfunction. Resonances between thermal ions and electrons and the expected level of ITER toroidal rotation were found to be important to stability, as were alpha particles, which are not in rotational resonance. MISK calculations show that without alpha particles, ITER is projected to be unstable to the RWM, but the expected level of alphas is calculated to provide a sufficient level of stability.
Type:
Dataset
Issue Date:
October 2017

11. Application of IR imaging for free-surface velocity measurement in liquid-metal systems

Author(s):
Hvasta, M.H.; Kolemen, E.; Fisher, A.
Abstract:
Measuring free-surface, liquid-metal flow velocity is challenging to do in a reliable and accurate manner. This paper presents a non-invasive, easily-calibrated method of measuring the surface velocities of open-channel liquid-metal flows using an IR camera. Unlike other spatially-limited methods, this IR camera particle tracking technique provides full field-of-view data that can be used to better understand open-channel flows and determine surface boundary conditions. This method could be implemented and automated for a wide range of liquid-metal experiments, even if they operate at high-temperatures or within strong magnetic fields.
Type:
Dataset
Issue Date:
January 2017

12. Application of IR imaging for free-surface velocity measurement in liquid-metal systems

Author(s):
Hvasta, M. G.; Kolemen, E.; Fisher, A.
Abstract:
Measuring free-surface, liquid-metal flow velocity is challenging to do in a reliable and accurate manner. This paper presents a non-invasive, easily-calibrated method of measuring the surface velocities of open-channel liquid-metal flows using an IR camera. Unlike other spatially-limited methods, this IR camera particle tracking technique provides full field-of-view data that can be used to better understand open-channel flows and determine surface boundary conditions. This method could be implemented and automated for a wide range of liquid-metal experiments, even if they operate at high-temperatures or within strong magnetic fields.
Type:
Dataset
Issue Date:
January 2017

13. Application of transient CHI plasma startup to future ST and AT devices

Author(s):
Hammond, K.C.; Raman, R.; Jardin, S.C.
Abstract:
Employment of non-inductive plasma start-up techniques would considerably simplify the design of a spherical tokamak fusion reactor. Transient coaxial helicity injection (CHI) is a promising method, expected to scale favorably to next-step reactors. However, the implications of reactor-relevant parameters on the initial breakdown phase for CHI have not yet been considered. Here, we evaluate CHI breakdown in reactor-like configurations using an extension of the Townsend avalanche theory. We find that a CHI electrode concept in which the outer vessel wall is biased to achieve breakdown, while previously successful on NSTX and HIT-II, may exhibit a severe weakness when scaled up to a reactor. On the other hand, concepts which employ localized biasing electrodes such as those used in QUEST would avoid this issue. Assuming that breakdown can be successfully attained, we then apply scaling relationships to predict plasma parameters attainable in the transient CHI discharge. Assuming the use of 1 Wb of injector flux, we find that plasma currents of 1 MA should be achievable. Furthermore, these plasmas are expected to Ohmically self-heat with more than 1 MW of power as they decay, facilitating efficient hand-off to steady-state heating sources. These optimistic scalings are supported by TSC simulations.
Type:
Dataset
Issue Date:
February 2019

14. Blob Structure and Motion in the Edge and SOL of NSTX

Author(s):
Zweben, S.J.; J.R. Myra; W.M. Davis; D.A. D'Ippolito; T.K. Gray; S.M. Kaye; B.P. LeBlanc; R.J. Maqueda; D.A. Russell; D.P. Stotler
Abstract:
Blob analysis dataset
Type:
Dataset
Issue Date:
January 2016

15. Blob-hole correlation model for edge turbulence and comparisons with NSTX GPI data

Author(s):
Myra, J.R.; Zweben, S.J.; Russell, D.A.
Abstract:
Gas puff imaging (GPI) observations made in NSTX [Zweben S J, et al., 2017 Phys. Plasmas 24 102509] have revealed two-point spatial correlations of edge and scrape-off layer turbulence in the plane perpendicular to the magnetic field. A common feature is the occurrence of dipole-like patterns with significant regions of negative correlation. In this paper, we explore the possibility that these dipole patterns may be due to blob-hole pairs. Statistical methods are applied to determine the two-point spatial correlation that results from a model of blob-hole pair formation. It is shown that the model produces dipole correlation patterns that are qualitatively similar to the GPI data in several respects. Effects of the reference location (confined surfaces or scrape-off layer), a superimposed random background, hole velocity and lifetime, and background sheared flows are explored and discussed with respect to experimental observations. Additional analysis of the experimental GPI dataset is performed to further test this blob-hole correlation model. A time delay two-point spatial correlation study did not reveal inward propagation of the negative correlation structures that were postulated to correspond to holes in the data nor did it suggest that the negative correlation structures are due to neutral shadowing. However, tracing of the highest and lowest values (extrema) of the normalized GPI fluctuations shows strong evidence for mean inward propagation of minima and outward propagation of maxima, in qualitative agreement with theoretical expectations. Other properties of the experimentally observed extrema are discussed.
Type:
Dataset
Issue Date:
July 2018

16. Calibrationless rotating Lorentz-force flowmeters for low flow rate applications

Author(s):
Hvasta, M. G.; Dudt, D.; Fisher, A. E.; Kolemen, E.
Abstract:
A 'weighted magnetic bearing' has been developed to improve the performance of rotating Lorentz-force flowmeters (RLFFs). Experiments have shown that the new bearing reduces frictional losses within a double-sided, disc-style RLFF to negligible levels. Operating such an RLFF under 'frictionless' conditions provides two major benefits. First, the steady-state velocity of the RLFF magnets matches the average velocity of the flowing liquid at low flow rates. This enables an RLFF to make accurate volumetric flow measurements without any calibration or prior knowledge of the fluid properties. Second, due to minimized frictional losses, an RLFF is able to measure low flow rates that cannot be detected when conventional, high-friction bearings are used. This paper provides a brief background on RLFFs, gives a detailed description of weighted magnetic bearings, and compares experimental RLFF data to measurements taken with a commercially available flowmeter.
Type:
Dataset
Issue Date:
29 May 2018

17. Collisional dependence of Alfven mode saturation in tokamaks

Author(s):
Zhou, M.; White, R.
Abstract:
Saturation of \alfven modes driven unstable by a distribution of high energy particles as a function of collisionality is investigated with a guiding center code, using numerical eigenfunctions produced by linear theory and numerical high energy particle distributions. The most important resonance is found and it is shown that when the resonance domain is bounded, not allowing particles to collisionlessly escape, the saturation amplitude is given by the balance of the resonance mixing time with the time for nearby particles to collisionally diffuse across the resonance width. Saturation amplitudes are in agreement with theoretical predictions as long as the mode amplitude is not so large that it produces stochastic loss from the resonance domain.
Type:
Dataset
Issue Date:
December 2016

18. Comment on ‘Numerical modeling of tokamak breakdown phase driven by pure Ohmic heating under ideal conditions’

Author(s):
Yoo, Min-Gu; Na, Yong-Su
Abstract:
In this comment, we point out possible critical numerical flaws of recent particle simulation studies (Jiang et al 2016 Nucl. Fusion 56 126017, Peng et al 2018 Nucl. Fusion 58 026007) on the electrical gas breakdown in a simple one-dimensional periodic slab geometry. We show that their observations on the effects of the ambipolar electric fields during the breakdown, such as the sudden reversal of the ion flow direction, could not be real physical phenomena but resulting from numerical artifacts violating the momentum conservation law. We show that an incomplete implementation of the direct-implicit scheme can cause the artificial electric fields and plasma transports resulting in fallacies in simulation results. We also discuss that their simple plasma model without considering poloidal magnetic fields seriously mislead the physical mechanism of the electrical gas breakdown because it cannot reflect important dominant plasma dynamics in the poloidal plane (Yoo et al 2018 Nat. Commun. 9 3523).
Type:
Dataset
Issue Date:
June 2019

19. Compact and multi-view solid state neutral particle analyzer arrays on National Spherical Torus Experiment-Upgrade

Author(s):
Liu, D.; Heidbrink, W.W.; Tritz, K.; Fredrickson, E.D.; Hao, G.Z.; Zhu, Y.B.
Abstract:
A compact and multi-view Solid State Neutral Particle Analyzer (SSNPA) diagnostic based on silicon photodiode arrays has been successfully tested on the National Spherical Torus Experiment-Upgrade (NSTX-U). The SSNPA diagnostic provides spatially, temporally, and pitch-angle resolved measurements of fast-ion distribution by detecting fast neutral flux resulting from charge exchange (CX) reactions. The system consists of three 16-channel subsystems: t-SSNPA viewing the plasma mid-radius and neutral beam (NB) line #2 tangentially, r-SSNPA viewing the plasma core and NB line #1 radially and p-SSNPA with no intersection with any NB lines. Due to the setup geometry, the active CX signals of t-SSNPA and r-SSNPA are mainly sensitive to passing and trapped particles respectively. In addition, both t-SSNPA and r-SSNPA utilize three vertically stacked arrays with different filter thickness to obtain coarse energy information. The experimental data show that all channels are operational. The signal to noise ratio is typically larger than 10 and the main noise is x-ray induced signal. The active and passive CX signals are clearly observed on t-SSNPA and r-SSNPA during NB modulation. The SSNPA data also indicate significant losses of passing particles during sawteeth, while trapped particles are weakly affected. Fluctuations up to 120 kHz, have been observed on SSNPA, and they are strongly correlated with magnetohydrodynamics (MHD) instabilities.
Type:
Dataset
Issue Date:
November 2016

20. Compact steady-state tokamak performance dependence on magnet and core physics limits

Author(s):
Menard, J.E.
Abstract:
Compact tokamak fusion reactors utilizing advanced high-temperature superconducting magnets for the toroidal field coils have received considerable recent attention due to the promise of more compact devices and more economical fusion energy development. Facilities with combined Fusion Nuclear Science (FNS) and Pilot Plant missions to provide both the nuclear environment needed to develop fusion materials and components while also potentially achieving sufficient fusion performance to generate modest net electrical power are considered. The performance of the tokamak fusion system is assessed using a range of core physics and toroidal field magnet performance constraints to better understand which parameters most strongly influence the achievable fusion performance.
Type:
Dataset
Issue Date:
December 2018

Limit your search

View distribution

Current results range from 2015 to 2023

View larger »