Search

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

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

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

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

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

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

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

109. An extended MHD study of the 16 October 2015 MMS diffusion region crossing

Abstract:
The Magnetospheric Multiscale (MMS) mission has given us unprecedented access to high cadence particle and field data of magnetic reconnection at Earth's magnetopause. MMS first passed very near an X-line on 16 October 2015, the Burch event, and has since observed multiple X-line crossings. Subsequent 3D particle-in-cell (PIC) modeling efforts of and comparison with the Burch event have revealed a host of novel physical insights concerning magnetic reconnection, turbulence induced particle mixing, and secondary instabilities. In this study, we employ the Gkeyll simulation framework to study the Burch event with different classes of extended, multi-fluid magnetohydrodynamics (MHD), including models that incorporate important kinetic effects, such as the electron pressure tensor, with physics-based closure relations designed to capture linear Landau damping. Such fluid modeling approaches are able to capture different levels of kinetic physics in global simulations and are generally less costly than fully kinetic PIC. We focus on the additional physics one can capture with increasing levels of fluid closure refinement via comparison with MMS data and existing PIC simulations. In particular, we find that the ten-moment model well captures the agyrotropic structure of the pressure tensor in the vicinity of the X-line and the magnitude of anisotropic electron heating observed in MMS and PIC simulations. However, the ten-moment model has difficulty resolving the lower hybrid drift instability, which has been observed to plays a fundamental role in heating and mixing electrons in the current layer.
Type:
Dataset

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

111. Calibrated VIC Model Parameters over CONUS

Author(s):
Yang, Yuan; Pan, Ming; Beck, Hylke; Fisher, Colby; Beighley, R. Edward; Kao, Shih-Chieh; Hong, Yang; Wood, Eric
Abstract:
Conventional basin-by-basin approaches to calibrate hydrologic models are limited to gauged basins and typically result in spatially discontinuous parameter fields. Moreover, the consequent low calibration density in space falls seriously behind the need from present-day applications like high resolution river hydrodynamic modeling. In this study we calibrated three key parameters of the Variable Infiltration Capacity (VIC) model at every 1/8° grid-cell using machine learning-based maps of four streamflow characteristics for the conterminous United States (CONUS), with a total of 52,663 grid-cells. This new calibration approach, as an alternative to parameter regionalization, applied to ungauged regions too. A key difference made here is that we tried to regionalize physical variables (streamflow characteristics) instead of model parameters whose behavior may often be less well understood. The resulting parameter fields no longer presented any spatial discontinuities and the patterns corresponded well with climate characteristics, such as aridity and runoff ratio. The calibrated parameters were evaluated against observed streamflow from 704/648 (calibration/validation period) small-to-medium-sized catchments used to derive the streamflow characteristics, 3941/3809 (calibration/validation period) small-to-medium-sized catchments not used to derive the streamflow characteristics) as well as five large basins. Comparisons indicated marked improvements in bias and Nash-Sutcliffe efficiency. Model performance was still poor in arid and semiarid regions, which is mostly due to both model structural and forcing deficiencies. Although the performance gain was limited by the relative small number of parameters to calibrate, the study and results here served as a proof-of-concept for a new promising approach for fine-scale hydrologic model calibrations.
Type:
Dataset
Issue Date:
4 September 2019

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

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

114. Dust and Starlight Maps for Galaxies in the KINGFISH Sample

Author(s):
Aniano, G.; Draine, B.T.; Hunt, L.K.; Sandstrom, K.; Calzetti, D.; Kennicutt, R.C.; Dale, D.A.; Galametz, M.; Gordon, K.D.; Leroy, A.K.; Smith, J.-D.T.; Roussel, H.; Sauvage, M.; Walter, F.; Armus, L.; Bolatto, A.D.; Boquien, M.; Crocker, A.; De Looze, I.; Donovan Meyer, J.; Helou, G.; Hinz, J.; Johnson, B.D.; Koda, J.; Miller, A.; Montiel, E.; Murphy, E.J.; Relano, M.; Rix, H.-W.; Schinnerer, E.; Skibba, R.; Wolfire, M.G.; Engelbracht, C.W.
Abstract:
Dust and starlight have been modeled for the KINGFISH project galaxies. For each pixel in each galaxy, we estimate: (1) dust surface density; (2) q_PAH, the dust mass fraction in PAHs; (3) distribution of starlight intensities heating the dust; (4) luminosity emitted by the dust; and (5) dust luminosity from regions with high starlight intensity. The modeling is as described in the paper "Modeling Dust and Starlight in Galaxies Observed by Spitzer and Herschel: The KINGFISH Sample", by G. Aniano, B.T. Draine, L.K. Hunt, K. Sandstrom, D. Calzetti, R.C. Kennicutt, D.A, Dale, and 26 other authors, accepted for publication in The Astrophysical Journal.
Type:
Dataset and Image

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

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

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

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

119. Instability of shear flows in spatially periodic domains

Abstract:
Force-driven parallel shear flow in a spatially periodic domain is shown to be linearly unstable with respect to both the Reynolds number and the domain aspect ratio. This finding is confirmed by computer simulations, and a simple expression is derived to determine stable flow conditions. Periodic extensions of Couette and Poiseuille flows are unstable at Reynolds numbers two orders of magnitude smaller than their aperiodic equivalents because the periodic boundaries impose fundamentally different constraints. This instability has important implications for designing computational models of nonlinear dynamic processes with periodicity.
Type:
Dataset

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

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

122. Methane emissions from oil and gas platforms in the North Sea.

Author(s):
Riddick, Stuart; Mauzerall, Denise
Abstract:
Since 1850 the concentration of atmospheric methane (CH4), a potent greenhouse gas, has more than doubled. Recent studies suggest that emission inventories may be missing sources and underestimating emissions. To investigate whether offshore oil and gas platforms leak CH4 during normal operation, we measured CH4 mole fractions around eight oil and gas production platforms in the North Sea which were neither flaring gas nor off-loading oil. We use the measurements from summer 2017, along with meteorological data, in a Gaussian plume model to estimate CH4 emissions from each platform. We find CH4 mole fractions of between 11 and 370 ppb above background concentrations downwind of the platforms measured, corresponding to a median CH4 emission of 6.8 g CH4 s-1 for each platform, with a range of 2.9 to 22.3 g CH4 s-1. When matched to production records, during our measurements individual platforms lost between 0.04% and 1.4% of gas produced with a median loss of 0.23%. When the measured platforms are considered collectively, (i.e. the sum of platforms’ emission fluxes weighted by the sum of the platforms’ production), we estimate the CH4 loss to be 0.19% of gas production. These estimates are substantially higher than the emissions most recently reported to the National Atmospheric Emission Inventory (NAEI) for total CH4 loss from United Kingdom platforms in the North Sea. The NAEI reports CH4 losses from the offshore oil and gas platforms we measured to be 0.13% of gas production, with most of their emissions coming from gas flaring and offshore oil loading, neither of which were taking place at the time of our measurements. All oil and gas platforms we observed were found to leak CH4 during normal operation and much of this leakage has not been included in UK emission inventories. Further research is required to accurately determine total CH4 leakage from all offshore oil and gas operations and to properly include the leakage in national and international emission inventories.
Type:
Dataset
Issue Date:
February 2019

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

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

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

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

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

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

129. Sowing the Seeds for More Usable Web Archives: A Usability Study of Archive-It

Author(s):
Abrams, Samantha; Antracoli, Alexis; Appel, Rachel; Caust-Ellenbogen, Celia; Dennison, Sarah; Duncan, Sumitra; Ramsay, Stefanie
Abstract:
In 2017, seven members of the Archive-It Mid-Atlantic Users Group (AITMA) conducted a study of 14 subjects representative of their stakeholder populations to assess the usability of Archive-It, a web archiving subscription service of the Internet Archive. While Archive-It is the most widely-used tool for web archiving, little is known about how users interact with the service. This study intended to teach us what users expect from web archives, which exist as another form of archival material. End-user subjects executed four search tasks using the public Archive-It interface and the Wayback Machine to access archived information on websites from the facilitators’ own harvested collections and provide feedback about their experiences. The tasks were designed to have straightforward pass or fail outcomes, and the facilitators took notes on the subjects’ behavior and commentary during the sessions. Overall, participants reported mildly positive impressions of Archive-It public user interface based on their session. The study identified several key areas of improvement for the Archive-It service pertaining to metadata options, terminology display, indexing of dates, and the site’s search box.
Type:
Dataset
Issue Date:
2019

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

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

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

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

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

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

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

138. Collapse of reacted fracture surface decreases permeability and frictional strength

Author(s):
Peters, Catherine A.; Spokas, Kasparas
Abstract:
Geochemical and geomechanical perturbations of the subsurface caused by the injection of fluids present the risk of leakage and seismicity. This study investigated how flow of acidic fluids affects hydraulic and frictional properties of fractures using experiments with 3.8 cm-long specimens of Eagle Ford shale, a laminated shale with carbonate-rich strata. In low-pressure flow cells, one set of samples was exposed to an acidic brine and another set was exposed to a neutral brine. X-ray computed tomography and x-ray fluorescence analysis revealed that samples exposed to the acidic brine were calcite-depleted and had developed a porous altered layer, while the other set showed little evidence of alteration. After reaction, samples were compacted and sheared in a triaxial cell that supplied normal stress and differential pore pressure at prescribed sliding velocities, independently measuring friction and permeability. During the initial compaction, the porous altered layer collapsed into fine particles that filled the fracture aperture. This effectively impeded flow and sealed the fracture, resulting in a decrease in fracture permeability by 1 to 2 orders of magnitude relative to the compressed unaltered fractures. During shear, the collapsed layer of fine-grained particles prevented the formation of interlocking micro-asperities resulting in lower frictional strength. With regard to subsurface risks, this study showcases how coupled geochemical and geomechanical processes could favorably seal fractures to inhibit leakage, but also could increase the likelihood of induced seismicity. These findings have important implications for geological carbon sequestration, pressurized fluid energy storage, geothermal energy, and other subsurface technologies.
Type:
Dataset
Issue Date:
2018

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

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

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

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

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

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

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

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

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

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

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