Search

Start Over Domain Natural Sciences

101. First impurity powder injection experiments in LHD

Author(s):
Nespoli F., Ashikawa N., Gilson E.P., Lunsford R., Masuzaki S., Shoji M., Oishi T., Suzuki C., Nagy A., Mollen A., Pablant N.A., Ida K., Yoshinuma M., Tamura N., Gates D.A., Morisaki T., and the LHD experiment group
Abstract:
Injection of impurities in the form of sub-millimeter powder grains is performed for the first time in the Large Helical Device (LHD) plasma, employing the Impurity Powder Dropper (IPD) [A. Nagy et al., RSI 2018], developed and built by PPPL. Controlled amounts of boron (B) and boron nitride (BN) powder are injected into the helical plasma. Visible camera imaging, UV and charge exchange spectroscopy measurements show that the injected impurities effectively penetrate into the plasma in two different magnetic configurations.The prompt effects of the impurities on the plasma are characterized as the injection rate is scanned. The injected impurities provide a supplemental electron source, causing the plasma density to increase, together with the radiated power. Beneficial effects on the confined plasma temperature are observed at low plasma densities, due to an increased efficiency in NBI power absorption. For $n_{e,av}<10^{19}m^{-3}$ the powder grains penetrate deeper into the plasma, as they can be less effectively deflected by the plasma flow in the divertor leg, which they have to cross first as they are injected from the top of the machine.In this case, the created B ions are observed to move outwards from UV spectroscopy and charge exchange measurements, due to the outwards direction of the radial electric field. This makes low density plasmas a better candidate for powder boronization techniques.
Type:
Dataset
Issue Date:
November 2020

102. Geometric concepts for stellarator permanent magnet arrays

Author(s):
Hammond, K. C.; Zhu, C.; Brown, T.; Corrigan, K.; Gates, D. A.; Sibilia, M.
Abstract:
The development of stellarators that use permanent magnet arrays to shape their confining magnetic fields has been a topic of recent interest, but the requirements for how such magnets must be shaped, manufactured, and assembled remain to be determined. To address these open questions, we have performed a study of geometric concepts for magnet arrays with the aid of the newly developed MAGPIE code. A proposed experiment similar to the National Compact Stellarator Experiment (NCSX) is used as a test case. Two classes of magnet geometry are explored: curved bricks that conform to a regular grid in cylindrical coordinates, and hexahedra that conform to the toroidal plasma geometry. In addition, we test constraints on the magnet polarization. While magnet configurations constrained to be polarized normally to a toroidal surface around the plasma are unable to meet the required magnetic field parameters when subject to physical limitations on the strength of present-day magnets, configurations with unconstrained polarizations are shown to satisfy the physics requirements for a targeted plasma.
Type:
Dataset
Issue Date:
July 2020

103. Gyrokinetic understanding of the edge pedestal transport driven by resonant magnetic perturbations in a realistic divertor geometry

Author(s):
Hager, R.; Chang, C. S.; Ferraro, N. M.; Nazikian R.
Abstract:
Self-consistent simulations of neoclassical and electrostatic turbulent transport in a DIII-D H-mode edge plasma under resonant magnetic perturbations (RMPs) have been performed using the global total-f gyrokinetic particle-in-cell code XGC, in order to study density-pump out and electron heat confinement.The RMP field is imported from the extended magneto-hydrodynamics (MHD) code M3D-C1, taking into account the linear two-fluid plasma response.With both neoclassical and turbulence physics considered together, the XGC simulation reproduces two key features of experimentally observed edge transport under RMPs: increased radial particle transport in the pedestal region that is sufficient to account for the experimental pump-out rate, and suppression of the electron heat flux in the steepest part of the edge pedestal.In the simulation, the density fluctuation amplitude of modes moving in the electron diamagnetic direction increases due to interaction with RMPs in the pedestal shoulder and outward, while the electron temperature fluctuation amplitude decreases.
Type:
Dataset
Issue Date:
June 2020

104. MHD-blob correlations in NSTX

Author(s):
Zweben SJ; Fredrickson ED; Myra JR; Podesta M; Scotti F
Abstract:
This paper describes a study of the cross-correlations between edge fluctuations as seen in the gas puff imaging (GPI) diagnostic and low frequency coherent magnetic fluctuations (MHD) in H-mode plasmas in NSTX. The main new result was that large blobs in the SOL were significantly correlated with MHD activity the 3-6 kHz range in 21 of the 223 shots examined. There were also many other shots in which fluctuations in the GPI signal level and its peak radius Rpeak were correlated with MHD activity, but without any significant correlation of the MHD with large blobs. The structure and motion of the MHD is compared with that of the correlated blobs, and some possible theoretical mechanisms for the MHD-blob correlation are discussed.
Type:
Dataset
Issue Date:
May 2020

105. Machine Learning Characterization of Alfvénic and Sub-Alfvénic Chirping and Correlation With Fast-Ion Loss at NSTX

Author(s):
Woods, B. J. Q.; Duarte, V. N.; Fredrickson, E. D.; Gorelenkov, N. N.; Podestà, M.; Vann, R. G. L.
Abstract:
Abrupt large events in the Alfvenic and sub-Alfvenic frequency bands in tokamaks are typically correlated with increased fast-ion loss. Here, machine learning is used to speed up the laborious process of characterizing the behavior of magnetic perturbations from corresponding frequency spectrograms that are typically identified by humans. The analysis allows for comparison between different mode character (such as quiescent, fixed frequency, and chirping, avalanching) and plasma parameters obtained from the TRANSP code, such as the ratio of the neutral beam injection (NBI) velocity and the Alfven velocity (v_inj./v_A), the q-profile, and the ratio of the neutral beam beta and the total plasma beta (beta_beam,i / beta). In agreement with the previous work by Fredrickson et al., we find a correlation between beta_beam,i and mode character. In addition, previously unknown correlations are found between moments of the spectrograms and mode character. Character transition from quiescent to nonquiescent behavior for magnetic fluctuations in the 50200-kHz frequency band is observed along the boundary v_phi ~ (1/4)(v_inj. - 3v_A), where v_phi is the rotation velocity.
Type:
Dataset
Issue Date:
December 2019

106. March Mathness: Effects of basketball on the brain

Author(s):
Antony, James; McDougle, Sam
Abstract:
Surprise signals a discrepancy between past and current beliefs. It is theorized to be linked to affective experiences, the creation of particularly resilient memories, and segmentation of the flow of experience into discrete perceived events. However, the ability to precisely measure naturalistic surprise has remained elusive. We used advanced basketball analytics to derive a quantitative measure of surprise and characterized its behavioral, physiological, and neural correlates in human subjects observing basketball games. We found that surprise was associated with segmentation of ongoing experiences, as reflected by subjectively perceived event boundaries and shifts in neocortical patterns underlying belief states. Interestingly, these effects differed by whether surprising moments contradicted or bolstered current predominant beliefs. Surprise also positively correlated with pupil dilation, activation in subcortical regions associated with dopamine, game enjoyment, and long-term memory. These investigations support key predictions from event segmentation theory and extend theoretical conceptualizations of surprise to real-world contexts.
Type:
Dataset
Issue Date:
2020

107. Modeling of resistive plasma response in toroidal geometry using an asymptotic matching approach

Author(s):
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.
Type:
Dataset
Issue Date:
October 2020

108. Simulating pitch angle scattering using an explicitly solvable energy-conserving algorithm

Author(s):
Zhang, Xin; Fu, Yichen; Qin, Hong
Abstract:
Particle distribution functions evolving under the Lorentz operator can be simulated with the Langevin equation for pitch angle scattering. This approach is frequently used in particle based Monte-Carlo simulations of plasma collisions, among others. However, most numerical treatments do not guarantee energy conservation, which may lead to unphysical artifacts such as numerical heating and spectra distortions. We present a novel structure-preserving numerical algorithm for the Langevin equation for pitch angle scattering. Similar to the well-known Boris algorithm, the proposed numerical scheme takes advantage of the structure-preserving properties of the Cayley transform when calculating the velocity-space rotations. The resulting algorithm is explicitly solvable, while preserving the norm of velocities down to machine precision. We demonstrate that the method has the same order of numerical convergence as the traditional stochastic Euler-Maruyama method.
Type:
Dataset
Issue Date:
September 2020

109. Solitary zonal structures in subcritical drift waves: a minimum model

Author(s):
Yao Zhou; Hongxuan Zhu; I. Y. Dodin
Abstract:
{\rtf1\ansi\ansicpg1252\cocoartf1561\cocoasubrtf610{\fonttbl\f0\fswiss\fcharset0 Helvetica;}{\colortbl;\red255\green255\blue255;\red0\green0\blue0;}{\*\expandedcolortbl;;\cssrgb\c0\c0\c0;}\margl1440\margr1440\vieww10800\viewh8400\viewkind0\pard\tx887\tx1775\tx2662\tx3550\tx4438\tx5325\tx6213\tx7101\tx7988\tx8876\tx9764\tx10651\tx11539\tx12427\tx13314\tx14202\tx15090\tx15977\tx16865\tx17753\tx18640\tx19528\tx20416\tx21303\tx22191\tx23079\tx23966\tx24854\tx25742\tx26629\tx27517\tx28405\tx29292\tx30180\tx31067\tx31955\tx32843\tx33730\tx34618\tx35506\tx36393\tx37281\tx38169\tx39056\tx39944\tx40832\tx41719\tx42607\tx43495\tx44382\tx45270\tx46158\tx47045\tx47933\tx48821\tx49708\tx50596\tx51484\tx52371\tx53259\tx54147\tx55034\tx55922\tx56810\tx57697\tx58585\tx59472\tx60360\tx61248\tx62135\tx63023\tx63911\tx64798\tx65686\tx66574\tx67461\tx68349\tx69237\tx70124\tx71012\tx71900\tx72787\tx73675\tx74563\tx75450\tx76338\tx77226\tx78113\tx79001\tx79889\tx80776\tx81664\tx82552\tx83439\tx84327\tx85215\tx86102\tx86990\tx87877\tx88765\slleading20\pardirnatural\partightenfactor0\f0\fs38 \cf2 Solitary zonal structures have recently been identified in gyrokinetic simulations of subcritical drift-wave (DW) turbulence with background shear flows. However, the nature of these structures has not been fully understood yet. Here, we show that similar structures can be obtained within a reduced model, which complements the modified Hasegawa\'97Mima equation with a generic primary instability and a background shear flow. We also find that these structures can be qualitatively reproduced in the modified Hasegawa\'97Wakatani equation, which subsumes the reduced model as a limit. In particular, we illustrate that in both cases, the solitary zonal structures approximately satisfy the same \'93equation of state\'94, which is a local relation connecting the DW envelope with the zonal-flow velocity. Due to this generality, our reduced model can be considered as a minimum model for solitary zonal structures in subcritical DWs. }
Type:
Dataset
Issue Date:
March 2020

110. Statistical properties of magnetic structures and energy dissipation during turbulent reconnection in the Earth's magnetotail

Author(s):
Bergstedt, K.; Ji, H.; Jara-Almonte, J.; Yoo, J.; Ergun, R. E.; Chen, L.-J.
Abstract:
We present the first statistical study of magnetic structures and associated energy dissipation observed during a single period of turbulent magnetic reconnection, by using the in situ measurements of the Magnetospheric Multiscale mission in the Earth's magnetotail on 26 July 2017. The structures are selected by identifying a bipolar signature in the magnetic field and categorized as plasmoids or current sheets via an automated algorithm which examines current density and plasma flow. The size of the plasmoids forms a decaying exponential distribution ranging from subelectron up to ion scales. The presence of substantial number of current sheets is consistent with a physical picture of dynamic production and merging of plasmoids during turbulent reconnection. The magnetic structures are locations of significant energy dissipation via electric field parallel to the local magnetic field, while dissipation via perpendicular electric field dominates outside of the structures. Significant energy also returns from particles to fields.
Type:
Dataset
Issue Date:
14 September 2020

111. Theory of the tertiary instability and the Dimits shift from reduced drift-wave models

Author(s):
Zhu, Hongxuan; Zhou, Yao; Dodin, I. Y.
Abstract:
Tertiary modes in electrostatic drift-wave turbulence are localized near extrema of the zonal velocity $U(x)$ with respect to the radial coordinate $x$. We argue that these modes can be described as quantum harmonic oscillators with complex frequencies, so their spectrum can be readily calculated. The corresponding growth rate $\gamma_{\rm TI}$ is derived within the modified Hasegawa--Wakatani model. We show that $\gamma_{\rm TI}$ equals the primary-instability growth rate plus a term that depends on the local $U''$; hence, the instability threshold is shifted compared to that in homogeneous turbulence. This provides a generic explanation of the well-known yet elusive Dimits shift, which we find explicitly in the Terry--Horton limit. Linearly unstable tertiary modes either saturate due to the evolution of the zonal density or generate radially propagating structures when the shear $|U'|$ is sufficiently weakened by viscosity. The Dimits regime ends when such structures are generated continuously.
Type:
Dataset
Issue Date:
January 2020

112. Theory of the tertiary instability and the Dimits shift within a scalar model

Author(s):
Zhu, Hongxuan; Zhou Yao; Dodin, I.Y.
Abstract:
The Dimits shift is the shift between the threshold of the drift-wave primary instability and the actual onset of turbulent transport in magnetized plasma. It is generally attributed to the suppression of turbulence by zonal flows, but developing a more detailed understanding calls for consideration of specific reduced models. The modified Terry--Horton system has been proposed by St-Onge [J. Plasma Phys. {\bf 83}, 905830504 (2017)] as a minimal model capturing the Dimits shift. Here, we use this model to develop an analytic theory of the Dimits shift and a related theory of the tertiary instability of zonal flows. We show that tertiary modes are localized near extrema of the zonal velocity $U(x)$, where $x$ is the radial coordinate. By approximating $U(x)$ with a parabola, we derive the tertiary-instability growth rate using two different methods and show that the tertiary instability is essentially the primary drift-wave instability modified by the local $U''$. Then, depending on $U''$, the tertiary instability can be suppressed or unleashed. The former corresponds to the case when zonal flows are strong enough to suppress turbulence (Dimits regime), while the latter corresponds to the case when zonal flows are unstable and turbulence develops. This understanding is different from the traditional paradigm that turbulence is controlled by the flow shear $U'$. Our analytic predictions are in agreement with direct numerical simulations of the modified Terry--Horton system.
Type:
Dataset
Issue Date:
June 2020

113. Toward fusion plasma scenario planning for NSTX-U using machine-learning-accelerated models

Author(s):
Mark D. Boyer
Abstract:
One of the most promising devices for realizing power production through nuclear fusion is the tokamak. To maximize performance, it is preferable that tokamak reactors achieve advanced operating scenarios characterized by good plasma confinement, improved magnetohydrodynamic (MHD) stability, and a largely non-inductively driven plasma current. Such scenarios could enable steady-state reactor operation with high \emph{fusion gain} --- the ratio of produced fusion power to the external power provided through the plasma boundary. Precise and robust control of the evolution of the plasma boundary shape as well as the spatial distribution of the plasma current, density, temperature, and rotation will be essential to achieving and maintaining such scenarios. The complexity of the evolution of tokamak plasmas, arising due to nonlinearities and coupling between various parameters, motivates the use of model-based control algorithms that can account for the system dynamics. In this work, a learning-based accelerated model trained on data from the National Spherical Torus Experiment Upgrade (NSTX-U) is employed to develop planning and control strategies for regulating the density and temperature profile evolution around desired trajectories. The proposed model combines empirical scaling laws developed across multiple devices with neural networks trained on empirical data from NSTX-U and a database of first-principles-based computationally intensive simulations. The reduced execution time of the accelerated model will enable practical application of optimization algorithms and reinforcement learning approaches for scenario planning and control development. An initial demonstration of applying optimization approaches to the learning-based model is presented, including a strategy for mitigating the effect of leaving the finite validity range of the accelerated model. The approach shows promise for actuator planning between experiments and in real-time.
Type:
Dataset
Issue Date:
May 2020

114. Vertical forces during VDEs in an ITER plasma and the role of halo currents

Author(s):
Clasuer, C; Jardin, S; Ferraro, N
Abstract:
Vertical displacement events (VDEs) can occur in elongated tokamaks causing large currents to flow in the vessel and other adjacent metallic structures. To better understand the potential magnitude of the associated forces and the role of the so called ``halo currents'' on them, we have used the M3D-C1 code to simulate potential VDEs in ITER. We used actual values for the vessel resistivity and pre-quench temperatures and, unlike most of the previous studies, the halo region is naturally formed by triggering the thermal quench with an increase in the plasma thermal conductivity. We used the 2D non-linear version of the code and vary the post-thermal quench thermal conductivity profile as well as the boundary temperature in order to generate a wide range of possible cases that could occur in the experiment. We also show that, for a similar condition, increasing the halo current does not increase the total force on the wall since it is offset by a decrease in the toroidal contribution.
Type:
Dataset
Issue Date:
February 2020

115. Visual Analogy Extrapolation Challenge (VAEC)

Author(s):
Webb, Taylor; Dulberg, Zachary; Frankland, Steven; Petrov, Alexander; O'Reilly, Randall; Cohen, Jonathan
Abstract:
Extrapolation -- the ability to make inferences that go beyond the scope of one's experiences -- is a hallmark of human intelligence. By contrast, the generalization exhibited by contemporary neural network algorithms is largely limited to interpolation between data points in their training corpora. In this paper, we consider the challenge of learning representations that support extrapolation. We introduce a novel visual analogy benchmark that allows the graded evaluation of extrapolation as a function of distance from the convex domain defined by the training data. We also introduce a simple technique, context normalization, that encourages representations that emphasize the relations between objects. We find that this technique enables a significant improvement in the ability to extrapolate, considerably outperforming a number of competitive techniques.
Type:
Dataset and Image
Issue Date:
2020

116. 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

117. 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

118. Blob wakes in NSTX

Author(s):
Zweben SJ, Myra JR, Diallo A, Russell DA, Scotti F, Stotler DP
Abstract:
Transient small-scale structures were identified in the wake of blobs movingpoloidally through the SOL of high-powered H-mode plasmas in NSTX, using the gaspuff imaging (GPI) diagnostic. These blob wakes had a poloidal wavelength in therange 3.5 cm, which is significantly smaller than the average blob scale of~12 cm, and the wakes had a poloidal velocity of 1.5 km/sec in theelectron diamagnetic direction, which is opposite to the blob poloidal velocity inthese shots. These wakes were radially localized 0-4 cm outside the separatrix andoccurred within ~50 microsec after the passage of a blob through the GPI field of view.The clearest wakes were seen when the GPI viewing angle was well aligned with thelocal B field line, as expected for such small-scale structures given the diagnosticgeometry. A plausible theoretical interpretation of the wakes is discussed: theobserved wakes share some features of drift waves and/or drift-Alfven waves whichcould be excited
Type:
Dataset
Issue Date:
July 2019

119. CKavity Library: Next-Generation Sequencing

Abstract:
Protein sequence space is vast; nature uses only an infinitesimal fraction of possible sequences to sustain life. Are there solutions to biological problems other than those provided by nature? Can we create artificial proteins that sustain life? To investigate this question, the Hecht lab has created combinatorial collections, or libraries, of novel sequences with no homology to those found in living organisms. These libraries were subjected to screens and selections, leading to the identification of sequences with roles in catalysis, modulating gene regulation, and metal homeostasis. However, the resulting functional proteins formed dynamic rather than well-ordered structures. This impeded structural characterization and made it difficult to ascertain a mechanism of action. To address this, Christina Karas's thesis work focuses on developing a new model of libraries based on the de novo protein S-824, a four-helix bundle with a very stable three-dimensional structure. The first part of this research focused on mutagenesis of S-824 and characterization of the resulting proteins, revealing that this scaffold tolerates amino acid substitutions, including buried polar residues and the removal of hydrophobic side chains to create a putative cavity. Distinct from previous libraries, Karas targeted variability to a specific region of the protein, seeking to create a cavity and potential active site. The second part of this work details the design and creation of a library encoding 1.7 x 10^6 unique proteins, assembled from degenerate oligonucleotides. The third and fourth parts of this work cover the screening effort for a range of activities, both in vitro and in vivo. I found that this collection binds heme readily, leading to abundant peroxidase activity. Hits for lipase and phosphatase activity were also detected. This work details the development of a new strategy for creating de novo sequences geared toward function rather than structure.
Type:
Dataset

120. Climate Impacts from Large Volcanic Eruptions in a High-resolution Climate Model: the Importance of Forcing Structure

Author(s):
Yang, Wenchang; Vecchi, Gabriel; Fueglistaler, Stephan; Horowitz, Larry; Luet, David; Muñoz, Ángel; Paynter, David; Underwood, Seth
Abstract:
Explosive volcanic eruptions have large climate impacts, and can serve as observable tests of the climatic response to radiative forcing. Using a high resolution climate model, we contrast the climate responses to Pinatubo, with symmetric forcing, and those to Santa Maria and Agung, which had meridionally asymmetric forcing. Although Pinatubo had larger global-mean forcing, asymmetric forcing strongly shifts the latitude of tropical rainfall features, leading to larger local precipitation/TC changes. For example, North Atlantic TC activity over is enhanced/reduced by SH-forcing (Agung)/NH-forcing (Santa Maria), but changes little in response to the Pinatubo forcing. Moreover, the transient climate sensitivity estimated from the response to Santa Maria is 20% larger than that from Pinatubo or Agung. This spread in climatic impacts of volcanoes needs to be considered when evaluating the role of volcanoes in global and regional climate, and serves to contextualize the well-observed response to Pinatubo.
Type:
Dataset
Issue Date:
2019

121. 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

122. Competitive learning modulates memory consolidation during sleep

Author(s):
Antony, James W.; Cheng, Larry Y.; Brooks, Paula P.; Paller, Ken A.; Norman, Kenneth A.
Abstract:
Competition between memories can cause weakening of those memories. Here we investigated memory competition during sleep in human participants by presenting auditory cues that had been linked to two distinct picture-location pairs during wake. We manipulated competition during learning by requiring participants to rehearse picture-location pairs associated with the same sound either competitively (choosing to rehearse one over the other, leading to greater competition) or separately; we hypothesized that greater competition during learning would lead to greater competition when memories were cued during sleep. With separate-pair learning, we found that cueing benefited spatial retention. With competitive-pair learning, no benefit of cueing was observed on retention, but cueing impaired retention of well-learned pairs (where we expected strong competition). During sleep, post-cue beta power (16–30 Hz) indexed competition and predicted forgetting, whereas sigma power (11–16 Hz) predicted subsequent retention. Taken together, these findings show that competition between memories during learning can modulate how they are consolidated during sleep.
Type:
Dataset
Issue Date:
November 2018

123. Design and simulation of the snowflake divertor control for NSTX-U

Author(s):
Vail, P. J.; Boyer, M. D.; Welander, A. S.; Kolemen, E.; U.S. Department of Energy contract number DE-AC02-09CH11466
Abstract:
This paper presents the development of a physics-based multiple-input-multiple-output algorithm for real-time feedback control of snowflake divertor (SFD) configurations on the National Spherical Torus eXperiment Upgrade (NSTX-U). A model of the SFD configuration response to applied voltages on the divertor control coils is first derived and then used, in conjunction with multivariable control synthesis techniques, to design an optimal state feedback controller for the configuration. To demonstrate the capabilities of the controller, a nonlinear simulator for axisymmetric shape control was developed for NSTX-U which simultaneously evolves the currents in poloidal field coils based upon a set of feedback-computed voltage commands, calculates the induced currents in passive conducting structures, and updates the plasma equilibrium by solving the free-boundary Grad-Shafranov problem. Closed-loop simulations demonstrate that the algorithm enables controlled operations in a variety of SFD configurations and provides capabilities for accurate tracking of time-dependent target trajectories for the divertor geometry. In particular, simulation results suggest that a time-varying controller which can properly account for the evolving SFD dynamical response is not only desirable but necessary for achieving acceptable control performance. The algorithm presented in this paper has been implemented in the NSTX-U Plasma Control System in preparation for future control and divertor physics experiments.
Type:
Dataset
Issue Date:
April 2019

124. Electron inertial effects on linearly polarized electromagnetic ion cyclotron waves at Earth's magnetosphere

Author(s):
Kim, Eun-Hwa; Johnson, Jay; Lee, Dong-Hun
Abstract:
We discuss a role of the electron inertial effect on linearly polarized electromagnetic ion cyclotron (EMIC) waves at Earth. The linearly polarized EMIC waves have been previously suggested to be generated via mode conversion from the fast compressional wave at the ion-ion hybrid (IIH) resonance. When the electron inertial effects are neglected, the wave normal angle of the mode-converted IIH waves is 90 degrees because the wavevector perpendicular to the magnetic field becomes infinite at the IIH resonance. When the electron inertial effect is considered, the mode-converted IIH waves can propagate across the magnetic field lines and the wavelength perpendicular to the magnetic field approaches the electron inertial length scale near the Buchsbaum resonance. These waves are referred to as electron inertial waves. Due to the electron inertial effect, the perpendicular wavenumber to the ambient magnetic field near the IIH resonance remains finite and the wave normal angle is less than 90 degrees. The wave normal angle where the maximum absorption occurs in a dipole magnetic field is 30-80 degrees, which is consistent with the observed values near the magnetic equator. Therefore, the numerical results suggest that the linearly polarized EMIC wave generated via mode conversion near the IIH resonance can be detected in between the Buchsbaum and the IIH resonance frequencies, and these waves can have normal angle less than 90 degrees.
Type:
Dataset
Issue Date:
April 2019

125. Elemental and topographical imaging of microscopic variations in deposition on NSTX-U and DIII-D samples2

Author(s):
C.H. Skinner, C.P. Chrobak, R. Kaita, B.E.Koel
Abstract:
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.
Type:
Dataset
Issue Date:
April 2019

126. Formation of solitary zonal structures via the modulational instability of drift waves

Author(s):
Zhou, Yao; Zhu, Hongxuan; Dodin, I. Y.
Abstract:
{\rtf1\ansi\ansicpg1252\cocoartf1561\cocoasubrtf600 {\fonttbl\f0\fswiss\fcharset0 Helvetica;} {\colortbl;\red255\green255\blue255;\red0\green0\blue0;} {\*\expandedcolortbl;;\cssrgb\c0\c0\c0;} \margl1440\margr1440\vieww10800\viewh8400\viewkind0 \pard\tx887\tx1775\tx2662\tx3550\tx4438\tx5325\tx6213\tx7101\tx7988\tx8876\tx9764\tx10651\tx11539\tx12427\tx13314\tx14202\tx15090\tx15977\tx16865\tx17753\tx18640\tx19528\tx20416\tx21303\tx22191\tx23079\tx23966\tx24854\tx25742\tx26629\tx27517\tx28405\tx29292\tx30180\tx31067\tx31955\tx32843\tx33730\tx34618\tx35506\tx36393\tx37281\tx38169\tx39056\tx39944\tx40832\tx41719\tx42607\tx43495\tx44382\tx45270\tx46158\tx47045\tx47933\tx48821\tx49708\tx50596\tx51484\tx52371\tx53259\tx54147\tx55034\tx55922\tx56810\tx57697\tx58585\tx59472\tx60360\tx61248\tx62135\tx63023\tx63911\tx64798\tx65686\tx66574\tx67461\tx68349\tx69237\tx70124\tx71012\tx71900\tx72787\tx73675\tx74563\tx75450\tx76338\tx77226\tx78113\tx79001\tx79889\tx80776\tx81664\tx82552\tx83439\tx84327\tx85215\tx86102\tx86990\tx87877\tx88765\slleading20\pardirnatural\partightenfactor0 \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.\ }
Type:
Dataset
Issue Date:
June 2019

127. Global modeling of wall material migration following boronization in NSTX-U

Author(s):
Nichols, J.H.; Jaworski, M.A.; Skinner, C.H.; Bedoya, F.; Scotti, F.; Soukhanovskii, V.A.; Schmid, K.
Abstract:
Boronization is commonly utilized in tokamaks to suppress intrinsic impurities, most notably oxygen from residual water vapor. However, this is a temporary solution, as oxygen levels typically return to pre-boronization levels following repeated plasma exposure. The global impurity migration model WallDYN has been applied to the post-boronization surface impurity evolution in NSTX-U. A “Thin Film Model” has been incorporated into WallDYN to handle spatially inhomogeneous conditioning films of varying thicknesses, together with an empirical boron conditioning model for the NSTX-U glow discharge boronization process. The model qualitatively reproduces the spatial distribution of boron in the NSTX-U vessel, the spatially-resolved divertor emission pattern, and the increase in oxygen levels following boronization. The simulations suggest that oxygen is primarily sourced from wall locations without heavy plasma flux or significant boron deposition, namely the lower and upper passive plates and the lower private flux zone.
Type:
Dataset
Issue Date:
March 2019

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

Author(s):
Guttenfelder, W.; Kaye, S.M.; Kreite, D.M.; Bell, R.E.; Diallo, A.; LeBlanc, B.P.; McKee, G.R.; Podesta, M.; Sabbagh, S.A.; Smith, D.R.
Abstract:
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 $\rho$=0.6-0.8, although $E\timesB$ shearing rates are larger than the linear growth rates over much of that region. Deeper in the core ($\rho$=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 $\rho$=0.76, whereas electromagnetic MTM transport is important at $\rho$=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. However, in the region between ($\rho$=0.56, 0.66), ion-scale simulations are strongly suppressed by the locally large $E\timesB$ 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 $\rho_*$=$\rho_i$/a at low aspect ratio), including electromagnetic and possibly multi-scale effects, is also discussed.
Type:
Dataset
Issue Date:
March 2019

129. Intrinsic Rotation in Axisymmetric Devices

Author(s):
T Stoltzfus-Dueck
Abstract:
Toroidal rotation is critical for fusion in tokamaks, since it stabilizes instabilities that can otherwise cause disruptions or degrade confinement. Unlike present-day devices, ITER might not have enough neutral-beam torque to easily avoid these instabilities. We must therefore understand how the plasma rotates intrinsically, that is, without applied torque. Experimentally, torque-free plasmas indeed rotate, with profiles that are often non-flat and even non-monotonic. The rotation depends on many plasma parameters including collisionality and plasma current, and exhibits sudden bifurcations (rotation reversals) at critical parameter values.Since toroidal angular momentum is conserved in axisymmetric systems, and since experimentally inferred momentum transport is much too large to be neoclassical, theoretical work has focused on rotation drive by nondiffusive turbulent momentum fluxes. In the edge, intrinsic rotation relaxes to a steady state in which the total momentum outflux from the plasma vanishes. Ion drift orbits, scrape-off-layer flows, separatrix geometry, and turbulence intensity gradient all play a role. In the core, nondiffusive and viscous momentum fluxes balance to set the rotation gradient at each flux surface. Although many mechanisms have been proposed for the nondiffusive fluxes, most are treated in one of two distinct but related gyrokinetic formulations. In a radially local fluxtube, appropriate for rho star <<1, the lowest-order gyrokinetic formulations exhibit a symmetry that prohibits nondiffusive momentum flux for nonrotating plasmas in an up- down symmetric magnetic geometry with no ExB shear. Many symmetry-breaking mechanisms have been identified, but none have yet been conclusively demonstrated to drive a strong enough flux to explain commonly observed experimental rotation profiles. Radially global gyrokinetic simulations naturally include many symmetry-breaking mechanisms, and have shown cases with experimentally relevant levels of nondiffusive flux. These promising early results motivate further work to analyze, verify, and validate.This article provides a pedagogical introduction to intrinsic rotation in axisymmetric devices. Intended for both newcomers to the topic and experienced practitioners, the article reviews a broad range of topics including experimental and theoretical results for both edge and core rotation, while maintaining a focus on the underlying concepts.
Type:
Dataset
Issue Date:
November 2019

130. Machine Learning Characterization of Alfvénic and Sub-Alfvénic Chirping and Correlation With Fast-Ion Loss at NSTX

Author(s):
Woods, B. J. Q.; Duarte, V. N.; Fredrickson, E. D.; Gorelenkov, N. N.; Podestà, M.; Vann, R. G. L.
Abstract:
Abrupt large events in the Alfvenic and sub-Alfvenic frequency bands in tokamaks are typically correlated with increased fast-ion loss. Here, machine learning is used to speed up the laborious process of characterizing the behavior of magnetic perturbations from corresponding frequency spectrograms that are typically identified by humans. The analysis allows for comparison between different mode character (such as quiescent, fixed frequency, and chirping, avalanching) and plasma parameters obtained from the TRANSP code, such as the ratio of the neutral beam injection (NBI) velocity and the Alfven velocity (v_inj./v_A), the q-profile, and the ratio of the neutral beam beta and the total plasma beta (beta_beam,i / beta). In agreement with the previous work by Fredrickson et al., we find a correlation between beta_beam,i and mode character. In addition, previously unknown correlations are found between moments of the spectrograms and mode character. Character transition from quiescent to nonquiescent behavior for magnetic fluctuations in the 50200-kHz frequency band is observed along the boundary v_phi ~ (1/4)(v_inj. - 3v_A), where v_phi is the rotation velocity.
Type:
Dataset
Issue Date:
December 2019

131. Modelling of Ablatant Deposition from Electromagnetically Driven Radiative Pellets for Disruption Mitigation Studies

Author(s):
Lunsford, Robert; Raman, Roger; Brooks, Arthur; Ellis, Robert A.; Lay, W-S;
Abstract:
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
Type:
Dataset
Issue Date:
June 2019

132. Nonlinear saturation and oscillations of collisionless zonal flows

Author(s):
Zhu, Hongxuan; Zhou, Yao; Dodin, I. Y.
Abstract:
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.
Type:
Dataset
Issue Date:
May 2019

133. Physics design for a lithium vapor box divertor experiment on Magnum-PSI

Author(s):
Schwartz, Jacob; Emdee, Eric; Goldston, Robert; Jaworski, Michael
Abstract:
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.
Type:
Dataset
Issue Date:
30 January 2019

134. Quasioptical modeling of wave beams with and without mode conversion: II. Numerical simulations of single-mode beams

Author(s):
K. Yanagihara, I. Y. Dodin, and S. Kubo
Abstract:
This work continues a series of papers where we propose an algorithm for quasioptical modeling of electromagnetic beams with and without mode conversion. The general theory was reported in the first paper of this series, where a parabolic partial differential equation was derived for the field envelope that may contain one or multiple modes with close group velocities. Here, we present a corresponding code PARADE (PAraxial RAy DEscription) and its test applications to single-mode beams in vacuum and also in inhomogeneous magnetized plasma. The numerical results are compared, respectively, with analytic formulas from Gaussian-beam optics and also with cold-plasma ray tracing. Quasioptical simulations of mode-converting beams are reported in the next, third paper of this series.
Type:
Dataset
Issue Date:
August 2019

135. Quasioptical modeling of wave beams with and without mode conversion: III. Numerical simulations of mode-converting beams

Author(s):
K. Yanagihara, I. Y. Dodin, and S. Kubo
Abstract:
This work continues a series of papers where we propose an algorithm for quasioptical modeling of electromagnetic beams with and without mode conversion. The general theory was reported in the first paper of this series, where a parabolic partial differential equation was derived for the field envelope that may contain one or multiple modes with close group velocities. In the second paper, we presented a corresponding code PARADE (PAraxial RAy DEscription) and its test applications to single-mode beams. Here, we report quasioptical simulations of mode-converting beams for the first time. We also demonstrate that PARADE can model splitting of two-mode beams. The numerical results produced by PARADE show good agreement with those of one-dimensional full-wave simulations and also with conventional ray tracing (to the extent that one-dimensional and ray-tracing simulations are applicable).
Type:
Dataset
Issue Date:
August 2019

136. Real-time capable modeling of neutral beam injection on NSTX-U using neural networks

Author(s):
Boyer, M.D.; Kaye, S.; Erickson, K.
Abstract:
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.
Type:
Dataset
Issue Date:
February 2019

137. Reduced Model for Direct Induction Startup Scenario Development on MAST-U and NSTX-U

Author(s):
Battaglia, D.J.; Thornton, A.J.; Gerhardt, S.P.; Kirk, A.; Kogan, L; Menard, J.E.
Abstract:
A reduced semi-empirical model using time-dependent axisymmetric vacuum field calculations is used to develop the prefill and feed-forward coil current targets required for reliable direct induction (DI) startup on the new MA-class spherical tokamaks, MAST-U and NSTX-U. The calculations are constrained by operational limits unique to each device, such as the geometry of the conductive elements and active coils, power supply specifications and coil heating and stress limits. The calculations are also constrained by semi-empirical models for sufficient breakdown, current drive, equilibrium and stability of the plasma developed from a shared database. A large database of DI startup on NSTX and NSTX-U is leveraged to quantify the requirements for achieving a reliable breakdown (Ip ~ 20 kA). It is observed that without pre-ionization, STs access the large E/P regime at modest loop voltage (Vloop) where the electrons in the weakly ionized plasma are continually accelerating along the open field lines. This ensures a rapid (order millisecond) breakdown of the neutral gas, even without pre-ionization or high-quality field nulls. The timescale of the initial increase in Ip on NSTX is reproduced in the reduced model provided a mechanism for impeding the applied electric field is included. Most discharges that fail in the startup phase are due to an inconsistency in the evolution of the plasma current (Ip) and equilibrium field or loss of vertical stability during the burn-through phase. The requirements for the self-consistent evolution of the fields in the weakly and full-ionized plasma states are derived from demonstrated DI startup on NSTX, NSTX-U and MAST. The predictive calculations completed for MAST-U and NSTX-U illustrate that the maximum Ip ramp rate (dIp/dt) in the early startup phase is limited by the voltage limits on the poloidal field coils on MAST-U and passive vertical stability on NSTX-U.
Type:
Dataset
Issue Date:
August 2019

138. Sleep spindle refractoriness segregates periods of memory reactivation

Author(s):
Antony, James W.; Piloto, Luis; Wang, Margaret; Brooks, Paula P.; Norman, Kenneth A.; Paller, Ken A.
Abstract:
The stability of long-term memories is enhanced by reactivation during sleep. Correlative evidence has linked memory reactivation with thalamocortical sleep spindles, although their functional role is not fully understood. Our initial study replicated this correlation and also demonstrated a novel rhythmicity to spindles, such that a spindle is more likely to occur approximately 3–6 s following a prior spindle. We leveraged this rhythmicity to test the role of spindles in memory by using real-time spindle tracking to present cues within versus just after the presumptive refractory period; as predicted, cues presented just after the refractory period led to better memory. Our findings demonstrate a precise temporal link between sleep spindles and memory reactivation. Moreover, they reveal a previously undescribed neural mechanism whereby spindles may segment sleep into two distinct substates: prime opportunities for reactivation and gaps that segregate reactivation events.
Type:
Dataset
Issue Date:
4 June 2018

139. Sound velocities in shock-synthesized stishovite to 72 GPa

Author(s):
Berryman, Eleanor J.; Winey, J. M.; Gupta, Yogendra M.; Duffy, Thomas S.
Abstract:
Stishovite (rutile-type SiO2) is the archetype of dense silicates and may occur in post-garnet eclogitic rocks at lower-mantle conditions. Sound velocities in stishovite are fundamental to understanding its mechanical and thermodynamic behavior at high pressure and temperature. Here, we use plate-impact experiments combined with velocity interferometry to determine the stress, density, and longitudinal sound speed in stishovite formed during shock compression of fused silica at 44 GPa and above. The measured sound speeds range from 12.3(8) km/s at 43.8(8) GPa to 9.8(4) km/s at 72.7(11) GPa. The decrease observed at 64 GPa reacts a decrease in the shear modulus of stishovite, likely due to the onset of melting. By 72 GPa, the measured sound speed agrees with the theoretical bulk sound speed indicating loss of all shear stiffness due to complete melting. Our sound velocity results provide direct evidence for shock-induced melting, in agreement with previous pyrometry data.
Type:
Dataset
Issue Date:
2019

140. Spontaneous multi-keV electron generation in a low-RF-power axisymmetric mirror machine

Author(s):
Swanson, C.;Cohen, S.A.
Abstract:
Title: Spontaneous multi-keV electron generation in a low-RF-power axisymmetric mirror machine Abstract: X-ray emission shows the existence of multi-keV electrons in low-temperature, low-power, capacitively-coupled RF-heated magnetic-mirror plasmas that also contain a warm (300 eV) minority electron population. Though these warm electrons are initially passing particles, we suggest that collisionless scattering -- mu non-conservation in the static vacuum field -- is responsible for a minority of them to persist in the mirror cell for thousands of transits during which time a fraction are energized to a characteristic temperature of 3 keV, with some electrons reaching energies above 30 keV. A heuristic model of the heating by a Fermi-acceleration-like mechanism is presented, with mu non-conservation in the static vacuum field as an essential feature.
Type:
Dataset
Issue Date:
May 2019

141. Study of liquid metal surface wave damping in the presence of magnetic fields and electrical currents

Author(s):
Fisher, A.E.; Hvasta, M.G.; Kolemen, E.
Abstract:
Experiments and predictions of surface wave damping in liquid metal due to a surface aligned magnetic field and externally regulated j × B force are presented. Fast-flowing, liquid-metal plasma facing components (LM-PFCs) are a proposed alternative to solid PFCs that are unable to handle the high heat flux, thermal stresses, and radiation damage in a tokamak. The significant technical challenges associated with LM-PFCs compared to solid PFCs are justified by greater heat flux management, self-healing properties, and reduced particle recycling. However, undesirable engineering challenges such as evaporation and splashing of the liquid metal introduce excessive impurities into the plasma and degrade plasma performance. Evaporation may be avoided through high-speed flow that limits temperature rise of the liquid metal by reducing heat flux exposure time, but as flow speed increases the surface may become more turbulent and prone to splashing and uneven surfaces. Wave damping is one mechanism that reduces surface disturbance and thus the chances of liquid metal impurity introduction into the plasma. Experiments on the Liquid Metal eXperiment Upgrade (LMX-U) examined damping under the influence of transverse magnetic fields and vertically directed Lorentz force.
Type:
Dataset
Issue Date:
March 2019

142. Study of the impact of pre- and real-time deposition of lithium on plasma performance on NSTX

Author(s):
Canal, G.P.; Maingi, R.; Evans, T.E.; Kaye, S.M.; Mansfield, D.K.
Abstract:
The efficiency of two lithium (Li) injection methods used on the National Spherical Torus Experiment (NSTX) are compared in terms of the amount of Li used to produce equivalent plasma performance improvements, namely Li evaporation over the divertor plates, prior to the initiation of the discharge, and real-time Li injection directly into the plasma scrape-off layer during the discharge. The measurements show that the real-time method can affect the energy confinement and edge stability of NSTX plasmas in a more efficient way than the Li evaporation method as it requires only a fraction of the amount of Li used by the evaporation method to produce similar improvements.
Type:
Dataset
Issue Date:
January 2019

143. TRANSP-based Optimization Towards Tokamak Scenario Development

Author(s):
W.P. Wehner, E. Schuster, M.D. Boyer, F. Poli
Abstract:
An optimization approach that incorporates the predictive transport code TRANSP is proposed for tokamak scenario development. Optimization methods are often employed to develop open-loop control strategies to aid access to high performance tokamak scenarios. In general, the optimization approaches use control-oriented models, i.e. models that are reduced in complexity and prediction accuracy as compared to physics-oriented transport codes such as TRANSP. In the presented approach, an optimization procedure using the TRANSP code to simulate the tokamak plasma is considered for improved predictive capabilities. As a test case, the neutral beam injection (NBI) power is optimized to develop a control strategy that maximizes the non-inductive current fraction during the ramp-up phase for NSTX-U. Simulation studies towards the achievement of non-inductive ramp up in NSTX-U have already been carried out with the TRANSP code. The optimization-based approach proposed in this work is used to maximize the non-inductive current fraction during ramp-up in NSTX-U, demonstrating that the scenario development task can be automated. An additional test case considers optimization of the current ramp rate in DIII-D for obtaining a stationary plasma characterized by a flat loop voltage profile in the flattop phase.
Type:
Dataset
Issue Date:
April 2019

144. TRANSP-based closed-loop simulations of current profile optimal regulation in NSTX-Upgrade

Author(s):
Ilhan, Z.O.; Boyer. M.D.; Schuster, E.
Abstract:
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.
Type:
Dataset
Issue Date:
March 2019

146. Verification of the global gyrokinetic stellarator code XGC-S for linear ion temperature gradient driven modes

Author(s):
Cole M; Hager R; Moritaka T; Dominski J; Kleiber R; Ku S; Lazerson S; Riemann J; Chang C
Abstract:
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.
Type:
Dataset
Issue Date:
August 2019

147. 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

148. 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

149. 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

150. Demonstrating electromagnetic control of free-surface, liquid-metal flows relevant to fusion reactors

Author(s):
Hvasta, M. G.; Kolemen, E.; Fisher, A. E.; Ji, H.
Abstract:
Plasma-facing components (PFC's) made from solid materials may not be able to withstand the large heat and particle fluxes that will be produced within next-generation fusion reactors. To address the shortcomings of solid PFC's, a variety of liquid-metal (LM) PFC concepts have been proposed. Many of the suggested LM-PFC designs rely on electromagnetic restraint (Lorentz force) to keep free-surface, liquid-metal flows adhered to the interior surfaces of a fusion reactor. However, there is very little, if any, experimental data demonstrating that free-surface, LM-PFC's can actually be electromagnetically controlled. Therefore, in this study, electrical currents were injected into a free-surface liquid-metal that was flowing through a uniform magnetic field. The resultant Lorentz force generated within the liquid-metal affected the velocity and depth of the flow in a controllable manner that closely matched theoretical predictions. These results show the promise of electromagnetic control for LM-PFC's and suggest that electromagnetic control could be further developed to adjust liquid-metal nozzle output, prevent splashing within a tokamak, and alter heat transfer properties for a wide-range of liquid-metal systems.
Type:
Dataset
Issue Date:
January 2018

151. Design and measurement methods for a lithium vapor box similarity experiment

Author(s):
Schwartz, J. A.; Emdee, E. D.; Jaworski, M. A; Goldston, R. J.
Abstract:
The lithium vapor box divertor is a concept for handling the extreme divertor heat fluxes in magnetic fusion devices. In a baffled slot divertor, plasma interacts with a dense cloud of Li vapor which radiates and cools the plasma, leading to recombination and detachment. Before testing on a tokamak the concept should be validated: we plan to study detachment and heat redistribution by a Li vapor cloud in laboratory experiments. Mass changes and temperatures are measured to validate a Direct Simulation Monte Carlo model of neutral Li. The initial experiment involves a 5 cm diameter steel box containing 10g of Li held at 650 degrees C as vapor flows out a wide nozzle into a similarly-sized box at a lower temperature. Diagnosis is made challenging by the required material compatibility with lithium vapor. Vapor pressure is a steep function of temperature, so to validate mass flow models to within 10%, absolute temperature to within 4.5K is required. The apparatus is designed to be used with an analytical balance to determine mass transport. Details of the apparatus and methods of temperature and mass flow measurements are presented.
Type:
Dataset
Issue Date:
August 2018

152. Detection of an electron beam in a high density plasma via an electrostatic probe

Author(s):
Majeski, Stephen; Yoo, Jongsoo; Zweben, Stewart; Yamada, Masaaki
Abstract:
An electron beam is detected by a 1D floating potential probe array in a relatively high density (10e12 − 10e13 cm−3) and low temperature (∼ 5 eV) plasma of the Magnetic Reconnection Experiment (MRX). Clear perturbations in the floating potential profile by the electron beam are observed. Based on the floating potential profile and a current balance equation to the probe array tips, the effective width of the electron beam is determined, from which we determine the radial and toroidal beam current density profiles. After the profile of the electron beam is specified from the measured beam current, we demonstrate the consistency of the current balance equation and the location of the perturbation is also in agreement with field line mapping. No significant broadening of the electron beam is observed after the beam propagates for tens of centimeters through the high density plasma. These results prove that the field line mapping is, in principle, possible in high density plasmas.
Type:
Dataset
Issue Date:
2018

153. ELM frequency enhancement and discharge modification through lithium granule injection into EAST H-modes

Author(s):
Lunsford; Hsu, J.S.; Sun, Z.; Maingi, R.; Mansfield, D.K.; Xu, W.; Zuo, G.Z.; Huang, M.; Diallo, A.; Osborne, T.; Tritz, K.; Canik, J.; Meng, X.C.; Gong, X.Z.; Wan, B.N.; Li, J.G.; EAST Team
Abstract:
The injection of impurity granules into fusion research discharges can serve as a catalyst for ELM events. For sufficiently low ELM frequencies, and granule sizes above a threshold, this can result in full control of the ELM cycle, referred to as ELM pacing. For this research, we extend the investigation to conditions where the natural ELM frequency is too high for ELM pacing to be realized. Utilizing multiple sizes of lithium granules and classifying their effects by granule size, we demonstrate that ELM mitigation through frequency multiplication can be used at ELM triggering rates that nominally make ELM pacing unrealizable. We find that above a size threshold, injected granules promptly trigger ELMs and commensurately enhance the ELM frequency . Below this threshold size, injection of an individual granule does not always lead to the prompt triggering of an ELM; however, collective ablation in the edge pedestal region does enhance the ELM frequency. Specifically, Li granules too small to individually trigger ELMs were injected into EAST H-mode discharges at frequencies up to 2.3 kHz; collectively the granules were observed to enhance the natural ELM frequency up to 620 Hz, resulting in a ~2.4x multiplication of the natural ELM frequency and a 50% decrease of the ELM size.
Type:
Dataset
Issue Date:
October 2018

154. ELM-free and inter-ELM divertor heat flux broadening induced by Edge Harmonics Oscillation in NSTX

Author(s):
Gan, K.; Ahn, J.-W.; Gray, T.K.; Zweben, S.J.; Fredrickson, E.D.; Scotti, F.; Maingi, R.; Park, J.-K.; Canal, G.P.; Soukhanovskii, V.A.; McLean, A.G.; Wirth, B.D.
Abstract:
A new n=1 dominated Edge Harmonic Oscillation (EHO) has been found in NSTX. The new EHO, rotating toroidally in the counter-current direction and the opposite direction of the neutral beam, was observed during certain inter-ELM and ELM-free periods of H-mode operation. This EHO is associated with a significant broadening of the integral heat flux width (?_int) by up to 150%, and a decrease in the divertor peak heat flux by >60%. An EHO induced filament was also observed by the gas puff imaging diagnostic. The toroidal rotating filaments could change the edge magnetic topology resulting in toroidal rotating strike point splitting and heat flux broadening. Experimental result of the counter current rotation of strike points splitting is consistent with the counter-current EHO.
Type:
Dataset
Issue Date:
December 2017

155. Effects of Axial Boundary Conductivity on a Free Stewartson-Shercliff Layer

Author(s):
Caspary, Kyle J.; Choi, Dahan; Ebrahimi, Fatima; Gilson, Erik P.; Goodman, Jeremy; Ji, Hantao
Abstract:
The effects of axial boundary conductivity on the formation and stability of a magnetized free Stewartson-Shercliff layer (SSL) in a short Taylor-Couette device are reported. As the axial field increases with insulating endcaps, hydrodynamic Kelvin-Helmholtz-type instabilities set in at the SSLs of the conducting fluid, resulting in a much reduced flow shear. With conducting endcaps, SSLs respond to an axial field weaker by the square root of the conductivity ratio of endcaps to fluid. Flow shear continuously builds up as the axial field increases despite the local violation of the Rayleigh criterion, leading to a large number of hydrodynamically unstable modes. Numerical simulations of both the mean flow and the instabilities are in agreement with the experimental results.
Type:
Dataset
Issue Date:
2018

156. Electromagnetic Particle Injector for Fast Time Response Disruption Mitigation in Tokamaks

Author(s):
Raman, R.; Lay, W.-S.; Jarboe, T.R.; Menard, J.E.; Ono, M.
Abstract:
A novel, rapid time-response, disruption mitigation system referred to as the Electromagnetic Particle Injector (EPI) is described. This method can accurately deliver the radiative payload to the plasma center on a <10 ms time scale, much faster, and deeper, than what can be achieved using conventional methods. The EPI system accelerates a sabot electromagnetically. The sabot is a metallic capsule that can be accelerated to desired velocities by an electromagnetic impeller. At the end of its acceleration, within 2 ms, the sabot will release a radiative payload, which is composed of low-z granules, or a shell pellet containing smaller pellets. The primary advantage of the EPI concept over gas propelled systems is its potential to meet short warning time scales, while accurately delivering the required particle size and materials at the velocities needed for achieving the required penetration depth in high power ITER-scale discharges for thermal and runaway current disruption mitigation. The present experimental tests from a prototype system have demonstrated the acceleration of a 3.2 g sabot to over 150 m/s within 1.5 ms, consistent with the calculations, giving some degree of confidence that larger ITER-scale injector can be developed.
Type:
Dataset
Issue Date:
November 2018

157. Electron-density-sensitive Line Ratios of Fe XIII--XVI from Laboratory Sources Compared to CHIANTI

Author(s):
Weller, M.E.; Beiersdorfer, P.; Soukhanovskii, V.A.; Scotti, F.; LeBlanc, B.P.
Abstract:
We present electron-density-sensitive line ratios for Fe XIII–XVI measured in the spectral wavelength range of 200–440 Å and an electron density range of (1–4)×10^13 cm^−3. The results provide a test at the high-density limit of density-sensitive line ratios useful for astrophysical studies. The measurements were performed on the National Spherical Torus Experiment-Upgrade, where electron densities were measured independently by the laser Thomson scattering diagnostic. Spectra were collected with a flat-field grazing-incidence spectrometer, which provided a spectral resolution of up to 0.3 Å, i.e., high resolution across the broad wavelength range. The response of the instrument was relatively calibrated using spectroscopic techniques in order to improve accuracy. The line ratios are compared to other laboratory sources and the latest version of CHIANTI (8.0.2), and an agreement within 30% is found.
Type:
Dataset
Issue Date:
February 2018

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

Author(s):
Skinner, C.H.; Chrobak, C.P.; Kaita, R.; Koel, B.E.
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.
Type:
Dataset
Issue Date:
October 2018

159. Experimental demonstration of hydraulic jump control in liquid metal channel flow using Lorentz force

Author(s):
Fisher, Adam; Kolemen, Egemen; Hvasta, Mike
Abstract:
In this paper, hydraulic jump control using electromagnetic force in a liquid metal flow is presented. The control methods used give insight into the hydraulic jump behavior in the presence of magnetic fields and electrical currents. Flowing liquid metals is a proposed solution to heat flux challenges posed in fusion reactors, specifically the tokamak. Unfortunately, thin, fast-flowing liquid metal divertor concepts for fusion reactors are susceptible to hydraulic jumps that drastically reduce the liquid metal flow speed, leading to potential problems such as excessive evaporation, unsteady power removal, and possible plasma disruption. Highly electrically conductive flows within the magnetic fields do not exhibit traditional hydraulic jump behavior. There is very little research investigating the use of externally injected electrical currents and magnetic fields to control liquid metal hydraulic jumps. By using externally injected electrical currents and a magnetic field, a Lorentz force (also referred to as j × B force) may be generated to control the liquid metal jump behavior. In this work, a free-surface liquid metal—GaInSn eutectic or “galinstan”—flow through an electrically insulating rectangular duct was investigated. It was shown that applying a Lorentz force has a repeatable and predictable impact on the hydraulic jump, which can be used for liquid metal control within next-generation fusion reactors.
Type:
Dataset
Issue Date:
2018

160. Fast animal pose estimation using deep neural networks

Author(s):
Pereira, Talmo D.; Aldarondo, Diego E.; Willmore, Lindsay; Kislin, Mikhail; Wang, Samuel S.-H.; Murthy, Mala; Shaevitz, Joshua W.
Abstract:
Recent work quantifying postural dynamics has attempted to define the repertoire of behaviors performed by an animal. However, a major drawback to these techniques has been their reliance on dimensionality reduction of images which destroys information about which parts of the body are used in each behavior. To address this issue, we introduce a deep learning-based method for pose estimation, LEAP (LEAP Estimates Animal Pose). LEAP automatically predicts the positions of animal body parts using a deep convolutional neural network with as little as 10 frames of labeled data for training. This framework consists of a graphical interface for interactive labeling of body parts and software for training the network and fast prediction on new data (1 hr to train, 185 Hz predictions). We validate LEAP using videos of freely behaving fruit flies (Drosophila melanogaster) and track 32 distinct points on the body to fully describe the pose of the head, body, wings, and legs with an error rate of <3% of the animal's body length. We recapitulate a number of reported findings on insect gait dynamics and show LEAP's applicability as the first step in unsupervised behavioral classification. Finally, we extend the method to more challenging imaging situations (pairs of flies moving on a mesh-like background) and movies from freely moving mice (Mus musculus) where we track the full conformation of the head, body, and limbs.
Type:
Dataset
Issue Date:
30 May 2018

162. Modeling of a Laser-Induced Rydberg Spectroscopy diagnostic for Direct Measurement of the Local Electric Field in the Edge Region of NSTX/NSTX-U

Author(s):
Reymond, L.; Diallo, A.; Vekselman, V.
Abstract:
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}.
Type:
Dataset
Issue Date:
July 2018

163. On the scattering correction of fast-ion D-alpha signal on NSTX-U

Author(s):
Hao, G.Z; Heidbrink, W.W.; Liu, D.; Stagner, L.; Podesta, M.; Bortolon, A.
Abstract:
Analysis of fast-ion D-alpha (FIDA) data on National Spherical Torus Experiment-Upgrade (NSTX-U) shows that the cold Dα line contaminates the FIDA baseline. The scattered light is comparable to the FIDA emission. A scattering correction is required to extract the FIDA signal. Two methods that relate the scattered light contamination to the intensity of the cold Dα line are employed. One method uses laboratory measurements with a calibration lamp; the other method uses data acquired during plasma operation and singular value decomposition analysis. After correction, both the FIDA spectra and spatial profile are in better agreement with theoretical predictions.
Type:
Dataset
Issue Date:
June 2018

164. Plasma boundary shape control and real-time equilibrium reconstruction on NSTX-U

Author(s):
Boyer, M.; Battaglia, D.; Mueller, D.; Eidietis, N.; Erickson, K.; Ferron, J.; Gates, D.; Gerhardt, S.; Johnson, R.; Kolemen, E.; Menard, J.; Myers, C.; Sabbagh, S.; Scotti, F.; Vail, P.
Abstract:
The upgrade to the National Spherical Torus eXperiment (NSTX-U) included two main improvements: a larger center-stack, enabling higher toroidal field and longer pulse duration, and the addition of three new tangentially aimed neutral beam sources, which increase available heating and current drive, and allow for flexibility in shaping power, torque, current, and particle deposition profiles. To best use these new capabilities and meet the high-performance operational goals of NSTX-U, major upgrades to the NSTX-U Control System (NCS) hardware and software have been made. Several control algorithms, including those used for real-time equilibrium reconstruction and shape control, have been upgraded to improve and extend plasma control capabilities. As part of the commissioning phase of first plasma operations, the shape control system was tuned to control the boundary in both inner-wall limited and diverted discharges. It has been used to accurately track the requested evolution of the boundary (including the size of the inner gap between the plasma and central solenoid, which is a challenge for the ST configuration), X-point locations, and strike point locations, enabling repeatable discharge evolutions for scenario development and diagnostic commissioning.
Type:
Dataset
Issue Date:
March 2018

165. Quantitative imaging of carbon dimer precursor for nanomaterial synthesis in the carbon arc

Author(s):
Vekselman, V.; Khrabry, A.; Kaganovich, I.; Stratton, B.; Selinsky, R. S.; Raitses, Y.
Abstract:
Delineating the dominant processes responsible for nanomaterial synthesis in a plasma environment requires measurements of the precursor species contributing to the growth of nanostructures. We performed comprehensive measurements of spatial and temporal profiles of carbon dimers in sub-atmospheric-pressure carbon arc by laser-induced fluorescence. Measured spatial profiles of carbon dimers coincide with the growth region of carbon nanotubes (Fang et al 2016 Carbon 107 273-80) and vary depending on the arc operation mode, which is determined by the discharge current and the ablation rate of the graphite anode. The carbon dimer density profile exhibits large spatial and time variations due to motion of the arc core. A comparison of the experimental data with the 2D simulation results of self-consistent arc modeling shows a good agreement. The model predicts well the main processes determining spatial profiles of carbon dimers.
Type:
Dataset
Issue Date:
January 2018

166. Reductions in Retrieval Competition Predict the Benefit of Repeated Testing

Author(s):
Rafidi, Nicole S; Hulbert, Justin C; Brooks, Paula P; Norman, Kenneth A
Abstract:
Repeated testing (as opposed to repeated study) leads to improved long-term memory retention, but the mechanism underlying this improvement remains controversial. In this work, we test the hypothesis that retrieval practice benefits subsequent recall by reducing competition from related memories. This hypothesis implies that the degree of reduction in competition between retrieval practice attempts should predict subsequent memory for the practiced items. To test this prediction, we collected electroencephalography (EEG) data across two sessions. In the first session, participants practiced selectively retrieving exemplars from superordinate semantic categories (high competition), as well as retrieving the names of the superordinate categories from exemplars (low competition). In the second session, participants repeatedly studied and were then tested on Swahili-English vocabulary. One week after session two, participants were again tested on the vocabulary. We trained a within-subject classifier on the data from session one to distinguish high and low competition states. We then used this classifier to measure competition across multiple retrieval practice attempts in the second session. The degree to which competition decreased for a given vocabulary word predicted whether that item was subsequently remembered in the third session. These results are consistent with the hypothesis that repeated testing improves retention by reducing competition.
Type:
Dataset
Issue Date:
April 2018

167. Regarding the optimization of O1-mode ECRH and the feasibility of EBW startup on NSTX-U

Author(s):
Lopez, N; Poli, F
Abstract:
Recently published scenarios for fully non-inductive startup and operation on the National Spherical Torus eXperiment Upgrade (NSTX-U) (Menard et al 2012 Nucl. Fusion 52 083015) show Electron Cyclotron Resonance Heating (ECRH) as an important component in preparing a target plasma for efficient High Harmonic Fast Wave and Neutral Beam heating. The modeling of the propagation and absorption of EC waves in the evolving plasma is required to define the most effective window of operation, and to optimize the launcher geometry for maximal heating and current drive during this window. Here, we extend a previous optimization of O1-mode ECRH on NSTX-U to account for the full time-dependent performance of the ECRH using simulations performed with TRANSP. We find that the evolution of the density profile has a prominent role in the optimization by defining the time window of operation, which in certain cases may be a more important metric to compare launcher performance than the average power absorption. This feature cannot be captured by analysis on static profiles, and should be accounted for when optimizing ECRH on any device that operates near the cutoff density. Additionally, the utility of the electron Bernstein wave (EBW) in driving current and generating closed flux surfaces in the early startup phase has been demonstrated on a number of devices. Using standalone GENRAY simulations, we find that efficient EBW current drive is possible on NSTX-U if the injection angle is shifted below the midplane and aimed towards the top half of the vacuum vessel. However, collisional damping of the EBW is projected to be significant, in some cases accounting for up to 97% of the absorbed EBW power
Type:
Dataset
Issue Date:
June 2018

168. Regional hydroclimatic variability due to contemporary deforestation in southern Amazonia and associated boundary layer characteristics

Author(s):
Khanna, Jaya; Medvigy, David; Fisch, Gilberto; Neves, Theomar Trindade de Araújo Tiburtino
Abstract:
Amazonian deforestation causes systematic changes in regional dry season precipitation. Some of these changes at contemporary large scales (a few hundreds of kilometers) of deforestation have been associated with a ‘dynamical mesoscale circulation’, induced by the replacement of rough forest with smooth pasture. In terms of decadal averages, this dynamical mechanism yields increased precipitation in downwind regions and decreased precipitation in upwind regions of deforested areas. Daily, seasonal, and interannual variations in this phenomenon may exist, but have not yet been identified or explained. This study uses observations and numerical simulations to develop relationships between the dynamical mechanism and the local- and continental-scale atmospheric conditions across a range of time scales. It is found that the strength of the dynamical mechanism is primarily controlled by the regional-scale thermal and dynamical conditions of the boundary layer, and not by the continental- and global-scale atmospheric state. Lifting condensation level and wind speed within the boundary layer have large and positive correlations with the strength of the dynamical mechanism. The strength of these relationships depends on time scale and is strongest over the seasonal cycle. Overall, the dynamical mechanism is found to be strongest during times when the atmosphere is relatively stable. Hence, for contemporary large scales of deforestation this phenomenon is found to be the prevalent convective triggering mechanism during the dry and parts of transition seasons (especially during the dry-to-wet transition), significantly affecting the hydroclimate during this period.
Type:
Dataset and Software
Issue Date:
2018

169. Scenario Development During Commissioning Operations on the National Spherical Torus Experiment Upgrade

Author(s):
Battaglia, D.J.; Boyer, M.D.; Gerhardt, S.; Mueller, D.; Myers, C.E.; Guttenfelder, W.; Menard, J.E.; Sabbagh, S.A.; Scotti, F.; Bedoya, F.; Bell, R.E.; Berkery, J.W.; Diallo, A.; Ferraro, N.; Jaworski, M.A.; Kaye, S.M.; LeBlanc, B.P.; Ono, M.; Park, J.-K.; Podesta, M.; Raman, R.; Soukhanovskii, V.A.
Abstract:
The National Spherical Torus Experiment Upgrade (NSTX-U) will advance the physics basis required for achieving steady-state, high-beta, and high-confinement conditions in a tokamak by accessing high toroidal field (1 T) and plasma current (1.0 - 2.0 MA) in a low aspect ratio geometry (A = 1.6 - 1.8) with flexible auxiliary heating systems (12 MW NBI, 6 MW HHFW). This paper describes progress in the development of L- and H-mode discharge scenarios and the commissioning of operational tools in the first ten weeks of operation that enable the scientific mission of NSTX-U. Vacuum field calculations completed prior to operations supported the rapid development and optimization of inductive breakdown at different values of ohmic solenoid current. The toroidal magnetic field (B_T0 = 0.65 T) exceeded the maximum values achieved on NSTX and novel long-pulse L-mode discharges with regular sawtooth activity exceeded the longest pulses produced on NSTX (tpulse > 1.8s). The increased flux of the central solenoid facilitated the development of stationary L-mode discharges over a range of density and plasma current (Ip). H-mode discharges achieved similar levels of stored energy, confinement (H98y,2 > 1) and stability (beta_N/beta_N-nowall > 1) compared to NSTX discharges for Ip < 1 MA. High-performance H-mode scenarios require an L-H transition early in the Ip ramp-up phase in order to obtain low internal inductance (li) throughout the discharge, which is conducive to maintaining vertical stability at high elongation (kappa > 2.2) and achieving long periods of MHD quiescent operations. The rapid progress in developing L- and H-mode scenarios in support of the scientific program was enabled by advances in real-time plasma control, efficient error field identification and correction, effective conditioning of the graphite wall and excellent diagnostic availability.
Type:
Dataset
Issue Date:
April 2018

170. Two-dimensional full-wave simulations of waves in space and tokamak plasmas

Author(s):
Kim, E.-W.; Bertelli, N.; Johnson, J.R.; Valeo, E.; Hosea, J.; Perkins, R.
Abstract:
We illustrate the capabilities of a recently developed two-dimensional full wave code (FW2D) in space and tokamak plasmas by adopting various values of density, magnetic field configuration and strength as well as boundary shape. As example, we first showed fast compressional wave propagation in the inner magnetosphere is dramatically modified by a plasmaspheric plume at Earth's magnetosphere. The results show that wave energy is trapped in the plume showing a leaky eigenmode-like structure with plume, which is similar to the detected magnetosonic waves. We also performed simulations of high harmonic fast waves in the scrape-off layer (SOL) plasmas of the National Spherical Torus eXperiment (NSTX)/NSTX-Upgrade. Comparison the results with previous full-wave simulations show that although the FW2D code uses a cold plasma approximation, the electric field and the fraction of the power losses in the SOL plasmas show excellent consistency and agreement with the previous full wave simulations performed by the AORSA code.
Type:
Dataset
Issue Date:
October 2018

171. Whistler wave generation by anisotropic tail electrons during asymmetric magnetic reconnection in space and laboratory

Author(s):
Yoo, Jongsoo; Jara-almonte, J.; Yerger, Evan; Wang, Shan; Qian, Tony; Le, Ari; Ji, Hantao; Yamada, Masaaki; Fox, William; Kim, Eun-Hwa; Chen, Li-Jen; Gershman, Daniel
Abstract:
Whistler wave generation near the magnetospheric separatrix during reconnection at the dayside magnetopause is studied with data from the Magnetospheric Multiscale (MMS) mission. The dispersion relation of the whistler mode is measured for the first time near the reconnection region in space, which shows that whistler waves propagate nearly parallel to the magnetic field line. A linear analysis indicates that the whistler waves are generated by temperature anisotropy in the electron tail population. This is caused by loss of electrons with a high velocity parallel to the magnetic field to the exhaust region. There is a positive correlation between activities of whistler waves and the lower-hybrid drift instability (LHDI) both in laboratory and space, indicating the enhanced transport by LHDI may be responsible for the loss of electrons with a high parallel velocity.
Type:
Dataset
Issue Date:
August 2018

172. 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

173. 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

174. 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

175. 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

176. Amplification of local changes along the timescale processing hierarchy

Author(s):
Yeshurun, Yaara; Nguyen, Mai; Hasson, Uri
Abstract:
Small changes in word choice can lead to dramatically different interpretations of narratives. How does the brain accumulate and integrate such local changes to construct unique neural representations for different stories? In this study we created two distinct narratives by changing only a few words in each sentence (e.g. “he” to “she” or “sobbing” to “laughing”) while preserving the grammatical structure across stories. We then measured changes in neural responses between the two stories. We found that the differences in neural responses between the two stories gradually increased along the hierarchy of processing timescales. For areas with short integration windows, such as early auditory cortex, the differences in neural responses between the two stories were relatively small. In contrast, in areas with the longest integration windows at the top of the hierarchy, such as the precuneus, temporal parietal junction, and medial frontal cortices, there were large differences in neural responses between stories. Furthermore, this gradual increase in neural difference between the stories was highly correlated with an area’s ability to integrate information over time. Amplification of neural differences did not occur when changes in words did not alter the interpretation of the story (e.g. “sobbing” to “crying”). Our results demonstrate how subtle differences in words are gradually accumulated and amplified along the cortical hierarchy as the brain constructs a narrative over time.
Type:
Dataset
Issue Date:
August 2017

177. 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

178. 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

179. 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

180. Application of Townsend avalanche theory to tokamak startup by coaxial helicity injection

Author(s):
Hammond, K.C.; Raman, R.; Volpe, F.A.
Abstract:
Townsend avalanche theory is employed to model and interpret plasma initiation in NSTX by Ohmic heating and coaxial helicity injection (CHI). The model is informed by spatially resolved vacuum calculations of electric field and magnetic field line connection length in the poloidal cross-section. The model is shown to explain observations of Ohmic startup including the duration and location of breakdown. Adapting the model to discharges initiated by CHI offers insight into the causes of upper divertor (absorber) arcs in cases where the discharge fails to initiate in the lower divertor gap. Finally, upper and lower limits are established for vessel gas fill based on requirements for breakdown and radiation. It is predicted that CHI experiments on NSTX-U should be able to use as much as four times the amount of prefill gas employed in CHI experiments in NSTX. This should provide greater flexibility for plasma start-up, as the injector flux is projected to be increased in NSTX-U.
Type:
Dataset
Issue Date:
September 2017

181. Computation of Alfvèn eigenmode stability and saturation through a reduced fast ion transport model in the TRANSP tokamak transport code

Author(s):
Podesta, M.; Gorelenkova, M.; Gorelenkov, N.N.; White, R.B.
Abstract:
Alfvénic instabilities (AEs) are well known as a potential cause of enhanced fast ion transport in fusion devices. Given a specific plasma scenario, quantitative predictions of (i) expected unstable AE spectrum and (ii) resulting fast ion transport are required to prevent or mitigate the AE- induced degradation in fusion performance. Reduced models are becoming an attractive tool to analyze existing scenarios as well as for scenario prediction in time-dependent simulations. In this work, a neutral beam heated NSTX discharge is used as reference to illustrate the potential of a reduced fast ion transport model, known as kick model, that has been recently implemented for interpretive and predictive analysis within the framework of the time-dependent tokamak transport code TRANSP. Predictive capabilities for AE stability and saturation amplitude are first assessed, based on given thermal plasma profiles only. Predictions are then compared to experimental results, and the interpretive capabilities of the model further discussed. Overall, the reduced model captures the main properties of the instabilities and associated effects on the fast ion population. Additional information from the actual experiment enables further tuning of the model’s parameters to achieve a close match with measurements.
Type:
Dataset
Issue Date:
September 2017

182. Data for Nature Climate Change article 'Regional dry-season climate changes due to three decades of Amazonian deforestation'

Author(s):
Khanna, Jaya; Medvigy, David; Fueglistaler, Stephan; Walko, Robert
Abstract:
More than 20% Amazon rainforest has been cleared in the past three decades triggering important hydroclimatic changes. Small-scale (~few kilometers) deforestation in the 1980s has caused thermally-triggered atmospheric circulations that increase regional cloudiness and precipitation frequency. However, these circulations are predicted to diminish as deforestation increases. Here we use multi-decadal satellite records and numerical model simulations to show a regime shift in the regional hydroclimate accompanying increasing deforestation in Rondônia, Brazil. Compared to the 1980s, present-day deforested areas in downwind western Rondônia are found to be wetter than upwind eastern deforested areas during the local dry season. The resultant precipitation change in the two regions is approximately ±25% of the deforested area mean. Meso-resolution simulations robustly reproduce this transition when forced with increasing deforestation alone, showing a negligible role of large-scale climate variability. Furthermore, deforestation-induced surface roughness reduction is found to play an essential role in the present-day dry season hydroclimate. Our study illustrates the strong scale-sensitivity of the climatic response to Amazonian deforestation and suggests that deforestation is sufficiently advanced to have caused a shift from a thermally- to a dynamically-driven hydroclimatic regime.
Type:
Dataset and Software
Issue Date:
2017

183. Density perturbation mode structure of high frequency compressional and global Alfvén eigenmodes in the National Spherical Torus Experiment using a novel reflectometer analysis technique

Author(s):
Crocker, N.A.; Kubota, S.; Peebles, W.A.; Rhodes, T.L.; Fredrickson, E.D.; Belova, E.; Diallo, A.; LeBlanc, B.P.; Sabbagh, S.A.
Abstract:
Reflectometry measurements of compressional (CAE) and global (GAE) Alfvén eigenmodes are analyzed to obtain the amplitude and spatial structure of the density perturbations associated with the modes. A novel analysis technique developed for this purpose is presented. The analysis also naturally yields the amplitude and spatial structure of the density contour radial displacement, which is found to be 2–4 times larger than the value estimated directly from the reflectometer measurements using the much simpler ‘mirror approximation’. The modes were driven by beam ions in a high power (6 MW) neutral beam heated H-mode discharge (#141398) in the National Spherical Torus Experiment. The results of the analysis are used to assess the contribution of the modes to core energy transport and ion heating. The total displacement amplitude of the modes, which is shown to be larger than previously estimated (Crocker et al 2013 Nucl. Fusion 53 43017), is compared to the predicted threshold (Gorelenkov et al 2010 Nucl. Fusion 50 84012) for the anomalously high heat diffusion inferred from transport modeling in similar NSTX discharges. The results of the analysis also have strong implications for the energy transport via coupling of CAEs to kinetic Alfvén waves seen in simulations with the Hybrid MHD code (Belova et al 2015 Phys. Rev. Lett. 115 15001). Finally, the amplitudes of the observed CAEs fall well below the threshold for causing significant ion heating by stochastic velocity space diffusion (Gates et al 2001 Phys. Rev. Lett. 87 205003).
Type:
Dataset
Issue Date:
November 2017

184. Design of Faraday cup ion detectors built by thin film deposition

Author(s):
Szalkowski, G.A.; Darrow, D.S.; Cecil, F.E.
Abstract:
Thin film Faraday cup detectors can provide measurements of fast ion loss from magnetically confined fusion plasmas. These multilayer detectors can resolve the energy distribution of the lost ions in addition to giving the total loss rate. Prior detectors were assembled from discrete foils and insulating sheets. Outlined here is a design methodology for creating detectors using thin film deposition that are suited to particular scientific goals. The intention is to use detectors created by this method on JET and NSTX-U. The detectors will consist of alternating layers of aluminum and silicon dioxide, with layer thicknesses chosen to isolate energies of interest. Thin film deposition offers the advantage of relatively simple and more mechanically robust construction compared to other methods, as well as allowing precise control of film thickness. Furthermore, this depositional fabrication technique places the layers in intimate thermal contact, providing for three-dimensional conduction and dissipation of the ion-produced heating in the layers, rather than the essentially two-dimensional heat conduction in the discrete foil stack implementation.
Type:
Dataset
Issue Date:
January 2017

185. Drift kinetic effects on the plasma response in high beta spherical tokamak experiments

Author(s):
Wang, Z.R.; Park, J.-K.; Menard, J.E.; Liu, Y.Q.; Kaye, S.M.; Gerhardt, S.
Abstract:
High $\beta$ plasma response to the rotating n=1 external magnetic perturbations is numerically studied and compared with National Spherical Torus eXperiment (NSTX). The hybrid magnetohydrodynamic(MHD)-kinetic modeling shows the drift kinetic effects are important to resolve the disagreement of plasma response between the ideal MHD prediction and the NSTX experimental observation when plasma pressure reaches and exceeds the no-wall limit [F. Troyon et al., Plasma Phys. Control. Fusion \textbf{26}, 209 (1984)]. Since the external rotating fields and high plasma rotation are presented in NSTX experiments, the importance of resistive wall effect and plasma rotation on determining the plasma response is also identified, where the resistive wall suppresses the plasma response through the wall eddy current. The inertial energy, due to plasma rotation, destabilizes the plasma. The complexity of plasma response, in this study, indicates that MHD modeling, including comprehensive physics e.g. the drift kinetic effects, resistive wall and plasma rotation, is essential to reliably predict the plasma behavior in high beta spherical tokamak device.
Type:
Dataset
Issue Date:
September 2017

186. ELM elimination with Li powder injection in EAST discharges using the tungsten upper divertor

Author(s):
Maingi, R.; Hu, J.S.; Sun, Z.; Tritz, K.; Zuo, G.Z.; Xu, W.; Huang, M.; Meng, X.C.; Canik, J.M.; Diallo, A.; Lunsford, R.; Mansfield, D.K.; Osborne, T.H.; Gong, X.Z.; Wang, Y.F.; Li, Y.Y.
Abstract:
We report the first successful use of lithium (Li) to eliminate edge-localized modes (ELMs) with tungsten divertor plasma-facing components in the EAST device. Li powder injected into the scrape-off layer of the tungsten upper divertor successfully eliminated ELMs for 3-5 sec in EAST. The ELM elimination became progressively more effective in consecutive discharges at constant lithium delivery rates, and the divertor D-alpha baseline emission was reduced, both signatures of improved wall conditioning. A modest decrease in stored energy and normalized energy confinement was also observed, but the confinement relative to H98 remained well above 1, extending the previous ELM elimination results via Li injection into the lower carbon divertor in EAST [J.S. Hu et al., Phys. Rev. Lett. 114 (2015) 055001]. These results can be compared with recent observations with lithium pellets in ASDEX-Upgrade that failed to mitigate ELMs [P.T. Lang et al., Nucl. Fusion 57 (2017) 016030], highlighting one comparative advantage of continuous powder injection for real-time ELM elimination.
Type:
Dataset
Issue Date:
December 2017

187. Electron heating and energy inventory during asymmetric reconnection in a laboratory plasma

Author(s):
Yoo, Jongsoo; Na, Byungkeun; Jara-Almonte, Jonathan; Yamada, Maasaki; Ji, Hantao; Roytershteyn, V.; Argall, M. R.; Fox, W.; Chen, Li-Jen
Abstract:
Electron heating and the energy inventory during asymmetric reconnection are studied in the laboratory plasma with a density ratio of about 8 across the current sheet. Features of asymmetric reconnection such as the large density gradients near the low-density-side separatrices, asymmetric in-plane electric field, and bipolar out-of-plane magnetic field are observed. Unlike the symmetric case, electrons are also heated near the low-density-side separatrices. The measured parallel electric field may explain the observed electron heating. Although large fluctuations driven by lower-hybrid drift instabilities are also observed near the low-density-side separatrices, laboratory measurements and numerical simulations reported here suggest that they do not play a major role in electron energization. The average electron temperature increase in the exhaust region is proportional to the incoming magnetic energy per an electron/ion pair but exceeds scalings of the previous space observations. This discrepancy is explained by differences in the boundary condition and system size. The profile of electron energy gain from the electric field shows that there is additional electron energy gain associated with the electron diamagnetic current besides a large energy gain near the X-line. This additional energy gain increases electron enthalpy, not the electron temperature. Finally, a quantitative analysis of the energy inventory during asymmetric reconnection is conducted. Unlike the symmetric case where the ion energy gain is about twice more than the electron energy gain, electrons and ions obtain a similar amount of energy during asymmetric reconnection.
Type:
Dataset
Issue Date:
August 2017

188. Energetic-particle-modified global Alfven eigenmodes

Author(s):
Lestz, J.B.; Belova, E.V.; Gorelenkov, N.N.
Abstract:
Fully self-consistent hybrid MHD/particle simulations reveal strong energetic particle modifications to sub-cyclotron global Alfven eigenmodes (GAE) in low-aspect ratio, NSTX-like conditions. Key parameters defining the fast ion distribution function -- the normalized injection velocity v_0/v_A and central pitch -- are varied in order to study their influence on the characteristics of the excited modes. It is found that the frequency of the most unstable mode changes significantly and continuously with beam parameters, in accordance with the Doppler-shifted cyclotron resonances which drive the modes, and depending most substantially on v_0/v_A. This unexpected result is present for both counter-propagating GAEs, which are routinely excited in NSTX, and high frequency co-GAEs, which have not been previously studied. Large changes in frequency without clear corresponding changes in mode structure could indicate the existence of a new energetic particle mode, referred to here as an energetic-particle-modified GAE (EP-GAE). Additional simulations conducted for a fixed MHD equilibrium demonstrate that the GAE frequency shift cannot be explained by the equilibrium changes due to energetic particle effects.
Type:
Dataset
Issue Date:
December 2017

189. Energy Exchange Dynamics across L-H Transitions in NSTX

Author(s):
Diallo, A.; Banerjee, S.; Zweben, S.; Stoltzfus-Dueck, T.
Abstract:
We studied the energy exchange dynamics across the low-to-high-confinement (L-H) in NSTX discharges using the gas-puff imaging (GPI) diagnotic. The investigation focused on the energy exchange between flows and turbulence, to help clarify the mechanism of the L-H transition. We apply this study to three type 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-programming decomposition), velocity fields of a 24 $\times$ 30 cm GPI view during the L-H transition were obtained with good spatial ($\sim$1 cm) and temporal ($\sim$2.5 $\mu$s) resolutions. Analysis using these velocity fields shows that the production term is systematically negative just prior to the L-H transition indicating transfer from mean flows to turbulence, which is inconsistent with the predator-prey paradigm. Moreover, using the inferred absolute value of the production term, an estimate of the L-H transition duration is found to be 25 ms, which is much larger than the measured duration. 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 turbulence depletion mechanism may not be playing an important role in the transition to the H-mode. Although the Reynolds work is 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.
Type:
Dataset
Issue Date:
June 2017

190. Experimental calibration procedures for rotating Lorentz-force flowmeters

Author(s):
Hvasta, M. G.; Slighton, N. T.; Kolemen, E.; Fisher, A. E.
Abstract:
Rotating Lorentz-force flowmeters are a novel and useful technology with a range of applications in a variety of different industries. However, calibrating these flowmeters can be challenging, time-consuming, and expensive. In this paper, simple calibration procedures for rotating Lorentz-force flowmeters are presented. These procedures eliminate the need for expensive equipment, numerical modeling, redundant flowmeters, and system down-time. The calibration processes are explained in a step-by-step manner and compared to experimental results.
Type:
Dataset
Issue Date:
2017

191. Feedback control design for non-inductively sustained scenarios in NSTX-U using TRANSP

Author(s):
Boyer, M.D.; Andre, R.G.; Gates, D.A.; Gerhardt, S.P.; Menard, J.E.; Poli, F.M.
Abstract:
This paper examines a method for real-time control of non-inductively sustained scenarios in NSTX-U by using TRANSP, a time-dependent integrated modeling code for prediction and interpretive analysis of tokamak experimental data, as a simulator. The actuators considered for control in this work are the six neutral beam sources and the plasma boundary shape. To understand the response of the plasma current, stored energy, and central safety factor to these actuators and to enable systematic design of control algorithms, simulations were run in which the actuators were modulated and a linearized dynamic response model was generated. A multi-variable model-based control scheme that accounts for the coupling and slow dynamics of the system while mitigating the effect of actuator limitations was designed and simulated. Simulations show that modest changes in the outer gap and heating power can improve the response time of the system, reject perturbations, and track target values of the controlled values.
Type:
Dataset
Issue Date:
June 2017

192. Full-wave simulations of ICRF heating regimes in toroidal plasmas with non-Maxwellian distribution functions

Author(s):
Bertelli, N; Valeo, E.J.; Green, D.L.; Gorelenkova, M.; Phillips, C.K.; Podesta, M.; Lee, J.P.; Wright, J.C.; Jaeger, E.
Abstract:
At the power levels required for significant heating and current drive in magnetically-confined toroidal plasma, modification of the particle distribution function from a Maxwellian shape is likely [T. H. Stix, Nucl. Fusion, 15 737 (1975)], with consequent changes in wave propagation and in the location and amount of absorption. In order to study these effects computationally, both the finite-Larmor-radius and the high-harmonic fast wave (HHFW), versions of the full-wave, hot-plasma toroidal simulation code TORIC [M. Brambilla, Plasma Phys. Control. Fusion 41, 1 (1999) and M. Brambilla, Plasma Phys. Control. Fusion 44, 2423 (2002)], have been extended to allow the prescription of arbitrary velocity distributions of the form f(v||, v_perp, psi , theta). For hydrogen (H) minority heating of a deuterium (D) plasma with anisotropic Maxwellian H distributions, the fractional H absorption varies significantly with changes in parallel temperature but is essentially independent of perpendicular temperature. On the other hand, for HHFW regime with anisotropic Maxwellian fast ion distribution, the fractional beam ion absorption varies mainly with changes in the perpendicular temperature. The evaluation of the wave-field and power absorption, through the full wave solver, with the ion distribution function provided by either aMonte-Carlo particle and Fokker-Planck codes is also examined for Alcator C-Mod and NSTX plasmas. Non-Maxwellian effects generally tends to increase the absorption with respect to the equivalent Maxwellian distribution.
Type:
Dataset
Issue Date:
May 2017

193. Gas Puff Imaging Diagnostics of Edge Plasma Turbulence in Magnetic Fusion Devices

Author(s):
Zweben, S.J.; Terry, J.L.; Stotler, D.P.; Maqueda, R.J.
Abstract:
Gas puff imaging (GPI) is a diagnostic of plasma turbulence which uses a puff of neutral gas at the plasma edge to increase the local visible light emission for improved space-time resolution of plasma fluctuations. This paper reviews gas puff imaging diagnostics of edge plasma turbulence in magnetic fusion research, with a focus on the instrumentation, diagnostic cross-checks, and interpretation issues. The gas puff imaging hardware, optics, and detectors are described for about 10 GPI systems implemented over the past ~15 years. Comparison of GPI results with other edge turbulence diagnostic results are described and many common features are observed. Several issues in the interpretation of GPI measurements are discussed, and potential improvements in hardware and modeling are suggested.
Type:
Dataset
Issue Date:
April 2017

194. Injected mass deposition thresholds for lithium granule instigated triggering of edge localized modes on EAST

Author(s):
Lunsford, R.; Sun, Z.; Maingi, R.; Hu, J.S.; Mansfield, D.; Xu, W.; Zuo, G.Z.; Diallo, A.; Osborne, T.; Tritz, K.; Canik, J.; Huang, M.; Meng, X.C.; Gong, X.Z.; Wan, B.N.; Li, J.G.
Abstract:
The ability of an injected lithium granule to promptly trigger an edge localized mode (ELM) has been established in multiple experiments. By horizontally injecting granules ranging in diameter from 200 microns to 1mm in diameter into the low field side of EAST H-mode discharges we have determined that granules with diameter > 600 microns are successful in triggering ELMs more than 95% of the time. It was also demonstrated that below 600 microns the triggering efficiency decreased roughly with granule size. Granules were radially injected from the outer midplane with velocities ~ 80 m/s into EAST upper single null discharges with an ITER like tungsten monoblock divertor. These granules were individually tracked throughout their injection cycle in order to determine their efficacy at triggering an ELM. For those granules of sufficient size, ELM triggering was a prompt response to granule injection. By simulating the granule injection with an experimentally benchmarked neutral gas shielding (NGS) model, the ablatant mass deposition required to promptly trigger an ELM is calculated and the fractional mass deposition is determined.
Type:
Dataset
Issue Date:
December 2017

195. Kinetic neoclassical calculations of impurity radiation profiles

Author(s):
Stotler, D.P.; Battaglia, D.J.; Hager, R.; Kim, K.; Koskela, T.; Park, G.; Reinke, M.L.
Abstract:
Modifications of the drift-kinetic transport code XGC0 to include the transport, ionization, and recombination of individual charge states, as well as the associated radiation, are described. The code is first applied to a simulation of an NSTX H-mode discharge with carbon impurity to demonstrate the approach to coronal equilibrium. The effects of neoclassical phenomena on the radiated power profile are examined sequentially through the activation of individual physics modules in the code. Orbit squeezing and the neoclassical inward pinch result in increased radiation for temperatures above a few hundred eV and changes to the ratios of charge state emissions at a given electron temperature. Analogous simulations with a neon impurity yield qualitatively similar results.
Type:
Dataset
Issue Date:
2017

196. M3D-C1 simulations of the plasma response to RMPs in NSTX-U single-null and snowflake divertor configurations

Author(s):
Canal, G.P.; Ferraro, N.M.; Evans, T.E.; Osborne, T.H.; Menard, J.E.; Ahn, J.-W.; Maingi, R.; Wingen, A.; Ciro, D.; Frerichs, H.; Schmitz, O.; Soukhanovskii, V.; Waters, I.
Abstract:
Non-axisymmetric control coils and the so-called snowflake divertor configuration are two potential solutions proposed to solve two separate outstanding issues on the path towards self-sustained burning plasma operations, namely the transient energy bursts caused by edge localized modes and the steady state heat exhaust problem. In a reactor, these two proposed solutions would have to operate simultaneously and it is, therefore, important to investigate their compatibility and to identify possible conflicts that could prevent them from operating simultaneously. In this work, single- and two-fluid resistive magnetohydrodynamic calculations are used to investigate the effect of externally applied magnetic perturbations on the snowflake divertor configuration. The calculations are based on simulated NSTX-U plasmas and the results show that additional and longer magnetic lobes are created in the null-point region of the snowflake configuration, compared to those in the conventional single-null. The intersection of these longer and additional lobes with the divertor plates are expected to cause more striations in the particle and heat flux target profiles. In addition, the results indicate that the size of the magnetic lobes, in both single-null and snowflake configurations, are more sensitive to resonant magnetic perturbations than to non-resonant magnetic perturbations. The results also suggest that lower values of current in non-axisymmetric control coils would be required to suppress edge localized modes in plasmas with the snowflake configuration.
Type:
Dataset
Issue Date:
July 2017

197. Modeling of Lithium Granule Injection in NSTX with M3D-C1

Author(s):
Fil, A.; Kolemen, E.; Bortolon, A.; Ferraro, N.; Jardin, S.; Parks, P.B.; Lunsford, R.; Maingi, R.
Abstract:
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.
Type:
Dataset
Issue Date:
January 2017

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

Author(s):
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.
Type:
Dataset
Issue Date:
July 2017

200. Neutral recycling effects on ITG turbulence

Author(s):
Stotler, D.P.; Lang, J.; Chang, C.S.; Churchill, R.M.; Ku, S.-H.
Abstract:
The effects of recycled neutral atoms on tokamak ion temperature gradient (ITG) driven turbulence have been investigated in a steep edge pedestal, magnetic separatrix configuration, with the full-f edge gryokinetic code XGC1. Ion temperature gradient turbulence is the most fundamental and robust edge plasma instability, having a long radial correlation length and an ability to impact other forms of pedestal turbulence. The neutral atoms enhance the ITG turbulence, first, by increasing the ion temperature gradient in the pedestal via the cooling effects of charge exchange and, second, by a relative reduction in the ExB shearing rate.
Type:
Dataset
Issue Date:
August 2017