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

A novel conditional formulation of the Vlasov–Ampère equations: a conservative, positivity, asymptotic and Gauss law preserving scheme

We propose a novel reformulation of the Vlasov–Ampère equations for plasmas that reveals discrete symmetries that enables simultaneous conservation of mass, momentum and energy; preservation of Gauss’s law; positivity of the distribution function; and consistency with quasi-neutral asymptotics. The approach employs variable and coordinate transformations to yield a coupled system comprising a modified Vlasov equation and associated moment–field equations. The modified Vlasov equation advances a conditional distribution function that excludes mass, momentum and energy densities, which are instead evolved through moment equations enforcing the relevant symmetries, conservation laws and involution constraints. This reformulation aligns naturally with a recent slow-manifold reduction technique, which separates fast electron time scales and simplifies the treatment of the quasi-neutral limit within the reduced moment–field subsystem. Using this framework, we develop a numerical method for the reduced 1D1V subsystem that, for the first time in the literature, satisfies all key physical constraints while maintaining a quasi-neutral asymptotic behaviour. The advantages of the method are demonstrated on canonical electrostatic test problems, including the multiscale ion acoustic shock wave.

1D1V↗
DOE OSTI2025

Cosmic Reionization on Computers: Statistical Properties of the Distributions of Mean Opacities

Quasar absorption lines provide a unique window to the relationship between galaxies and the intergalactic medium during the Epoch of Reionization. In particular, high redshift quasars enable measurements of the neutral hydrogen content of the universe. However, the limited sample size of observed quasar spectra, particularly at the highest redshifts, hampers our ability to fully characterize the intergalactic medium during this epoch from observations alone. In this work, we characterize the distributions of mean opacities of the intergalactic medium in simulations from the Cosmic Reionization on Computers (CROC) project. We find that the distribution of mean opacities along sightlines follows a non-trivial distribution that cannot be easily approximated by a known distribution. When comparing the cumulative distribution function of mean opacities measurements in subsamples of sample sizes similar to observational measurements from the literature, we find consistency between CROC and observations at redshifts $z\lesssim 5.7$. However, at higher redshifts ($z\gtrsim5.7$), the cumulative distribution function of mean opacities from CROC is notably narrower than those from observed quasar sightlines implying that observations probe a systematically more opaque intergalactic medium at higher redshifts than the intergalactic medium in CROC boxes at these same redshifts. This is consistent with previous analyses that indicate that the universe is reionized too early in CROC simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

The Effect of Spatially Varying Collision Frequency on the Development of the Rayleigh–Taylor Instability

Abstract The Rayleigh–Taylor (RT) instability is ubiquitously observed, yet has traditionally been studied using ideal fluid models. Collisionality can vary strongly across the fluid interface, and previous work demonstrates the necessity of kinetic models to completely capture dynamics in certain collisional regimes. Where previous kinetic simulations used spatially and temporally constant collision frequency, this work presents five-dimensional (two spatial, three velocity dimensions) continuum-kinetic simulations of the RT instability using a more realistic spatially varying collision frequency. Three cases of collisional variation are explored for two Atwood numbers: low to intermediate, intermediate to high, and low to high. The low-to-intermediate case exhibits no RT instability growth, while the intermediate-to-high case is similar to a fluid-limit kinetic case with interface widening biased toward the lower-collisionality region. A novel contribution of this work is the low-to-high collisionality case that shows significantly altered instability growth through an upward movement of the interface and damped spike growth due to increased free-streaming particle diffusion in the lower region. Contributions to the energy flux from the non-Maxwellian portions of the distribution function are not accessible to fluid models and are greatest in magnitude in the spike and regions of low collisionality. Increasing the Atwood number results in greater RT instability growth and reduced upward interface movement. Deviation of the distribution function from Maxwellian is inversely proportional to collision frequency and concentrated around the fluid interface. The linear phase of RT instability growth is well described by theoretical linear growth rates accounting for viscosity and diffusion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Photoproduction of three jets in the CGC: gluon TMDs and dilute limit

We study the process γA → \( q\overline{q}g \) + X in the Color Glass Condensate (CGC) effective theory. After obtaining the cross section, we consider two kinematic limits which are encompassed in our result. In the so-called correlation limit, the vector sum of the transverse momenta of the three outgoing particles is small with respect to the individual transverse momenta; the cross section then simplifies considerably and can be written in a factorized form, sensitive to both the unpolarized and linearly-polarized Weizsäcker-Williams transverse momentum dependent gluon distribution function (gluon TMD). The second limit of the CGC cross section that we consider is the dilute limit, which we obtain after performing a weak-field expansion; we recover a typical linear-regime expression, involving a single unintegrated gluon distribution function. Using numerical simulations of the small-x QCD evolution of the TMDs, we investigate the rapidity dependence of the cross section in the correlation limit.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

A Power Outage Data Informed Resilience Assessment Framework

Catastrophic impacts to power systems due to disruptive events have increased significantly during the last decade. These events highlight the need to develop approaches to assess the resilience of power systems against extreme events. However, the availability of data that capture power system performance during and after disruptive events is scarce. This paper proposes an assessment framework to evaluate the performance aspects of the grid system during extreme outage events using the Environment for Analysis of Geo-Located Energy Information (EAGLE-I) data. EAGLE-I includes information related to the number of impacted customers, duration, and location of power outages in the United States. Statistical analyses were conducted to extract resilient-based outage data and derive probability distribution functions of their impact and recovery characteristics. A list of extreme events is identified based on few predetermined threshold values. Metrics from other power outage assessments were used to measure the characteristics of each event, including impact rate and duration, recovery rate and duration, and impact level. A probability distribution function is obtained for each metric. The obtained results provide a representation of national grid performance during extreme events, which can be applied as a framework to evaluate various resilience enhancement techniques.

EAGLE-I↗
DOE OSTI2022

Atomic structure of an FeCrMoCBY metallic glass revealed by high energy x-ray diffraction

Amorphous bulk metallic glasses with the composition Fe 48 Cr 15 Mo 14 C 15 B 6 Y 2 have been of interest due to their special mechanical and electronic properties, including corrosion resistance, high yield-strength, large elasticity, catalytic performance, and soft ferromagnetism. Here, we apply a reverse Monte Carlo technique to unravel the atomic structure of these glasses. The pair-distribution functions for various atomic pairs are computed based on the high-energy x-ray diffraction data we have taken from an amorphous sample. Monte Carlo cycles are used to move the atomic positions until the model reproduces the experimental pair-distribution function. The resulting fitted model is consistent with our ab initio simulations of the metallic glass. Our study contributes to the understanding of functional properties of Fe-based bulk metallic glasses driven by disorder effects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Local atomic configurations in mechanically alloyed amorphous (FeCoNi) 70 Ti 10 B 20 powders

Here, the atomic structure of amorphous (FeCoNi) 70 Ti 10 B 20 alloy synthesized by mechanical alloying was investigated using high energy synchrotron X-ray diffraction and inverse Monte Carlo simulations of pair distribution functions. Empirical potential structure refinement indicates a chemical short-range order at the length scales of 2.1–2.5 Å via local atomic arrangements forming deformed bcc-like clusters. The structural model obtained was described by bond lengths, coordination numbers, and bond angle distribution functions determined for the first neighbor atoms by x-ray scattering supplemented with 3D Monte Carlo simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Lower hybrid drift instability in nonthermal plasmas

Lower hybrid drift instability (LHDI) is driven by the cross-field current and operates in the vicinity of the lower-hybrid frequency, between the ion- and electron-gyro frequencies, and with wavelengths between the electron and ion thermal gyro radii. The free energy source that drives this instability resides in the density gradient associated with an inhomogeneous plasma. The existing literature on LHDI assumes that the charged particle distribution function is given by a Maxwellian form, but the space plasma is pervasively observed to feature nonthermal characteristics. This paper extends the theory of LHDI to nonthermal plasmas. The generalized theory of LHDI is, thus, applicable to various space plasma environments characterized by nonthermal plasma velocity distribution functions.

Arya, Neetasha (ORCID:0009000236075605)↗
DOE OSTI2024

Particle Acceleration in Relativistic Alfvénic Turbulence

Abstract Strong magnetically dominated Alfvénic turbulence is an efficient engine of nonthermal particle acceleration in a relativistic collisionless plasma. We argue that in the limit of strong magnetization, the type of energy distribution attained by accelerated particles depends on the relative strengths of turbulent fluctuations δ B 0 and the guide field B 0 . If δ B 0 ≪ B 0 , the particle magnetic moments are conserved, and the acceleration is provided by magnetic curvature drifts. Curvature acceleration energizes particles in the direction parallel to the magnetic field lines, resulting in log-normal tails of particle energy distribution functions. Conversely, if δ B 0 ≳ B 0 , interactions of energetic particles with intense turbulent structures can scatter particles, creating a population with large pitch angles. In this case, magnetic mirror effects become important, and turbulent acceleration leads to power-law tails of the energy distribution functions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Physical implications of the extrapolation and statistical bootstrap of nucleon structure function ratios F 2 n F 2 p for mirror nuclei He 3 and H 3

A nuclear physics example of statistical bootstrap is used on the MARATHON nucleon structure function ratio data in the quark momentum fraction regions x B → 0 and x B → 1. The extrapolated F 2 ratio as quark momentum fraction x B → 1 is $\frac{F^n_2}{F^p_2}$ → 0.4 ± 0.05 and this value is compared to theoretical predictions. The extrapolated ratio when x B → 0 favors the simple model of isospin symmetry with the complete dominance of sea quarks at low momentum fraction. At high- x B , the proton quark distribution function ratio d/u is derived from the F 2 ratio and found to be d/u → 1/6. Our extrapolated values for both the $\frac{F^n_2}{F^p_2}$ ratio and the d/u parton distribution function ratio are within uncertainties of perturbative QCD values from quark counting, helicity conservation arguments, and a Dyson-Schwinger equation with a contact interaction model. In addition, it is possible to match the statistical bootstrap value to theoretical predictions by allowing two compatible models to act simultaneously in the nucleon wave function. Finally, one such example is nucleon wave functions composed of a linear combination of a quark-diquark state and a three-valence quark correlated state with coefficients that combine to give the extrapolated F 2 ratio at x B = 1.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Physics-Informed Gaussian Process Inference of Liquid Structure from Scattering Data

We present a nonparametric Bayesian framework to infer radial distribution functions from experimental scattering measurements with uncertainty quantification using nonstationary Gaussian processes. The Gaussian process prior mean and kernel functions are designed to mitigate well-known numerical challenges with the Fourier transform, including discrete measurement binning and detector windowing, while encoding fundamental yet minimal physical knowledge of the liquid structure. We demonstrate uncertainty propagation of the Gaussian process posterior to unmeasured quantities of interest. Experimental radial distribution functions of liquid argon and water with uncertainty quantification are provided as both a proof of principle for the method and a benchmark for molecular models.

Chemical structure↗
DOE OSTI2023

Inclusive-photon production and its dependence on photon isolation in $pp$ collisions at $\sqrt{s}$ = 13 TeV using 139 fb -1 of ATLAS data

Measurements of differential cross sections are presented for inclusive isolated-photon production in $pp$ collisions at a centre-of-mass energy of 13 TeV provided by the LHC and using 139 -1 fb of data recorded by the ATLAS experiment. The cross sections are measured as functions of the photon transverse energy in different regions of photon pseudorapidity. The photons are required to be isolated by means of a fixed-cone method with two different cone radii. The dependence of the inclusive-photon production on the photon isolation is investigated by measuring the fiducial cross sections as functions of the isolation-cone radius and the ratios of the differential cross sections with different radii in different regions of photon pseudorapidity. The results presented in this paper constitute an improvement with respect to those published by ATLAS earlier: the measurements are provided for different isolation radii and with a more granular segmentation in photon pseudorapidity that can be exploited in improving the determination of the proton parton distribution functions. These improvements provide a more in-depth test of the theoretical predictions. Next-to-leading-order QCD predictions from JETPHOX and SHERPA and next-to-next-to-leading-order QCD predictions from NNLOJET are compared to the measurements, using several parameterisations of the proton parton distribution functions. The measured cross sections are well described by the fixed-order QCD predictions within the experimental and theoretical uncertainties in most of the investigated phase-space region.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

A Search for Correlations between Turbulence and Star Formation in LITTLE THINGS Dwarf Irregular Galaxies

Turbulence has the potential for creating gas density enhancements that initiate cloud and star formation (SF), and it can be generated locally by SF. To study the connection between turbulence and SF, we looked for relationships between SF traced by FUV images, and gas turbulence traced by kinetic energy density (KED) and velocity dispersion (v {sub disp}) in the LITTLE THINGS sample of nearby dIrr galaxies. We performed 2D cross-correlations between FUV and KED images, measured cross-correlations in annuli to produce correlation coefficients as a function of radius, and determined the cumulative distribution function of the cross-correlation value. We also plotted on a pixel-by-pixel basis the locally excess KED, v {sub disp}, and H i mass surface density, Σ{sub HI}, as determined from the respective values with the radial profiles subtracted, versus the excess SF rate density Σ{sub SFR}, for all regions with positive excess Σ{sub SFR}. We found that Σ{sub SFR} and KED are poorly correlated. The excess KED associated with SF implies a ∼0.5% efficiency for supernova energy to pump local H i turbulence on the scale of the resolution here, which is a factor of ∼2 too small for all of the turbulence on a galactic scale. The excess v {sub disp} in SF regions is also small, only ∼0.37 km s{sup −1}. The local excess in Σ{sub HI} corresponding to an excess in Σ{sub SFR} is consistent with a H i consumption time of ∼1.6 Gyr in the inner parts of the galaxies. The similarity between this timescale and the consumption time for CO implies that CO-dark molecular gas has comparable mass to H i in the inner disks.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Statistical description of coalescing magnetic islands via magnetic reconnection

The physical picture of interacting magnetic islands provides a useful paradigm for certain plasma dynamics in a variety of physical environments, such as the solar corona, the heliosheath and the Earth's magnetosphere. In this work, we derive an island kinetic equation to describe the evolution of the island distribution function (in area and in flux of islands) subject to a collisional integral designed to account for the role of magnetic reconnection during island mergers. This equation is used to study the inverse transfer of magnetic energy through the coalescence of magnetic islands in two dimensions. We solve our island kinetic equation numerically for three different types of initial distribution: Dirac delta, Gaussian and power-law distributions. The time evolution of several key quantities is found to agree well with our analytical predictions: magnetic energy decays as $\tilde {t}^{-1}$ , the number of islands decreases as $\tilde {t}^{-1}$ and the averaged area of islands grows as $\tilde {t}$ , where $\tilde {t}$ is the time normalised to the characteristic reconnection time scale of islands. General properties of the distribution function and the magnetic energy spectrum are also studied. Finally, we discuss the underlying connection of our island-merger models to the (self-similar) decay of magnetohydrodynamic turbulence.

Physics↗
DOE OSTI2020

Augmenting machine learning of energy landscapes with local structural information

We present a machine learning approach for accurately predicting formation energies of binary compounds in the context of crystal structure predictions. The success of any machine learning model depends significantly on the choice of representation used to encode the relevant physical information into machine-learnable data. We test different representation schemes based on partial radial and angular distribution functions (RDF+ADF) on Al–Ni and Cd–Te structures generated using our genetic algorithm for structure prediction. We observe a remarkable improvement in predictive accuracy upon transitioning from global to atom-centered representations, resulting in a threefold decrease in prediction errors. We show that a support vector regression model using a combination of atomic radial and angular distribution functions performs best at the formation energy prediction task, providing small root mean squared errors of 3.9 meV/atom and 10.9 meV/atom for Al–Ni and Cd–Te, respectively. We test the performance of our models against common traditional descriptors and find that RDF- and ADF-based representations significantly outperform many of those in the prediction of formation energies. The high accuracy of predictions makes our machine learning models great candidates for the exploration of energy landscapes.

Honrao, Shreyas J. (ORCID:0000000166155257)↗
DOE OSTI2020

Fast resonance decays in nuclear collisions

Here, we present a new method to calculate the final decay spectrum of direct resonance decays directly from hydrodynamic fields on a freeze-out surface. The method is based on identifying components of the final particle spectrum that transform in an irreducible way under rotations in the fluid-restframe. Corresponding distribution functions can be pre-computed including all resonance decays. Just a few of easily tabulated scalar functions then determine the Lorentz invariant decay spectrum from each space-time point, and simple integrals of these scalar functions over the freeze-out surface determine the final decay products. This by-passes numerically costly event-by-event calculations of the intermediate resonances. The method is of considerable practical use for making realistic data to model comparisons of the identified particle yields and flow harmonics, and for studying the viscous corrections to the freeze-out distribution function.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

A small-box approach to the local crystal structure of Y 3 NbO 7

The local crystal structure of Y 3 NbO 7 was probed using neutron total scattering. Five different structural models were fit to the experimental pair distribution function of this material up to 15 Å. Fits of defect fluorite (Fm$\overline{3}$m), pyrochlore (Fd$\overline{3}$m), and orthorhombic C weberite models (C222 1 ) were attempted based on prior structural investigations of Y 3 NbO 7 and related fluorite-based materials. Orthorhombic P weberite (P222 1 ) and monoclinic weberite models (P112 1 ) were derived via symmetry reduction of the orthorhombic C weberite structure. All models except monoclinic weberite failed to reproduce the experimental pair distribution function. Misfits were clearly visible for peaks arising from atom pairs within the first coordination shell (defect fluorite and pyrochlore) and beyond (orthorhombic C and P weberite). By contrast, monoclinic weberite (γ ≈ 91.8°) provided an excellent fit, with no major misfits observed. This unit cell features distorted YO 8 dodecahedra, YO 7 monocapped octahedra, and NbO 6 octahedra as building blocks. Comparison between the topologies of the cationic and anionic substructures in monoclinic weberite and its parent orthorhombic C structure revealed displacements of the metal atoms along the [100] direction and a less compact spatial distribution of oxygen atoms along the a axis. In conclusion, the monoclinic weberite model presented herein may serve as a starting point to revisit structural modeling of fluorite-based materials with challenging local structures.

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

Molecular Origins of the Barriers to Proton Transport in Acidic Aqueous Solutions

The self-consistent iterative multistate empirical valence bond (SCI-MS-EVB) method is used to analyze the structure, thermodynamics, and dynamics of hydrochloric acid solutions. The reorientation timescales of irreversible proton transport are elucidated by simulating 0.43 M, 0.85 M, 1.68 M, and 3.26 M HCl solutions at 270 K, 285 K, 300 K, 315 K, and 330 K. The results indicate increased counter-ion pairing with increasing concentration, which manifests itself via a reduced hydronium oxygen–chloride (O*-Cl) structuring in the radial distribution functions. Increasing ionic concentration also reduces the diffusion of the hydrated excess protons, principally by reducing the contribution of the Grotthuss proton hopping (shuttling) mechanism to the overall diffusion process. In agreement with prior experimental findings, a decrease in the activation energy of reorientation timescales was also observed, which is explicitly explained using activated rate theory and an energy-entropy decomposition of the state-averaged radial distribution functions. Furthermore, these results provide atomistic verification of suggestions from recent two-dimensional infrared spectroscopy experiments that chloride anions (as opposed to hydrated excess protons) create entropic barriers to proton transport.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗