Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 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↗
DOE OSTI2026

Tensor network simulations of quasi-GPDs in the massive Schwinger model

Generalized parton distribution functions (GPDs) are off-diagonal light-cone matrix elements that encode the internal structure of hadrons in terms of quark and gluon degrees of freedom. In this work, we present the first nonperturbative study of quasi-GPDs in the massive Schwinger model, quantum electrodynamics in 1+1 dimensions (QED 2 ), within the Hamiltonian formulation of lattice field theory. Quasidistributions are spatial correlation functions of boosted states, which approach the relevant light-cone distributions in the luminal limit. Using tensor networks, we prepare the first excited state in the strongly coupled regime and boost it to close to the light-cone on lattices of up to 400 lattice sites. We compute both quasiparton distribution functions and, for the first time, quasi-GPDs, and study their convergence for increasingly boosted states. In addition, we perform analytic calculations of GPDs in the two-particle Fock-space approximation and in the Reggeized limit, providing qualitative benchmarks for the tensor network results. Our analysis establishes computational benchmarks for accessing partonic observables in low-dimensional gauge theories, offering a starting point for future extensions to higher dimensions, non-Abelian theories, and quantum simulations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Angular momentum eigenstates of the isotropic 3-D harmonic oscillator: Phase-space distributions and coalescence probabilities

The isotropic 3-dimensional harmonic oscillator potential can serve as an approximate description of many systems in atomic, solid state, nuclear, and particle physics. In particular, the question of 2 particles binding (or coalescing) into angular momentum eigenstates in such a potential has interesting applications. Here we compute the probabilities for coalescence of two distinguishable, non-relativistic particles into such a bound state, where the initial particles are represented by generic wave packets of given average positions and momenta. We use a phase-space formulation and hence need the Wigner distribution functions of angular momentum eigenstates in isotropic 3-dimensional harmonic oscillators. These distribution functions have been discussed in the literature before but we utilize an alternative approach to obtain these functions. Along the way, we derive a general formula that expands angular momentum eigenstates in terms of products of 1-dimensional harmonic oscillator eigenstates.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Atomic structure and grain shape evolution of nanodiamond during annealing in oxidizing atmosphere from neutron diffraction and MD simulations

In this work, real space pair-distribution function (PDF) analysis in conjunction with molecular dynamics (MD) simulations was applied to characterize crystal structure, shape and surface structure of diamond nanoparticles annealed in air up to 880 °C. MD simulated models of individual nanodiamond particles of various shapes and sizes were used to calculate theoretical radial distribution functions G(r) for powders with different grain size distributions. The calculated curves, showing distinct differences for the explored parameter space, were compared to the experimental G(r) functions obtained from neutron scattering for annealed detonated nanodiamond samples in a search for a model best describing the actual powders. It was found that as-synthesized nanodiamond grains are terminated by equally abundant (100), (110) and (111) crystallographic faces and their structure does not change upon annealing in air up to 280 °C. The shape and size of grains annealed above 480 °C changes due to surface etching. The etching rate differs between crystals faces and depends on temperature, resulting in the observed differences in grain shape. Below 680 °C the share of the (110) facets gradually increases and at 680 °C the grains become rhombic dodecahedra terminated solely by (110) surfaces. Further increase of the annealing temperature promotes formation of (111) surfaces and at 880 °C the preferred grain shape becomes an octahedron terminated by a variant of (111) surfaces possessing three dangling sp3 bonds per atom. At lower temperatures all grains become smaller and the size distribution moves towards the lower values. Furthermore, as the annealing temperature increases, the smallest grains disappear entirely and the average grain size grows back.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Insights on the structure and properties of sodium iron phosphate glasses from molecular dynamics simulations

Iron phosphate glasses are promising nuclear waste forms while more detailed understanding of their structures and structure-property relations are still needed to better design waste glass compositions. In this work we report studies of three series of sodium iron phosphate (NFP) glasses: 60P 2 O 5 -(40-x)Fe 2 O 3 -xNa 2 O (x = 0→40), (100–2x)P 2 O 5 -xFe 2 O 3 -xNa 2 O (x = 5→17.5) and one with different iron redox ratio, to understand the composition as well as the iron redox effects on the structure and properties of these glasses using molecular dynamics simulations with effective two-body and three-body potentials. Structural analyses, including pair distribution function, bond angle distribution, Q n distribution, and polyhedral connectivity, were performed to obtain in-depth information on short-range and medium-range structural features. The P-O pair distributions showed a first peak splitting with phosphorus-bridging and non-bridging oxygen contributions. This and the average P-O and other cation-oxygen bond distances are in excellent agreement with experiments. The coordination number of P 5+ remained four while that of Fe 3+ increased from 4.30 to 4.72 with decreasing Fe/Na ratio. Polyhedral linkage analysis showed [PO 4 ] units linked with [PO 4 ] and [FeO x ] through corner-sharing while the [PO 4 ]-[FeO x ] linkages become dominant for compositions with Fe 2 O 3 larger than 15 mol%. The effect of iron redox ratio on the structure of NFP glasses was also studied and it was found that bond lengths and coordination numbers were not strongly affected, while the reduction of iron introduced higher network distortions, as evident by O-P-O bond angle and Q n distribution. The glass transition temperature (T g ) showed a monotonic increase with Fe 2 O 3 in the first series, in good agreement with experiments, while those of the second series showed a maximum at P 2 O 5 = 82 mol%. Here, calculated elastic moduli were found to increase with Fe 2 O 3 in the first glass series, which was be explained by the increase of network connectivity, while those of the second series decrease with Fe 2 O 3 due to decrease of P 2 O 5 .

36 MATERIALS SCIENCE↗
DOE OSTI2022

An unconditionally stable, time-implicit algorithm for solving the one-dimensional Vlasov–Poisson system

The development of an implicit, unconditionally stable, numerical method for solving the Vlasov–Poisson system in one dimension using a phase-space grid is presented. The algorithm uses the Crank–Nicolson discretization scheme and operator splitting allowing for direct solution of the finite difference equations. This method exactly conserves particle number, enstrophy and momentum. A variant of the algorithm which does not use splitting also exactly conserves energy but requires the use of iterative solvers. This algorithm has no dissipation and thus fine-scale variations can lead to oscillations and the production of negative values of the distribution function. We find that overall, the effects of negative values of the distribution function are relatively benign. We consider a variety of test cases that have been used extensively in the literature where numerical results can be compared with analytical solutions or growth rates. We examine higher-order differencing and construct higher-order temporal updates using standard composition methods.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Estimating Electron Temperature and Density Using Van Allen Probe Data: Typical Behavior of Energetic Electrons in the Inner Magnetosphere

Abstract The Earth's inner magnetosphere contains multiple electron populations influenced by different factors. The cold electrons of the plasmasphere, warm plasma that contributes to the ring current, and the relativistic plasma of the radiation belts often seem to behave independently. Using omni‐directional flux and energy measurements from the HOPE and Magnetic Electron Ion Spectrometer instruments aboard the Van Allen Probes, we provide a detailed density and temperature description of the inner magnetosphere, offering a comprehensive statistical analysis of the entire Van Allen Probe era. While number density and temperature data at geosynchronous orbit are available, this study focuses on the warm plasma in the inner magnetosphere . Values of density and temperature are extracted by fitting energy and phase space density to obtain the distribution function. The fitted distributions are related to the zeroth and second moments to estimate the number density and temperature. Analysis has indicated that a two Maxwellian fit is sufficient over a wide range of and that there are two independent plasma populations. The more energetic population has a median number density of approximately and a temperature of around 130 keV, with a temperature peak observed between L * = 4 and L * = 4.5. This population is relatively uniform in magnetic local time (MLT). In contrast, the less energetic warm electron population has a median number density of about and a temperature of 7.4 keV. Strong statistical trends in density and temperature across both L * and MLT are presented, along with potential sources driving these variations.

58 GEOSCIENCES↗
DOE OSTI2024

Temperature-induced gradual polymorphic structural crossover in liquid indium, tin, and antimony

Temperature-induced polymorphic structural crossover in metallic liquid indium (In), tin (Sn), and antimony (Sb) is investigated by in situ high-energy x-ray diffraction and ab initio molecular dynamics simulations. The results demonstrate the existence of a temperature-induced reversible and gradual structural crossover in these three liquids, reflecting from both “static” structures, e.g., peak positions of the structure factor and pair distribution function, bond angle distribution, number of tetrahedra, cluster connection, free volume, and dynamical behaviors, including the relaxation time, diffusion coefficient, and non-Gaussian parameter. Further, the main difference from liquid In, Sn, to Sb is that not only the large-sized free volume but also the anharmonicity and the dynamical heterogeneity increase significantly with the appearance of strong covalent bonds. These findings will deepen the understanding of liquid structure and the liquid-liquid transition in metallic liquids.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Measurements of jet cross-section ratios in 13 TeV proton-proton collisions with ATLAS

Measurements of jet cross-section ratios between inclusive bins of jet multiplicity are performed in 140 fb -1 of proton-proton collisions with √s = 13 TeV center-of-mass energy, recorded with the ATLAS detector at CERN’s Large Hadron Collider. These ratios are constructed from double-differential cross-section measurements that are made in bins of jet multiplicity and other observables that are sensitive the energy scale and angular distribution of radiation due to the strong interaction in the final state. Additionally, the scalar sum of the two leading jets’ transverse momenta is measured triple differentially, in bins of the third jet’s transverse momentum and of jet multiplicity. These measurements are unfolded to account for acceptance and detector-related effects. The measured distributions are used to construct ratios of the inclusive jet-multiplicity bins, which have been shown to be sensitive to the strong coupling Formula Presented while being less sensitive than other observables to systematic uncertainties and parton distribution functions. The measured distributions are compared with state-of-the-art QCD calculations, including next-to-next-to-leading-order predictions for two- and three-jet events. These predictions are generally found to model the data well and perform best in bins with a modest requirement on the third jet’s transverse momentum. Significant differences between data and Monte Carlo predictions are observed in events with large rapidity gaps and invariant masses of the leading jet pair. Studies leading to reduced jet energy scale uncertainties significantly improve the precision of this work and are documented herein.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Landau Modes are Eigenmodes of Stellar Systems in the Limit of Zero Collisions

Abstract We consider the spectrum of eigenmodes in a stellar system dominated by gravitational forces in the limit of zero collisions. We show analytically and numerically using the Lenard–Bernstein collision operator that the Landau modes, which are not true eigenmodes in a strictly collisionless system (except for the Jeans unstable mode), become part of the true eigenmode spectrum in the limit of zero collisions. Under these conditions, the continuous spectrum of true eigenmodes in a collisionless system, also known as the Case–van Kampen modes, is eliminated. Furthermore, because the background distribution function in a weakly collisional system can exhibit significant deviations from a Maxwellian distribution function over long times, we show that the spectrum of Landau modes can change drastically even in the presence of slight deviations from a Maxwellian, primarily through the appearance of weakly damped modes that may be otherwise heavily damped for a Maxwellian distribution. Our results provide important insights for developing statistical theories to describe thermal fluctuations in a stellar system, which are currently a subject of great interest for N -body simulations as well as observations of gravitational systems.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Double-sided electron energy analyzer for measurement of non-Maxwellian electron energy distributions

Here, a double-sided electron energy analyzer is developed for studies of magnetic reconnection. It can measure electron energy distribution functions along two directions opposite to each other at the same time. Each side is composed of a floating reference grid, an energy selector grid, and a collector plate. The voltage of the selector grid is swept from –40 to 0 V with respect to the reference grid with a frequency of 1 MHz. This fast sweeping is required to resolve sub-Alfvénic changes in plasma quantities of the Magnetic Reconnection Experiment, where the typical Alfvénic time is a few microseconds. The reliability of the energy analyzer is checked in Maxwellian plasmas away from the reconnection region. In this case, the electron temperature computed from the electron energy distribution function agrees with measurements of a reference triple Langmuir probe. When it is located near the reconnection region, the temperatures of the tail electron population from both sides, facing into the electron flow and facing away from it, exceed the bulk electron temperature measured by the Langmuir probe by a factor of about 2.

47 OTHER INSTRUMENTATION↗