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}.
Lampert,Mate; Diallo,Ahmed; Myra,James R.; Zweben, Stewart J.
Abstract:
Edge localized modes (ELMs) are routinely observed in H-mode plasma regimes of the National Spherical Torus Experiment (NSTX). Due to the explosive nature of the instability, only diagnostics with high temporal and spatial resolution could provide a detailed insight into the dynamics associated with the ELMs. Gas-puff imaging (GPI) at NSTX provides 2D measurements of the magnetic field aligned fluctuations (e.g. ELM filaments) in the scrape-off layer and the at the plasma edge with 2.5 us temporal and 10 mm optical resolution.A novel analysis technique was developed to estimate the frame-by-frame velocities and the spatial parameters of the dominant structures associated with the ELMs. The analysis was applied to single ELM events to characterize the ELM crash dynamics, and then extended to a database of 169 ELM events.Statistical analysis was performed in order to find the characterizing dynamics of the ELM crash. The results show that on average an ELM crash consists of a filament with a circular cross-section which is propelled outwards with a characterizing peak radial velocity of ~3.3 km/s. The radial velocity was found to be linearly dependent on the distance of the filament from the separatrix, which has never been seen before. The ELM filament is characterized by propagation in the ion-diamagnetic direction poloidally with a peak velocity of 11.4 km/s. The ELM crash lasts for approximately 100us until the radial propulsion settles back to the pre-ELM level. The experimental findings were compared with analytical theory. Two possible mechanisms were identified for explaining the observations: the curvature interchange model and the current-filament interaction model.
The engineering limits of plasma facing components (PFCs) constrain the allowable operational space of tokamaks. Poorly managed heat fluxes that push the PFCs beyond their limits not only degrade core plasma performance via elevated impurities, but can also result in PFC failure due to thermal stresses or melting. Simple axisymmetric assumptions fail to capture the complex interaction between 3D PFC geometry and 2D or 3D plasmas. This results in fusion systems that must either operate with increased risk or reduce PFC loads, potentially through lower core plasma performance, to maintain a nominal safety factor. High precision 3D heat flux predictions are necessary to accurately ascertain the state of a PFC given the evolution of the magnetic equilibrium. A new code, the Heat flux Engineering Analysis Toolkit (HEAT), has been developed to provide high precision 3D predictions and analysis for PFCs. HEAT couples many otherwise disparate computational tools together into a single open source python package. Magnetic equilibrium, engineering CAD, finite volume solvers, scrape off layer plasma physics, visualization, high performace computing, and more, are connected in a single web-based user interface. Linux users may use HEAT without any software prerequisites via an appImage. This manuscript introduces HEAT, discusses the software architecture, presents first HEAT results, and outlines physics modules in development.
The dielectric function for "Astrodust" grain material is provided for different assumed values of the dust grain shape (spheroid axis ratio) and porosity (vacuum fraction), and fraction of the interstellar iron present as metallic inclusions. For each case, the dielectric function is obtained by requiring that the grains reproduce the observed infrared opacity, and match to a physically reasonable dielectric function at 1 micron, and extending to X-ray energies. The derived dielectric functions satisfy the Kramers-Kronig relations. Dielectric functions are provided from 1 Angstrom to 5 cm (12.4 keV to 2.59e-5 eV).
For each dielectric function, we also calculate absorption and scattering corss sections for spheroidal grains, for three orientations of the grain relative to incident linearly-polarized light, for wavelengths from the Lyman limit (0.0912 micron) to the microwave (4 cm), and grain "effective radii" a_eff from 3.162A to 5.012 micron.
Notterman, Daniel A; Schneper, Lisa M; Drake, Amanda; Piyasena, Chinthika
Abstract:
This entry contains the data used in the PLOS ONE publication entitled, "Characteristics of salivary telomere length shortening in preterm infants" by Schneper et al. The objective of the study was to examine the association between gestational age, telomere length (TL) and rate of shortening in newborns. Genomic DNA was isolated from buccal samples of 39 term infants at birth and one year and 32 preterm infants at birth, term-adjusted age (40 weeks post-conception) and age one-year corrected for gestational duration. Telomere length was measured by quantitative real-time PCR. Demographic and clinical data were collected during clinic or research visits and from hospital records. Socioeconomic status was estimated using the deprivation category (DEPCAT) scores derived from the Carstairs score of the subject's postal code.
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.
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.
The history of organismal evolution, seawater chemistry, and paleoclimate is recorded in layers of carbonate sedimentary rock. Meter-scale cyclic stacking patterns in these carbonates often are interpreted as representing sea level change. A reliable sedimentary proxy for eustasy would be profoundly useful for reconstructing paleoclimate, since sea level responds to changes in temperature and ice volume. However, the translation from water depth to carbonate layering has proven difficult, with recent surveys of modern shallow water platforms revealing little correlation between carbonate facies (i.e., grain size, sedimentary bed forms, ecology) and water depth. We train a convolutional neural network with satellite imagery and new field observations from a 3,000 km2 region northwest of Andros Island (Bahamas) to generate a facies map with 5 m resolution. Leveraging a newly-published bathymetry for the same region, we test the hypothesis that one can extract a signal of water depth change, not simply from individual facies, but from sequences of facies transitions analogous to vertically stacked carbonate strata. Our Hidden Markov Model (HMM) can distinguish relative sea level fall from random variability with ∼90% accuracy. Finally, since shallowing-upward patterns can result from local (autogenic) processes in addition to forced mechanisms such as eustasy, we search for statistical tools to diagnose the presence or absence of external forcings on relative sea level. With a new data-driven forward model that simulates how modern facies mosaics evolve to stack strata, we show how different sea level forcings generate characteristic patterns of cycle thicknesses in shallow carbonates, providing a new tool for quantitative reconstruction of ancient sea level conditions from the geologic record.
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.
Schwartz, Jacob; Emdee, Eric; Goldston, Robert; Jaworski, Michael
Abstract:
The lithium vapor box divertor is a potential solution for power exhaust in toroidal confinement devices. The divertor plasma interacts with a localized, dense cloud of lithium vapor, leading to volumetric radiation, cooling, recombination, and detachment. To minimize contamination of the core plasma, lithium vapor is condensed on cool (300°C to 400°C) baffles upstream of the detachment point. Before implementing this in a toroidal plasma device with a slot divertor geometry, we consider an experiment with a scaled baffled-pipe geometry in the high-power linear plasma device Magnum-PSI. Three 15 cm-scale open cylinders joined by 6 cm diameter ‘nozzles’ are positioned on the plasma beam axis upstream of a target. The central box may be loaded with several tens of grams of lithium, which can be evaporated at 650°C to produce a vapor predicted, using a simple plasma-neutral interaction model, to be dense enough to cause volumetric detachment in the plasma. The power delivered to the target and box walls as measured by increases in their temperatures after a 10 s plasma pulse can be compared to determine the effectiveness of the vapor in detaching the plasma. Direct Simulation Monte Carlo simulations are performed to estimate the flow rates of lithium vapor between the boxes and to estimate the trapping of H2 delivered by the plasma in the boxes, which could inadvertently lead to detachment. Details of the geometry, simulations, and possible diagnostic techniques are presented.
These data include 39 structured interview transcripts. Each case is someone who worked at the time for Uber, UberEats, Lyft, and/or Amazon Flex (Amazon’s contractor delivery service). These data were collected between July and September 2019. All but one of the interviews occurred over the phone. My questions are focused on the structure of their gig work jobs and the technology they used at work or expected to use at work in the future. I included a description of the data, the recruitment methods, and the discussion guide in this ReadMe file.
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.
Amazonian deforestation causes systematic changes in regional dry season precipitation. Some of these changes at contemporary large scales (a few hundreds of kilometers) of deforestation have been associated with a ‘dynamical mesoscale circulation’, induced by the replacement of rough forest with smooth pasture. In terms of decadal averages, this dynamical mechanism yields increased precipitation in downwind regions and decreased precipitation in upwind regions of deforested areas. Daily, seasonal, and interannual variations in this phenomenon may exist, but have not yet been identified or explained. This study uses observations and numerical simulations to develop relationships between the dynamical mechanism and the local- and continental-scale atmospheric conditions across a range of time scales. It is found that the strength of the dynamical mechanism is primarily controlled by the regional-scale thermal and dynamical conditions of the boundary layer, and not by the continental- and global-scale atmospheric state. Lifting condensation level and wind speed within the boundary layer have large and positive correlations with the strength of the dynamical mechanism. The strength of these relationships depends on time scale and is strongest over the seasonal cycle. Overall, the dynamical mechanism is found to be strongest during times when the atmosphere is relatively stable. Hence, for contemporary large scales of deforestation this phenomenon is found to be the prevalent convective triggering mechanism during the dry and parts of transition seasons (especially during the dry-to-wet transition), significantly affecting the hydroclimate during this period.
Measuring free-surface, liquid-metal flow velocity is challenging to do in a reliable and accurate manner. This paper presents a non-invasive, easily-calibrated method of measuring the surface velocities of open-channel liquid-metal flows using an IR camera. Unlike other spatially-limited methods, this IR camera particle tracking technique provides full field-of-view data that can be used to better understand open-channel flows and determine surface boundary conditions. This method could be implemented and automated for a wide range of liquid-metal experiments, even if they operate at high-temperatures or within strong magnetic fields.
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.
Hvasta, M. G.; Slighton, N. T.; Kolemen, E.; Fisher, A. E.
Abstract:
Rotating Lorentz-force flowmeters are a novel and useful technology with a range of
applications in a variety of different industries. However, calibrating these flowmeters can
be challenging, time-consuming, and expensive. In this paper, simple calibration procedures
for rotating Lorentz-force flowmeters are presented. These procedures eliminate the need for
expensive equipment, numerical modeling, redundant flowmeters, and system down-time.
The calibration processes are explained in a step-by-step manner and compared to experimental results.
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.
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.
Large-eddy simulations were employed over half-ice and half-water surfaces, with varying surface temperatures, wind speeds, directions, as to test if the atmospheric interaction with the heterogeneous surface can be predicted via a heterogeneity Richardson number. This dataset was used to determine that surface heat fluxes over ice, water, and the aggregate surface seem to be captured reasonably well by the wind direction and the heterogeneity Richardson number, but the mean wind and turbulent kinetic energy (TKE) profiles were not, suggesting that not only the difference in stability between the two surface, but also the individual stabilities over each surface influence the dynamics.
Large-eddy simulations were employed over five different sea ice patterns, with a constant ice fraction, to test if the overlying atmospheric boundary layer (ABL) dynamics and thermodynamics differs. The results of these simulations were used to determine that there were differences in vertical heat flux, momentum flux, and horizontal wind speed, and that more surface information is needed to predict the ABL over the sea ice surface. To see what other surface information is needed, twenty-two landscape metrics were calculated over forty-four different maps at differing resolutions, using the FRAGSTATs program. The results of that analysis are available in a .csv file in this dataset.
Myers, Clayton; Yamada, Masaaki; Ji, Hantao; Yoo, Jongsoo; Jara-Almonte, Jonathan; Fox, William
Abstract:
The loss-of-equilibrium is a solar eruption mechanism whereby a sudden breakdown of the magnetohydrodynamic force balance in the Sun's corona ejects a massive burst of particles and energy into the heliosphere. Predicting a loss-of-equilibrium, which has more recently been formulated as the torus instability, relies on a detailed understanding of the various forces that hold the pre-eruption magnetic flux rope in equilibrium. Traditionally, idealized analytical force expressions are used to derive simplified eruption criteria that can be compared to solar observations and modeling. What is missing, however, is a validation that these idealized analytical force expressions can be applied to the line-tied, low-aspect-ratio conditions of the corona. In this paper, we address this shortcoming by using a laboratory experiment to study the forces that act on long-lived, arched, line-tied magnetic flux ropes. Three key force terms are evaluated over a wide range of experimental conditions: (1) the upward hoop force; (2) the downward strapping force; and (3) the downward toroidal field tension force. First, the laboratory force measurements show that, on average, the three aforementioned force terms cancel to produce a balanced line-tied equilibrium. This finding validates the laboratory force measurement techniques developed here, which were recently used to identify a dynamic toroidal field tension force that can prevent flux rope eruptions [Myers et al., Nature 528, 526 (2015)]. The verification of magnetic force balance also confirms the low-beta assumption that the plasma thermal pressure is negligible in these experiments. Next, the measured force terms are directly compared to their corresponding analytical expressions. While the measured and analytical forces are found to be well correlated, the low-aspect-ratio, line-tied conditions in the experiment are found to both reduce the measured hoop force and increase the measured tension force with respect to analytical expectations. These two co-directed effects combine to generate laboratory flux rope equilibria at lower altitudes than are predicted analytically. Such considerations are expected to modify the loss-of-equilibrium eruption criteria for analogous flux ropes in the solar corona.
Myers, Clayton; Yamada, Masaaki; Ji, Hantao; Yoo, Jongsoo; Jara-Almonte, Jonathan; Fox, William
Abstract:
Solar eruptions are often driven by magnetohydrodynamic instabilities such as the torus and kink instabilities that act on line-tied magnetic flux ropes. Recent laboratory experiments designed to study these eruptive instabilities have demonstrated the key role of both dynamic (Myers et al 2015 Nature 528, 526) and quasi-static (Myers et al 2016 Phys. Plasmas, in press) magnetic tension forces in contributing to the equilibrium and stability of line-tied magnetic flux ropes. In this paper, we synthesize these laboratory results and explore the relationship between the dynamic and quasi-static tension forces. While the quasi-static tension force is found to contribute to the flux rope equilibrium in a number of regimes, the dynamic tension force is substantial mostly in the so-called failed torus regime where magnetic self-organization events prevent the flux rope from erupting.
One aspect of the interaction between fast ions and magnetohydrodynamic (MHD) instabilities is the fast ion transport. Coupled kink and tearing MHD instabilities have also been reported to cause fast ion transport. Recently, the ''kick" model has been developed to compute the evolution of the fast ion distribution from the neutral beam injection using instabilities as phase space resonance sources. The goal of this paper is to utilize the kick model to understand the physics of fast ion transport caused by the coupled kink and tearing modes. Soft X-ray diagnostics are used to identify the mode parameters in NSTX. The comparison of neutron rates measured and computed from time-dependent TRANSP simulation with the kick model shows the coupling of kink and tearing mode is important in determination of the fast ion transport. The numerical scan of the mode parameters shows that the relative phase of the kink and tearing modes and the overlapping of kink and tearing mode resonances in the phase space can affect the fast ion transport, suggesting that the synergy of the coupled modes may be causing the fast ion transpor
The growth of magnetic islands in NSTX is modeled successfully, with the consideration of passing fast ions. It is shown that a good quantitative agreement between simulation and experimental measurement can be achieved when the uncompensated cross-field current induced by passing fast ions is included in the island growth model. The fast ion parameters,
along with other equilibrium parameters, are obtained self-consistently using the TRANSP code with the assumptions of the ‘kick’ model (Podestà et al 2017 Plasma Phys. Control. Fusion 59 095008). The results show that fast ions can contribute to overcoming the stabilizing effect of polarization current for magnetic island growth.
Particle distribution functions evolving under the Lorentz operator can be simulated with the Langevin equation for pitch angle scattering. This approach is frequently used in particle based Monte-Carlo simulations of plasma collisions, among others. However, most numerical treatments do not guarantee energy conservation, which may lead to unphysical artifacts such as numerical heating and spectra distortions. We present a novel structure-preserving numerical algorithm for the Langevin equation for pitch angle scattering. Similar to the well-known Boris algorithm, the proposed numerical scheme takes advantage of the structure-preserving properties of the Cayley transform when calculating the velocity-space rotations. The resulting algorithm is explicitly solvable, while preserving the norm of velocities down to machine precision. We demonstrate that the method has the same order of numerical convergence as the traditional stochastic Euler-Maruyama method.
