Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “KINETIC THEORY”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Plasma hydrodynamics from mean force kinetic theory

Mean force kinetic theory is used to evaluate the electrical conductivity, thermal conductivity, electrothermal coefficient, thermoelectric coefficient, and shear viscosity of a two-component plasma. The results are compared with molecular dynamics simulations. These simulations are made possible by assuming a repulsive Coulomb force for all interactions. Good agreement is found for all coefficients up to a Coulomb coupling strength of Γ ≈ 20. This is over 100-times larger than the coupling strength at which traditional theories break down. It is concluded that mean force kinetic theory is a promising framework for extending hydrodynamics to dense plasmas.

Electrical conductivity↗
DOE OSTI2024

Relativistic Guiding-Center Motion: Action Principle, Kinetic Theory and Hydrodynamics

We treat the guiding-center dynamics in a varying external Maxwell field using a relativistically covariant action principle which reproduces the known Vandervoort expression for the drift velocity and extends it to curved spacetime. We derive the corresponding kinetic theory and ideal hydrodynamic theory. In contrast to conventional five-equation hydrodynamics, the guiding-center hydrodynamics needs only three equations due to a constraint on the motion across magnetic field. Furthermore, we argue that such a hydrodynamics is applicable to strongly coupled plasmas where kinetic theory fails.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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 OSTI2025

Thermal relaxation and the complete set of second-order transport coefficients for the unitary Fermi gas from kinetic theory

We compute the complete set of second-order transport coefficients of the unitary Fermi gas, a dilute gas of spin-1/2 particles interacting via an 𝑠 -wave interaction tuned to infinite scattering length. The calculation is based on kinetic theory and the Chapman-Enskog method at second order in the Knudsen expansion. We take into account the exact two-body collision integral. We extend previous results on second-order coefficients related to shear stress by including terms related to heat flow and gradients of the fugacity. We confirm that the thermal relaxation time is given by the simple estimate 𝜏 𝜅 = 𝜅⁢𝑚/(𝑐 𝑃 ⁢𝑇) even if the full collision kernel is taken into account. Furthermore, 𝜅 is the thermal conductivity, 𝑚 is the mass of the particles, 𝑐𝑃 is the specific heat at constant pressure, and 𝑇 is the temperature.

Kinetic theory↗
DOE OSTI2024

Formulation of a one-dimensional electrostatic plasma model for testing the validity of kinetic theory

Here, we present a one-dimensional (1-D) model composed of aligned, electrostatically interacting charged disks, conceived to address in a computable model the validity of the Bogoliubov assumption on the decay of particle correlations in the Born–Bogoliubov–Green–Kirkwood–Yvon hierarchy. This assumption is a basic premise of plasma kinetic theory. The disk model exhibits spatially 1-D features at short distances, but retains 3-D features at large distances. Here the collective dynamics of this model plasma is investigated by solving the corresponding Vlasov equation. In addition, the implementation of the model for the numerical validation of the Bogoliubov assumption is formulated.

1-D plasma model↗
DOE OSTI2026

Expansion-driven Weibel instability in magnetohydrodynamics: Linear theory for static systems

The Weibel dispersion relation is obtained in the magnetohydrodynamic (MHD) approximation by including a tensor expression for temperature. MHD gives an upper cutoff wavenumber for Weibel growth identical to kinetic theory if the electron drift velocity is included in the evolution of the electron temperature. MHD overestimates growth rates compared to kinetic theory and gives maximum growth at larger wavenumbers, but it only leads to divergent results as temperature anisotropy tends to infinity. Thermal conduction in MHD lowers the growth rates and shifts maximum growth to smaller wavenumbers, but flux-limited or nonlocal thermal conduction is found to significantly limit these effects. The results are used as a first step in evaluating the potential of MHD to simulate expansion-driven Weibel instability, which occurs due to the temperature tensor only cooling in the direction of expansion. Finally, the small temperature anisotropy found in expansion-driven Weibel instability means that MHD could be an adequate model because growth rates are much less than the electron plasma frequency, and wavelengths are much greater than the Debye length.

Davies, J. R. [Univ. of Rochester, NY (United Stat↗
DOE OSTI2025

A geometric perspective on kinetic matter–radiation interaction and moment systems

Here, in this paper, we provide a geometric perspective on the kinetic interaction of matter and radiation, based on a pair bracket approach. We discuss the interaction of kinetic theories via dissipative brackets, with our fundamental example being the coupling of matter, described by the Boltzmann equation, and radiation, described by the radiation transport equation. We explore the transition from kinetic systems to their corresponding moment systems, provide a Hamiltonian description of such moment systems, and give a geometric interpretation of the moment closure problem for kinetic theories. As an application, we discuss in detail diffusion radiation hydrodynamics as an example of a pair bracket formulation on a space of moments corresponding to kinetic matter–radiation interaction. Additionally, using the variable moment closure framework of Burby [J. W. Burby, Variable-moment fluid closures with Hamiltonian structure, Sci. Rep. 13 (2023) 18286, doi:10.1038/s41598-023-45416-5], we show how to construct Hamiltonian moment closures for kinetic transport equations with arbitrary Hamiltonian. Using this general construction, we derive novel Hamiltonian moment closures for pure radiation transport.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Experimental and computational study of phase space dynamics in strongly coupled plasmas with steep density gradients

Understanding how plasmas thermalize when density gradients are steep remains a fundamental challenge in plasma physics, with direct implications for fusion experiments and astrophysical phenomena. Standard hydrodynamic models break down in these regimes, and kinetic theories make predictions that have never been directly tested. Here, we present the first detailed phase-space measurements of a strongly coupled plasma as it evolves from sharp density gradients to thermal equilibrium. Using laser-induced fluorescence imaging of an ultracold calcium plasma, we track the complete ion distribution function f(x,v,t). We discover that commonly used kinetic models (Bhatnagar–Gross–Krook and Lenard–Bernstein) overpredict thermalization rates, even while correctly capturing the initial counterstreaming plasma formation. Our measurements reveal that the initial ion acceleration response scales linearly with electron temperature, and that the simulations underpredict the initial ion response. In our geometry we demonstrate the formation of well-controlled counterpropagating plasma beams. This experimental platform enables precision tests of kinetic theories and opens new possibilities for studying plasma stopping power and flow-induced instabilities in strongly coupled systems.

Bergeson, Scott (ORCID:0000000231249226)↗
DOE OSTI2026

A conservative discontinuous Galerkin algorithm for particle kinetics on smooth manifolds

A novel, conservative discontinuous Galerkin algorithm is presented for particle kinetics on manifolds. The motion of particles on the manifold is represented using both canonical and non-canonical Hamiltonian formulations. Our schemes apply to both formulations, but the canonical formulation results in a particularly efficient scheme that also conserves particle density and energy exactly. The collisionless update is coupled to a Bhatnagar-Gross-Krook (BGK) collision operator that provides a simplified model for relaxation to local thermodynamic equilibrium. An iterative scheme is constructed to ensure collisional invariants (density, momentum and energy) are preserved numerically. Rotation of the manifold is incorporated by modifying the Hamiltonian while ensuring a canonical formulation. Several test problems, including a kinetic version of the classical Sod shock problem, Kelvin-Helmholtz instability on the surfaces of a sphere and a hyperboloid, with and without rotations, are presented. A prospectus for further development of this approach to simulation of kinetic theory in general relativity is presented.

Discontinuous Galerkin↗
DOE OSTI2026

Collisional stopping power of ions in warm dense matter

Here, a model for the collisional stopping of ions on free electrons in warm dense matter is developed and explored. It is based on plasma kinetic theory, but with modifications to address the warm dense matter regime. Specifically, it uses the Boltzmann-Uehling-Uhlenbeck kinetic equation to incorporate effects of Fermi degeneracy of electrons. The cross section is computed from quantum scattering of electrons and ions occurring via the potential of mean force derived from an average atom model, which incorporates effects of strong Coulomb correlations. Predictions from this model show comparable accuracy to results from time-dependent density functional theory calculations for deuterium near solid density and a temperature of several electronvolts, at a fraction of the computational cost. Further, the model captures the transition of a plasma from the classical limit to the degenerate limit, including qualitative behaviors of solid state theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Perfect spinfluid: A divergence-type approach

We present a new formulation of nondissipative relativistic spin hydrodynamics that incorporates spin degrees of freedom into the divergence-type theory framework. Due to the divergence-type structure, it is straightforward to enforce nonlinear causality and symmetric hyperbolicity of the equations of motion, ensuring local well-posedness of the initial-value problem and stability of the theory. Furthermore, in a specific realization based on spin kinetic theory, we prove that the equations of motion remain nonlinearly causal and symmetric-hyperbolic to all orders in the spin potential, provided a specific thermodynamic constraint is satisfied. Here, this framework can be applied for numerical simulations to study the dynamics of spin-polarized fluids, such as the quark-gluon plasma in heavy-ion collisions.

Chirality↗
DOE OSTI2025

Challenges in integrating dissolved organic matter chemodiversity into kinetic models of soil respiration

The chemodiversity of dissolved organic matter (DOM) in soil has been proposed to influence the microbial metabolism and fate of belowground organic carbon (C). However, integrating DOM chemistry into soil C cycle models to improve predictions of C stocks and fluxes—beyond simply considering DOM pool size—remains a challenge. While recent research suggests that incorporating DOM chemodiversity into models can improve predictions of microbial respiration, there is still a lack of mechanistic understanding describing how DOM chemodiversity affects microbial metabolism and soil respiration. Here, we evaluated whether DOM chemodiversity was a determinant of soil respiration using paired measurements of high-resolution DOM chemistry, obtained from Fourier transform ion cyclotron resonance mass spectrometry (FTICR-MS), and potential soil respiration rates from across the United States (U.S.), all data provided by the Molecular Observation Network. Our objectives were to (1) assess statistical relationships between DOM chemodiversity and microbial respiration, and (2) evaluate the ability of kinetic models to leverage DOM chemistry to explain empirical relationships found in statistical models. Statistical regressions revealed that DOM chemodiversity (alpha diversity) was nonlinearly related to potential soil respiration rates, both independently and through its interactions with DOM and total C concentrations. In soils with relatively high DOM but low total C concentrations, potential soil respiration rates were negatively correlated with DOM alpha diversity, whereas in soils with relatively low DOM and high total C concentrations showed the opposite trend. However, when metabolic transition theory kinetic models were modified to include chemodiversity, their performance was comparable to traditional Monod kinetics approaches, which simulate respiration rates as a function of DOM concentration. The inability to account for nonlinearities in DOM chemodiversity–respiration relationships highlight an opportunity to advance substrate uptake kinetics by establishing causal links between DOM chemodiversity, microbial metabolism trade-offs, and potential interactions under varied environmental conditions.

Bioenergetic model↗
DOE OSTI2025

Asymptotic-preserving dynamical low-rank method for the stiff nonlinear Boltzmann equation

In kinetic theory, numerically solving the full Boltzmann equation is extremely expensive. This is because the Boltzmann collision operator involves a high-dimensional, nonlinear integral that must be evaluated at each spatial grid point and every time step. The challenge becomes even more pronounced in the fluid (strong collisionality) regime, where the collision operator exhibits strong stiffness, causing explicit time integrators to impose severe stability restrictions. In this paper, we propose addressing this problem through a dynamical low-rank (DLR) approximation. The resulting algorithm requires evaluating the Boltzmann collision operator only r 2 times, where r, the rank of the approximation, is much smaller than the number of spatial grid points. We propose a novel DLR integrator, called the XL integrator, which reduces the number of steps compared to the available alternatives (such as the projector splitting or basis update & Galerkin (BUG) integrator). For a class of problems including the Boltzmann collision operator which enjoys a separation property between physical and velocity space, we further propose a specialized version of the XL integrator, called the sXL integrator. This version requires solving only one differential equation to update the low-rank factors. Furthermore, the proposed low-rank schemes are asymptotic-preserving, meaning they can capture the asymptotic fluid limit in the case of strong collisionality. Our numerical experiments demonstrate the efficiency and accuracy of the proposed methods across a wide range of regimes, from non-stiff (kinetic) to stiff (fluid).

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Phase transformation kinetics model for metals

We develop a new model for phase transformation kinetics in metals by generalizing the Levitas–Preston (LP) phase field model of martensite phase transformations (see Levitas and Preston (2002a,b) and Levitas et al. (2003)) to arbitrary pressure. Furthermore, we account for and track: the interface speed of the pressure-driven phase transformation, properties of critical nuclei, as well as nucleation at grain sites and on dislocations and homogeneous nucleation. The volume fraction evolution of each phase is described by employing KJMA (Kolmogorov, 1937; Johnson and Mehl, 1939; Avrami, 1939, 1940, 1941) kinetic theory. We then test our new model for iron under ramp loading conditions and compare our predictions for the α → ϵ iron phase transition to experimental data of Smith et al. (2013). In conclusion, more than one combination of material and model parameters (such as dislocation density and interface speed) led to good agreement of our simulations to the experimental data, thus highlighting the importance of having accurate microstructure data for the sample under consideration.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Stochastic fluctuations and the relaxation time in transient relativistic fluids

We argue that the ratio between the shear viscosity and the shear relaxation time, η/τ π , should be defined as a thermodynamic quantity obtained from the equal-time symmetric correlator of the shear-stress tensor. In kinetic theory, we show that this ratio does not depend on the type of interaction. Similarly, an exact expression for this ratio is obtained for holographic gauge theories. We also determine how stochastic fluctuations change η/τ π in transient relativistic hydrodynamics and show that thermal fluctuations do not spoil causality and stability.

Denicol, Gabriel S. [Universidade Federal Fluminen↗
DOE OSTI2026

Machine-learned quantum molecular dynamics calculations of warm dense equation of state and ionic transport coefficients of deuterated water

White dwarf models require accurate equations of state and ionic transport coefficients in the warm dense matter regime, where kinetic theory models and tabulated equations of state are often inaccurate. In this work, spectral-partitioned density functional theory and machine-learned interatomic potentials are combined to perform large-scale, first-principles quantum molecular dynamics simulations of deuterated water (D 2 O) near the principal Hugoniot. This approach retains Kohn-Sham accuracy while achieving orders-of-magnitude speedup, yielding converged equation of state and transport properties over a broad pressure and temperature range. The results reveal the thermodynamic conditions under which ionic transport models for interdiffusivity and shear viscosity converge and identify those in closest agreement with density functional theory benchmarks at temperatures in the warm dense matter regime. The present framework extends first-principles transport calculations to higher temperatures than previously achieved, and provides an efficient, scalable, and general approach for studying transport properties in complex multicomponent mixtures.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Origin of anomalous magnetotransport in kagome superconductors 𝐴⁢V 3 ⁢Sb 5 (𝐴 = K,Rb,Cs)

Multiple anomalous features in electronic spectra of metals with a kagome lattice structure—van Hove singularities, Dirac points, and flat bands—imply that materials containing this structural motif may lie at a nexus of topological and correlated electron physics. Due to the prospects of such exceptional electronic behavior, the recent discovery of superconductivity coexisting with charge-density wave (CDW) order in the layered kagome metals 𝐴⁢V 3 ⁢Sb 5 (𝐴 = K,Rb,Cs) has attracted considerable attention. Notably, these archetypal kagome metals express unconventional magnetotransport behavior, including an unexpected linear-in-𝐻 diagonal resistivity at low fields, and an even more peculiar, nonmonotonic sign-changing behavior of the Hall resistivity, which has been speculated to arise from a chiral CDW. We argue here that this unusual magnetotransport derives not from such unconventional phenomena, but rather from the unique fermiology of the 𝐴⁢V 3 ⁢Sb 5 materials. Specifically, it is caused by a large, concave hexagonal Fermi surface sheet formed in the close proximity to the van Hove singularities, which is backfolded into a small hexagonal sheet and two large triangular sheets in the CDW state. We introduce and analyze a model of the electronic structure of these Fermi surface sheets that allows for a full analytical treatment within Boltzmann kinetic theory and that enables semi-quantitative fits of our transport data. Specifically, we find that the anomalous magnetotransport behavior is caused by the confluence of strong reduction of the Fermi velocity near the van Hove singularities located near the vertices of the hexagonal sheet and sharp corners in Fermi surface generated by the CDW reconstruction. In conclusion, our analytical approach not only explains the anomalous magnetotransport in the kagome superconductors but also can be extended to a variety of metallic systems hosting singular features in their Fermi surfaces.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗