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.
The ability of an injected lithium granule to promptly trigger an edge localized mode (ELM) has been established in multiple experiments. By horizontally injecting granules ranging in diameter from 200 microns to 1mm in diameter into the low field side of EAST H-mode discharges we have determined that granules with diameter > 600 microns are successful in triggering ELMs more than 95% of the time. It was also demonstrated that below 600 microns the triggering efficiency decreased roughly with granule size. Granules were radially injected from the outer midplane with velocities ~ 80 m/s into EAST upper single null discharges with an ITER like tungsten monoblock divertor. These granules were individually tracked throughout their injection cycle in order to determine their efficacy at triggering an ELM. For those granules of sufficient size, ELM triggering was a prompt response to granule injection. By simulating the granule injection with an experimentally benchmarked neutral gas shielding (NGS) model, the ablatant mass deposition required to promptly trigger an ELM is calculated and the fractional mass deposition is determined.
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.
Monitoring the attention of others is fundamental to social cognition. Most of the literature on the topic assumes that our social cognitive machinery is tuned specifically to the gaze direction of others as a proxy for attention. This standard assumption reduces attention to an externally visible parameter. Here we show that this assumption is wrong and a deeper, more meaningful representation is involved. We presented subjects with two cues about the attentional state of a face: direction of gaze and emotional expression. We tested whether people relied predominantly on one cue, the other, or both. If the traditional view is correct, then the gaze cue should dominate. Instead, people employed a variety of strategies, some relying on gaze, some on expression, and some on an integration of cues. We also assessed people’s social cognitive ability using two, independent, standard tests. If the traditional view is correct, then social cognitive ability, as assessed by the independent tests, should correlate with the degree to which people successfully use the gaze cue to judge the attention state of the face. Instead, social cognitive ability correlated best with the degree to which people successfully integrated the cues together, instead of with the use of any one specific cue. The results suggest a rethink of a fundamental component of social cognition: monitoring the attention of others involves constructing a deep model that is informed by a combination of cues. Attention is a rich process and monitoring the attention of others involves a similarly rich representation.
A matrix inversion technique is derived to calculate local ion temperature from line-integrated measurements of an extended emission source in an axisymmetric plasma which exactly corrects for both toroidal velocity and radial velocity components. Local emissivity and toroidal velocity can be directly recovered from line-integrated spectroscopic measurements, but an independent measurement of the radial velocity is necessary to complete the temperature inversion. The extension of this technique to handle the radial velocity is relevant for magnetic reconnection and merging compression devices where temperature inversion from spectroscopic measurements is desired. A simulation demonstrates the effects of radial velocity on the determination of ion temperature.
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.
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.
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.
Non-axisymmetric control coils and the so-called snowflake divertor configuration are two potential solutions proposed to solve two separate outstanding issues on the path towards self-sustained burning plasma operations, namely the transient energy bursts caused by edge localized modes and the steady state heat exhaust problem. In a reactor, these two proposed solutions would have to operate simultaneously and it is, therefore, important to investigate their compatibility and to identify possible conflicts that could prevent them from operating simultaneously. In this work, single- and two-fluid resistive magnetohydrodynamic calculations are used to investigate the effect of externally applied magnetic perturbations on the snowflake divertor configuration. The calculations are based on simulated NSTX-U plasmas and the results show that additional and longer magnetic lobes are created in the null-point region of the snowflake configuration, compared to those in the conventional single-null. The intersection of these longer and additional lobes with the divertor plates are expected to cause more striations in the particle and heat flux target profiles. In addition, the results indicate that the size of the magnetic lobes, in both single-null and snowflake configurations, are more sensitive to resonant magnetic perturbations than to non-resonant magnetic perturbations. The results also suggest that lower values of current in non-axisymmetric control coils would be required to suppress edge localized modes in plasmas with the snowflake configuration.
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.
This is the dataset for the plots presented in the article "CO2-leakage-driven diffusiophoresis causes spontaneous accumulation of charged materials in channel flow."
Natural gas vehicles (NGVs) have been promoted in China to mitigate air pollution, yet our measurements and analyses show that NGV growth in China may have significant negative impacts on climate change. We conducted real-world vehicle emission measurements in China and found high methane emissions from heavy-duty NGVs (90% higher than current emission limits). These emissions have been ignored in previous emission estimates, leading to biased results. Applying our observations to life-cycle analyses, we found that switching to NGVs from conventional vehicles in China has led to a net increase in greenhouse gas (GHG) emissions since 2000. With scenario analyses, we also show that the next decade will be critical for China to reverse the trend with the upcoming China VI standard for heavy-duty vehicles. Implementing and enforcing the China VI standard is challenging, and the method demonstrated here can provide critical information regarding the fleet-level CH4 emissions from NGVs.
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.
Kramer, G. J.; Bortolon, A.; Ferraro, N. M.; Spong, D. A.; Crocker, N. A.; Darrow, D. S.; Fredrickson, E. D.; Kubota, S.; Park, J.-K.; Podesta, M.; Heidbrink, W. W.
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}.
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.
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
These GROMACS trajectories show the existence of a critical point in deeply supercooled WAIL water. Also included is the code necessary to reproduce the figures in the corresponding paper from these trajectories. From this data the critical temperature, pressure, and density of the model can be found, and critical fluctuations in the deeply supercooled liquid can be directly observed (in a computer-simulation sense).
Movies of relativistic reconnection and particle acceleration in relativistic reconnection accompanying the article "Relativistic Reconnection: an Efficient Source of Nonthermal Particles" by Lorenzo Sironi and Anatoly Spitkovsky.
Magnetic reconnection in partially ionized plasmas is a ubiquitous and important phenomena in both laboratory and astrophysical systems. Here, simulations of partially ionized magnetic reconnection with well-matched initial conditions are performed using both multi-fluid and fully-kinetic approaches. Despite similar initial conditions, the time-dependent evolution differs between the two models. In multi-fluid models, the reconnection rate locally obeys either a decoupled Sweet-Parker scaling, where neutrals are unimportant, or a fully coupled Sweet-Parker scaling, where neutrals and ions are strongly coupled, depending on the resistivity. In contrast, kinetic models show a faster reconnection rate that is proportional to the fully-coupled, bulk Alfv\'en speed, $v_A^\star$. These differences are interpreted as the result of operating in different collisional regimes. Multi-fluid simulations are found to maintain $\nu_{ni}L/v_A^\star \gtrsim 1$, where $\nu_{ni}$ is the neutral-ion collision frequency and $L$ is the time-dependent current sheet half-length. This strongly couples neutrals to the reconnection outflow, while kinetic simulations evolve to allow $\nu_{ni}L/v_A^\star < 1$, decoupling neutrals from the reconnection outflow. Differences in the way reconnection is triggered may explain these discrepancies.
