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 289 records · Page 16

On the Use of the Log-Normal Particle Size Distribution to Characterize Global Rain

Although most parameterizations of the drop size distributions (DSD) use the gamma function, there are several advantages to the log-normal form, particularly if we want to characterize the large scale space-time variability of the DSD and rain rate. The advantages of the distribution are twofold: the logarithm of any moment can be expressed as a linear combination of the individual parameters of the distribution; the parameters of the distribution are approximately normally distributed. Since all radar and rainfall-related parameters can be written approximately as a moment of the DSD, the first property allows us to express the logarithm of any radar/rainfall variable as a linear combination of the individual DSD parameters. Another consequence is that any power law relationship between rain rate, reflectivity factor, specific attenuation or water content can be expressed in terms of the covariance matrix of the DSD parameters. The joint-normal property of the DSD parameters has applications to the description of the space-time variation of rainfall in the sense that any radar-rainfall quantity can be specified by the covariance matrix associated with the DSD parameters at two arbitrary space-time points. As such, the parameterization provides a means by which we can use the spaceborne radar-derived DSD parameters to specify in part the covariance matrices globally. However, since satellite observations have coarse temporal sampling, the specification of the temporal covariance must be derived from ancillary measurements and models. Work is presently underway to determine whether the use of instantaneous rain rate data from the TRMM Precipitation Radar can provide good estimates of the spatial correlation in rain rate from data collected in 5(sup 0)x 5(sup 0) x 1 month space-time boxes. To characterize the temporal characteristics of the DSD parameters, disdrometer data are being used from the Wallops Flight Facility site where as many as 4 disdrometers have been used to acquire data over a 2 km path. These data should help quantify the temporal form of the covariance matrix at this site.

Meneghini, Robert↗

Techniques for estimating Space Station aerodynamic characteristics

A method was devised and calculations were performed to determine the effects of reflected molecules on the aerodynamic force and moment coefficients for a body in free molecule flow. A procedure was developed for determining the velocity and temperature distributions of molecules reflected from a surface of arbitrary momentum and energy accommodation. A system of equations, based on momentum and energy balances for the surface, incident, and reflected molecules, was solved by a numerical optimization technique. The minimization of a 'cost' function, developed from the set of equations, resulted in the determination of the defining properties of the flow reflected from the arbitrary surface. The properties used to define both the incident and reflected flows were: average temperature of the molecules in the flow, angle of the flow with respect to a vector normal to the surface, and the molecular speed ratio. The properties of the reflected flow were used to calculate the contribution of multiply reflected molecules to the force and moments on a test body in the flow. The test configuration consisted of two flat plates joined along one edge at a right angle to each other. When force and moment coefficients of this 90 deg concave wedge were compared to results that did not include multiple reflections, it was found that multiple reflections could nearly double lift and drag coefficients, with nearly a 50 percent increase in pitching moment for cases with specular or nearly specular accommodation. The cases of diffuse or nearly diffuse accommodation often had minor reductions in axial and normal forces when multiple reflections were included. There were several cases of intermediate accommodation where the addition of multiple reflection effects more than tripled the lift coefficient over the convex technique.

Thomas, Richard E.↗

Nonlocal chiral contributions to generalized parton distributions of the proton at nonzero skewness

We compute the one-loop contributions to spin-averaged generalized parton distributions (GPDs) in the proton from pseudoscalar mesons with intermediate octet and decuplet baryon states at nonzero skewness. Our framework is based on nonlocal covariant chiral effective theory, with ultraviolet divergences regularized by introducing a relativistic regulator derived consistently from the nonlocal Lagrangian. Using the splitting functions calculated from the nonlocal Lagrangian, we find the nonzero skewness GPDs from meson loops by convoluting with the phenomenological pion GPD and the generalized distribution amplitude, and verify that these satisfy the correct polynomiality properties. We also compute the lowest two moments of GPDs to quantify the meson loop effects on the Dirac, Pauli, and gravitational form factors of the proton.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A model for the three-dimensional spacecraft control laboratory experiment

A model for the three-dimensional Spacecraft Control Laboratory Experiment (SCOLE) is developed. The objective behind this method of modelling is to utilize the basic partial differential equations of motion for this distributed parameter system and not to use the modal expansion in developing the model. The final model obtained is in terms of a transfer function matrix which relates the flexible mast parameters like displacement, slope, shear stress, etc. to external forces and moments.

Kakad, Yogendra P.↗

Multipole correction of atomic monopole models of molecular charge distribution. I. Peptides

The defects in atomic monopole models of molecular charge distribution have been analyzed for several model-blocked peptides and compared with accurate quantum chemical values. The results indicate that the angular characteristics of the molecular electrostatic potential around functional groups capable of forming hydrogen bonds can be considerably distorted within various models relying upon isotropic atomic charges only. It is shown that these defects can be corrected by augmenting the atomic point charge models by cumulative atomic multipole moments (CAMMs). Alternatively, sets of off-center atomic point charges could be automatically derived from respective multipoles, providing approximately equivalent corrections. For the first time, correlated atomic multipoles have been calculated for N-acetyl, N'-methylamide-blocked derivatives of glycine, alanine, cysteine, threonine, leucine, lysine, and serine using the MP2 method. The role of the correlation effects in the peptide molecular charge distribution are discussed.

NASA Discipline Exobiology↗

Weibel instability beyond bi-Maxwellian anisotropy

The shape of the anisotropic velocity distribution function, beyond the realm of strict Maxwellians can play a significant role in determining the evolution of the Weibel instability dictating the dynamics of self-generated magnetic fields. For non-Maxwellian distribution functions, we show that the direction of the maximum growth rate wave vector changes with shape. In this work, we investigate different laser-plasma interaction model distributions which show that their Weibel generated magnetic fields may require closer scrutiny beyond the second moment (temperature) anisotropy ratio characterization.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Digital Twin + AI: Control Room of the Future [Slides]

The control room functions as the central brain of the grid, essential for balancing supply and demand and ensuring moment-to-moment grid reliability. Like the human brain, which processes sensory data to make decisions, control room operators analyze operational data from power generation, transmission, and distribution to make informed decisions. Currently, decision-making primarily rests with operators due to hardware and software limitations. However, with technological advancements, Digital Twins and AI are becoming high interest points in the control room's decision-making pilot programs. NREL is developing a comprehensive decision-making platform that integrates Digital Twins, AI, and advanced visualization techniques. As this integration progresses, the role of Digital Twins will evolve from conducting automated simulations to serving as a Trustworthy AI enabler, offering verification and validation of AI-generated response for power systems or providing physics-aware synthetic data of AI pre-training.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Radar Cross Section of Orbital Debris Objects

This discussion is concerned with the radar-data analysis and usage involved in the building of model orbital debris (OD) populations in the near-Earth environment, focusing on radar cross section (RCS). While varying with radar wavelength, physical dimension, material composition, overall shape and structure, the RCS of an irregular object is also strongly dependent on its spatial orientation. The historical records of observed RCSs for cataloged OD objects in the Space Surveillance Network are usually distributed over an RCS range, forming respective characteristic patterns. The National Aeronautics and Space Administration (NASA) Size Estimation Model provides an empirical probability-density function of RCS as a function of “effective diameter” (or “characteristic length”), which makes it feasible to predict possible RCS distributions for a given model OD population and to link data with model from a statistical perspective. The discussion also includes application of the widely used method of moments (MoM) and the Generalized Multi-particle Mie-solution (GMM) in the prediction of the RCS of arbitrarily shaped objects. Theoretical calculation results for an aluminum cube are compared with corresponding experimental measurements.

Xu, Yu-Lin↗

A Variable Eddington Factor Model for Thermal Radiative Transfer with Closure Based on Data-Driven Shape Function

Here, a new variable Eddington factor (VEF) model is presented for nonlinear problems of thermal radiative transfer (TRT). The VEF model is data-driven and acts on known (a-priori) radiation-diffusion solutions for material temperatures in the TRT problem. A linear auxiliary problem is constructed for the radiative transfer equation (RTE) whose emission source and opacities are evaluated at these known material temperatures. The solution to this RTE approximates the specific intensity distribution in phase-space and time. It is applied as a shape function to define the Eddington tensor for the presented VEF model. The shape function computed via the auxiliary RTE problem will capture some degree of transport effects within the TRT problem. The VEF moment equations closed with this approximate Eddington tensor will thus carry with them these captured transport effects. In this study, the temperature data comes from multigroup P 1 , P 1/3 , and flux-limited diffusion radiative transfer models. The proposed VEF model can be interpreted as a transport-corrected diffusion reduced-order model. Numerical results are presented on the Fleck-Cummings test problem which models a supersonic wavefront of radiation. The VEF model is shown to improve accuracy by 1–2 orders of magnitude compared to the considered radiation-diffusion model solutions to the TRT problem.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Scalar Dissipation Modeling for Passive and Active Scalars: a priori Study Using Direct Numerical Simulation

Transitional databases from Direct Numerical Simulation (DNS) of three-dimensional mixing layers for single-phase flows and two-phase flows with evaporation are analyzed and used to examine the typical hypothesis that the scalar dissipation Probability Distribution Function (PDF) may be modeled as a Gaussian. The databases encompass a single-component fuel and four multicomponent fuels, two initial Reynolds numbers (Re), two mass loadings for two-phase flows and two free-stream gas temperatures. Using the DNS calculated moments of the scalar-dissipation PDF, it is shown, consistent with existing experimental information on single-phase flows, that the Gaussian is a modest approximation of the DNS-extracted PDF, particularly poor in the range of the high scalar-dissipation values, which are significant for turbulent reaction rate modeling in non-premixed flows using flamelet models. With the same DNS calculated moments of the scalar-dissipation PDF and making a change of variables, a model of this PDF is proposed in the form of the (beta)-PDF which is shown to approximate much better the DNS-extracted PDF, particularly in the regime of the high scalar-dissipation values. Several types of statistical measures are calculated over the ensemble of the fourteen databases. For each statistical measure, the proposed (beta)-PDF model is shown to be much superior to the Gaussian in approximating the DNS-extracted PDF. Additionally, the agreement between the DNS-extracted PDF and the (beta)-PDF even improves when the comparison is performed for higher initial Re layers, whereas the comparison with the Gaussian is independent of the initial Re values. For two-phase flows, the comparison between the DNS-extracted PDF and the (beta)-PDF also improves with increasing free-stream gas temperature and mass loading. The higher fidelity approximation of the DNS-extracted PDF by the (beta)-PDF with increasing Re, gas temperature and mass loading bodes well for turbulent reaction rate modeling.

scalar dissipation modeling↗

Moments of dipole oscillator-strength distribution for the helium sequence

The moments S(mu) for mu at least -6 but no more than 2 and L(mu) for mu = 0, 1, and 2 are calculated for the helium sequence for atomic numbers (Z) up to 30 under a screened hydrogenic model. The model describes the atom by single-particle hydrogenic wave functions and treats the initial and the final state as characterized by two different effective charge parameters Zi and Zf, respectively. The differential oscillator strength of the screened hydrogenic model is asymptotically expanded. Assuming the value of 287.6 for the coefficient of the term epsilon to the -7/2 for helium atoms, the parameter Zf is determined for the helium sequence.

Khan, F.↗

A deterministic verification strategy for electrostatic particle-in-cell algorithms in arbitrary spatial dimensions using the method of manufactured solutions

As simulations of kinetic plasmas continue to increase in scope and complexity, a rigorous and straightforward method for verifying particle-in-cell (PIC) implementations is necessary to ensure their correctness. Here, in this paper, we present a deterministic method for the rigorous verification of multidimensional, multispecies, electrostatic particle-in-cell codes based on the method of manufactured solutions. Specifically, we prove that rigorous verification is possible through the exclusive examination of errors of grid quantities (i.e., moments and/or fields), allowing for a very light-weight and non-intrusive implementation in existing PIC codes. This is a marked improvement over earlier PIC verification approaches (only demonstrated with one species in 1D-1V), which rely on the comparison of cumulative distribution functions, and are invasive on the code base, introduce additional stochastic noise, are significantly more computationally expensive, and lack rigorous convergence properties. Interestingly, we show that different grid quantities feature different rates of convergence with the number of particles and mesh size, impacting the verification process. These theoretical results are confirmed numerically with a multi-species 2D-2V particle-in-cell code, including the ability of the method to detect order reduction due to an incorrect implementation.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Statistical comparison of galaxy formation models - The bispectrum

To discriminate between galaxy-formation models that give similar two-point correlation functions, higher moments, specifically the three-point function in transform space or bispectrum, are looked to for an ensemble of numerical simulations. Different behaviors of this statistic are found that can be used to distinguish between classes of models. Up to the uncertainty in identifying the galaxy distribution with the mass distribution in the models, the results seem to prefer the isothermal or the cold-particle models. No model, however, agrees in all detail with the observational bispectrum of the Shane-Wirtanen galaxy catalog.

Fry, J. N.↗

A New Method for Interfacing Unsuited Subjects to Overhead Suspension Partial Gravity Simulators

The purpose of performing unsuited testing as part of a reduced gravity extravehicular (EVA) suited human performance research program is to define baseline performance. These results are then coupled with suited test results to evaluate how the suit system affects human performance at reduced gravity. The primary drawback to this approach is that previous studies used notably different systems to interface suited and unsuited subjects to overhead-suspension, partial-gravity simulators. A spreader bar (SB) assembly previously used for unsuited tests allowed limited pitch and roll of the subject, whereas the gimbal for suited tests allowed more pitch and roll, although the mass distribution led to large moments of inertia in the yaw axis. It is hypothesized that use of the same methods for offload of both unsuited and suited subjects is needed to make meaningful comparisons. A new gimbal (GIM) was designed with the idea that it could function with both suited and unsuited subjects. GIM was designed to minimize mass and moments of inertia and to be adjustable to co-locate the 3 axes of rotation with the subject s center of gravity. OBJECTIVE: To evaluate human performance differences between SB and GIM. METHODS: Ten unsuited subjects were off-loaded to 1/6-g using both interfaces. Subjects completed tasks including overground and treadmill ambulation, picking up objects, shoveling, postural stability, range of motion testing, and recovery from the kneeling and prone positions. Metabolic, biomechanical, and/or subjective data were collected based on task. RESULTS: Initial analyses suggest that subjects completed all tasks with lower levels of compensation and a more terrestrial approach to movement when suspended via GIM. With SB, subjects were not able to fall or get into a prone position and had increased difficulty both retrieving objects off the floor and with overground ambulation, especially at gait initiation, because they were unable to bend their torso. GIM shows promise as a new method.

Norcross, Jason R.↗

A New Method for Interfacing Unsuited Subjects to Overhead Suspension Partial Gravity Simulators

The purpose of performing unsuited testing as part of a reduced gravity extravehicular (EVA) suited human performance research program is to define baseline performance. These results are then coupled with suited test results to evaluate how the suit system affects human performance at reduced gravity. The primary drawback to this approach is that previous studies used notably different systems to interface suited and unsuited subjects to overhead-suspension, partial-gravity simulators. A spreader bar (SB) assembly previously used for unsuited tests allowed limited pitch and roll of the subject, whereas the gimbal for suited tests allowed more pitch and roll, although the mass distribution led to large moments of inertia in the yaw axis. It is hypothesized that use of the same methods for offload of both unsuited and suited subjects is needed to make meaningful comparisons. A new gimbal (GIM) was designed with the idea that it could function with both suited and unsuited subjects. GIM was designed to minimize mass and moments of inertia and to be adjustable to co-locate the 3 axes of rotation with the subject s center of gravity. OBJECTIVE: To evaluate human performance differences between SB and GIM. METHODS: Ten unsuited subjects were offloaded to 1/6-g using both interfaces. Subjects completed tasks including overground and treadmill ambulation, picking up objects, shoveling, postural stability, range of motion testing, and recovery from the kneeling and prone positions. Metabolic, biomechanical, and/or subjective data were collected based on task. RESULTS: Initial analyses suggest that subjects completed all tasks with lower levels of compensation and a more terrestrial approach to movement when suspended via GIM. With SB, subjects were not able to fall or get into a prone position and had increased difficulty both retrieving objects off the floor and with overground ambulation, especially at gait initiation, because they were unable to bend their torso. GIM shows promise as a new method.

Norcross, Jason↗

Learning PDFs through interpretable latent representations in Mellin space

Representing the parton distribution functions (PDFs) of the proton and other hadrons through flexible, high-fidelity parametrizations has been a long-standing goal of particle physics phenomenology. This is particularly true since the chosen parametrization methodology can play an influential role in the ultimate PDF uncertainties as extracted in QCD global analyses; these, in turn, are often determinative of the reach of experiments at the LHC and other facilities to nonstandard physics, including at large 𝑥, where parametrization effects can be significant. In this study, we explore a series of encoder-decoder machine-learning (ML) models with various neural-network topologies as efficient means of reconstructing PDFs from meaningful information stored in an interpretable latent space. Given recent effort to pioneer synergies between QCD analyses and lattice-gauge calculations, we formulate a latent representation based on the behavior of PDFs in Mellin space, i.e., their integrated moments, and test the ability of various models to decode PDFs from this information faithfully. We introduce a numerical package, PDFdecoder, which implements several encoder-decoder models to reconstruct PDFs with high fidelity and use this end-to-end tool to explore how such neural-network-based models might connect PDF parametrizations to underlying properties like their Mellin moments. We additionally dissect patterns of learned correlations between encoded Mellin moments and reconstructed PDFs that suggest opportunities for further improvements to ML-based approaches to PDF parametrizations and uncertainty quantification.

Machine learning↗

Constitutive relations associated with the Mott-Smith distribution function

It is shown that the distribution function assumed by Mott-Smith determines a unique relation between heat flux, stress, and fluid velocity given by q = (3/2) tau u - i.e., it provides a constitutive relation for heat flux, and it also determines a simple expression for the ratio of third-order central moments. These expressions allow the equation of transfer for (C sub x) squared to be cast in a form that yields a nonlinear constitutive relation for stress. The results obtained from the Mott-Smith ansatz are compared with the theory of Baganoff and Nathenson (1970) and results from a numerical solution of the Boltzmann equation for shock-wave structure obtained by Hicks et al. (1972).

Nathenson, M.↗