Search

1. 2-km resolution oversampled NH3 data

Author(s):
Wang, Rui; Guo, Xuehui; Pan, Da; Kelly, James; Bash, Jesse; Sun, Kang; Paulot, Fabien; Clarisse, Lieven; Van Damme, Martin; Whitburn, Simon; Coheur, Pierre-François; Clerbaux, Cathy; Zondlo, Mark
Abstract:
Monthly, high resolution (~2 km) ammonia (NH3) column maps from the Infrared Atmospheric Sounding Interferometer (IASI) were developed across the contiguous United States and adjacent areas. Ammonia hotspots (95th percentile of the column distribution) were highly localized with a characteristic length scale of 12 km and median area of 152 km2. Five seasonality classes were identified with k-means++ clustering. The Midwest and eastern United States had a broad, spring maximum of NH3 (67% of hotspots in this cluster). The western United States, in contrast, showed a narrower mid-summer peak (32% of hotspots). IASI spatiotemporal clustering was consistent with those from the Ammonia Monitoring Network. CMAQ and GFDL-AM3 modeled NH3 columns have some success replicating the seasonal patterns but did not capture the regional differences. The high spatial-resolution monthly NH3 maps serve as a constraint for model simulations and as a guide for the placement of future, ground-based network sites.
Type:
Dataset
Issue Date:
December 2020

3. 60 degree Inclined Pulsar

Abstract:
Current sheet and open field lines with footpoints near the edge of the polar cap. The magnetic axis is inclined relative to the rotation axis by 60 degrees. Red field lines originate on the north polar cap and green field lines in the right panel originate on the south polar cap. Purple and grey colors indicate positive and negative net local charge density in the current sheet, which is shown between 1.2-2 light cylinder radii.
Type:
moving image
Issue Date:
21 November 2013

4. 90 degree Inclined Pulsar

Abstract:
Current sheet and open field lines with footpoints near the edge of the polar cap. The magnetic axis is inclined relative to the rotation axis by 90 degrees. Red field lines originate on the north polar cap and green field lines in the right panel originate on the south polar cap. Purple and grey colors indicate positive and negative net local charge density in the current sheet, which is shown between 1.2-2 light cylinder radii.
Type:
moving image
Issue Date:
21 November 2013

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

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

7. A dual-mechanism antibiotic kills Gram-negative bacteria and avoids drug resistance

Author(s):
Martin, James K; Sheehan, Joseph P; Bratton, Benjamin P; Moore, Gabriel M; Mateus, André; Li, Sophia Hsin-Jung; Kim, Hahn; Rabinowitz, Joshua D; Typas, Athanasios; Savitski, Mikhail M; Wilson, Maxwell Z; Gitai, Zemer
Abstract:
The rise of antibiotic resistance and declining discovery of new antibiotics have created a global health crisis. Of particular concern, no new antibiotic classes have been approved for treating Gram-negative pathogens in decades. Here, we characterize a compound, SCH-79797, that kills both Gram-negative and Gram-positive bacteria through a unique dual-targeting mechanism of action (MoA) with undetectably-low resistance frequencies. To characterize its MoA, we combined quantitative imaging, proteomic, genetic, metabolomic, and cell-based assays. This pipeline demonstrates that SCH-79797 has two independent cellular targets, folate metabolism and bacterial membrane integrity, and outperforms combination treatments in killing MRSA persisters. Building on the molecular core of SCH-79797, we developed a derivative, Irresistin-16, with increased potency and showed its efficacy against Neisseria gonorrheae in a mouse vaginal infection model. This promising antibiotic lead suggests that combining multiple MoAs onto a single chemical scaffold may be an underappreciated approach to targeting challenging bacterial pathogens.
Type:
Dataset
Issue Date:
20 May 2020

8. A dynamic magnetic tension force as the cause of failed solar eruptions

Abstract:
Coronal mass ejections are solar eruptions driven by a sudden release of magnetic energy stored in the Sun’s corona. In many cases, this magnetic energy is stored in long-lived, arched structures called magnetic flux ropes. When a flux rope destabilizes, it can either erupt and produce a coronal mass ejection or fail and collapse back towards the Sun. The prevailing belief is that the outcome of a given event is determined by a magnetohydrodynamic force imbalance called the torus instability. This belief is challenged, however, by observations indicating that torus-unstable flux ropes sometimes fail to erupt. This contradiction has not yet been resolved because of a lack of coronal magnetic field measurements and the limitations of idealized numerical modelling. Here we report the results of a laboratory experiment that reveal a previously unknown eruption criterion below which torus-unstable flux ropes fail to erupt. We find that such ‘failed torus’ events occur when the guide magnetic field (that is, the ambient field that runs toroidally along the flux rope) is strong enough to prevent the flux rope from kinking. Under these conditions, the guide field interacts with electric currents in the flux rope to produce a dynamic toroidal field tension force that halts the eruption. This magnetic tension force is missing from existing eruption models, which is why such models cannot explain or predict failed torus events.
Type:
Dataset
Issue Date:
11 December 2015

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

10. A novel scheme for error field correction in permanent magnet stellarators

Author(s):
Rutkowski, Adam; Hammond, Kenneth; Zhu, Caoxiang; Gates, David; Chambliss, Amelia
Abstract:
Stellarators offer a promising path towards fusion reactors, but their design and construction are complicated by stringent tolerance requirements on highly complex 3D coils. A potential way to simplify the engineering requirements for stellarators is to use simple planar toroidal field coils along with permanent magnet arrays to generate shaping fields. In order to ensure sufficient field accuracy while minimizing engineering complexity and system cost, new techniques are required to correct the field produced by the permanent magnet arrays to within requirements set by plasma physics. This work describes a novel correction method developed for this purpose. This analysis is applied to the design of a quasi-axisymmetric stellarator that employs a combination of permanent magnets and planar toroidal field coils to generate its magnetic field. Analysis techniques and initial results using the method for error correction on a proposed permanent magnet stellarator are shown, and it is demonstrated that the method successfully meets the design requirements of the project.
Type:
Dataset
Issue Date:
7 December 2022

11. A scalable real-time framework for Thomson scattering analysis: Application to NSTX-U

Author(s):
F. M. Laggner, A. Diallo, B. P. LeBlanc, R. Rozenblat, G. Tchilinguirian, E.Kolemen, the NSTX-U team
Abstract:
A detailed description of a prototype setup for real-time (rt) Thomson scattering (TS) analysis is presented and implemented in the multi-point Thomson scattering (MPTS) diagnostic system at the National Spherical Torus Experiment Upgrade(NSTX-U). The data acquisition hardware was upgraded with rt capable electronics (rt-analog digital converters (ADCs) and a rt server) that allow for fast digitization of the laser pulse signal of eight radial MPTS channels. In addition, a new TS spectrum analysis software for a rapid calculation of electron temperature (Te) and electron density (ne) was developed. Testing of the rt hardware and data analysis soft-ware was successfully completed and benchmarked against the standard, post-shot evaluation. Timing tests were performed showing that the end-to-end processing time was reproducibly below 17 ms for the duration of at least 5 s, meeting a 60 Hz deadline by the laser pulse repetition rate over the length of a NSTX-U discharge. The presented rt framework is designed to be scalable in system size, i.e. incorporation of additional radial channels by solely adding additional rt capable hardware. Furthermore, it is scalable in its operation duration and was continuously run for up to 30 min, making it an attractive solution for machines with long discharge duration such as advanced, non-inductive tokamaks or stellarators.
Type:
Dataset
Issue Date:
March 2019

12. A software package for plasma facing component analysis and design: the Heat flux Engineering Analysis Toolkit (HEAT)

Author(s):
Looby, Tom; Reinke, Matthew; Wingen, Andreas; Menard, Jonathan; Gerhardt, Stefan; Gray, Travis; Donovan, David; Unterberg, Ezekial; Klabacha, Jonathan; Messineo, Mike
Abstract:
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.
Type:
Dataset
Issue Date:
March 2021

13. A thermodynamic phase transition in magnetic reconnection

Author(s):
Jara-Almonte, Jonathan; Hantao, Ji
Abstract:
Data supporting the manuscript "A thermodynamic phase transition in magnetic reconnection" published in Physical Review Letters.
Type:
Dataset
Issue Date:
7 July 2021

14. Accessing Real-Life Episodic Information from Minutes versus Hours Earlier Modulates Hippocampal and High-Order Cortical Dynamics

Author(s):
Chen, Janice; Honey, Christopher; Simony, Erez; Arcaro, Michael; Norman, Kenneth; Hasson, Uri
Abstract:
It is well known that formation of new episodic memories depends on hippocampus, but in real-life settings (e.g., conversation), hippocampal amnesics can utilize information from several minutes earlier. What neural systems outside hippocampus might support this minutes-long retention? In this study, subjects viewed an audiovisual movie continuously for 25 min; another group viewed the movie in 2 parts separated by a 1-day delay. Understanding Part 2 depended on retrieving information from Part 1, and thus hippocampus was required in the day-delay condition. But is hippocampus equally recruited to access the same information from minutes earlier? We show that accessing memories from a few minutes prior elicited less interaction between hippocampus and default mode network (DMN) cortical regions than accessing day-old memories of identical events, suggesting that recent information was available with less reliance on hippocampal retrieval. Moreover, the 2 groups evinced reliable but distinct DMN activity timecourses, reflecting differences in information carried in these regions when Part 1 was recent versus distant. The timecourses converged after 4 min, suggesting a time frame over which the continuous-viewing group may have relied less on hippocampal retrieval. We propose that cortical default mode regions can intrinsically retain real-life episodic information for several minutes.
Type:
Dataset
Issue Date:
3 August 2015

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

16. Aligned Nonspinning Dipoles

Abstract:
Magnetic field lines and current sheets for an orbiting neutron star binary with the magnetic moments of both stars aligned with the rotation axis. The stars are not spinning, i.e., R_{LC,∗} = ∞. Fields are by and large confined to the half of the magnetosphere closer to their source star. This movie shows the corotating field pattern as the orbit progresses.
Type:
moving image
Issue Date:
21 November 2013

17. Aligned Rapidly Spinning Dipoles

Abstract:
Magnetic field lines and current sheets for an orbiting neutron star binary with the magnetic moments of both stars aligned with the rotation axis. The stars are spinning rapidly at ∼ ms periods, with R_{LC,∗}/R_∗ = 2.7. Stellar spin winds fields backwards toroidally, and they can propagate to the far side of the magnetosphere closer to the opposing star. This movie shows the corotating field pattern as the orbit progresses.
Type:
moving image
Issue Date:
21 November 2013

18. Alternative energy vehicle deployment coupled with decarbonizing electricity generation amplifies climate, air quality and health co-benefits in China

Author(s):
Peng, Liqun; Liu, Feiqi; Zhou, Mi; Li, Mingwei; Zhang, Qiang; Mauzerall, Denise L.
Abstract:
China is the world's largest carbon emitter and suffers from severe air pollution. About one million deaths in China were attributable to air pollution in 2017. Alternative energy vehicles (AEVs), e.g. electric, hydrogen fuel cell, and natural gas vehicles, can help achieve both carbon emission mitigation and air quality improvement. However, climate, air quality and health co-benefit of AEVs powered by deeply decarbonized electricity generation remain poorly quantified. Here, we conduct a quantitative integrated assessment of the air quality, health, carbon emission mitigation and economic benefits of AEV deployment as the electricity grid decarbonizes in China. We find population-weighted annual PM2.5 and summer O3 concentration can decrease as large as 5.7μgm−3 and 4.9ppb. Annual avoided premature mortalities and years of life lost resulting from improved ambient air pollution can be as large as ~329,000 persons and ~1,611,000 years. We thus show that maximizing climate, air quality and health benefits of AEV deployment in China requires rapid decarbonization of the power system.
Type:
Dataset
Issue Date:
17 June 2021

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

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

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

21. Analytic stability boundaries for compressional and global Alfven eigenmodes driven by fast ions. I. Interaction via ordinary and anomalous cyclotron resonances.

Author(s):
Lestz J.B., Gorelenkov N.N., Belova E.V., Tang S.X., Crocker N.A.
Abstract:
Conditions for net fast ion drive are derived for beam-driven, sub-cyclotron compressional (CAE) and global (GAE) Alfven eigenmodes, such as those routinely observed in spherical tokamaks such as NSTX(-U) and MAST. Both co- and counter-propagating CAEs and GAEs are investigated, driven by the ordinary and anomalous Doppler-shifted cyclotron resonance with fast ions. Whereas prior results were restricted to vanishingly narrow distributions in velocity space, broad parameter regimes are identified in this work which enable an analytic treatment for realistic fast ion distributions generated by neutral beam injection. The simple, approximate conditions derived in these regimes for beam distributions of realistic width compare well to the numerical evaluation of the full analytic expressions for fast ion drive. Moreover, previous results in the very narrow beam case are corrected and generalized to retain all terms in omega/omega_{ci} and k_{||}/kperp, which are often assumed to be small parameters but can significantly modify the conditions of drive and damping when they are non-negligible. Favorable agreement is demonstrated between the approximate stability criterion, simulation results, and a large database of NSTX observations of cntr-GAEs.
Type:
Dataset
Issue Date:
September 2019

22. Analytic stability boundaries for compressional and global Alfven eigenmodes driven by fast ions. II. Interaction via Landau resonance.

Author(s):
Lestz, J.B.; Gorelenkov, N.N.; Belova, E.V.; Tang, S.X.; Crocker, N.A.
Abstract:
Conditions for net fast ion drive are derived for beam-driven, co-propagating, sub-cyclotron compressional (CAE) and global (GAE) Alfven eigenmodes driven by the Landau resonance with super-Alfvenic fast ions. Approximations applicable to realistic neutral beam distributions and mode characteristics observed in spherical tokamaks enable the derivation of marginal stability conditions for these modes. Such conditions successfully reproduce the stability boundaries found from numerical integration of the exact expression for local fast ion drive/damping. Coupling between the CAE and GAE branches of the dispersion due to finite \omega/\omega_{ci} and k_\parallel/k_\perp is retained and found to be responsible for the existence of the GAE instability via this resonance. Encouraging agreement is demonstrated between the approximate stability criterion, simulation results, and a database of NSTX observations of co-CAEs.
Type:
Dataset
Issue Date:
January 2020

23. Antialigned Nonspinning Dipoles

Abstract:
Magnetic field lines and current sheets for an orbiting neutron star binary with the magnetic moment of one star aligned with the rotation axis, and the magnetic moment of the other star tilted and antialigned with the rotation axis. The stars are not spinning, i.e., R_{LC,∗} = ∞. Fields from each star encircle the other star and force fields coming off the second star backwards toroidally. This movie shows the corotating field pattern as the orbit progresses.
Type:
moving image
Issue Date:
21 November 2013

24. Antialigned Rapidly Spinning Dipoles

Abstract:
Magnetic field lines and current sheets for an orbiting neutron star binary with the magnetic moment of one star aligned with the rotation axis, and the magnetic moment of the other star tilted and antialigned with the rotation axis. The stars are spinning rapidly at ∼ ms periods, with R_{LC,∗} /R_∗ = 2.7. Stellar spin winds fields backwards toroidally. This movie shows the corotating field pattern as the orbit progresses.
Type:
moving image
Issue Date:
21 November 2013

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

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

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

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

29. Application of transient CHI plasma startup to future ST and AT devices

Author(s):
Hammond, K.C.; Raman, R.; Jardin, S.C.
Abstract:
Employment of non-inductive plasma start-up techniques would considerably simplify the design of a spherical tokamak fusion reactor. Transient coaxial helicity injection (CHI) is a promising method, expected to scale favorably to next-step reactors. However, the implications of reactor-relevant parameters on the initial breakdown phase for CHI have not yet been considered. Here, we evaluate CHI breakdown in reactor-like configurations using an extension of the Townsend avalanche theory. We find that a CHI electrode concept in which the outer vessel wall is biased to achieve breakdown, while previously successful on NSTX and HIT-II, may exhibit a severe weakness when scaled up to a reactor. On the other hand, concepts which employ localized biasing electrodes such as those used in QUEST would avoid this issue. Assuming that breakdown can be successfully attained, we then apply scaling relationships to predict plasma parameters attainable in the transient CHI discharge. Assuming the use of 1 Wb of injector flux, we find that plasma currents of 1 MA should be achievable. Furthermore, these plasmas are expected to Ohmically self-heat with more than 1 MW of power as they decay, facilitating efficient hand-off to steady-state heating sources. These optimistic scalings are supported by TSC simulations.
Type:
Dataset
Issue Date:
February 2019

30. Attention and awareness in the dorsal attention network

Author(s):
Wilterson, Andrew; Nastase, Samuel; Bio, Branden; Guterstam, Arvid; Graziano, Michael
Abstract:
The attention schema theory (AST) posits a specific relationship between subjective awareness and attention, in which awareness is the control model that the brain uses to aid in the endogenous control of attention. We proposed that the right temporoparietal junction (TPJ) is involved in that interaction between awareness and attention. In previous experiments, we developed a behavioral paradigm in human subjects to manipulate awareness and attention. The paradigm involved a visual cue that could be used to guide a shift of attention to a target stimulus. In task 1, subjects were aware of the visual cue, and their endogenous control mechanism was able to use the cue to help control attention. In task 2, subjects were unaware of the visual cue, and their endogenous control mechanism was no longer able to use it to control attention, even though the cue still had a measurable effect on other aspects of behavior. Here we tested the two tasks while scanning brain activity in human volunteers. We predicted that the right TPJ would be active in relation to the cue in task 1, but not in task 2. This prediction was confirmed. The right TPJ was active in relation to the cue in task 1; it was not measurably active in task 2; the difference was significant. In our interpretation, the right TPJ is involved in a complex interaction in which awareness aids in the control of attention.
Type:
Dataset
Issue Date:
2020

31. Attentional Modulation of Brain Responses to Primary Appetitive and Aversive Stimuli

Author(s):
Cara L. Buck; Jonathan D. Cohen; Field, Brent; Daniel Kahneman; Samuel M. McClure; Leigh E. Nystrom
Abstract:
Studies of subjective well-being have conventionally relied upon self-report, which directs subjects’ attention to their emotional experiences. This method presumes that attention itself does not influence emotional processes, which could bias sampling. We tested whether attention influences experienced utility (the moment-by-moment experience of pleasure) by using functional magnetic resonance imaging (fMRI) to measure the activity of brain systems thought to represent hedonic value while manipulating attentional load. Subjects received appetitive or aversive solutions orally while alternatively executing a low or high attentional load task. Brain regions associated with hedonic processing, including the ventral striatum, showed a response to both juice and quinine. This response decreased during the high-load task relative to the low-load task. Thus, attentional allocation may influence experienced utility by modulating (either directly or indirectly) the activity of brain mechanisms thought to represent hedonic value.
Type:
Dataset, Software, and text
Issue Date:
11 February 2015

33. Basin-Scale Leakage Risks from Geologic Carbon Sequestration: Impact on Carbon Capture and Storage Energy Market Competitiveness

Abstract:
This three-year project, performed by Princeton University in partnership with the University of Minnesota and Brookhaven National Laboratory, examined geologic carbon sequestration in regard to CO2 leakage and potential subsurface liabilities. The research resulted in basin-scale analyses of CO2 and brine leakage in light of uncertainties in the characteristics of leakage processes, and generated frameworks to monetize the risks of leakage interference with competing subsurface resources. The geographic focus was the Michigan sedimentary basin, for which a 3D topographical model was constructed to represent the hydrostratigraphy. Specifically for Ottawa County, a statistical analysis of the hydraulic properties of underlying sedimentary formations was conducted. For plausible scenarios of injection into the Mt. Simon sandstone, leakage rates were estimated and fluxes into shallow drinking-water aquifers were found to be less than natural analogs of CO2 fluxes. We developed the Leakage Impact Valuation (LIV) model in which we identified stakeholders and estimated costs associated with leakage events. It was found that costs could be incurred even in the absence of legal action or other subsurface interference because there are substantial costs of finding and fixing the leak and from injection interruption. We developed a model framework called RISCS, which can be used to predict monetized risk of interference with subsurface resources by combining basin-scale leakage predictions with the LIV method. The project has also developed a cost calculator called the Economic and Policy Drivers Module (EPDM), which comprehensively calculates the costs of carbon sequestration and leakage, and can be used to examine major drivers for subsurface leakage liabilities in relation to specific injection scenarios and leakage events. Finally, we examined the competitiveness of CCS in the energy market. This analysis, though qualitative, shows that financial incentives, such as a carbon tax, are needed for coal combustion with CCS to gain market share. In another part of the project we studied the role of geochemical reactions in affecting the probability of CO2 leakage. A basin-scale simulation tool was modified to account for changes in leakage rates due to permeability alterations, based on simplified mathematical rules for the important geochemical reactions between acidified brines and caprock minerals. In studies of reactive flows in fractured caprocks, we examined the potential for permeability increases, and the extent to which existing reactive transport models would or would not be able to predict it. Using caprock specimens from the Eau Claire and Amherstburg, we found that substantial increases in permeability are possible for caprocks that have significant carbonate content, but minimal alteration is expected otherwise. We also found that while the permeability increase may be substantial, it is much less than what would be predicted from hydrodynamic models based on mechanical aperture alone because the roughness that is generated tends to inhibit flow.
Type:
text
Issue Date:
5 June 2014

34. Blob Structure and Motion in the Edge and SOL of NSTX

Author(s):
Zweben, S.J.; J.R. Myra; W.M. Davis; D.A. D'Ippolito; T.K. Gray; S.M. Kaye; B.P. LeBlanc; R.J. Maqueda; D.A. Russell; D.P. Stotler
Abstract:
Blob analysis dataset
Type:
Dataset
Issue Date:
January 2016

35. Blob wakes in NSTX

Author(s):
Zweben SJ, Myra JR, Diallo A, Russell DA, Scotti F, Stotler DP
Abstract:
Transient small-scale structures were identified in the wake of blobs movingpoloidally through the SOL of high-powered H-mode plasmas in NSTX, using the gaspuff imaging (GPI) diagnostic. These blob wakes had a poloidal wavelength in therange 3.5 cm, which is significantly smaller than the average blob scale of~12 cm, and the wakes had a poloidal velocity of 1.5 km/sec in theelectron diamagnetic direction, which is opposite to the blob poloidal velocity inthese shots. These wakes were radially localized 0-4 cm outside the separatrix andoccurred within ~50 microsec after the passage of a blob through the GPI field of view.The clearest wakes were seen when the GPI viewing angle was well aligned with thelocal B field line, as expected for such small-scale structures given the diagnosticgeometry. A plausible theoretical interpretation of the wakes is discussed: theobserved wakes share some features of drift waves and/or drift-Alfven waves whichcould be excited
Type:
Dataset
Issue Date:
July 2019

36. Blob-hole correlation model for edge turbulence and comparisons with NSTX GPI data

Author(s):
Myra, J.R.; Zweben, S.J.; Russell, D.A.
Abstract:
Gas puff imaging (GPI) observations made in NSTX [Zweben S J, et al., 2017 Phys. Plasmas 24 102509] have revealed two-point spatial correlations of edge and scrape-off layer turbulence in the plane perpendicular to the magnetic field. A common feature is the occurrence of dipole-like patterns with significant regions of negative correlation. In this paper, we explore the possibility that these dipole patterns may be due to blob-hole pairs. Statistical methods are applied to determine the two-point spatial correlation that results from a model of blob-hole pair formation. It is shown that the model produces dipole correlation patterns that are qualitatively similar to the GPI data in several respects. Effects of the reference location (confined surfaces or scrape-off layer), a superimposed random background, hole velocity and lifetime, and background sheared flows are explored and discussed with respect to experimental observations. Additional analysis of the experimental GPI dataset is performed to further test this blob-hole correlation model. A time delay two-point spatial correlation study did not reveal inward propagation of the negative correlation structures that were postulated to correspond to holes in the data nor did it suggest that the negative correlation structures are due to neutral shadowing. However, tracing of the highest and lowest values (extrema) of the normalized GPI fluctuations shows strong evidence for mean inward propagation of minima and outward propagation of maxima, in qualitative agreement with theoretical expectations. Other properties of the experimentally observed extrema are discussed.
Type:
Dataset
Issue Date:
July 2018

37. Bubble pinch-off in turbulence

Abstract:
Experimental data, simulation code, and Python scripts to reproduce the data presented in "Bubble pinch-off in turbulence"
Type:
Dataset and Software

38. CKavity Library: Next-Generation Sequencing

Abstract:
Protein sequence space is vast; nature uses only an infinitesimal fraction of possible sequences to sustain life. Are there solutions to biological problems other than those provided by nature? Can we create artificial proteins that sustain life? To investigate this question, the Hecht lab has created combinatorial collections, or libraries, of novel sequences with no homology to those found in living organisms. These libraries were subjected to screens and selections, leading to the identification of sequences with roles in catalysis, modulating gene regulation, and metal homeostasis. However, the resulting functional proteins formed dynamic rather than well-ordered structures. This impeded structural characterization and made it difficult to ascertain a mechanism of action. To address this, Christina Karas's thesis work focuses on developing a new model of libraries based on the de novo protein S-824, a four-helix bundle with a very stable three-dimensional structure. The first part of this research focused on mutagenesis of S-824 and characterization of the resulting proteins, revealing that this scaffold tolerates amino acid substitutions, including buried polar residues and the removal of hydrophobic side chains to create a putative cavity. Distinct from previous libraries, Karas targeted variability to a specific region of the protein, seeking to create a cavity and potential active site. The second part of this work details the design and creation of a library encoding 1.7 x 10^6 unique proteins, assembled from degenerate oligonucleotides. The third and fourth parts of this work cover the screening effort for a range of activities, both in vitro and in vivo. I found that this collection binds heme readily, leading to abundant peroxidase activity. Hits for lipase and phosphatase activity were also detected. This work details the development of a new strategy for creating de novo sequences geared toward function rather than structure.
Type:
Dataset

39. CLEVR-Matrices

Author(s):
Mondal, Shanka Subhra; Webb, Taylor; Cohen, Jonathan
Abstract:
A dataset of Raven’s Progressive Matrices (RPM)-like problems using realistically rendered 3D shapes, based on source code from CLEVR (a popular visual-question-answering dataset) (Johnson, J., Hariharan, B., Van Der Maaten, L., Fei-Fei, L., Lawrence Zitnick, C., & Girshick, R. (2017). Clevr: A diagnostic dataset for compositional language and elementary visual reasoning. In Proceedings of the IEEE conference on computer vision and pattern recognition (pp. 2901-2910)).
Type:
Dataset
Issue Date:
1 March 2023

40. CPCB continuous monitoring air quality data for India, from 2015-2019

Author(s):
Sharma, Disha; Mauzerall, Denise
Abstract:
The dataset is a compilation of real time ground observations of criteria pollutants monitored at the Central Pollution Control Board (CPCB) continuous stations in India, from 2015-2019. Pollutants included are PM2.5, PM10, SO2, NO2 and O3 and are archived at every hour for all stations across India.
Type:
Dataset
Issue Date:
15 October 2021

42. Calibrated VIC Model Parameters over CONUS

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

43. Calibrationless rotating Lorentz-force flowmeters for low flow rate applications

Author(s):
Hvasta, M. G.; Dudt, D.; Fisher, A. E.; Kolemen, E.
Abstract:
A 'weighted magnetic bearing' has been developed to improve the performance of rotating Lorentz-force flowmeters (RLFFs). Experiments have shown that the new bearing reduces frictional losses within a double-sided, disc-style RLFF to negligible levels. Operating such an RLFF under 'frictionless' conditions provides two major benefits. First, the steady-state velocity of the RLFF magnets matches the average velocity of the flowing liquid at low flow rates. This enables an RLFF to make accurate volumetric flow measurements without any calibration or prior knowledge of the fluid properties. Second, due to minimized frictional losses, an RLFF is able to measure low flow rates that cannot be detected when conventional, high-friction bearings are used. This paper provides a brief background on RLFFs, gives a detailed description of weighted magnetic bearings, and compares experimental RLFF data to measurements taken with a commercially available flowmeter.
Type:
Dataset
Issue Date:
29 May 2018

44. Chinese overseas development financing of electric power generation: A comparative analysis

Author(s):
Chen, Xu; Gallagher, Kevin P.; Mauzerall, Denise L.
Abstract:
Global power generation must rapidly decarbonize by mid-century to meet the goal of stabilizing global warming below 2 degree Celsius. To meet this objective, multilateral development banks (MDBs) have gradually reduced fossil fuel and increased renewable energy financing. Meanwhile, globally active national development finance institutions (DFIs) from Japan and South Korea have continued to finance overseas coal plants. Less is known about the increasingly active Chinese DFIs. Here we construct a new dataset of China’s policy banks’ overseas power generation financing and compare their technology choices and impact on generation capacity with MDBs and Japanese and South Korean DFIs. We find Chinese DFI power financing since 2000 has dramatically increased, surpassing other East Asian national DFIs and the major MDBs’ collective public sector power financing in 2013. As most Chinese DFI financing is currently in coal, decarbonization of their power investments will be critical in reducing future carbon emissions from recipient countries.
Type:
Dataset
Issue Date:
23 October 2020

45. Chiral Tetramer Molecular Model

Author(s):
Petsev, Nikolai D.; Stillinger, Frank H.; Debenedetti, Pablo G.
Abstract:
Source code for our energy-conserving reformulation of the 4-site molecular model for chiral phenomena originally introduced by Latinwo et al. [F. Latinwo, F. H. Stillinger, and P. G. Debenedetti, Molecular Model for Chirality Phenomena, J. Chem. Phys. 145, 154503 (2016)]. The reformulation includes an additional 8-body force that arises from an explicit configuration-dependent term in the potential energy function, resulting in a coarse-grained energy-conserving force field for molecular dynamics simulations of chirality phenomena. In this model, the coarse-grained interaction energy between two tetramers depends on their respective chiralities, and is controlled by a parameter λ, where favors local configurations involving tetramers of opposite chirality, and gives energetic preference to configurations involving tetramers of the same chirality. The source code is for use with the LAMMPS simulation package.
Type:
Software
Issue Date:
3 August 2021

46. Climate Impacts from Large Volcanic Eruptions in a High-resolution Climate Model: the Importance of Forcing Structure

Author(s):
Yang, Wenchang; Vecchi, Gabriel; Fueglistaler, Stephan; Horowitz, Larry; Luet, David; Muñoz, Ángel; Paynter, David; Underwood, Seth
Abstract:
Explosive volcanic eruptions have large climate impacts, and can serve as observable tests of the climatic response to radiative forcing. Using a high resolution climate model, we contrast the climate responses to Pinatubo, with symmetric forcing, and those to Santa Maria and Agung, which had meridionally asymmetric forcing. Although Pinatubo had larger global-mean forcing, asymmetric forcing strongly shifts the latitude of tropical rainfall features, leading to larger local precipitation/TC changes. For example, North Atlantic TC activity over is enhanced/reduced by SH-forcing (Agung)/NH-forcing (Santa Maria), but changes little in response to the Pinatubo forcing. Moreover, the transient climate sensitivity estimated from the response to Santa Maria is 20% larger than that from Pinatubo or Agung. This spread in climatic impacts of volcanoes needs to be considered when evaluating the role of volcanoes in global and regional climate, and serves to contextualize the well-observed response to Pinatubo.
Type:
Dataset
Issue Date:
2019

47. Code and data from "Comparative genomic analysis reveals varying levels of mammalian adaptation to coronavirus infections"

Author(s):
King, Sean
Abstract:
Severe acute respiratory coronavirus 2 (SARS-CoV-2), the causative agent of COVID-19, is of zoonotic origin. Evolutionary analyses assessing whether coronaviruses similar to SARS-CoV-2 infected ancestral species of modern-day animal hosts could be useful in identifying additional reservoirs of potentially dangerous coronaviruses. We reasoned that if a clade of species has been repeatedly exposed to a virus, then their proteins relevant for viral entry may exhibit adaptations that affect host susceptibility or response. We perform comparative analyses across the mammalian phylogeny of angiotensin-converting enzyme 2 (ACE2), the cellular receptor for SARS-CoV-2, in order to uncover evidence for selection acting at its binding interface with the SARS-CoV-2 spike protein. We uncover that in rodents there is evidence for adaptive amino acid substitutions at positions comprising the ACE2-spike interaction interface, whereas the variation within ACE2 proteins in primates and some other mammalian clades is not consistent with evolutionary adaptations. We also analyze aminopeptidase N (APN), the receptor for the human coronavirus 229E, a virus that causes the common cold, and find evidence for adaptation in primates. Altogether, our results suggest that the rodent and primate lineages may have had ancient exposures to viruses similar to SARS-CoV-2 and HCoV-229E, respectively. Included in this repository are the instructions and corresponding code required to build the dataset and run the analysis in the manuscript.
Type:
Dataset
Issue Date:
28 September 2021

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

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

49. Collisional dependence of Alfven mode saturation in tokamaks

Author(s):
Zhou, M.; White, R.
Abstract:
Saturation of \alfven modes driven unstable by a distribution of high energy particles as a function of collisionality is investigated with a guiding center code, using numerical eigenfunctions produced by linear theory and numerical high energy particle distributions. The most important resonance is found and it is shown that when the resonance domain is bounded, not allowing particles to collisionlessly escape, the saturation amplitude is given by the balance of the resonance mixing time with the time for nearby particles to collisionally diffuse across the resonance width. Saturation amplitudes are in agreement with theoretical predictions as long as the mode amplitude is not so large that it produces stochastic loss from the resonance domain.
Type:
Dataset
Issue Date:
December 2016

50. Comment on ‘Numerical modeling of tokamak breakdown phase driven by pure Ohmic heating under ideal conditions’

Author(s):
Yoo, Min-Gu; Na, Yong-Su
Abstract:
In this comment, we point out possible critical numerical flaws of recent particle simulation studies (Jiang et al 2016 Nucl. Fusion 56 126017, Peng et al 2018 Nucl. Fusion 58 026007) on the electrical gas breakdown in a simple one-dimensional periodic slab geometry. We show that their observations on the effects of the ambipolar electric fields during the breakdown, such as the sudden reversal of the ion flow direction, could not be real physical phenomena but resulting from numerical artifacts violating the momentum conservation law. We show that an incomplete implementation of the direct-implicit scheme can cause the artificial electric fields and plasma transports resulting in fallacies in simulation results. We also discuss that their simple plasma model without considering poloidal magnetic fields seriously mislead the physical mechanism of the electrical gas breakdown because it cannot reflect important dominant plasma dynamics in the poloidal plane (Yoo et al 2018 Nat. Commun. 9 3523).
Type:
Dataset
Issue Date:
June 2019