Search NASA⌕ Search

SEARCH · Search NASA

Results for “distribution function”

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 OSTI2023

Magnetic structures and excitations in sawtooth olivine chalcogenides Mn 2 SiX 4 (X = S, Se)

The Mn lattice in olivine chalcogenide Mn 2 SiX 4 (X = S, Se) compounds forms a sawtooth, which is of special interest in magnetism owing to the possibility of realizing flat bands in magnon spectra, a key component in magnonics. In this work, we investigate the Mn 2 SiX 4 olivines using magnetic susceptibility, and X-ray and neutron diffraction. We have determined the average and local crystal structures of Mn 2 SiS 4 and Mn 2 SiSe 4 using synchrotron X-ray, neutron diffraction, and X-ray total scattering data followed by Rietveld and pair distribution function analyses. It is found from the pair distribution function analysis that the Mn triangle that constitutes the sawtooth is isosceles in Mn 2 SiS 4 and Mn 2 SiSe 4 . The temperature evolution of magnetic susceptibility of Mn 2 SiS 4 and Mn 2 SiSe 4 shows anomalies below 83 K and 70 K, respectively, associated with magnetic ordering. From the neutron powder diffraction measurements the magnetic space groups of Mn 2 SiS 4 and Mn 2 SiSe 4 are found to be Pnma and Pnm'a', respectively. Here, we find that the Mn spins adopt a ferromagnetic alignment on the sawtooth in both Mn 2 SiS 4 and Mn 2 SiSe 4 but along different crystallographic directions for the S and the Se compounds. From the temperature evolution of Mn magnetic moments obtained from refining neutron diffraction data, the transition temperatures are accurately determined as T N (S) = 83(2) K and T N (Se) = 70.0(5) K. Broad diffuse magnetic peaks are observed in both the compounds, and are prominently seen close to T N , suggesting the presence of a short-range magnetic order. The magnetic excitations studied using inelastic neutron scattering reveal a magnon excitation with an energy corresponding to approximately 4.5 meV in both S and Se compounds. Spin correlations are observed to persist up to 125 K much above the ordering temperature and we suggest the possibility of short-range spin correlations responsible for this.

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

Electrostatic upper-hybrid mode instability driven by a ring electron distribution

Quasi electrostatic fluctuations in the upper-hybrid frequency range are commonly detected in the planetary magnetospheric environment. The origin of such phenomena may relate to the instability driven by a loss-cone feature associated with the electrons populating the dipole-like magnetic field. The present paper carries out a one-dimensional electrostatic particle-in-cell simulation accompanied by a reduced quasilinear kinetic theoretical analysis to investigate the dynamics of the upper-hybrid mode instability driven by an initial ring electron distribution function, which is a form of loss-cone distribution. A favorable comparison is found between the two approaches, which shows that the reduced quasilinear theory, which is grounded in the concept of a model of the particle distribution function that is assumed to maintain a fixed mathematical form except that the macroscopic parameters that define the distribution are allowed to evolve in time, can be an effective tool in the study of plasma instabilities, especially if it is guided by and validated against the more rigorous simulation result.

Physics↗
DOE OSTI2025

Elongated particles in flow: commentary on small-angle scattering investigations

Here, this work thoroughly examines several analytical tools, each possessing a different level of mathematical intricacy, for the purpose of characterizing the orientation distribution function of elongated objects under flow. Our investigation places an emphasis on connecting the orientation distribution to the small-angle scattering spectra measured experimentally. The diverse range of mathematical approaches investigated herein provide insights into the flow behavior of elongated particles from different perspectives and serve as powerful tools for elucidating the complex interplay between flow dynamics and the orientation distribution function.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Whistler waves generated by nongyrotropic and gyrotropic electron beams during asymmetric guide field reconnection

Using a two-dimensional particle-in-cell simulation of asymmetric reconnection with a guide field whose strength is 0.3 times the reconnecting magnetic field, we study electron distribution functions and wave intensities in the diffusion region, focusing on the electron diffusion region (EDR). Wave activities with frequencies below the electron cyclotron frequency are observed, and these are whistler waves propagating almost anti-parallel to the magnetic field. The waves are concentrated near the magnetospheric separatrix away from the X line, but the wave activity also spreads through the EDR near the X line. The reconnection outflows are asymmetric in the outflow direction in the magnetospheric side, and the wave intensity is stronger in the side of the faster electron outflow. We study the whistler waves using the fast Fourier transform, analyses of electron velocity distribution functions, and the dispersion solver calculation. Along the magnetospheric separatrix in the stronger outflow side, highly anisotropic electron beams exist with super-Alfvénic drift speeds. The dispersion analysis shows that there are two modes: a temperature anisotropy mode and a beam mode. Outside the EDR, the whistler wave intensity is highest near the separatrix, but the wave intensity decreases if we move away from the separatrix toward the magnetic neutral line because of the increase in the electron population near zero parallel velocity. In the EDR, in the velocity plane perpendicular to the magnetic field, ring/crescent electron distribution functions are observed. Near the X-line, the wave power is enhanced where nongyrotropic electrons contribute to increase the perpendicular temperature anisotropy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ∼1/Np with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly non-equilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Hybrid particle-spectral method for kinetic plasma simulations

A hybrid model for numerical solutions of the Vlasov–Poisson equations is presented, which blends spectral and particle approaches. The model splits the distribution function for plasma species into both spectral and particle representations in the velocity space to combine the advantages of each approach. The spectral representation leverages asymmetrically weighted Hermite basis, whereas the particle representation leverages the particle-in-cell method. Configuration phase space is decomposed with the Fourier method, which is well suited for periodic problems. We derive conservation equations for mass, momentum, and energy for the proposed combined method. It is shown that the coupling error between the two methods is absent in the semi-discrete setting (not taking into account time discretization). Finally, numerical test cases are presented simulating a weak electron beam interaction with plasma, leading to beam–plasma instability. The initially localized electron beam evolved into a highly non-equilibrium distribution function in the velocity space. A small growth rate and the resonance nature of instability make it difficult to obtain accurate solutions for purely particle methods due to noise, which falls as ~1/$\sqrt{N_p}$ with a number of particles. At the same time, purely spectral methods may require a large number of modes to capture the highly nonequilibrium state of the evolved beam. We show that the hybrid method is well suited for such problems: it reproduces the linear stage as well as nonlinear dynamics with sufficient accuracy using a highly non-equilibrium distribution function.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Unraveling the convoluted and dynamic interphasial mechanisms on Li metal anode

Accurate understanding of the chemistry of solid-electrolyte interphase (SEI) is key to developing new electrolytes for high-energy batteries using lithium metal (Li-0) anodes(1). SEI is generally believed to be formed by the reactions between Li-0 and electrolyte(2,3). However, our new study shows this is not the whole story. Through synchrotron-based X-ray diffraction and pair distribution function analysis, we reveal a much more convoluted formation mechanism of SEI, which receives considerable contributions from electrolyte, cathode, moisture and native surface species on Li-0, with highly dynamic nature during cycling. Using isotope labelling, we traced the origin of LiH to electrolyte solvent, moisture and a new source: the native surface species (LiOH) on pristine Li-0. When lithium accessibility is very limited as in the case of anode-free cells, LiOH develops into plate-shaped large crystals during cycling. Alternatively, when the lithium source is abundant, as in the case of Li||NMC811 cells, LiOH reacts with Li-0 to form LiH and Li2O. While the desired anion-derived LiF-rich SEI is typically found in the concentrated electrolytes or their derivatives, we found it can also be formed in low-concentration electrolyte via the crosstalk effect, emphasizing the importance of formation cycle protocol and opening up opportunities for low-cost electrolyte development.

Polzin, Bryant J.↗
DOE OSTI2021

Physics-informed neural networks for solving forward and inverse flow problems via the Boltzmann-BGK formulation

We report the Boltzmann equation with the Bhatnagar-Gross-Krook collision model (Boltzmann-BGK equation) has been employed to describe multiscale flows, i.e., from the hydrodynamic limit to free molecular flow. In this study we employ physics-informed neural networks (PINNs) to solve forward and inverse problems via the Boltzmann-BGK formulation (PINN-BGK), enabling PINNs to model flows in both the continuum and rarefied regimes. In particular, the PINN-BGK is composed of three sub-networks, i.e., the first for approximating the equilibrium distribution function, the second for approximating the non-equilibrium distribution function, and the third one for encoding the Boltzmann-BGK equation as well as the corresponding boundary/initial conditions. By minimizing the residuals of the governing equations and the mismatch between the predicted and provided boundary/initial conditions, we can approximate the Boltzmann-BGK equation for both continuous and rarefied flows. For forward problems, the PINN-BGK is utilized to solve various benchmark flows given boundary/initial conditions, e.g. Kovasznay flow, Taylor-Green flow, cavity flow, and micro Couette flow for Knudsen number up to 5. For inverse problems, we focus on rarefied flows in which accurate boundary conditions are difficult to obtain. We employ the PINN-BGK to infer the flow field in the entire computational domain given a limited number of interior scattered measurements on the velocity without using the (unknown) boundary conditions. Results for the two-dimensional micro Couette and micro cavity flows with Knudsen numbers ranging from 0.1 to 10 indicate that the PINN-BGK can infer the velocity field in the entire domain with good accuracy. Finally, we also present some results on using transfer learning to accelerate the training process. Specifically, we can obtain a three-fold speedup comparted to the standard training process (e.g., Adam plus L-BFG-B) for the two-dimensional flow problems considered in our work.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Cross-beam energy transfer saturation by ion trapping-induced detuning

The performance of direct-drive inertial confinement fusion implosions relies critically on the coupling of laser energy to the target plasma. Cross beam energy transfer (CBET), the resonant exchange of energy between intersecting laser beams mediated by ponderomotively driven ion-acoustic waves (IAW), inhibits this coupling by scattering light into unwanted directions. Furthermore, the variety of beam intersection angles and varying plasma conditions in an implosion results in IAWs with a range of phase velocities. Here we show that CBET saturates through a resonance detuning that depends on the IAW phase velocity and that results from trapping-induced modifications to the ion distribution functions. For smaller phase velocities, the modifications to the distribution functions can rapidly thermalize in the presence of mid-Z ions, leading to a blueshift in the resonant frequency. For larger phase velocities, the modifications can persist, leading to a redshift in the resonant frequency. Ultimately, these results may reveal pathways towards CBET mitigation and inform reduced models for radiation hydrodynamics codes to improve their predictive capability.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

2021 IC Project Report Viewgraph [Slides]

Here we show that cross-beam energy transfer (CBET) saturates through a resonance detuning that depends on the IAW phase velocity and that results from trapping-induced modifications to the ion distribution functions. For smaller phase velocities, the modifications to the distribution functions can rapidly thermalize in the presence of mid-Z ions, leading to a blueshift in the resonant frequency. For larger velocities, the modifications can persist, leading to a redshift in the resonant frequency.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Post-Newtonian kinetic theory

A kinetic theory for relativistic gases in the presence of gravitational fields is developed in the second post-Newtonian approximation. The corresponding Boltzmann equation is determined from the evolution of the one-particle distribution function with respect to the proper time along the world line of the particle. From the knowledge of the equilibrium Maxwell–Jüttner distribution function in the second post-Newtonian approximation the components of the particle four-flow and energy–momentum tensor are obtained. The Eulerian hydrodynamic equations for the mass density, mass-energy density and momentum density in the second post-Newtonian approximation are determined from the Boltzmann equation. It is shown that the combination of the hydrodynamic equations of mass and mass-energy densities leads to the hydrodynamic equation for the internal energy density in the first post-Newtonian approximation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Alkali hydroxide (LiOH, NaOH, KOH) in water: Structural and vibrational properties, including neutron scattering results

Structural and vibrational properties of aqueous solutions of alkali hydroxides (LiOH, NaOH, and KOH) are computed using quantum molecular dynamics simulations for solute concentrations ranging between 1 and 10M. Element-resolved partial radial distribution functions, neutron and x-ray structure factors, and angular distribution functions are computed for the three hydroxide solutions as a function of concentration. The vibrational spectra and frequency-dependent conductivity are computed from the Fourier transforms of velocity autocorrelation and current autocorrelation functions. Our results for the structure are validated with the available neutron data for 17M concentration of NaOH in water [Semrouni et al., Phys. Chem. Chem. Phys. 21, 6828 (2019)]. We found that the larger ionic radius [rLi+

Chemistry↗
DOE OSTI2020

Progress in modelling fast-ion D-alpha spectra and neutral particle analyzer fluxes using FIDASIM

FIDASIM is a code that models signals produced by charge-exchange reactions between neutrals and ions (both fast and thermal) in magnetically confined plasmas. With the ion distribution function as input, the code predicts the efflux to a neutral particle analyzer diagnostic and the photon radiance of Balmer-alpha light to a fast-ion D α diagnostic, in addition to many other related quantities. A new, parallelized version of the Monte Carlo code FIDASIM has been developed in Fortran90 that is substantially faster than the original interactive data language version. Modified algorithms include more accurate treatments of the time dependent collisional-radiative equations that describe neutral energy levels, of the cloud of ‘halo’ neutrals that surround the injected neutral beam, and of finite Larmor radius effects. Enhanced physics capabilities include modelling ‘passive’ signals from cold edge neutrals, the ability to treat general three-dimensional magnetic confinement configurations, and calculations of diagnostic-specific weight functions that enable tomographic reconstructions of the fast-ion distribution function. Neutral beam attenuation, beam emission, and fast-ion birth profiles are also modelled. Finally, the new algorithms have been successfully validated against experimental data and new features have been tested through benchmarks between two independently developed versions of the code.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Isotope effects in molecular structures and electronic properties of liquid water via deep potential molecular dynamics based on the SCAN functional

Here, Feynman path-integral deep potential molecular dynamics (PI-DPMD) calculations have been employed to study both light (H 2 O) and heavy water (D 2 O) within the isothermalisobaric ensemble. In particular, the deep neural network is trained based on ab initio data obtained from the strongly constrained and appropriately normed (SCAN) exchange-correlation functional. Because of the lighter mass of hydrogen than deuteron, the properties of light water are more influenced by nuclear quantum effect than those of heavy water. Clear isotope effects are observed and analyzed in terms of hydrogen-bond structure and electronic properties of water that are closely associated with experimental observables. The molecular structures of both liquid H 2 O and D 2 O agree well with the data extracted from scattering experiments. The delicate isotope effects on radial distribution functions and angular distribution functions are well reproduced as well. Our approach demonstrates that deep neural network combined with SCAN functional based ab initio molecular dynamics provides an accurate theoretical tool for modeling water and its isotope effects.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Establishing criteria for the transition from kinetic to fluid modeling in hollow cathode analysis

Hollow cathodes for plasma switch applications are investigated via 2D3V particle-in-cell simulations of the channel and plume region. The kinetic nature of the plasma within the channel is dependent on the thermalization rate of electrons, emitted from the insert. When Coulomb collisions occur at a much greater rate than ionization or excitation collisions, the electron energy distribution function rapidly relaxes to a Maxwellian and the plasma within the channel can be described accurately via a fluid model. In contrast, if inelastic processes are much faster than Coulomb collisions, then the electron energy distribution function in the channel exhibits a notable high-energy tail, and a kinetic treatment is required. This criterion is applied to hollow cathodes from the literature, revealing that a fluid approach is suitable for most electric propulsion applications, whereas a kinetic treatment can be more critical to accurate modeling of plasma switches.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Light mesons with one dynamical gluon on the light front

We obtain the light meson mass spectroscopy from the light-front quantum chromodynamics (QCD) Hamiltonian, determined for their constituent quark-antiquark and quark-antiquark-gluon Fock components, together with a three-dimensional confinement. The eigenvectors of the light-front effective Hamiltonian provide a good quality description of the pion electromagnetic form factor, decay constant, and the valence quark distribution functions following QCD scale evolution. We also show that the pion's gluon densities can be probed through the pion-nucleus induced $J/ψ$ production data. Our pion parton distribution functions provide excellent agreement with $J/ψ$ production data from widely different experimental conditions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Molecular dynamics based study on the effects of cation size on the local structure and diffusion in polymerized ionic liquids

We have used coarse-grained molecular dynamics simulations to understand the effect of cation size on ion diffusion in polymerized ionic liquids at temperatures well above the glass transition temperature of the polymers. We investigated dependencies of the diffusion constant on the cation radius and static dielectric constant and interpreted these results in terms of underlying structural changes, and decoupling phenomenon. We have found non-monotonic effects of the cation radius on the diffusion constant, with a maximum at an intermediate radius resulting from two different size dependent effects. Changes in the radial distribution functions characterizing spatial distribution of cations with respect to anions result in a monotonic increase in the coordination number. This yields an increase in the effective interaction energy as a function of cation radius, which can explain the decrease in the diffusion constant. This increase is counteracted at small cation radii by a corresponding increase in decoupling of cation and anion dynamics. In addition, the diffusion constant is found to increase on increasing the static dielectric constant for all radii of the cations. These results highlight the significant role of decoupling in designing polymerized ionic liquids with an enhanced diffusion constant of cations.

Gillespie, Colin [ORNL] (ORCID:0000000320287080)↗
DOE OSTI2024

Spectral analysis and kinetic modeling of radioluminescence in air and nitrogen

Here, in this article we present a quantitative analysis of the second positive system of molecular nitrogen and the first negative system of the molecular nitrogen cation excited in the presence of ionizing radiation. Optical emission spectra of atmospheric air and nitrogen surrounding 210 Po sources were measured from 250 to 400 nm. Multi-Boltzmann and non-Boltzmann vibrational distribution spectral models were used to determine the vibrational temperature and vibrational distribution function of the emitting N 2 (C 3 Π u ) and N 2 + (B 2 Σ + u ) states. A zero-dimensional kinetic model, based on the electron energy distribution function (EEDF) and steady-state excitation and de-excitation of N 2 (X 1 Σ + g ), N 2 + (B 2 Σ + u ), N 2 + (X 2 Σ + g ), N 4 + , O 2 + , and N 2 (C 3 Π u , v), was developed for the prediction of the relative spectral intensity of both the N 2 + (B 2 Σ + u → X 2 Σ + g ) emission band and the vibrational bands of N 2 (C 3 Π u → B 3 Π g ) for comparison with the experimental data.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗