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 739 records · Page 41

Scale Dependence of Multiplier Distributions for Particle Concentration, Enstrophy, and Dissipation in the Inertial Range of Homogeneous Turbulence

Turbulent flows preferentially concentrate inertial particles depending on their stopping time or Stokes number, which can lead to significant spatial variations in the particle concentration. Cascade models are one way to describe this process in statistical terms. Here, we use a direct numerical simulation (DNS) dataset of homogeneous, isotropic turbulence to determine probability distribution functions (PDFs) for cascade multipliers, which determine the ratio by which a property is partitioned into subvolumes as an eddy is envisioned to decay into smaller eddies. We present a technique for correcting effects of small particle numbers in the statistics. We determine multiplier PDFs for particle number, flow dissipation, and enstrophy, all of which are shown to be scale dependent. However, the particle multiplier PDFs collapse when scaled with an appropriately defined local Stokes number. As anticipated from earlier works, dissipation and enstrophy multiplier PDFs reach an asymptote for sufficiently small spatial scales. From the DNS measurements, we derive a cascade model that is used it to make predictions for the radial distribution function (RDF) for arbitrarily high Reynolds numbers, Re, finding good agreement with the asymptotic, infinite Re inertial range theory of Zaichik and Alipchenkov [New J. Phys. 11, 103018 (2009)]. We discuss implications of these results for the statistical modeling of the turbulent clustering process in the inertial range for high Reynolds numbers inaccessible to numerical simulations.

turbulent flows↗

The energy associated with MHD waves generation in the solar wind plasma

Gyrotropic symmetry is usually assumed in measurements of electron distribution functions in the heliosphere. This prevents the calculation of a net current perpendicular to the magnetic field lines. Previous theoretical results derived by one of the authors for a collisionless plasma with isotropic electrons in a strong magnetic field have shown that the excitation of MHD modes becomes possible when the external perpendicular current is non-zero. We consider then that any anisotropic electron population can be thought of as 'external', interacting with the remaining plasma through the self-consistent electromagnetic field. From this point of view any perpendicular current may be due to the anisotropic electrons, or to an external source like a stream, or to both. As perpendicular currents cannot be derived from the measured distribution functions, we resort to Ampere's law and experimental data of magnetic field fluctuations. The transfer of energy between MHD modes and external currents is then discussed.

delaTorre, A.↗

Review of Current Software for Analyzing Total X-ray Scattering Data from Liquids

The popularity of the pair distribution function (PDF) analysis of X-ray total scattering data has steadily grown as access to ex situ synchrotron data has expanded. Due to the broadening of the PDF user community, there is a growing demand for software that can be used to extract PDFs and is accessible to non-expert users. While user-friendly options have been developed over the past decade for fast, streamlined data analysis, care must be taken in both processing the data and understanding any limitations, especially in the case of liquids. In this review, the same scattering data are analyzed using different total X-ray scattering software, in order to compare the accuracy of the extracted structure factors and associated pair distribution functions. The goal is to assess the best practices for extracting the most accurate liquid data for each software package. The importance of absolute normalization and the application of the most appropriate corrections are emphasized via quantitative comparisons between liquid sulfur and water. Additionally, an awareness of the competing conventions used to define the PDF in crystallography and liquids/glasses is crucial for both the downstream analyses of the data and a comparison with the previous results in the literature.

Gallington, Leighanne C. (ORCID:0000000203837522)↗

Joint Factorization of QCD and QED Radiation in Lepton-Hadron Scattering

The factorization theorem plays an important role in the analysis of high energy quantum chromodynamic (QCD) processes, separating the nonperturbative hadronic interaction into the universal parton distribution functions (PDFs) and fragmentation functions (FFs) and the process-dependent interactions into short distance perturbative calculations, with any interference power suppressed. With a virtual photon exchange, lepton-hadron deep inelastic scattering (DIS) provides an electromagnetic hard probe for the partonic structure of colliding hadrons and has played an important role in the development of QCD factorization. However, the collision induced QED radiation can change the momentum of the exchanged but unobserved virtual photon, making the photon-hadron frame, where the factorization formalism for DIS and semi-inclusive DIS (SIDIS) was derived, ill defined. A new analogous factorization approach has been introduced to separate the leading power process-independent QED radiative contributions to the single photon exchange by introducing lepton distribution functions (LDFs) and lepton fragmentation functions (LFFs), while process-dependent effects are perturbatively calculated with large logarithms removed [J. High Energ. Phys. 2021, 157 (2021)]. These LDFs and LFFs are considered global, as they appear in many different interactions, such as e+e-, DIS and SIDIS, so data from experiments can be used to fit and describe these functions across a wide range of lepton scattering. In this work, I will apply this new hybrid factorization approach to lepton-hadron DIS and SIDIS. For DIS, I derive the NLO short distance perturbative contribution to the cross section and demonstrate the effects the QED radiation has on the cross section using this approach using the CTEQ parameterization for the QCD functions. As part of the SIDIS analysis, I study the cross-section in two different kinematic regions: (1) the scattered lepton and observed hadron are not near back-to-back, and (2) they are close to back-to-back, where collinear QCD factorization works for (1) and TMD QCD factorization for (2) while collinear QED factorization works for both. As part of this work, I show the effects on the SIDIS cross section using fixed order calculations for the unpolarized structure function by first showing the effect of the radiative corrections on the main kinematic variables, especially how the internal transverse momentum is significantly correlated to the external angular dependence, and then the unpolarized structure function (or cross section) with matching between the descriptions for low and high transverse momentum. This work will impact the calculations for predictions for data from COMPASS and various Jefferson Lab experiments.

Cammarota, Justin [Univ. of Kentucky, Lexington, K↗

Wide-Field Hubble Space Telescope Observations of the Globular Cluster System in NGC 1399*

We present a comprehensive high spatial resolution imaging study of globular clusters (GCs) in NGC 1399, thecentral giant elliptical cD galaxy in the Fornax galaxy cluster, conducted with the Advanced Camera for Surveys(ACS) aboard theHubble Space Telescope(HST).Using a novel technique to construct drizzled point-spreadfunction libraries for HSTACS data, we accurately determine the fidelity of GC structural parameter measurementsfrom detailed artificial star cluster experiments and show the superior robustness of the GC half-light radius,rh,compared with other GC structural parameters, such as King core and tidal radius. The measurement ofrhfor themajor fraction of the NGC 1399 GC system reveals a trend of increasingrhversus galactocentric distance,Rgal,out to about 10 kpc and a flat relation beyond. This trend is very similar for blue and red GCs, which are found tohave a mean size ratio ofrh,redrh,blue0.820.11 at all galactocentric radii from the core regions of the galaxyout to40 kpc. This suggests that the size difference between blue and red GCs is due to internal mechanismsrelated to the evolution of their constituent stellar populations. Modeling the mass density profile of NGC 1399shows that additional external dynamical mechanisms are required to limit the GC size in the galaxy halo regionstorh2 pc. We suggest that this may be realized by an exotic GC orbit distribution function, an extended darkmatter halo, andor tidal stress induced by the increased stochasticity in the dwarf halo substructure at largergalactocentric distances. We compare our results with the GCrhdistribution functions in various galaxies and findthat the fraction of extended GCs withrh5 pc is systematically larger in late-type galaxies compared with GCsystems in early-type galaxies. This is likely due to the dynamically more violent evolution of early-type galaxies.We match our GCrhmeasurements with radial velocity data from the literature and split the resulting sample at themedianrhvalue into compact and extended GCs. We find that compact GCs show a significantly smaller line-of-sight velocity dispersion,cmp22525 km s1, than their extended counterparts,ext31721 km s1.Considering the weaker statistical correlation in the GCrhcolor and the GCrhRgalrelations, the more significantGC sizedynamics relation appears to be astrophysically more relevant and hints at the dominant influence of theGC orbit distribution function on the evolution of GC structural parameters.

Hubble↗

Scaling and adiabaticity in a rapidly expanding gluon plasma

In this work we aim to gain qualitative insight on the far-from-equilibrium behavior of the gluon plasma produced in the early stages of a heavy-ion collision. It was recently discovered [1] that the distribution functions of quarks and gluons in QCD effective kinetic theory (EKT) exhibit self-similar “scaling” evolution with time-dependent scaling exponents long before those exponents reach their pre-hydrodynamic fixed-point values. In this work we shed light on the origin of this time-dependent scaling phenomenon in the small-angle approximation to the Boltzmann equation. We first solve the Boltzmann equation numerically and find that time-dependent scaling is a feature of this kinetic theory, and that it captures key qualitative features of the scaling of hard gluons in QCD EKT. We then proceed to study scaling analytically and semi-analytically in this equation. We find that an appropriate momentum rescaling allows the scaling distribution to be identified as the instantaneous ground state of the operator describing the evolution of the distribution function, and the approach to the scaling function is described by the decay of the excited states. That is to say, there is a frame in which the system evolves adiabatically. Furthermore, from the conditions for adiabaticity we can derive evolution equations for the time-dependent scaling exponents. In addition to the known free-streaming and BMSS fixed points, we identify a new “dilute” fixed point when the number density becomes small before hydrodynamization. Corrections to the fixed point exponents in the small-angle approximation agree quantitatively with those found previously in QCD EKT and arise from the evolution of the ratio between hard and soft scales.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The effect of wave-particle interactions on the polar winds O(+)

The escape of the polar wind plasma is an important element in the ionosphere-magnetosphere coupling. Both theory and observations indicate that the wave-particle interactions (WPI) play a significant role in the dynamics of ion outflow along open geomagnetic field lines. A Monte Carlo simulation was developed in order to include the effect of the WPI in addition to the factors that are traditionally included in the 'classical' polar wind (i.e. gravity, electrostatic field, and divergence of geomagnetic field lines). The ion distribution function (f(sub j)), as well as the profiles of its moments (density, drift velocity, temperature, etc.) were found for different levels of WPI, that is, for different values of the normalized diffusion rate in the velocity space (D(tilde) (sub j perpendicular). Although the model included O(=), H(+) and eletrons, we presented only the results related to the O(+) ion. We found that (1) both the density and drift velocity of O(+) increased with the WPI strength, and consequently, the O(+) escape flux was enhanced by a factor of up to 10(exp 5), (2) The O(+) ions could be energized up to a few electron volts; (3) for moderate and high levels of WPI D(tilde) (sub perpendicular) (O(+) greater than (tilde) 1, the distribution function f(O(+)) displayed very pronounced conic features at altitudes around 3 R(sub e). Finally, the interplay between the downward body force, the upward mirror force, and the perpendicular heating resulted in the formation of the 'pressure cooker' effect. This phenomena explained some interesting features of our solution, such as, the peak in the O(=) temperature, and the formation of 'ears' and conics for f(O(+)) around 2.5 R(sub e).

Barakat, A. R.↗

Utilizing the maximum likelihood estimator for flow analysis

We explore the possibility of evaluating flow harmonics by employing the maximum likelihood estimator (MLE). For a given finite multiplicity, the MLE simultaneously furnishes estimations for all the parameters of the underlying distribution function while efficiently suppressing the variance of measures. Also, the method provides a means to assess a specific class of mixed harmonics, which is not straightforwardly feasible by the approaches primarily based on particle correlations. The results are analyzed using the Wald, likelihood ratio, and score tests of hypotheses. Besides, the resultant flow harmonics obtained using MLE are compared with those derived using particle correlations and event plane methods. Here, the dependencies of extracted flow harmonics on the multiplicity of individual events and the total number of events are analyzed. It is shown that the proposed approach works efficiently to deal with the deficiency in detector acceptability. Moreover, we elaborate on a fictitious scenario where the event plane is not a well-defined quantity in the distribution function. For the latter case, the MLE is shown to largely perform better than the two-particle correlation estimator. In this regard, one concludes that the MLE furnishes a meaningful alternative to the existing approaches for flow analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

A review of the kinetic detail required for accurate predictions of normal shock waves

Several aspects of the kinetic models used in the collision phase of Monte Carlo direct simulations have been studied. Accurate molecular velocity distribution function predictions require a significantly increased number of computational cells in one maximum slope shock thickness, compared to predictions of macroscopic properties. The shape of the highly repulsive portion of the interatomic potential for argon is not well modeled by conventional interatomic potentials; this portion of the potential controls high Mach number shock thickness predictions, indicating that the specification of the energetic repulsive portion of interatomic or intermolecular potentials must be chosen with care for correct modeling of nonequilibrium flows at high temperatures. It has been shown for inverse power potentials that the assumption of variable hard sphere scattering provides accurate predictions of the macroscopic properties in shock waves, by comparison with simulations in which differential scattering is employed in the collision phase. On the other hand, velocity distribution functions are not well predicted by the variable hard sphere scattering model for softer potentials at higher Mach numbers.

Muntz, E. P.↗

Method of Distributions for Two‐Phase Flow in Heterogeneous Porous Media

Abstract Multiscale heterogeneity and insufficient characterization data for a specific subsurface formation of interest render predictions of multi‐phase fluid flow in geologic formations highly uncertain. Quantification of the uncertainty propagation from the geomodel to the fluid‐flow response is typically done within a probabilistic framework. This task is computationally demanding due to, for example, the slow convergence of Monte Carlo simulations (MCS), especially when computing the tails of a distribution that are necessary for risk assessment and decision‐making under uncertainty. The frozen streamlines method (FROST) accelerates probabilistic predictions of immiscible two‐phase fluid flow problems; however, FROST relies on MCS to compute the travel‐time distribution, which is then used to perform the transport (phase saturation) computations. To alleviate this computational bottleneck, we replace MCS with a deterministic equation for the cumulative distribution function (CDF) of travel time. The resulting CDF‐FROST approach yields the CDF of the saturation field without resorting to sampling‐based strategies. Our numerical experiments demonstrate the high accuracy of CDF‐FROST in computing the CDFs of both saturation and travel time. For the same accuracy, it is about 5 and 10 times faster than FROST and MCS, respectively.

Yang, Hyung Jun↗

The effect of wave-particle interactions on the polar wind: Preliminary results

A Monte Carlo simulation was developed in order to study the effects of wave-particle interactions (WPI) on the plasma outflow in the polar wind. The simulation also considered the other mechanisms included in the `classical' polar wind studies such as gravity, the polarization electrostatic field, and the divergence of geomagnetic field lines. Although the plasma consisted of electrons, H(+) and O(+) ions, we emphasized the behavior of H(+) in this preliminary study. The ion distribution function, as well as the profiles of its moments (density, drift velocity, temperatures, etc.) were found for different levels of WPI, that is, for different values of normalized diffusion rates in the velocity space (normalized D(sub perpendicular) H(+)). We found that as the WPI strength increases: (1) the ion drift velocity increases and its density decreases; (2) the perpendicular temperature T(sub perpendicular) (H(+)) increases; (3) the parallel temperature T(sub parallel) (H(+)) first decreases and then increases due to the balance between the parallel adiabatic cooling and the transfer of the energy from the perpendicular to the parallel direction; and (4) the temperature anisotropy ((T(sub parallel) (H(+))/T(sub perpendicular) (H(+)) is reduced and even reversed in some cases. For strong WPI (normalized D(sub perpendicular) (H(+) much greater than 1), the ion distribution function shows weak conic features at high altitudes.

Barakat, A. R.↗

Progress on the National Solar Radiation Data Base (NSRDB): A New DNI Computation

This study introduces a new technique to compute direct normal irradiance (DNI) for improving the National Solar Radiation Data Base (NSRDB). A finite-surface integration algorithm is developed to compute solar radiation in differential solid angles and efficiently infer its contribution to a surface perpendicular to the solar direction. A lookup table of cloud bi-directional transmittance distribution function (BTDF) is developed by use of the discrete ordinates radiative transfer (DISORT) model for possible solar and observing directions and various cloud optical and microphysical properties. In each solar incident direction, DNI is given by the cloud BTDFs from approximately 200 differential solid angles. The simulated DNI is calibrated and evaluated using surface observations at the National Renewable Energy Laboratory's (NREL's) Solar Energy Research Laboratory (SRRL) and the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) facility.

41 EE - Solar Energy Technologies Office (EE-4S)↗

Inertial and anisotropic pressure effects on cross-field electron transport in low-temperature magnetized plasmas

In this paper, a one-dimensional (1D) particle-in-cell Monte Carlo collision (PIC-MCC) model is developed to investigate the effects of anisotropic pressure and inertial terms due to non-Maxwellian velocity distribution functions on cross-field electron transport. The conservation of momentum is evaluated by taking the moments of the first-principles gas-kinetic equation. A steady-state discharge is obtained without any low-frequency ionization oscillations by considering an anomalous electron scattering profile. Further, the results obtained from the 1D PIC-MCC model are compared with fluid models, including the quasi-neutral drift-diffusion (DD), non-neutral DD, and full fluid moment models. The discharge current obtained from the PIC-MCC model is in good agreement with the fluid models. The cross-field electron transport due to the inertial terms, i.e. the gradient of axial and azimuthal drift, is evaluated. Moreover, PIC-MCC simulation results show non-zero, anisotropic, off-diagonal pressure tensor terms due to asymmetric non-Maxwellian electron velocity distribution function, potentially contributing to cross-field electron transport.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Electron collection by a highly positive satellite in the ionosphere: Test aprticle simulation

Collection of electrons by a satellite at a relatively high potential in low Earth orbit, including the effects of the satellite's orbital motion, remains a poorly understood phenomenon. Using a test particle simulation in which charged particle motion is tracked in prescribed electric fields, we calculate here the current collection and the nature of the energy distribution function of the electrons collected by the satellite, including the effects of the satellite orbital motion. Calculations of the collected current without the orbital motion show an excellent agreement with the current from the Parker-Murphy model, but with the orbital motion of the satellite the current collection is enhanced, the degree of enhancement depending on the size of the satellite sheath extending along the magnetic field line. In the latter case, the flow of electrons around the satellite shows some interesting behaviors including azimuthal E x B drift around the satellite, axial trapping along the magnetic field, and formation of field-aligned flow of electrons in the wake region. The total energy of the collected electrons is approximately equal e phi(sub 0), where phi(sub 0) is the satellite potential, but the partition of the energy into components parallel (W(sub parallel)) and perpendicular (W(sub perpendicular)) to the magnetic field shows interesting features. For the magnetic field along the polar axis Z, the energy distributions near the poles are perfectly field aligned, that is, the parallel energy (W(sub parallel) = e phi(sub 0). The perpendicular component (W(sub perpendicular)) progressively increases toward the 'equator' of the spherical satellite. For a sufficiently large bias potential phi(sub 0), the equatorial energy distributions f(W(sub parallel)) and f(W(sub perpendicular)) show fine structures having multiple peaks. These features of the distribution functions result from the chaotic motion of electrons and the associated coupling between the electron motions parallel and perpendicular to the magnetic field in strongly nonuniform electric fields near the satellite.

Singh, Nagendra↗

Langmuir Probe and Mass Spectroscopic Measurements in Inductively Coupled CF4 Plasmas

Abstract Electron and ion energy distribution functions and other plasma parameters such as plasma potential (V(sub p)) , electron temperature (T(sub e)), and electron and ion number densities (n (sub e) and n(sub i)) in low pressure CF4 plasmas have been measured. The experiments were conducted in a GEC cell using an inductively coupled plasma (ICP) device powered by a 13.56 MHz radio-frequency (rf) power source. The measurements were made at 300 W of input rf power at 10, 30 and 50 mTorr gas pressures. Langmuir probe measurements suggest that n(sub e), n(sub i) and V(sub p) remain constant over 60% of the central electrode area, beyond which they decrease. Within the limits of experimental error (+/- 0.25 eV), T(sub e) remains nearly constant over the electrode area. T(sub e) and V(sub p) increase with a decrease in pressure. n(sub e) and n(sub i) are not affected as significantly as T(sub e) or V(sub p) by variation in the gas pressure. The electron energy distribution function (EEDF) measurements indicate a highly non-Maxwellian plasma. CF3+ is the most dominant ion product of the plasma, followed by CF2+ and CF+. The concentrations of CF2+ and CF+ are much larger than that is possible from direct electron impact ionization of the parent gas. The cross-section data suggest that the direct electron impact ionization of fragment neutrals and negative ion production by electron attachment may be responsible for increase of the minor ions.

Rao, M. V. V. S.↗

Acceleration of interstellar pickup ions in the disturbed solar wind observed on Ulysses

Acceleration of interstellar pickup H(+) and He(+) as well as of solar wind protons and alpha particles has been observed on Ulysses during the passage of a corotating interaction region (CIR) at approximately 4.5 AU. Injection efficiencies for both the high thermal speed interstellar pickup ions (H(+) and He(+)) and the low thermal speed solar wind ions (H(+) and He(++) are derived using velocity distribution functions of protons, pickup He(+) and alpha particles from less than 1 to 60 keV/e and of ions (principally protons) above approximately 60 keV. The observed spatial variations of the few keV and the few hundred keV accelerated pickup protons across the forward shock of CIR indicate a two stage acceleration mechanism. Thermal ions are first accelerated to speeds of 3 to 4 times the solar wind speed inside the CIR, presumably by some statistical mechanism, before reaching higher energies by a shock acceleration process. Our results also indicate that (1) the injection efficiencies for pickup ions are almost 100 times higher than they are for solar wind ions, (2) pickup H(+) and He(+) are the two most abundant suprathermal ion species and they carry a large fraction of the particle thermal pressure, (3) the injection efficiency is highest for protons, lowest for He(+), and intermediate for alpha particles, (4) both H(+) and He(+) have identical spectral shapes above the cutoff speed for pickup ions, and (5) the solar wind frame velocity distribution function of protons has the form F(w) = F(sub o)w(sup -4) for 1 is less than w is less than approximately 5, where w is the ion speed divided by the solar wind speed. Above w approximately 5-10 the proton spectrum becomes steeper. These results have important implications concerning acceleration of ions by shocks and CIRs, acceleration of anomalous cosmic rays, and particle dynamics in the outer heliosphere.

Gloeckler, G.↗

Theory of flux anisotropies in a guiding center plasma

The one particle distribution function f on the scale of the bounce motion of particles in a magnetic field B is considered. The Vlasov equation is expanded through O(epsilon) in the adiabatic parameter which is the ratio of particle gyroradius to scale length of the magnetic field. Because f is directly proportional to particle flux differential in kinetic energy and solid angle, f is in principle measurable in space experiments, and the analysis is tailored to be explicitly applicable to space problems. To O(1), f is gyrotropic; its first velocity moment is (if non-vanishing) parallel to B, and hence macroscopic parallel flow is included in this term. The O(epsilon) contribution is non-gyrotropic and macroscopic flow parallel to B plus additional parallel flow results from these terms. The degree of non-gyrotropy and the amount of cross-field macroscopic flow depend on the perpendicular component of the electric field, on curvature and shear in the magnetic field, and on the spatial gradient, pitch angle derivative, and speed derivative of the lowest order distribution function.

Birmingham, T. J.↗

Evolution and interplay of lithium metal interphase components revealed by experimental and theoretical studies

ABSTRACT: Lithium metal batteries (LMB) have high energy densities and are crucial for clean energy solutions. The characterization of lithium metal interphase is fundamentally and practically important but technically challenging. Taking advantage of synchrotron x-ray which has the unique capability of analyzing crystalline/amorphous phases quantitative-ly with statistical significance, we study the composition and dynamics of LMB interphase for a newly developed im-portant LMB electrolyte that is based on fluorinated ether. Pair distribution function analysis revealed the sequential role of anion and solvent in interphase formation during cycling. The relative ratio between Li2O and LiF first increases and then decreases during cycling, suggesting suppressed Li2O formation in both initial and long extended cycles. Theoretical studies revealed that in initial cycles, this is due to the energy barrier in many-electron transfer. In long extended cycles, the anion decomposition product Li2O encourages solvent decomposition by facilitating solvent adsorption on Li2O which is followed by concurrent depletion of both. This work highlights the important role of Li2O in transitioning from anion-derived interphase to a solvent-derived one.

Tan, Sha↗