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 145 records · Page 8

Dynamics of rapidly spinning blob-filaments: Fluid theory with a parallel kinetic extension

Blob-filaments (or simply “blobs”) are coherent structures formed by turbulence and sustained by nonlinear processes in the edge and scrape-off layer (SOL) of tokamaks and other magnetically confined plasmas. The dynamics of these blob-filaments, in particular, their radial motion, can influence the scrape-off layer width and plasma interactions with both the divertor target and with the main chamber walls. Motivated by recent results from the XGC1 gyrokinetic simulation code reported on elsewhere [J. Cheng et al., Nucl. Fusion 63, 086015 (2023)], a theory of rapidly spinning blob-filaments has been developed for this work. The theory treats blob-filaments in the closed flux surface region or the region that is disconnected from sheaths in the SOL. It extends previous work by treating blob spin, arising from partially or fully adiabatic electrons, as the leading-order effect and retaining inertial (ion charge polarization) physics in next order. Spin helps to maintain blob coherency and affects the blob's propagation speed. Dipole charge polarization, treated perturbatively, gives rise to blob-filaments with relatively slow radial velocity, comparable to that observed in the simulations. The theory also treats the interaction of rapidly spinning blob-filaments with a zonal flow layer. It is shown analytically that the flow layer can act like a transport barrier for these structures. Finally, parallel electron kinetic effects are incorporated into the theory. Various asymptotic parameter regimes are discussed, and asymptotic expressions for the radial and poloidal motion of the blob-filaments are obtained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Comparing ab initio and quantum-kinetic approaches to electron transport in warm dense matter

Accurate knowledge of the electronic transport properties of warm dense matter is one of the main concerns of research in high-energy-density physics. Three modern approaches with vastly different levels of fidelity are reviewed and compared: the Kubo–Greenwood (KG) approach based on density-functional-theory molecular dynamics simulations (QMD), quantum kinetic theory based on average-atom models, and time-dependent density functional theory. Throughout, emphasis is placed on the connection between static properties of the electrons (e.g., density of states) and transport properties. Overall, it is found that whenever the conduction electrons can be modeled as being nearly free, fair to excellent agreement is found between QMD and kinetic theory approaches. Such a circumstance is required for modeling warm dense matter as a plasma of ions and free electrons, which is assumed in most kinetic theory approaches. The sensitivity of transport properties to the electronic structure is further highlighted by comparing different exchange–correlation approximations in QMD and KG calculations. It is found that the inclusion of exact exchange via thermal hybrid functionals can make a pronounced impact on electrical and thermal conduction in warm dense matter. We also investigate dynamic screening physics via kinetic theory and time-dependent density functional theory calculations of the mean free path of an electron in a hot dense plasma. In sum, we identify three axes along which to make progress in predicting electron transport in warm dense matter.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electron-scale temperature gradients in kinetic equilibrium: MMS observations and Vlasov–Maxwell solutions

For an ideal gas, the notion of thermodynamic equilibrium requires that the system’s temperature be spatially uniform. We demonstrate that for a collisionless plasma, kinetic equilibrium can be achieved even in the presence of a temperature gradient. Motivated by recent MMS magnetopause observations of thin current layers exhibiting electron-scale temperature gradients, here we present an exact solution to the Vlasov-Maxwell system which accounts for the observed electron temperature gradient and reproduces much of the plasma’s spatial variation along one dimension. Furthermore, the current layer in the Vlasov-Maxwell model is self-consistently supported by a prominent crescent shaped electron velocity distribution that naturally arises due to the temperature transition, which agrees with the electron distribution structures measured by MMS during two asymmetric current layers encountered at the magnetopause.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quantum turbulence, superfluidity, non-Markovian dynamics, and wave function thermalization

While quantum turbulence has been addressed both experimentally (predominantly for superfluid He 4 and He 3 ) and theoretically, the dynamics of various ensembles of quantized vortices has been followed in time only until the vortices have decayed into phonons. How this “thermalization” is achieved is still an unaddressed and thus an unelucidated question. The unitary Fermi gas (UFG) is a unique quantum system, which has no classical counterpart and is of relevance to neutron stars, cold atoms, condensed-matter and nuclear many-body systems. The non-Markovian evolution of an isolated UFG is put in evidence and its entire nonequilibrium evolution can be studied theoretically within a unified theoretical framework. The initial lattice of quantum vortices and antivortices evolves through a couple of vortex tangles and excitation of Kelvin waves, where vortices cross and reconnect, until very slowly thermalization sets in. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Collision integrals within the Chapman–Enskog theory for a generalized Lennard-Jones potential

Here we report the values of the collision integrals, needed for the calculation of the macroscopic transport properties such as viscosity (η) and diffusion coefficient (D) of gases within the Chapman–Enskog kinetic gas theory, for a generalized Lennard-Jones potential (gLJ), a more general potential with an adjustable long range 1/r dependence that can describe a wide range of intermolecular interactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Theory, observations, and simulations of kinetic entropy in a magnetotail electron diffusion region

In this work, we examine velocity-space kinetic entropy, a spatially local measure of entropy for systems out of thermal equilibrium, during an encounter of an electron diffusion region at a magnetic reconnection site in Earth’s magnetotail by the Magnetospheric Multiscale (MMS) mission. We start by generalizing the theory of kinetic entropy to the case of non-uniform velocity space grids and transforming the equations into spheri- cal energy coordinates useful to experimental plasma detectors. The theory is then applied to MMS data and compared to particle-in-cell simulations of reconnection. We demonstrate that the entropy-based non-Maxwellianity measure from the MMS data is of sufficiently high precision to reliably identify non-Maxwellian distributions and therefore the measurements when kinetic effects are most significant. By comparing two different non-Maxwellian measures, we show that total entropy density suffers from “information loss” because it lacks a dependence on the velocity space grid, and so has lost information about how well a distribution function is resolved. Local velocity-space kinetic entropy density recovers this information. We quantify information loss and argue that the considerations needed to minimize it are crucial for instruments designed to measure distribution functions in situ.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Efficient Nonthermal Ion and Electron Acceleration Enabled by the Flux-Rope Kink Instability in 3D Nonrelativistic Magnetic Reconnection

The relaxation of field-line tension during magnetic reconnection gives rise to a universal Fermi acceleration process involving the curvature drift of particles. However, the efficiency of this mechanism is limited by the trapping of energetic particles within flux ropes. Using 3D fully kinetic simulations, we demonstrate that the flux-rope kink instability leads to strong field-line chaos in weak-guide-field regimes where the Fermi mechanism is most efficient, thus allowing particles to transport out of flux ropes and undergo further acceleration. As a consequence, both ions and electrons develop clear power-law energy spectra that contain a significant fraction of the released energy. The low-energy bounds are determined by the injection physics, while the high-energy cutoffs are limited only by the system size. Furthermore, these results have strong relevance to observations of nonthermal particle acceleration in space and astrophysics.

79 ASTRONOMY AND ASTROPHYSICS↗

Kinetic entropy-based measures of distribution function non-Maxwellianity: theory and simulations

We investigate kinetic entropy-based measures of the non-Maxwellianity of distribution functions in plasmas, i.e. entropy-based measures of the departure of a local distribution function from an associated Maxwellian distribution function with the same density, bulk flow and temperature as the local distribution. First, we consider a form previously employed by Kaufmann & Paterson (J. Geophys. Res., vol. 114, 2009, A00D04), assessing its properties and deriving equivalent forms. To provide a quantitative understanding of it, we derive analytical expressions for three common non-Maxwellian plasma distribution functions. We show that there are undesirable features of this non-Maxwellianity measure including that it can diverge in various physical limits and elucidate the reason for the divergence. We then introduce a new kinetic entropy-based non-Maxwellianity measure based on the velocity-space kinetic entropy density, which has a meaningful physical interpretation and does not diverge. We use collisionless particle-in-cell simulations of two-dimensional anti-parallel magnetic reconnection to assess the kinetic entropy-based non-Maxwellianity measures. We show that regions of non-zero non-Maxwellianity are linked to kinetic processes occurring during magnetic reconnection. We also show the simulated non-Maxwellianity agrees reasonably well with predictions for distributions resembling those calculated analytically. These results can be important for applications, as non-Maxwellianity can be used to identify regions of kinetic-scale physics or increased dissipation in plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Quasicrystal stability and nucleation kinetics from density functional theory

The aperiodic order of quasicrystals bridges the amorphous and crystalline regime, so it has remained unclear whether quasicrystals are metastable or stable phases of matter. Density functional theory is often used to evaluate thermodynamic stability, but quasicrystals are long-range aperiodic and their energies cannot be calculated using conventional ab initio methods. Here, in this work, we perform first-principles calculations on quasicrystal nanoparticles of increasing size, from which we can directly extrapolate their bulk and surface energies. Using this technique, we determine with high confidence that the icosahedral quasicrystals ScZn 7.33 and YbCd 5.7 are ground-state phases, thus revealing that translational symmetry is not a necessary condition for the zero-temperature stability of inorganic solids. Although we found the ScZn 7.33 quasicrystal to be thermodynamically stable, we show on a mixed thermodynamic and kinetic phase diagram that its solidification from the melt is limited by nucleation, which illustrates why even stable materials may be kinetically challenging to grow. Our techniques broadly open the door to first-principles investigations into the structure–bonding–stability relationships of aperiodic materials.

density functional theory↗

Computational Fluid Dynamics analysis of a jet-induced attrition

The flow dynamics, collisional behavior and attrition rates of a grid jet attrition unit are computationally investigated using a kinetic theory-based approach. Profiles of the flow statistical properties are presented and discussed to describe the jet dynamics. Attrition behavior and rates are analyzed from the local reconstruction of the collision energy spectra. The commonly used jet-induced attrition mechanisms which involve the entrainment of the particles into the jet followed by their acceleration and collision at the jet tip is reproduced by the CFD calculations in this work. Here, supported by the local attrition rates, it is shown that the shear-based collisions around the jets, due to counter-current particle axial motions, are an additional jet-induced mechanism generating more fines than at the jet tips.

42 ENGINEERING↗

Diffusion of NiH on Ni(111) and of Ni with Adsorbed H on Ni(111) in the Context of Ni Coarsening in Solid Oxide Cells

In the Ni-based hydrogen electrode of a solid oxide cell (SOC), Ni coarsening is an important degradation mechanism and could be enhanced by NiH diffusing on the Ni particle surfaces. Here, in this work, the average lifetime and diffusion distance of NiH on Ni(111) are computed using density-functional theory and kinetic Monte Carlo methods. It is found that NiH is extremely short-lived and, thus, cannot promote coarsening in the SOC. Also, the diffusion of Ni on Ni(111) is shown to be at most slightly accelerated by H along NiH dissociation paths involving the movement of Ni with a nearby H. Based on this result, coarsening is not likely to be dramatically accelerated by the diffusion of Ni with H on Ni(111). However, support is provided for the experimental procedure of measuring the product of surface coverage and single-species diffusivity of Ni on Ni(111) in a hydrogen-rich atmosphere.

Ni coarsening↗

Combined density functional theory/kinetic Monte Carlo investigation of surface morphology during cycling of Li-Cu electrodes

In this work, we employed a combined density functional theory and kinetic Monte Carlo approach in order to investigate the nanoscale behavior of lithium stripping and plating on lithium-copper nanoalloys. Density functional theory calculations were used to develop an atom-by-atom understanding of the nucleation behavior of lithium on copper surfaces with different surface morphologies. We found that nanometric pits resulted in significantly more favorable lithium deposition when compared to pristine or lightly-defective Cu surfaces, which is in agreement with previous experimental studies. We then used density functional theory and constrained ab initio molecular dynamics to develop interaction and reaction parameters for a larger-scale kinetic Monte Carlo model in order to provide insight into the stripping and plating behavior of nanostructured Li-Cu nanoalloy anodes. Substantial losses in capacity and strong indications of morphological change were observed for a wide range of Cu:Li ratios over 100 cycles. We observed a complex relationship between the nanoparticle composition and the properties of the formed solid-electrolyte interphase.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Density functional theory calculations of the thermodynamic and kinetic properties of point defects in β -U

Density functional theory (DFT) calculations of the thermodynamic and kinetic properties of point defects in the β phase of uranium are reported. Defect energies and entropies were calculated using 2 x 2 x 2 supercells and the Generalized Gradient Approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) for the exchange-correlation potential. Due to computational cost, calculations of the vibrational properties governing entropies were performed by only displacing atoms within (roughly) the 3rd or 4th nearest neighbor shell of the defect, which implicitly assumes that atoms beyond this distance are unaffected by the defect. Migration barriers were estimated by nudged elastic band (NEB) calculations. The low symmetry of the β-U phase (the unit cell is tetragonal and contains 30 atoms) results in many point defect configurations and even more migration pathways. A connectivity map, starting from the most stable point defects, was developed in order to identify the rate-limiting step controlling the net diffusion rate in each crystallographic direction. The uranium self-diffusivity tensor was calculated by combining the defect formation energies, entropies, migration barriers and attempt frequencies. The fastest diffusion rate was determined to be a vacancy mechanism in the z crystallographic direction. The predicted uranium self-diffusivity for this mechanism agrees well with available experimental data. Finally, the diffusion mechanisms and rates identified in this study will inform fuel performance models of swelling and gas evolution.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Theoretically Informed Kinetics (ThInK): Establishing a modern C 0 -C 3 mechanism for combustion modeling

In contrast to the adage “Models are to be used, not believed”, combustion kinetics models have been intended to be predictive in nature. Theoretical chemical kinetics is now understood to provide a firm foundation for the reaction parameters, thereby facilitating predictive simulations of chemical reactivity, even in regimes that are poorly characterized by chemical kinetic and/or combustion experiments. In this study, we describe a theory-informed kinetics model (ThInK) for small molecule combustion chemistry (H 2 and C 1 – C 3 species) that is based on the prodigious use of theoretical predictions for reaction rate coefficients, thermochemistry, and transport parameters. The distinct features of this kinetics model, which was developed over the course of several decades, are illustrated through simulations of flame propagation and auto-ignition.

Energy transfer↗

First-Principles Insights into the Thermocatalytic Cracking of Ammonia-Hydrogen Blends on Fe(110). 2. Kinetics

Ammonia (NH 3 ) is an energy-rich molecule that is routinely synthesized from nitrogen (N 2 ) and hydrogen (H 2 ). NH 3 ’s more favorable physical properties compared to H 2 suggests it may offer a way to more conveniently store, transport, and, when needed, extract H 2 via thermal decomposition. However, the high kinetic barrier and endoergicity to decompose to H 2 and N 2 require high temperatures. The standard reaction free energy indicates nearly 100% thermodynamic conversion to the diatomic molecules only at ~673 K and higher. However, even at these temperatures, a catalyst, e.g., iron (Fe), is needed for favorable kinetic conversion. Here, in this study, we explore via density functional theory the kinetics of NH 3 decomposition on the most stable facet of body-centered cubic Fe, namely, (110), under typical high-temperature and finite-pressure operando conditions. We predict coverage-dependent energetics of elementary surface reactions, often neglected in atomic-scale modeling. From these models, we find the recombinative desorption of adsorbed N as N 2 is rate-determining at 573.15–773.15 K and even at an extreme case of 1173.15 K. From microkinetic modeling, we find that the steady-state turnover frequencies (TOFs) for N 2 and H 2 generation rates (r$_{H_2}$) depend exponentially on temperature. The catalyst achieves a steady-state TOF of 36.4 s –1 and an r$_{H_2}$ of 0.107 μmol cm –2 s –1 for a feed of 1.8 bar NH 3 with 0.2 bar H 2 at 1173.15 K. However, at 773.15 K, with the same feed composition and velocity, the steady-state TOF and r$_{H_2}$ decrease to 0.14 s –1 and 4.10 × 10 –4 μmol cm –2 s –1 , respectively, as the process is significantly hindered by slow N 2 desorption. Although at first glance counterintuitive, our simulations suggest that surface modifications that reduce Fe’s reactivity toward NH x species should enhance its overall NH 3 decomposition activity.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗