Search NASA⌕ Search

SEARCH · Search NASA

Results for “MOMENT DISTRIBUTION”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Convex relaxation for Fokker–Planck equation

We propose an approach to directly estimate the moments or marginals for a high-dimensional equilibrium distribution in statistical mechanics by solving the high-dimensional Fokker–Planck equation in terms of low-order cluster moments or marginals. With this approach, we bypass the exponential complexity of estimating the full high-dimensional distribution and directly solve the simplified partial differential equations for low-order moments/marginals. Moreover, the proposed moment/marginal relaxation is fully convex and can be solved via off-the-shelf solvers. We further propose a time-dependent version of the convex programs to study non-equilibrium dynamics. In a specific setting, we show the proposed method can recover a mean-field-type equilibrium density. Numerical results are provided to demonstrate the performance of the proposed algorithm for high-dimensional systems.

Chen, Yian↗
DOE OSTI2021

New method for detecting fast neutrino flavor conversions in core-collapse supernova models with two-moment neutrino transport

Fast-pairwise neutrino oscillations potentially affect many aspects of core-collapse supernovae (CCSNe): the explosion mechanism, neutrino signals, and nucleosynthesis in the ejecta. This particular mode of collective neutrino oscillations has a deep connection to the angular structure of neutrinos in momentum space; for instance, the appearance of electron-neutrino lepton number (ELN) angular crossings in momentum space is a good indicator of the occurrence of flavor conversions. However, many multidimensional CCSN simulations are carried out with approximate neutrino transport (such as two-moment methods), which limits access to the angular distributions of neutrinos, i.e., inhibits ELNcrossing searches. In this paper, we develop a new method of searching for ELN crossing in these CCSN simulations. The required data is the zeroth and first angular moments of neutrinos and the matter profile, all of which are available in CCSN models with the two-moment method. One of the novelties of our new method is the use of a ray-tracing neutrino transport to determine ELNs in the direction of the stellar center. It is designed to compensate for shortcomings of the crossing searches with only two angular moments. We assess the capability of the method by carrying out a detailed comparison with results of full Boltzmann neutrino transport in 1D and 2D CCSN models. In this work, we find that the ray-tracing neutrino transport improves the accuracy of crossing searches; indeed, the appearance/disappearance of the crossings is accurately detected even in the region of forward-peaked angular distributions. The new method is computationally cheap and has the benefit of efficient parallelization; hence, it will be useful for ELN-crossing searches in any CCSN models that employ two-moment neutrino transport.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Population balance modeling of polyurethane foam formation with pressure‐dependent growth kernel

Abstract Polyurethane foams are widely used materials often chosen for their useful characteristics such as low thermal conductivity, ease of application, and high strength‐to‐weight ratios. Computational models are needed to predict the dynamics of the flow and expansion, and the resulting material properties, to improve manufacturing processes. In this paper, a model for PMDI, a water‐blown polyurethane foam, is presented. By extending a kinetics‐based approach by adding bubble‐scale information via a population balance equation (PBE) using the quadrature method of moments, we can track bubble size distributions during foaming. We present results from a three‐dimensional computational fluid dynamics model using arbitrary Lagrangian–Eulerian interface tracking implemented in finite element software. The model compares favorably with experimental data, including dynamics, bubble distributions measured by both camera and diffusion wave spectroscopy, and post‐test bubble size from scanning electron microscopy and density measurements from x‐ray computed tomography.

Ortiz, Weston↗
DOE OSTI2022

Kinetic Landau-fluid closures of non-Maxwellian distributions

New kinetic Landau-fluid closures, based on the cutoff Maxwellian distribution, are derived. A special static case is considered (the frequency ω=0). In the strongly collisional regime, our model reduces to Braginskii's heat flux model, and the transport is local. In the weak collisional regime, our model indicates that the heat flux is non-local and recovers the Hammett–Perkins model while the value of the cutoff velocity approaches to infinity. We compare the thermal transport coefficient χ of Maxwellian, cutoff Maxwellian and super-Gaussian distribution. The results show that the reduction of the high-speed tail particles leads to the corresponding reduction of the thermal transport coefficient χ across the entire range of collisionality, more reduction of the free streaming transport toward the weak collisional regime. In the collisionless limit, χ approaches to zero for the cutoff Maxwellian and the super-Gaussian distribution but remains finite for Maxwellian distribution. χ is complex if the cutoff Maxwellian distribution is asymmetric. The Im(χ) approaches to different convergent values in both collisionless and strongly collisional limit, respectively. It yields an additional streaming heat flux in comparison with the symmetric cutoff Maxwellian distribution. Furthermore, due to the asymmetric distribution, there is a background heat flux q0 though there is no perturbation. Furthermore, the derived Landau-fluid closures are general for fluid moment models, and applicable for the cutoff Maxwellian distribution in an open magnetic field line region, such as the scape-off-layer of Tokamak plasmas, in the thermal quench plasmas during a tokamak disruption, and the super-Gaussian electron distribution function due to inverse bremsstrahlung heating in laser-plasma studies.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

A physics-constrained deep learning treatment of runaway electron dynamics

An adjoint formulation leveraging a physics-informed neural network (PINN) is employed to advance the density moment of a runaway electron (RE) distribution forward in time. A distinguishing feature of this approach is that once the adjoint problem is solved, its solution can be used to project the RE density forward in time for an arbitrary initial momentum space distribution of REs. Furthermore, by employing a PINN, a parametric solution to the adjoint problem can be learned. Thus, once trained, this adjoint-deep learning framework is able to efficiently project the RE density forward in time across various plasma conditions while still including a fully kinetic description of RE dynamics. As an example application, the temporal evolution of the density of primary electrons is studied, with particular emphasis on evaluating the decay of a RE population when below threshold. Predictions from the adjoint-deep learning framework are found to be in good agreement with a traditional relativistic electron Fokker–Planck solver, for several distinct initial conditions, and across an array of physics parameters. Once trained, the PINN thus provides a means of generating RE density time histories with exceptionally low online execution time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Equation of state within the EPOS3 model

Transitions between different states of matter and their thermodynamic properties are described by the Equation of State (EoS). A universal representation of the EoS of quantum chromodynamics for the wide range of phase diagrams has yet to be determined. The expectation of the systems to undergo various types of transitions depending on the temperature (T), the chemical potential (μ B ), and other thermodynamic features make solving that puzzle challenging. Furthermore, it needs to be apparent which experimentally measurable observables could provide helpful information for determining EoS. The application of different EoS for hydrodynamical evolution was introduced in the EPOS3 generator, which allows one to study its changing effect on the experimental observables. The family of EoS proposed by the BEST Collaboration was implemented. The critical point location and the strength of criticality variations were investigated with particle yield, transverse momentum spectra, flow, and moments of the net-proton distributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2023

Dose-Weighted Radiographic Source-Spot Size

The spatial resolution of a flash radiograph is determined by the time-integrated size of the radiation source. Since the radiation pulse includes the beam head and tail, which may be substantially larger than the beam body, the radiation source is larger than the focused flattop of the electron beam pulse. Moreover, the radiation dose rate is weighted by both beam current and electron energy, and it is the radiation dose that forms the radiograph image. Therefore, in order to estimate radiographic resolution from computer simulations of a focused electron beam, one should use a dose weighted moment analysis of the radiation source distribution. This approach is used to show that, in the absence of beam effects other than space-charge and emittance, radiographic resolution for Scorpius exceeds requirements, regardless of reasonable timing options for current and gap-voltage pulses.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Examination of New Theory for Neutron Multiplicity Counting of Non-Point-Like Sources of Special Nuclear Material

The purpose of a nondestructive assay is to accurately verify the declared mass of special nuclear material (SNM) samples in a limited amount of time. One measurement modality is neutron multiplicity counting (NMC), which relates time-correlated neutron detection rates to the mass of SNM present. Traditional theory assumes point-like sources, which can be ill-posed for kilogram-quantity, bulk samples. Recent theory was developed to improve NMC accuracy for non-point-like samples. This work preliminarily examines the theory with measured data. The OSCAR prototype (a 3-by-4 array of 5.08-cm-thick, 5.08-cm-diameter trans-stilbene organic scintillators) measures configurations of 29.41-49.00 kg of highly enriched uranium (93wt% 235 U). Adjusted factorial moments for the emitted neutron multiplicity distributions are used to account for the shape and type of SNM being measured.

252Cf↗
DOE OSTI2026

First simultaneous global QCD analysis of kaon and pion parton distributions with lattice QCD constraints

We perform the first simultaneous global QCD analysis of pion and kaon parton distribution functions (PDFs) constrained by pion- and kaon-induced Drell-Yan (DY) and leading neutron electroproduction data together with lattice QCD data on pion and kaon PDF moments. The analysis indicates a softer valence u ¯ distribution in the K − than in the π − and a significantly more peaked valence s -quark density in K − compared with the u ¯ . The effective exponent governing the high- x behavior of the PDF is found to be larger for u ¯ in the kaon, β u ¯ K − = 1.6 ( 2 ) , than in the pion, β u ¯ π − = 1.16 ( 4 ) , in the range 0.7 ≤ x ≤ 0.95 . From the gluon momentum fractions, we find the pion’s gluon content accounts for ≈ 1 / 3 of the mass budget of the pion at μ = 2 GeV , but only ≈ 1 / 4 for the kaon.

Barry, P. C. [Argonne National Laboratory] (ORCID:↗
DOE OSTI2020

Transverse-energy-energy correlations in deep inelastic scattering

Event shape observables have been widely used for precision QCD studies at various lepton and hadron colliders. We present the most accurate calculation of the transverse-energy-energy correlation event shape variable in deep-inelastic scattering. In the framework of soft-collinear effective theory the cross section is factorized as the convolution of the hard function, beam function, jet function and soft function in the back-to-back limit. A close connection to TMD factorization is established, as the beam function when combined with part of the soft function is identical to the conventional TMD parton distribution function, and the jet function is the second moment of the TMD fragmentation function matching coefficient. We validate our framework by comparing the obtained LO and NLO leading singular distributions to the full QCD calculations in the back-to-back limit. We report the resummed transverse-energy-energy correlation distributions up to N 3 LL accuracy matched with the NLO cross section for the production of a lepton and two jets. Our work provides a new way to precisely study TMD physics at the future Electron-Ion Collider.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

The heliospheric ambipolar potential inferred from sunward-propagating halo electrons

ABSTRACT We provide evidence that the sunward-propagating half of the solar wind electron halo distribution evolves without scattering in the inner heliosphere. We assume the particles conserve their total energy and magnetic moment, and perform a ‘Liouville mapping’ on electron pitch angle distributions measured by the Parker Solar Probe SPAN-E instrument. Namely, we show that the distributions are consistent with Liouville’s theorem if an appropriate interplanetary potential is chosen. This potential, an outcome of our fitting method, is compared against the radial profiles of proton bulk flow energy. We find that the inferred potential is responsible for nearly 100 per cent of the proton acceleration in the solar wind at heliocentric distances 0.18-0.79 AU. These observations combine to form a coherent physical picture: the same interplanetary potential accounts for the acceleration of the solar wind protons as well as the evolution of the electron halo. In this picture the halo is formed from a sunward-propagating population that originates somewhere in the outer heliosphere by a yet-unknown mechanism.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

A Bayesian Approach for Statistical–Physical Bulk Parameterization of Rain Microphysics. Part I: Scheme Description

A novel framework is proposed for the bulk parameterization of rain microphysics: the Bayesian Observationally Constrained Statistical–Physical Scheme (BOSS). It is designed to facilitate direct constraint by observations using Bayesian inference. BOSS combines existing process-level microphysical knowledge with flexible process rate formulations and parameters constrained by observations within a Bayesian framework. Furthermore, using a raindrop size distribution (DSD) normalization method that relates DSD moments to one another via generalized power series, generalized multivariate power expressions are derived for the microphysical process rates as functions of a set of prognostic DSD moments. The scheme is flexible and can utilize any number and combination of prognostic moments and any number of terms in the process rate formulations. This means that both uncertainty in parameter values and structural uncertainty associated with the process rate formulations can be investigated systematically, which is not possible using traditional schemes. In this paper, BOSS is compared to two- and three-moment versions of a traditional bulk rain microphysics scheme (denoted as MORR). It is demonstrated that some process formulations in MORR are analytically equivalent to the generalized power expressions in BOSS using one or two terms, while others are not. BOSS is able to replicate the behavior of MORR in idealized one-dimensional rainshaft tests, but with a much more flexible and systematic design. Part II of this study describes the application of BOSS to derive rain microphysical process rates and posterior parameter distributions in Bayesian experiments using Markov chain Monte Carlo sampling constrained by synthetic observations.

58 GEOSCIENCES↗
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 OSTI2020

Distributionally robust facility location problem under decision-dependent stochastic demand

While the traditional facility location problem considers exogenous demand, in some applications, locations of facilities could affect the willingness of customers to use certain types of services, e.g., carsharing, and therefore they also affect realizations of random demand. Moreover, a decision maker may not know the exact distribution of such endogenous demand and how it is affected by location choices. In this paper, we consider a distributionally robust facility location problem, in which we interpret the moments of stochastic demand as functions of facility-location decisions. We reformulate a two-stage decision-dependent distributionally robust optimization model as a monolithic formulation, and then derive exact mixed-integer linear programming reformulation as well as valid inequalities when the means and variances of demand are piecewise linear functions of location solutions. We conduct extensive computational studies, in which we compare our model with a decision-dependent deterministic model, as well as stochastic programming and distributionally robust models without the decision-dependent assumption. Here, the results show superior performance of our approach with remarkable improvement in profit and quality of service under various settings, in addition to computational speed-ups given by formulation enhancements. These results draw attention to the need of considering the impact of location decisions on customer demand within this strategic-level planning problem.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Magnetic oxygen in transition metal oxides: A case study of Ba 2 CoO 4

Transition metal oxides (TMOs) exhibit exotic magnetic properties in both naturally formed and artificially structured materials, often difficult to understand in conventional wisdom. Magnetic insulator Ba 2 CoO 4 has mystified the community, because the CoO 4 tetrahedron is completely isolated with the nearest Co atoms far apart (~5 Å), making it impossible to account for long-range magnetic ordering seen experimentally using only Co. By theoretically investigating magnetism and relating our findings to experimental observations in bulk Ba 2 CoO 4 , we illustrate for the first time that the magnetic moment on oxygen atoms are the origin of the unexpected long-range magnetic ordering and low magnetic dimensionality. We find that the magnetic moment is not only localized on Co atoms, as assumed in all conventional data analysis, but also distributed on its tetrahedrally-coordinated O atoms. The total magnetic moment of the CoO 4 building block is 4.63μ B with the magnetic moment on Co being only 3.08μ B . Therefore, the magnetic building block is CoO 4 not Co. Furthermore, our first principles calculations are capable of explaining the origin of the unique magnetic response, including the presence of long-range magnetic ordering with two-dimensional character, and a one-dimensional magnetoelastic behavior. Having oxygen contribute to the magnetic moment will undoubtedly be identified as a universal property of magnetic TMOs, which will require a fresh look at conventional models of magnetism in TMOs.

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

Small-𝑥 gluon GPD constrained from deeply virtual 𝐽/𝜓 production and gluon PDF through universal-moment parametrization

We phenomenologically constrain the small-𝑥 and small-𝜉 gluon generalized parton distributions (GPDs) with the deeply virtual 𝐽/𝜓 production (DV⁢𝐽/𝜓⁢P) in the framework of GPDs through universal moment parameterization (GUMP). We use a hybrid cross-section formula combining collinear factorization to the next-to-leading order (NLO) accuracy of the strong coupling 𝛼 𝑠 , with corrections from nonrelativistic QCD to account for the power corrections due to the heavy 𝐽/𝜓 mass. We reach reasonable fit to the measured differential cross sections of DV⁢𝐽/𝜓⁢P by H1 at the Hadron-Electron Ring Accelerator as well as forward gluon parton distribution functions (PDFs) from JAM22 global analysis. We find that both NLO and nonrelativistic corrections are significant for heavy vector meson productions. Of course, the gluon GPD we obtain still contains considerable freedom in need of inputs from other constraints, particularly in the distribution-amplitude-like region.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2025

Hot or Not? An Evaluation of Methods for Identifying Hot Moments of Nitrous Oxide Emissions From Soils

Abstract Effectively quantifying hot moments of nitrous oxide (N 2 O) emissions from agricultural soils is critical for managing this potent greenhouse gas. However, we are challenged by a lack of standard approaches for identifying hot moments, including (a) determining thresholds above which emissions are considered hot moments, and (b) considering seasonal variation in the magnitude and frequency distribution of net N 2 O fluxes. We used one year of hourly N 2 O flux measurements from 16 autochambers that varied in flux magnitude and frequency distribution in a conventionally tilled maize field in central Illinois, USA, to compare three approaches to identify hot moment thresholds: standard deviations (SD) above the mean, 1.5x the interquartile range (IQR), and isolation forest (IF) identification of anomalous values. We also compared these approaches on seasonally subdivided data (early, late, and non‐growing seasons) versus the whole year. Our analyses revealed that 1.5x IQR method best identified N 2 O hot moments. In contrast, using 2 or 4 SD both yielded hot moment threshold values too high, and IF yielded threshold values too low, leading to missed N 2 O hot moments or low net N 2 O fluxes mischaracterized as hot moments, respectively. Furthermore, seasonally subdividing the data set not only facilitated identification of smaller hot moments in the late‐ and non‐growing seasons when N 2 O hot moments were generally smaller but it also increased hot moment threshold values in the early growing season when N 2 O hot moments were larger. Consequently, of the methods evaluated here, we recommend using the 1.5x IQR method on whole year data sets to identify N 2 O hot moments.

Stuchiner, Emily R. [Institute for Sustainability,↗
DOE OSTI2026

Long-Range Fit: A Software Package for the Representation and Study of Long-Range Molecular Interactions

Describing intermolecular forces is fundamental to modeling and predicting the behavior of molecular systems. In particular, long-range molecular interactions─with electrostatic, induction, and dispersion as the main components─play a critical role, especially for low-temperature and low-density regimes. Long-range interactions are often described through perturbation theory, representing the electronic charge distribution via a multipolar series of the moments and polarizability tensors corresponding to each molecule. However, while the theory is well established, obtaining the resulting analytical expressions (and their practical implementation) constitutes a highly complex and system-dependent task. To address this challenge, we developed long-range-fit (LRF), an interactive and user-friendly software package designed to automate the generation and fitting of long-range interaction terms for arbitrary molecules in nondegenerate (ground or excited) electronic states. We have derived and implemented all terms up to 15th order, without approximations, via a spherical tensor representation, with symmetry adaptation to all molecular point-group symmetries. The resulting potential energy surface is compatible with most representations of the close interaction region.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗