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 253 records · Page 14

Integral-moment analysis of the BATSE gamma-ray burst intensity distribution

We have applied the technique of integral-moment analysis to the intensity distribution of the first 260 gamma-ray bursts observed by the Burst and Transient Source Experiment (BATSE) on the Compton Gamma Ray Observatory. This technique provides direct measurement of properties such as the mean, variance, and skewness of the convolved luminosity-number density distribution, as well as associated uncertainties. Using this method, one obtains insight into the nature of the source distributions unavailable through computation of traditional single parameters such as V/V(sub max)). If the luminosity function of the gamma-ray bursts is strongly peaked, giving bursts only a narrow range of luminosities, these results are then direct probes of the radial distribution of sources, regardless of whether the bursts are a local phenomenon, are distributed in a galactic halo, or are at cosmological distances. Accordingly, an integral-moment analysis of the intensity distribution of the gamma-ray bursts provides for the most complete analytic description of the source distribution available from the data, and offers the most comprehensive test of the compatibility of a given hypothesized distribution with observation.

Horack, John M.↗

Parallel expansion of a fuel pellet plasmoid

The problem of the assimilation of a cryogenic fuel pellet injected into a hot plasma is considered. Due to the transparency to ambient particles of the plasmoid, the localised region of high-density plasma created by ionisation of the ablated pellet material, electrons reach a ‘quasiequilibrium’ (QE) state which is characterised by a steady-state on the fastest collisional time scale. The simplified electron kinetic equation of the QE state is solved. Taking a velocity moment of the higher-order electron kinetic equation, which is valid on the expansion time scale, permits a fluid closure, yielding an evolution equation for the macroscopic parameters describing the QE distribution function. In contrast to the Braginskii equations, the closure does not require that electrons have a short mean free path compared with the size of density perturbations, and permits an anisotropic and highly non-Maxwellian distribution function. As the QE distribution function accounts for both trapped and passing electrons, the self-consistent electric potential that causes the expansion can be properly described, in contrast to earlier models of pellet plasmoid expansion with an unbounded potential. The plasmoid expansion is simulated using both a Vlasov model and a cold-fluid model for the ions. During the expansion plasmoid ions and electrons obtain nearly equal amounts of energy; as hot ambient electrons provide this energy in the form of collisional heating of plasmoid electrons, the expansion of a pellet plasmoid is expected to be a potent mechanism for the transfer of energy from electrons to ions on a time scale shorter than that of ion–electron thermalisation.

Physics↗

Relative advantages of thin-layer Navier-Stokes and interactive boundary-layer procedures

Numerical procedures for solving the thin-shear-layer Navier-Stokes equations and for the interaction of solutions to inviscid and boundary-layer equations are described and evaluated. To allow appraisal of the numerical and fluid dynamic abilities of the two schemes, they have been applied to one airfoil as a function of angle of attack at two slightly different Reynolds numbers. The NACA 0012 airfoil has been chosen because it allows comparison with measured lift, drag, and moment and with surface-pressure distributions. Calculations have been performed with algebraic eddy-viscosity formulations, and they include consideration of transition. The results are presented in a form that allows easy appraisal of the accuracy of both procedures and of the relative costs. The interactive procedure is computationally efficient but restrictive relative to the thin-layer Navier-Stokes procedure. The latter procedure does a better job of predicting drag than does the former. In both procedures, the location of transition is crucial for accurate or detailed computations, particularly at high angles of attack. When the upstream influence of pressure field through the shear layer is important, the thin-layer Navier-Stokes procedure has an edge over the interactive procedure.

Mehta, U.↗

A comparison of interactive boundary-layer and thin-layer Navier-Stokes procedures

Numerical procedures for solving the thin-shear-layer Navier-Stokes equations and for the interaction of solutions to inviscid and boundary-layer equations are described and evaluated. To allow appraisal of the numerical and fluid dynamic abilities of the two schemes, they have been applied to one airfoil as a function of angle of attack at two slightly different Reynolds numbers. The NACA 0012 airfoil has been chosen because it allows comparison with measured lift, drag, and moment and with surface-pressure distributions. Calculations have been performed with algebraic eddy-viscosity formulations, and they include consideration of transition. The results are presented in a form that allows easy appraisal of the accuracy of both procedures and of the relative costs. The interactive procedure is computationally efficient but restrictive relative to the thin-layer Navier-Stokes procedure. The latter procedure does a better job of predicting drag than does the former. In both procedures, the location of transition is crucial for accurate or detailed computations, particularly at high angles of attack. When the upstream influence of pressure field through the shear layer is important, the thin-layer Navier-Stokes procedure has an edge over the interactive procedure.

Mehta, U.↗

Probabilistic assessment of space trusses subjected to combined mechanical and thermal loads

A three-bay, space, cantilever truss is probabilistically evaluated to quantify the range of uncertainties of buckling loads and member forces due to nonuniform thermal loads, applied loads and moments (mechanical loads), and combination of both. The truss members are assumed to be made from Aluminum tubes or high modulus graphite-fiber/intermediate modulus epoxy-matrix composite tubes. Cumulative distribution function results show that certain combinations of thermal loads with mechanical loads reduce the probabilistic buckling loads and increase the magnitude of the member axial forces for the aluminum truss. The same trend is observed for the composite truss as well, as however, the thermal effects on the probabilistic buckling loads and member axial forces are not as substantial as that for an aluminum truss. This can be attributed to the large differences in the values of coefficient of thermal expansion. Finally, the sensitivities associated with the uncertainties in the structural, material, and load variables (primitive variables) are investigated. They show that buckling loads and member axial forces are most sensitive to the uncertainties in spacial (geometry) variables.

Pai, Shantaram S.↗

Lattice Observables for Parton Distributions

In recent years there has been a momentous surge of interest in the determination of the parton content of nucleons using lattice QCD. Scalar field theory is a convenient toy model to understand the relation between Euclidean correlators and Parton Distribution Functions (PDFs). In this talk we review the factorization formulae that allow the extraction of PDFs from the pseudo-PDFs and quasi-PDFs computed on the lattice. The scalar theory provides a simplified setting where the main physical properties are easily analysed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Modified Spectral Fatigue Methods for S-N Curves With MIL-HDBK-5J Coefficients

The rainflow method is used for counting fatigue cycles from a stress response time history, where the fatigue cycles are stress-reversals. The rainflow method allows the application of Palmgren-Miner's rule in order to assess the fatigue life of a structure subject to complex loading. The fatigue damage may also be calculated from a stress response power spectral density (PSD) using the semi-empirical Dirlik, Single Moment, Zhao-Baker and other spectral methods. These methods effectively assume that the PSD has a corresponding time history which is stationary with a normal distribution. This paper shows how the probability density function for rainflow stress cycles can be extracted from each of the spectral methods. This extraction allows for the application of the MIL-HDBK-5J fatigue coefficients in the cumulative damage summation. A numerical example is given in this paper for the stress response of a beam undergoing random base excitation, where the excitation is applied separately by a time history and by its corresponding PSD. The fatigue calculation is performed in the time domain, as well as in the frequency domain via the modified spectral methods. The result comparison shows that the modified spectral methods give comparable results to the time domain rainflow counting method.

Irvine, Tom↗

Direct Numerical Simulation of Transitional Multicomponent-Species Gaseous and Multicomponent-Liquid Drop-Laden Mixing

A model of multicomponent-liquid (MC-liquid) drop evaporation in a three-dimensional mixing layer is here exercised at larger Reynolds numbers than in a previous study, and transitional states are obtained. The gas phase is followed in an Eulerian frame and the multitude of drops is described in a Lagrangian frame. Complete coupling between phases is included with source terms in the gas conservation equations accounting for the drop/flow interaction in terms of drop drag, drop heating and species evaporation. The liquid composition, initially specified as a single-Gamma (SG) probability distribution function (PDF) depending on the molar mass is allowed to evolve into a linear combination of two SGPDFs, called the double-Gamma PDF (DGPDF). The compositions of liquid and vapor emanating from the drops are calculated through four moments of the DGPDFs, which are drop-specific and location-specific, respectively. The mixing layer is initially excited to promote the double pairing of its four initial spanwise vortices into an ultimate vortex in which small scales proliferate. Simulations are performed for four liquids of different compositions and the effect of the initial mass loading and initial free-stream gas temperature are explored. For reference, Simulations are also performed for gaseous multicomponent mixing layers for which the effect of Reynolds number is investigated. The results encompass examination of the global layer characteristics, flow visualizations and homogeneous-plane statistics at transition. Comparisons are performed with previous pre-transitional MC-liquid simulations and with transitional single-component (SC) liquid studies. It is found that MCC flows at transition, the classical energy cascade is of similar strength, but that the smallest scales contain orders of magnitude less energy than SC flows, which is confirmed by the larger viscous dissipation in the former case. Contrasting to pre-transitional MC flows, the vorticity and drop organization depend on the initial gas temperature, this being due to the drop/turbulence coupling. The vapor-composition mean molar mass and standard deviation distributions strongly correlate with the initial liquid-composition PDF; such a correlation only exists for the magnitude of the mean but not for that of the standard deviation. Unlike in pre-transitional situations, regions of large composition standard deviation no longer necessarily coincide with regions of large mean molar mass. The kinetic energy, rotational and composition characteristics, and dissipation are liquid specific and the variation among liquids is amplified with increasing free-stream gas temperature. Eulerian and Lagrangian statistics of gas-phase quantities show that the different. Observation framework may affect the perception of the flow characteristics. The gas composition, of which the first four moments are calculated, is shown to be close to, but distinct from a SGPDF. The PDF of the scalar dissipation rate is calculated for drop-laden layers and is shown to depart more significantly from the typically assumed Gaussian in gaseous flows than experimentally measured gaseous scalar dissipation rates, this being attributed to the increased heterogeneity due to drop/flow interactions.

tubulent multicompenent-species↗

Single-Molecule Kinetics of Styrene Hydrogenation on Silica-Supported Vanadium: The Role of Disorder for Single-Atom Catalysts

We report a theoretical approach for the study of supported atom catalysis is developed based on recent advances in the study of single-molecule kinetics. This view is particularly useful in exhibiting the role of disorder in single-atom and single-site catalysts on amorphous supports. The distribution of passage times (or waiting times) through a complex catalytic network originating from a set of coupled active sites is described by a probability distribution function (PDF), f(t), that reflects the local environment of the reaction center. An efficient algorithm is developed based on the linear algebra of the Markov transition matrix that produces f(t) or its moments. The kinetics of the hydrogenation reaction of styrene on an organovanadium(III) catalyst supported on amorphous silica is studied. A kinetic model consisting of three intertwined catalytic cycles emanating from three chemically distinct active sites is proposed to describe the chemistry. Density functional theory (DFT) calculations are employed to determine the free energy barriers of the reactions, which are used to construct the rate coefficient matrix. The disorder induced by the amorphous support material is divided into a low-dimensional short-range component reflecting the covalent structures near the reaction center and a weaker long-range component modeling the bulk randomness. The results are computed and analyzed for a wide range of concentration values and disorder scenarios. The unusual structure in the f(t) PDF is found to occur for certain cases that reveal the contribution of multiple catalytic pathways acting in concert.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Gluon Sivers function from forward exclusive χ c 1 photoproduction on unpolarized protons

Exclusive production of a χ c 1 axial vector quarkonia in photon-proton scattering at high energies requires a C -odd t -channel exchange. In the limit of vanishing momentum transfer this occurs either via the exchange of a photon, the Primakoff process, where the spin of the proton does not change. For axial-vector meson production, as a consequence of the Landau-Yang theorem, the Primakoff cross section is finite as t → 0 . Alternatively, a C -odd spin dependent Odderon can be exchanged, which involves a spin flip of the proton. The resulting cross section is related to the square of the collinear trigluon correlator or the k ⊥ -moment of the gluon Sivers function. Using two models for the gluon Sivers function from the literature we compute the ratio of Sivers to Primakoff cross sections and the angular coefficient λ θ governing the angular distribution of the χ c 1 → J / ψ + γ decay as functions of x . We point out that these observables constrain the magnitude of the gluon Sivers function at small x which could be accessed in electron-proton scattering and ultraperipheral proton-proton and nucleus-proton collisions. Published by the American Physical Society 2025

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Particle Acceleration in Relativistic Alfvénic Turbulence

Abstract Strong magnetically dominated Alfvénic turbulence is an efficient engine of nonthermal particle acceleration in a relativistic collisionless plasma. We argue that in the limit of strong magnetization, the type of energy distribution attained by accelerated particles depends on the relative strengths of turbulent fluctuations δ B 0 and the guide field B 0 . If δ B 0 ≪ B 0 , the particle magnetic moments are conserved, and the acceleration is provided by magnetic curvature drifts. Curvature acceleration energizes particles in the direction parallel to the magnetic field lines, resulting in log-normal tails of particle energy distribution functions. Conversely, if δ B 0 ≳ B 0 , interactions of energetic particles with intense turbulent structures can scatter particles, creating a population with large pitch angles. In this case, magnetic mirror effects become important, and turbulent acceleration leads to power-law tails of the energy distribution functions.

79 ASTRONOMY AND ASTROPHYSICS↗

Nonuniform thermal environmental effects on space truss structural reliability

A three-bay, space, cantilever truss is probabilistically evaluated to quantify the range of uncertainties of buckling loads and member forces due to nonuniform thermal loads, applied loads and moments (mechanical loads), and combination of both. The truss members are assumed to be made from (1) aluminum tubes or (2) high modulus graphite-fiber/intermediate modulus epoxy-matrix composite tubes. Cumulative distribution function results show that certain combinations of thermal loads with mechanical loads reduce the probabilistic buckling loads and increase the magnitude of the member axial forces for aluminum truss. However, the reverse is true for composite truss due to very low coefficient of thermal expansion and higher stiffness of the composite. Finally, the sensitivities associated with the uncertainties in the structural, material, and load variables (primitive variables) are investigated. They show that buckling loads and member axial forces are most sensitive to the uncertainties in spacial (geometry) variables.

Pai, Shantaram S.↗

Probabilistic structural analysis of space truss structures for nonuniform thermal environmental effects

A three-bay, space, cantilever truss is probabilistically evaluated to quantify the range of uncertainties of vibration frequencies and truss free end vertical displacements due to nonuniform thermal loads, applied loads and moments (mechanical loads), and combination of both. The truss members are assumed to be made from (1) aluminum tubes or (2) high modulus graphite-fiber/intermediate modulus epoxy-matrix composite tubes. Cumulative distribution function results show that certain combination of thermal loads with mechanical loads reduce the probabilistic vibration frequencies and increase the magnitude of the truss free end vertical displacements for the aluminum truss. However, for the composite truss, the thermal effects on the probabilistic vibration frequencies are not as substantial as that for an aluminum truss. This can be attributed to the large differences in the values of coefficient of thermal expansion. Finally, the sensitivities associated with the uncertainties in the structural, material and load variables (primative variables) are investigated. Results indicate that vibration frequencies and truss free end vertical displacements are sensitive to the uncertainties in spacial (geometry) variables as well as member modulus.

Pai, Shantaram S.↗

Saturn's Magnetosphere and Properties of Upstream Flow at Titan: Preliminary Results

Using Cassini Plasma Spectrometer (CAPS) Ion Mass Spectrometer (IMS) measurements, we present the ion fluid properties and its ion composition of the upstream flow for Titan's interaction with Saturn's magnetosphere. A 3D ion moments algorithm is used which is essentially model independent with only requirement is that ion flow is within the CAPS IMS 2(pi) steradian field-of-view (FOV) and that the ion 'velocity distribution function (VDF) be gyrotropic. These results cover the period from TA flyby (2004 day 300) to T22 flyby (2006 363). Cassini's in situ measurements of Saturn's magnetic field show it is stretched out into a magnetodisc configuration for Saturn Local Times (SLT) centered about midnight local time. Under those circumstances the field is confined near the equatorial plane with Titan either above or below the magnetosphere current sheet. Similar to Jupiter's outer magnetosphere where a magnetodisc configuration applies, one expects the heavy ions within Saturn's outer magnetosphere to be confined within a few degrees of the current sheet while at higher magnetic latitudes protons should dominate. We show that when Cassini is between dusk-midnight-dawn local time and spacecraft is not within the current sheet that light ions (H, 142) tend to dominate the ion composition for the upstream flow. If true, one may expect the interaction between Saturn's magnetosphere, locally devoid of heavy ions and Titan's upper atmosphere and exosphere to be significantly different from that for Voyager 1, TA and TB when heavy ions were present in the upstream flow. We also present observational evidence for Saturn's magnetosphere interaction with Titan's extended H and H2 corona which can extend approx. 1 Rs from Titan.

Sittler, E. C., Jr.↗

Gaussian Kernel Methods for Seismic Fragility and Risk Assessment of Mid-Rise Buildings

Seismic fragility functions can be evaluated using the cloud analysis method with linear regression which makes three fundamental assumptions about the relation between structural response and seismic intensity: log-linear median relationship, constant standard deviation, and Gaussian distributed errors. While cloud analysis with linear regression is a popular method, the degree to which these individual and compounded assumptions affect the fragility and the risk of mid-rise buildings needs to be systematically studied. This paper conducts such a study considering three building archetypes that make up a bulk of the building stock: RC moment frame, steel moment frame, and wood shear wall. Gaussian kernel methods are employed to capture the data-driven variations in the median structural response and standard deviation and the distributions of residuals with the intensity level. With reference to the Gaussian kernels approach, it is found that while the linear regression assumptions may not affect the fragility functions of lower damage states, this conclusion does not hold for the higher damage states (such as the Complete state). In addition, the effects of linear regression assumptions on the seismic risk are evaluated. For predicting the demand hazard, it is found that the linear regression assumptions can impact the computed risk for larger structural response values. However, for predicting the loss hazard with downtime as the decision variable, linear regression can be considered adequate for all practical purposes.

58 GEOSCIENCES↗

Thermophysical properties of FLiBe using moment tensor potentials

Fluoride salts are prospective materials for applications in some next-generation nuclear reactors and their thermophysical properties at various conditions are of interest. Experimental measurement of the properties of these salts is often difficult and, in some cases, unfeasible due to challenges from high temperatures, impurity control, and corrosivity. Therefore, accurate theoretical methods are needed for fluoride salt property prediction. In this work, we used moment tensor potentials (MTP) to approximate the potential energy surface of eutectic FLiBe (66.6% LiF – 33.3% BeF2) predicted by the ab initio (DFT-D3) method. Here, we then used the developed potential and molecular dynamics to obtain several thermophysical properties of FLiBe, including radial distribution functions, density, self-diffusion coefficients, thermal expansion, specific heat capacity, bulk modulus, viscosity, and thermal conductivity. Our results show that the MTP potential approximates the potential energy surface accurately and the overall approach yields very good agreement with experimental values. The converged fitting can be obtained with less than 600 configurations generated from DFT calculations, which data can be generated in just 1200 core hours on today's typical processors. The MTP potential is faster than many machine learning potentials and about one order of magnitude slower than widely used empirical molten salt potentials such as Tosi/Fumi.

36 MATERIALS SCIENCE↗

On isostatic geoid anomalies

In regions of slowly varying lateral density changes, the gravity and geoid anomalies may be expressed as power series expansions in topography. Geoid anomalies in isostatically compensated regions can be directly related to the local dipole moment of the density-depth distribution. This relationship is used to obtain theoretical geoid anomalies for different models of isostatic compensation. The classical Pratt and Airy models give geoid height-elevation relationships differing in functional form but predicting geoid anomalies of comparable magnitude. The thermal cooling model explaining ocean floor subsidence away from mid-ocean ridges predicts a linear age-geoid height relationship of 0.16 m/m.y. Geos 3 altimetry profiles were examined to test these theoretical relationships. A profile over the mid-Atlantic ridge is closely matched by the geoid curve derived from the thermal cooling model. The observed geoid anomaly over the Atlantic margin of North America can be explained by Airy compensation. The relation between geoid anomaly and bathymetry across the Bermuda Swell is consistent with Pratt compensation with a 100-km depth of compensation.

Haxby, W. F.↗

Feautrier solution of the electron transport equation

An efficient, accurate, and stable solution to the electron transport problem is presented. The two-stream approximation is discussed in detail, but the solution is developed for the general multistream case. Applied to the photoelectron problem, it is shown that accurate escape fluxes are obtained using only 10 altitudes distributed over the region 100 - 1000 km. The local equilibrium, two-stream, and multistream solutions to the photoelectron problem are compared. In addition, a recursion relation for the moments of the elastic angular scattering phase function is presented.

Link, R.↗