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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
NASA NTRS1993

A multi-dimensional kinetic-based upwind solver for the Euler equations

A multidimensional kinetic fluctuation-splitting scheme has been developed for the Euler equations. The scheme is based on an N-scheme discretization of the Boltzmann equation at the kinetic level for triangulated Cartesian meshes with a diagonal-adaptive strategy. The resulting Euler scheme is a cell-vertex fluctuation-splitting scheme where fluctuations in the conserved-variable vector Q are obtained as moments of the fluctuation in the Maxwellian velocity distribution function at the kinetic level. Encouraging preliminary results have been obtained for perfect gases on Cartesian meshes with first-order spatial accuracy. The present approach represents an improvement to the well-established dimensionally-split upwind schemes.

Eppard, W. M.↗
NASA NTRS1994

Resolution requirements for velocity gradients in turbulence

Since high resolution numerical simulations of turbulent flows, or experiments at high Reynolds numbers, represent a substantial investment in resources, the estimation of the minimum resolution required for the study of a given property has been the subject of continued interest. Different properties require, in general, different resolutions, and the present paper is dedicated to the requirements for the measurement of the probability distribution functions of the velocity gradients and, in particular, of their low order moments. The deviation of these quantities from the values corresponding to a Gaussian distribution was one of the first indications of the presence of Reynolds number-dependent intermittency and has been the object of recent interest as numerical simulations have become able to explore the distribution of gradients in the low Reynolds number range, while new experiments have extended the range to increasingly high Reynolds numbers. We will use progressive filtering of the results of numerical simulations of isotropic turbulence as a model for the effect of a sensor of finite size. The numerical issues will be addressed first to ensure that the simulations are fully resolved from the point of view of the velocity gradients. This will also give us an estimate for the numerical resolution required for the different quantities.

Jimenez, Javier↗
DOE OSTI2022

Synthesis of observed and simulated rain microphysics to inform a new Bayesian statistical framework for microphysical parameterization in climate models (Final Scientific Report)

This project aimed to investigate a new approach to bulk microphysics schemes. We identified the need to move beyond the fixed structural assumptions and approximation of existing schemes. For example, most schemes assume some functional form for the rain drop size distribution (e.g. a gamma or exponential distribution). Most schemes then evolve some number of statistical moments of that distribution via various processes, such as evaporation, sedimentation, collision-coalescence, and collisional breakup. These microphysical processes, in turn, are typically some fixed functional form derived from either some other (more detailed) model, or using some single-particle rates that are then integrated over the size distribution. While these process rate functions may have some free parameters to adjust, in other cases doing so is impossible (e.g. when the process rates are analytical or piecewise solutions to some target function). Our goal was to build and test a scheme that assumed no size distribution form, and used a series of power laws as the basis for the process rates. Moments of the size distribution would be predicted, but no underlying size distribution would be specified. Any number or choice of size distribution moments could be used, and any number of power laws could be employed to model the process rates. Thus, our approach is flexible, and can seamlessly scale across levels of complexity, as demanded by the data. Bayesian inference then would provide the formalism to estimate the model parameters and structure, allowing for robust uncertainty quantification. We proposed testing this framework in idealized simulations using bin microphysical schemes as a data source. We also proposed using real data to inform our microphysics scheme, and also that we would integrate our scheme into WRF.

58 GEOSCIENCES↗
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↗
NASA NTRS1986

On the Maxwellian distribution, symmetric form, and entropy conservation for the Euler equations

The Euler equations of gas dynamics have some very interesting properties in that the flux vector is a homogeneous function of the unknowns and the equations can be cast in symmetric hyperbolic form and satisfy the entropy conservation. The Euler equations are the moments of the Boltzmann equation of the kinetic theory of gases when the velocity distribution function is a Maxwellian. The present paper shows the relationship between the symmetrizability and the Maxwellian velocity distribution. The entropy conservation is in terms of the H-function, which is a slight modification of the H-function first introduced by Boltzmann in his famous H-theorem. In view of the H-theorem, it is suggested that the development of total H-diminishing (THD) numerical methods may be more profitable than the usual total variation diminishing (TVD) methods for obtaining wiggle-free solutions.

Deshpande, S. M.↗
NASA NTRS1991

Statistics of primordial density perturbations from discrete seed masses

The statistics of density perturbations for general distributions of seed masses with arbitrary matter accretion is examined. Formal expressions for the power spectrum, the N-point correlation functions, and the density distribution function are derived. These results are applied to the case of uncorrelated seed masses, and power spectra are derived for accretion of both hot and cold dark matter plus baryons. The reduced moments (cumulants) of the density distribution are computed and used to obtain a series expansion for the density distribution function. Analytic results are obtained for the density distribution function in the case of a distribution of seed masses with a spherical top-hat accretion pattern. More generally, the formalism makes it possible to give a complete characterization of the statistical properties of any random field generated from a discrete linear superposition of kernels. In particular, the results can be applied to density fields derived by smoothing a discrete set of points with a window function.

Scherrer, Robert J.↗
NASA NTRS1985

TRO-2D - A code for rational transonic aerodynamic optimization

Features and sample applications of the transonic rational optimization (TRO-2D) code are outlined. TRO-2D includes the airfoil analysis code FLO-36, the CONMIN optimization code and a rational approach to defining aero-function shapes for geometry modification. The program is part of an effort to develop an aerodynamically smart optimizer that will simplify and shorten the design process. The user has a selection of drag minimization and associated minimum lift, moment, and the pressure distribution, a choice among 14 resident aero-function shapes, and options on aerodynamic and geometric constraints. Design variables such as the angle of attack, leading edge radius and camber, shock strength and movement, supersonic pressure plateau control, etc., are discussed. The results of calculations of a reduced leading edge camber transonic airfoil and an airfoil with a natural laminar flow are provided, showing that only four design variables need be specified to obtain satisfactory results.

Davis, W. H., Jr.↗
NASA NTRS1935

Principal Effects of Axial Load on Moment-Distribution Analysis of Rigid Structures

This thesis presents the method of moment distribution modified to include the effect of axial load upon the bending moments. This modification makes it possible to analyze accurately complex structures, such as rigid fuselage trusses, that heretofore had to be analyzed by approximate formulas and empirical rules. The method is simple enough to be practicable even for complex structures, and it gives a means of analysis for continuous beams that is simpler than the extended three-moment equation now in common use. When the effect of axial load is included, it is found that the basic principles of moment distribution remain unchanged, the only difference being that the factors used, instead of being constants for a given member, become functions of the axial load. Formulas have been developed for these factors, and curves plotted so that their applications requires no more work than moment distribution without axial load. Simple problems have been included to illustrate the use of the curves.

James, Benjamin Wylie↗
DOE OSTI2024

Advanced Closures for Nonlinear Fluid Simulations of Plasmas (Final Report)

Advanced Closures for Nonlinear Fluid Simulations of Plasmas E. Held, Utah State University (Principal Investigator) J.-Y. Ji, Utah State University (Co-Investigator) Utah State University (USU) continued to develop and apply closures for hybrid fluid/kinetic simulations of magnetized fusion plasmas. The closures developed by USU's fusion theory and computation group use the drift ordering and are referred to as "continuum" and "general moment" approaches. They differ in their representation for the distribution functions which probabilistically describe the density of particles in 5 dimensional phase space. The two chosen velocity variables, needed in studies of kinetic physics for magnetized plasmas, are the speed and the component of a particle's velocity along the magnetic field. In magnetized plasmas, the charged particle trajectories are tied closely to the magnetic field lines. Both approaches account for binary interactions between particles using the Coulomb collision operator. Theoretical and computational development along both lines provided avenues for verification between the two approaches. This funding allowed USU to further develop physically accurate and numerically efficient models that advance our understanding of kinetic physics in experiments like the International Thermonuclear Experimental Reactor (ITER). The ITER device will likely prove the physical practicality of thermonuclear fusion by producing 500 megawatts of clean, abundant fusion power. This work helped to realize this goal by advancing our simulation capability using hybrid fluid/kinetic models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1986

Spatial analysis of storm depths from an Arizona raingage network

Eight years of summer rainstorm observations are analyzed by a dense network of 93 raingages operated by the U.S. Department of Agriculture, Agricultural Research Service, in the 150 km Walnut Gulch experimental catchment near Tucson, Arizona. Storms are defined by the total depths collected at each raingage during the noon-to-noon period for which there was depth recorded at any of the gages. For each of the resulting 428 storm days, the gage depths are interpolated onto a dense grid and the resulting random field analyzed to obtain moments, isohyetal plots, spatial correlation function, variance function, and the spatial distribution of storm depth.

Fennessey, N. M.↗
NASA NTRS1987

Application of spatial Poisson process models to air mass thunderstorm rainfall

Eight years of summer storm rainfall observations from 93 stations in and around the 154 sq km Walnut Gulch catchment of the Agricultural Research Service, U.S. Department of Agriculture, in Arizona are processed to yield the total station depths of 428 storms. Statistical analysis of these random fields yields the first two moments, the spatial correlation and variance functions, and the spatial distribution of total rainfall for each storm. The absolute and relative worth of three Poisson models are evaluated by comparing their prediction of the spatial distribution of storm rainfall with observations from the second half of the sample. The effect of interstorm parameter variation is examined.

Eagleson, P. S.↗
NASA NTRS2006

What does Reflection from Cloud Sides tell us about Vertical Distribution of Cloud Droplet Sizes?

Cloud development, the onset of precipitation and the effect of aerosol on clouds depend on the structure of the cloud profiles of droplet size and phase. Aircraft measurements of cloud profiles are limited in their temporal and spatial extent. Satellites were used to observe cloud tops not cloud profiles with vertical profiles of precipitation-sized droplets anticipated from CloudSat. The recently proposed CLAIM-3D satellite mission (cloud aerosol interaction mission in 3-D) suggests to measure profiles of cloud microphysical properties by retrieving them from the solar and infrared radiation reflected or emitted from cloud sides. Inversion of measurements from the cloud sides requires rigorous understanding of the 3-dimentional(3-D) properties of clouds. Here we discuss the reflected sunlight from the cloud sides and top at two wavelengths: one nonabsorbing to solar radiation (0.67 microns) and one with liquid water efficient absorption of solar radiation (2.1 microns). In contrast to the plane-parallel approximation, a conventional approach to all current operational retrievals, 3-D radiative transfer is used for interpreting the observed reflectances. General properties of the radiation reflected from the sides of an isolated cloud are discussed. As a proof of concept, the paper shows a few examples of radiation reflected from cloud fields generated by a simple stochastic cloud model with the prescribed vertically resolved microphysics. To retrieve the information about droplet sizes, we propose to use the probability density function of the droplet size distribution and its first two moments instead of the assumption about fixed values of the droplet effective radius. The retrieval algorithm is based on the Bayesian theorem that combines prior information about cloud structure and microphysics with radiative transfer calculations.

Marshak, A.↗
NASA NTRS2005

What Does Reflection from Cloud Sides Tell Us About Vertical Distribution of Cloud Droplet Sizes?

Cloud development, the onset of precipitation and the effect of aerosol on clouds depend on the structure of the cloud profiles of droplet size and phase. Aircraft measurements of cloud profiles are limited in their temporal and spatial extent. Satellites were used to observe cloud tops not cloud profiles with vertical profiles of precipitation-sized droplets anticipated from Cloudsat. The recently proposed CLAIM-3D satellite mission (cloud aerosol interaction mission in 3D) suggests to measure profiles of cloud microphysical properties by retrieving them from the solar and infrared radiation reflected or emitted from cloud sides. Inversion of measurements from the cloud sides requires rigorous understanding of the 3-dimensional (3D) properties of clouds. Here we discuss the reflected sunlight from the cloud sides and top at two wavelengths: one nonabsorbing to solar radiation (0.67 micrometers) and one with liquid water efficient absorption of solar radiation (2.1 micrometers). In contrast to the plane-parallel approximation, a conventional approach to all current operational retrievals, 3D radiative transfer is used for interpreting the observed reflectances. General properties of the radiation reflected from the sides of an isolated cloud are discussed. As a proof of concept, the paper shows a few examples of radiation reflected from cloud fields generated by a simple stochastic cloud model with the prescribed vertically resolved microphysics. To retrieve the information about droplet sizes, we propose to use the probability density function of the droplet size distribution and its first two moments instead of the assumption about fixed values of the droplet effective radius. The retrieval algorithm is based on the Bayesian theorem that combines prior information about cloud structure and microphysics with radiative transfer calculations.

Marshak, Alexander↗
NASA NTRS1990

A dipole moment of the microwave background as a cosmological effect

A spherically symmetrical Tolman-Bondi cosmological model is presented in which the curvature of space and the entropy variety with distance from the center. The dipole and quadrupole moments in the distribution of the microwave background radiation are calculated as a function of cosmic time and position of an observer, assuming that the distance to the horizon is much smaller than any characteristic scale in the model. The quadrupole moment is found to be affected mostly by the gradient in the curvature of space while the dipole moment is dominated by the gradient of entropy. The results indicate that the observed dipole in the microwave background may be cosmological in origin. Observational tests of this argument are suggested.

Paczynski, Bohdan↗
DOE OSTI2025

First Global Extraction of Generalized Parton Distributions from Experiment and Lattice Data with Next-to-Leading-Order Accuracy

We report the first global extraction of generalized parton distributions, GUMP 1.0, by combining deeply virtual Compton scattering and 𝜌-meson production data from Jefferson Lab and the Hadron-Electron Ring Accelerator with global fits of parton distribution functions, charge form factors, and lattice quantum chromodynamics simulations. Using a conformal moment space parametrization, we achieve a unified description across low- and high-𝑥 regions at next to leading order accuracy in perturbative corrections. The results provide state-of-the-art generalized parton distributions consistent with almost all known facts, enabling three-dimensional nucleon imaging in impact parameter space and, at the same time, establishing a benchmark for future theoretical and experimental studies of the nucleon structure.

Form factors↗
DOE OSTI2021

Linear Combination of Atomic Dipoles to Calculate the Bond and Molecular Dipole Moments of Molecules and Molecular Liquids

Here we report a linear combination of atomic dipole (LCAD) method for calculating the bond dipole moments of molecules. We show that the LCAD method reproduces the known molecular dipole moments of small to large molecules with a small error with respect to experimental and benchmark ab initio calculations, and molecular dipole distributions of bulk water that agree with maximally localized Wannier functions. The bond dipole moments derived from LCAD are also chemically interpretable in terms of the trend in bond ionicity in going from neutral to charged molecules. Moreover, the LCAD method accurately captures the influence of electric fields, supported by the correct trend in the change of the dipole moment under a uniform external electric field. The better grounding of bond dipole calculations indicates that it should also serve as a useful approach to bond dipole-field models used in catalysis or to reconstruct the small dipole of a H-terminated graphene flake.

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

Ten-moment fluid model for low-temperature magnetized plasmas

In this paper, a one-dimensional 10-moment multi-fluid plasma model is developed and applied to low-temperature magnetized plasmas. The 10-moment model solves for six anisotropic pressure terms, in addition to density and three components of fluid momentum, which allows the model to capture finite kinetic effects. The results are benchmarked against a 5-moment model, which assumes that the gas constituents follow a Maxwellian velocity distribution function (VDF), and a particle-in-cell Monte Carlo collision model that allows for arbitrary non-Maxwellian VDFs. The models are compared in a one-dimensional, low-temperature, partially magnetized plasma test case. The 10-moment results accurately reproduce the anisotropic temperature profile in low-temperature magnetized plasmas, where shear gradients exist due to the E×B drift. We discuss the mechanisms by which the anisotropic pressure can be generated in low-temperature magnetized plasmas. In addition, the importance of a self-consistent heat flux closure to the 10-moment model is studied, showing consistency with other models only when the assumptions of the underlying model are met. The 10-moment model allows for study of electron inertia effects and non-Maxwellian VDFs without the need for kinetic methods that are more computationally expensive.

Kuldinow, Derek Amur (ORCID:0000000319730196)↗
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.↗