SEARCH · Engineering Papers

Results for “Probability 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.

SourceResearch connectionFrom yearThrough year

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

At least 397 records · Page 22

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Entropy production in model colloidal suspensions under shear via the fluctuation theorem

Dissipative systems often exhibit novel and unexpected properties. This is, for instance, the case of simple liquids, which, when subjected to shear and after reaching a steady state, can exhibit a negative entropy production over finite length scales and timescales. This result, among others, is captured by nonequilibrium relations known as fluctuation theorems. Using nonequilibrium molecular dynamics simulations, we examine how, by fine-tuning the properties of the components of a complex fluid, we can steer the nonequilibrium response of the fluid. More specifically, we show how we control the nonequilibrium probability distribution for the shear stress and, in turn, how often states with a negative entropy production can occur. To achieve this, we start by characterizing how the size for the liquid matrix impacts the probability of observing negative entropy states, as well as the timescale over which these can be observed. We then measure how the addition of larger particles to this liquid matrix, i.e., simulating a model colloidal suspension, results in an increase in the occurrence of such states. As a result, this suggests how modifications in the composition of the mixture and in the properties of its components lead to an increase in the probability of observing states of negative entropy production and, thus, for the system to run in reverse.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Bayesian reconstruction of anisotropic flow fluctuations at fixed impact parameter

The cumulants of the distribution of anisotropic flow are measured accurately in Pb+Pb collisions at the LHC as a function of centrality classifiers (charged multiplicity and/or transverse energy). Using Bayesian inference, we reconstruct from these measurements the probability distribution of anisotropic flow in the ``theorists' frame'' where the impact parameter has a fixed magnitude and orientation, up to ∼70% centrality. The variation of flow fluctuations with impact parameter displays direct evidence of viscous damping, which is larger for higher Fourier harmonics, in line with expectations from hydrodynamics. We use intensive measures of non-Gaussian flow fluctuations, which have reduced dependence on centrality. Here, we infer from ATLAS data the magnitude of these intensive non-Gaussianities in each Fourier harmonic. They provide data-driven estimates of response coefficients to initial anisotropies, without resorting to any specific microscopic model of initial conditions. These estimates agree with viscous hydrodynamic calculations.

Bayesian methods↗

DOE OSTI2022

Statistics of Green's functions on a disordered Cayley tree and the validity of forward scattering approximation

The accuracy of the forward scattering approximation for two-point Green's functions of the Anderson localization model on the Cayley tree is studied. A relationship between the moments of the Green's function and the largest eigenvalue of the linearized transfer-matrix equation is proved in the framework of the supersymmetric functional-integral method. The new large-disorder approximation for this eigenvalue is derived and its accuracy is established. Using this approximation the probability distribution of the two-point Green's function is found and compared with that in the forward scattering approximation (FSA). It is shown that FSA overestimates the role of resonances and thus the probability for the Green's function to be significantly larger than its typical value. The error of FSA increases with increasing the distance between points in a two-point Green's function.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2021

Bayesian estimations of orientation distribution functions from small-angle scattering enable direct prediction of mechanical stress in anisotropic materials

Properties of soft materials are influenced by their anisotropic structuring under nonequilibrium fields. Although anisotropic structure-property relationships have been extensively explored theoretically, comparison to experiments requires determination of the microstructural orientation probability distribution function (OPDF) of microstructural elements. Small angle scattering (SAS) measurements encode information about the OPDF, but tools to navigate this connection are incomplete. Here, we develop and validate an explicit framework to link arbitrary OPDFs to SAS measurements. Specifically, we propose, validate, and apply a method, maximum a posteriori scattering inference (MAPSI), whereby the OPDF may be obtained from SAS measurements using a Bayesian estimation method. Using this method, we obtain estimates of the full 3D OPDF for two model semidilute fd-virus (rodlike) dispersions at concentrations that are approximately equal to and twice the overlap concentration. From the OPDF, we calculate its second and fourth moments and compare these to predictions for a dilute suspension of rigid rods and to a recent theory for semidilute suspensions. Finally, we use both the theoretical and measured moments to calculate the stress, both for dilute and semidilute suspensions. These predictions are not only compared to each other, but also to measured values of the shear stress, and point to new insights into the behavior of suspensions of highly elongated particles in the transition between dilute and semidilute behavior. We also use this new framework to provide perspective on the connection between scalar parameterizations of scattering and the OPDF that have frequently been used in the past. Here, the new tools developed in this work provide an unprecedented path toward experimental validation of dynamical theories of rodlike colloids and polymers, and for measurement of nonequilibrium structures and stresses of other complex fluids and soft materials with SAS.

36 MATERIALS SCIENCE↗

NASA NTRS1990

BATSE software for the analysis of the gamma ray burst spatial distribution

The Burst and Transient Source Experiment (BATSE) on the Gamma Ray Observatory (GRO) is designed to study astronomical gamma ray sources and to provide better positional, spectral, and time resolution about these objects than has previously been possible from one experiment. The procedure to be used in the analysis of the gamma ray burst spatial distribution is presented. Data is input from BATSE via the Gamma Ray Burst Catalog (listing individual burst positions, flux values, and associated errors) and the Sky Sensitivity Map (which summarizes observational selection effects in table format). A FORTRAN program generates Monte Carlo burst catalogs, which are models to be compared to the actual distribution. The Monte Carlo models are then filtered through the Sky Sensitivity Map so that they suffer from the same selection effects as the actual catalog data. Additionally, each burst position is converted into a probability distribution to mimic BATSE positional sensitivity. The Burst Catalog, Monte Carlo burst catalog, and Sky Sensitivity Map are then passed onto an IDL program that compares the catalogs for statistical significance. The Sky Sensitivity Map is used to estimate how often each sky area is observed above the minimum flux level in question. Each burst found in this sky area is then weighted according to the frequency with which this sky area is observed. The catalogs are then compared via tests of homogeneity (based on their radial distributions) and isotropy (based upon their angular distributions). The results of the statistical comparisons along with graphs and charts of the summaries, are output from the IDL program for study.

Hakkila, Jon↗

NASA NTRSDate not supplied

Inter-Calibrating the CERES SW Instrument Fluxes Using Deep Convective Clouds(DCC) Within the Same Sun-Synchronous Orbit

The NASA CERES EBAF product provides the scientific community observed TOA fluxes to monitor the Earth’s energy imbalance and to validate climate models. To provide a seamless EBAF 24-year record, the CERES instrument calibration must be stable and consistent across satellite records. Once the Terra and Aqua spacecraft are decommissioned at the end of 2025, the CERES project will rely on CERES instrument observations on the SNPP and NOAA-20 1:30 PM sun-synchronous orbits. The SNPP and NOAA-20 satellite orbit placement will prevent any time matched observations for inter-calibration efforts. The future Libera instrument, which will continue the CERES record and scheduled for launch in 2028 onboard the NOAA-22 satellite, will also need to be inter-calibrated without the aid of time-matched observations. Deep convective clouds (DCC) are the most Lambertian, brightest, tropical Earth invariant targets located at the tropopause making them ideal to radiometrically scale the CERES SW observed radiance to a common reference. The SNPP, NOAA-20 and future NOAA-22 satellites are in the same 16-day repeatable orbits, allowing the DCC targets to be observed with the same angular configuration. The empirically derived SNPP CERES SW channel Bidirectional Reflectance Distribution Function (BRDF) can be easily be applied to the NOAA-20 CERES SW channel radiances since they observe the nearly the same DCC systems. The DCC BRDF corrected radiances are analyzed collectively into radiance probability distribution functions (PDF). By comparing the PDF statistics, the SNPP and NOAA-20 CERES SW channel can be inter-calibrated. The inter-calibration coefficients are validated with the CERES instrument calibration team coefficients to optimize the methodology.

David Doelling↗

DOE OSTI2020

Nearest neighbour distributions: New statistical measures for cosmological clustering

ABSTRACT The use of summary statistics beyond the two-point correlation function to analyse the non-Gaussian clustering on small scales, and thereby, increasing the sensitivity to the underlying cosmological parameters, is an active field of research in cosmology. In this paper, we explore a set of new summary statistics – the k-Nearest Neighbour Cumulative Distribution Functions (kNN-CDF). This is the empirical cumulative distribution function of distances from a set of volume-filling, Poisson distributed random points to the k-nearest data points, and is sensitive to all connected N-point correlations in the data. The kNN-CDF can be used to measure counts in cell, void probability distributions, and higher N-point correlation functions, all using the same formalism exploiting fast searches with spatial tree data structures. We demonstrate how it can be computed efficiently from various data sets – both discrete points, and the generalization for continuous fields. We use data from a large suite of N-body simulations to explore the sensitivity of this new statistic to various cosmological parameters, compared to the two-point correlation function, while using the same range of scales. We demonstrate that the use of kNN-CDF improves the constraints on the cosmological parameters by more than a factor of 2 when applied to the clustering of dark matter in the range of scales between 10 and $40\, h^{-1}\, {\rm Mpc}$. We also show that relative improvement is even greater when applied on the same scales to the clustering of haloes in the simulations at a fixed number density, both in real space, as well as in redshift space. Since the kNN-CDF are sensitive to all higher order connected correlation functions in the data, the gains over traditional two-point analyses are expected to grow as progressively smaller scales are included in the analysis of cosmological data, provided the higher order correlation functions are sensitive to cosmology on the scales of interest.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS2012

Global Warming Induced Changes in Rainfall Characteristics in IPCC AR5 Models

Changes in rainfall characteristic induced by global warming are examined from outputs of IPCC AR5 models. Different scenarios of climate warming including a high emissions scenario (RCP 8.5), a medium mitigation scenario (RCP 4.5), and 1% per year CO2 increase are compared to 20th century simulations (historical). Results show that even though the spatial distribution of monthly rainfall anomalies vary greatly among models, the ensemble mean from a sizable sample (about 10) of AR5 models show a robust signal attributable to GHG warming featuring a shift in the global rainfall probability distribution function (PDF) with significant increase (>100%) in very heavy rain, reduction (10-20% ) in moderate rain and increase in light to very light rains. Changes in extreme rainfall as a function of seasons and latitudes are also examined, and are similar to the non-seasonal stratified data, but with more specific spatial dependence. These results are consistent from TRMM and GPCP rainfall observations suggesting that extreme rainfall events are occurring more frequently with wet areas getting wetter and dry-area-getting drier in a GHG induced warmer climate.

Lau, William K. M.↗

DOE OSTI2025

LTAU-FF: Loss Trajectory Analysis for Uncertainty in atomistic Force Fields

Model ensembles are effective tools for estimating prediction uncertainty in deep learning atomistic force fields. However, their widespread adoption is hindered by high computational costs and overconfident error estimates. In this work, we address these challenges by leveraging distributions of per-sample errors obtained during training and employing a distance-based similarity search in the model latent space. Our method, which we call LTAU (Loss Trajectory Analysis for Uncertainty), efficiently estimates the full probability distribution function of errors for any test point using the logged training errors, achieving speeds that are 2–3 orders of magnitudes faster than typical ensemble methods and allowing it to be used for tasks where training or evaluating multiple models would be infeasible. We apply LTAU towards estimating parametric uncertainty in atomistic force fields (LTAU-FF), demonstrating that it produces well-calibrated confidence intervals and predicts errors that correlate strongly with the true errors for data near the training domain. Furthermore, we show that the errors predicted by LTAU-FF can be used in practical applications for detecting out-of-domain data, tuning model performance, and predicting failure during simulations. We believe that LTAU will be a valuable tool for uncertainty quantification in atomistic force fields and is a promising method that should be further explored in other domains of machine learning.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1979

Simulation of the pressure field near a jet by randomly distributed vortex rings

Fluctuations of the pressure field in the vicinity of a jet are simulated numerically by a flow model consisting of axially symmetric vortex rings with viscous cores submerged in a uniform stream. The time interval between the shedding of successive vortices is taken to be a random variable with a probability distribution chosen to match that from experiments. It is found that up to 5 diameters downstream of the jet exit, statistics of the computed pressure field are in good agreement with experimental results. Statistical comparisons are provided for the overall sound pressure level, the peak amplitude, and the Strouhal number based on the peak frequency of the pressure signals.

Fung, Y. T.↗

NASA NTRS1989

Infrared backscattering

All particles in the atmosphere are not spherical. Moreover, the scattering properties of randomly oriented nonspherical particles are not equivalent to those of spherical particles no matter how the term equivalent is defined. This is especially true for scattering in the backward direction and at the infrared wavelengths at which some atmospheric particles have strong absorption bands. Thus calculations based on Mie theory of infrared backscattering by dry or insoluble atmospheric particles are suspect. To support this assertion, it was noted that peaks in laboratory-measured infrared backscattering spectra show appreciable shifts compared with those calculated using Mie theory. One example is ammonium sulfate. Some success was had in modeling backscattering spectra of ammonium sulfate particles using a simple statistical theory called the continuous distribution of ellipsoids (CDE) theory. In this theory, the scattering properties of an ensemble are calculated. Recently a modified version of this theory was applied to measured spectra of scattering by kaolin particles. The particles were platelike, so the probability distribution of ellipsoidal shapes was chosen to reflect this. As with ammonium sulfate, the wavelength of measured peak backscattering is shifted longward of that predicted by Mie theory.

Bohren, Craig F.↗

NASA NTRS2002

Risk-Based Probabilistic Approach to Aeropropulsion System Assessment

In an era of shrinking development budgets and resources, where there is also an emphasis on reducing the product development cycle, the role of system assessment, performed in the early stages of an engine development program, becomes very critical to the successful development of new aeropropulsion systems. A reliable system assessment not only helps to identify the best propulsion system concept among several candidates, it can also identify which technologies are worth pursuing. This is particularly important for advanced aeropropulsion technology development programs, which require an enormous amount of resources. In the current practice of deterministic, or point-design, approaches, the uncertainties of design variables are either unaccounted for or accounted for by safety factors. This could often result in an assessment with unknown and unquantifiable reliability. Consequently, it would fail to provide additional insight into the risks associated with the new technologies, which are often needed by decision makers to determine the feasibility and return-on-investment of a new aircraft engine. In this work, an alternative approach based on the probabilistic method was described for a comprehensive assessment of an aeropropulsion system. The statistical approach quantifies the design uncertainties inherent in a new aeropropulsion system and their influences on engine performance. Because of this, it enhances the reliability of a system assessment. A technical assessment of a wave-rotor-enhanced gas turbine engine was performed to demonstrate the methodology. The assessment used probability distributions to account for the uncertainties that occur in component efficiencies and flows and in mechanical design variables. The approach taken in this effort was to integrate the thermodynamic cycle analysis embedded in the computer code NEPP (NASA Engine Performance Program) and the engine weight analysis embedded in the computer code WATE (Weight Analysis of Turbine Engines) with the fast probability integration technique (FPI). FPI was developed by Southwest Research Institute under contract with the NASA Glenn Research Center. The results were plotted in the form of cumulative distribution functions and sensitivity analyses and were compared with results from the traditional deterministic approach. The comparison showed that the probabilistic approach provides a more realistic and systematic way to assess an aeropropulsion system. The current work addressed the application of the probabilistic approach to assess specific fuel consumption, engine thrust, and weight. Similarly, the approach can be used to assess other aspects of aeropropulsion system performance, such as cost, acoustic noise, and emissions. Additional information is included in the original extended abstract.

Tong, Michael T.↗

DOE OSTI2023

Bayesian weighting of climate models based on climate sensitivity

Using climate model ensembles containing members that exhibit very high climate sensitivities to increasing CO 2 concentrations can result in biased projections. Various methods have been proposed to ameliorate this ‘hot model’ problem, such as model emulators or model culling. Here, we utilize Bayesian Model Averaging as a framework to address this problem without resorting to outright rejection of models from the ensemble. Taking advantage of multiple lines of evidence used to construct the best estimate of the earth’s climate sensitivity, the Bayesian Model Averaging framework produces an unbiased posterior probability distribution of model weights. The updated multi-model ensemble projects end-of-century global mean surface temperature increases of 2 °C for a low emissions scenario (SSP1-2.6) and 5 °C for a high emissions scenario (SSP5-8.5). These estimates are lower than those produced using a simple multi-model mean for the CMIP6 ensemble. The results are also similar to results from a model culling approach, but retain some weight on low-probability models, allowing for consideration of the possibility that the true value could lie at the extremes of the assessed distribution. Our results showcase Bayesian Model Averaging as a path forward to project future climate change that is commensurate with the available scientific evidence.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2024

Path integrals, complex probabilities and the discrete Weyl representation

Abstract A discrete formulation of the real-time path integral as the expectation value of a functional of paths with respect to a complex probability on a sample space of discrete valued paths is explored. The formulation in terms of complex probabilities is motivated by a recent reinterpretation of the real-time path integral as the expectation value of a potential functional with respect to a complex probability distribution on cylinder sets of paths. The discrete formulation in this work is based on a discrete version of the Weyl algebra that can be applied to any observable with a finite number of outcomes. The origin of the complex probability in this work is the completeness relation. In the discrete formulation the complex probability exactly factors into products of conditional probabilities and exact unitarity is maintained at each level of approximation. The approximation of infinite dimensional quantum systems by discrete systems is discussed. The method is illustrated by applying it to scattering theory and quantum field theory. The implications of these applications for quantum computing is discussed.

Physics↗

DOE OSTI2025

Unavoidable multilevel biaxial symmetry breaking in chiral hybrid liquid crystals

Chiral nematic or cholesteric liquid crystals (LCs) are fluid mesophases with long-ranged orientational order featuring a quasilayered periodicity imparted by a helical director configuration but lacking long-range positional order. Doping molecular cholesteric LCs with strongly anisotropic uniaxial colloidal particles adds another level of complexity because of the interplay between weak surface-anchoring boundary conditions and bulk-based elastic distortions near the particle-LC interface. Using cylindrical colloidal disks and rods with different geometric shapes and surface conditions, we demonstrate that these colloidal inclusions generically exhibit biaxial orientational probability distributions which may impart anomalously strong local biaxiality onto the hybrid cholesteric LC structure. Unlike nonchiral hybrid molecular-colloidal LCs, where biaxial order emerges only at critical colloid volume fractions exceeding some uniaxial-biaxial transition value, the orientational probability of the colloidal inclusions immersed in chiral nematic hosts is unambiguously biaxial even at infinite dilution. We demonstrate that the colloids induce local biaxial perturbations within the molecular orientational order of the LC host medium which strongly enhances the weak but native biaxial order of chiral nematic LC induced by the chiral symmetry breaking of the director field. With the help of analytical modeling and computer simulations based on the Landau–de Gennes free energy of the host LC around the colloids, we rationalize the observed multilevel biaxial order and conclude that it is not only unavoidable but also strongly enhanced compared to both achiral hybrid LCs and purely molecular cholesteric LCs.

biaxial liquid crystals↗

DOE OSTI2026

Wavelet flow for extragalactic foreground simulations

Extragalactic foregrounds in cosmic microwave background (CMB) observations are both a source of cosmological and astrophysical information and a nuisance to the CMB. Effective field-level modeling that captures their non-Gaussian statistical distributions is increasingly important for optimal information extraction, particularly given the low-noise observations from current and upcoming experiments. Here, we explore the use of Wavelet Flow (WF) models to tackle the novel task of modeling the field-level probability distributions of multi-component CMB secondaries and foregrounds. Specifically, we jointly train correlated CMB lensing convergence (κ) and cosmic infrared background (CIB) maps with a WF model and obtain a network that statistically recovers the input to high accuracy — the trained network generates samples of κ and CIB fields whose average power spectra are within a few percent of the inputs across all scales, and whose Minkowski functionals are similarly accurate compared to the inputs. Leveraging the multiscale architecture of these models, we fine-tune both the model parameters and the priors at each scale independently, optimizing performance across different resolutions. These results demonstrate that WF models can accurately simulate correlated components of CMB secondaries, supporting improved analysis of cosmological data. Our code and trained models can be found on this GitHub repo.

cosmological simulations↗

DOE OSTI2025

Evolution of in-medium baryon-baryon scattering cross sections and stiffness of dense nuclear matter from Bayesian analyses of FOPI proton-flow excitation functions

Within a Bayesian statistical framework using a Gaussian process (GP) emulator for an isospin-dependent Boltzmann-Uehling-Uhlenbeck (IBUU) transport model simulator of heavy-ion reactions with momentum-independent Skyrme interactions, we infer from the proton directed and elliptical flow in mid-central Au + Au reactions, measured at beam energies from 150 to 1200 MeV/nucleon by the FOPI Collaboration, the posterior probability distribution functions (PDFs) of the in-medium baryon-baryon scattering cross section (BBSCS) modification factor 𝑋 (with respect to their free-space values) and the stiffness parameter 𝐾 of dense nuclear matter. Here, we find that the most probable value of 𝑋 evolves from around 0.7 to 1.0 as the beam energy 𝐸 beam /𝐴 increases. On the other hand, the posterior PDF⁡(𝐾) may have dual peaks having roughly the same height or extended shoulders at high 𝐾 values. More quantitatively, the posterior PDF⁡(𝐾) changes from having a major peak around 220 MeV, characterizing a soft EOS in the reaction at 𝐸 beam /𝐴 = 150 MeV, to one that peaks around 320 MeV, indicating a stiff EOS in the reactions at 𝐸 beam /𝐴 higher than about 600 MeV. The transition from soft to stiff happens in mid-central Au + Au reactions at beam energies around 250 MeV/nucleon in which 𝐾 = 220 MeV and 𝐾 = 320 MeV are approximately equally probable. Altogether, the FOPI proton-flow excitation function data indicate a gradual hardening of hot and dense nuclear matter as its density and temperature increase in reactions with higher beam energies.

nuclear matter↗

NASA NTRS1979

Laser beam propagation in atmospheric turbulence

The optical effects of atmospheric turbulence on the propagation of low power laser beams are reviewed in this paper. The optical effects are produced by the temperature fluctuations which result in fluctuations of the refractive index of air. The commonly-used models of index-of-refraction fluctuations are presented. Laser beams experience fluctuations of beam size, beam position, and intensity distribution within the beam due to refractive turbulence. Some of the observed effects are qualitatively explained by treating the turbulent atmosphere as a collection of moving gaseous lenses of various sizes. Analytical results and experimental verifications of the variance, covariance and probability distribution of intensity fluctuations in weak turbulence are presented. For stronger turbulence, a saturation of the optical scintillations is observed. The saturation of scintillations involves a progressive break-up of the beam into multiple patches; the beam loses some of its lateral coherence. Heterodyne systems operating in a turbulent atmosphere experience a loss of heterodyne signal due to the destruction of coherence.

Murty, S. S. R.↗