Search NASA⌕ Search

SEARCH · Search NASA

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Quantum kinetic modeling of KEEN waves in a warm-dense regime

We report the first fully kinetic, quantum study of kinetic electrostatic electron nonlinear (KEEN) waves, showing that quantum diffraction systematically erodes the classical trapping mechanism, narrows harmonic locking to the fundamental, and hastens post-drive decay. Electrons are evolved with a second-order Strang-split 1D1V Wigner–Poisson solver that couples conservative semi-Lagrangian WENO advection to an analytic Fourier space update for the non-local Wigner term, while ions remain classical. We focus on collisionless dynamics in a weakly coupled regime, providing a controlled baseline before collisional extensions. Short, frequency-tuned ponderomotive pulses drive KEEN formation in a uniform Maxwellian plasma; as the dimensionless quantum parameter H rises from the classical limit to values relevant to warm-dense matter, doped semiconductors, and 2D electron systems, the drive threshold increases, higher harmonics are damped, trapped electron vortices diffuse, and the subplasma electrostatic energy relaxes to a lower stationary level, as confirmed by continuous wavelet analysis. These microscopic changes carry macroscopic weight. Ignition-scale capsules now compress matter to regimes where the electron de Broglie wavelength rivals the Debye length, making classical kinetic descriptions insufficient. By extending KEEN physics into this quantum domain, our results offer a potential diagnostic of non-equilibrium electron dynamics for next-generation inertial-confinement designs and high-energy-density platforms, indicating that predictive fusion modeling may benefit from the integration of kinetic fidelity with quantum effects.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Molecular dynamics simulation of hydrodynamic transport coefficients in plasmas

Molecular dynamics simulations are used to calculate transport coefficients in a two-component plasma interacting through a repulsive Coulomb potential. The thermal conductivity, electrical conductivity, electrothermal coefficient, thermoelectric coefficient, and shear viscosity are computed using the Green–Kubo formalism over a broad range of Coulomb coupling strength, 0.01 ≤ Γ ≤ 140. Emphasis is placed on testing standard results of the Chapman–Enskog solution in the weakly coupled regime (Γ ≪ 1) using these first-principles simulations. As expected, the results show good agreement for Γ ≲ 0.1. However, this agreement is only possible if careful attention is paid to the definitions of linear constitutive relations in each of the theoretical models, a point that is often overlooked. For example, the standard Green–Kubo expression for thermal conductivity is a linear combination of thermal conductivity, electrothermal, and thermoelectric coefficients computed in the Chapman–Enskog formalism. Meaningful results for electrical conductivity are obtained over the full range of coupling strengths explored, but it is shown that potential and virial components of the other transport coefficients diverge in the strongly coupled regime (Γ ≫ 1). In this regime, only the kinetic components of the transport coefficients are meaningful for a classical plasma.

Electrical conductivity↗
DOE OSTI2022

Cross-beam energy transfer saturation: ion heating and pump depletion

Cross-beam energy transfer (CBET) was measured in two regimes where the energy transfer saturation mechanism was determined by the plasma and laser beam conditions. Here, linear kinetic CBET theory agrees well with the measured energy transfer in all experimental configurations and at all probe beam intensities when accounting for pump depletion and the plasma conditions measured using Thomson-scattering. Simultaneous CBET and Thomson-scattering measurements enabled uncertainties in the plasma conditions to be isolated from CBET theory, which allowed the saturation mechanisms to be identified. In the perpendicular-beam configuration the saturation mode was through ion heating, which resulted from ion trapping in the driven waves and subsequent ion-ion collisional heating. In the co-propagating beam configuration there was minimal ion heating and the saturation mode was through pump depletion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A kinetic-based regularization method for data science applications

We propose a physics-based regularization technique for function learning, inspired by statistical mechanics. By drawing an analogy between optimizing the parameters of an interpolator and minimizing the energy of a system, we introduce corrections that impose constraints on the lower-order moments of the data distribution. This minimizes the discrepancy between the discrete and continuum representations of the data, in turn allowing to access more favorable energy landscapes, thus improving the accuracy of the interpolator. Our approach improves performance in both interpolation and regression tasks, even in high-dimensional spaces. Unlike traditional methods, it does not require empirical parameter tuning, making it particularly effective for handling noisy data. We also show that thanks to its local nature, the method offers computational and memory efficiency advantages over Radial Basis Function interpolators, especially for large datasets.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Emergent viscous hydrodynamics from a single quantum particle

We investigate an explicit example of how spatial decoherence can lead to hydrodynamic behavior in the late-time, long-wavelength regime of open quantum systems. We focus on the case of a single nonrelativistic quantum particle linearly coupled to a thermal bath of noninteracting harmonic oscillators at temperature T , a la Caldeira and Leggett. Taking advantage of decoherence in the position representation, we expand the reduced density matrix in powers of the off-diagonal spatial components, so that high-order terms are suppressed at late times. Truncating the resulting power series at second order leads to a set of dissipative transient hydrodynamic equations similar to the nonrelativistic limit of equations widely used in simulations of the quark-gluon plasma formed in ultrarelativistic heavy-ion collisions. Transport coefficients are directly determined by the damping constant γ , which quantifies the influence of the environment. The asymptotic limit of our hydrodynamic equations reduces to the celebrated Navier-Stokes equations for a compressible fluid in the presence of a drag force. Furthermore, our results shed new light on the onset of hydrodynamic behavior in open quantum systems where a system with few degrees of freedom is coupled to a large thermal environment.

Hydrodynamics↗
DOE OSTI2025

Relaxation time approximation for a multispecies relativistic gas

We generalize a recent prescription for the relaxation time approximation for the relativistic Boltzmann equation for systems with multiple particle species at finite temperature. This is performed by adding counter-terms to the traditional Anderson-Witting ansatz for each particle species. Our approach allows for the use of momentum-dependent relaxation times and the obedience of local conservation laws regardless of the definition of the local equilibrium state. As an application, we derive the first order Chapman-Enskog corrections to the equilibrium distribution and display results for the hadron-resonance gas. We also demonstrate that our collision term ansatz obeys the second law of thermodynamics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Neutrino flavor instabilities in neutron star mergers with moment transport: Slow, fast, and collisional modes

Determining where, when, and how neutrino flavor oscillations must be included in large-scale simulations of hot and dense astrophysical environments is an enduring challenge that must be tackled to obtain accurate predictions. Here, using an angular moment-based linear stability analysis framework, we examine the different kinds of flavor instabilities that can take place in the context of the postprocessing of a neutron star merger simulation, with a particular focus on the collisional flavor instability and a careful assessment of several commonly used approximations. First, neglecting anisotropies of the neutrino field, we investigate the extent to which commonly used monoenergetic growth rates reproduce the results obtained from a full multienergy treatment. Contrary to the large discrepancies found in core-collapse supernova environments, we propose a simple combination of energy-averaged estimates that reproduces the multienergy growth rates in our representative simulation snapshot. We then quantify the impact of additional physical effects, including nuclear many-body corrections, scattering opacities, and the inclusion of the vacuum term in the neutrino Hamiltonian. Finally, we include the neutrino distribution anisotropies, which allows us to explore, for the first time in a multienergy setting, the interplay between collisional, fast, and slow modes in a moment-based neutron star merger simulation. We find that, despite a dominance of the fast instability in most of the simulation volume, certain regions exhibit only a collisional instability, while others, especially at large distances, exhibit a slow instability that is largely underestimated if anisotropic effects are neglected.

neutrino oscillations↗
DOE OSTI2021

Conformal maps of viscous electron flow in the Gurzhi crossover

Herein we investigate the impact of geometric constriction on the viscous flow of electron liquid through quantum point contacts. We provide analysis on the electric potential distribution given the setup of a slit configuration and use the method of conformal mapping to obtain analytical results. The potential profile can be tested and contrasted experimentally with the scanning tunneling potentiometry technique. We discuss intricate physics that underlies the Gurzhi effect, i.e., the enhancement of conductivity in the viscous flow, and compare results for different boundary conditions. In addition, we calculate the temperature dependence of the momentum relaxation time as a result of impurity assisted quasiballistic interference effects and discuss various correlational corrections that lead to the violation of Matthiessen's rule in the hydrodynamic regime. We caution that spatially inhomogeneous profiles of current in the Gurzhi crossover between Ohmic and Stokes flows might also appear in the nonhydrodynamic limit where nonlocality plays an important role. This conclusion is corroborated by calculation of dispersive conductivity in the weakly impure limit.

2-dimensional systems↗
DOE OSTI2021

Prehydrodynamic evolution and its signatures in final-state heavy-ion observables

We investigate the effects of prehydrodynamic evolution on final-state observables in heavy-ion collisions using state-of-the art event simulations coupled to different prehydrodynamic scenarios, which include the recently developed effective kinetic transport theory evolution model KøMPøST. Differential flow observables are found to be mostly insensitive to the details of prehydrodynamic evolution. The main effect we observe is in the p T spectra, particularly the mean transverse momentum. However, at least part of this effect is a consequence of the underlying conformal invariance assumption currently present in such approaches, which is known to be violated in the temperature regime probed in heavy-ion collisions. This assumption of early time conformal invariance leads to an artificially large out-of-equilibrium bulk pressure when switching from (conformal) prehydrodynamic evolution to hydrodynamics (using the nonconformal QCD equation of state), which in turn increases the transverse momentum. Our study indicates that a consistent treatment of prehydrodynamic evolution in heavy-ion collisions requires the use of nonconformal models of early time dynamics.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Modified equilibrium distributions for Cooper-Frye particlization

We introduce a positive-definite single-particle distribution that is suitable for describing the transition from a macroscopic hydrodynamic to a microscopic kinetic description during the late stages of heavy-ion collisions in the presence of moderately large viscous corrections. The modified equilibrium distribution function can be constructed with hydrodynamic input from either relativistic viscous fluid dynamics or anisotropic fluid dynamics. We test the modified equilibrium distribution's hydrodynamic output for a stationary hadron resonance gas subject to either shear stress, bulk pressure, or baryon diffusion current at a given freeze-out temperature and baryon chemical potential. While it does not reproduce all components of the net baryon current and energy-momentum tensor exactly, it significantly improves upon the customary linearized approximations for the nonequilibrium correction δf n which typically lead to unphysical negative distribution functions at large particle momenta. A comparison of particle spectra and p T –differential elliptic flow coefficients from the Cooper–Frye formula computed with the modified equilibrium distribution and with linearized δf n corrections is presented, for two different (2 + 1)–dimensional hypersurfaces corresponding to central and noncentral Pb + Pb collisions at the Large Hadron Collider (LHC).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Building a testable shear viscosity across the QCD phase diagram

Current experiments at the Relativistic Heavy Ion Collider (RHIC) are probing finite baryon densities where the shear viscosity to enthalpy ratio ηT/w of the Quark Gluon Plasma remains unknown. Here we use the Hadron Resonance Gas (HRG) model with the most up-to-date hadron list to calculate ηT/w at low temperatures and at finite baryon densities ρ B . We then match ηT/w to a QCD-based shear viscosity calculation within the deconfined phase to create a table across {T,μ B } for different cross-over and critical point scenarios at a specified location. We find that these new ηT/w(T,μ B ) values would require initial conditions at significantly larger ρB, compared to ideal hydrodynamic trajectories, in order to reach the same freeze-out point.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Nonequilibrium effects and transverse spherocity in ultrarelativistic proton-nucleus collisions

We investigate the effects of nonequilibrium dynamics in small colliding systems by comparing a nonequilibrium transport approach, Parton-Hadron-String Dynamics (PHSD), with a (2+1)-dimensional viscous hydrodynamic model, VISHNew. Focusing on p + Pb collisions at energies available at the Large Hadron Collider, we extract the initial conditions for the hydrodynamic model from PHSD, in order to reduce the impact of the early out-of-equilibrium dynamics and focus on the traces of nonequilibrium in the ensuing medium evolution. We find that in the transport approach quantities like energy density and bulk viscous pressure are highly inhomogeneous on the transverse plane during the whole evolution, whereas the hydrodynamic simulations dissolve more efficiently the initial spatial irregularities, still keeping a high degree of inhomogeneity due to the smaller size and lifetime of the medium produced in p + Pb collisions with respect to heavy-ion reactions. As a first step that will help to identify the impact of these nonequilibrium effects on final observables in proton-nucleus collisions, we perform an analysis of the transverse spherocity, an event-shape observable able to distinguish between jetty and isotropic configurations of transverse momenta. We found that the spherocity distribution in PHSD is slightly shifted towards the isotropic limit with respect to the hydrodynamic result. Even though this dissimilarity is partially due to the difference in the final charged particle production, it mainly comes from the different description within the two frameworks of the medium produced in small colliding systems. This finding supports the idea that multidifferential measurements, such as those based on event categorization according to multiplicity and spherocity, are useful to study final-state observables in ultrarelativistic proton-nucleus collisions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Phase transition amplification of proton number fluctuations in nuclear collisions from a transport model approach

The time evolution of particle number fluctuations in nuclear collisions at intermediate energies (E lab = 1.23–10⁢A GeV) is studied by means of the UrQMD-3.5 transport model. The transport description incorporates baryonic interactions through a density-dependent potential. This allows for an implementation of a first-order phase transition including a mechanically unstable region at large baryon density. The scaled variance of the baryon and proton number distributions is calculated in the central cubic spatial volume of the collisions at different times. A significant enhancement of fluctuations associated with the unstable region is observed. This enhancement persists to late times, reflecting a memory effect for the fluctuations. The presence of the phase transition has a much smaller influence on the observable event-by-event fluctuations of protons in momentum space.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Neutrino fast flavor oscillations with moments: Linear stability analysis and application to neutron star mergers

Providing an accurate modeling of neutrino physics in dense astrophysical environments such as binary neutron star mergers presents a challenge for hydrodynamic simulations. Nevertheless, understanding how flavor transformation can occur and affect the dynamics, the mass ejection, and the nucleosynthesis will need to be achieved in the future. Computationally expensive, large-scale simulations frequently evolve the first classical angular moments of the neutrino distributions. By promoting these quantities to matrices in flavor space, we develop a linear stability analysis of fast flavor oscillations using only the first two “quantum” moments, which notably requires generalizing the classical closure relations that appropriately truncate the hierarchy of moment equations in order to treat quantum flavor coherence. After showing the efficiency of this method on a well-understood test situation, we perform a systematic search of the occurrence of fast flavor instabilities in a neutron star merger simulation. Here, we discuss the successes and shortcomings of moment linear stability analysis, as this framework provides a time-efficient way to design and study better closure prescriptions in the future.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Bhatnagar-Gross-Krook subgrid model for neutrino quantum kinetics

Here, we present a new subgrid model for neutrino quantum kinetics, which is primarily designed to incorporate effects of collective neutrino oscillations into neutrino-radiation-hydrodynamic simulations for core-collapse supernovae and mergers of compact objects. We approximate the neutrino oscillation term in a quantum kinetic equation using the Bhatnagar-Gross-Krook (BGK) relaxation-time prescription, and the transport equation is directly applicable for classical neutrino transport schemes. The BGK model is motivated by recent theoretical indications that nonlinear phases of collective neutrino oscillations settle into quasisteady structures. We explicitly provide basic equations of the BGK subgrid model for both multiangle and moment-based neutrino transport to facilitate the implementation of the subgrid model in the existing neutrino transport schemes. We also show the capability of our BGK subgrid model by comparing it to fully quantum kinetic simulations for fast neutrino-flavor conversion. We find that the overall properties can be well reproduced in the subgrid model; the error of angular-averaged survival probability of neutrinos is within ~20%. By identifying the source of error, we also discuss perspectives to improve the accuracy of the subgrid model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Generation of arbitrary probability distributions from chaotic dynamics

A framework is developed to derive nonlinear dynamical systems that chaotically generate arbitrary multivariate probability distributions with smooth, deterministic trajectories. Herein, the ideas are used to extend the Nosé-Hoover thermostat methodology to three-dimensional cases where the momentum distribution is nonisotropic and/or non-Gaussian. Toy models that generate several well-known distributions in physics are given as pedagogical examples.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗