Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “variance components”

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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Spatio-Temporal Modeling of Field Surveyed Backsheet Degradation

Assessing photovoltaic module backsheet durability is critical to increasing module lifetime. Lab based accelerating testing has recently failed to predict large scale failures of widely adopted polymeric materials. Field surveyed data is critical to assess the performance of component lifetime. Using a documented field survey protocol, 13 field surveys where conducted. Each measurement is encoded with it's spatial location in respect to the other modules. By combining field survey data on degradation predictors with real time satellite weather data, data-driven predictive models of backsheet degradation were trained. LOESS models were constructed to investigate the spatial dependence of measurements. It was found that micro-climatic effects like tree-lines, ground surface changes, and elevation changes effected the magnitude and variance of the measurements. A GAM model was created to predict the value of degradation based on measured predictors. The model includes variables on the climate of the system and the location of each measurement in the PV mounting structure. The model performed well with an adj.R 2 of 0.95 for yellowness index prediction. The model was cross-validated using k-folds.

backsheet↗
DOE OSTI2022

Co-optimized machine-learned manifold models for large eddy simulation of turbulent combustion

Many modeling approaches in large eddy simulation (LES) of turbulent combustion employ a projection of the thermochemical state onto a low-dimensional manifold within state space to reduce the number of transported variables and hence computational cost. Flamelet-generated manifolds (FGM) is an example of a well-established, physics-based approach, but increasingly, principal component analysis (PCA) is being used as a data-driven method for generating manifold models. For both approaches, the nonlinear relationship between the location on the predefined manifold and the outputs of interest, such as reaction rates, can be tabulated or encoded in a neural network. This work proposes a new approach for manifold modeling that extends these existing approaches. A modified neural network structure simultaneously encodes the definition of the manifold variables, the nonlinear mapping, and the subfilter closure for LES. This allows all three of these aspects of the model to be co-optimized, generating a model from any source of combustion thermochemical state data. The manifold parameterizing variables are constrained to be linear combinations of species, as in FGM and PCA-based models, to aid in interpretability and implementation. For LES, subfilter variances of the manifold variables are also included as inputs. Two types of a priori analysis are performed to evaluate the new approach. In the first, the model is trained on data from one-dimensional premixed flames. In this case, the approach recovers the behavior of flamelet-based manifold approaches, and in fact slightly improves performance by identifying an optimized progress variable. The approach is also applied to data from direct numerical simulations of spherical ignition kernels in isotropic turbulence. For any specified manifold dimensionality, the new approach provides substantially lower prediction errors than a PCA-based model developed from the same data set. Additionally, the LES formulation of the new approach can provide accurate predictions for filtered reaction rates across a variety of filter widths.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1986

Viscoelastic deformation near active plate boundaries

Model deformations near the active plate boundaries of Western North America using space-based geodetic measurements as constraints are discussed. The first six months of this project were spent gaining familarity with space-based measurements, accessing the Crustal Dynamics Data Information Computer, and building time independent deformation models. The initial goal was to see how well the simplest elastic models can reproduce very long base interferometry (VLBI) baseline data. From the Crustal Dynamics Data Information Service, a total of 18 VLBI baselines are available which have been surveyed on four or more occasions. These data were fed into weighted and unweighted inversions to obtain baseline closure rates. Four of the better quality lines are illustrated. The deformation model assumes that the observed baseline rates result from a combination of rigid plate tectonic motions plus a component resulting from elastic strain build up due to a failure of the plate boundary to slip at the full plate tectonic rate. The elastic deformation resulting from the locked plate boundary is meant to portray interseismic strain accumulation. During and shortly after a large interplate earthquake, these strains are largely released, and points near the fault which were previously retarded suddenly catch up to the positions predicted by rigid plate models. Researchers judge the quality of fit by the sum squares of weighted residuals, termed total variance. The observed baseline closures have a total variance of 99 (cm/y)squared. When the RM2 velocities are assumed to model the data, the total variance increases to 154 (cm/y)squared.

Ward, S. N.↗
NASA NTRS2006

Observations of Intermediate-Scale Diurnal Waves in the Equatorial Mesosphere and Lower Thermosphere

The purpose of this study is to seek observational evidence for diurnal, vertically propagating inertia-gravity waves (IGWs) in the mesosphere and lower thermosphere (MLT). Numerical modeling studies indicate that vertically propagating IGWs excited by tropical heating provide a significant source of momentum for the semiannual oscillation (SAO) in the equatorial zonal mean winds. The power spectrum of these waves has a strong diurnal component. We analyze global patterns of ascending and descending node differences in MLT satellite temperatures, which are assumed in this study to be proxies for waves of diurnal period. Juxtaposed upon the planetary-scale features are localized variations with longitudinal wavelengths ranging from 25 deg to 50 deg, and with vertical wavelengths between 13 and 25 km. These intermediate-scale variations are spatially coherent, and persist for several weeks. Their amplitudes generally increase with altitude, while their phase structures suggest both eastward and westward propagation. The variance of wave numbers 9-17 in the MLT is examined in relation to the underlying stratospheric zonal mean zonal gradient winds. Stronger variances generally coincide with periods where the underlying zonal mean winds are relatively light, or unidirectional. The weak inverse association between variance strength and zonal wind magnitude is suggestive of a wave filtering mechanism.

Lieberman, R. S.↗
NASA NTRS2012

Comparison of Gravity Wave Temperature Variances from Ray-Based Spectral Parameterization of Convective Gravity Wave Drag with AIRS Observations

The realism of ray-based spectral parameterization of convective gravity wave drag, which considers the updated moving speed of the convective source and multiple wave propagation directions, is tested against the Atmospheric Infrared Sounder (AIRS) onboard the Aqua satellite. Offline parameterization calculations are performed using the global reanalysis data for January and July 2005, and gravity wave temperature variances (GWTVs) are calculated at z = 2.5 hPa (unfiltered GWTV). AIRS-filtered GWTV, which is directly compared with AIRS, is calculated by applying the AIRS visibility function to the unfiltered GWTV. A comparison between the parameterization calculations and AIRS observations shows that the spatial distribution of the AIRS-filtered GWTV agrees well with that of the AIRS GWTV. However, the magnitude of the AIRS-filtered GWTV is smaller than that of the AIRS GWTV. When an additional cloud top gravity wave momentum flux spectrum with longer horizontal wavelength components that were obtained from the mesoscale simulations is included in the parameterization, both the magnitude and spatial distribution of the AIRS-filtered GWTVs from the parameterization are in good agreement with those of the AIRS GWTVs. The AIRS GWTV can be reproduced reasonably well by the parameterization not only with multiple wave propagation directions but also with two wave propagation directions of 45 degrees (northeast-southwest) and 135 degrees (northwest-southeast), which are optimally chosen for computational efficiency.

convective gravity wave drag↗
DOE OSTI2024

Multilevel Monte Carlo Estimators For Derivative-Free Optimization Under Uncertainty

Optimization is a key tool for scientific and engineering applications; however, in the presence of models affected by uncertainty, the optimization formulation needs to be extended to consider statistics of the quantity of interest. Optimization under uncertainty (OUU) deals with this endeavor and requires uncertainty quantification analyses at several design locations; i.e., its overall computational cost is proportional to the cost of performing a forward uncertainty analysis at each design location. An OUU workflow has two main components: an inner loop strategy for the computation of statistics of the quantity of interest, and an outer loop optimization strategy tasked with finding the optimal design, given a merit function based on the inner loop statistics. Here, in this work, we propose to alleviate the cost of the inner loop uncertainty analysis by leveraging the so-called multilevel Monte Carlo (MLMC) method, which is able to allocate resources over multiple models with varying accuracy and cost. The resource allocation problem in MLMC is formulated by minimizing the computational cost given a target variance for the estimator. We consider MLMC estimators for statistics usually employed in OUU workflows and solve the corresponding allocation problem. For the outer loop, we consider a derivative-free optimization strategy implemented in the SNOWPAC library; our novel strategy is implemented and released in the Dakota software toolkit. We discuss several numerical test cases to showcase the features and performance of our approach with respect to its Monte Carlo single fidelity counterpart.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Accurate determination of genotypic variance of cell wall characteristics of a Populus trichocarpa pedigree using high-throughput pyrolysis-molecular beam mass spectrometry

Abstract Background Pyrolysis-molecular beam mass spectrometry (py-MBMS) analysis of a pedigree of Populus trichocarpa was performed to study the phenotypic plasticity and heritability of lignin content and lignin monomer composition. Instrumental and microspatial environmental variability were observed in the spectral features and corrected to reveal underlying genetic variance of biomass composition. Results Lignin-derived ions (including m/z 124, 154, 168, 194, 210 and others) were highly impacted by microspatial environmental variation which demonstrates phenotypic plasticity of lignin composition in Populus trichocarpa biomass. Broad-sense heritability of lignin composition after correcting for microspatial and instrumental variation was determined to be H 2 = 0.56 based on py-MBMS ions known to derive from lignin. Heritability of lignin monomeric syringyl/guaiacyl ratio ( S / G ) was H 2 = 0.81. Broad-sense heritability was also high (up to H 2 = 0.79) for ions derived from other components of the biomass including phenolics (e.g., salicylates) and C5 sugars (e.g., xylose). Lignin and phenolic ion abundances were primarily driven by maternal effects, and paternal effects were either similar or stronger for the most heritable carbohydrate-derived ions. Conclusions We have shown that many biopolymer-derived ions from py-MBMS show substantial phenotypic plasticity in response to microenvironmental variation in plantations. Nevertheless, broad-sense heritability for biomass composition can be quite high after correcting for spatial environmental variation. This work outlines the importance in accounting for instrumental and microspatial environmental variation in biomass composition data for applications in heritability measurements and genomic selection for breeding poplar for renewable fuels and materials.

09 BIOMASS FUELS↗
DOE OSTI2021

Flux Ropes, Turbulence, and Collisionless Perpendicular Shock Waves: High Plasma Beta Case

With the onset of solar maximum and the expected increased prevalence of interplanetary shock waves, Parker Solar Probe is likely to observe numerous shocks in the next few years. An outstanding question that has received surprisingly little attention has been how turbulence interacts with collisionless shock waves. Turbulence in the supersonic solar wind is described frequently as a superposition of a majority 2D and a minority slab component. Here we formulate a collisional perpendicular shock-turbulence transmission problem in a way that enables investigation of the interaction and transmission of quasi-perpendicular fluctuations such as magnetic flux ropes/islands and vortices as well as entropy and acoustic modes in the large plasma beta regime. We focus on the transmission of an upstream spectrum of these modes, finding that the downstream spectral amplitude is typically increased significantly (a factor of 10 or more), and that the upstream spectral index of the inertial range, and indeed the general spectral shape, is unchanged for the downstream magnetic variance, kinetic energy, and density variance. A comparison of the theoretically predicted downstream magnetic variance, kinetic energy, and density variance spectra with those observed at 1, 5, and 84 au by Wind, Ulysses, and Voyager 2 shows excellent agreement. The overall theoretically predicted characteristics of the transmission of turbulence across shocks observed in the solar wind appear to be largely consistent with recent observational studies by Pitňa et al. and Borovsky.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1995

Statistically generated weighted curve fit of residual functions for modal analysis of structures

A statistically generated weighting function for a second-order polynomial curve fit of residual functions has been developed. The residual flexibility test method, from which a residual function is generated, is a procedure for modal testing large structures in an external constraint-free environment to measure the effects of higher order modes and interface stiffness. This test method is applicable to structures with distinct degree-of-freedom interfaces to other system components. A theoretical residual function in the displacement/force domain has the characteristics of a relatively flat line in the lower frequencies and a slight upward curvature in the higher frequency range. In the test residual function, the above-mentioned characteristics can be seen in the data, but due to the present limitations in the modal parameter evaluation (natural frequencies and mode shapes) of test data, the residual function has regions of ragged data. A second order polynomial curve fit is required to obtain the residual flexibility term. A weighting function of the data is generated by examining the variances between neighboring data points. From a weighted second-order polynomial curve fit, an accurate residual flexibility value can be obtained. The residual flexibility value and free-free modes from testing are used to improve a mathematical model of the structure. The residual flexibility modal test method is applied to a straight beam with a trunnion appendage and a space shuttle payload pallet simulator.

Bookout, P. S.↗
NASA NTRS1979

Polar motion from laser range measurements of GEOS-3

Using two-day arcs of GEOS-3 laser data, simultaneous solutions for pole position components, x sub p and y sub p, and orbit elements have been obtained for the period spanning 3 February to 6 March 1976 using three NASA Goddard Space Flight Center laser stations located near Washington, D.C. (STALAS) and on the islands of Bermuda and Grand Turk. The results are in general agreement with the BIH results. However, because of the locations of the laser sites, the x sub p solution is weaker than the y sub p solution. The x sub p and y sub p estimates were smoothed with a straight line by weighted least squares using the variance associated with the pole estimates as weights in order to reflect the effect of widely different data distributions. The smoothed y sub p differs by one meter with respect to the BIH smoothed values and the smoothed x sub p differs by about two meters. Spectral analysis of the results has identified frequencies associated with the orbital motion indicating the need for further improvements in the model of the physical system.

Schutz, B. E.↗
DOE OSTI2025

Geothermal well testing pressure prediction by using a hybrid transformer model system: FORGE well use case

Geothermal has huge potential to become an indispensable component in achieving the goal of sustainable energy economy, given its capability to provide consistent baseload power to the electric grid. Injection tests are crucial in geothermal energy system as they naturally help to evaluate reservoir properties, understand fluid flow and even enhance reservoir performance. In this research, we developed a hybrid model system that integrates machine learning (ML) regression, a physics-based mathematical model, and transformer deep learning. Trained and validated using FORGE injection test dataset, this system can forecast the pressure variations both upward and downward over time. The pressure prediction achieved prediction accuracy within 3-6% variance of true pressure values. The system can significantly save time and reduce costs by testing only a few cycles and then using model predictions for further analysis, instead of conducting additional real injection cycle tests. The developed model system also holds promise for designing injection test processes and maintaining well production in geothermal energy. Presented at the IMAGE ‘25 Conference led by Shell.

FORGE↗
DOE OSTI2025

Computationally efficient subglacial drainage modelling using Gaussian process emulators: GlaDS-GP v1.0

Subglacial drainage models represent water flow at the ice–bed interface through coupled distributed and channelized systems to determine water pressure, discharge, and drainage system geometry. While they are used to understand processes such as the relationship between surface melt and ice flow, the number of uncertain model parameters and the computational cost of running models makes it difficult to adequately explore the high-dimensional parameter space and evaluate uncertainty in model predictions. Here, we develop Gaussian process (GP) emulators that make fast predictions with associated uncertainty of subglacial drainage model outputs. Using a truncated principal component (PC) basis representation, we construct a GP emulator for diurnally averaged subglacial water pressure. We also explore emulation of scalar variables describing drainage efficiency and configuration. We train the emulators using ensembles of up to 512 simulations varying eight parameters of the Glacier Drainage System (GlaDS) model on a synthetic domain intended to represent an ice-sheet margin. The emulators make predictions ∼ 1000 times faster than GlaDS simulations, with errors <3 % for the water pressure field and ∼ 5 %–9 % for drainage efficiency and configuration. We apply the emulators to explore the eight-dimensional parameter space by computing variance-based parameter sensitivity indices, finding that three parameters (ice flow coefficient, bed bump aspect ratio, and the subglacial cavity system conductivity) explain 90 % of the variance in modelled water pressure in response to parameter changes. The GP emulator approach described here is well suited to integrating observational data with models to make calibrated, credible predictions of subglacial drainage.

58 GEOSCIENCES↗
NASA NTRS1989

Probabilistic modeling for simulation of aerodynamic uncertainties in propulsion systems

The numerical simulation of the probabilistic aerothermodynamic response of propulsion system components to randomness in their environment was explored. The reusable rocket engine turbopumps were selected as an example because of the severe cryogenic environment in which they operate. The thermal and combustion instabilities, coupled with the engine thrust requirements from start up to shut down, lead to randomness in the flow variables and uncertainties in the aerodynamic loading. The probabilistic modeling of the turbopumps aerodynamic response was accomplished using the panel method coupled with fast probability integration methods. The aerodynamic response in the form of probabilistic rotor blades and splitter loading were predicted and the results presented for specified flow coefficient and rotor preswirl variance. Possible future applications of the aerothermodynamic probabilistic modeling in engine transient simulation, condition monitoring and engine life prediction are briefly discussed.

Hamed, Awatef↗
NASA NTRS1990

Probabilistic modeling for simulation of aerodynamic uncertainties in propulsion systems

The numerical simulation of the probabilistic aerothermodynamic response of propulsion system components to randomness in their environment was explored. The reusable rocket engine turbopumps were selected as an example because of the severe cryogenic environment in which they operate. The thermal and combustion instabilities, coupled with the engine thrust requirements from start up to shut down, lead to randomness in the flow variables and uncertainties in the aerosynamic loading. The probabilistic modeling of the turbopumps' aerodynamic response was accomplished using the panel method coupled with fast probability integration methods. The aerodynamic response in the form of probabilistic rotor blades and splitter loading were predicted and the results presented for specified flow coefficient and rotor preswirl variance. Possible future applications of the aerothermodynamic probabilistic modeling in engine transient simulation, condition monitoring, and engine life prediction are briefly discussed.

Hamed, A.↗
DOE OSTI2025

Bayesian calibration of irradiated graphite property models under high temperatures

Graphite under high temperatures and irradiation is central to advanced reactors. We develop a Bayesian calibration framework for graphite property models that explicitly represents model-data mismatch via a Gaussian-process discrepancy. The approach propagates uncertainty from parameters, experimental noise, and model form, with a hierarchical variance structure to capture group and cross-group noise. Using two predictive models across five grades (IG-110, NBG-18, PCEA, NBG-17, 2114) and four properties-irradiation-induced dimension change, creep, Young’s modulus change ratio, and coefficient of thermal expansion change ratio-we obtain average predictive-error reductions of 54%, 65%, 17%, and 17% when discrepancy is included. We illustrate engineering impact with a multiphysics model of a very-high-temperature reactor prismatic reflector brick, analyzing stresses under high fluence and temperature. Accounting for model discrepancy markedly improves predictive accuracy and provides a robust basis for reliable graphite component design in advanced reactors.

36 - MATERIALS SCIENCE↗
DOE OSTI2022

Energy Exascale Earth System Model v2.0.2

Second patch release of v2.0.0 Changes since v2.0.1 [Important change] Add and update SSP370 and SSP585 cases, add tests, fix use-case files, Change ocean and sea-ice IC for ARRM60to10. [EAM] allow thetaxx as a CAM_TARGET, enable northamericax4v1pg2_WC14to60E2r3 for AMIP, fix ndrop initialization, allow up to 15 history files. [EAMxx] allow use of readnl [EAM-MMF] allow transient SST case for C++ MMF, remove specific task/thread count for ESMT test, modify the variance transport diagnostic, reorg tests, [HOMME] Allow optional sponge layer. [ELM] add 2 land-atm compsets, update FATES to API 17.0.0, allow FATES sp mode, allow ELM harvest to drive FATES harvest, fix memleaks, reduce test build times, modified parameters for miscanthus and switchgrass based on calibration [MPAS-Ocean] turn on ocean BGC in BGC cases, fix interface locations for 60L PHC grid, add mode spec to conservation check streams, allow oRRS18to6v3 grid to run with JRA, update ocean and sea ice ICs for ARRM60to10 (needed for v2), GPU port of thickness tendency, add CMPASO-JRA1p4, [MPAS framework] Add new reproducible global sum module for MPAS components, Add mostRecentAccessTime attribute to streams [MPAS-landice] Update MALI version and Greenland mesh [MPAS-seaice] add single-cell test case, fix BGC restart, Adds omp critical directives for ice warnings seen on cori. [CIME] Fix component namelist creation bugs when NINST>1, add OpenACC, OpenMP and CXX GPU tests, fix nonBFB tests, stop using config_compilers, reduce ELM test build times with shared executable, update OpenMPI on Chrysalis. fix baseline handling, provenance handling with update to cime6.0.33, also use component-specific config_pes files and add the ones from master [Machine updates] cori modules after maint, Chrysalis to OpenMPI-4.1.3 [run_e3sm] replace default case name and group. [Externals] update SCORPIO to 1.3.2

ECP↗
DOE OSTI2022

Development of Short-Term Forecasting Models Using Plant Asset Data and Feature Selection

Nuclear power plants collect and store large volumes of heterogeneous data from various components and systems. With recent advances in machine learning (ML) techniques, these data can be leveraged to develop diagnostic and short-term forecasting models to better predict future equipment condition. Maintenance operations can then be planned in advance whenever degraded performance is predicted, thus resulting in fewer unplanned outages and the optimization of maintenance activities. This enables lower maintenance costs and improves the overall economics of nuclear power. This paper focuses on developing a short-term forecasting process that leverages a feature selection process to distill large volumes of heterogeneous data and predict specific equipment parameters. A variety of feature selection methods, including Shapley Additive Explanations (SHAP) and variance inflation factor (VIF), were used to select the optimal features as inputs for three ML methods: long short-term memory (LSTM) networks, support vector regression (SVR), and random forest (RF). Each combination of model and input features was used to predict a pump bearing temperature both 1 and 24 hours in advance, based on actual plant system data. The optimal inputs for the LSTM and SVR were selected using the SHAP values, while the optimal input for the RF consisted solely of the response variable itself. Each model produced similar 1-hour-ahead predictions, with root mean square errors (RMSEs) of roughly 0.006. For the 24-hour-ahead predictions, differences could be seen between LSTM, SVR, and RF, as reflected by model performances of 0.036 +- 0.014, 0.0026 +- 0, and 0.063 +- 0.004 RMSE, respectively. As big data and continuous online monitoring become more widely available, the proposed feature selection process can be used for many applications beyond the prediction of process parameters within nuclear infrastructure.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Analytics-at-scale of Sensor Data for Digital Monitoring in Nuclear Plants (3 rd Annual Report)

Nuclear power plants collect and store large volumes of heterogeneous data from various components and systems. With recent advances in machine learning (ML) techniques, these data can be leveraged to develop diagnostic and short-term forecasting models to better predict future equipment condition. Maintenance operations can then be planned in advance whenever degraded performance is predicted, thus resulting in fewer unplanned outages and the optimization of maintenance activities. This enables lower maintenance costs and improves the overall economics of nuclear power. This report primarily focuses on developing a short-term forecasting process that leverages a feature selection process to distill large volumes of heterogeneous data and predict specific equipment parameters. A variety of feature selection methods, including Shapley Additive Explanations (SHAP) and variance inflation factor (VIF), were used to select the optimal features as inputs for three ML methods: long short-term memory (LSTM) networks, support vector regression (SVR), and random forest (RF). Each combination of model and input features was used to predict a pump bearing temperature both 1 and 24 hours in advance, based on actual plant system data. The optimal inputs for the LSTM and SVR were selected using the SHAP values, while the optimal input for the RF consisted solely of the response variable itself. Each model produced similar 1-hour-ahead predictions, with root mean square errors (RMSEs) of roughly 0.006. For the 24-hour-ahead predictions, differences could be seen between LSTM, SVR, and RF, as reflected by model performances of 0.036 ± 0.014, 0.0026 ± 0, and 0.063 ± 0.004 RMSE, respectively. As big data and continuous online monitoring become more widely available, the proposed feature selection process can be used for many applications beyond the prediction of process parameters within nuclear infrastructure. This report summarizes the Fiscal Year 2021 research progress encompassing the (1) data cleaning and feature selection necessary for ML applications; (2) development of short-term forecasting models to predict future plant process parameters for both single and multiple time steps ahead; and (3) validation of the feature selection methods and short-term forecasting models given new data from different systems.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗