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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

Theory of a cylindrical probe in a collisionless magnetoplasma

A theory is presented for a cylindrical electrostatic probe in a collisionless plasma in the case where the probe axis is inclined at an angle to a uniform magnetic field. The theory is applicable to electron collection, and under more restrictive conditions, to ion collection. For a probe at space potential, the theory is exact in the limit where probe radius is much less than Debye length. At attracting probe potentials, the theory yields an upper bound and an adiabatic limit for current collection. At repelling probe potentials, it provides a lower bound. The theory is valid if the ratios of probe radius to Debye length and probe radius to mean gyroradius are not simultaneously large enough to produce extrema in the probe sheath potential. The numerical current calculations are based on the approximation that particle orbits are helices near the probe, together with the use of kinetic theory to relate velocity distributions near the probe to those far from it. Probe characteristics are presented for inclination angles from 0 to 90 deg and for probe-radius mean-gyroradius ratios from 0.1 to infinity. For an angle of 0 deg, the end-effect current is calculated separately.

Laframboise, J. G.↗
DOE OSTI2022

Electrochemically-Induced Phase Transformations in Battery Storage Compounds (Final Technical Report)

Compounds of interest for ion storage in advanced batteries frequently exhibit phase transformations as the working ion concentration varies. Under large electrochemical driving forces inherent to practical use, systems are often driven far from equilibrium. This program combines experiments and theory to understand the phase transition behavior of ion insertion compounds when electrochemically driven far from equilibrium. As model systems, we focus on alkaline metal phosphates AMPO4 (A = alkali; M = first row transition metal) of olivine structure, which are both technologically interesting and ideally suited for fundamental study due to the ability to systematically tune transformation strain, and along with it, the phase transformation pathway. Behavior in compositions having large transformation strains (~15 vol%) requiring plasticity for strain accommodation is emphasized. Experimental techniques include operando characterization of structure while simultaneously varying electrokinetic parameters, and high resolution microscopy of nanoscale and interfacial phenomena. Phase-field modeling is used to model the thermodynamics and kinetics of competing transformation pathways, extended to include the effects of plasticity, and integrated with porous electrode kinetic theory to treat multi-particle effects. Success in this project will lead to an ability to design ion storage compounds with predictable transformation pathways, electrochemical kinetics, capacity utilization, and durability. New technologically important compounds may also be discovered.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Electron kinetics in a positive column of AC discharges in a dynamic regime

Abstract We have performed hybrid kinetic-fluid simulations of a positive column in alternating current (AC) argon discharges over a range of driving frequencies f and gas pressure p for the conditions when the spatial nonlocality of the electron energy distribution function (EEDF) is substantial. Our simulations confirmed that the most efficient conditions of plasma maintenance are observed in the dynamic regime when time modulations of mean electron energy (temperature) are substantial. The minimal values of the root mean square electric field and the electron temperature have been observed at f/p values of about 3 kHz Torr −1 in a tube of radius R = 1 cm. The ionization rate and plasma density reached maximal values under these conditions. The numerical solution of a kinetic equation allowed accounting for the kinetic effects associated with spatial and temporal nonlocality of the EEDF. Using the kinetic energy of electrons as an independent variable, we solved an anisotropic tensor diffusion equation in phase space. We clarified the role of different flux components during electron diffusion in phase space over surfaces of constant total energy. We have shown that the kinetic theory uncovers a more exciting and rich physics than the classical ambipolar diffusion (Schottky) model. Non-monotonic radial distributions of excitation rates, metastable densities, and plasma density have been observed in our simulations at pR > 6 Torr cm. The predicted off-axis plasma density peak in the dynamic regime has never been observed in experiments so far. We hope our results stimulate further experimental studies of the AC positive column. The kinetic analysis could help uncover new physics even for such a well-known plasma object as a positive column in noble gases.

Physics↗
DOE OSTI2020

A general metriplectic framework with application to dissipative extended magnetohydrodynamics

General equations for conservative yet dissipative (entropy producing) extended magnetohydrodynamics are derived from two-fluid theory. Keeping all terms generates unusual cross-effects, such as thermophoresis and a current viscosity that mixes with the usual velocity viscosity. While the Poisson bracket of the ideal version of this model has already been discovered, we determine its metriplectic counterpart that describes the dissipation. This is done using a new and general thermodynamic point of view to derive dissipative brackets, a means of derivation that is natural for understanding and creating dissipative dynamics without appealing to underlying kinetic theory orderings. Finally, the formalism is used to study dissipation in the Lagrangian variable picture where, in the context of extended magnetohydrodynamics, non-local dissipative brackets naturally emerge.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1993

An upwind, kinetic flux-vector splitting method for flows in chemical and thermal non-equilibrium

We have developed new upwind kinetic difference schemes for flows with non-equilibrium thermodynamics and chemistry. These schemes are derived from the Boltzmann equation with the resulting Euler schemes developed as moments of the discretized Boltzmann scheme with a locally Maxwellian velocity distribution. Splitting the velocity distribution at the Boltzmann level is seen to result in a flux-split Euler scheme and is called Kinetic Flux Vector Splitting (KFVS). Extensions to flows with finite-rate chemistry and vibrational relaxation is accomplished utilizing nonequilibrium kinetic theory. Computational examples are presented comparing KFVS with the schemes of Van Leer and Roe for a quasi-one-dimensional flow through a supersonic diffuser, inviscid flow through two-dimensional inlet, and viscous flow over a cone at zero angle-of-attack. Calculations are also shown for the transonic flow over a bump in a channel and the transonic flow over an NACA 0012 airfoil. The results show that even though the KFVS scheme is a Riemann solver at the kinetic level, its behavior at the Euler level is more similar to the existing flux-vector splitting algorithms than to the flux-difference splitting scheme of Roe.

Eppard, W. M.↗
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↗
NASA NTRS1987

Cellular automatons applied to gas dynamic problems

This paper compares the results of a relatively new computational fluid dynamics method, cellular automatons, with experimental data and analytical results. This technique has been shown to qualitatively predict fluidlike behavior; however, there have been few published comparisons with experiment or other theories. Comparisons are made for a one-dimensional supersonic piston problem, Stokes first problem, and the flow past a normal flat plate. These comparisons are used to assess the ability of the method to accurately model fluid dynamic behavior and to point out its limitations. Reasonable results were obtained for all three test cases, but the fundamental limitations of cellular automatons are numerous. It may be misleading, at this time, to say that cellular automatons are a computationally efficient technique. Other methods, based on continuum or kinetic theory, would also be very efficient if as little of the physics were included.

Long, Lyle N.↗
DOE OSTI2023

Expanding-box Quasilinear Model of the Solar Wind

The expanding-box model of the solar wind has been adopted in the literature within the context of magnetohydrodynamics, hybrid, and full particle-in-cell simulations to investigate the dynamic evolution of the solar wind. The present paper extends such a method to the framework of self-consistent quasilinear kinetic theory. It is shown that the expanding-box quasilinear methodology is largely equivalent to the inhomogeneous steady-state quasilinear model discussed earlier in the literature, but a distinction regarding the description of wave dynamics between the two approaches is also found. The expanding-box quasilinear formalism is further extended to include the effects of a spiraling solar-wind magnetic field as well as collisional age effects. The present finding shows that the expanding-box quasilinear approach and the steady-state global-kinetic models may be employed interchangeably in order to address other more complex problems associated with the solar-wind dynamics.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1979

A fluid description of plasma double-layers

The space-charge double-layer that forms between two plasmas with different densities and thermal energies was investigated using three progressively realistic models which are treated by fluid theory, and take into account four species of particles: electrons and ions reflected by the double-layer, and electrons and ions transmitted through it. The two plasmas are assumed to be cold, and the self-consistent potential, electric field and space-charge distributions within the double-layer are determined. The effects of thermal velocities are taken into account for the reflected particles, and the modifications to the cold plasma solutions are established. Further modifications due to thermal velocities of the transmitted particles are examined. The applicability of a one dimensional fluid description, rather than plasma kinetic theory, is discussed. Theoretical predictions are compared with double layer potentials and lengths deduced from laboratory and space plasma experiments.

Levine, J. S.↗
DOE OSTI2023

Recent progress on in-medium properties of heavy mesons from finite-temperature EFTs

Mesons with heavy flavor content are an exceptional probe of the hot QCD medium produced in heavy-ion collisions. In the past few years, significant progress has been made toward describing the modification of the properties of heavy mesons in the hadronic phase at finite temperature. Ground-state and excited-state thermal spectral properties can be computed within a self-consistent many-body approach that employs appropriate hadron-hadron effective interactions, providing a unique opportunity to confront hadronic Effective Field Theory predictions with recent and forthcoming lattice QCD simulations and experimental data. In this article, we revisit the application of the imaginary-time formalism to extend the calculation of unitarized scattering amplitudes from the vacuum to finite temperature. These methods allow us to obtain the ground-state thermal spectral functions. The thermal properties of the excited states that are dynamically generated within the molecular picture are also directly accessible. We present here the results of this approach for the open-charm and open-bottom sectors. We also analyze how the heavy-flavor transport properties, which are strongly correlated to experimental observables in heavy-ion collisions, are modified in hot matter. In particular, transport coefficients can be computed using an off-shell kinetic theory that is fully consistent with the effective theory describing the scattering processes. The results of this procedure for both charm and bottom transport coefficients are briefly discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
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 OSTI2021

Non-Contact Mass Density and Thermal Conductivity Measurements of Organic Thin Films Using Frequency–Domain Thermoreflectance

Thin-film organic materials are broadly used to study amorphous stabilization of active pharmaceuticals, control explosive detonation phenomena, and introduce insulation in novel thermal barriers. Their synthesis, however, introduces defects and thickness variations that warrant careful characterization of local thermophysical properties such as thermal conductivity and mass density. In this study, wide bandwidth (200 Hz to 20 MHz) frequency–domain thermoreflectance (FDTR) is demonstrated to simultaneously extract the thermal conductivity and mass density of 1 μm physical vapor-deposited indomethacin films on Si and SiO 2 substrates, as well as 10 and 100 μm films on Si. By assuming a bulk specific heat capacity, mass densities are determined with FDTR measurements of volumetric heat capacity and are in good agreement with the literature, as well as models based upon a dependence on porosity and the kinetic theory for phonons. Lastly, it is found that for broad-band FDTR measurements, insulating substrates provide improved fidelity for the extraction of thermal conductivity and volumetric heat capacity in organic thin films. Overall, this work demonstrates the potential for FDTR as a non-contact method to determine microscale mass density variations across the surface and thickness of organic thin films.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Spin relaxation rate for heavy quarks in weakly coupled QCD plasma

We compute the relaxation rate of the spin density of heavy quarks in a perturbative QCD plasma to leading-log order in the coupling constant g. The spin relaxation rate Γ s in spin hydrodynamics is shown to be Γ s ~ g 4 log(1/g)T(T/M) 2 in the heavy-quark limit T/M $\ll$ 1, which is smaller than the relaxation rate of other non- hydrodynamic modes by additional powers of T/M. We demonstrate three different methods to evaluate the spin relaxation rate: 1) the Green-Kubo formula in the spin hydrodynamic regime, 2) the spin density correlation function in the strict hydrodynamic limit, and 3) quantum kinetic theory of the spin distribution function in momentum space. We highlight the interesting differences between these methods, while they are ultimately connected to each other by the underlying Ward-Takahashi identity for the non-conserved spin density.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Thermal diffusivity recovery and defect annealing kinetics of self-ion implanted tungsten probed by insitu transient grating spectroscopy

Tungsten is a promising candidate material for plasma-facing armour components in future fusion reactors. A key concern is irradiation-induced degradation of its normally excellent thermal transport properties. In this comprehensive study, thermal diffusivity degradation in ion-implanted tungsten and its evolution from room temperature (RT) to 1073 K is considered. Five samples were exposed to 20 MeV self-ions at RT to achieve damage levels ranging from 3.2 x 10 -4 to 3.2 displacements per atom (dpa). Transient grating spectroscopy with in-situ heating was then used to study thermal diffusivity evolution as a function of temperature. Using a kinetic theory model, an equivalent point defect density is estimated from the measured thermal diffusivity. The results showed a prominent recovery of thermal diffusivity between 450 K and 650 K, which coincides with the onset of mono-vacancy mobility. After 1073 K annealing samples with initial damage of 3.2 x 10 -3 dpa or less recover close to the pristine value of thermal diffusivity. For doses of 3.2 x 10 -2 dpa or higher, on the other hand, a residual reduction in thermal diffusivity remains even after 1073 K annealing. Transmission electron microscopy reveals that this is associated with extended, irradiation-induced dislocation structures that are retained after annealing. A sensitivity analysis shows that thermal diffusivity provides an efficient tool for the assessing total defect content in tungsten up to 1000 K.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Direct numerical simulation of a turbulent methane/air flame impinged by a sub-breakdown electric field

Here, the effects of an external electric field on a turbulent methane/air diffusion flame are analyzed in this work using direct numerical simulations. The analyzed configuration consists of a temporally evolving mixing layer of air and a mixture of methane and nitrogen at 1 atm that is impinged by an electric field in the direction normal to the mean mixing plane. The combustion and chemi-ionization reactions involved in the flow are modeled using finite-rate chemistry and a reduced reaction scheme consisting of 26 species and 134 reactions. The mass diffusion and ion-wind effects are modeled using a detailed description of the diffusion coefficients and electric mobilities based on kinetic theory. The presented calculations show that the turbulence generated within the mixing layer is mostly unaffected by the applied electric field for the configuration under exam. In fact, the electric body force is developed away from the mixing region where the flow is uniform. Conversely, the turbulence is able to introduce very high intermittency in the electrically charged species concentration and, consequently, in the electric body force. Such intermittency will constitute a challenge in the future formulation of combustion models that take into account ion-wind effects.

09 BIOMASS FUELS↗
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 OSTI2021

Probing Early-Time Dynamics and Quark-Gluon Plasma Transport Properties with Photons and Hadrons

In this work, we use a hybrid model which relies on QCD effective kinetic theory - KøMPøST - to dynamically bridge the gap between IP-Glasma initial states and viscous hydrodynamics, for the theoretical interpretation of heavy-ion collision results. The hydrodynamic phase is then followed by dynamical freeze-out handled by UrQMD. It is observed that the new pre-equilibrium/pre-hydro phase does influence the extraction of transport coefficients. We also show how the early dynamics is reflected in photon observables.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗