Search NASA⌕ Search

SEARCH · Search NASA

Results for “Kinetics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 811 records · Page 45
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Learning heterogeneous reaction kinetics from X-ray videos pixel by pixel

Reaction rates at spatially heterogeneous, unstable interfaces are notoriously difficult to quantify, yet are essential in engineering many chemical systems, such as batteries and electrocatalysts. Experimental characterizations of such materials by operando microscopy produce rich image datasets, but data-driven methods to learn physics from these images are still lacking because of the complex coupling of reaction kinetics, surface chemistry and phase separation. Here we show that heterogeneous reaction kinetics can be learned from in situ scanning transmission X-ray microscopy (STXM) images of carbon-coated lithium iron phosphate (LFP) nanoparticles. Combining a large dataset of STXM images with a thermodynamically consistent electrochemical phase-field model, partial differential equation (PDE)-constrained optimization and uncertainty quantification, we extract the free-energy landscape and reaction kinetics and verify their consistency with theoretical models. We also simultaneously learn the spatial heterogeneity of the reaction rate, which closely matches the carbon-coating thickness profiles obtained through Auger electron microscopy (AEM). Across 180,000 image pixels, the mean discrepancy with the learned model is remarkably small (<7%) and comparable with experimental noise. Our results open the possibility of learning nonequilibrium material properties beyond the reach of traditional experimental methods and offer a new non-destructive technique for characterizing and optimizing heterogeneous reactive surfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Crystallographic and kinetic analyses of human IPMK reveal disordered domains modulate ATP binding and kinase activity

Inositol polyphosphate multikinase (IPMK) is a member of the IPK-superfamily of kinases, catalyzing phosphorylation of several soluble inositols and the signaling phospholipid PI(4,5)P 2 (PIP 2 ). IPMK also has critical non-catalytic roles in p53, mTOR/Raptor, TRAF6 and AMPK signaling mediated partly by two disordered domains. Although IPMK non-catalytic functions are well established, it is less clear if the disordered domains are important for IPMK kinase activity or ATP binding. Here, kinetic and structural analyses of an engineered human IPMK lacking all disordered domains (ΔIPMK) are presented. Although the K M for PIP 2 is identical between ΔIPMK and wild type, ΔIPMK has a 1.8-fold increase in k cat for PIP 2 , indicating the native IPMK disordered domains decrease IPMK activity in vitro . The 2.5 Å crystal structure of ΔIPMK is reported, confirming the conserved ATP-grasp fold. A comparison with other IPK-superfamily structures revealed a putative “ATP-clamp” in the disordered N-terminus, we predicted would stabilize ATP binding. Consistent with this observation, removal of the ATP clamp sequence increases the K M for ATP 4.9-fold, indicating the N-terminus enhances ATP binding to IPMK. Together, these structural and kinetic studies suggest in addition to mediating protein-protein interactions, the disordered domains of IPMK impart modulatory capacity to IPMK kinase activity through multiple kinetic mechanisms.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Sulfate adenylyl transferase kinetics and mechanisms of metabolic inhibitors of microbial sulfate respiration

Abstract Sulfate analog oxyanions that function as selective metabolic inhibitors of dissimilatory sulfate reducing microorganisms (SRM) are widely used in ecological studies and industrial applications. As such, it is important to understand the mode of action and mechanisms of tolerance or adaptation to these compounds. Different oxyanions vary widely in their inhibitory potency and mechanism of inhibition, but current evidence suggests that the sulfate adenylyl transferase/ATP sulfurylase (Sat) enzyme is an important target. We heterologously expressed and purified the Sat from the model SRM, Desulfovibrio alaskensis G20. With this enzyme we determined the turnover kinetics (kcat, KM) for alternative substrates (molybdate, selenate, arsenate, monofluorophosphate, and chromate) and inhibition constants (KI) for competitive inhibitors (perchlorate, chlorate, and nitrate). These measurements enable the first quantitative comparisons of these compounds as substrates or inhibitors of a purified Sat from a respiratory sulfate reducer. We compare predicted half-maximal inhibitory concentrations (IC50) based on Sat kinetics with measured IC50 values against D. alaskensis G20 growth and discuss our results in light of known mechanisms of sensitivity or resistance to oxyanions. This analysis helps with the interpretation of recent adaptive laboratory evolution studies and illustrates the value of interpreting gene–microbe–environment interactions through the lens of enzyme kinetics.

Carlson, Hans K. (ORCID:0000000215835313)↗
DOE OSTI2018

Elementary kinetics of nitrogen electroreduction to ammonia on late transition metals

Developing small scale nitrogen to ammonia conversion devices could greatly help sustainable agriculture. The lack of effective cathode materials for reduction of nitrogen to ammonia puts a major limitation on developing small scale electrochemical ammonia production devices that can operate at low pressures and temperatures. In this study, the electrochemical nitrogen reduction reaction (NRR) mechanism is investigated over late transition metals. Activation barriers for possible rate limiting steps are calculated using density functional theory (DFT) methods and are converted to potential dependent electrochemical barriers. Associative and dissociative paths are evaluated and the associative path shows lower barriers on all metals considered at NRR potentials. Brønsted–Evans–Polanyi (BEP) relationships are evaluated for rate limiting steps and a “kinetic volcano” is demonstrated for catalyst optimization. Rhodium (Rh) and iron (Fe) appear to have the lowest kinetic barriers to convert N 2 to NH 3 . A large over-potential, however, is required to convert nitrogen to ammonia over all the surfaces considered. Finally, the hydrogen evolution reaction (HER) has a lower activation barrier compared to NRR, demonstrating the kinetic selectivity challenge. ; A

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Hydration structure and water exchange kinetics at xenotime–water interfaces: implications for rare earth minerals separation

Hydration of surface ions gives rise to structural heterogeneity and variable exchange kinetics of water at complex mineral–water interfaces. In this work, we employ ab initio molecular dynamics (AIMD) simulations and water adsorption calorimetry to examine the aqueous interfaces of xenotime, a phosphate mineral that contains predominantly Y 3+ and heavy rare earth elements. Consistent with natural crystal morphology, xenotime is predicted to have a tetragonal prismatic shape, dominated by the {100} surface. Hydration of this surface induces multilayer interfacial water structures with distinct OH orientations, which agrees with recent crystal truncation rod measurements. The exchange kinetics between two adjacent water layers exhibits a wide range of underlying timescales (5–180 picoseconds), dictated by ion–water electrostatics. Adsorption of a bidentate hydroxamate ligand reveals that {100} xenotime surface can only accommodate monodentate coordination with water exchange kinetics strongly depending on specific ligand orientation, prompting us to reconsider traditional strategies for selective separation of rare-earth minerals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

One-electron redox kinetics of aqueous transition metal couples Zn 2+/+ , Co 2+/+ , and Ni 2+/+ using pulse radiolysis

The one-electron redox potentials for aqueous metal couples Co 2+/+ and Ni 2+/+ have been investigated by using pulse radiolysis using their reactions with a series of reference compounds to establish the most positive upper limits of E 0 . Experiments with Zn + were also carried out to confirm the characteristic shape of the expected reduction kinetics. Both formate ions and t-BuOH were employed to scavenge ˙OH radicals and ˙H atoms. Kinetics and fitted first and second order reaction rates have been reported for reactions with methyl viologen, fluorescein, Ru(NH 3 ) 6 3+ , Co(en) 3 2+ , Co(sepulcrate) 3+ , Ru(bpy) 3 2+ , Cr(bpy)33+, and Ni(Me 6 [14]4,11-dieneN4) 2+ . Previous work demonstrated that both Co 2+ and Ni 2+ can be reduced by CO 2 ˙- radicals, giving a negative E 0 limit of -1.9 V vs. SHE. A definite reaction of Ni + with fluorescein di-anions provides a new upper limit of the Ni 2+/+ couple as -0.906 V vs. SHE. The reaction of Co + with Ru(bpy) 3 2+ has been confirmed, giving E 0 = -1.3 V vs. SHE as a rigorous upper limit of the Co 2+/+ couple. In the case of Co 2+/+ , kinetics were complicated by a self-catalyzed metal clustering phenomenon. Initiation rate constants of this process have also been reported.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

Kinetic origins of the metastable zone width in the manganese oxide Pourbaix diagram

Pourbaix diagrams show solid–aqueous phase stability as a function of pH and redox potential and can be a valuable tool to guide the hydrothermal synthesis of transition metal oxides. However, the Pourbaix diagram is based on thermodynamics, and nucleation kinetics are not readily apparent in this framework. Here, we conduct a combined experimental and theoretical study to measure the onset of MnO 2 precipitation from an MnO 4 - (aq) solution at various pH conditions, which we then compare against Pourbaix diagram phase boundaries. Using a combination of in situ X-ray absorption spectroscopy and X-ray wide-angle scattering, we directly observe the transformation kinetics from the tetrahedral MnO 4 - (aq) ion to octahedrally coordinated Mn pre-nuclei, as well as its initiation into crystalline δ'-MnO 2 at various pH. The kinetics of octahedral-Mn precursor availability is observed to govern induction times of crystalline MnO 2 nucleation, which results in a metastable zone width in the Mn–H 2 O Pourbaix diagram. These results suggest that synthesis conditions often need to be prepared far beyond the phase boundaries in a Pourbaix diagram to initiate crystallization within a reasonable timescale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Stereoisomer-dependent unimolecular kinetics of 2,4-dimethyloxetanyl peroxy radicals

2,4,dimethyloxetane is an important cyclic ether intermediate that is produced from hydroperoxyalkyl (QOOH) radicals in the low-temperature combustion of n-pentane. However, the reaction mechanisms and rates of consumption pathways remain unclear. In the present work, the pressure- and temperature-dependent kinetics of seven cyclic ether peroxy radicals, which stem from 2,4,dimethyloxetane via H-abstraction and O 2 addition, were determined. The automated kinetic workflow code, KinBot, was used to model the complexity of the chemistry in a stereochemically resolved manner and solve the resulting master equations from 300–1000 K and from 0.01–100 atm. The main conclusions from the calculations include (i) diastereomeric cyclic ether peroxy radicals show significantly different reactivities, (ii) the stereochemistry of the peroxy radical determines which QOOH isomerization steps are possible, (iii) conventional QOOH decomposition pathways, such as cyclic ether formation and HO 2 elimination, compete with ring-opening reactions, which primarily produce OH radicals, the outcome of which is sensitive to stereochemistry. Ring-opening reactions lead to unique products, such as unsaturated, acyclic peroxy radicals, that form direct connections with species present in other chemical kinetics mechanisms through "cross-over" reactions that may complicate the interpretation of experimental results from combustion of n-pentane and, by extension, other alkanes. For example, one cross-over reaction involving 1-hydroperoxy-4-pentanone-2-yl produces 2-(hydroperoxymethyl)-3-butanone-1-yl, which is an iso-pentane-derived ketohydroperoxide (KHP). At atmospheric pressure, the rate of chemical reactions of all seven peroxy radicals compete with that of collisional stabilization, resulting in well-skipping reactions. However, at 100 atm, only one out of seven peroxy radicals undergoes significant well-skipping reactions. Here, the rates produced from the master equation calculations provide the first foundation for the development of detailed sub-mechanisms for cyclic ether intermediates. In addition, analysis of the complex reaction mechanisms of 2,4-dimethyloxetane-derived peroxy radicals provides insights into the effects of stereoisomers on reaction pathways and product yields.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Kinetic processes of phase separation and aggregation behaviors in slot-die processed high efficiency Y6-based organic solar cells

Morphological optimization has proved to be one crucial factor contributing to the 19% efficiency of Y6-based organic solar cells (OSCs). Although the relationship between component miscibility and film morphology has been established, it has not been clarified how the film formation processes proceed, especially for methods compatible with large-area fabrication, i.e. slot-die coating, which restricts the further optimization of the morphology. Herein, we comprehensively investigate the effect of miscibility and the film-formation kinetic process on film morphology during processing with different solvents. A highest power conversion efficiency (PCE) of 17.38% can be obtained in a D18:Y6 device processed with slot-die coating in an open-air environment. However, due to their relatively poor miscibility, incorporating a trace amount of chlorobenzene (CB) into CF can sufficiently promote Y6 aggregation, leading to increased phase separation and thus the average PCE drops to 15.16%. In contrast, a PM6:Y6 blend shows insensitive changes in the kinetic process as well as the final morphology, and thus in the comparable PCEs when cast with different solvents, which is ascribed to their relatively good miscibility. As a result, this work provides scientific guidelines for device optimization by combining the intermediate gap from miscibility to kinetic process and their impact on final morphology.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Quantifying how the cis/trans ratio of N,N -dimethyl-3,5-dimethylpiperidinium hydroxide impacts the growth kinetics, composition and local structure of SSZ-39

This work integrates experiments and computational methods to quantify how the cis/trans ratio of the OSDA used in SSZ-39 synthesis impacts the crystallization kinetics, material properties, and final product composition. The crystallization kinetics increase by 30% when increasing the trans isomer content from 14% to 80%. Per prior work, in all cases based on the synthesis gel composition and product yield aluminum is the limiting reagent, and the absence of any amorphous material detected in the time resolved PXRD studies leads us to conclude that FAU dissolution is the rate limiting step in the formation of SSZ-39 in this synthesis protocol. The TGA and NMR results suggest that the trans isomer of OSDA is selectively incorporated into the product. The NMR binding studies, and corresponding DFT-based results show that the trans isomer binds to FAU more strongly than the cis isomer, providing one possible explanation for this enhancement in kinetics and preferential uptake of the trans isomer. The EDS analysis indicates that the Si/Al ratios are between 7.7 and 8.6 at low and high trans OSDA content, indicating zeolite composition is mildly sensitive to the trans isomer content. EDS results show this decrease in aluminum content leads to a corresponding decrease in sodium uptake. DFT-based calculations confirm OSDA–sodium interactions cannot explain any decrease in sodium uptake, reinforcing lower aluminum content as the cause of lower sodium uptake. Preliminary cobalt titration experiments show a surprisingly low cobalt uptake but also show a clear dependence of the cobalt uptake on the solution pH.

Cui, Zheng [Tulane University, New Orleans, LA (Un↗
DOE OSTI2025

Dynamics of inverse metal oxides on metal catalysts using spectro-kinetics: reversible Brønsted acid site formation and irreversible reduction

Brønsted acid sites (BASs) in inverse catalysts are vital for the selective hydrogenolysis of polyols, specifically cleaving secondary C–O bonds. These BASs form dynamically in situ in an H 2 environment. While H 2 enables rapid BAS generation on short timescales, it reduces the catalyst at prolonged exposures. The active center for BAS generation, the kinetics of BAS formation, its reverse decomposition, and the irreversible oxide reduction have lacked direct experimental evidence. Here, aided by advanced spectro-kinetic studies, we identify trimeric W 3 O x sites on Pt as the active centers for BAS generation, whereas isolated WO x species on SiO 2 act merely as spectator species, demonstrated using an inverse WO x /Pt catalyst as a representative system. A detailed kinetic profile capturing the dynamics of W 3 O x sites on Pt is also established. The rate constant for BAS formation is two orders of magnitude higher than for its decomposition, which is one order of magnitude faster than the irreversible site reduction. Co-fed H 2 O suppresses the site reduction by ∼50%. Furthermore, the H 2 partial pressure plays an important role. While lower gas-phase H 2 partial pressure does not influence the reversible BAS formation, it can significantly (∼3×) suppress catalyst reduction. Finally, these findings offer critical insights into optimizing reaction conditions through periodic H 2 pulsing, enhancing catalyst stability and performance in hydrogenolysis reactions.

Sourav, Sagar [Indian Inst. of Technology (IIT), M↗
DOE OSTI2025

Probing the proton exchange kinetics of BaZr 0.1 Ce 0.7 Y 0.1 Yb 0.1 O 3− δ ceramic electrolyte by operando diffuse reflectance infrared Fourier transform spectroscopy

Proton exchange kinetics plays an important role in governing the performance of intermediate-temperature protonic ceramic electrolysis cells (PCECs) for hydrogen production. Our understanding of the nature of the surface hydration reaction at the single-cell level, however, remains very limited, hampering further efficiency improvements. Here, in this study, we developed a custom operando diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) platform that operates under high temperature and steam conditions with applied bias. Quantitative investigations of surface H 2 O/D 2 O isotope exchange in a BaZr 0.1 Ce 0.7 Y 0.1 Yb 0.1 O 3−δ (BZCYYb1711) protonic electrolyte-based single cell were conducted under different applied voltages using this DRIFTS platform, to gain molecular-level insight into hydration kinetics. The findings show that the application of an external voltage significantly enhances the surface proton exchange rate, decreasing the apparent activation energy from 29.1 kJ mol −1 at open-circuit voltage (OCV) to 6.8 kJ mol −1 at 1.3 V. In addition, distinct voltage-induced spectral shifts in O–D vibrations point to dynamic changes in surface hydration. These findings demonstrate a sensitive spectroscopic platform for probing interfacial proton processes and reveal strong electrochemical control over surface proton kinetics, offering new opportunities for probing electrolyte hydration behavior in PCECs.

36 - MATERIALS SCIENCE↗
DOE OSTI2026

Solubility-limited depolymerization kinetics in the glycolysis of carbonyl-containing polymers

Chemical recycling of condensation polymers is often rationalized on the basis of the intrinsic reactivity of ester and carbonate functional groups. However, under heterogeneous conditions relevant to plastic waste processing and environmental degradation, bulk depolymerization rates often diverge from trends predicted by homogeneous chemistry. Here, we investigate how polymer–solvent compatibility, catalyst strength, and phase behavior govern the heterogeneous glycolysis of carbonyl-containing polymers. Using poly(ethylene terephthalate) (PET), glycol-modified PET (PETG), and bisphenol-A polycarbonate (PC) as model systems, we examine depolymerization kinetics at 180 °C with ethylene glycol and bisphenol A as diols under both amphoteric organosalt (TBD : MSA) and strong base (TBD) catalysis. Despite substantial differences in crystallinity and glycol uptake, PET and PETG depolymerize at comparable rates under organosalt catalysis, while PC depolymerizes significantly more slowly under identical conditions. Time-resolved molecular weight analysis and thermal characterization demonstrate that these rate differences do not arise from crystallinity, swelling, or inherent carbonyl reactivity, but instead reflect solubility-limited kinetics that constrain the transition from heterogeneous to homogeneous reaction regimes. When polymer solubility is low, depolymerization remains heterogeneous and slow; when solubility is enhanced—either through increased polymer–diol compatibility or stronger base catalysis—rapid homogeneous depolymerization is observed, reversing apparent reactivity trends. These results establish solubility and phase behavior as primary determinants of depolymerization kinetics in heterogeneous polymer recycling systems. By demonstrating how catalyst selection and solvent compatibility can expose or overcome solubility limitations, this work provides mechanistic insight to design more energy-efficient and selective chemical recycling processes. More broadly, these findings suggest that polymers with limited solvent or water compatibility may resist chemical degradation in the environment, favoring fragmentation and persistence as micro- and nanoplastics. Understanding solubility-controlled depolymerization offers a pathway toward more sustainable polymer design and end-of-life chemical recovery.

Watson-Sanders, Shelby [Department of Chemistry, U↗
DOE OSTI2017

Kinetic-scale flux rope reconnection in periodic and line-tied geometries

The collisionless reconnection of two parallel flux ropes driven by both the coalescence and kink instabilities is examined using fully kinetic simulations in periodic and line-tied geometries. The three-dimensional reconnection rate is computed from the maximum of the quasi-potential, Ξ≡−∫E·dℓ, where the integral of the electric field is taken along the magnetic field lines across the system. In periodic simulations in which the kink mode is nearly suppressed, reconnection is driven by the coalescence instability, and the peak rate is within 3%–8% of comparable 2D simulations. When a strong kink growth is observed, the peak reconnection rate drops by 10%–25%, and there is a larger drop for lower guide field. With line-tied boundary conditions, the kink instability plays a key role in allowing the flux ropes to interact and partially reconnect. In this limit, the field lines with maximum quasi-potential are associated with a quasi-separatrix layer, and the electric field along these special field lines is supported predominantly by the divergence of the electron pressure tensor. Both of these features, along with the observed reconnection rate, are consistent with recent laboratory experiments on kinetic-scale flux ropes. In kinetic simulations, the non-gyrotropic pressure tensor terms contribute significantly more to the reconnecting electric field than do the gyrotropic terms, while contributions from the electron inertia are significant for field lines adjacent to the quasi-separatrix layer.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Nonlinear saturation of the slab ITG instability and zonal flow generation with fully kinetic ions

Fully kinetic turbulence models are of interest for their potential to validate or replace gyrokinetic models in plasma regimes where the gyrokinetic expansion parameters are marginal. Here, we demonstrate fully kinetic ion capability by simulating the growth and nonlinear saturation of the ion-temperature-gradient instability in shearless slab geometry assuming adiabatic electrons and including zonal flow dynamics. The ion trajectories are integrated using the Lorentz force, and the cyclotron motion is fully resolved. Linear growth and nonlinear saturation characteristics show excellent agreement with analogous gyrokinetic simulations across a wide range of parameters. The fully kinetic simulation accurately reproduces the nonlinearly generated zonal flow. In conclusion, this work demonstrates nonlinear capability, resolution of weak gradient drive, and zonal flow physics, which are critical aspects of modeling plasma turbulence with full ion dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2018

Kinetic simulation of magnetic field generation and collisionless shock formation in expanding laboratory plasmas

Recent laboratory experiments with laser-produced plasmas have observed and studied a number of fundamental physical processes relevant to magnetized astrophysical plasmas, including magnetic reconnection, collisionless shocks, and magnetic field generation by Weibel instability, opening up new experimental platforms for laboratory astrophysics. We develop a fully kinetic simulation model for first-principles simulation of these systems including the dynamics of magnetic fields—magnetic field generation by the Biermann battery effect or Weibel instability; advection by the ion flow, Hall effect, and Nernst effect; and destruction of the field by dissipative mechanisms. Key dimensionless parameters describing the system are derived for scaling between kinetic simulation, recent experiments, and astrophysical plasmas. First, simulations are presented which model Biermann battery magnetic field generation in plasmas expanding from a thin target. Ablation of two neighboring plumes leads to the formation of a current sheet as the opposing Biermann-generated fields collide, modeling recent laser-driven magnetic reconnection experiments. Second, we simulate recent experiments on collisionless magnetized shock generation, by expanding a piston plasma into a pre-magnetized ambient plasma. For parameters considered, the Biermann effect generates additional magnetic fields in the curved shock front and thereby increases shock particle reflection. In conclusion, both cases show the importance of kinetic processes in the interaction of plasmas with magnetic fields and open opportunities to benchmark these important processes through comparison of theory and experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Nonlinear saturation of kinetic ballooning modes by zonal fields in toroidal plasmas

Kinetic ballooning modes (KBMs) are widely believed to play a critical role in disruptive dynamics as well as turbulent transport in magnetic fusion and space plasmas. While the nonlinear evolution of the ballooning modes has been proposed as a mechanism for “detonation” in various scenarios such as the edge localized modes in tokamaks, the role of the kinetic effects in such nonlinear dynamics remains largely unexplored. In this study, global gyrokinetic simulation results of KBM nonlinear behavior are presented. Instead of the finite-time singularity predicted by ideal magnetohydrodynamic theory, the kinetic instability is shown to develop into an intermediate nonlinear regime of exponential growth, followed by a nonlinear saturation regulated by spontaneously generated zonal fields. In the intermediate nonlinear regime, rapid growth of localized current sheets, which can induce magnetic reconnection, is observed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Customizable two-species kinetic equilibria for nonuniform low-beta plasmas

Here, two methods are developed for constructing self-consistent two-species kinetic equilibria for low-beta electrostatic plasmas, in which the magnetic field is uniform and fixed in time. The first method uses Taylor series approximations to construct distribution functions that can be specified analytically. The second method relies on numerically solving a nonlinear ordinary differential equation and produces exact—to numerical precision—equilibria. In both methods, the equilibrium distribution functions for ions and electrons are expressed in terms of constants of motion and satisfy the steady-state Vlasov-Poisson equation system. Provided that the ion drift speed does not exceed the ion thermal speed, the equilibria can be specified with customizable density and electrostatic potential profiles. The methods can thereby be tailored to different applications and are successfully applied to construct kinetic equilibria for cross-field plasmas with sheared flows, large density variations, and different levels of magnetization. The equilibria are used to initialize fourth-order finite-volume Vlasov-Poisson simulations in (x, v x , v y ) coordinates and the associated temporal evolution is used to assess the accuracy of each method. Lastly, the low-amplitude deviations observed in these simulations demonstrate that the kinetic equilibria are robust and that they provide a valuable means of studying the dynamics of nonuniform magnetized plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗