SEARCH · Search NASA

Results for “distribution function moments”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7

All sources NASA NTRS DOE OSTI

NASA NTRS1994

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.↗

NASA NTRS1957

Theoretical Calculation of the Power Spectra of the Rolling and Yawing Moments on a Wing in Random Turbulence

The correlation functions and power spectra of the rolling and yawing moments on an airplane wing due to the three components of continuous random turbulence are calculated. The rolling moments to the longitudinal (horizontal) and normal (vertical) components depend on the spanwise distributions of instantaneous gust intensity, which are taken into account by using the inherent properties of symmetry of isotropic turbulence. The results consist of expressions for correlation functions or spectra of the rolling moment in terms of the point correlation functions of the two components of turbulence. Specific numerical calculations are made for a pair of correlation functions given by simple analytic expressions which fit available experimental data quite well. Calculations are made for four lift distributions. Comparison is made with the results of previous analyses which assumed random turbulence along the flight path and linear variations of gust velocity across the span.

Eggleston, John M↗

NASA NTRS1994

Kinetic and thermodynamic properties of a convecting plasma in a two-dimensional dipole field

Charged particle guiding center motion is considered in the magnetic field of a two-dimensional ('line') dipole on which is superimposed a small, static, perpendicular electric field. The parallel equation of motion is that of a simple harmonic oscillator for cos theta, the cosine of magnetic colatitude theta. Equations for the perpendicular electric and magnetic drifts are derived as well as their bounce-averaged forms. The latter are solved to yield a bounce-averaged guiding center trajectory, which is the same as that obtained from conversation of magnetic moment mu, longitudinal invariant J, and total (kinetic plus electrostatic) energy K. The algebraic simplicity of the trajectory equations is also manifest in the forms of the invariants. An interesting result is that guiding centers drift in such a way that they preserve the values of their equatorial pitch angles and (equivalently) mirror latitudes. The most general Maxwellian form of the equilibrium one-particle distribution function f is constructed from the invariants, and spatially varying density and pressure moments, parallel and perpendicular to the magnetic field, are identified. Much of the paper deals with the more restricted problem in which f is specified as a bi-Maxwellian over a straight line of finite length in the equatorial plane of the dipole and perpendicular to field lines. This might be thought of as specifying a cross-tail ion injection source; our formalism then describes the subsequent spatial development. The distribution away from the source is a scaled bi-Maxwellian but one that is cut off at large and small kinetic energies, which depend on position. Density and pressure components are reduced from the values they would have if the total content of individual flux tubes convected intact. The equatorial and meridional variations of density and pressure components are examined and compared systematically for the isotropic and highly anisotropic situations. There appears to be little qualitative difference due to anisotropy. An anisotropy measure is defined, and its spatial variation determined as a signature of possible MHD instability. Extreme values are found, larger than at the source, but the plasma beta in such regions is probably so low as to render the effect inconsequential energetically. Finally, the possible consequence of 'nonadia- batic' pressure profiles on electrostatic interchanges is considered, and a boundary delineating stabilizing and destabilizing regions determined.

Huang, T. S.↗

DOE OSTI2021

Constructing angular distributions of neutrinos in core-collapse supernovae from zeroth and first moments calibrated by full Boltzmann neutrino transport

Two-moment neutrino transport methods have been widely used for developing theoretical models of core-collapse supernovae (CCSN), since they substantially reduce the computational burden inherent in the multidimensional neutrino-radiation hydrodynamical simulations. The approximation, however, comes at a price; the detailed structure of angular distribution of neutrinos is sacrificed, that is the main drawback of this approach. In this paper, we develop a novel method by which to construct angular distributions of neutrinos from the zeroth and first angular moments. In our method, the angular distribution is expressed with two quadratic functions of the neutrino angle in a piecewise fashion. We determine the best parameters in the fitting function by comparing to the neutrino data in a spherically symmetric CCSN model with full Boltzmann neutrino transport. We demonstrate the capability of our method by using our recent 2D CCSN model. In this work, we find that the essential features of the angular distributions can be well reconstructed, whereas the angular distributions of incoming neutrinos tend to have large errors that increase with flux factor ($\textit{κ}$). This issue originates from the insensitiveness of incoming neutrinos to $\textit{κ}$, that is an intrinsic limitation in moment methods. Based on the results of the demonstration, we assess the reliability of electron- neutrinos lepton number (ELN)-crossing searches with the two-moment neutrino transport. This analysis is complementary to our [L. Johns and H. Nagakura arXiv:2104.04106] that scrutinizes the limitation of crossing searches with a few moments. We find that the systematic errors of angular distributions for incoming neutrinos lead to misjudgements of the crossing at $\textit{κ}$ ≳ 0.5. This casts doubt on the results of ELN-crossing searches based on two-moment methods in some previous studies.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2022

Limitations of Separate Cloud and Rain Categories in Parameterizing Collision-Coalescence for Bulk Microphysics Schemes

Warm rain collision-coalescence has been persistently difficult to parameterize in bulk microphysics schemes. We use a flexible bulk microphysics scheme with bin scheme process parameterizations, called AMP, to investigate reasons for the difficulty. AMP is configured in a variety of ways to mimic bulk schemes and is compared to simulations with the bin scheme upon which AMP is built. We find that an important limitation in traditional bulk schemes is the use of separate cloud and rain categories. When the drop size distribution is instead represented by a continuous distribution, the simulation of cloud-to-rain conversion is substantially improved. We also find large sensitivity to the threshold size to distinguish cloud and rain in traditional schemes; substantial improvement is found by decreasing the threshold from 40 to 25 μm. Neither the use of an assumed functional form for the size distribution nor the choice of predicted distribution moments has a large impact on the ability of AMP to simulate rain production. When predicting four total moments of the liquid drop size distribution, either with a traditional two-category, two-moment scheme with a reduced size threshold, or a four-moment single-category scheme, errors in the evolution of mass and the cloud size distribution are similar, but the single-category scheme has a substantially better representation of the rain size distribution. Optimal moment combinations for the single-category approach are investigated and appear to be linked more to the information content they provide for constraining the size distributions than to their correlation with collision-coalescence rates.

54 ENVIRONMENTAL SCIENCES↗

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 OSTI2026

Gradient Flow for Parton Distribution Functions: First Application to the Pion

Parton distribution functions (PDFs) are central to precision QCD phenomenology. Their Mellin moments can be computed on the lattice, but direct determinations using local operators, besides ⟨𝑥⟩, face severe challenges from reduced hypercubic symmetry, limiting results to the lowest moments. A recently proposed method resolves these issues using gradient flow. We demonstrate the efficacy of this method by computing ratios of flavor nonsinglet pion PDF moments up to ⟨𝑥 5 ⟩ on four lattice spacings at 𝑚 𝜋 ≃ 411 MeV. The moments and reconstructed PDF agree quantitatively with recent phenomenological extractions.

Francis, Anthony [National Yang Ming Chiao Tung Un↗

NASA NTRS1972

Kinetic theory analysis of rarefied gas flow through finite length slots

An analytic study is made of the flow a rarefied monatomic gas through a two dimensional slot. The parameters of the problem are the ratios of downstream to upstream pressures, the Knudsen number at the high pressure end (based on slot half width) and the length to slot half width ratio. A moment method of solution is used by assuming a discontinuous distribution function consisting of four Maxwellians split equally in angular space. Numerical solutions are obtained for the resulting equations. The characteristics of the transition regime are portrayed. The solutions in the free molecule limit are systematically lower than the results obtained in that limit by more accurate numerical methods.

Raghuraman, P.↗

NASA NTRS1981

Two-stream Maxwellian kinetic theory of cloud droplet growth by condensation

A new growth rate formula (NGRF) is developed for the rate of growth of cloud droplets by condensation. The theory used is a modification of the Lees-Shankar theory in which the two-stream Maxwellian distribution function of Lees is used in Maxwell's method of moments to determine the transport of water vapor to and heat away from the droplet. Boundary conditions at the droplet are the usual conditions set in terms of accommodation coefficients, and the solution passes smoothly into diffusion flow in the far region. Comparisons are given between NGRF and the conventional formula showing close agreement (approximately 0.1%) for large radii with significant difference (approximately 5%) for small radii (not greater than 1 micron). Growth times for haze droplets in a Laktionov chamber are computed.

Robinson, N. F.↗

NASA NTRS2022

Limitations of Separate Cloud and Rain Categories in Parameterizing Collision-Coalescence for Bulk Microphysics Schemes

autoconversion↗

NASA NTRS2005

Plasma sheet turbulence observed by Cluster II

Cluster fluxgate magnetometer (FGM) and ion spectrometer (CIS) data are employed to analyze magnetic field fluctuations within the plasma sheet during passages through the magnetotail region in the summers of 2001 and 2002 and, in particular, to look for characteristics of magnetohydrodynamic (MHD) turbulence. Power spectral indices determined from power spectral density functions are on average larger than Kolmogorov's theoretical value for fluid turbulence as well as Kraichnan's theoretical value for MHD plasma turbulence. Probability distribution functions of the magnetic fluctuations show a scaling law over a large range of temporal scales with non-Gaussian distributions at small dissipative scales and inertial scales and more Gaussian distribution at large driving scales. Furthermore, a multifractal analysis of the magnetic field components shows scaling behavior in the inertial range of the fluctuations from about 20 s to 13 min for moments through the fifth order. Both the scaling behavior of the probability distribution functions and the multifractal structure function suggest that intermittent turbulence is present within the plasma sheet. The unique multispacecraft aspect and fortuitous spacecraft spacing allow us to examine the turbulent eddy scale sizes. Dynamic autocorrelation and cross correlation analysis of the magnetic field components allow us to determine that eddy scale sizes fit within the plasma sheet. These results suggest that magnetic field turbulence is occurring within the plasma sheet resulting in turbulent energy dissipation.

Weygand, James M.↗

NASA NTRS1991

Magnetic moment distribution of magnetic cataclysmic variables

A simulation study is made of the relative numbers of the AM Herculis binaries and the intermediate polars as a function of the orbital period using random variables subject to suitable constraints to describe the various parameters. It is shown that the observations can be matched by a single distribution in the magnetic moment equals 0.7 +/- 0.3. For such an ensemble, the intermediate polars are distributed in the log(Porb) - log(Ps) diagram about the critical disk line but with a larger scatter than observed.

Wu, Kinwah↗

DOE OSTI2021

Measurement of charged particle multiplicity distributions in DIS at HERA and its implication to entanglement entropy of partons

Charged particle multiplicity distributions in positron-proton deep inelastic scattering at a centre-of-mass energy √s = 319 GeV are measured. The data are collected with the H1 detector at HERA corresponding to an integrated luminosity of 136 pb -1 . Charged particle multiplicities are measured as a function of photon virtuality Q 2 , inelasticity y and pseudorapidity η in the laboratory and the hadronic centre-of-mass frames. Predictions from different Monte Carlo models are compared to the data. The first and second moments of the multiplicity distributions are determined and the KNO scaling behaviour is investigated. The multiplicity distributions as a function of Q 2 and the Bjorken variable x bj are converted to the hadron entropy S hadron , and predictions from a quantum entanglement model are tested.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2026

Rapidity-dependent spin decomposition of the nucleon

We revisit the two-dimensional Fourier transform of generalized parton distributions (GPDs) at nonzero skewness. At 𝜂 = 0 it reduces to the standard impact-parameter density, while at 𝜂 ≠ 0 it is an off-forward amplitude that we interpret as a genuine parton–nucleon correlation. Its overall strength (the transverse-plane integral of the density) is fixed by the GPD at the kinematic point 𝑡 =−𝑐 𝜂 =−4⁢𝜂 2 ⁢𝑚$^{2}_{𝑁}$/(1 − 𝜂 2 ) and decreases monotonically with the rapidity gap Δ⁢𝑦 = ln⁡[(1+𝜂)/(1−𝜂)] = 2 artanh⁡(𝜂). This rapidity dependence implies rapidity-modified Ji identities that connect helicity, orbital, and total angular momenta of the correlation in closed form. To quantify these effects, we construct leading-twist quark and gluon GPDs in a string-based conformal framework: conformal moments are parametrized by linear open- and closed-string Regge trajectories with slopes constrained by parton distribution functions (PDFs), hadron/glueball spectroscopy, and form-factor data, and GPDs are reconstructed over the full (𝑥,𝜂,𝑡) domain by Mellin-Barnes inversion with next-to-leading order evolution. We find qualitative agreement (and fair quantitative agreement within quoted uncertainties) for several moments and selected nonsinglet 𝑥-space channels at 𝜇 = 2 GeV when compared with lattice QCD, while we also identify channels with visible tension and discuss likely sources (PDF priors and 𝑡-slope systematics).

Gauge-gravity dualities↗

DOE OSTI2023

Generalized parton distributions through universal moment parameterization: non-zero skewness case

We present the first global analysis of generalized parton distributions (GPDs) combing lattice quantum chromodynamics (QCD) calculations and experiment measurements including global parton distribution functions (PDFs), form factors (FFs) and deeply virtual Compton scattering (DVCS) measurements. Following the previous work where we parameterize GPDs in terms of their moments, we extend the framework to allow for the global analysis at non-zero skewness. Together with the constraints at zero skewness, we fit GPDs to global DVCS measurements from both the recent JLab and the earlier Hadron-Electron Ring Accelerator (HERA) experiments with two active quark flavors and leading order QCD evolution. With certain choices of empirical constraints, both sea and valence quark distributions are extracted with the combined inputs, and we present the quark distributions in the proton correspondingly. We also discuss how to extend the framework to accommodate more off-forward constraints beyond the small ξ expansion, especially the lattice calculated GPDs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS2004

Interplanetary Radiation and Internal Charging Environment Models for Solar Sails

A Solar Sail Radiation Environment (SSRE) model has been developed for characterizing the radiation dose and internal charging environments in the solar wind. The SSRE model defines the 0.01 keV to 1 MeV charged particle environment for use in testing the radiation dose vulnerability of candidate solar sail materials and for use in evaluating the internal charging effects in the interplanetary environment. Solar wind and energetic particle instruments aboard the Ulysses spacecraft provide the particle data used to derive the environments for the high inclination 0.5 AU Solar Polar Imager mission and the 1.0 AU L1 solar sail missions. Ulysses is the only spacecraft to sample high latitude solar wind environments far from the ecliptic plane and is therefore uniquely capable of providing the information necessary for defining radiation environments for the Solar Polar Imager spacecraft. Cold plasma moments are used to derive differential flux spectra based on Kappa distribution functions. Energetic particle flux measurements are used to constrain the high energy, non-thermal tails of the distribution functions providing a comprehensive electron, proton, and helium spectra from less than 0.01 keV to a few MeV.

Minow, Joseph I.↗

DOE OSTI2026

Proton isovector helicity PDF at NNLO and the twist-3 moment $\tilde{d}$ 2 from lattice QCD at physical quark masses

We present a lattice quantum chromodynamics calculation of the 𝑥-dependent isovector quark helicity parton distribution function (PDF) of the proton in the large momentum effective theory (LaMET) framework. Through operator product expansion (OPE) we also extract the $\tilde{d}$ 2 moment of the twist-3 PDF 𝑔 𝑇 ⁡(𝑥) for the first time in the $\overline{MS}$ scheme, which is proportional to the average color Lorentz force experienced by the quark in the proton. This calculation is performed on a lattice of spacing 𝑎 =0.076 fm at physical quark masses. The quasi-PDF matrix elements are measured in proton states boosted to momenta 𝑃 𝑧 ={0,0.25,1.02,1.53} GeV. We first extract the lowest few helicity PDF moments from the renormalization-group (RG) invariant ratios of the matrix elements with OPE. Combined with the matrix elements relevant for 𝑔 𝑇 ⁡(𝑥), we obtain $\tilde{d}$$^{u-d}_2$⁡(2 GeV) =0.0024⁢(46) at next-to-leading order in $\overline{MS}$. Then, the helicity quasi-PDF matrix elements are renormalized in the hybrid scheme with linear renormalon resummation and Fourier transformed to the 𝑥-space after an asymptotic extrapolation. The quasi-PDF is perturbatively matched to the $\overline{MS}$ PDF with RG and threshold resummations at next-to-leading power and next-to-next-to-leading logarithmic accuracies. After resummations, we determine the PDF in the region 𝑥 ∈[0.25,0.75]. The end-point regions are then parametrized, combined with the LaMET prediction at moderate 𝑥, and fitted to the short-distance matrix elements in coordinate space.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Coexistence of ferromagnetism and antiferromagnetic dimers in topological insulators

The addition of magnetic impurities in topological insulators (TIs) can drive ferromagnetic order that leads to quantum anomalous Hall transport well below the Curie temperature. The fragility of the quantized regime has been ascribed to the random nature of the magnetic moment distribution. Here, we refine this hypothesis by using inelastic neutron scattering and density-functional theory calculations to show that two antagonistic components define the magnetism in Mn-substituted SnTe, thereby limiting the effectiveness of dilute magnetic TIs. One component is strongly bound antiferromagnetic dimers that compete with ferromagnetic order. In conclusion, the other component consists of undimerized moments where ferromagnetic order develops via long-range interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗