Search NASA⌕ Search

SEARCH · Search NASA

Results for “stochastic 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.

152 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Nonlocal suppression of Biermann battery magnetic-field generation for arbitrary atomic numbers and magnetization

The Biermann battery term of magnetohydrodynamics (MHD) generates a magnetic field where electron density gradients and electron temperature gradients are perpendicular to one another. Kinetic simulations and experiments have shown that the rate of magnetic-field generation is lower than Biermann when the electron mean free path becomes comparable to or greater than the temperature gradient scale length, known as the nonlocal regime. We investigate the nonlocal suppression of the Biermann term using simplified Fokker–Planck simulations covering a wide range of parameters. We provide the first fit for nonlocal Biermann suppression that has physically accurate behavior for small and large values of a suitable nonlocality parameter, valid for an arbitrary atomic number, and that includes the effect of magnetization on nonlocality. Lastly, the fit is intended to provide an approximate method to account for reduced magnetic-field generation in MHD codes and theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A physics-constrained deep learning treatment of runaway electron dynamics

An adjoint formulation leveraging a physics-informed neural network (PINN) is employed to advance the density moment of a runaway electron (RE) distribution forward in time. A distinguishing feature of this approach is that once the adjoint problem is solved, its solution can be used to project the RE density forward in time for an arbitrary initial momentum space distribution of REs. Furthermore, by employing a PINN, a parametric solution to the adjoint problem can be learned. Thus, once trained, this adjoint-deep learning framework is able to efficiently project the RE density forward in time across various plasma conditions while still including a fully kinetic description of RE dynamics. As an example application, the temporal evolution of the density of primary electrons is studied, with particular emphasis on evaluating the decay of a RE population when below threshold. Predictions from the adjoint-deep learning framework are found to be in good agreement with a traditional relativistic electron Fokker–Planck solver, for several distinct initial conditions, and across an array of physics parameters. Once trained, the PINN thus provides a means of generating RE density time histories with exceptionally low online execution time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Unveiling the catalytic mechanism of GTP hydrolysis in microtubules

Microtubules (MTs) are large cytoskeletal polymers, composed of αβ-tubulin heterodimers, capable of stochastically converting from polymerizing to depolymerizing states and vice versa. Depolymerization is coupled with hydrolysis of guanosine triphosphate (GTP) within β-tubulin. Hydrolysis is favored in the MT lattice compared to a free heterodimer with an experimentally observed rate increase of 500- to 700-fold, corresponding to an energetic barrier lowering of 3.8 to 4.0 kcal/mol. Mutagenesis studies have implicated α-tubulin residues, α:E254 and α:D251, as catalytic residues completing the β-tubulin active site of the lower heterodimer in the MT lattice. The mechanism for GTP hydrolysis in the free heterodimer, however, is not understood. Additionally, there has been debate concerning whether the GTP-state lattice is expanded or compacted relative to the GDP state and whether a “compacted” GDP-state lattice is required for hydrolysis. In this work, extensive quantum mechanics/molecular mechanics simulations with transition-tempered metadynamics free-energy sampling of compacted and expanded interdimer complexes, as well as a free heterodimer, have been carried out to provide clear insight into the GTP hydrolysis mechanism. α:E254 was found to be the catalytic residue in a compacted lattice, while in the expanded lattice, disruption of a key salt bridge interaction renders α:E254 less effective. The simulations reveal a barrier decrease of 3.8 ± 0.5 kcal/mol for the compacted lattice compared to a free heterodimer, in good agreement with experimental kinetic measurements. Additionally, the expanded lattice barrier was found to be 6.3 ± 0.5 kcal/mol higher than compacted, demonstrating that GTP hydrolysis is variable with lattice state and slower at the MT tip.

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

Non-local transport in radiation-hydrodynamics codes for ICF by efficient coupling to an external Vlasov–Fokker–Planck code

Accurately incorporating non-local transport into radiation-hydrodynamics codes, and indeed any fluid system, has long been elusive. To date, a simplified and accurate theory that can be easily integrated has not been available. This limitation affects modeling in inertial confinement fusion (ICF) and magnetic confinement fusion systems, among others, where non-local transport is well-known to be present. Here, we present a coupling methodology between a full Vlasov–Fokker–Planck (VFP) electron kinetic code and radiation-hydrodynamics (rad-hydro) codes. The VFP code is used to adjust native electron transport in the rad-hydro code, thus enabling improved transport without the need to integrate a full electron VFP solver into the rad-hydro code. This approach necessitates only occasional invocation of the VFP code, reducing computational intensity compared to following the dynamic evolution entirely with the VFP code on fluid time scales. We illustrate that the methodology is more accurate than other simplified methods in thermal decay systems relevant to ICF and can replicate standard theoretical results with high accuracy.

Electronic transport↗
DOE OSTI2020

Magnetization around mix jets entering inertial confinement fusion fuel

Engineering features are known to cause jets of ablator materials to enter the fuel hot-spot in inertial confinement fusion implosions. The Biermann battery mechanism wraps them in a self-generated magnetic field. Here we show that higher-Z jets have an additional thermoelectric magnetic source term that is not present for hydrogen jets, verified here through a kinetic simulation. It has similar magnitude to the Biermann term. We then include this in an extended magneto-hydrodynamic approach to post process an xRAGE radiation-hydrodynamic implosion simulation. The simulation includes an accurate model for the capsule fill tube, producing a dense carbon jet that becomes wrapped in a 4000 T magnetic field. A simple spherical carbon mix model shows that this insulates the electron heat conduction enough to cause contraction of the jet to an optically thick equilibrium. The denser magnetized jet hydrodynamics could change its core penetration and therefore the final mix mass, which is known to be well correlated with fusion yield degradation. Fully exploring this will require self-consistent magneto-hydrodynamic simulations. Experimental signatures of this self-magnetization may emerge in the high energy neutron spectrum.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Flame acceleration and transition to detonation in a pre-/main-chamber combustion system

Numerical simulations are performed to study the mechanism of deflagration to detonation transition (DDT) in a pre-/main-chamber combustion system with a stoichiometric ethylene–oxygen mixture. A Godunov algorithm, fifth-order in space, and third-order in time, is used to solve the fully compressible Navier–Stokes equations on a dynamically adapting mesh. A single-step, calibrated chemical diffusive model described by Arrhenius kinetics is used for energy release and conservation between the fuel and the product. The two-dimensional simulation shows that a laminar flame grows in the pre-chamber and then develops into a jet flame as it passes through the orifice. A strong shock forms immediately ahead of the flame, reflecting off the walls and interacting with the flame front. The shock–flame interactions are crucial for the development of flame instabilities, which trigger the subsequent flame development. The DDT arises due to a shock-focusing mechanism, where multiple shocks collide at the flame front. A chemical explosive mode analysis (CEMA) criterion is developed to study the DDT ignition mode. Preliminary one-dimensional computations for a laminar propagating flame, a fast flame deflagration, and a Chapman–Jouguet detonation are conducted to demonstrate the validity of CEMA on the chemical-diffusive model, as well as to determine the proper conditioning value for CEMA diagnostic. The two-dimensional analysis with CEMA indicates that the DDT initiated by the shock-focusing mechanism can form a strong thermal expansion region at the flame front that features large positive eigenvalues for the chemical explosive mode and dominance of the local autoignition mode. Thus, the CEMA criterion proposed in this study provides a robust diagnostic for identifying autoignition-supported DDT, of which the emergence of excessive local autoignition mode is found to be a precursor. The effect of grid size, initial temperature, and orifice size are then evaluated, and results show that although the close-chamber DDT is highly stochastic, the detonation initiation mechanism remains robust.

42 ENGINEERING↗
DOE OSTI2026

Generalized master equation for particle transport in binary random media with renewal statistics

Particle transport in binary stochastic mixtures is classically modeled assuming Markovian or exponential mixing statistics but in many applications material memory invalidates the Markov assumption. For non-Markovian mixing characterized by alternating renewal processes, a transport-theoretic framework is presented that provides an exact description of transport in nonscattering random binary media with general non-exponential statistics. Our approach is to Markovianize the problem by augmenting the {material type, particle flux} state space with the age or distance from the last interface. A Chapman-Kolmogorov equation is formulated for the joint probability density of the material type, particle flux, and age, and subsequently reduced to a generalized Master equation (GME) in differential form. This constitutes the primary result of this work. A state-updating Monte Carlo algorithm consistent with the GME is developed and benchmarked against analytical solutions for multiple chord-length laws. For purely absorbing renewal statistical media, the GME reproduces analytical benchmarks for the equilibrium age distribution, interior mean/variance of material-conditioned fluxes, and boundary transmittance. Simulations further demonstrate that a Markov (exponential) approximation of non-exponential statistics can introduce large errors in transmittance and interior flux profiles. Lastly, the reintroduction of memory due to scattering is briefly addressed through heuristic considerations.

Fluctuations & noise↗
DOE OSTI2025

Excitation of whistler and slow-X waves by runaway electrons in a collisional plasma

Runaway electrons are known to provide robust ideal or collisionless kinetic drive for plasma wave instabilities in both the whistler and slow-X branches, via the anomalous Doppler-shifted cyclotron resonances. In a cold and dense post-thermal-quench plasma, collisional damping of the plasma waves can compete with the collisionless drive. Previous studies have found that, due to their higher wavelength and frequency, slow-X waves suffer stronger collisional damping than the whistlers, while the ideal growth rate of slow-X modes is higher. Here, we study runaway avalanche distributions that maintain the same eigen distribution and increase only in magnitude over time. The distributions are computed from the relativistic Fokker–Planck–Boltzmann solver, upon which a linear dispersion analysis is performed to search for the most unstable or least damped slow-X and whistler modes. Taking into account the effect of plasma density, plasma temperature, and effective charge number, we find that the slow-X modes tend to be excited before the whistlers in a runaway current ramp-up. Furthermore, even when the runaway current density is sufficiently high that both branches are excited, the most unstable slow-X mode has a much higher growth rate than the most unstable whistler mode. The qualitative and quantitative trends uncovered in the current study indicate that even though past experiments and modeling efforts have concentrated on whistler modes, there is a compelling case that slow-X modes should also be a key area of focus in the runaway self-mediation through wave instabilities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗