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 1,009 records · Page 56
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

A mechanistic damage model for solid oxide fuel cell ceramic materials - Part I: Constitutive modeling

A multiscale mechanistic damage model for solid oxide fuel cell (SOFC) ceramics that combines micromechanics of stiffness reduction due to material porosity change and microcracking with a continuum damage mechanics (CDM) description for the evolution of microcracks up to fracture is developed in this work. The model also accounts for volumetric swelling of the anode due to redox cycling. Porosity change and swelling during redox combined with normal operating thermomechanical loads could increase the stress, strain and damage distributions in SOFC stacks leading them to failure. The proposed model involves three governing parameters: (i) the porosity in terms of the pore or void volume fraction, (ii) material swelling magnitude prescribed at a given loading step, and (iii) the damage variable that describes microcracking caused by thermomechanical loads and swelling. Pores and microcracks in the ceramics are modeled as randomly distributed ellipsoidal inclusions with negligible stiffness by an Eshelby-Mori-Tanaka formulation combined with an inclusion orientation distribution method. Microcracking damage evolution is described by a CDM formulation. Volumetric swelling is treated in a similar way to thermal expansion in the constitutive relations. After validation and numerical checks, the damage model is used to analyze a simplified SOFC positive-electrode/electrolyte/negative-electrode (PEN) structure subjected a redox cycle in addition to thermomechanical loads experienced during normal operation of a SOFC stack.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2021

Dependency of active pressure and equation of state on stiffness of wall

Autonomous motion and motility are hallmarks of active matter. Active agents, such as biological cells and synthetic colloidal particles, consume internal energy or extract energy from the environment to generate self-propulsion and locomotion. These systems are persistently out of equilibrium due to continuous energy consumption. It is known that pressure is not always a state function for generic active matter. Torque interaction between active constituents and confinement renders the pressure of the system a boundary-dependent property. The mechanical pressure of anisotropic active particles depends on their microscopic interactions with a solid wall. Using self-propelled dumbbells confined by solid walls as a model system, we perform numerical simulations to explore how variations in the wall stiffness influence the mechanical pressure of dry active matter. In contrast to previous findings, we find that mechanical pressure can be independent of the interaction of anisotropic active particles with walls, even in the presence of intrinsic torque interaction. Particularly, the dependency of pressure on the wall stiffness vanishes when the stiffness is above a critical level. In such a limit, the dynamics of dumbbells near the walls are randomized due to the large torque experienced by the dumbbells, leading to the recovery of pressure as a state variable of density.

42 ENGINEERING↗
NASA NTRS1989

Stochastic modeling of crack initiation and short-crack growth under creep and creep-fatigue conditions

A simplified stochastic model is proposed for crack initiation and short-crack growth under creep and creep-fatigue conditions. Material inhomogeneity provides the random nature of crack initiation and early growth. In the model, the influence of microstructure is introduced by the variability of: (1) damage accumulation along grain boundaries, (2) critical damage required for crack initiation or growth, and (3) the grain-boundary length. The probabilities of crack initiation and growth are derived by using convolution integrals. The model is calibrated and used to predict the crack density and crack-growth rate of short cracks of 304 stainless steel under creep and creep-fatigue conditions. The mean-crack initiation lives are predicted to be within an average deviation of about 10 percent from the experimental results. The predicted cumulative distributions of crack-growth rate follow the experimental data closely. The applicability of the simplified stochastic model is discussed and the future research direction is outlined.

Kitamura, Takayuki↗
NASA NTRS1992

Stochastic modeling of crack initiation and short-crack growth under creep and creep-fatigue conditions

A simplified stochastic model is proposed for crack initiation and short-crack growth under creep and creep-fatigue conditions. Material inhomogeneity provides the random nature of crack initiation and early growth. In the model, the influence of microstructure is introduced by the variability of: (1) damage accumulation along grain boundaries, (2) critical damage required for crack initiation or growth, and (3) the grain-boundary length. The probabilities of crack initiation and growth are derived by using convolution integrals. The model is calibrated and used to predict the crack density and crack-growth rate of short cracks of 304 stainless steel under creep and creep-fatigue conditions. The mean-crack initiation lives are predicted to be within an average deviation of about 10 percent from the experimental results. The predicted comulative distributions of crack-growth rate follow the experimental data closely. The applicability of the simplified stochastic model is discussed and the future research direction is outlined.

Kitamura, Takayuki↗
NASA NTRS2001

Spatial scale of motion segmentation from speed cues

For the accurate perception of multiple, potentially overlapping, surfaces or objects, the visual system must distinguish different local motion vectors and selectively integrate similar motion vectors over space to segment the retinal image properly. We recently showed that large differences in speed are required to yield a percept of motion transparency. In the present study, to investigate the spatial scale of motion segmentation from speed cues alone, we measured the speed-segmentation threshold (the minimum speed difference required for 75% performance accuracy) for 'corrugated' random-dot patterns, i.e. patterns in which dots with two different speeds were alternately placed in adjacent bars of variable width. In a first experiment, we found that, at large bar widths, a smaller speed difference was required to segment and perceive the corrugated pattern of moving dots, while at small bar-widths, a larger speed difference was required to segment the two speeds and perceive two transparent surfaces of moving dots. Both the perceptual and segmentation performance transitions occurred at a bar width of around 0.4 degrees. In a second experiment, speed-segmentation thresholds were found to increase sharply when dots with different speeds were paired within a local pooling area. The critical pairing distance was about 0.2 degrees in the fovea and increased linearly with stimulus eccentricity. However, across the range of eccentricities tested (up to 15 degrees ), the critical pairing distance did not change much and remained close to the receptive field size of neurons within the primate primary visual cortex. In a third experiment, increasing dot density changed the relationship between speed-segmentation thresholds and bar width. Thresholds decreased for large bar widths, but increased for small bar widths. All of these results are well fit by a simple stochastic model, which estimates the probabilities of having identical or different motion vectors within a local pooling area whose size is the same as that of primate V1 neurons. Altogether, these results demonstrate that speed-based segmentation can function well, even at small spatial scales (i.e. high-spatial frequencies of spatial corrugation) and thereby emphasizes the critical role of a local pooling process early in the cortical motion-processing pathway.

NASA Discipline Neuroscience↗
NASA NTRS2013

A Review of Spectral Methods for Variable Amplitude Fatigue Prediction and New Results

A comprehensive review of the available methods for estimating fatigue damage from variable amplitude loading is presented. The dependence of fatigue damage accumulation on power spectral density (psd) is investigated for random processes relevant to real structures such as in offshore or aerospace applications. Beginning with the Rayleigh (or narrow band) approximation, attempts at improved approximations or corrections to the Rayleigh approximation are examined by comparison to rainflow analysis of time histories simulated from psd functions representative of simple theoretical and real world applications. Spectral methods investigated include corrections by Wirsching and Light, Ortiz and Chen, the Dirlik formula, and the Single-Moment method, among other more recent proposed methods. Good agreement is obtained between the spectral methods and the time-domain rainflow identification for most cases, with some limitations. Guidelines are given for using the several spectral methods to increase confidence in the damage estimate.

Larsen, Curtis E.↗
DOE OSTI2015

Detectability of historical trends in station‐based precipitation characteristics over the continental United States

Abstract The goal of this paper is to detect secular trends in observed, station‐based precipitation variations and extreme event occurrences over the United States. By definition, detectable trends are those that are unlikely to have arisen from internal variability alone. To represent this internal variability, we use station‐specific, seasonally varying, daily time scale stationary stochastic weather models—models in which the simulated interannual‐to‐multidecadal precipitation variance is purely the result of the random evolution of daily precipitation events within a given time period—to first estimate the trends in various means and extremes that can occur even with fixed, climatological daily precipitation characteristics. Detection of secular trends in the observed climate—whether naturally or anthropogenically induced—can then be defined relative to this stochastic variability, i.e., as trends in the means and/or extremes that are likely to have occurred only through a change in the underlying precipitation characteristics. The derived results have two important ramifications. First, they identify “hot spot” regions in which trends in precipitation characteristics are already emerging from within the envelope of stochastic variability, including (but not limited to) positive trends in annual occurrence across most of the U.S. and positive trends in annual intensity and heavy‐event accumulations across the Interior Plains and around the Great Lakes. Further, they identify “sentinel” metrics which show the greatest detectability—e.g., annual precipitation occurrence and intensity—as well as those which show the least detectability and hence are unlikely to produce any statistically meaningful signals—e.g., seasonal total precipitation and extremes.

Anderson, Bruce T.↗
DOE OSTI2021

Learning-based CO 2 concentration prediction: Application to indoor air quality control using demand-controlled ventilation

There have been increasing concerns over the air quality inside buildings as high levels of bio-effluents can cause nausea, dizziness, headaches, and fatigue to the people working in those spaces. First published in 2004 as Standard 62.1, ASHRAE Standard 62.2-2019 requires highly occupied spaces to implement heating, ventilation, and air conditioning (HVAC) that can dilute contaminants produced by occupants. In this regard, occupant-centric ventilation control has been regarded as an effective practice to maintain a satisfactory indoor air quality (IAQ) when dealing with highly variable occupancy environments. However, few established models in current literature and practice consider dynamic occupancy behavior and adaptive IAQ control. To address this gap, a dynamic indoor CO2 model is constructed using machine learning algorithms to forecast CO2concentrations across a range of forecasting horizons. Herein, we tuned and compared six state-of-the-algorithms—including Support Vector Machine, Ada Boost, Random Forest, Gradient Boosting, Logistic Regression, and Multilayer Perceptron. The algorithms’ performances are validated using CO 2 and historical meteorological data collected from a campus classroom with a variable occupancy rate. Simulation results showed that Multilayer Perceptron can strongly predict the volatile CO 2 behavior and also outperforms other algorithms in terms of accuracy. Furthermore, a control strategy capable of modeling and detecting dynamic patterns of CO 2 level is utilized to modulate the ventilation rate in real-time and also reduce the energy consumption. The proposed controller reduced the HVAC fan’s energy consumption by 51.4% and provide ventilation as needed per the ASHRAE standards.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Surrogate Hessian accelerated structural optimization for stochastic electronic structure theories

In this work, we present an efficient energy-based method for structural optimization with stochastic electronic structure theories, such as diffusion quantum Monte Carlo (DMC). This method is based on robust line-search energy minimization in reduced parameter space, exploiting approximate but accurate Hessian information from a surrogate theory, such as density functional theory. The surrogate theory is also used to characterize the potential energy surface, allowing for simple but reliable ways to maximize statistical efficiency while retaining controllable accuracy. We demonstrate the method by finding the minimum DMC energy structures of the selected flake-like aromatic molecules, such as benzene, coronene, and ovalene, represented by 2, 6, and 19 structural parameters, respectively. In each case, the energy minimum is found within two parallel line-search iterations. The method is near-optimal for a line-search technique and suitable for a broad range of applications. It is easily generalized to any electronic structure method where forces and stresses are still under active development and implementation, such as diffusion Monte Carlo, auxiliary-field Monte Carlo, and stochastic configuration interaction, as well as deterministic approaches such as the random-phase approximation. Accurate and efficient means of geometry optimization could shed light on a broad class of materials and molecules, showing high sensitivity of induced properties to structural variables.

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

Radiation Transport in Random Media With Large Fluctuations

Strongly heterogeneous media arise in several applications that include radiation shields, nuclear fuel, BWR moderators, clouds, planetary and stellar atmospheres, turbulent gases and plasmas. Neutral and charged particle transport computations in such media have relied heavily on formulating transport equations with spatially random coefficients (physical data) and developing solution methods to deal with the additional stochastic dimensions. Attempts at developing approximate closures that yield only low order statistical information (e.g., mean and variance of the flux) have proved to be highly restrictive under real physics conditions or rely on techniques that require fluctuation amplitudes to be small for robustness. Recently, stochastic spectral methods such as polynomial chaos and stochastic collocation have been developed for aleatoric uncertainty quantification and sensitivity analysis, and successfully applied in radiation transport work. Advances in UQ techniques have tended to focus on efficiently handling large numbers of uncertain variables but the rigorous stochastic basis of the approach also promotes its use in situations where stochasticity is due to spatial heterogeneity and the associated uncertainty is large. Here we apply these techniques to radiation transport in media with spatially randomly varying cross sections without restriction on fluctuation amplitudes. Specifically, we represent the cross sections as a lognormal spatial random process with specified mean, variance and covariance function and use a Karhunen-Loeve (KL) decomposition to generate cross section realizations that are strictly positive. Woodcock Monte Carlo (WMC) is then used to simulate transport using random sampling of cross sections and deterministic sampling based on a stochastic collocation technique. Numerical results for the mean and variance in the scalar flux and leakage currents are obtained for weak and strong random variations. We focus in this work on demonstrating the approach and defer application to specific problems to a future investigation. We have demonstrated that a combination of stochastic spectral representation and Woodcock Monte Carlo simulation allows efficient computation of radiation transport in strongly random media. The use of a lognormal distribution to represent the fluctuations in medium properties facilitates enforcement of strict positivity of cross section realizations, a condition that has challenged standard approaches that assume Gaussian fluctuations. The fact that the lognormal random process is a memoryless nonlinear transformation of a Gaussian process further enables efficient and accurate reconstruction of the random cross section realizations using a Karhunen-Loeve representation of the cross section. Finally, the robustness of this model enables fluctuations of arbitrarily large amplitude to be studied, using both random sampling Monte Carlo and more efficient stochastic collocation techniques. The present work is being extended to allow reconstruction of probability density functions of output variables using polynomial chaos representation, and to accommodate general cross section covariance functions. Also, techniques are being explored to more judiciously select ceiling cross sections to enhance the efficiency of Woodcock Monte Carlo for transport in random media. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Multiparameter Methods are the New Field in Statistics

During the last decade statistical problems of a new kind have appeared more and more often, in which the dimensionality of observations is large, and the sample size is small or comparable with the dimensionality of observations. These problems may be characterized as multiparameter problems. The theory of multiparameter analysis, proposed by A.N.Kolmogorov and supported by the studies on the spectral theory of random matrices of V. A.Marchenko, L. A.Pastur, V. L.Girko, and others, revealed some specific phenomena appearing in statistics with a large number of weakly dependent variables. Particularly, there are stable relations between the principal parts of parameter set functions and the set of observed variables, which may be used to improve the statistical analysis. The use of these phenomena allows one to develop a thorough, always stable, and approximately optimal (irrespective of samples) versions of mostly used statistical procedures. The multiparameter versions of the discriminant analysis are of special interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Diffusion in random networks

The ensemble averaging technique is applied to model mass transport by diffusion in random networks. The system consists of an ensemble of random networks, where each network is made of pockets connected by tortuous channels. Inside a channel, fluid transport is assumed to be governed by the one-dimensional diffusion equation. Mass balance leads to an integro-differential equation for the pocket mass density. The so-called dual-porosity model is found to be equivalent to the leading order approximation of the integration kernel when the diffusion time scale inside the channels is small compared to the macroscopic time scale. As a test problem, we consider the one-dimensional mass diffusion in a semi-infinite domain. Because of the required time to establish the linear concentration profile inside a channel, for early times the similarity variable is xt $-$1/4 rather than xt $-$1/2 as in the traditional theory. We found this early time similarity can be explained by random walk theory through the network.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Multitask graph neural networks for elastoplastic response prediction in dual-phase polycrystals

Microstructure-sensitive prediction of elastoplastic response remains a recurring bottleneck in multiscale damage and fatigue modeling, where large ensembles of statistically distinct polycrystals are required to quantify variability and extreme-value behavior. In this work, we develop a multitask graph neural network (GNN) surrogate that maps dual-phase ferrite–martensite polycrystal microstructures to Statistical Volume Element (SVE)-level elastoplastic Quantities of Interest (QoIs). Each SVE is represented as a grain-adjacency graph, with node features encoding phase, geometry, and crystallographic orientation, and edge features encoding relative misorientation. A message-passing graph convolution generates node embeddings, which are pooled into a graph representation and passed to a multitask regression head that jointly predicts 10 scalar QoIs and vector-valued stress–strain responses in orthogonal loading directions across multiple martensite volume fractions and SVE sizes. Results show high accuracy for scalar QoIs and strong agreement for full stress–strain trajectories, with population envelopes reproducing both median behavior and finite-SVE variability across compositions and partition scales. A unified model trained on pooled volume-fraction data preserves most within-regime accuracy relative to regime-specific models while also capturing the broader cross-regime variation reflected in the pooled test set. Distributional comparisons further demonstrate that the surrogate preserves heterogeneity under SVE partitioning, enabling statistically consistent block-wise random-field construction for mesoscale analyses. Overall, the proposed grain-graph surrogate provides a practical pathway to accelerate ensemble-based studies of SVE-level constitutive variability in dual-phase polycrystals.

Crystal plasticity↗
NASA NTRS2022

Estimating Species-Specific Leaf Area Index and Basal Area Using Optical and SAR Remote Sensing Data in Acadian Mixed Spruce-Fir Forests, USA

This study combined Sentinel-1 synthetic aperture radar (SAR), Sentinel-2 multispectral, and site variable datasets to model leaf area index (LAI) and basal area per ha (BAPH) of two economically important tree species in Northeast, USA; red spruce (Picea rubens Sarg.; RS), and balsam fir (Abies balsamea (L.) Mill.; BF). We used Random Forest (RF), and Multi-Layer Perceptron (MLP) algorithms for LAI and BAPH modeling. The results showed that RF outperformed MLP by reducing the normalized root mean square error (nRMSE) by 0.01 and 0.06 for LAI and BAPH, respectively. The final variables selected for modeling of both LAI and BAPH indicated the superiority of Sentinel-2 variables over the Sentinel-1 SAR with minor contributions of site variables (mainly elevation). The red-edge spectral vegetation indices played a significant role in both LAI and BAPH estimation. We attained the lowest nRMSEs of 0.12, and 0.16 for the final LAI model of RS, and BF, respectively using Sentinel-2 and site variables. The lowest nRMSE for both RS and BF BAPH models was 0.12. As RS and BF are the primary host species for a cyclically occurring and most destructive pest of the region, eastern spruce budworm (Choristoneura fumiferana; SBW), these estimations will be useful to evaluate SBW dynamics in the region.

Forest inventory↗
DOE OSTI2015

SDSS J14584479+3720215: A BENCHMARK JHK{sub S} BLAZAR LIGHT CURVE FROM THE 2MASS CALIBRATION SCANS

Active galactic nuclei (AGNs) are well-known to exhibit flux variability across a wide range of wavelength regimes, but the precise origin of the variability at different wavelengths remains unclear. To investigate the relatively unexplored near-IR (NIR) variability of the most luminous AGNs, we conduct a search for variability using well sampled JHK{sub s}-band light curves from the Two Micron All Sky Survey (2MASS) calibration fields. Our sample includes 27 known quasars with an average of 924 epochs of observation over three years, as well as one spectroscopically confirmed blazar (SDSS J14584479+3720215) with 1972 epochs of data. This is the best-sampled NIR photometric blazar light curve to date, and it exhibits correlated, stochastic variability that we characterize with continuous auto-regressive moving average (CARMA) models. None of the other 26 known quasars had detectable variability in the 2MASS bands above the photometric uncertainty. A blind search of the 2MASS calibration field light curves for active galactic nucleus (AGN) candidates based on fitting CARMA(1,0) models (damped-random walk) uncovered only seven candidates. All seven were young stellar objects within the ρ Ophiuchus star forming region, five with previous X-ray detections. A significant γ-ray detection (5σ) for the known blazar using 4.5 yr of Fermi photon data is also found. We suggest that strong NIR variability of blazars, such as seen for SDSS J14584479+3720215, can be used as an efficient method of identifying previously unidentified γ-ray blazars, with low contamination from other AGNs.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Modeling Stochastic Variability in Multiband Time-series Data

In preparation for the era of time-domain astronomy with upcoming large-scale surveys, we propose a state-space representation of a multivariate damped random walk process as a tool to analyze irregularly-spaced multifilter light curves with heteroscedastic measurement errors. We adopt a computationally efficient and scalable Kalman filtering approach to evaluate the likelihood function, leading to maximum O(k 3 n) complexity, where k is the number of available bands and n is the number of unique observation times across the k bands. This is a significant computational advantage over a commonly used univariate Gaussian process that can stack up all multiband light curves in one vector with maximum O(k 3 n 3 ) complexity. Using such efficient likelihood computation, we provide both maximum likelihood estimates and Bayesian posterior samples of the model parameters. Three numerical illustrations are presented: (i) analyzing simulated five-band light curves for a comparison with independent single-band fits; (ii) analyzing five-band light curves of a quasar obtained from the Sloan Digital Sky Survey Stripe 82 to estimate short-term variability and timescale; (iii) analyzing gravitationally lensed g- and r-band light curves of Q0957+561 to infer the time delay. Two R packages, Rdrw and timedelay, are publicly available to fit the proposed models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Performance of reverse osmosis membrane with large feed pressure fluctuations from a wave-driven desalination system

Wave-driven desalination systems are proposed water treatment systems that involve reverse osmosis of seawater powered directly by wave motion. Such a configuration would result in drastic feed pressure fluctuations. For a technology conventionally operated with a constant feed condition, the effect of these variable pressures on membrane integrity and performance is unknown. Here, experiments were conducted with spiral wound membranes coupled to a system capable of producing feed pressure fluctuations of more than 400 psi. Feed composition included 5, 20, and 35 g/L NaCl, and a synthetic seawater at normal and 1.5x concentration. The variable feed conditions included sine-like pressure waves swings of 200-500 and 500-900 psi with frequencies of 1.25, 7.5, and 12 waves/min, and a model-generated random waveform. Between each wave experiment we performed membrane integrity tests at 650 psi and 25 g/L NaCl feed, which showed a 7.4% drop in the membrane's water permeability coefficient, an 18.4% flux decline, and more than 99% salt rejection over 1770 h of cumulative experimental time. Analysis of permeate samples showed high salt rejection. In general, variable feed pressure had no significant deleterious effect on membrane integrity or performance.

16 TIDAL AND WAVE POWER↗
DOE OSTI2023

Seasonal drivers of dissolved oxygen across a tidal creek–marsh interface revealed by machine learning

Abstract Dissolved oxygen (DO) is a key biogeochemical control in coastal systems, and its concentration and drivers vary markedly through time and space. This makes it difficult to accurately represent coastal DO and associated biogeochemical processes in models, limiting our ability to predict how these systems will respond to global change. We obtained high‐frequency (5‐min) in situ measurements of DO collected at three locations across the interface of a tidal creek and coastal marsh in the Pacific Northwest, USA. Random Forest machine learning models quantified the importance of three categories of environmental drivers (Aquatic, Climatic, and Terrestrial) of DO variability across the creek–marsh interface. We selected two 4‐month datasets representing Summer and Winter seasonal periods to test two hypotheses on the dominant drivers of DO at the coastal interface. We found that the Terrestrial driver—characterized by long periods of anaerobic conditions and episodic pulses in DO after floods—was most important during the Winter, whereas the Aquatic driver—characterized by variability over tidal, diel, and lunar cycles—was most important during the Summer. We explored how future climate change scenarios could alter the drivers of DO variability using a cumulative sums driver–response framework. Our results suggest that under climate change, Aquatic and Climatic drivers may increase in importance during the Summer, potentially linked to changing metabolic regimes and sea level, with Terrestrial driver importance potentially increasing during the Winter. Our approach highlights useful methods for understanding the spatiotemporal complexity of oxygen across coastal interfaces and quantifying the relative importance of distinct environmental drivers.

54 ENVIRONMENTAL SCIENCES↗