Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “kinetic”

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 955 records · Page 53
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Calculation of self-diffusion coefficients in supercritical carbon dioxide using mean force kinetic theory

This paper presents an application of mean force kinetic theory (MFT) to the calculation of the self-diffusivity of CO 2 in the supercritical fluid regime. Two modifications to the typical application of MFT are employed to allow its application to a system of molecular species. Therefore, the first is the assumption that the inter-particle potential of mean force can be obtained from the molecule center-of-mass pair correlation function, which in the case of CO 2 is the C–C pair correlation function. The second is a new definition of the Enskog factor that describes the effect of correlations at the surface of the collision volume. The new definition retains the physical picture that this quantity represents a local density increase, resulting from particle correlations, relative to that in the zero density homogeneous fluid limit. These calculations are facilitated by the calculation of pair correlation functions from molecular dynamics (MD) simulations using the FEPM2 molecular CO2 model. The self-diffusivity calculated from theory is in good agreement with that from MD simulations up to and slightly beyond the density at the location of the Frenkel line. The calculation is compared with and is found to perform similarly well to other commonly used models but has a greater potential for application to systems of mixed species and to systems of particles with long range interatomic potentials due to electrostatic interactions.

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

Kinetic Landau-fluid closures of non-Maxwellian distributions

New kinetic Landau-fluid closures, based on the cutoff Maxwellian distribution, are derived. A special static case is considered (the frequency ω=0). In the strongly collisional regime, our model reduces to Braginskii's heat flux model, and the transport is local. In the weak collisional regime, our model indicates that the heat flux is non-local and recovers the Hammett–Perkins model while the value of the cutoff velocity approaches to infinity. We compare the thermal transport coefficient χ of Maxwellian, cutoff Maxwellian and super-Gaussian distribution. The results show that the reduction of the high-speed tail particles leads to the corresponding reduction of the thermal transport coefficient χ across the entire range of collisionality, more reduction of the free streaming transport toward the weak collisional regime. In the collisionless limit, χ approaches to zero for the cutoff Maxwellian and the super-Gaussian distribution but remains finite for Maxwellian distribution. χ is complex if the cutoff Maxwellian distribution is asymmetric. The Im(χ) approaches to different convergent values in both collisionless and strongly collisional limit, respectively. It yields an additional streaming heat flux in comparison with the symmetric cutoff Maxwellian distribution. Furthermore, due to the asymmetric distribution, there is a background heat flux q0 though there is no perturbation. Furthermore, the derived Landau-fluid closures are general for fluid moment models, and applicable for the cutoff Maxwellian distribution in an open magnetic field line region, such as the scape-off-layer of Tokamak plasmas, in the thermal quench plasmas during a tokamak disruption, and the super-Gaussian electron distribution function due to inverse bremsstrahlung heating in laser-plasma studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Continuum kinetic investigation of the impact of bias potentials in the current saturation regime on sheath formation

In this work, we examine sheath formation in the presence of bias potentials in the current saturation regime for pulsed power fusion experiments. It is important to understand how the particle and heat fluxes at the wall may impact the wall material and affect electrode degradation. Here, simulations are performed using the 1X-1V Boltzmann–Poisson system for a proton–electron plasma in the presence of bias potentials ranging from 0 to 10 kV. The results indicate that the sheath near the anode remains generally the same as that of a classical sheath without the presence of a bias potential. However, the sheath near the cathode becomes more prominent with a larger potential drop, a significant decrease in the electron density, and larger sheath lengths. The spatially constant current density increases to a saturation value with the increasing bias potential. For high bias potentials, the material choice needs to consider that the anode faces significantly larger particle and heat fluxes compared to the cathode. In general, the results trend with theory with differences attributed to the simplified assumptions in the theory and the kinetic effects considered in the simulations. Due to the significant computational cost of a well resolved 1X-2V simulation, only one such simulation is performed for the 5 kV case showing higher current.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

An extended Vlasov–Fokker–Planck approach for kinetic simulations of laser plasmas

Vlasov–Fokker–Planck simulation codes occupy an important niche in modeling laser-produced plasmas, since they are well suited to studying the effect of collisions on electron kinetic phenomena, especially energy transport. One of the most important elements of energy transport is the absorption of laser light by the plasma; however, simulating this in detail requires resolving oscillations of the laser light, whose characteristic timescale is orders of magnitude shorter than the simulation time needed to study transport physics. For this reason, most Vlasov–Fokker–Planck codes used to study electron transport in laser plasmas rely on simplified models of the laser–plasma coupling. Their underlying assumptions nominally preclude their use for modeling laser light having short-scale structure in space or time, such as broadband lasers. In this work, we derive a more general computational framework suitable for arbitrarily structured laser fields. Furthermore, our approach is based on an extended set of Vlasov–Fokker–Planck equations that separately solve for the low- and high-frequency plasma response. We implement these extended Vlasov–Fokker–Planck equations in the spherical harmonic code K2 and demonstrate the performance of the method on several laser absorption test problems, with particular attention to the judicious selection of time steps, time integrators, and spherical harmonic truncation, according to the intensity and spectrum of the laser light under consideration. Comparison with the widely used Langdon absorption operator shows the Langdon operator performs remarkably well for predicting laser heating in the simple cases considered here, even in situations that would seem to violate its underlying assumptions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Probabilistic and maximum entropy modeling of chemical reaction systems: Characteristics and comparisons to mass action kinetic models

We demonstrate and characterize a first-principles approach to modeling the mass action dynamics of metabolism. Starting from a basic definition of entropy expressed as a multinomial probability density using Boltzmann probabilities with standard chemical potentials, we derive and compare the free energy dissipation and the entropy production rates. We express the relation between entropy production and the chemical master equation for modeling metabolism, which unifies chemical kinetics and chemical thermodynamics. Because prediction uncertainty with respect to parameter variability is frequently a concern with mass action models utilizing rate constants, we compare and contrast the maximum entropy model, which has its own set of rate parameters, to a population of standard mass action models in which the rate constants are randomly chosen. We show that a maximum entropy model is characterized by a high probability of free energy dissipation rate and likewise entropy production rate, relative to other models. We then characterize the variability of the maximum entropy model predictions with respect to uncertainties in parameters (standard free energies of formation) and with respect to ionic strengths typically found in a cell.

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

Similarity for downscaled kinetic simulations of electrostatic plasmas: Reconciling the large system size with small Debye length

A simple similarity has been proposed for kinetic (e.g., particle-in-cell) simulations of plasma transport that can effectively address the long-standing challenge of reconciling the tiny Debye length with the vast system size. This applies to both transport in unmagnetized plasma and parallel transport in magnetized plasmas, where the characteristics length scales are given by the Debye length, collisional mean free paths, and the system or gradient lengths. The controlled scaled variables are the configuration space, x/L, and an artificial Coulomb Logarithm, L ln Λ, for collisions, while the scaled time, t/L, and electric field, LE, are automatic outcomes. The similarity properties are examined, demonstrating that the macroscopic transport physics is preserved through a similarity transformation while keeping the microscopic physics at its original scale of Debye length. To showcase the utility of this approach, two examples of 1D plasma transport problems were simulated using the VPIC code: the plasma thermal quench in tokamaks [Li et al., Nuclear Fusion 63, 066030 (2023)] and the plasma sheath in the high-recycling regime [Li et al., Physics of Plasmas 30, 063505 (2023)].

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Challenges in the study of chemistry and photochemistry at air–water interfaces: Toward in situ monitoring of reaction kinetics with spectroscopic techniques

Air–water interfaces, including those on droplet surfaces, have been the subject of many recent experiments due to their propensity for unique chemistry. Here, in this study, an overview of some recent advancements in understanding interfacial reaction kinetics is provided, highlighting non-surface-specific methods—such as mass spectrometry—compared with the advantages of surface-specific techniques—such as reflection–absorption spectroscopy, sum frequency generation, and photoelectron spectroscopy. This Perspective discusses the information depth and time constraints of common surface-specific spectroscopic methods that need to be addressed to monitor interfacial chemistry in situ and uses a few key examples from the literature as case studies. It concludes by advocating for the continued development of advanced spectroscopic methods to further investigate interfacial chemistry, underscoring the need for interdisciplinary collaboration to bridge the gap between molecular-level insights and macroscopic observations in future research.

Air-water interface↗
DOE OSTI2021

Kinetic proofreading of lipochitooligosaccharides determines signal activation of symbiotic plant receptors

Plant cell surface receptors perceive carbohydrate signaling molecules and hereby establish communication with surrounding microbes. Genetic studies have identified two different classes of lysin motif receptor kinases as gatekeepers that together trigger the symbiotic pathway in plants; however, no structural or functional data of the perception mechanisms switching these receptors from resting state into activation is known. In this study, we use structural biology, biochemical, and genetic approaches to demonstrate how the NFP/NFR5 class of lipochitooligosaccharide (LCO) receptors discriminate bacterial symbionts based on a kinetic proofreading mechanism that controls receptor activation and signaling specificity. We show that the LCO binding site can be engineered to support symbiotic functions, which greatly advance future opportunities for receptor engineering in legumes and nonlegumes.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

In vivo kinetics of SARS-CoV-2 infection and its relationship with a person’s infectiousness

Significance Quantifying the kinetics of SARS-CoV-2 infection and individual infectiousness is important for understanding SARS-CoV-2 transmission and evaluating intervention strategies. Here, we developed within-host models of SARS-CoV-2 infection, and by fitting them to clinical data, we estimated key within-host viral dynamic parameters. We also developed a mechanistic model for viral transmission and show that the logarithm of the viral load in the upper respiratory tract serves as an appropriate surrogate for a person’s infectiousness. Using data on how viral load changes during infection, we further evaluated the effectiveness of PCR and antigen-based testing strategies for averting transmission and identifying infected individuals.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2023

Two-Phase Turbulent Kinetic Energy Budget Computation in Co-Current Taylor Bubble Flow

Simulation of two-phase flows is relevant for reactor design and safety at normal operation or during accident scenarios. Often, the two-phase flow is in a regime in which slugs are formed or where the flow stratifies. Modeling such situations using standard single-phase Reynolds-averaged Navier-Stokes (RANS) turbulence models fails due to an overestimation of the eddy viscosity at the resolved two-phase interface. To solve this, an ad hoc turbulence damping term has been proposed in the literature that reduces the turbulence production locally at a two-phase interface, analogously to turbulence wall functions. However, this approach must be tailored to the specific setting and does not consider physical contributions such as surface tension or flow topology. Therefore, the problem of two-phase interfacial turbulence must be studied more in-depth. In this work, we consider co-current turbulent Taylor bubble flow using high-fidelity numerical simulation. The Basilisk code is used to simulate a Taylor bubble rising in a vertical pipe. By simulating the bubble in a moving frame of reference, we may study the turbulent kinetic energy (TKE) budgets ahead of the bubble, in its wake, and across the interface. The implementation of the TKE budget computation and the underlying averaging techniques are validated for the single-phase region ahead of the Taylor bubble using reference direct numerical simulation data. The analysis of the TKE budgets in the setting of Taylor bubble flow allows for the study of how turbulence behaves due to the presence of a two-phase interface and, in turn, supports the improvement of two-phase RANS models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

A physics-constrained neural ordinary differential equations approach for robust learning of stiff chemical kinetics

The high computational cost associated with solving for detailed chemistry poses a significant challenge for predictive computational fluid dynamics (CFD) simulations of turbulent reacting flows. While deep learning techniques have been explored to develop faster surrogate models, they often fail to integrate reliably with CFD solvers. This instability arises because traditional deep learning approaches optimize for training error without ensuring compatibility with ordinary differential equation (ODE) solvers, resulting in accumulation of errors over time. Recently, neuralODE (NODE) based approaches have been shown to be a promising technique to emulate and accelerate detailed chemistry computations. Here, in the present work, we extend this NODE framework for stiff chemical kinetics by incorporating mass conservation constraints directly into the loss function during training. This ensures that the total mass as well as the individual elemental species masses are conserved in an a-posteriori manner. Proof-of-concept studies are performed with the novel physics-constrained NODE (PC-NODE) approach for homogeneous autoignition of hydrogen-air mixture over a range of composition and thermodynamic conditions. It is demonstrated that the PC-NODE framework not only improves the physical consistency of the resulting data-driven model with respect to mass conservation criteria, but also improves training efficiency. PC-NODE is shown to achieve 2–100× speedup relative to the hydrogen-air detailed chemical mechanism depending on the type of the ODE solver (implicit or explicit) used during autoregressive inference tests. Lastly, a-posteriori studies are performed wherein the trained PC-NODE model is coupled with a CFD solver. It is shown that higher accuracy is achieved with PC-NODE relative to the purely data-driven NODE approach. Moreover, PC-NODE also exhibits robustness and generalizability to unseen initial conditions from within (interpolative capability) as well as outside (extrapolative capability) the training regime.

computational combustion↗
DOE OSTI2023

Sulfidation kinetics of titanium and Ti-6Al-4V with elemental sulfur

The sulfidation kinetics of titanium and Ti-6Al-4V (Ti64) using elemental sulfur were studied by measuring the amount of metal consumed instead of the change in mass or scale thickness. Here this metal consumption approach combines the parabolic oxidation law, Arrhenius equation, and Pilling-Bedworth ratio to determine the apparent activation energy and the pre-exponential factor of consumption of Ti64 during sulfidation. The estimations of apparent activation energy and pre-exponential factor were 110.4 kJ·mol -1 and 1.3×10 -9 m 2 s -1 for the titanium-sulfur, and 116.2 kJ·mol -1 and 1.2×10 -9 m2s -1 for the Ti64-sulfur system. Finally, a material removal predictor was developed to predict the amount of Ti64 consumed for a given sulfidation temperature and time.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Sheath collapse at critical shallow angle due to kinetic effects

The Debye sheath is known to vanish completely in magnetised plasmas for a sufficiently small electron gyroradius and small angle between the magnetic field and the wall. This angle depends on the current onto the wall. When the Debye sheath vanishes, there is still a potential drop between the wall and the plasma across the magnetic presheath. The magnetic field angle corresponding to the predicted sheath collapse is shown to be much smaller than previous estimates, scaling with the electron-ion mass ratio and not with the square root of the mass ratio. This is shown to be a consequence of the kinetic electron and finite ion orbit width effects, which are not captured by fluid models. The wall potential with respect to the bulk plasma at which the Debye sheath vanishes is calculated. Above this wall potential, it is possible that the Debye sheath will invert.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Drift kinetic theory of the NTM magnetic islands in a finite beta general geometry tokamak plasma

In (Imada et al 2019 Nucl. Fusion 59 046016 and references therein) a new 4D drift kinetic nonlinear theory, valid in the limit of a low beta, small inverse aspect ratio, circular cross section, toroidal geometry, to describe the plasma response to the neoclassical tearing mode (NTM) magnetic perturbation is derived. In (Dudkovskaia et al 2021 Plasma Phys. Control. Fusion 63 054001) this theory is reduced in a low collisionality limit, which allows a dimensionality reduction to a 3D problem to efficiently resolve the collisional dissipation layer in the vicinity of the trapped-passing boundary. (Dudkovskaia et al 2021 Plasma Phys. Control. Fusion 63 054001) adopts an improved model for the magnetic drift frequency, which reduces the threshold magnetic island half-width from $8.73\rho_{b i}$, where $\rho_{b i}$ is the trapped ion banana orbit width, to $1.46 \rho_{b i}$, making it in closer agreement with experimental observations for the large aspect ratio tokamak equilibrium. In the present paper, the theory is extended to a high beta, arbitrary tokamak geometry to capture the plasma shaping effects on the NTM threshold physics with the focus on the non-zero triangularity discharges that are known to have a strong impact on the plasma MHD stability. First, it is found that the higher triangularity plasma is more prone to NTMs which is in agreement with the 2/1 tearing mode onset relative frequency measurements in DIII-D. Second, the NTM threshold dependence on the tokamak inverse aspect ratio obtained in (Dudkovskaia et al 2021 Plasma Phys. Control. Fusion 63 054001) is refined and extended to a finite aspect ratio limit. Third, the NTM threshold dependence on poloidal beta is obtained and benchmarked against the EAST threshold island width measurements.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Viral Load Kinetics of Severe Acute Respiratory Syndrome Coronavirus 2 in Hospitalized Individuals With Coronavirus Disease 2019

Abstract Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) kinetics remain understudied, including the impact of remdesivir. In hospitalized individuals, peak sputum viral load occurred in week 2 of symptoms, whereas viremia peaked within 1 week of symptom-onset, suggesting early systemic seeding of SARS-CoV-2. Remdesivir treatment was associated with faster viral decay.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2024

Modeling stochastic fluctuations in relativistic kinetic theory

Using the information current, we develop a Lorentz-covariant framework for modeling equilibrium fluctuations in relativistic kinetic theory in the grand-canonical ensemble. The resulting stochastic theory is proven to be causal and covariantly stable, and its predictions do not depend on the choice of spacetime foliation used to define the grand-canonical probabilities. As expected, in a box containing N > 5 particles, Boltzmann’s molecular chaos postulate is broken with (almost exact) probability N -1/2 , leading to a breakdown of the Boltzmann equation in small systems. Here, we also verify that, in ultrarelativistic gases, transient hydrodynamics already accounts for at least 80% of the equilibrium fluctuations of the stress-energy tensor at a given time. Finally, we compute the correlators at nonequal times for two selected collision kernels: that of a chemically active diluted solution, and that of ultrarelativistic scalar particles self-interacting via a quartic potential. For the former, we compute the density-density correlators analytically in real space, and dehydrodynamization of the stochastic theory is proven to occur whenever the mean free path diverges at high energy.

Astronomy & Astrophysics↗
DOE OSTI2020

Measurement of Kinetic-Scale Current Filamentation Dynamics and Associated Magnetic Fields in Interpenetrating Plasmas

Here, we present the first local, quantitative measurements of ion current filamentation and magnetic field amplification in interpenetrating plasmas, characterizing the dynamics of the ion Weibel instability. The interaction of a pair of laser-generated, counterpropagating, collisionless, supersonic plasma flows is probed using optical Thomson scattering (TS). Analysis of the TS ion-feature revealed anticorrelated modulations in the density of the two ion streams at the spatial scale of the ion skin depth c/ω pi = 120 μm, and a correlated modulation in the plasma current. The inferred current profile implies a magnetic field amplitude ~30 ± 6 T , corresponding to ~1% of the flow kinetic energy, indicating that magnetic trapping is the dominant saturation mechanism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Stochastic Kinetics of Nanocatalytic Systems

Catalytic reaction events occurring on the surface of a nanoparticle constitute a complex stochastic process. Although advances in modern single-molecule experiments enable direct measurements of individual catalytic turnover events occurring on a segment of a single nanoparticle, we do not yet know how to measure the number of catalytic sites in each segment or how the catalytic turnover counting statistics and the catalytic turnover time distribution are related to the microscopic dynamics of catalytic reactions. Here, we address these issues by presenting a stochastic kinetics for nanoparticle catalytic systems. We propose a new experimental measure of the number of catalytic sites in terms of the mean and variance of the catalytic event count. By considering three types of nanocatalytic systems, we investigate how the mean, the variance, and the distribution of the catalytic turnover time depend on the catalytic reaction dynamics, the heterogeneity of catalytic activity, and communication among catalytic sites. This work enables accurate quantitative analyses of single-molecule experiments for nanocatalytic systems and enzymes with multiple catalytic sites.

Statistical Physics↗