Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 739 records · Page 41

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↗

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↗

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↗

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↗

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↗

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.↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

Identifying Transient Candidates in the Dark Energy Survey Using Convolutional Neural Networks

The ability to discover new transient candidates via image differencing without direct human intervention is an important task in observational astronomy. For these kind of image classification problems, machine learning techniques such as Convolutional Neural Networks (CNNs) have shown remarkable success. In this work, we present the results of an automated transient candidate identification on images with CNNs for an extant data set from the Dark Energy Survey Supernova program, whose main focus was on using Type Ia supernovae for cosmology. By performing an architecture search of CNNs, we identify networks that efficiently select non-artifacts (e.g., supernovae, variable stars, AGN, etc.) from artifacts (image defects, mis-subtractions, etc.), achieving the efficiency of previous work performed with random Forests, without the need to expend any effort in feature identification. The CNNs also help us identify a subset of mislabeled images. Performing a relabeling of the images in this subset, the resulting classification with CNNs is significantly better than previous results, lowering the false positive rate by 27% at a fixed missed detection rate of 0.05.

79 ASTRONOMY AND ASTROPHYSICS↗

Understanding Growth Dynamics and Yield Prediction of Sorghum Using High Temporal Resolution UAV Imagery Time Series and Machine Learning

Unmanned aerial vehicles (UAV) carrying multispectral cameras are increasingly being used for high-throughput phenotyping (HTP) of above-ground traits of crops to study genetic diversity, resource use efficiency and responses to abiotic or biotic stresses. There is significant unexplored potential for repeated data collection through a field season to reveal information on the rates of growth and provide predictions of the final yield. Generating such information early in the season would create opportunities for more efficient in-depth phenotyping and germplasm selection. This study tested the use of high-resolution time-series imagery (5 or 10 sampling dates) to understand the relationships between growth dynamics, temporal resolution and end-of-season above-ground biomass (AGB) in 869 diverse accessions of highly productive (mean AGB = 23.4 Mg/Ha), photoperiod sensitive sorghum. Canopy surface height (CSM), ground cover (GC), and five common spectral indices were considered as features of the crop phenotype. Spline curve fitting was used to integrate data from single flights into continuous time courses. Random Forest was used to predict end-of-season AGB from aerial imagery, and to identify the most informative variables driving predictions. Improved prediction of end-of-season AGB (RMSE reduction of 0.24 Mg/Ha) was achieved earlier in the growing season (10 to 20 days) by leveraging early- and mid-season measurement of the rate of change of geometric and spectral features. Early in the season, dynamic traits describing the rates of change of CSM and GC predicted end-of-season AGB best. Late in the season, CSM on a given date was the most influential predictor of end-of-season AGB. The power to predict end-of-season AGB was greatest at 50 days after planting, accounting for 63% of variance across this very diverse germplasm collection with modest error (RMSE 1.8 Mg/ha). End-of-season AGB could be predicted equally well when spline fitting was performed on data collected from five flights versus 10 flights over the growing season. This demonstrates a more valuable and efficient approach to using UAVs for HTP, while also proposing strategies to add further value.

54 ENVIRONMENTAL SCIENCES↗

The Feinheit method: A phase-independent formulation of the period-luminosity relation for Cephids

The tight correlation between color and luminosity during the cyclical variation of a Cepheid is calibrated and shown to be a direct means by which most of the phase-dependent temperature-induced amplitude of the variable can be transformed away. The resulting Feinheit function, F = V - alpha (B-V) gives rise to a random-phase period-luminosity relation with an rms scatter of less than 0.20 mag. A comparison shows that the blue Feinheit method for distance determinations method from single-phase observations is as accurate as near-infrared photometry but has the added advantage of being able to use panoramic detectors for the data acquisition. The Feinheit function is identified as the surface area variation of the Cepheid during its cycle.

Madore, B. F.↗

Turbulence structure in microburst phenomena

The efect on turbulence of a variable mean wind along the flight path of an aircraft is modeled and analyzed. It is found that the effect of a variable head or tail wind alters the magnitude of the length-scale of sensed microburst turbulence, rendering turbulence more random than usually encountered in the upper atmosphere. This, coupled with accompanying aerodynamic lift loss experienced during the headwind-to-tailwind swing, is what collectively creates the hazardous environment for a microburst-encountering aircraft attempting to land during a thunderstorm.

Trevino, George↗

Orbital magnetoconductance in the variable-range-hopping regime

The orbital magnetoconductance (MC) in the variable-range-hopping (VRH) regime is evaluated by use of a model which approximately takes into account the interference among random paths in the hopping process. Instead of logarithmic averaging the MC is obtained by the critical percolating resistor method. The small-field MC is quadratic in H; it is positive deep in the VRH regime and changes sign when the zero-field conductivity is high enough. This behavior (except for the sign change) and the relevant magnetic field scale are in agreement with recent experiments. The calculated MC is always positive for strong fields and is predicted to saturate at sufficiently large fields.

Sivan, U.↗