Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “statistical modeling”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Bayesian Spatial Models for Projecting Corn Yields

Climate change is predicted to impact corn yields. Previous studies analyzing these impacts differ in data and modeling approaches and, consequently, corn yield projections. We analyze the impacts of climate change on corn yields using two statistical models with different approaches for dealing with county-level effects. The first model, which is novel to modeling corn yields, uses a computationally efficient spatial basis function approach. We use a Bayesian framework to incorporate both parametric and climate model structural uncertainty. We find that the statistical models have similar predictive abilities, but the spatial basis function model is faster and hence potentially a useful tool for crop yield projections. We also explore how different gridded temperature datasets affect the statistical model fit and performance. Compared to the dataset with only weather station data, we find that the dataset composed of satellite and weather station data results in a model with a magnified relationship between temperature and corn yields. For all statistical models, we observe a relationship between temperature and corn yields that is broadly similar to previous studies. We use downscaled and bias-corrected CMIP5 climate model projections to obtain detrended corn yield projections for 2020–2049 and 2069–2098. In both periods, we project a decrease in the mean corn yield production, reinforcing the findings of other studies. However, the magnitude of the decrease and the associated uncertainties we obtain differ from previous studies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Attribution of heterogeneous stress distributions in low-grain polycrystals under conditions leading to damage

In high-purity polycrystalline metallic materials, voids tend to favor grain boundaries as nucleation sites due to the elevated stress states produced by granular interactions and the weakened grain boundary from the relative atomic disorder. To quantify the key factors of this elevated stress state, simple compression of a small multi-grain cylinder of body-centered cubic tantalum was simulated using a single crystal plasticity model that incorporates non-Schmid effects. Four increasingly complex synthetic microstructures were created to tractably incorporate grain boundary interactions, and a statistically significant number of combinations were performed by varying the initial crystallographic orientations of the microstructure. Most of these simulations produce the maximum von Mises stress on a grain boundary and less frequently at the multi-grain junctions. To build a statistical model for the maximum von Mises stress at the grain boundary, physically based features that could contribute to the elevated stress state were selected. Then, a learning algorithm based on information theory was used to identify which of these features contributed the most information to the data set. The identified features include a grain’s propensity to accommodate both elastic and plastic deformations and their directional components. The misalignment of the direction of each grain’s mechanical response was found to be strongly correlated to the magnitude of the stress near the grain boundary. For all of the synthetic microstructures, the statistical models produce a residual distribution that is nearly Gaussian with a variance of, at most, 10% of the prior distribution. The successful performance of the statistical model implies the correct identification of the physical features that cause severe stress localization in polycrystalline materials. The statistical models constructed here can be used to formulate a physically motivated void nucleation model which is sensitive to a microstructure’s propensity to produce elevated stress states. As a result, these statistical models also enable the design of material microstructures, in which the crystallographic orientation is chosen to resist void nucleation.

36 MATERIALS SCIENCE↗
NASA NTRS1984

Simulation and assimilation of satellite altimeter data at the oceanic mesoscale

An improved "objective analysis' technique is used along with an altimeter signal statistical model, an altimeter noise statistical model, an orbital model, and synoptic surface current maps in the POLYMODE-SDE area, to evaluate the performance of various observational strategies in catching the mesoscale variability at mid-latitudes. In particular, simulated repetitive nominal orbits of ERS-1, TOPEX, and SPOT/POSEIDON are examined. Results show the critical importance of existence of a subcycle, scanning in either direction. Moreover, long repeat cycles ( 20 days) and short cross-track distances ( 300 km) seem preferable, since they match mesoscale statistics. Another goal of the study is to prepare and discuss sea-surface height (SSH) assimilation in quasigeostrophic models. Restored SSH maps are shown to meet that purpose, if an efficient extrapolation method or deep in-situ data (floats) are used on the vertical to start and update the model.

Demay, P.↗
DOE OSTI2020

Damage and Failure in a Statistical Crack Model

In this paper I will take a close look at a statistical crack model (SCM) as is used in engineering computer codes to simulate fracture at high strain rates. My general goal is to understand the macroscopic behavior effected by the microphysical processes incorporated into an SCM. More specifically, I will assess the importance of including local interactions between cracks into the growth laws of an SCM. My strategy will be to construct a numerical laboratory that represents a single computational cell containing a realization of a statistical distribution of cracks. The cracks will evolve by the microphysical models of the SCM, leading to quantifiable damage and failure of the computational cell. I will use the numerical data generated by randomly generated ensembles of the fracture process to establish scaling laws that will modify and simplify the implementation of the SCM into large scale engineering codes.

42 ENGINEERING↗
NASA NTRS2012

Using Statistical Multivariable Models to Understand the Relationship Between Interplanetary Coronal Mass Ejecta and Magnetic Flux Ropes

In-situ measurements of interplanetary coronal mass ejections (ICMEs) display a wide range of properties. A distinct subset, "magnetic clouds" (MCs), are readily identifiable by a smooth rotation in an enhanced magnetic field, together with an unusually low solar wind proton temperature. In this study, we analyze Ulysses spacecraft measurements to systematically investigate five possible explanations for why some ICMEs are observed to be MCs and others are not: i) An observational selection effect; that is, all ICMEs do in fact contain MCs, but the trajectory of the spacecraft through the ICME determines whether the MC is actually encountered; ii) interactions of an erupting flux rope (PR) with itself or between neighboring FRs, which produce complex structures in which the coherent magnetic structure has been destroyed; iii) an evolutionary process, such as relaxation to a low plasma-beta state that leads to the formation of an MC; iv) the existence of two (or more) intrinsic initiation mechanisms, some of which produce MCs and some that do not; or v) MCs are just an easily identifiable limit in an otherwise corntinuous spectrum of structures. We apply quantitative statistical models to assess these ideas. In particular, we use the Akaike information criterion (AIC) to rank the candidate models and a Gaussian mixture model (GMM) to uncover any intrinsic clustering of the data. Using a logistic regression, we find that plasma-beta, CME width, and the ratio O(sup 7) / O(sup 6) are the most significant predictor variables for the presence of an MC. Moreover, the propensity for an event to be identified as an MC decreases with heliocentric distance. These results tend to refute ideas ii) and iii). GMM clustering analysis further identifies three distinct groups of ICMEs; two of which match (at the 86% level) with events independently identified as MCs, and a third that matches with non-MCs (68 % overlap), Thus, idea v) is not supported. Choosing between ideas i) and iv) is more challenging, since they may effectively be indistinguishable from one another by a single in-situ spacecraft. We offer some suggestions on how future studies may address this.

Riley, P.↗
DOE OSTI2018

Sensitivity Analysis and Emulation for Functional Data using Bayesian Adaptive Splines

When a computer code is used to simulate a complex system, one of the fundamental tasks is to assess the sensitivity of the simulator to the different input parameters. In the case of computationally expensive simulators, this is often accomplished via a surrogate statistical model, a statistical output emulator. An effective emulator is one that provides good approximations to the computer code output for wide ranges of input values. In addition, an emulator should be able to handle large dimensional simulation output for a relevant number of inputs; it should flexibly capture heterogeneities in the variability of the response surface; it should be fast to evaluate for arbitrary combinations of input parameters, and it should provide an accurate quantification of the emulation uncertainty. In this paper we discuss the Bayesian approach to multivariate adaptive regression splines (BMARS) as an emulator for a computer model that outputs curves. We introduce modifications to traditional BMARS approaches that allow for fitting large amounts of data and allow for more efficient MCMC sampling. We emphasize the ease with which sensitivity analysis can be performed in this situation. We present a sensitivity analysis of a computer model of the deformation of a protective plate used in pressure-driven experiments. Our example serves as an illustration of the ability of BMARS emulators to fulfill all the necessities of computability, flexibility and reliable calculation on relevant measures of sensitivity.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Predicting electromagnetic ion cyclotron wave amplitude from unstable ring current plasma conditions

Electromagnetic ion cyclotron (EMIC) waves in the Earth's inner magnetosphere are enhanced fluctuations driven unstable by ring current ion temperature anisotropy. EMIC waves can resonate with relativistic electrons and play an important role in precipitation of MeV radiation belt electrons. In this study, we investigate the excitation and saturation of EMIC instability in a homogeneous plasma using both linear theory and nonlinear hybrid simulations. We have explored a four-dimensional parameter space, carried out a large number of simulations, and derived a scaling formula that relates the saturation EMIC wave amplitude to initial plasma conditions. Finally, such scaling can be used in conjunction with ring current models like ring current-atmosphere interactions model with self-consistent magnetic field to provide global dynamic EMIC wave maps that will be more accurate inputs for radiation belt modeling than statistical models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Particle Size Distribution Dynamics Can Help Constrain the Phase State of Secondary Organic Aerosol

Particle phase state is a property of atmospheric aerosols that has important implications for the formation, evolution, and gas/particle partitioning of secondary organic aerosol (SOA). In this work, we use a size-resolved chemistry and microphysics model (SOM-TOMAS), updated to include an explicit treatment of particle phase state, to constrain the bulk diffusion coefficient (Db) of SOA produced from ??- pinene ozonolysis. By leveraging data from laboratory experiments performed in the absence of a seed and under dry conditions, we find that the Db for SOA can be constrained (1-7 ×10-15 cm2 s-1 in these experiments) by simultaneously reproducing the time-varying SOA mass concentrations and the evolution of the particle size distribution. Another version of our model that used the predicted SOA composition to calculate the glass transition temperature, viscosity, and, ultimately, Db (~10-15 cm2 s-1) of the SOA was able to reproduce the mass and size distribution measurements when we included oligomer formation (oligomers accounted for about a fifth of the SOA mass). Our work highlights the potential of a size resolved SOA model to constrain the particle phase state of SOA by utilizing historical measurements of the evolution of the particle size distribution.

He, Yicong↗
DOE OSTI2017

Ensemble Kalman Filter for Dynamic State Estimation of Power Grids Stochastically Driven by Time-correlated Mechanical Input Power

State and parameter estimation of power transmission networks is important for monitoring power grid operating conditions and analyzing transient stability. Wind power generation depends on fluctuating input power levels, which are correlated in time and contribute to uncertainty in turbine dynamical models. The ensemble Kalman filter (EnKF), a standard state estimation technique, uses a deterministic forecast and does not explicitly model time-correlated noise in parameters such as mechanical input power. However, this uncertainty affects the probability of fault-induced transient instability and increased prediction bias. Here a novel approach is to model input power noise with time-correlated stochastic fluctuations, and integrate them with the network dynamics during the forecast. While the EnKF has been used to calibrate constant parameters in turbine dynamical models, the calibration of a statistical model for a time-correlated parameter has not been investigated. In this study, twin experiments on a standard transmission network test case are used to validate our time-correlated noise model framework for state estimation of unsteady operating conditions and transient stability analysis, and a methodology is proposed for the inference of the mechanical input power time-correlation length parameter using time-series data from PMUs monitoring power dynamics at generator buses.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2019

Degradation Modeling and Reliability Assessment of Capacitors

The degradation of capacitors under accelerated stress conditions occur in a monotonic and non-linear fashion. Several efforts have been made to model the degradation behavior of capacitor considering either physics-of-failure models or statistical models and subsequently estimate its reliability and lifetime parameters. But most of these models fail to reflect the physical properties of the degradation path, which varies according to several intrinsic and extrinsic factors. These factors introduce random and temporal uncertainty among the population of capacitors. The gamma stochastic process can model both type of uncertainties among the population of capacitors. In this paper, we model the capacitor degradation by non-homogeneous gamma stochastic process in which both the model parameters (shape and scale) are dependent on stress variables. The model parameters are estimated using the maximum likelihood estimation approach.

42 ENGINEERING↗
DOE OSTI2015

Wind Power Curve Modeling in Simple and Complex Terrain

Our previous work on wind power curve modeling using statistical models focused on a location with a moderately complex terrain in the Altamont Pass region in northern California (CA). The work described here is the follow-up to that work, but at a location with a simple terrain in northern Oklahoma (OK). The goal of the present analysis was to determine the gain in predictive ability afforded by adding information beyond the hub-height wind speed, such as wind speeds at other heights, as well as other atmospheric variables, to the power prediction model at this new location and compare the results to those obtained at the CA site in the previous study. While we reach some of the same conclusions at both sites, many results reported for the CA site do not hold at the OK site. In particular, using the entire vertical profile of wind speeds improves the accuracy of wind power prediction relative to using the hub-height wind speed alone at both sites. However, in contrast to the CA site, the rotor equivalent wind speed (REWS) performs almost as well as the entire profile at the OK site. Another difference is that at the CA site, adding wind veer as a predictor significantly improved the power prediction accuracy. The same was true for that site when air density was added to the model separately instead of using the standard air density adjustment. At the OK site, these additional variables result in no significant benefit for the prediction accuracy.

17 WIND ENERGY↗
DOE OSTI2017

A Hybrid Multi-Scale Model of Crystal Plasticity for Handling Stress Concentrations

Microstructural effects become important at regions of stress concentrators such as notches, cracks and contact surfaces. A multiscale model is presented that efficiently captures microstructural details at such critical regions. The approach is based on a multiresolution mesh that includes an explicit microstructure representation at critical regions where stresses are localized. At regions farther away from the stress concentration, a reduced order model that statistically captures the effect of the microstructure is employed. The statistical model is based on a finite element representation of the orientation distribution function (ODF). As an illustrative example, we have applied the multiscaling method to compute the stress intensity factor K I around the crack tip in a wedge-opening load specimen. The approach is verified with an analytical solution within linear elasticity approximation and is then extended to allow modeling of microstructural effects on crack tip plasticity.

36 MATERIALS SCIENCE↗
NASA NTRS1981

Applications systems verification and transfer project. Volume 4: Operational applications of satellite snow cover observations. Colorado Field Test Center

The study was conducted on six watersheds ranging in size from 277 km to 3460 km in the Rio Grande and Arkansas River basins of southwestern Colorado. Six years of satellite data in the period 1973-78 were analyzed and snowcover maps prepared for all available image dates. Seven snowmapping techniques were explored; the photointerpretative method was selected as the most accurate. Three schemes to forecast snowmelt runoff employing satellite snowcover observations were investigated. They included a conceptual hydrologic model, a statistical model, and a graphical method. A reduction of 10% in the current average forecast error is estimated when snowcover data in snowmelt runoff forecasting is shown to be extremely promising. Inability to obtain repetitive coverage due to the 18 day cycle of LANDSAT, the occurrence of cloud cover and slow image delivery are obstacles to the immediate implementation of satellite derived snowcover in operational streamflow forecasting programs.

Shafer, B. A.↗
NASA NTRS1994

Empirical convection models for northward IMF

It is clear that polar cap convection during times of northward Interplanetary Magnetic Field (IMF) is more structured and of lower mean speed than at times of southward IMF. This, coupled with the fact that the polar cap is smaller, means that empirical models are more difficult to construct with certainty. It is also clear that sunward flow deep in the polar cap is often observed, but its connection with the rest of the flow pattern is controversial. At present, empirical models are of three types: 'statistical' models wherein data from different days but with similar IMF conditions are averaged together; 'pattern recognition' models, which are built up by examining individually hundreds of passes to derive a 'typical' pattern which embodies features frequently observed; and 'assimilative' models, which use data of different types and from as many locations as possible, but all taken at the same time, in order to derive a snapshot (or series of snapshots) of the entire pattern. Each type of model has its own difficulties. Statistical models, by their very nature, smooth out flow features (e.g. the convection reversal, and the locus of sunward flow deep in the polar cap) which are not found at precisely the same invariant latitudes and magnetic local times on different days. Pattern recognition models are better at reproducing small-scale features, but the large-scale pattern can be a matter of interpretation. Assimilative models (such as AMIE) hold out the best hope for creating instantaneous, global convection patterns; however, the analysis technique tends to be most irregular (and least reliable) in the regions which are not well covered by in situ data. It appears that, at least at times, a four cell model with sunward flow at the highest and lowest latitudes, and antisunward flow in between, is consistent with the observations. At other times, the observations may be consistent with a two-cell convection pattern, but which includes significant meanders within the polar cap.

Moses, Julie J.↗
NASA NTRS1998

Interactions between CO2 chemoreflexes and arterial baroreflexes

We studied interactions between CO2 chemoreflexes and arterial baroreflexes in 10 supine healthy young men and women. We measured vagal carotid baroreceptor-cardiac reflexes and steady-state fast Fourier transform R-R interval and photoplethysmographic arterial pressure power spectra at three arterial pressure levels (nitroprusside, saline, and phenylephrine infusions) and three end-tidal CO2 levels (3, 4, and 5%, fixed-frequency, large-tidal-volume breathing, CO2 plus O2). Our study supports three principal conclusions. First, although low levels of CO2 chemoreceptor stimulation reduce R-R intervals and R-R interval variability, statistical modeling suggests that this effect is indirect rather than direct and is mediated by reductions of arterial pressure. Second, reductions of R-R intervals during hypocapnia reflect simple shifting of vagally mediated carotid baroreflex responses on the R-R interval axis rather than changes of baroreflex gain, range, or operational point. Third, the influence of CO2 chemoreceptor stimulation on arterial pressure (and, derivatively, on R-R intervals and R-R interval variability) depends critically on baseline arterial pressure levels: chemoreceptor effects are smaller when pressure is low and larger when arterial pressure is high.

NASA Discipline Cardiopulmonary↗
NASA NTRS1991

Evaluation of an efficient statistical transport model for turbulent droplet dispersion in dilute combusting sprays

Evaluation of an efficient statistical transport model for turbulent droplet dispersion is made for a dilute spray of methanol droplets injected into a turbulent, axisymmnetric methane-fuelled diffusion flame burning in stagnant air. In the dispersion model, a computational parcel representing a group of physical particles (droplets) is considered to have a normal probability density function (pdf) in space. The mean is determined by Lagrangian tracking through a sequence of stochastically generated turbulent eddies and the variance is evaluated from a statistical formulation based on the linearized particle equations of motion. The basic validity of this model is established through comparison with available experimental data and with theoretical predictions using a conventional stochastic direct modeling approach. The conclusion of the evaluation is that the proposed dispersion model compares favorably with experimental data and provides a valid technique for simulating turbulent combusting sprays with significant computational savings over conventional methods.

Litchford, Ron J.↗
DOE OSTI2014

Automatic segmentation of head and neck CT images for radiotherapy treatment planning using multiple atlases, statistical appearance models, and geodesic active contours

Purpose: Accurate delineation of organs at risk (OARs) is a precondition for intensity modulated radiation therapy. However, manual delineation of OARs is time consuming and prone to high interobserver variability. Because of image artifacts and low image contrast between different structures, however, the number of available approaches for autosegmentation of structures in the head-neck area is still rather low. In this project, a new approach for automated segmentation of head-neck CT images that combine the robustness of multiatlas-based segmentation with the flexibility of geodesic active contours and the prior knowledge provided by statistical appearance models is presented. Methods: The presented approach is using an atlas-based segmentation approach in combination with label fusion in order to initialize a segmentation pipeline that is based on using statistical appearance models and geodesic active contours. An anatomically correct approximation of the segmentation result provided by atlas-based segmentation acts as a starting point for an iterative refinement of this approximation. The final segmentation result is based on using model to image registration and geodesic active contours, which are mutually influencing each other. Results: 18 CT images in combination with manually segmented labels of parotid glands and brainstem were used in a leave-one-out cross validation scheme in order to evaluate the presented approach. For this purpose, 50 different statistical appearance models have been created and used for segmentation. Dice coefficient (DC), mean absolute distance and max. Hausdorff distance between the autosegmentation results and expert segmentations were calculated. An average Dice coefficient of DC = 0.81 (right parotid gland), DC = 0.84 (left parotid gland), and DC = 0.86 (brainstem) could be achieved. Conclusions: The presented framework provides accurate segmentation results for three important structures in the head neck area. Compared to a segmentation approach based on using multiple atlases in combination with label fusion, the proposed hybrid approach provided more accurate results within a clinically acceptable amount of time.

62 RADIOLOGY AND NUCLEAR MEDICINE↗