Runaway Regime Research Articles

Generation of Runaway Electrons near Micro-Inhomogeneities on the Cathode Surface in Subnanosecond Self-Sustained Discharges in a Wide Range of High Pressures

The results of numerical 3D modeling of the development of an electron avalanche initiated by a field emission electron in a small-sized region of an amplified electric field near the microinhomogeneities at the cathode have been presented. The simulation has been carried out for the discharge gaps with an initially homogeneous distribution of the electric field with a reduced intensity significantly lower than that required by the electron runaway criterion. The possibility of the transition of the field emission electrons initiating avalanches and the electrons in these avalanches into runaway regime has been investigated. The microinhomogeneities in the form of a cone, metal droplets, and boundaries between pores or microcraters have been considered. The calculations were carried out for nitrogen in the pressure range from atmospheric to 40 atm. It has been shown that the initial energy obtained near the microinhomogeneity can significantly facilitate the transition of the electron into the runaway mode. And the electron will continue to run away in a discharge gap electric field weak according to the runaway criterion. It has been shown that this effect is especially noticeable at gas pressures above 10 atm. A comparative analysis of the simulation results with the experimental data obtained by us on the switching characteristics of a discharge gap filled with nitrogen when exposed to voltage pulses with subnanosecond fronts of different steepness has been carried out. This made it possible to divide the ranges of experimental conditions into those when only the amplification of the electric field near the microinhomogeneities is sufficient for the runaway of electrons and when the electric field of an avalanche of critical or close to critical size is additionally necessary for the runaway.

Dynamics and products of potato crop residue conversion under a pyrolytic runaway regime - Influences of feedstock variability

Dynamics and products of potato crop residue conversion under a pyrolytic runaway regime - Influences of feedstock variability

The maximum accretion rate of a protoplanet: how fast can runaway be?

ABSTRACT The hunt is on for dozens of protoplanets hypothesized to reside in protoplanetary discs with imaged gaps. How bright these planets are, and what they will grow to become, depend on their accretion rates, which may be in the runaway regime. Using 3D global simulations, we calculate maximum gas accretion rates for planet masses Mp from 1$\, \mathrm{ M}_{{\oplus }}$ to $10\, \mathrm{ M}_{\rm J}$. When the planet is small enough that its sphere of influence is fully embedded in the disc, with a Bondi radius rBondi smaller than the disc’s scale height Hp – such planets have thermal mass parameters qth ≡ (Mp/M⋆)/(Hp/Rp)3 ≲ 0.3, for host stellar mass M⋆ and orbital radius Rp – the maximum accretion rate follows a Bondi scaling, with $\max \dot{M}_{\rm p} \propto \rho _{\rm g}M_{\rm p}^2 / (H_{\rm p}/R_{\rm p})^3$ for ambient disc density ρg. For more massive planets with 0.3 ≲ qth ≲ 10, the Hill sphere replaces the Bondi sphere as the gravitational sphere of influence, and $\max \dot{M}_{\rm p} \propto \rho _{\rm g}M_{\rm p}^1$, with no dependence on Hp/Rp. In the strongly superthermal limit when qth ≳ 10, the Hill sphere pops well out of the disc, and $\max \dot{M}_{\rm p} \propto \rho _{\rm g}M_{\rm p}^{2/3} (H_{\rm p}/R_{\rm p})^1$. Applied to the two confirmed protoplanets PDS 70b and c, our numerically calibrated maximum accretion rates imply that their Jupiter-like masses may increase by up to a factor of ∼2 before their parent disc dissipates.

Computational modelling and statistical evaluation of thermal runaway safety regime response on lithium-ion battery with different cathodic chemistry and varying ambient condition

As one of the most promising energy storage mediums, Lithium-ion batteries (LIBs) have attracted extensive research interest. A major challenge associated with LIB application is thermal runaway, which can be triggered under abused conditions and impose direct threats to life and properties. Thermal runaway is a major source of fire and explosions in the battery energy storage industry. In this paper, the thermal runaway numerical model is used to explain the evolution of thermal runaway trigger point over different cathodic chemistry under the influence of varying ambient temperature. The test is performed using nail penetration assisted thermal abuse method. The model is based on the conservation of energy Equation, with the empirical corelation of the thermal abuse equations which is used to simulate the source term. Additionally, experimental literature is used to verify the simulation results extensively. The simulation results are utilized to monitor internal short circuit circumstances that lead to the thermal runaway trigger point, and a comparison scaled map is created between various chemistry to show which chemistry triggers prematurely between LiNiMnCoO2 (NMC), LiFePO4 (LFP), and LiCoO2 (LCO). The existing critical uncertainty, such as nail diameter, heating profile, variation in operating ambient temperature, battery jelly roll thickness and effect of convective parameters is investigated in detailed and decision tree is prepared using ANOVA statistical method.

Chapman–Jouguet deflagration criteria and compressibility dynamics of turbulent fast flames for turbulence-induced deflagration-to-detonation transition

This work characterizes the compressibility dynamics in turbulent fast flames for a range of turbulent flame speeds. These turbulent fast flames experience increased effects of compressibility through the formation of strong shocks and may develop a runaway acceleration combined with a pressure buildup that leads to turbulence induced deflagration-to-detonation transition (tDDT). Simultaneous high-speed particle image velocimetry, OH* chemiluminescence, schlieren, and pressure measurements are used to examine the reacting flow field and flame dynamics. We examine flames with turbulent flame speeds ranging from 100 to 600 m/s. At lower turbulent flame speeds, the flame is not able to produce favorable background conditions for deflagration-to-detonation transition (DDT) onset, and thus flame compressibility and turbulence amplification are less dominant, resulting in a weaker acoustic coupling between the flame and compressed region. As the turbulent burning velocities exceed the Chapman–Jouguet deflagration speed, favorable background conditions are produced, as we observe flame-generated shocks and flame-generated turbulence with higher turbulent velocities and larger turbulent scales. At this regime, the flame is categorized to be at the runaway transition regime that leads to tDDT.

Stability of the Néel quantum critical point in the presence of Dirac fermions

We investigate the stability of the N\'eel quantum critical point of two-dimensional quantum antiferromagnets, described by a nonlinear $\ensuremath{\sigma}$ model, in the presence of a Kondo coupling to ${N}_{f}$ flavors of two-component Dirac fermion fields. The long-wavelength order parameter fluctuations are subject to Landau damping by electronic particle-hole fluctuations. Using the momentum-shell renormalization group (RG), we demonstrate that the Landau damping is weakly irrelevant at the N\'eel quantum critical point, despite the fact that the corresponding self-energy correction dominates over the quadratic gradient terms in the IR limit. In the ordered phase, the Landau damping increases under the RG, indicative of damped spin-wave excitations. Although the Kondo coupling is weakly relevant, sufficiently strong Landau damping renders the N\'eel quantum critical point quasistable for ${N}_{f}\ensuremath{\ge}4$ and thermodynamically stable for ${N}_{f}<4$. In the latter case, we identify a multicritical point which describes the transition between the N\'eel critical and Kondo runaway regimes. The symmetry breaking at this fixed point results in the opening of a gap in the Dirac fermion spectrum. Approaching the multicritical point from the disordered phase, the fermionic quasiparticle residue vanishes, giving rise to non-Fermi-liquid behavior.

Критерии убегания электронов в газовом диоде с игольчатым катодом

The conditions of electron runaway in a gas diode with a cathode in the form of a needle are studied theoretically. It is shown that the runaway conditions are qualitatively different for needles with relatively large and small tip radii, i.e., in fact, for different degrees of electric field inhomogeneity. In a weakly inhomogeneous field, the transition of electrons to the runaway regime is determined by the local distribution of the field near the place of their start – the tip of the needle. In a strongly inhomogeneous field, the runaway condition has a nonlocal character: it is determined by the behavior of electrons in the near-anode region. This difference leads to a nonmonotonic dependence of the threshold runaway voltage on the tip radius.

Global model on oscillation discharge characteristics during deep oscillation magnetron sputtering of Cr target

A time-dependent global model is developed for deep oscillation magnetron sputtering (DOMS) with a Cr target to characterize the plasma using the working Ar gas at a pressure of 0.8 Pa. An input of the global model is the target voltage and current waveforms at charging voltages from 260 to 400 V. The plasma densities vary synchronously with the oscillations of target power. During the on-time of a micropulse, the gas atoms are ionized first to ignite the plasma and subsequently the metal ions follow to sustain the discharge. During the off-time, the gas ions dominate the afterglow. The DOMS possesses a characteristic of alternating gas/metal discharge in the time domain. The peak values of plasma densities increase linearly with charging voltages, which is mainly attributed to the higher Cr+ ion density. The discharge transits from gas dominated to metal dominated at high charging voltage. The working gas rarefaction in DOMS discharge is enhanced at higher charging voltages. At charging voltages higher than 360 V, the metal self-sputtering comes into the runaway regime temporarily as indicated by the self-sputtering parameters exceeding unity, generating the dense and metal-rich plasma. The metal self-sputtering is promoted by gas rarefaction. The electrons tend to ionize the metal atoms after the gas atoms are depleted. The decreased grain size and elevated nano-hardness of the Cr thin films are explained by structure transition from zone I to zone T due to the efficient metal ion bombardment to the growing thin films as the charging voltage increases.

Spin impurities, Wilson lines and semiclassics

We consider line defects with large quantum numbers in conformal field theories. First, we consider spin impurities, both for a free scalar triplet and in the Wilson-Fisher O(3) model. For the free scalar triplet, we find a rich phase diagram that includes a perturbative fixed point, a new nonperturbative fixed point, and runaway regimes. To obtain these results, we develop a new semiclassical approach. For the Wilson-Fisher model, we propose an alternative description, which becomes weakly coupled in the large spin limit. This allows us to chart the phase diagram and obtain numerous rigorous predictions for large spin impurities in 2 + 1 dimensional magnets. Finally, we also study 1/2-BPS Wilson lines in large representations of the gauge group in rank-1 mathcal{N} = 2 superconformal field theories. We contrast the results with the qualitative behavior of large spin impurities in magnets.

Spontaneous runaway of fast turbulent flames for turbulence-induced deflagration-to-detonation transition

One of the fundamental mechanisms for detonation initiation is deflagration-to-detonation transition (DDT). This research experimentally explores the runaway condition for highly turbulent fast flames before DDT, which are characterized by extremely high turbulent flame speeds. Such fast turbulent flames experience increased effects of compressibility and may develop a runaway acceleration combined with a pressure buildup that leads to a turbulence-induced DDT (tDDT) mechanism that has been recently reported. The flame dynamics and the associated reacting flow field are characterized using simultaneous high-speed particle image velocimetry, OH* chemiluminescence, pressure measurements, and schlieren imaging. We study the flow-field conditions for runaway acceleration of fast turbulent flames and effects of compressibility on the evolution of these flames. The locally measured turbulent flame speed is found to be greater than that of a Chapman–Jouguet deflagration speed, which places the flame in the runaway transition regime that would eventually lead to a detonation.

Role of Suprathermal Runaway Electrons Returning to the Acceleration Region in Solar Flares

Abstract During solar flares, a large flux of energetic electrons propagate from the tops of reconnecting magnetic flux tubes toward the lower atmosphere. Over the course of the electrons’ transport, a co-spatial counterstreaming return current is induced, thereby balancing the current density. In response to the return-current electric field, a fraction of the ambient electrons will be accelerated into the runaway regime. However, models describing the accelerated electron beam/return-current system have not previously taken these suprathermal runaway electrons into account self-consistently. We develop a model in which an accelerated electron beam drives a steady-state, sub-Dreicer co-spatial return-current electric field, which locally balances the direct beam current and freely accelerates a fraction of background (return-current) electrons. The model is self-consistent, i.e., the electric field induced by the coevolution of the direct beam and the runaway current is considered. We find that (1) the return-current electric field can return a significant number of suprathermal electrons to the acceleration region, where they can be further accelerated to higher energies, runaway electrons can be a few tens of percent of the return-current flux returning to the nonthermal beam’s acceleration region, (2) the energy gain of the suprathermal electrons can be up to 10–35 keV, (3) the heating rate in the corona can be reduced by a factor of 3 for medium range injected fluxes in comparison to models which neglect the runaway component. The results depend on the injected beam flux density, the temperature, and density of the background plasma.

Steady state by recycling prevents premature collapse of protoplanetary atmospheres

Context. In recent years, space missions such as Kepler and TESS have discovered many close-in planets with significant atmospheres consisting of hydrogen and helium: mini-Neptunes. This indicates that these planets formed early in gas-rich disks while avoiding the runaway gas accretion that would otherwise have turned them into hot-Jupiters. A solution is to invoke a long Kelvin-Helmholtz contraction (or cooling) timescale, but it has also been suggested that thermodynamical cooling can be prevented by hydrodynamical planet atmosphere-disk recycling. Aims. We investigate the efficacy of the recycling hypothesis in preventing the collapse of the atmosphere, check for the existence of a steady state configuration, and determine the final atmospheric mass to core mass ratio. Methods. We use three-dimensional radiation-hydrodynamic simulations to model the formation of planetary proto-atmospheres. Equations are solved in a local frame centered on the planet. Results. Ignoring small oscillations that average to zero over time, the simulations converge to a steady state where the velocity field of the gas becomes constant in time. In a steady state, the energy loss by radiative cooling is fully compensated by the recycling of the low entropy gas in the planetary atmosphere with high entropy gas from the circumstellar disk. Conclusions. For close-in planets, recycling naturally halts the cooling of planetary proto-atmospheres, preventing them from contracting toward the runaway regime and collapsing into gas giants.

The Effect of Land Albedo on the Climate of Land-dominated Planets in the TRAPPIST-1 System

Abstract Variations in the reflective properties of the bulk material that comprises the surface of land-dominated planets will affect the planetary energy balance by interacting differently with incident radiation from the host star. Furthermore, low-mass cool stars, such as nearby M8V dwarf TRAPPIST-1, emit a significant fraction of their flux in longer wavelengths relative to the Sun in regions where terrestrial materials may exhibit additional variability in albedo. Using the Community Earth System Model, we investigate the effect of the composition of the land surface and its albedo on planetary climate in the context of spatially homogeneous, entirely land-covered planets with dry atmospheres at the orbital separation of TRAPPIST-1d, TRAPPIST-1e, and TRAPPIST-1f. We use empirically derived spectra of four terrestrial compositional endmembers (granite, calcite, aridisol, and dune sand) and a composite spectrum of TRAPPIST-1 for these simulations and compare these model outputs to an aquaplanet and several Sol-spectrum control cases. We report a difference of approximately 50 K in global mean surface temperature, variations in atmospheric rotational features, and a reduction in cross-equatorial heat transport between scenarios in which materials with higher albedo in the infrared (calcite and dune sand) were used and those with more absorptive crustal material, such as granite or dry soils. An aquaplanet TRAPPIST-1d scenario results in an unstable runaway greenhouse regime. Therefore, we demonstrate that determining the composition and albedo of continental landmasses is crucial for making accurate determinations of the climate of terrestrial exoplanets.

Modeling the spatial dynamics of marsh ponds in New England salt marshes

Ponds are common features on salt marshes, yet it is unclear how they affect large-scale marsh evolution. We developed a spatially explicit model that combines cellular automata for pond formation, expansion, and drainage, and partial differential equations for elevation dynamics. We use the mesotidal Barnstable marsh (MA, USA) as a case study, for which we measured pond expansion rate by remote sensing analysis over a 41-year time span. We estimated pond formation rate by comparing observed and modeled pond size distribution, and predicted pond deepening by comparing modeled and measured pond depth. The Barnstable marsh is currently in the pond recovery regime, i.e.,every pond revegetates and recovers the necessary elevation to support plant growth after re-connecting to the channel network. This pond dynamic creates an equivalent (i.e.,spatially and temporally averaged over the whole marsh) 0.5–2 mm/yr elevation loss that needs to be supplemented by excess vertical accretion. We explore how the pond regime would change with decreased sediment supply and increased relative sea-level rise (RSLR) rate, focusing on the case in which the vegetated marsh keeps pace with RSLR. When the RSLR rate remains below the minimum unvegetated deposition rate, the pond dynamics is nearly unaltered and ponds always occupy ~10% of the marsh area. However, when RSLR rate exceeds this threshold, the ponds in the marsh interior – which receive the least amount of suspended sediment – do not recover after drainage. These ponds transition to mudflats and permanently occupy up to 30% of the marsh area depending on RSLR rate. For marshes with a small tidal range, such as the microtidal Sage Lot Pond marsh on the opposite side of the peninsula from Barnstable marsh, high RSLR rates could bring every portion of the marsh into the pond runaway regime, with the whole marsh eventually converting into mudflats. In this regime, the existing marsh would disappear within centuries to millennia depending on the RSLR rate. Because of the spatial and temporal components of marsh evolution, a single RSLR threshold value applied across the entire marsh landscape provides a limited description of the marsh vulnerability to RSLR.

Self-gravitating planetary envelopes and the core-nucleated instability

Abstract Planet formation scenarios can be constrained by the ratio of the gaseous envelope mass relative to the solid core mass in the observed exoplanet populations. One-dimensional calculations find a critical (maximal) core mass for quasi-static envelopes to exist, suggesting that envelopes around more massive cores should collapse due to a ‘core-nucleated’ instability. We study self-gravitating planetary envelopes via hydrodynamic simulations, progressively increasing the dimensionality of the problem. We characterize the core-nucleated instability and its non-linear evolution into runaway gas accretion in one-dimensional spherical envelopes. We show that rotationally supported envelopes can enter a runaway accretion regime via polar shocks in a two-dimensional axisymmetric model. This picture remains valid for high-mass cores in three dimensions, where the gas gravity mainly adds up to the core gravity and enhances the mass accretion rate of the planet in time. We relate the core-nucleated instability to the absence of equilibrium connecting the planet to its parent disc and discuss its relevance for massive planet formation.

Venus as a Laboratory for Exoplanetary Science

AbstractThe current goals of the astrobiology community are focused on developing a framework for the detection of biosignatures, or evidence thereof, on objects inside and outside of our solar system. A fundamental aspect of understanding the limits of habitable environments (surface liquid water) and detectable signatures thereof is the study of where the boundaries of such environments can occur. Such studies provide the basis for understanding how a once inhabitable planet might come to be uninhabitable. The archetype of such a planet is arguably Earth's sibling planet, Venus. Given the need to define the conditions that can rule out bio‐related signatures of exoplanets, Venus provides a unique opportunity to explore the processes that led to a completely uninhabitable environment by our current definition of the term. Here we review the current state of knowledge regarding Venus, particularly in the context of remote‐sensing techniques that are being or will be employed in the search for and characterization of exoplanets. We discuss candidate Venus analogs identified by the Kepler and TESS exoplanet missions and provide an update to exoplanet demographics that can be placed in the potential runaway greenhouse regime where Venus analogs are thought to reside. We list several major outstanding questions regarding the Venus environment and the relevance of those questions to understanding the atmospheres and interior structure of exoplanets. Finally, we outline the path toward a deeper analysis of our sibling planet and the synergy to exoplanetary science.

Cold Electron Runaway Below the Friction Curve

AbstractCold electron runaway means that free electrons in gases are accelerated by electric fields from eV energies to energies above tens of keV where they can be accelerated further. To run away, the electrons need to overcome a barrier at intermediate energies where they can lose much energy in collisions. When they have reached the runaway regime, they can produce high‐energy radiation by bremsstrahlung that can be detected as (terrestrial) gamma ray flashes. When can thermal electrons from active discharges like streamers and leaders reach the runaway regime? The deterministic approach to this question is based on an energy‐dependent electron friction that has to be overcome by electric acceleration. Taking the stochastic nature of the electron molecule collisions into account, we find (1) that the classical friction curve in the energy regime up to 1 keV does not characterize the mean electron energy, but rather, it seems to approximate the upper limit of the electron energy distribution and (2) that electrons can “tunnel” through the barrier when the field is close to 3 MV/m, below the so‐called cold runaway threshold (or critical field) of approximately 26 MV/m in air at standard temperature and pressure. (3) This is only true in the bulk perspective where the electron liberation and attachment in a given electric field is taken into account in the continuously refreshing electron ensemble. In a flux simulation that follows individual electrons as long as they are free, electron attachment reduces electron runaway very strongly in air, differently to what we observe in nitrogen.

Numerical Analysis of Electron Runaway in the Presence of Enhanced Field in the Vicinity of a Microtip

Transition of field-emission electrons to the runaway regime in the region of enhanced electric field determined by the configuration of a microtip on a cathode is studied at several pressures of gas medium. The problem is solved using simulation of electron motion in the presence of nonuniform electric field with the aid of the Monte Carlo procedure in the 2D configuration. Nitrogen is used as a working gas. Passage through a relatively small region of the enhanced field in the vicinity of the microtip may substantially facilitate electron escape to the runaway regime, especially, at pressures of greater than 10 atm. In our opinion, the resulting runaway electrons may provide preionization of gas medium and formation of the initial stage of a 3D discharge.

Three regimes of high-voltage breakdown in helium

A quasi-analytical model was developed to map out the low-pressure (left-hand) branch of the Paschen curve at very high voltage when electrons are in the runaway regime and charge exchange/ionization avalanche sustained by ions and fast neutral atoms becomes important. The model was applied to helium gas between parallel-plate electrodes, at potentials ranging in magnitude between 10 and 1000 kV. The respective value of reduced electric field E/n varied in the range of 50−6000 kTd (1 kTd = 10−18 Vm2), with reduced density nd (where n is the gas density and d is the inter-electrode distance) on the order of 1020 m−2. Three regimes of the breakdown have been identified according to the relative share of impact ionization by electrons, by ions, and by fast neutrals. The analytically derived Paschen curve is compared to those obtained with a detailed particle-in-cell/Monte Carlo simulation, and also through experimental measurements (L Xu, A V Khrabrov, I D Kaganovich and T J Sommerer 2017 Phys. Plasmas 24 093511). The model provides accurate predictions for E/n up to ∼103 kTd, constrained by availability and quality of required input data.

Analysis of the Pyrolytic Runaway Dynamics during Agricultural Waste Conversion

Potato plant waste pyrolysis is investigated by means of a cylindrical packed-bed reactor exposed along the lateral surface to heat fluxes of 20.6–33.2 kW/m2 (corresponding to heating temperatures, Ts, around 530–770 K). After an abrupt transition from torrefaction, a pyrolytic runaway regime is established (570 K ≤ Ts ≤ 620 K), characterized by temperature overshoots up to 300 K and actual conversion times only 4–15% of the preheating times. A uniform radial temperature profile is first associated with a trigger point positioned at the bed center and the propagation of a progressively thick thermal wave toward the lateral surface and then with one or more trigger points and the almost simultaneous conversion of the entire bed. Subsequently the trigger point moves to the lateral surface and a thermal wave propagates toward the bed center. The latter heating dynamics are also qualitatively preserved for Ts > 620 K when, in the presence of increasing radial gradients, a self-controlled pyrolysis regime is established with the reaction heat mainly exploited to preheat the bed.