Search NASA⌕ Search

SEARCH · Search NASA

Results for “random variables”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Higher-Order Statistical Moments of Predicted Sonic Boom Waveforms Through Turbulence

Sonic booms generated by supersonic aircraft are affected by turbulence in the atmospheric boundary layer through which they propagate. Turbulence effects lead to random variability of the sonic boom waveforms measured on the ground, complicating the prediction of such waveforms. As an initial effort to predict the waveform variability, the solution of the coherent or mean sonic boom waveform has previously been formulated by the author. The current paper extends the formulation to derive an expression of the second-order statistical moment necessary to calculate the variance of the spectral amplitudes. Since the derivation uses a full wave equation, results using the derived expression will be compared with those obtained using a parabolic approximation. In order to fully quantify the uncertainty of our predictions, the higher order statistical moments are also formulated and are shown to be approximately zero. Consequently, the probability density function (pdf) of the spectral amplitudes is predicted to be Gaussian. The formulation is further extended to determine the pdf of the loudness of sonic booms and quantify the uncertainties associated with the loudness prediction. Results from the formulation are compared with available flight test data.

sonic booms↗
NASA NTRSDate not supplied

Evaluation of Finite Impulse Response Filters for Turbulence Effects on Sonic Booms

Turbulence effects on sonic booms lead to random variability of sonic boom waveforms measured on the ground, complicating the prediction of such waveforms using sonic boom propagation codes that do not account for turbulence effects. The NASA PCBoom software is one such code shown to accurately predict sonic booms above the atmospheric boundary layer but not those on the ground. Analyses of measured sonic booms show that, on average, the turbulence effects result in sonic boom loudness reduction, on average, that increases with propagation distance and turbulence strength. To efficiently account for such effects, a signal processing-based approach has been developed at NASA using finite impulse response filters derived from predicted waveforms by solving a nonlinear parabolic equation. Analyses indicate that the mean loudness reduction obtained using the filters does not increase with propagation distance or with increasing turbulence strength, in disagreement with flight test data. An alternative approach to the filters is currently under development employing the multiple scattering theory (MST) of wave propagation. Available results from the physics-based MST approach are used to evaluate the predictive capability of the filters, suggesting that the filters can severely underestimate the prediction of the MST approach.

sonic booms↗
NASA NTRSDate not supplied

Evaluation of Finite Impulse Response Filters for Turbulence Effects on Sonic Booms

Turbulence effects on sonic booms lead to random variability of sonic boom waveforms measured on the ground, complicating the prediction of such waveforms using sonic boom propagation codes that do not account for turbulence effects. The NASA PCBoom software is one such code shown to accurately predict sonic booms above the atmospheric boundary layer but not those on the ground. Analyses of measured sonic booms show that, on average, the turbulence effects result in sonic boom loudness reduction, on average, that increases with propagation distance and turbulence strength. To efficiently account for such effects, a signal processing-based approach has been developed at NASA using finite impulse response filters derived from predicted waveforms by solving a nonlinear parabolic equation. Analyses indicate that the mean loudness reduction obtained using the filters does not increase with propagation distance or with increasing turbulence strength, in disagreement with flight test data. An alternative approach to the filters is currently under development employing the multiple scattering theory (MST) of wave propagation. Available results from the physics-based MST approach are used to evaluate the predictive capability of the filters, suggesting that the filters can severely underestimate the prediction of the MST approach.

sonic booms↗
DOE OSTI2013

Classifying with confidence from incomplete information.

For this paper, we consider the problem of classifying a test sample given incomplete information. This problem arises naturally when data about a test sample is collected over time, or when costs must be incurred to compute the classification features. For example, in a distributed sensor network only a fraction of the sensors may have reported measurements at a certain time, and additional time, power, and bandwidth is needed to collect the complete data to classify. A practical goal is to assign a class label as soon as enough data is available to make a good decision. We formalize this goal through the notion of reliability—the probability that a label assigned given incomplete data would be the same as the label assigned given the complete data, and we propose a method to classify incomplete data only if some reliability threshold is met. Our approach models the complete data as a random variable whose distribution is dependent on the current incomplete data and the (complete) training data. The method differs from standard imputation strategies in that our focus is on determining the reliability of the classification decision, rather than just the class label. We show that the method provides useful reliability estimates of the correctness of the imputed class labels on a set of experiments on time-series data sets, where the goal is to classify the time-series as early as possible while still guaranteeing that the reliability threshold is met.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2017

Delineating Facies Spatial Distribution by Integrating Ensemble Data Assimilation and Indicator Geostatistics with Level Set Transformation.

A new approach is developed to delineate the spatial distribution of discrete facies (geological units that have unique distributions of hydraulic, physical, and/or chemical properties) conditioned not only on direct data (measurements directly related to facies properties, e.g., grain size distribution obtained from borehole samples) but also on indirect data (observations indirectly related to facies distribution, e.g., hydraulic head and tracer concentration). Our method integrates for the first time ensemble data assimilation with traditional transition probability-based geostatistics. The concept of level set is introduced to build shape parameterization that allows transformation between discrete facies indicators and continuous random variables. The spatial structure of different facies is simulated by indicator models using conditioning points selected adaptively during the iterative process of data assimilation. To evaluate the new method, a two-dimensional semi-synthetic example is designed to estimate the spatial distribution and permeability of two distinct facies from transient head data induced by pumping tests. The example demonstrates that our new method adequately captures the spatial pattern of facies distribution by imposing spatial continuity through conditioning points. The new method also reproduces the overall response in hydraulic head field with better accuracy compared to data assimilation with no constraints on spatial continuity on facies.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Karhunen-Loeve expansion analysis of uncertainties in cloud microphysical property retrievals

Here, this study proposes a methodology of quantifying uncertainties for cloud retrievals on model resolution to facilitate the comparison with model outputs. Primary component analysis is applied to reduce the dimension of random variables (up to a factor of 50) and reveal the cross correlations in the input data, making large sampling computationally feasible and uncertainty quantification accurate and reliable. Our approach has the capability of parameterizing input uncertainties and attributing the uncertainties in the retrieval output to each individual source, which allows sensitivity analysis of cloud retrieval algorithms and provides directions for improving observation instruments as well as strategies. We applied the method to characterize uncertainties in cloud ice water content (IWC) retrieved from the DOE Atmospheric Radiation Measurement (ARM) programs baseline cloud microphysical retrieval algorithm (MICROBASE). We test it with a selected ice cloud case observed on 9 March 2000 at the ARM Southern Great Plains site during its 2000 cloud intensive observing period. The test results indicate that (1) uncertainties in the output retrieved by MICROBASE are comparable amongst different retrievals; (2) The mean values obtained by our UQ method are closer to the aircraft data with less errors compared to the direct ensemble average; (3) Ice water path (IWP) generally incurred larger uncertainty in optically thin ice clouds and there was more variability in vertical in the retrieved IWC; and (4) Uncertainties in the output are mainly due to the interactions among different modes of ARM radar pro files.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Bayesian State-Space Modeling Framework for Understanding and Predicting Golden Eagle Movements Using Telemetry Data

Predicting raptor movements through a wind power plant under given atmospheric and topographical conditions is a crucial first step in the overall goal of quantifying the risk of turbine-related collisions and mortalities. Extracting behavioral traits of golden eagles (Aquila chrysaetos) from telemetry data requires the fusion of noisy and sparse movement data (location, heading, velocity) with a stochastic mathematical representation of the eagles' decision-making processes. In this study, we framed this problem in a Bayesian state-space framework where both observations and decision-making are assumed to be stochastic processes connected through hidden states (mode of flight, intent), and the unknown model parameters are assumed to be random variables that are calibrated using the available telemetry data. This framework allowed for rigorous consideration of underlying uncertainties while allowing for both data and prior biological knowledge to contribute to a probabilistic and predictive agent-based movement model. We implemented and applied the Bayesian framework to understand movement behavior of 23 GPS-tagged golden eagles travelling in the western US for years 2019 and 2020. Our preliminary findings show that the Bayesian state-space framework provides a robust inverse modeling apparatus to decode eagle behavioral characteristics from telemetry data. This study was primarily aimed at verifying and validating the framework with selected golden eagle tracks (both long- and short-ranged), with future research aimed at extending the framework to include multi-mode flight, consideration of atmospheric data and uplift mechanisms, eagle-to-eagle interaction, and eagle-to-turbine interaction.

Bayesian modeling↗
DOE OSTI2018

Hanford Direct Feed Low Activity Waste Program Throughput Optimization Through Modeling - 18587

The Direct Feed Low Activity Waste (DFLAW) program at Hanford allows for the early treatment of tank wastes years earlier than previously planned. Operational Research (OR) modeling is being used to provide invaluable insight into how complex and integrated processes will operate. The DFLAW OR models are used to help optimize the design and operation of the integrated system. The DFLAW OR models include the processes and real-life complexities involved in the day-to-day operation of the systems and equipment, and can predict key performance characteristics including throughput, mission timescale, equipment and labor utilization, lag storage and capacity requirements. The DFLAW OR models also have the in-built capability to simulate random events such as equipment failures. Adding random variability to models introduces a real-life element to model results. For example, OR modeling breakdown analysis of equipment failure and restoration times may inform specific improvements to a design and plant operation, including areas where the provision of spare (duty/standby) equipment would prove beneficial; the optimum number and capacity of lag storage tanks (for waste intermediate storage, reagents, etc.); operational issues, including scheduling constraints, labor resource constraints, effect of shift schedules (e.g. days only, extended shifts, 24/7),, maintenance strategies, including planned maintenance, equipment upgrades, on-hand spare part stockholding, etc. At Hanford, standalone OR models to date, have been used to help identify key operational issues and improve the design and operation of individual systems within the Tank Farms, Low Activity Waste Pretreatment System (LAWPS), Low Activity Waste (LAW) Vitrification and Effluent Management Facility (EMF) and Integrated Disposal Facility (IDF), logistics and infrastructure, 222-S laboratory [References 1 and 2] including: - Identification of critical spares including transfer pumps; - Underpin size and number of lag storage tanks to maintain continuous LAW feed to the LAW melters; - Underpin sample turnaround times at the 222-S Laboratory (e.g. feed prequalification and evaporator campaign samples) and determine how targets can be achieved without impacting waste feed delivery to LAWPS and the LAW Vitrification Facility; - Determine operational strategy required at IDF including shift schedules and sizes and the number of ILAW transporters that can support DFLAW throughput; - Identification of lag storage management strategies for glass formers, reagents, other consumables due to delivery distance and delays due to bad weather; and - Truck delivery badging needs to support 24/7 operation and management of shift in/out times required to reduce congestion. An Integrated DFLAW OR model has been developed and combines the standalone OR models such that integrated throughput and availability can be assessed from a program perspective. The Integrated OR model will provide a more detailed understanding of the interactions between Tank Farms, LAWPS, LAW VIT and IDF such that key bottleneck areas can be identified and ensure improvements are focused. Examples of future improvements that may be considered include identification of critical spare parts, required on-site maintenance capabilities, shift scheduling, and craft and labor availability. The methods, results and insights derived by DFLAW OR modeling efforts are expected to provide useful applications to other projects that involve combinations of mechanical handling, chemical processing, waste treatment/packaging or environmental remediation operations. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

Unbiased Ratio Estimator for k-Eigenvalue in Monte Carlo Criticality Calculations

In criticality calculations, the effective multiplication factor (k{sub eff}) is an important parameter for reactor analysis. It can be interpreted as the ratio of the neutron production rate to the loss rate in a reactor core. In Monte Carlo (MC) criticality calculations, k{sub eff} is obtained for each active cycle i as: k-i = P{sub i}/L{sub i} for I = 1 to I where P{sub i} and L{sub i} are the production rate and the loss rate at cycle i, respectively, and they are normalized by the nominal source size for cycle. Since all neutron particles in a cycle will be absorbed or leaks out at the end, the loss rate, L{sub i}, can be set as the source rate, which is the unity. Thus, in MCNP5, k{sub i} is simply calculated with P{sub i} by using three conventional estimators; collision, absorption, and track-length estimators. In domain decomposed MC criticality calculations based on the fission and surface source (FSS) iteration method, k{sub i} is calculated by a ratio of P{sub i} and L{sub i}, where both P{sub i} and L{sub i} are random variables, sampled by using the same particle histories. When the numerator and the denominator are correlated, the variance of the ratio can be smaller than that of the numerator. Then, the final estimate of k{sub eff} is provided as the sample mean of k{sub i}, which will be referred as 'simple' ratio estimator in this paper. In Ref., it is shown that this ratio estimator converges faster and fluctuates less than the conventional estimator. A similar ratio estimator is also used in Serpent 2, as an implicit estimate of k{sub eff}. However, this type of 'simple' ratio estimator is biased, which means that the expectation of the ratio is different from the ratio of the expectation of the numerator and that of the denominator. In this paper, an unbiased ratio estimator which was originally derived in the cytometry analysis is introduced to estimate k{sub eff} and is tested on a continuous-energy UO{sub 2} pin-cell problem. In this test problem, it is shown that the variance of the unbiased ratio estimator is smaller than those of the conventional estimators provided by MCNP5, while the simple ratio estimator can be significantly biased when the number of histories per cycle is too small. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Fault Graph (fg)

This software package offers a novel fault graph modeling language and a simulation tool for executing models specified in the language. The modeling language has three primary features that distinguish it from similar reliability analysis tools. These are (1) a random variable modeling several distinct outcomes of a single fault; (2) chains of faults in which one fault triggers another; and (3) time to fail sampled from probability distributions including positive normal, exponential, Weibull with a minimum, or immediate.

Nutaro, James↗
DOE OSTI2013

IS QUASAR OPTICAL VARIABILITY A DAMPED RANDOM WALK?

The damped random walk (DRW) model is increasingly used to model the variability in quasar optical light curves, but it is still uncertain whether the DRW model provides an adequate description of quasar optical variability across all timescales. Using a sample of OGLE quasar light curves, we consider four modifications to the DRW model by introducing additional parameters into the covariance function to search for deviations from the DRW model on both short and long timescales. We find good agreement with the DRW model on timescales that are well sampled by the data (from a month to a few years), possibly with some intrinsic scatter in the additional parameters, but this conclusion depends on the statistical test employed and is sensitive to whether the estimates of the photometric errors are correct to within {approx}10%. On very short timescales (below a few months), we see some evidence of the existence of a cutoff below which the correlation is stronger than the DRW model, echoing the recent finding of Mushotzky et al. using quasar light curves from Kepler. On very long timescales (>a few years), the light curves do not constrain models well, but are consistent with the DRW model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Impedance Response Influenced by Variability in the Random Distribution of Physical Properties of Coated Materials in Two-Dimensional Space

Heterogeneous physical characteristics of a system featuring a single-layer film on a metallic surface have been explored via its impedance response. The Nyquist plot showed a distorted semicircle, indicative of the system's unique distribution characteristics. Utilizing a copula-based probability method, a two-dimensional deterministic impedance model was successfully integrated, accounting for spatial physical properties such as permittivity and electrical conductivity. This strategy enabled in-depth exploration and mechanistic quantification of a broad spectrum of properties. A quantitative understanding of impedance signal alterations, characterized by normally or log-normally correlated variables, was achieved through the variation in aspect ratio and characteristic frequency of the impedance spectra. Log-normally distributed electrical properties provided a superior representation of the distorted impedance spectra. As coefficient of variation (CV) values fluctuated, the aspect ratio and characteristic frequency showed heightened sensitivity to log-normal permittivity compared to log-normal electrical conductivity. Notably, a marked positive linear correlation between electrical properties resulted in an impedance response that approximated perfect semicircular spectra. The variability in the electrical properties' distribution was demonstrated by considering the correlation coefficient between electrical conductivity and the z-direction position. Furthermore, the highest aspect ratio of the impedance spectra was observed when the electrical conductivity was randomly distributed across the z-direction space.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Neural chaos: A spectral stochastic neural operator

Building surrogate models for operators with uncertainty quantification capabilities is essential for many engineering applications where randomness–such as variability in material properties, boundary conditions, and initial conditions–is unavoidable. Polynomial Chaos Expansion (PCE) is widely recognized as a go-to method for constructing stochastic surrogates in both intrusive and non-intrusive ways, and it has recently been used in the context of operator learning. However, its application becomes challenging for complex or high-dimensional processes, as achieving accuracy requires higher-order polynomials, which can increase computational demand and/or the risk of overfitting. Furthermore, PCE requires specialized treatments to manage random variables that are not independent, and these treatments may be problem-dependent or may fail with increasing complexity. Here, in this work, we adopt the same formalism as the spectral expansion used in PCE; however, we replace the classical polynomial basis functions with neural network (NN) basis functions to leverage their expressivity. To achieve this, we propose an algorithm that identifies NN-parameterized basis functions in a purely data-driven manner, without any prior assumptions about the joint distribution of the random variables involved, whether independent or dependent, or about their marginal distributions. The proposed algorithm identifies each NN-parameterized basis function sequentially, ensuring they are orthogonal with respect to the data distribution. The basis functions are constructed directly on the joint stochastic variables without requiring a tensor product structure or assuming independence of the random variables. This approach may offer greater flexibility for complex stochastic models, while simplifying implementation compared to the tensor product structures typically used in PCE to handle random vectors. This is particularly advantageous given the current state of open-source packages, where building and training neural networks can be done with just a few lines of code and extensive community support. We demonstrate the effectiveness of the proposed scheme through several numerical examples of varying complexity and provide comparisons with classical PCE.

Polynomial chaos expansion↗
DOE OSTI2015

Verifying and Validating Simulation Models

This presentation is a high-level discussion of the Verification and Validation (V&V) of computational models. Definitions of V&V are given to emphasize that “validation” is never performed in a vacuum; it accounts, instead, for the current state-of-knowledge in the discipline considered. In particular comparisons between physical measurements and numerical predictions should account for their respective sources of uncertainty. The differences between error (bias), aleatoric uncertainty (randomness) and epistemic uncertainty (ignorance, lack-of- knowledge) are briefly discussed. Four types of uncertainty in physics and engineering are discussed: 1) experimental variability, 2) variability and randomness, 3) numerical uncertainty and 4) model-form uncertainty. Statistical sampling methods are available to propagate, and analyze, variability and randomness. Numerical uncertainty originates from the truncation error introduced by the discretization of partial differential equations in time and space. Model-form uncertainty is introduced by assumptions often formulated to render a complex problem more tractable and amenable to modeling and simulation. The discussion concludes with high-level guidance to assess the “credibility” of numerical simulations, which stems from the level of rigor with which these various sources of uncertainty are assessed and quantified.

Mathematics & Computing(97)↗
NASA NTRS1985

The effect of material heterogeneity and random loading on the mechanics of fatigue crack growth

This paper reviews experimental work on the influence of variable amplitude or random loads on the mechanics and micromechanisms of fatigue crack growth. Implications are discussed in terms of the crack driving force, local plasticity, crack closure, crack blunting, and microstructure. Due to heterogeneity in the material's microstructure, the crack growth rate varies with crack tip position. Using the weakest link theory, an expression for crack growth rate is obtained as the expectation of a random variable. This expression is used to predict the crack growth rates for aluminum alloys, a titanium alloy, and a nickel steel in the mid-range region. It is observed, using the present theory, that the crack growth rate obeys the power law for small stress intensity factor range, and that the power is a function of a material constant.

Srivatsan, T. S.↗