F. M. Laggner, A. Diallo, B. P. LeBlanc, R. Rozenblat, G. Tchilinguirian, E.Kolemen, the NSTX-U team
A detailed description of a prototype setup for real-time (rt) Thomson scattering (TS) analysis is presented and implemented in the multi-point Thomson scattering (MPTS) diagnostic system at the National Spherical Torus Experiment Upgrade(NSTX-U). The data acquisition hardware was upgraded with rt capable electronics (rt-analog digital converters (ADCs) and a rt server) that allow for fast digitization of the laser pulse signal of eight radial MPTS channels. In addition, a new TS spectrum analysis software for a rapid calculation of electron temperature (Te) and electron density (ne) was developed. Testing of the rt hardware and data analysis soft-ware was successfully completed and benchmarked against the standard, post-shot evaluation. Timing tests were performed showing that the end-to-end processing time was reproducibly below 17 ms for the duration of at least 5 s, meeting a 60 Hz deadline by the laser pulse repetition rate over the length of a NSTX-U discharge. The presented rt framework is designed to be scalable in system size, i.e. incorporation of additional radial channels by solely adding additional rt capable hardware. Furthermore, it is scalable in its operation duration and was continuously run for up to 30 min, making it an attractive solution for machines with long discharge duration such as advanced, non-inductive tokamaks or stellarators.
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.
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.
Zweben SJ, Myra JR, Diallo A, Russell DA, Scotti F, Stotler DP
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
Yang, Yuan; Pan, Ming; Beck, Hylke; Fisher, Colby; Beighley, R. Edward; Kao, Shih-Chieh; Hong, Yang; Wood, Eric
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.
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.
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).
Antony, James W.; Cheng, Larry Y.; Brooks, Paula P.; Paller, Ken A.; Norman, Kenneth A.
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.
Vail, P. J.; Boyer, M. D.; Welander, A. S.; Kolemen, E.; U.S. Department of Energy contract number DE-AC02-09CH11466
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.
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.
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.
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.
\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.\
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.
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.
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.
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.
The data are 4554 light curves derived from images taken of the globular cluster M4 by the Kepler space telescope during the K2 portion of its mission, specifically during Campaign 2 of that mission, which occurred in 2014. A total of 3856 images were taken over approximately three months at a cadence of approximately half an hour. The purpose of these observations was to find stars and other objects that vary in brightness over time --- variable stars. Also included is a table with associated information for each of the 4554 objects and their light curves.