Search NASA⌕ Search

SEARCH · Search NASA

Results for “index models”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Investigating the Global Biogeophysical Impact of Area and Mass Based Wood Harvest in a Vegetation Demography Model

Wood harvesting alters land surface properties and energy redistribution, but there is a lack of studies estimating these changes on a global scale. We coupled a vegetation demographic model, the Functionally Assembled Terrestrial Ecosystem Simulator, with the E3SM land model to perform offline model simulation to investigate the land biogeophysical responses, including canopy coverage, leaf area index, albedo, surface roughness length, and energy fluxes, to historical wood harvest on the global scale. In this study, we found 50% less harvested carbon (C) when choosing the area-based harvest rate as driving data that has not been spatially harmonized, compared to reharmonized mass-based harvesting. By considering the uncertainty from reconstruction of historical wood harvest time series and the choice of wood harvest approach in the model, continuous wood harvest (1850–2015) results in 5%–10% of canopy coverage loss, contributing 0.5%–1% increase of albedo over disturbed land, which is much stronger than a non-demographic land surface model. Changes in energy flux from the wood harvest are negligible (<1%), but the responses of land surface properties vary (up to 30%) due to differences in model structure between the single canopy, sun-shade leaf model and vegetation demographic model.

Shu, Shijie [Lawrence Berkeley National Laboratory↗
DOE OSTI2022

Fine resolution remote sensing spectra improves estimates of gross primary production of croplands

Gross primary production (GPP) is a fundamental measure of the terrestrial carbon cycle critical to our understanding of ecosystem function under the changing climate and land use. Remote sensing enables access to continuous spatial coverage, but remains challenged in heterogeneous croplands. Coarse resolution products, like MOD17A (500 m), may aggregate fragmented land cover types commonly found in heavily managed landscapes and misrepresent their respective contribution to carbon production. Consequently, this study demonstrates the capability of fine-resolution imagery (20-30 m) and available red-edge vegetation indices to characterize GPP across seven Midwest cropping systems. Four sites were established on a 22-year-old USDA Conservation Reserve Program (CRP); and the other three on land conventionally farmed with corn-soybean-wheat rotation (AGR). We compare in situ GPP estimates from eddy-covariance towers with ten satellite models: eight variants of the vegetation photosynthesis models (VPM), of which five include a red-edge vegetation index, as well as conventional products Landsat CONUS GPP (30 m) and MOD17A2H V6 (500 m). Daily and cumulative fine-resolution imagery integrated within VPM generally agreed with tower-based GPP in heterogeneous landscapes more than those from MODIS 500 m VPM or conventional GPP products from MOD17AH V6 or Landsat 8 CONUS. Here, replacing EVI2 with red-edge indices NDRE2, NDRE1, and MTCI in Sentinel 2 VPMs notably improved explanation of variance and estimation of cumulative GPP. While existing methods using MODIS- and Landsat-derived GPP are important baselines for regional and global studies, future research may benefit from the higher spatial, temporal, and radiometric resolution.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Exact wave solver for nonparaxial laser beam propagation

Simulations of inertial confinement fusion (ICF) experiments require high-fidelity models for laser beam propagation in a nonuniform plasma with varying index of refraction. We describe a new numerical wave solver that is applicable to centimeter-scale length plasmas encountered in indirect drive ICF applications. The one-way Helmholtz equation (OHE) generalizes the time-harmonic paraxial wave equation to large angles. Here, we present a methodology to numerically evaluate the exact solution to the OHE. This solution is computed by analytically advancing eigenfunctions of the one-way Helmholtz operator along a propagation direction and is applicable to any given index of a refraction profile. We compare our exact method with a commonly used approximate split-step technique for solving the OHE. As a test problem, we consider nonparaxial propagation of Gaussian and speckled beams in a plasma density channel with internal reflection. We find that the split-step approach incurs significant errors compared to the exact solution computed using the novel algorithm.

Belyaev, Mikhail A. (ORCID:0000000224908887)↗
DOE OSTI2020

Functional Photoresists for Energy Applications. Final Report

Monolithic ultralow-density porous bulk materials have recently attracted much interest due to many emerging applications in the areas of catalysis, energy storage and conversion, and thermal insulation. They are also important components of high energy density (HED) and inertial confinement fusion (ICF) targets. However, despite tremendous progress that has been made in the synthesis of porous materials, deterministic and independent control over microscopic architecture, density and composition remain key issues, and their integration in high precision devices requires cost and time-intensive mechanical machining that not only reduces reproducibility by generating debris but also limits the complexity of the 3D shapes that can be realized. In this project, we overcame these limitations by developing a universal templating capability that provides deterministic and independent control over density, composition, architecture, and macroscopic sample shape. This was achieved by developing the technology to 1) 3D print ultrahigh resolution, ultra-high precision polymeric micro-lattice templates, 2) coat these templates with the desired materials, and 3) removing the template (Fig. 1a). Atomic layer deposition (ALD) provides the atomic scale coating thickness accuracy required for precisely controlling density. While this templating approach had been demonstrated in prior work, limitations in suitable photoresists, 3D print technologies, print design, and template removal techniques did not allow the fabrication of millimeter-sized high-precision parts with sub-micron resolution. To enable this technology, we developed 1) two-photon polymerization (TPP) print designs that enable the fabrication of millimeter-sized, mechanically robust polymeric templates with sub-micron resolution and 2) a continuum level TPP printing simulation capability for additional print design guidance; 3) atomistic models to study photoresist polymerization kinetics and network topography, 4) refractive index matched polymeric and preceramic TPP photoresists, and 5) functional TPP photoresists including porous voxel structures and self-immolative polymer photoresist chemistries; and 6) damage free template removal techniques that enable the fabrication of defect-free high-precision low-density foam components. We also developed a templating approach for pure carbon microlattices with a unique tube-in-tube ligament morphology. As a test platform, we pursued the fabrication of foam liners that promise to further increase the neutron yield in indirect drive ICF experiments by improving implosion symmetry control and coupling between the laser and the deuterium-tritium fuel. This application requires fabrication and integration of a ultra-high precision, millimeter-sized, thin-walled (200-400 micrometer thick), low-density (10-30 mg/cc), high atomic number (high Z) cylindrical foam tube into the gold hohlraum of an indirect drive ICF target (Fig. 1b). While our hohlraum liner test case will mainly find application in HED and ICF experiments, the underlying science will also directly apply to previously developed nanoparticle and additive manufacturing technologies and will advance those techniques as well.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Short-Term Rainfall Prediction Based on Radar Echo Using an Improved Self-Attention PredRNN Deep Learning Model

Accurate short-term precipitation forecast is extremely important for urban flood warning and natural disaster prevention. In this paper, we present an innovative deep learning model named ISA-PredRNN (improved self-attention PredRNN) for precipitation nowcasting based on radar echoes on the basis of the advanced PredRNN-V2. We introduce the self-attention mechanism and the long-term memory state into the model and design a new set of gating mechanisms. To better capture different intensities of precipitation, the loss function with weights was designed. We further train the model using a combination of reverse scheduled sampling and scheduled sampling to learn the long-term dynamics from the radar echo sequences. Experimental results show that the new model (ISA-PredRNN) can effectively extract the spatiotemporal features of radar echo maps and obtain radar echo prediction results with a small gap from the ground truths. From the comparison with the other six models, the new ISA-PredRNN model has the most accurate prediction results with a critical success index (CSI) of 0.7001, 0.5812 and 0.3052 under the radar echo thresholds of 10 dBZ, 20 dBZ and 30 dBZ, respectively.

Wu, Dali (ORCID:0000000231177074)↗
DOE OSTI2020

A numerical study on the sooting tendencies of bio-derived fuels for spark-ignition engines

While more efficient vehicles have been developed over the decades, their performance is limited by the properties of existing fuels. In response, the Co-Optimization of Fuels and Engines Initiative (Co-Optima) under the U.S. Department of Energy (DOE) has developed a rigorous screening approach to identify the most promising biomass-derived blendstocks that are suitable for advanced spark-ignition (SI) engines. A detailed kinetic model has also been developed to predict the combustion properties of the selected blendstocks. This kinetic model has been designed mainly targeting ignition and flame propagation properties, and it has not previously been validated for soot formation. In this work, we numerically predicted the sooting tendency of 20 Co-Optima SI blendstocks using the Co-Optima kinetic model. The sooting tendencies are determined quantitatively using the Yield Sooting Index (YSI) methodology. As shown in Fig. 1, the predicted YSIs show good agreement with measurements, except for 2,5-dimethylfuran. We also quantify the sensitivity of the YSI predictions to aromatic growth reactions and fuel decomposition reactions of the test fuels. It is found that PAH growth reactions have little impact on YSIs for all 20 fuels under investigation. On the other hand, fuel decomposition reactions have a significant influence on the YSI of 2,5-dimethylfuran. Perturbing their reaction rates by a factor of two is shown to lead to 11 % YSI prediction uncertainty.

09 BIOMASS FUELS↗
DOE OSTI2020

Predicting the Cetane Number, Yield Sooting Index, Kinematic Viscosity, and Cloud Point for Catalytically Upgraded Pyrolysis Oil Using Artificial Neural Networks

Abstract The conversion of biomass using fast pyrolysis has the potential to be significantly less expensive at scale compared to alternative methods such as fermentation and gasification. Selective upgrading of the products of fast pyrolysis through chemical catalysis produces compounds with lower oxygen content and lower acidity; however, identifying the specific catalytic pathways for producing viable fuels and fuel additives often requires a trial-and-error approach. Specifically, key properties of the compounds must be experimentally tested to evaluate the viability of the resultant compounds. The present work proposes predictive models constructed with artificial neural networks (ANNs) for cetane number (CN), yield sooting index (YSI), kinematic viscosity (KV), and cloud point (CP), with blind test set median absolute errors of 5.14 cetane units, 3.36 yield sooting index units, 0.07 millimeters squared per second, and 4.89 degrees Celsius, respectively. Furthermore, the cetane number, yield sooting index, kinematic viscosity, and cloud point were predicted for over three hundred expected products from the catalytic upgrading of pyrolysis oil. It was discovered that 130 of these compounds have predicted cetane numbers greater than 40, with four of these compounds possessing predicted yield sooting index values significantly less than that of diesel fuel and predicted viscosities and cloud points comparable to that of diesel fuel.

09 BIOMASS FUELS↗
DOE OSTI2025

Clumping index estimation with 30°-tilted cameras in row crops: Evaluation of methods and segment size effects

The clumping index (CI) quantifies the spatial distribution of foliage elements and is essential for accurately estimating the plant area index (PAI), canopy radiative transfer, and photosynthesis. Traditionally, the finite-length averaging method (LX), the gap size distribution method (CC), and a combined approach of CC and LX (CLX) have been applied to instruments like TRAC and digital hemispherical photography to estimate CI. However, a comprehensive evaluation of these methods in row crops remains limited, especially regarding the influence of segment size on CI. Meanwhile, digital cameras offer a cost-effective and user-friendly solution for canopy measurements in row crops, yet their application in this context remains underexplored. In this study, we employed a new approach using a 30°-tilted digital camera to estimate CI in corn and soybean fields, applying the LX, CC, and CLX methods. We systematically assessed the performance of these three methods by combining field measurements in real-world fields with simulations using the LESS 3D radiative transfer model. Our results showed that CLX applied to the whole image and 45° segment offered accurate estimation of CI (bias within ±0.1, RMSE < 0.2) and PAI (bias within ±0.4, RMSE < 1) in real-world fields and LESS simulations. The accuracy of the LX method was highly sensitive to segment size, with the best performance observed at the 15° segment (PAI bias within ±0.4). In contrast, the CC method remained stable across different segment sizes, and its performance was generally comparable to that of LX, except at the 15° segment. Across view zenith angles, CI derived from CC generally showed a continuous increase, while those from LX and CLX followed a rising trend at small zenith angles but began to decline at 68°, likely due to an increasing proportion of no-gap segments. Seasonally, LX tended to show decreasing CI during early growth stages but increased as the canopy matured, whereas CC and CLX showed gradually increasing CI before plateauing at peak PAI. The 30°-tilted camera effectively captured CI variations across different angles and growth stages, making it a practical and robust instrument for row crop canopy structure analysis. Furthermore, applying these CI methods to digital cameras offers a low-cost and accessible CI estimation alternative, improving canopy structure monitoring accuracy in row crops.

30°-tilted camera↗
DOE OSTI2022

Predicting Switchgrass Biomass Yields Using a Spectral Vegetation Index Derived from Multispectral Satellite Imagery

Successful scaling of perennial bioenergy crop production requires a landscape design that optimizes the benefits of finite lands for people, communities, and environments. Utilizing marginal areas is the key to sustainable bioenergy crop production (Ssegane et al., 2015, 2016). Marginal areas are often small-sized lands and unevenly distributed across the agricultural landscape (Ssegane et al., 2016); thus, a systematic, semi-automated remote sensing method is needed as an effective means of estimating bioenergy crop yields across landscapes. Argonne National Laboratory (Argonne) is currently developing a tool, Scaling Up Perennial Bioenergy Economics and Ecosystem Services Tool (SUPERBEEST), to identify marginal agricultural lands and quantify environmental and economic effects of perennial bioenergy crop production systems. The tool aims to provide users a path to foster the sustainable and productive integration of bioenergy crops in the Midwestern agricultural landscape. Reliable, cost-effective, and timely estimation of bioenergy crop yields using remote sensing would help calculate and track the success of integrated bioenergy crops in the landscape for those communities. Argonne previously conducted feasibility studies for estimating biomass yields for bioenergy feedstock, corn and perennial grass using spectral vegetation indices (SVIs)1 derived from optical imagery (Hamada et al., 2015, 2021). In both studies, SVIs, more specifically those sensitive to plant chlorophyll or nitrogen contents, showed potential for estimating or predicting biomass yields with a correlation of determination (R2) ranging from 0.54 to 0.96, indicating a value for further investigation as a viable means of quantifying bioenergy feedstock production across large landscapes. Thus, the goal of this study is to evaluate the feasibility of use of SVIs as a means of estimating or predicting switchgrass biomass yields at harvest using publicly available multispectral satellite imagery. The feasibility analysis was performed using four study areas of mature switchgrass located in Virginia. Objectives are to (1) examine the SVIs and establish their relationships with switchgrass biomass yields at harvest, (2) develop a parsimonious image processing model for predicting at-harvest yields by applying the relationships with the most promising spectral index and (3) map switchgrass yields predicted by the image processing model across the study sites. The calibration to field data will rely on switchgrass biomass yields determined by the baling method, representing a potential challenge to the analysis but an important practical aspect for future applications. With this research design, the study aimed to gain insights into enabling remote sensing-based estimation of bioenergy crop yields in a reliable, cost- effective, and timely manner across large, heterogeneous landscapes.

09 BIOMASS FUELS↗
DOE OSTI2022

Data used in “Enguehard et al. 2022, Machine-Learning Functional Zonation Approach for Characterizing Terrestrial–Aquatic Interfaces: Application to Lake Erie”

The package contains the data layers used in “Enguehard et al. 2022, Machine-Learning Functional Zonation Approach for Characterizing Terrestrial–Aquatic Interfaces: Application to Lake Erie”. Spatial data layers include: topography, wetland vegetation cover, time series of Landsat’s enhanced vegetation index (EVI) between 1990 and 2020. The study aims to characterize coastal wetlands with particular focus on the co-variability between plant dynamics, topography, soil, and other environmental factors. We proposed a functional zonation approach based on machine learning clustering to identify the spatial regions, i.e., zones that capture these co-varied properties. This approach was applied to publicly available datasets along Lake Erie, in the Great Lakes Region

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Microreactor Core Transportation Cask Model Description for Criticality Safety Validation Basis Assessment (Rev. 2)

Criticality safety analyses are completed on transportation casks used for microreactor whole core shipment to provide examples of models and analyses to industry, regulators, and nuclear community at large to be used in verification and validation analyses of similar applications. The microreactors considered are based on a Gas-Cooled Microreactor (GCMR) and a Heat-Pipe Microreactor (HPMR), both utilize HALEU fuel in the form of TRISO particles and various other design options considered in industry microreactor designs. Variant design options of GCMR and HPMR were also investigated to provide a wider application range for each technology. Criticality safety analyses for the GCMR and HPMR packages were performed using the CSAS6 sequence of SCALE 6.3.2 with the ENDF/B-VII.1-based continuous energy neutron libraries. Different scenarios were investigated, including normal operation and water flooded conditions to represent nominal and hypothetical accident scenarios. Sensitivity and similarity analyses are also performed using the TSUNAMI sequence of SCALE 6.3.2, and the similarity analysis uses all the experiments from the ICSBEP Handbook with High Enriched Uranium (HEU), Intermediate and mixed Enriched Uranium (IEU) and Low Enriched Uranium (LEU) systems. Many ICSBEP experiments were found marginally similar, with similarity index (ck) values greater than 0.8, to the GCMR and HPMR cask models, especially in flooded conditions that may be more constraining due to reduced reactivity margins. Some experiments using TRISO fuel particles and graphite moderator, including the recent THETA and Deimos experiments from LANL, together with some IRPhEP experiments, were also used for the similarity analyses. These TRISO-fueled and graphite moderated experiments behaved similarly to the GCMR models, with maximum ck value greater than 0.9, but they are less similar to the HPMR models. This observation indicates that additional critical experiments might be needed to further validate the criticality safety models, especially for heat-pipe based microreactor transport packages and for non-flooded transportation configurations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

Coupling localized Noah-MP-Crop model with the WRF model improved dynamic crop growth simulation across Northeast China

Croplands play a critical role in regulating the energy and moisture exchanges between the land surface and atmosphere. However, the interactions between cropland and climate are usually poorly represented due to a lack of detailed representation in crop types and field management. Here, we coupled the Noah-MP-Crop model with the state-of-the-art Weather Research and Forecasting (WRF) model to explore and evaluate the crop growth dynamics in response to climate variations across Northeast China. The default parameters of the crop model were not exactly suitable for the agricultural ecosystems in Northeast China. The detailed cropland distribution, and crop phenology parameters including growing degree days (GDD) and planting (harvesting) date were first created using multi-source remote sensing products and reanalysis data, and was then successfully used to simulate the growth and yield for corn and soybean and associated energy exchanges. We also optimized and calibrated other crop parameters using the time-series of the Moderate Resolution Imaging Spectroradiometer (MODIS) land surface products. The modified crop model substantially improved the simulation of crop growth, plant physiology, and biomass accumulation for both corn and soybean. Coupling the localized dynamic crop model into the WRF led to considerable decreases in the simulated mean-absolute-errors (MAEs) and biases of the leaf area index, evapotranspiration, and gross primary production compared with the MODIS observed values. Compared with the statistical yield from each province, the modified crop model underestimated the corn yield from 11.1% to 48.6%, whereas overestimated the soybean yield from 16.5% to 162.6%.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Surface Energies of LaMnO 3 High-Index Surfaces Obtained from Density-Functional Theory

Here, the first six cubic LaMnO 3 high-index surfaces, (210), (211), (221), (310), (311), and (320), are examined for the first time. The unrelaxed and relaxed surface energies of the surfaces are computed as a function of surface model thickness using unrelaxed and relaxed surface models, respectively. The (210), (320), (211), and (221) surfaces are found to be more stable than the low-index (011) and (111) surfaces. Helping to explain this result, the surface terminations of the (210), (320), (211), or (221) surface are seen to exhibit a rotational relaxation of the MnO 6 or oxygen octahedra extending from the surface into the bulk. By contrast, the relaxed surface energies of the (310) and (311) surfaces are concluded to be difficult to determine due to a structural transformation or phase change in the surface models of the surfaces away from cubic LaMnO 3 . The finding of a phase change is important. It indicates that a phase change in the surface models of a surface for another material could occur if the surface is relatively open and a different phase of the material is stable at low temperature. When modeling such a surface, two steps are suggested to be taken.

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

Kronecker-structured covariance models for multiway data

Many applications produce multiway data of exceedingly high dimension. Modeling such multi-way data is important in multichannel signal and video processing where sensors produce multi-indexed data, e.g. over spatial, frequency, and temporal dimensions. We will address the challenges of covariance representation of multiway data and review some of the progress in statistical modeling of multiway covariance over the past two decades, focusing on tensor-valued covariance models and their inference. We will illustrate through a space weather application: predicting the evolution of solar active regions over time.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Quantifying model prediction sensitivity to model-form uncertainty

Computational and mathematical models are essential to understanding complex systems and phenomena. However, when developing such models, limited knowledge and/or resources necessitates the use of simplifying assumptions. It is therefore crucial to quantify the impact of such simplifying assumptions on the reliability and accuracy of resulting model predictions. This work develops a first-of-its-kind approach to quantify the impact of physics modeling assumptions on predictions. Here, we leverage the emerging field of model-form uncertainty (MFU) representations, which are parameterized modifications to modeling assumptions, in combination with grouped Sobol’ indices to quantitatively measure an assumption’s importance. Specifically, we compute the grouped Sobol’ index for the MFU representation’s parameters as a single importance measure of the assumption for which the MFU representation characterizes uncertainty. To ensure this approach is robust to the subjective choice of how to parameterize a MFU representation, we establish bounds for the difference between sensitivity results for two different MFU representations based on differences in model prediction statistics. The capabilities associated with this approach are demonstrated on three exemplar problems: an upscaled subsurface contaminant transport problem, ablation modeling for hypersonic flight, and nuclear waste repository modeling. We found that our grouped approach is able to assess the impact of modeling assumptions on predictions and offers computational advantages over classical Sobol’ index computation while providing more interpretable results.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

High-resolution leaf area index maps generated from unoccupied aerial system, Teller Mile 27, Seward Peninsula, Alaska

Leaf area index (LAI), a measure of the amount of one-side leaf area per ground unit, is an important indicator of plant carbon, energy, and water cycle. In the heterogeneous Arctic landscapes, it has been challenging to accurately measure LAI across species and space needed for Earth system model validation. Here, we use multispectral unoccupied aerial systems (UASs) to scale up and map leaf area index (LAI) , in a low-Arctic tundra landscape on the Seward Peninsula, Alaska. We linked previous published LAI measurements with high-resolution, UAS-collected multispectral data collected over the region of Next Generation Ecosystem Experiments in the Arctic (NGEE Arctic)’s Teller Mile Maker 27 site in 2022 to develop random forest (RF) machine learning models to predict and map LAI. 100 RF models were developed to account for uncertainties in ground LAI plot measurements and process scaling. This dataset includes a raster (*.tif) map of the mean LAI value of the 100 RF models, a raster (*.tif) map of the standard deviation of the RF-modeled LAI data, and a user guide (*.pdf).

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Spatially Accelerated Winding Numbers for Curved Geometry

The generalized winding number (GWN) is a scalar field that supports robust containment queries on curved geometry, including non-watertight, overlapping, and nested boundary representations. While queries can be easily parallelized over samples, direct evaluation on parametric curves and surfaces remains costly for large and complex models. Fast, state-of-the-art GWN approaches leverage a spatial index to approximate the GWN, typically coupled with a Taylor expansion which approximates the GWN contribution for far clusters of geometric primitives. However, such methods operate only on discrete inputs such as triangle meshes and point clouds, and would introduce containment errors near boundaries if applied to curved input. We extend support for fast GWN evaluation over arbitrary collections of NURBS curves in 2D and trimmed NURBS patches in 3D via a Bounding Volume Hierarchy that stores efficiently precomputed moment data in the hierarchy nodes. When querying the hierarchy, approximations for far clusters are used alongside direct evaluation for nearby NURBS primitives, achieving sub-linear complexity while preserving the geometric features in the vicinity of the query point. Central to our performance improvements is an adaptive subdivision strategy for NURBS primitives during a preprocessing phase, creating better spatial partitions while retaining the same accuracy for containment decisions as a direct evaluation. We demonstrate the performance and accuracy of our approach across a large collection of 2D and 3D datasets.

Computer science↗
DOE OSTI2021

Adaptive resource allocation for surrogate modeling of systems comprised of multiple disciplines with varying fidelity

We present an adaptive algorithm for constructing surrogate models for integrated systems composed of a set of coupled components. With this goal we introduce ‘coupling’ variables with a priori unknown distributions that allow approximations of each component to be built independently. Once built, the surrogates of the components are combined and used to predict system-level quantities of interest (QoI) at a fraction of the cost of interrogating the full system model. We use a greedy experimental design procedure, based upon a modification of Multi-Index Stochastic Collocation (MISC), to minimize the error of the combined surrogate. This is achieved by refining each component surrogate in accordance with its relative contribution to error in the approximation of the system-level QoI. Our adaptation of MISC is a multi-fidelity procedure that can leverage ensembles of models of varying cost and accuracy, for one or more components, to produce estimates of system-level QoI. Several numerical examples demonstrate the efficacy of the proposed approach on systems involving feed-forward and feedback coupling. For a fixed computational budget, the proposed algorithm is able to produce approximations that are orders of magnitude more accurate than approximations that treat the integrated system as a black-box.

97 MATHEMATICS AND COMPUTING↗