SEARCH · Search NASA

Results for “statistical models”

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 1,315 records · Page 73

All sources NASA NTRS DOE OSTI

NASA NTRS2002

Scale Dependence of Statistics of Spatially Averaged Rain Rate Seen in TOGA COARE Comparison with Predictions from a Stochastic Model

A characteristic feature of rainfall statistics is that they in general depend on the space and time scales over which rain data are averaged. As a part of an earlier effort to determine the sampling error of satellite rain averages, a space-time model of rainfall statistics was developed to describe the statistics of gridded rain observed in GATE. The model allows one to compute the second moment statistics of space- and time-averaged rain rate which can be fitted to satellite or rain gauge data to determine the four model parameters appearing in the precipitation spectrum - an overall strength parameter, a characteristic length separating the long and short wavelength regimes and a characteristic relaxation time for decay of the autocorrelation of the instantaneous local rain rate and a certain 'fractal' power law exponent. For area-averaged instantaneous rain rate, this exponent governs the power law dependence of these statistics on the averaging length scale $L$ predicted by the model in the limit of small $L$. In particular, the variance of rain rate averaged over an $L \times L$ area exhibits a power law singularity as $L \rightarrow 0$. In the present work the model is used to investigate how the statistics of area-averaged rain rate over the tropical Western Pacific measured with ship borne radar during TOGA COARE (Tropical Ocean Global Atmosphere Coupled Ocean Atmospheric Response Experiment) and gridded on a 2 km grid depends on the size of the spatial averaging scale. Good agreement is found between the data and predictions from the model over a wide range of averaging length scales.

Kundu, Prasun K.↗

DOE OSTI2020

Light element ($$Z=1,2$$) production from spontaneous ternary fission of $$^{252}$$Cf

Abstract The yields of light elements ( $$Z=1,2$$ Z = 1 , 2 ) obtained from spontaneous ternary fission of $$^{252}$$ 252 Cf are treated within a nonequilibrium approach, and the contribution of unstable nuclei and excited bound states is taken into account. These light cluster yields may be used to probe dense matter, and to infer in-medium corrections such as Pauli blocking which is determined by the nucleon density. Continuum correlations are calculated from scattering phase shifts using the Beth-Uhlenbeck formula, and the effect of medium modification is estimated. The relevant distribution is reconstructed from the measured yields of isotopes. This describes the state of the nucleon system at scission and cluster formation, using only three Lagrange parameters which are the nonequilibrium counterparts of the temperature and chemical potentials, as defined in thermodynamic equilibrium. We concluded that a simple nuclear statistical equilibrium model neglecting continuum correlations and medium effects is not able to describe the measured distribution of H and He isotopes. Moreover, the freeze-out concept may serve as an important ingredient to the nonequilibrium approach using the relevant statistical operator concept.

Röpke, G.↗

NASA NTRS2023

Spacesuit Fit and Mobility Assessments by Digital Human Modeling

Spacesuit Fit and Mobility Assessments by Digital Human Modeling K. Han Kim (Leidos, Inc.) Elizabeth A. Benson (KBR, Inc.) Sudhakar L. Rajulu (NASA Johnson Space Center) Spacesuits are required to accommodate safe operations for astronauts across gender and a wide variety of body shapes and sizes. This goal has been of particular importance given the increasing diversity of NASA crewmembers for upcoming Missions. While testing with design prototypes is a critical step for spacesuit development, iterative mockup design, fabrication, and human subject tests can be extremely costly and time consuming. Moreover, testing with a limited subject pool has often raised questions for validity, as test subjects need to represent the entire astronaut population, not only of the past or current, but also the future. This study is aimed at demonstrating how digital human modeling (DHM) tools have been built and directly supported NASA spacesuit developments. With DHM, the computer aided design (CAD) model of a spacesuit was integrated with human body models. Two use scenarios are presented in this paper, namely fit and mobility. For fit assessments, the suit-to-body contact and compression patterns were estimated using 3D human body scan models virtually wearing a spacesuit model. A statistical fit classifier was made from the contact patterns and applied against a large database of body scans (N=2,500). With this technique, the NASA reference design spacesuit Exploration EVA Mobility Unit (xEMU)was verified to accommodate 90% of the astronaut-like population, with the critical dimensions covering 1stto 99thpercentiles of the target body measurements. A similar technique assessed XEMU mobility. The maximum reach envelopes were considered, within which the work objects and critical hardware interfaces should be located for safe and ergonomic operations. While the reach envelope geometry varies significantly with the suit wearer’s body size and strength, the existing test data did not include the subjects critically required to define suit mobility requirements, such as very small females or large males. Using the xEMU virtual model kinematically simulated and permuted for a hypothetical wearer, however, the existing data were statistically transformed and scaled. This method enabled for a parametric estimation of the reach envelopes from the 1stpercentile female or 99thpercentile male. The outcome was successfully incorporated for requirement developments. With the new DHM tools, human integration of the spacesuit was structurally simulated and predictively assessed, which would have been otherwise impossible. Also, the needs for iterative mockup and subject tests were significantly reduced, which resulted in time and cost saving. Additional work is in progress to integrate additional vehicles, tools, and hardware with DHM framework.

K Han Kim↗

DOE OSTI2020

Propagation of statistical uncertainties of Skyrme mass models to simulations of r -process nucleosynthesis

Uncertainties in nuclear models have a major impact on simulations that aim at understanding the origin of heavy elements in the universe through the rapid neutron capture process ( r process) of nucleosynthesis. Within the framework of the nuclear density functional theory, we use results of Bayesian statistical analysis to propagate uncertainties in the parameters of energy density functionals to the predicted r -process abundance pattern, by way not only of the nuclear masses but also through the influence of the masses on β -decay and neutron capture rates. We point out the importance of the nonequilibrium end stage of the r process in determining the width of the resulting abundance pattern uncertainty bands. We additionally make the first identifications of specific parameters of Skyrme-like energy density functionals which show tentative correlations with particular aspects of the r -process abundance pattern. While previous studies have explored the reduction in the abundance pattern uncertainties due to anticipated new measurements of neutron-rich nuclei, here we point out that an even larger reduction will occur when these new measurements are used to reduce the uncertainty of model predictions of masses, which are then propagated through to the abundance pattern. We make a quantitative prediction for how large this reduction will be.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2021

Windowed multipole representation of R -matrix cross sections

Nuclear cross sections are basic inputs to any nuclear computation. Campaigns of experiments are fitted with the parametric R-matrix model of quantum nuclear interactions, and the resulting cross sections are documented—both pointwise and as resonance parameters (with uncertainties)—in standard evaluated nuclear data libraries (ENDF, JEFF, BROND, JENDL, CENDL, TENDL): these constitute our common knowledge of fundamental low-energy nuclear cross sections. In the past decade, a collaborative effort has been deployed to establish a new nuclear cross-section library format—the Windowed Multipole Library—with the goal of considerably reducing the computational cost of cross-section calculations in nuclear transport simulations. This work lays the theoretical foundations underpinning these efforts. From general R-matrix scattering theory, we derive the windowed multipole representation of nuclear cross sections. Though physically and mathematically equivalent to R-matrix cross sections, the windowed multipole representation is particularly well suited for subsequent temperature treatment of angle-integrated cross sections, in particular Doppler broadening, which is the averaging of cross sections over the thermal motion of the target atoms. Doppler broadening is of critical importance in neutron transport applications, as it ensures the stability of many nuclear reactors (negative thermal reactivity). Yet, Doppler broadening of nuclear cross sections has been a considerable bottleneck for nuclear transport computations, often requiring memory-costly pretabulations. We show that the windowed multipole representation can perform accurate Doppler broadening analytically (up to the first reaction threshold), from which we derive cross-section temperature derivatives to any order—all computable on the fly (without precalculations stored in memory). Furthermore, we here establish a way of converting the R-matrix resonance parameters uncertainty (covariance matrices) into windowed multipole parameters uncertainty. We show that generating stochastic nuclear cross sections by sampling from the resulting windowed multipole covariance matrix can reproduce the cross-section uncertainty in the original nuclear data file. The windowed multipole representation is therefore a novel nuclear physics formalism able to generate Doppler broadened stochastic nuclear cross sections on the fly, unlocking breakthrough computational gains for nuclear computations. Through this foundational paper, we hope to make the windowed multipole representation accessible, reproducible, and usable for the nuclear physics community, as well as provide the theoretical basis for future research on expanding its capabilities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2025

Challenges in quantifying unparameterized spatial uncertainties in deep geologic repositories for nuclear waste

Spatially heterogeneous uncertainties are prevalent in geophysical modeling applications, such as probabilistic post-closure performance assessment (PA) of deep geologic repositories for nuclear waste. Such uncertainties are often highly influential to model outputs, so it is desirable to identify the most important mechanisms by which they influence model predictions. However, these uncertainties are often unparameterized in the sense that there is no set of parameters that can be specified to yield a particular realization of the uncertainty. Additionally, the uncertainty is not intrinsically endowed with a parameterization that captures a realization’s mechanistic influence on model outputs. Therefore, in this work we present a novel methodology to develop and assess a set of proxy variables that aim to represent this influence. We show how they can be computed, downselected, and incorporated into the construction of statistical surrogate models mapping model inputs to outputs. We present our methodology in the context of a motivating application problem in deep geologic repository PA and discuss the challenges in capturing the effects of these spatial heterogeneities in uncertainty analyses.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

NASA NTRS1988

Fluctuation statistics of outgoing longwave radiation in a general circulation model and in satellite data

Fluctuation statistics of Northern-Hemisphere-winter outgoing longwave radiation (OLR) obtained from data collected by Nimbus 7 and NOAA AVHRR are compared with those computed by a general circulation model (GCM). It was found that the GCM OLR has more variance than does the satellite, but that the GCM clear-sky OLR has less variance than the satellite OLR. The variance of the GCM outgoing full-sky OLR was found to be forced by clouds. While the standard deviation of the GCM full OLR is larger than that of the satellite OLR, the geographical patterns of both are similar. The autocorrelations of the modeled and the observed OLR also display similar geographical distributions, and they are in turn correlated with the dynamics.

Charlock, Thomas P.↗

DOE OSTI2022

Are Terrestrial Biosphere Models Fit for Simulating the Global Land Carbon Sink?

The Global Carbon Project estimates that the terrestrial biosphere has absorbed about one-third of anthropogenic CO 2 emissions during the 1959–2019 period. This sink-estimate is produced by an ensemble of terrestrial biosphere models and is consistent with the land uptake inferred from the residual of emissions and ocean uptake. The purpose of our study is to understand how well terrestrial biosphere models reproduce the processes that drive the terrestrial carbon sink. One challenge is to decide what level of agreement between model output and observation-based reference data is adequate considering that reference data are prone to uncertainties. To define such a level of agreement, we compute benchmark scores that quantify the similarity between independently derived reference data sets using multiple statistical metrics. Models are considered to perform well if their model scores reach benchmark scores. Our results show that reference data can differ considerably, causing benchmark scores to be low. Model scores are often of similar magnitude as benchmark scores, implying that model performance is reasonable given how different reference data are. While model performance is encouraging, ample potential for improvements remains, including a reduction in a positive leaf area index bias, improved representations of processes that govern soil organic carbon in high latitudes, and an assessment of causes that drive the inter-model spread of gross primary productivity in boreal regions and humid tropics. The success of future model development will increasingly depend on our capacity to reduce and account for observational uncertainties.

58 GEOSCIENCES↗

NASA NTRS2022

New Directions in Tropical Phenology

Earth’s most speciose biomes are in the tropics, yet tropical plant phenology remains poorly understood. Tropical phenological data are comparatively scarce and viewed through the lens of a ‘temperate phenological paradigm’ expecting phenological traits to respond to strong, predictably annual shifts in climate (e.g., between subfreezing and frost-free periods). Digitized herbarium data greatly expand existing phenological data for tropical plants; and circular data, statistics, and models are more appropriate for analyzing tropical (and temperate) phenological datasets. Phylogenetic information, which remains seldom applied in phenological investigations, provides new insights into phenological responses of large groups of related species to climate. Consistent combined use of herbarium data, circular statistical distributions, and robust phylogenies will rapidly advance our understanding of tropical – and temperate – phenology.

tropical phenology↗

DOE OSTI2024

Bayesian optimization algorithms for accelerator physics

Accelerator physics relies on numerical algorithms to solve optimization problems in online accelerator control and tasks such as experimental design and model calibration in simulations. The effectiveness of optimization algorithms in discovering ideal solutions for complex challenges with limited resources often determines the problem complexity these methods can address. The accelerator physics community has recognized the advantages of Bayesian optimization algorithms, which leverage statistical surrogate models of objective functions to effectively address complex optimization challenges, especially in the presence of noise during accelerator operation and in resource-intensive physics simulations. In this review article, we offer a conceptual overview of applying Bayesian optimization techniques toward solving optimization problems in accelerator physics. We begin by providing a straightforward explanation of the essential components that make up Bayesian optimization techniques. We then give an overview of current and previous work applying and modifying these techniques to solve accelerator physics challenges. Finally, we explore practical implementation strategies for Bayesian optimization algorithms to maximize their performance, enabling users to effectively address complex optimization challenges in real-time beam control and accelerator design. Published by the American Physical Society 2024

43 PARTICLE ACCELERATORS↗

DOE OSTI2022

Modeling Flow and Transport in Fractured Rock using Machine Learning [Slides]

Fractured systems in the subsurface play a role in many natural and engineered applications such as geologic carbon sequestration, hydraulic fracturing and underground nuclear test detection. Structural information (fracture size, orientation, etc.) plays a key role in governing the dominant physics for these systems but can only be known statistically. Traditional modeling approaches either ignore or idealize structural information at these larger scales because we lack a computational framework that utilizes it in its entirety. The work presented here integrates computational physics, machine learning and graph theory to make a paradigm shift from computationally intensive high-fidelity models to coarse-scale graphs without loss of critical structural information. We exploit the underlying discrete structure of fracture networks in systems considering flow through fractures and fracture propagation. We demonstrate that compact graph representations require significantly fewer degrees of freedom (dof) to capture micro-fracture information and further accelerate these models with Machine Learning.

58 GEOSCIENCES↗

NASA NTRS1975

A linear feedback guidance law for an Aeromaneuvering Orbit-to-Orbit Shuttle

The sensitivity of the trajectory of an Aeromaneuvering Orbit-to-Orbit Shuttle (AMOOS) to unpredictable variations in atmospheric density is shown and a compensating guidance law is developed. The sensitivity was determined using a varying atmosphere model whose statistical properties matched those of actual atmosphere data. This study showed that a compensating guidance law is crucial to the feasibility of AMOOS. Using an optimal regulator approach, a linear feedback law for bank angle was developed and checked out with the atmosphere model. This law uses the altitude profile of the trajectory to control the amount of energy dissipated within the atmosphere.

Rehder, J. J.↗

NASA NTRS1991

Dynamical evolution of topology of large-scale structure

The nonlinear effects of statistical biasing and gravitational evolution on the genus are studied. The biased galaxy subset is picked for the first time by actually identifying galaxy-sized peaks above a fixed threshold in the initial conditions, and their subsequent evolution is followed. It is found that in the standard cold dark matter (CDM) model the statistical biasing in the locations of galaxies produces asymmetry in the genus curve and coupling with gravitational evolution gives rise to a shift in the genus curve to the left in moderately nonlinear regimes. Gravitational evolution alone reduces the amplitude of the genus curve due to strong phase correlations in the density field and also produces asymmetry in the curve. Results on the genus of the mass density field for both CDM and hot dark matter models are consistent with previous work by Melott, Weinberg, and Gott (1987).

Park, Changbom↗

DOE OSTI2025

Rare events and Griffiths phases in topological quantum error correction

The performance of quantum error correcting (QEC) codes is often studied under the assumption of spatiotemporally uniform error rates. On the other hand, experimental implementations almost always produce heterogeneous error rates, in either space or time, as a result of effects such as imperfect fabrication and/or cosmic rays. It is therefore important to understand if and how their presence can affect the performance of QEC in qualitative ways. Here, in this work, we study the effects of nonuniform error rates in the representative examples of the 1D repetition code and the 2D toric code, focusing on when they have extended spatiotemporal correlations; these may arise, for instance, from rare events (such as cosmic rays) that temporarily elevate error rates over the entire code patch. These effects can be described in the corresponding statistical mechanics models for decoding, where long-range correlations in the error rates lead to extended rare regions of weaker coupling. For the 1D repetition code where the rare regions are linear, we find two distinct decodable phases: a conventional ordered phase in which logical failure rates decay exponentially with the code distance, and a rare-region dominated Griffiths phase in which failure rates are parametrically larger and decay as a stretched exponential. In particular, the latter phase is present when the error rates in the rare regions are above the bulk threshold. For the 2D toric code where the rare regions are planar, we find no decodable Griffiths phase: rare events which boost error rates above the bulk threshold lead to an asymptotic loss of threshold and failure to decode. Unpacking the failure mechanism implies that techniques for suppressing extended sequences of repeated rare events (which, without intervention, will be statistically present with high probability) will be crucial for QEC with the toric code.

classical statistical mechanics↗

DOE OSTI2025

On the Quantum Mechanics of Entropic Forces

It was conjectured 30 years ago that gravity could arise from the entropic rearrangement of information. We offer a set of microscopic quantum models which realize this idea in detail. In particular, we suggest a simple mechanism by which Newton’s law of gravity arises from extremization of the free energy of a collection of qubits or oscillators, rather than from the exchange of virtual quanta of a fundamental field. We give both a local and a nonlocal version of the construction and show how to distinguish a range of these entropic models from ordinary perturbative quantum gravity using existing observations and near-term experiments.

quantum fluctuations↗

DOE OSTI2023

Developing and applying quantifiable metrics for diagnostic and experiment design on Z

This project applies methods in Bayesian inference and modern statistical methods to quantify the value of new experimental data, in the form of new or modified diagnostic configurations and/or experiment designs. We demonstrate experiment design methods that can be used to identify the highest priority diagnostic improvements or experimental data to obtain in order to reduce uncertainties on critical inferred experimental quantities and select the best course of action to distinguish between competing physical models. Bayesian statistics and information theory provide the foundation for developing the necessary metrics, using two high impact experimental platforms on Z as exemplars to develop and illustrate the technique. We emphasize that the general methodology is extensible to new diagnostics (provided synthetic models are available), as well as additional platforms. We also discuss initial scoping of additional applications that began development in the last year of this LDRD.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2023

Optimization of the generator coordinate method with machine-learning techniques for nuclear spectra and neutrinoless double- β decay: Ridge regression for nuclei with axial deformation

The generator coordinate method (GCM) is an important tool of choice for modeling large-amplitude collective motion in atomic nuclei. The computational complexity of the GCM increases rapidly with the number of collective coordinates. It imposes a strong restriction on the applicability of the method. In this work, we propose a subspace-reduction algorithm that employs optimal statistical ML models as surrogates for exact quantum-number projection calculations for norm and Hamiltonian kernels. The model space of the original GCM is reduced to a subspace relevant for nuclear low energy spectra and the NME of ground state to ground state 0νββ decay based on the orthogonality condition (OC) and the energy-transition-orthogonality procedure (ENTROP), respectively. For simplicity, the polynomial ridge regression (RR) algorithm is used to learn the norm and Hamiltonian kernels of axially deformed configurations. The efficiency and accuracy of this algorithm are illustrated for 76 Ge and 76 Se by comparing results obtained using the optimal RR models to direct GCM calculations. The low-lying energy spectra of 76 Ge and 76 Se, as well as the 0νββ-decay NME between their ground states, are computed. Furthermore, the results show that the performance of the GCM+OC/ENTROP+RR is more robust than that of the GCM+RR alone, and the former can reproduce the results of the original GCM calculation accurately with a significantly reduced computational cost.

59 ≤ A ≤ 89↗

NASA NTRS1989

The discrimination of sea ice types using SAR backscatter statistics

X-band (HH) synthetic aperture radar (SAR) data of sea ice collected during the Marginal Ice Zone Experiment in March and April of 1987 was statistically analyzed with respect to discriminating open water, first-year ice, multiyear ice, and Odden. Odden are large expanses of nilas ice that rapidly form in the Greenland Sea and transform into pancake ice. A first-order statistical analysis indicated that mean versus variance can segment out open water and first-year ice, and skewness versus modified skewness can segment the Odden and multilayer categories. In additions to first-order statistics, a model has been generated for the distribution function of the SAR ice data. Segmentation of ice types was also attempted using textural measurements. In this case, the general co-occurency matrix was evaluated. The textural method did not generate better results than the first-order statistical approach.

Shuchman, Robert A.↗