Search NASA⌕ Search

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.

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

At least 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Plasma thermal transport with a generalized 8-moment distribution function

Moment equations that model plasma transport require an ansatz distribution function to close the system of equations. The resulting transport is sensitive to the specific closure used, and several options have been proposed in the literature. Two different 8-moment distribution functions can be generalized to form a single-parameter family of distribution functions. The transport coefficients resulting from this generalized distribution function can be expressed in terms of this free parameter. This provides the flexibility of matching the 8-moment model to some validating result at a given magnetization value, such as Braginskii’s transport, or the more recent results of Davies et al. [Physics of Plasma, 28, 012305 (2021)]. Here, this process can be thought of as solving for the 8-moment distribution function that matches the value of a transport coefficient given by a Chapman-Enskog expansion, while retaining the improved physical properties, such as finite propagation speeds and time dependence which belong to the hyperbolic moment models. Since the presented generalized distribution function only has a single free parameter, only a single transport coefficient can be matched at a time. However, this generalization process may be extended to provide multiple free parameters. The focus of this Brief Communication is on the dramatically improved thermal conductivity of the proposed model compared to the two base moment models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Moments of parton distribution functions of the pion from lattice QCD using gradient flow

We present a nonperturbative determination of the pion valence parton distribution function (PDF) moment ratios ⟨𝑥 𝑛−1 ⟩/⟨𝑥⟩ up to 𝑛 = 6, using the gradient flow in lattice quantum chromodynamics (QCD). As a testing ground, we employ SU(3) isosymmetric gauge configurations generated by the OpenLat initiative with a pseudoscalar mass of 𝑚 𝜋 ≃ 411 MeV. Our analysis uses four lattice spacings and a nonperturbatively improved action, enabling full control over the continuum extrapolation, and the limit of vanishing flow time, 𝑡 →0. The flowed ratios exhibit O(𝑎 2 ) scaling across the ensembles, and the continuum-extrapolated results, matched to the $\overline{MS}$ scheme at 𝜇 = 2 GeV using next-to-next-to-leading order matching coefficients, show only mild residual flow-time dependence. The resulting ratios, computed with a relatively small number of configurations, are consistent with phenomenological expectations for the pion’s valence distribution, with statistical uncertainties that are competitive with modern global fits. These findings demonstrate that the gradient flow provides an efficient and systematically improvable method to access partonic quantities from first principles. Future extensions of this work will target lighter pion masses toward the physical point, and applications to nucleon structure such as the proton PDFs and the gluon and sea-quark distributions.

lattice QCD↗
NASA NTRS2009

Spherical Harmonic Analysis of Particle Velocity Distribution Function: Comparison of Moments and Anisotropies using Cluster Data

This paper presents a spherical harmonic analysis of the plasma velocity distribution function using high-angular, energy, and time resolution Cluster data obtained from the PEACE spectrometer instrument to demonstrate how this analysis models the particle distribution function and its moments and anisotropies. The results show that spherical harmonic analysis produced a robust physical representation model of the velocity distribution function, resolving the main features of the measured distributions. From the spherical harmonic analysis, a minimum set of nine spectral coefficients was obtained from which the moment (up to the heat flux), anisotropy, and asymmetry calculations of the velocity distribution function were obtained. The spherical harmonic method provides a potentially effective "compression" technique that can be easily carried out onboard a spacecraft to determine the moments and anisotropies of the particle velocity distribution function for any species. These calculations were implemented using three different approaches, namely, the standard traditional integration, the spherical harmonic (SPH) spectral coefficients integration, and the singular value decomposition (SVD) on the spherical harmonic methods. A comparison among the various methods shows that both SPH and SVD approaches provide remarkable agreement with the standard moment integration method.

Gurgiolo, Chris↗
DOE OSTI2024

Moments of parton distribution functions of any order from lattice QCD

We describe a procedure to determine moments of parton distribution functions of any order in lattice quantum chromodynamics (QCD). The procedure is based on the gradient flow for fermion and gauge fields. The flowed matrix elements of twist-2 operators renormalize multiplicatively, and the matching with the physical matrix elements can be obtained using continuum symmetries and the irreducible representations of Euclidean 4-dimensional rotations. We calculate the matching coefficients at one-loop in perturbation theory for moments of any order in the flavor nonsinglet case. We also give specific examples of operators that could be used in lattice QCD computations. It turns out that it is possible to choose operators with identical Lorentz indices and still have a multiplicative matching. One can thus use twist-2 operators exclusively with temporal indices, thus substantially improving the signal-to-noise ratio in the computation of the hadronic matrix elements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1972

Approximation method for the kinetic Boltzmann equation

The further development of a method for approximating the Boltzmann equation is considered and a case of pseudo-Maxwellian molecules is treated in detail. A method of approximating the collision frequency is discussed along with a method for approximating the moments of the Boltzmann collision integral. Since the return collisions integral and the collision frequency are expressed through the distribution function moments, use of the proposed methods make it possible to reduce the Boltzmann equation to a series of approximating equations.

Shakhov, Y. M.↗
NASA NTRS1985

Analytical-numerical methods of calculations of energy and three-dimensional particle distributions in electromagnetic cascades

Analytical and numerical methods of calculation of the energy and three dimensional EPS characteristics are reported. The angular and lateral functions of electrons in EPS have been obtained by the Landau and small angle approximations A and B and compared with earlier data. A numerical method of solution of cascade equations for the EPS distribution function moments has been constructed. Considering the equilibrium rms angle as an example, errors appearing when approximating the elementary process cross sections by their asymptotic expressions are analyzed.

Ivanenko, I. P.↗
NASA NTRS1992

Reactive collision terms for fluid transport theory

A modified Bhatnagar-Gross-Krook (BGK) collision method is used for the derivation of reactive collision terms to allow for complete higher-moment approximations in fluid transport theory. The general reactive collision terms for the 8-, 10- and 13-moment approximation for binary reactions are derived. Then, reactive collision terms for specific chemistry are derived with the assumptions: (1) all reactants and products have 13-moment distribution functions; (2) simple, classical assumptions provide sufficient models for reaction dynamics; and (3) the reaction rate is energy-independent. The derived specific formulas reflect two extremes in reactive collision dynamics: direct and indirect reaction mechanisms. Finally, considerations for correct use of the reactive collision terms are discussed in the context of space plasma environments.

Eccles, J. V.↗
DOE OSTI2024

Gluon moment and parton distribution function of the pion from 𝑁 𝑓 = 2 +1 +1 lattice QCD

We present the first calculation of the pion gluon moment from lattice QCD in the continuum-physical limit. The calculation is done using clover fermions for the valence action with three pion masses, 220, 310 and 690 MeV, and three lattice spacings, 0.09, 0.12, and 0.15 fm, using ensembles generated by MILC Collaboration with 2+1+1 flavors of highly improved staggered quarks (HISQ). On the lattice, we nonperturbatively renormalize the gluon operator in RI/MOM scheme using the cluster-decomposition error reduction (CDER) technique to enhance the signal-to-noise ratio of the renormalization constant. We extrapolate the pion gluon moment to the continuum-physical limit and obtain ⟨𝑥⟩ 𝑔 = 0.394⁢(58) stat+NPR ⁢(39) mixing in the $\overline{MS}$ scheme at 2 GeV, with first error being the statistical error and uncertainties in nonperturbative renormalization, and the second being a systematic uncertainty estimating the effect of ignoring quark mixing. Our pion gluon momentum fraction has a central value lower than two recent single-ensemble lattice-QCD results near physical pion mass but is consistent with the recent global fits by JAM and xFitter and with most QCD-model estimates.

Astronomy & Astrophysics↗
NASA NTRS1971

Higher moment equations and the distribution function of the solar-wind plasma.

Study of the higher-moment equations for a collisionless fully ionized plasma. For a collisionless, heat-conducting plasma, the distribution function f is cylindrically symmetric about the direction of the magnetic field. It is shown that under a certain assumption the fourth moments of f can be expressed as simple functions of lower moments. Thus no higher-moment terms appear in the third-moment equations. The two third-moment equations, which are obtained in a simple form, join other lower-moment equations to form a closed set of moment equations. The new equations can be used to study the thermal anisotropy and the heat flux of the solar-wind proton. A special case of the cylindrically symmetric distribution function f is found to resemble the proton distribution function reconstructed from solar-wind data, and this resemblance justifies the assumption needed for decoupling the moment equations.

Whang, Y. C.↗
DOE OSTI2020

Moment preserving constrained resampling with applications to particle-in-cell methods

The Moment Preserving Constrained Resampling (MPCR) algorithm for particle resampling is introduced and applied to particle-in-cell (PIC) methods to increase simulation accuracy, reduce compute cost, and/or avoid numerical instabilities. The general algorithm partitions the system space into smaller subsets and resamples the distribution within each subset. Further, the algorithm is designed to conserve any number of particle and grid moments with a high degree of accuracy (i.e. machine accuracy). The effectiveness of MPCR is demonstrated with several numerical tests, including a use-case study in gyrokinetic fusion plasma simulations. Finally, the computational cost of MPCR is negligible compared to the cost of particle evolution in PIC methods, and the tests demonstrate that periodic particle resampling yields a significant improvement in the accuracy and stability of the results.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
NASA NTRS1980

Simulation of spacecraft charging environments by monoenergetic beams

To simulate the distributed spectra of plasmas in space, the differential current density spectrum of the plasma was divided into a number of energy bands and the beam energy and current were calculated for each band to provide a piecewise reproduction of the distributed spectrum. Beam energies and current densities were chosen to match the velocity moments of the plasma distribution function. The velocity moments are averages related to physical quantities such as particle density, flux, pressure, and energy flux, and have been used extensively to characterize the measured properties of plasmas in space. Combinations of one, two, and three beams were found to match two to six velocity moments of Maxwellian distributions. A computational model was used to compare the charging of a spacecraft by plasmas with distributed spectra and by monoenergetic beams.

Halverson, W.↗
DOE OSTI2024

Spatial imaging of proton via leading-twist nonskewed GPDs with basis light-front quantization

The internal image of the proton is unveiled by examining the generalized parton distributions (GPDs) at zero skewness, within the basis light-front quantized environment. Several distributions emerge when a quark is sampled with different currents depending upon the helicity arrangements of the active quark and the proton target. We investigate six of the eight leading-twist proton GPDs of the valence quarks, the helicity conserving distributions ( H , E , H ˜ ) and the helicity nonconserving ( H T , E T , H ˜ T ) distributions at skewness set to zero ( ζ = 0 ). We consider purely transverse momentum transfer and, hence, obtain results describe only the proton’s two-dimensional structure in the transverse plane. We present the Mellin moments of these distribution functions, where the first moment produces a form factor and the second Mellin moments help extract the information on partonic contributions to the hadronic angular momentum. We compare our results for the Mellin moments with those from lattice QCD and other approaches where available. We also present the GPDs in transverse position space. Published by the American Physical Society 2024

Astronomy & Astrophysics↗
DOE OSTI2024

Moments of nucleon unpolarized, polarized, and transversity parton distribution functions from lattice QCD at the physical point

The second Mellin moments ⟨ x ⟩ of the nucleon’s unpolarized, polarized, and transversity parton distribution functions are computed. Two lattice QCD ensembles at the physical pion mass are used: these were generated using a tree-level Symanzik-improved gauge action and 2 + 1 flavor tree-level improved Wilson Clover fermions coupling via 2-level HEX-smearing. The moments are extracted from forward matrix elements of local leading twist operators. We determine renomalization factors in RI-(S)MOM and match to MS ¯ at scale 2 GeV. Our findings show that operators that exhibit vanishing kinematics at zero momentum can have significantly reduced excited-state contamination. The resulting polarized moment is used to quantify the longitudinal contribution to the quark spin-orbit correlation. All our results agree within two sigma with previous lattice results. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2024

Seismic moment tensor classification using elliptical distribution functions on the hypersphere

Discrimination of underground explosions from naturally occurring earthquakes and other anthropogenic sources is one of the fundamental challenges of nuclear explosion monitoring. In an operational setting, the number of events that can be thoroughly investigated by analysts is limited by available resources. The capability to rapidly screen out events that can be robustly identified as not being explosions is, therefore, of great potential benefit. Nevertheless, possible mis-classification of explosions as earthquakes currently limits the use of screening methods for verification of test-ban treaties. Moment tensors provide a physics-based classification tool for the characterization of different seismic sources and have enabled the advent of new techniques for discriminating between earthquakes and explosions. Following normalization and projection of their six-degree vectors onto the hypersphere, existing screening approaches use spherically symmetric metrics to determine whether any new moment tensor may have been an explosion. Here, we show that populations of moment tensors for both earthquakes and explosions are anisotropically distributed on the hypersphere. Distributions possessing elliptical symmetry, such as the scaled von Mises–Fisher distribution, therefore provide a better description of these populations than the existing spherically symmetric models. We describe a method that uses these elliptical distributions in combination with a Bayesian classifier to achieve successful classification rates of 99 per cent for explosions and 98 per cent for earthquakes using existing catalogues of events from the western United States. The 1983 May 5 Crowdie underground nuclear test and 2018 July 20 DAG-1 deep-borehole chemical explosion are the only two explosions out of 140 that are incorrectly classified. Application of the method to the 2006–2017 nuclear tests in the Democratic People’s Republic of Korea yields 100 per cent identification rates and we provide a simple routine MTid for general usage. The approach provides a means to rapidly assess the likelihood of an event being an explosion and can be built into monitoring workflows that rely on simultaneously assessing multiple different discrimination metrics.

58 GEOSCIENCES↗