Search

Start Over Year Published 2017

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

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

2. Overview of NSTX Upgrade Initial Results and Modelling Highlights

Author(s):
Menard, J.E.; Allain, J.P.; Battaglia, D.J.; Bedoya, F.; Bell, R.E.; Belova, E.; Berkery, J.W.; Boyer, M.D.; Crocker, N.; Diallo, A.; Ebrahimi, F.; Ferrraro, N.; Fredrickson, E.; Frerichs, H.; Gerhardt, S.; Gorelenkov, N.; Guttenfelder, W.; Heidbrink, W.; Kaita, R.; Kaye, S.M.; Kriete, D.M.; Kubota, S.; LeBlanc, B.P.; Liu, D.; Lunsford, R.; Mueller, D.; Myers, C.E.; Ono, M.; Park, J.-K.; Podesta, M.; Raman, R.; Reinke, M.; Ren, Y.; Sabbagh, S.A.; Schmitz, O.; Scotti, F.; Sechrest, Y.; Skinner, C.H.; Smith, D.R.; Soukhanovskii, V.; Stoltzfus-Dueck, T.; Yuh, H.; Wang, Z.; Waters, I.; Ahn, J.-W.; Andre, R.; Barchfeld, R.; Beiersdorfer, P.; Bertelli, N.; Bhattacharjee, A.; Boyle, D.; Brennan, D.; Buttery, R.; Capece, A.; Canal, G.; Canik, J.; Chang, C.S.; Darrow, D.; Delgado-Aparicio, L.; Domier, C.; Ethier, S.; Evans, T.; Ferron, J.; Finkenthal, M.; Fonck, R.; Gan, K.; Gates, D.; Goumiri, I.; Gray, T.; Hosea, J.; Humphreys, D.; Jarboe, T.; Jardin, S.; Jaworski, M.A.; Koel, B.; Kolemen, E.; Ku, S.; LaHaye, R.J.; Levinton, F.; Luhmann Jr., N.; Maingi, R.; Maqueda, R.; McKee, G.; Meier, E.; Myra, J.; Perkins, R.; Poli, F.; Rhodes, T.; Riquezes, J.; Rowley, C.; Russell, D.; Schuster, E.; Stratton, B.; Stutman, D.; Taylor, G.; Tritz, K.; Wang, W.; Wirth, B.; Zweben, S.J.
Abstract:
The National Spherical Torus Experiment (NSTX) has undergone a major upgrade, and the NSTX Upgrade (NSTX-U) Project was completed in the summer of 2015. NSTX-U first plasma was subsequently achieved, diagnostic and control systems have been commissioned, H-Mode accessed, magnetic error fields identified and mitigated, and the first physics research campaign carried out. During 10 run weeks of operation, NSTX-U surpassed NSTX-record pulse-durations and toroidal fields, and high-performance ~1MA H-mode plasmas comparable to the best of NSTX have been sustained near and slightly above the n=1 no-wall stability limit and with H-mode confinement multiplier H98y2 above 1. Transport and turbulence studies in L-mode plasmas have identified the coexistence of at least two ion-gyro-scale turbulent micro-instabilities near the same radial location but propagating in opposite (i.e. ion and electron diamagnetic) directions. These modes have the characteristics of ion-temperature gradient and micro-tearing modes, respectively, and the role of these modes in contributing to thermal transport is under active investigation. The new second more tangential neutral beam injection was observed to significantly modify the stability of two types of Alfven Eigenmodes. Improvements in offline disruption forecasting were made in the areas of identification of rotating MHD modes and other macroscopic instabilities using the Disruption Event Characterization and Forecasting (DECAF) code. Lastly, the Materials Analysis and Particle Probe (MAPP) was utilized on NSTX-U for the first time and enabled assessments of the correlation between boronized wall conditions and plasma performance. These and other highlights from the first run campaign of NSTX-U are described.
Type:
Dataset
Issue Date:
October 2017

3. Advances in boronization on NSTX-Upgrade

Author(s):
Skinner, C.H.; Bedoya, F.; Scotti, F.; Allain, J.P.; Blanchard, W.; Cai, D.; Jaworski, M.; Koel, B.E.
Abstract:
Boronization has been effective in reducing plasma impurities and enabling access to higher density, higher confinement plasmas in many magnetic fusion devices. The National Spherical Torus eXperiment, NSTX, has recently undergone a major upgrade to NSTX-U in order to develop the physics basis for a ST-based Fusion Nuclear Science Facility (FNSF) with capability for double the toroidal field, plasma current, and NBI heating power and increased pulse duration from 1–1.5 s to 5–8 s. A new deuterated tri-methyl boron conditioning system was implemented together with a novel surface analysis diagnostic. We report on the spatial distribution of the boron deposition versus discharge pressure, gas injection and electrode location. The oxygen concentration of the plasma facing surface was measured by in-vacuo XPS and increased both with plasma exposure and with exposure to trace residual gases. This increase correlated with the rise of oxygen emission from the plasma.
Type:
Dataset
Issue Date:
January 2017

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

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

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

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

6. Stabilizing effect of resistivity towards ELM-free H-mode discharge in lithium-conditioned NSTX

Author(s):
Banerjee, D.; Zhu, P.; Maingi, R.
Abstract:
Linear stability analysis of the national spherical torus experiment (NSTX) Li-conditioned ELM-free H-mode equilibria is carried out in the context of the extended magneto-hydrodynamic (MHD) model in NIMROD. The purpose is to investigate the physical cause behind edge localized mode (ELM) suppression in experiment after the Li-coating of the divertor and the first wall of the NSTX tokamak. Besides ideal MHD modeling, including finite-Larmor radius effect and two-fluid Hall and electron diamagnetic drift contributions, a non-ideal resistivity model is employed, taking into account the increase of Z eff after Li-conditioning in ELM-free H-mode. Unlike an earlier conclusion from an eigenvalue code analysis of these equilibria, NIMROD results find that after reduced recycling from divertor plates, profile modification is necessary but insufficient to explain the mechanism behind complete ELMs suppression in ideal two-fluid MHD. After considering the higher plasma resistivity due to higher Z eff , the complete stabilization could be explained. A thorough analysis of both pre-lithium ELMy and with-lithium ELM-free cases using ideal and non-ideal MHD models is presented, after accurately including a vacuum-like cold halo region in NIMROD to investigate ELMs.
Type:
Dataset
Issue Date:
July 2017

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

Author(s):
Diallo, A.; Banerjee, S.; Zweben, S.; Stoltzfus-Dueck, T.
Abstract:
We studied the energy exchange dynamics across the low-to-high-confinement (L-H) in NSTX discharges using the gas-puff imaging (GPI) diagnotic. The investigation focused on the energy exchange between flows and turbulence, to help clarify the mechanism of the L-H transition. We apply this study to three type of heating schemes, including a total of 17 shots from the NSTX 2010 campaign run. Results show that the edge fluctuation characteristics (fluctuation levels, radial and poloidal correlation lengths) measured using GPI do not vary just prior to the H-mode transition, but change after the transition. Using a velocimetry approach (orthogonal-programming decomposition), velocity fields of a 24 $\times$ 30 cm GPI view during the L-H transition were obtained with good spatial ($\sim$1 cm) and temporal ($\sim$2.5 $\mu$s) resolutions. Analysis using these velocity fields shows that the production term is systematically negative just prior to the L-H transition indicating transfer from mean flows to turbulence, which is inconsistent with the predator-prey paradigm. Moreover, using the inferred absolute value of the production term, an estimate of the L-H transition duration is found to be 25 ms, which is much larger than the measured duration. These discrepancies are further reinforced by consideration of the ratio between the kinetic energy in the mean flow to the thermal free energy, which is estimated to be much less than 1, suggesting again that turbulence depletion mechanism may not be playing an important role in the transition to the H-mode. Although the Reynolds work is too small to directly deplete the turbulent free energy reservoir, order-of-magnitude analysis shows that the Reynolds stress may still make a non-negligible contribution to the observed poloidal flows.
Type:
Dataset
Issue Date:
June 2017

8. Kinetic neoclassical calculations of impurity radiation profiles

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

9. Suppression of Alfvén modes on NSTX-U with outboard beam injection

Author(s):
Fredrickson, E.D.; Belova, E.V.; Battaglia, D.J.; Bell, R.E.; Crocker, N.A.; Darrow, D.S.; Diallo, A.; Gerhardt, S.P.; Gorelenkov, N.N.; LeBlanc, B.P.; Podesta, M.
Abstract:
In this paper we present data from experiments on NSTX-U where it is shown for the first time that small amounts of high pitch-angle beam ions can strongly suppress the counter-propagating Global Alfvén Eigenmodes (GAE). GAE have been implicated in the redistribution of fast ions and modification of the electron power balance in previous experiments on NSTX. The ability to predict the stability of Alfvén modes, and developing methods to control them, is important for fusion reactor like the International Tokamak Experimental Reactor (ITER) which are heated by a large population of non-thermal, super-Alfvénic ions consisting of fusion generated alphas and beam ions injected for current profile control. We present a qualitative interpretation of these observations using an analytic model of the Doppler-shifted ion-cyclotron resonance drive responsible for GAE instability which has an important dependence on k⊥ρL. A quantitative analysis of this data with the HYM stability code predicts both the frequencies and instability of the GAE prior to, and suppression of the GAE after the injection of high pitch-angle beam ions.
Type:
Dataset
Issue Date:
June 2017

10. Plasma measurements of the Fe XVII L-shell emission and blending with F VIII and F IX

Author(s):
Beiersdorfer, P.; Lepson, J.K.; Gu, M.F.; Bitter, M.
Abstract:
We measured the L-shell emission spectrum of Fe XVII in a low-density, low-gradient magnetically confined laboratory plasma that contains predominantly C, O, Fe, and Ni as trace elements and find excellent agreement with the relative spectral emission obtained in solar and astrophysical observations. However, we obtained spectra that appear to have an usually large 1s^22s^22p^5_{1/2}3d_{3/2} --> 1s^22s^22p^6 Fe XVII resonance transition, commonly labeled 3C, from hot plasmas that also contain F. The wavelength of the Ly-alpha feature of F IX is coincident with the wavelength of the Fe XVII line 3C within one part in 538, and its flux, therefore, enhances the Fe XVII resonance line. Moreover, the resonance and forbidden lines of F VIII are close to the 3s --> 2p transitions in Fe XVII, and may further alter the inferred apparent Fe XVII line ratios, particularly in spectrometers with moderate spectral resolution. The enhanced emission of line 3C, thus, can serve as a new spectral diagnostic for the detection of fluorine in astrophysical plasmas.
Type:
Dataset
Issue Date:
November 2017

11. Noise correlations in the human brain and their impact on pattern classification

Author(s):
Bejjanki, Vikranth R.; da Silveira, Rava Azeredo; Cohen, Jonathan D.; Turk-Browne, Nicholas B.
Abstract:
Multivariate decoding methods, such as multivoxel pattern analysis (MVPA), are highly effective at extracting information from brain imaging data. Yet, the precise nature of the information that MVPA draws upon remains controversial. Most current theories emphasize the enhanced sensitivity imparted by aggregating across voxels that have mixed and weak selectivity. However, beyond the selectivity of individual voxels, neural variability is correlated across voxels, and such noise correlations may contribute importantly to accurate decoding. Indeed, a recent computational theory proposed that noise correlations enhance multivariate decoding from heterogeneous neural populations. Here we extend this theory from the scale of neurons to functional magnetic resonance imaging (fMRI) and show that noise correlations between heterogeneous populations of voxels (i.e., voxels selective for different stimulus variables) contribute to the success of MVPA. Specifically, decoding performance is enhanced when voxels with high vs. low noise correlations (measured during rest or in the background of the task) are selected during classifier training. Conversely, voxels that are strongly selective for one class in a GLM or that receive high classification weights in MVPA tend to exhibit high noise correlations with voxels selective for the other class being discriminated against. Furthermore, we use simulations to show that this is a general property of fMRI data and that selectivity and noise correlations can have distinguishable influences on decoding. Taken together, our findings demonstrate that if there is signal in the data, the resulting above-chance classification accuracy is modulated by the magnitude of noise correlations.
Type:
Dataset
Issue Date:
August 2017

12. A Reduced Resistive Wall Mode Kinetic Stability Model for Disruption Forecasting

Author(s):
Berkery, J.W.; Sabbagh, S.A.; Bell, R.E.; Gerhardt, S.P.; LeBlanc, B.P.
Abstract:
Kinetic modification of ideal stability theory from stabilizing resonances of mode-particle interaction has had success in explaining resistive wall mode (RWM) stability limits in tokamaks. With the goal of real-time stability forecasting, a reduced kinetic stability model has been implemented in the new Disruption Event Characterization and Forecasting (DECAF) code, which has been written to analyze disruptions in tokamaks. The reduced model incorporates parameterized models for ideal limits on beta, a ratio of plasma pressure to magnetic pressure, which are shown to be in good agreement with DCON code calculations. Increased beta between these ideal limits causes a shift in the unstable region of delta W_K space, where delta W_K is the change in potential energy due to kinetic effects that is solved for by the reduced model, such that it is possible for plasmas to be unstable at intermediate beta but stable at higher beta. Gaussian functions for delta W_K are defined as functions of E cross B frequency and collisionality, with parameters reflecting the experience of the National Spherical Torus Experiment (NSTX). The reduced model was tested on a database of discharges from NSTX and experimentally stable and unstable discharges were separated noticeably on a stability map in E cross B frequency, collisionality space. The reduced model only failed to predict an unstable RWM in 15.6% of cases with an experimentally unstable RWM and performed well on predicting stability for experimentally stable discharges as well.
Type:
Dataset
Issue Date:
May 2017

13. The Role of Recti ed Currents in Far-Field RF Sheaths and in SOL Losses of HHFW Power on NSTX

Author(s):
Perkins, R.J.; Hosea, J.C.; Jaworski, M.A.; Bell, R.E.; Bertelli, N.; Kramer, G.J.; Roquemore, L.; Taylor, G.; Wilson, J.R.
Abstract:
Radio-frequency (RF) rectification is an important sheath phenomenon for wave heating of plasma in fusion devices and is proposed to be the mechanism responsible for converting highharmonic fast-wave (HHFW) power in the National Spherical Torus eXperiment (NSTX) into a heat ux to the divertor. RF rectification has two aspects: current rectification and voltage recti- fication, and, while the latter is emphasized in many application, we demonstrate the importance of current rectification in analysis of the NSTX divertor during HHFW heating. When rectified currents are accounted for in first-principle models for the heat ux to the tiles, we predict a sizeable enhancement for the heat ux in the presence of an RF field: for one case studied, the predicted heat ux increases from 0:103 MW=m2 to 0:209 MW=m2. We also demonstrate how this rectification scales with injected HHFW power by tracking probe characteristics during a HHFW power ramp; the rectified current may be clamped at a certain level. This work is important for minimizing SOL losses of HHFW power in NSTX-U but may also have implications for near-field studies of ICRF antennae: ignoring rectified current may lead to underestimated heat uxes and overestimated rectified voltages.
Type:
Dataset
Issue Date:
May 2017

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

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

15. Energetic-particle-modified global Alfven eigenmodes

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

16. Nonlinear simulations of beam-driven Compressional Alfvén Eigenmodes in NSTX

Author(s):
Belova, E.V.; Gorelenkov, N.N.; Crocker, N.A.; Lestz, J.B.; Fredrickson, E.D.; Tang, S.; Tritz, K.
Abstract:
Results of 3D nonlinear simulations of neutral-beam-driven compressional Alfven eigenmodes (CAEs) in the National Spherical Torus Experiment (NSTX) are presented. Hybrid MHD-particle simulations for the H-mode NSTX discharge (shot 141398) using the HYM code show unstable CAE modes for a range of toroidal mode numbers, n=4-9, and frequencies below the ion cyclotron frequency. It is found that the essential feature of CAEs is their coupling to kinetic Alfven wave (KAW) that occurs on the high-field side at the Alfven resonance location. High-frequency Alfven eigenmodes are frequently observed in beam-heated NSTX plasmas, and have been linked to flattening of the electron temperature profiles at high beam power. Coupling between CAE and KAW suggests an energy channeling mechanism to explain these observations, in which beam-driven CAEs dissipate their energy at the resonance location, therefore significantly modifying the energy deposition profile. Nonlinear simulations demonstrate that CAEs can channel the energy of the beam ions from the injection region near the magnetic axis to the location of the resonant mode conversion at the edge of the beam density profile. A set of nonlinear simulations show that the CAE instability saturates due to nonlinear particle trapping, and a large fraction of beam energy can be transferred to several unstable CAEs of relatively large amplitudes and absorbed at the resonant location. Absorption rate shows a strong scaling with the beam power.
Type:
Dataset
Issue Date:
April 2017

17. Application of Benchmarked Kinetic Resistive Wall Mode Stability Codes to ITER, Including Additional Physics

Author(s):
Berkery, J.W.; Wang, Z.R.; Sabbagh, S.A.; Liu, Y.Q.; Betti, R.; Guazotto, L.
Abstract:
Leading resistive wall mode (RWM) stability codes MARS-K [Y. Liu, et al., Phys. Plasmas 15, 112503 (2008)] and MISK [B. Hu, et al., Phys. Plasmas 12, 057301 (2005)] have been previously benchmarked. The benchmarking has now been extended to include additional physics, and used to project the stability of ITER in a realistic operating space. Due to ITER's relatively low plasma rotation and collisionality, collisions and non-resonance rotational effects were both found to have little impact on stability, and these non-resonance rotational effects also will not self-consistently affect the ITER RWM eigenfunction. Resonances between thermal ions and electrons and the expected level of ITER toroidal rotation were found to be important to stability, as were alpha particles, which are not in rotational resonance. MISK calculations show that without alpha particles, ITER is projected to be unstable to the RWM, but the expected level of alphas is calculated to provide a sufficient level of stability.
Type:
Dataset
Issue Date:
October 2017

18. The Structured `Low Temperature' Phase of the Retinal Population Code

Author(s):
Ioffe, Mark Lev; Berry II, Michael J.
Abstract:
Recent advances in experimental techniques have allowed the simultaneous recordings of populations of hundreds of neurons, fostering a debate about the nature of the collective structure of population neural activity. Much of this debate has focused on the empirical findings of a phase transition in the parameter space of maximum entropy models describing the measured neural probability distributions, interpreting this phase transition to indicate a critical tuning of the neural code. Here, we instead focus on the possibility that this is a first-order phase transition which provides evidence that the real neural population is in a `structured', collective state. We show that this collective state is robust to changes in stimulus ensemble and adaptive state. We find that the pattern of pairwise correlations between neurons has a strength that is well within the strongly correlated regime and does not require fine tuning, suggesting that this state is generic for populations of 100+ neurons. We find a clear correspondence between the emergence of a phase transition, and the emergence of attractor-like structure in the inferred energy landscape. A collective state in the neural population, in which neural activity patterns naturally form clusters, provides a consistent interpretation for our results.
Type:
Dataset
Issue Date:
2017

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

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

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

Author(s):
Fil, A.; Kolemen, E.; Bortolon, A.; Ferraro, N.; Jardin, S.; Parks, P.B.; Lunsford, R.; Maingi, R.
Abstract:
In this paper we present initial simulations of pedestal control by Lithium Granule Injection (LGI) in NSTX. A model for small granule ablation has been implemented in the M3D-C1 code [1], allowing the simulation of realistic Lithium granule injections. 2D simulations in NSTX L-mode and H-mode plasmas are done and the effect of granule size, injection angle and velocity on the pedestal gradient increase are studied. For H-mode cases, the amplitude of the local pressure perturbation caused by the granules is highly dependent on the solid granule size. In our simulations, reducing the granule injection velocity allows one to inject more particles at the pedestal top.
Type:
Dataset
Issue Date:
January 2017

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

Author(s):
Lunsford, R.; Bortolon, A.; Roquemore, A.L.; Mansfield, D.K.; Jaworski, M.A.; Kaita, R.; Maingi, R.; Nagy, A.
Abstract:
By employing a neutral gas shielding (NGS) model to characterize impurity granule injection the pedestal atomic deposition for three different species of granule: lithium, boron, and carbon are determined. Utilizing the duration of ablation events recorded on experiments performed at DIII-D to calibrate the NGS model we are able to quantify the ablation rate and mass deposition location with respect to the plasma density profile. The species specific granule shielding constant is then used to model granule ablation within NSTX-U discharges. Simulations of 300, 500 and 700 micron diameter granules injected at 50 m/sec are presented for NSTX-U L-mode type plasmas as well as H-mode discharges with low natural ELM frequencies. Additionally, ablation calculations of 500 micron granules of each species are presented at velocities ranging from 50 � 150 m/sec. In H-mode type discharges these simulations show that the majority of the injected granule is ablated within or just past the steep gradient region of the discharge. At this radial position, the perturbation to the background plasma generated by the ablating granule can lead to conditions advantageous for the rapid triggering of an ELM crash event.
Type:
Dataset
Issue Date:
July 2017

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

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

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

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

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

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

26. Neutral recycling effects on ITG turbulence

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

27. Optimization of the angular orientation for a fast ion loss detector in a tokamak

Author(s):
Darrow, D.
Abstract:
A scintillator type fast ion loss detector measures the gyroradius and pitch angle distribution of superthermal ions escaping from a magnetically confined fusion plasma at a single location. Described here is a technique for optimizing the angular orientation of such a detector in an axisymmetric tokamak geometry in order to intercept losses over a useful and interesting ranges of pitch angle. The method consists of evaluating the detector acceptance as a function of the fast ion constants of motion, i.e. energy, canonical toroidal momentum, and magnetic moment. The detector acceptance can then be plotted in a plane of constant energy and compared with the relevant orbit class boundaries and fast ion source distributions. Knowledge of expected or interesting mechanisms of loss can further guide selection of the detector orientation. The example of a fast ion loss detector for the National Spherical Torus Experiment-Upgrade (NSTX-U) is considered.
Type:
Dataset
Issue Date:
January 2017

28. A multi-machine scaling of halo current rotation

Author(s):
Myers, C.E.; Eidietis, N.W.; Gerasimov, S.N.; Gerhardt, S.Pl.; Granetz, R.S.; Hender, T.C.; Pautasso, G.
Abstract:
Halo currents generated during unmitigated tokamak disruptions are known to develop rotating asymmetric features that are of great concern to ITER because they can dynamically amplify the mechanical stresses on the machine. This paper presents a multi-machine analysis of these phenomena. More specifically, data from C-Mod, NSTX, ASDEX Upgrade, DIII-D, and JET are used to develop empirical scalings of three key quantities: (1) the machine-specific minimum current quench time, tauCQ; (2) the halo current rotation duration, trot; and (3) the average halo current rotation frequency, <fh>. These data reveal that the normalized rotation duration, trot/tauCQ, and the average rotation velocity, <vh>, are surprisingly consistent from machine to machine. Furthermore, comparisons between carbon and metal wall machines show that metal walls have minimal impact on the behavior of rotating halo currents. Finally, upon projecting to ITER, the empirical scalings indicate that substantial halo current rotation above <fh> = 20 Hz is to be expected. More importantly, depending on the projected value of tauCQ in ITER, substantial rotation could also occur in the resonant frequency range of 6-20 Hz. As such, the possibility of damaging halo current rotation during unmitigated disruptions in ITER cannot be ruled out.
Type:
Dataset
Issue Date:
October 2017

29. Application of IR imaging for free-surface velocity measurement in liquid-metal systems

Author(s):
Hvasta, M. G.; Kolemen, E.; Fisher, A.
Abstract:
Measuring free-surface, liquid-metal flow velocity is challenging to do in a reliable and accurate manner. This paper presents a non-invasive, easily-calibrated method of measuring the surface velocities of open-channel liquid-metal flows using an IR camera. Unlike other spatially-limited methods, this IR camera particle tracking technique provides full field-of-view data that can be used to better understand open-channel flows and determine surface boundary conditions. This method could be implemented and automated for a wide range of liquid-metal experiments, even if they operate at high-temperatures or within strong magnetic fields.
Type:
Dataset
Issue Date:
January 2017

30. Application of IR imaging for free-surface velocity measurement in liquid-metal systems

Author(s):
Hvasta, M.H.; Kolemen, E.; Fisher, A.
Abstract:
Measuring free-surface, liquid-metal flow velocity is challenging to do in a reliable and accurate manner. This paper presents a non-invasive, easily-calibrated method of measuring the surface velocities of open-channel liquid-metal flows using an IR camera. Unlike other spatially-limited methods, this IR camera particle tracking technique provides full field-of-view data that can be used to better understand open-channel flows and determine surface boundary conditions. This method could be implemented and automated for a wide range of liquid-metal experiments, even if they operate at high-temperatures or within strong magnetic fields.
Type:
Dataset
Issue Date:
January 2017

31. Experimental calibration procedures for rotating Lorentz-force flowmeters

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

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

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

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

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

34. A Riccati Solution for the Ideal MHD Plasma Response with Applications to Real-time Stability Control

Author(s):
Glasser, A.; Kolemen, E.; Glasser, A.H.
Abstract:
To effectuate near real-time feedback control of ideal MHD instabilities in a tokamak geometry, a rapid solution for stability analysis is a prerequisite. Toward this end, we reformulate the δW stability method with a Hamilton-Jacobi theory, elucidating analytical and numerical features of the generic tokamak ideal MHD stability problem. The plasma response matrix is demonstrated to be the solution of an ideal MHD matrix Riccati differential equation (MRDE). Since Riccati equations are prevalent in the control theory literature, such a shift in perspective brings to bear a range of numerical methods that are well-suited to the robust, fast solution of control problems. We discuss the usefulness of Riccati techniques in solving the stiff ODEs often encountered in ideal MHD stability analyses-—for example, in tokamak edge and stellarator physics. We then demonstrate the applicability of such methods to an existing 2D ideal MHD stability code—DCON—enabling its parallel operation in near real-time. Output is shown to match with high accuracy, and with wall-clock time ≪ 1s. Such speed may help enable active feedback ideal MHD stability control, especially in tokamak plasmas whose ideal MHD equilibria evolve with inductive timescale τ > 1s-—as in ITER.
Type:
Dataset
Issue Date:
March 2018

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

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

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

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

37. Amplification of local changes along the timescale processing hierarchy

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

38. Sensitivity of WallDYN material migration modeling to uncertainties in mixed-material surface binding energies

Author(s):
Nichols, J.H.; Jaworski, M.A.; Schmid, K.
Abstract:
The WallDYN package has recently been applied to a number of tokamaks to self-consistently model the evolution of mixed-material plasma facing surfaces. A key component of the WallDYN model is the concentration-dependent surface sputtering rate, calculated using SDTRIM.SP. This modeled sputtering rate is strongly influenced by the surface binding energies (SBEs) of the constituent materials, which are well known for pure elements but often are poorly constrained for mixed-materials. This work examines the sensitivity of WallDYN surface evolution calculations to different models for mixed-material SBEs, focusing on the carbon/lithium/oxygen/deuterium system present in NSTX. A realistic plasma background is reconstructed from a high density, H-mode NSTX discharge, featuring an attached outer strike point with local density and temperature of 4e20 m^-3 and 4 eV, respectively. It is found that various mixed-material SBE models lead to significant qualitative and quantitative changes in the surface evolution profile at the outer divertor, with the highest leverage parameter being the C-Li binding model. Uncertainties of order 50%, appearing on time scales relevant to tokamak experiments, highlight the importance of choosing an appropriate mixed-material sputtering representation when modeling the surface evolution of plasma facing components. These results are generalized to other fusion-relevant materials with different ranges of SBEs.
Type:
Dataset
Issue Date:
2017

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

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

40. Two-dimensional turbulence cross-correlation functions in the edge of NSTX

Author(s):
Zweben, S.J.; Stotler, D.P.; Scotti, F.; Myra, J.R.
Abstract:
The 2-D radial vs. poloidal cross-correlation functions of edge plasma turbulence were measured near the outer midplane using the gas puff imaging (GPI) diagnostic on NSTX. These correlation functions were evaluated at radii r= 0 cm, ±3 cm, and ±6 cm from the separatrix and poloidal locations p=0 cm and ±7.5 cm from the GPI poloidal center line for 20 different shots. The ellipticity ε and tilt angle φ of the positive cross- correlation regions, and the minimum negative cross-correlation “cmin” and total negative over positive values “neg/pos” were evaluated for each of these cases. The average results over this data set were ε=2.2±0.9, φ=87±34o (i.e. poloidally oriented), cmin= -0.30±0.15, and neg/pos=0.25±0.24. Thus there was significant variation in these correlation results within this database, with dependences on the location within the image, the magnetic geometry, and the plasma parameters. Possible causes for this variation are discussed, including the misalignment of the GPI view with the local B field line, the magnetic shear of field lines in the edge, the poloidal flow shear of the turbulence, blob-hole correlations, and the neutral density ‘shadowing’ effect in GPI.
Type:
Dataset
Issue Date:
September 2017

41. Fracture aperture maps used to study reactive transport, channelization, and permeability evolution in carbonate rocks.

Author(s):
Peters, Catherine
Abstract:
Fractures in geological formations may enable migration of environmentally relevant fluids, as in leakage of CO2 through caprocks in geologic carbon sequestration. We investigated geochemically induced alterations of fracture geometry in Indiana Limestone specimens. Experiments were the first of their kind, with periodic high-resolution imaging using X-ray computed tomography (xCT) scanning while maintaining high pore pressure (100 bar). We studied two CO2-acidified brines having the same pH (3.3) and comparable thermodynamic disequilibrium but different equilibrated pressures of CO2 (PCO2 values of 12 and 77 bar). High-PCO2 brine has a faster calcite dissolution kinetic rate because of the accelerating effect of carbonic acid. Contrary to expectations, dissolution extents were comparable in the two experiments. However, progressive xCT images revealed extensive channelization for high PCO2, explained by strong positive feedback between ongoing flow and reaction. The pronounced channel increasingly directed flow to a small region of the fracture, which explains why the overall dissolution was lower than expected. Despite this, flow simulations revealed large increases in permeability in the high-PCO2 experiment. This study shows that the permeability evolution of dissolving fractures will be larger for faster-reacting fluids. The overall mechanism is not because more rock dissolves, as would be commonly assumed, but because of accelerated fracture channelization.
Type:
Dataset
Issue Date:
December 2017

Limit your search