Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Model calibration”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Multiscale Approach to Fast ModSim for Laser Processing of Metals for Future Nuclear Deterrence Environments

Predicting performance of parts produced using laser-metal processing remains an out- standing challenge. While many computational models exist, they are generally too computationally expensive to simulate the build of an engineering-scale part. This work develops a reduced order thermal model of a laser-metal system using analytical Green's function solutions to the linear heat equation, representing a step towards achieving a full part performance prediction in an "overnight" time frame. The developed model is able to calculate a thermal history for an example problem 72 times faster than a traditional FEM method. The model parameters are calibrated using a non-linear solution and microstructures and residual stresses calculated and compared to a non-linear case. The calibrated model shows promising agreement with a non-linear solution.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Simulation of Creep Deformation and Failure in Graded AM Microstructures

This report describes modeling tools and techniques developed to simulate the long-term material performance of 316H stainless steel manufactured using Laser Powder Bed Fusion (LPBF). A physics-based Crystal Plasticity Finite Element model is used to simulate creep in microstructures and to study the roles of grain morphology, porosity, and texture. We describe our modeling methodology, including an orientation-mapping technique to capture the spatially varying crystallographic orientation that results from the build conditions. Our study of microstructural features shows that AM microstructures produced by LPBF tend to creep faster in the build direction, while texture and grain boundaries strengthen the transverse directions. However, when grain-boundary porosity and the consequent cavity growth are included in the model, the transverse directions begin to creep faster. In examining texture, the results indicate that spatially varying orientation arising from the build conditions increases anisotropy in the material, making it critical to account for orientation gradients in the material to accurately model its mechanical behavior. We also describe a material-model calibration campaign in which we calibrated the constitutive model specifically for LPBF 316H stainless steel at 725℃ for both solution-annealed and as-built conditions. Finally, these tools and techniques are used to model creep in microstructures representing different regions of an LPBF material with graded microstructure, owing to intentional variation in processing conditions. The creep simulation results show good agreement with experimental data across all three microstructures, with future work planned to study rupture in the material.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Perspectives on the integration between first-principles and data-driven modeling

Efficiently embedding and/or integrating mechanistic information with data-driven models is essential if it is desired to simultaneously take advantage of both engineering principles and data-science. Further the opportunity for hybridization occurs in many scenarios, such as the development of a faster model of an accurate high-fidelity computer model; the correction of a mechanistic model that does not fully-capture the physical phenomena of the system; or the integration of a data-driven component approximating an unknown correlation within a mechanistic model. At the same time, different techniques have been proposed and applied in different literatures to achieve this hybridization, such as hybrid modeling, physics-informed Machine Learning (ML) and model calibration. In this paper we review the methods, challenges, applications and algorithms of these three research areas and discuss them in the context of the different hybridization scenarios. Moreover, we provide a comprehensive comparison of the hybridization techniques with respect to their differences and similarities, as well as advantages and limitations and future perspectives. Finally, we apply and illustrate hybrid modeling, physics-informed ML and model calibration via a chemical reactor case study.

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

Thermal performance and energy consumption validation of an occupied local government office building outfitted with ceiling tile phase change materials

Buildings present an opportunity for energy conservation and the modulation of peak energy demand through controlled Heating, Ventilation, and Air Conditioning (HVAC) energy use. The administrative and office building stock in the United States holds potential to achieve energy and demand savings through retrofits such as insulation, weatherization, and thermal energy storage. Specifically, there is a need to validate passive phase change material (PCM) applications in full scale in aging administrative buildings in the US to evaluate the energy benefits. Aim of this study was to conduct a whole building level thermal and energy validation of an operational building and explore an alternative method for evaluating energy efficiency. To accomplish this, the study employed PCMs in the drop ceiling and carry out an energy audit and on-site measurement of HVAC systems' energy demand and consumption. A full-scale EnergyPlus energy model, modeled by the authors, served as a baseline for evaluation. The results show that calibrated model's envelope temperature measures fall within the accepted errors. HVAC energy simulation results also fall within the accepted errors for monthly and hourly pre- and post- PCM retrofit electricity and natural gas data. The novelty of this study is that it employees energy scales per Heating Degree Hour and Cooling Degree Hour, in contrast to the commonly used Heating Degree Days and Cooling Degree Days as reported in the literature to analyze energy savings. These findings underscore the pivotal role of a calibrated model in assessing the efficacy of a singular energy measure, like a PCM-retrofitted ceiling, in an occupied office building.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2023

Inference of strength and phase transition kinetics in dynamically-compressed tin

Dynamic compression experiments in condensed matter are of interest in part because they provide opportunities to examine material response under extreme conditions; however, the inference of material behavior from dynamic experiments is challenging in the presence of phase transitions exhibiting kinetic processes. Demonstrating an approach to quantitative interpretation of such dynamic experiments, here we present a Bayesian model calibration of strength and phase transformation parameters to data drawn from pulsed power and gas gun shot experiments. The posterior predictions of the Bayesian model capture the experimental measurements and account for the various uncertainties in the experimental configurations. This holistic approach to model calibration utilizing multiple types of experimental data identifies important cross correlations among kinetics, strength, and the phase boundary. Improved insight into potential sources of current model form error is provided by comparing the differences between calibrations against different subsets of the experimental data.

42 ENGINEERING↗
DOE OSTI2022

Cross-correlation of Dark Energy Survey Year 3 lensing data with ACT and Planck thermal Sunyaev-Zel’dovich effect observations. I. Measurements, systematics tests, and feedback model constraints

Here we present a tomographic measurement of the cross-correlation between thermal Sunyaev-Zel’dovich (TSZ) maps from Planck and the Atacama Cosmology Telescope and weak galaxy lensing shears measured during the first three years of observations of the Dark Energy Survey. This correlation is sensitive to the thermal energy in baryons over a wide redshift range and is therefore a powerful probe of astrophysical feedback. We detect the correlation at a statistical significance of 21σ, the highest significance to date. We examine the TSZ maps for potential contaminants, including cosmic infrared background and radio sources, finding that cosmic infrared background has a substantial impact on our measurements and must be taken into account in our analysis. We use the cross-correlation measurements to test different feedback models. In particular, we model the TSZ using several different pressure profile models calibrated against hydrodynamical simulations. Our analysis marginalizes over redshift uncertainties, shear calibration biases, and intrinsic alignment effects. We also marginalize over Ω m and σ 8 using Planck or DES priors. We find that the data prefer the model with a low amplitude of the pressure profile at small scales, compatible with a scenario with strong active galactic nuclei feedback and ejection of gas from the inner part of the halos. When using a more flexible model for the shear profile, constraints are weaker, and the data cannot discriminate between different baryonic prescriptions.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Quantification of Lanthanides in a Molten Salt Reactor Surrogate Off-Gas Stream Using Laser-Induced Breakdown Spectroscopy

To enable the deployment of molten salt reactor technology, the development of off-gas treatment systems and advanced monitoring tools capable of operating with high temperatures and radiation fields while delivering near real-time information is necessary. Here, we aim to fulfill this requirement and proposes laser-induced breakdown spectroscopy (LIBS) for monitoring molten salt aerosol streams. A sheath gas measuring method was developed to protect optical elements from aerosol particles and to ensure a relatively constant aerosol stream for measurement. An aqueous system was studied to demonstrate the utility of LIBS for monitoring possible fission products in an aerosol stream: Gd, Nd, and Sm up to 2000 parts per million (ppm). A calibration model was built using partial least squares (PLS) regression with five, six, and nine latent variables for Gd, Nd, and Sm, respectively. This calibration model successfully estimated the concentrations of three test samples, which were validated with inductively charged plasma optical emission spectroscopy measurements at a 99.9% confidence interval. To enhance these models, a genetic algorithm was used to filter the spectra before entering the PLS model, thereby limiting the spectral features being regressed to those with greater correlations to concentration. This allowed for the number of latent variables used in the PLS models to be reduced to four, three, and three for Gd, Nd, and Sm, respectively. Finally, the genetic algorithm-filtered PLS models were used to predict the concentrations of the aerosol stream on a real-time dataset and resulted in a 73%, 18%, and 25% improvement in root mean squared error of prediction compared to the original PLS models developed.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

SIDDA: SInkhorn Dynamic Domain Adaptation for image classification with equivariant neural networks

Modern neural networks (NNs) often do not generalize well in the presence of a ‘covariate shift’; that is, in situations where the training and test data distributions differ, but the conditional distribution of classification labels given the data remains unchanged. In such cases, NN generalization can be reduced to a problem of learning more robust, domain-invariant features. Domain adaptation (DA) methods include a broad range of techniques aimed at achieving this; however, these methods have struggled with the need for extensive hyperparameter tuning, which then incurs significant computational costs. In this work, we introduce SInkhorn Dynamic Domain Adaptation (SIDDA), an out-of-the-box DA training algorithm built upon the Sinkhorn divergence, that can achieve effective domain alignment with minimal hyperparameter tuning and computational overhead. We demonstrate the efficacy of our method on multiple simulated and real datasets of varying complexity, including simple shapes, handwritten digits, real astronomical observations, and remote sensing data. These datasets exhibit covariate shifts due to noise, blurring, differences between telescopes, and variations in imaging wavelengths. SIDDA is compatible with a variety of NN architectures, and it works particularly well in improving classification accuracy and model calibration when paired with symmetry-aware equivariant NNs (ENNs). We find that SIDDA consistently enhances the generalization capabilities of NNs, achieving up to a ${\approx}40\%$ improvement in classification accuracy on unlabeled target data, while also providing a more modest performance gain of $\lesssim 1\%$ on labeled source data. We also study the efficacy of DA on ENNs with respect to the varying group orders of the dihedral group DN, and find that the model performance improves as the degree of equivariance increases. Finally, if SIDDA achieves proper domain alignment, it also enhances model calibration on both source and target data, with the most significant gains in the unlabeled target domain—achieving over an order of magnitude improvement in the expected calibration error and Brier score. SIDDA’s versatility across various NN models and datasets, combined with its automated approach to domain alignment, has the potential to significantly advance multi-dataset studies by enabling the development of highly generalizable models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Grain size estimation in fluvial gravel bars using uncrewed aerial vehicles: A comparison between methods based on imagery and topography

Abstract Grain size assessments are necessary for understanding the various geomorphological, hydrological and ecological processes that occur within rivers. Recent research has shown that the application of Structure‐from‐Motion (SfM) photogrammetry to imagery from uncrewed aerial vehicles (UAVs) shows promise for rapidly characterising grain sizes along rivers in comparison to traditional field‐based methods. Here, we evaluated the applicability of different methods for estimating grain sizes in gravel bars along a study reach in the Olentangy River in Columbus, Ohio. We collected imagery of these gravel bars with a UAV and processed those images with SfM photogrammetry software to produce three‐dimensional point clouds and orthomosaics. Our evaluation compared statistical models calibrated on topographic roughness, which was computed from the point clouds, and to those based on image texture, which was computed from the orthomosaics. Our results showed that statistical models calibrated on image texture were more accurate than those based on topographic roughness. This might be because of site‐specific patterns of grain size, shape and imbrication. Such patterns would have complicated the detection of topographic signatures associated with individual grains. Our work illustrates that UAV‐SfM approaches show potential to be used as an accessible method for characterising surface grain sizes along rivers at higher spatial and temporal resolutions than those provided by traditional methods.

Wong, Tyler↗
DOE OSTI2021

Comparison of Spectroscopic Techniques for Determining the Peroxide Value of 19 Classes of Naturally Aged, Plant-Based Edible Oils

The peroxide value of edible oils is a measure of the degree of oxidation, which directly relates to the freshness of the oil sample. Several studies previously reported in the literature have paired various spectroscopic techniques with multivariate analyses to rapidly determine peroxide values using field portable and process instrumentation; those efforts presented “best-case scenarios” with oils from narrowly defined training and test sets. The purpose of this paper is to evaluate the use of near- and mid-infrared absorption and Raman scattering spectroscopies on oil samples from different oil classes, including seasonal and vendor variations, to determine which measurement technique or combination thereof is best for predicting peroxide values. Following peroxide value assays of each oil class using an established titration-based method, global and global-subset calibration models were constructed from spectroscopic data collected on the 19 oil classes used in this study. Spectra from each optical technique were used to create partial least squares regression calibration models to predict the peroxide value of unknown oil samples. A global peroxide value model based on near-infrared (8 mm optical path length) oil measurements produced the lowest RMSEP (4.9), followed by 24 mm optical path length near infrared (5.1), Raman (6.9) and 50 µm optical path length mid-infrared (7.3). However, it was determined that the Raman RMSEP resulted from chance correlations. Global peroxide value models based on low-level fusion of the NIR (8 and 24 mm optical path length) data and all infrared data produced the same RMSEP of 5.1. Global subset models, based on any of the spectroscopies and olive oil training sets from any class (pure, extra light, extra virgin), all failed to extrapolate to the non-olive oils. However, the near-infrared global subset model built on extra virgin olive oil could extrapolate to test samples from other olive oil classes. This work demonstrates the difficulty of developing a truly global method for determining peroxide value of oils.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Utilizing 3D Mechanical Earth Models for Calibration and Validation in a Large-Scale Carbon Capture and Storage Project in North Dakota

Conference paper presented at 17th International Conference on Greenhouse Gas Control Technologies (GHGT-17), Calgary, Alberta, Canada, October 20–24, 2024. Three-dimensional (3D) mechanical earth models (MEMs) are pivotal in assessing geological sites for carbon dioxide (CO 2 ) storage and mitigating potential risks associated with storage and injection. The Energy & Environmental Research Center is exploring innovative carbon storage monitoring techniques at a CO 2 storage site near Beulah, North Dakota. These methods aim to provide lower-impact, faster feedback for commercial carbon capture and storage (CCS) projects. The research includes monitoring CO 2 injection using various approaches: 1) an automated, integrated modular monitoring station; 2) time-lapse electromagnetic (EM) field surveys; 3) drone-based surveillance; 4) time-lapse seismic techniques; and 5) advanced wellbore monitoring.

02 PETROLEUM↗
DOE OSTI2025

SIDDA: SInkhorn Dynamic Domain Adaptation

Modern neural networks (NNs) often do not generalize well in the presence of a "covariate shift"; that is, in situations where the training and test data distributions differ, but the conditional distribution of classification labels remains unchanged. In such cases, NN generalization can be reduced to a problem of learning more domain-invariant features. Domain adaptation (DA) methods include a range of techniques aimed at achieving this; however, these methods have struggled with the need for extensive hyperparameter tuning, which then incurs significant computational costs. In this work, we introduce SIDDA, an out-of-the-box DA training algorithm built upon the Sinkhorn divergence, that can achieve effective domain alignment with minimal hyperparameter tuning and computational overhead. We demonstrate the efficacy of our method on multiple simulated and real datasets of varying complexity, including simple shapes, handwritten digits, and real astronomical observations. SIDDA is compatible with a variety of NN architectures, and it works particularly well in improving classification accuracy and model calibration when paired with equivariant neural networks (ENNs). We find that SIDDA enhances the generalization capabilities of NNs, achieving up to a ≈40% improvement in classification accuracy on unlabeled target data. We also study the efficacy of DA on ENNs with respect to the varying group orders of the dihedral group DN, and find that the model performance improves as the degree of equivariance increases. Finally, we find that SIDDA enhances model calibration on both source and target data--achieving over an order of magnitude improvement in the ECE and Brier score. SIDDA's versatility, combined with its automated approach to domain alignment, has the potential to advance multi-dataset studies by enabling the development of highly generalizable models.

Pandya, Sneh [Northeastern U.]↗
DOE OSTI2022

Thermo-Hydro-Chemical Simulation of Mid-Ocean Ridge Hydrothermal Systems: Static 2D Models and Effects of Paleo-Seawater Chemistry

Decades of research have resulted in characterization of the ocean floor manifestations of mid-ocean ridge (MOR) hydrothermal systems, yet numerical models accounting for the connections between heat transfer, hydrology and geochemistry have been slow to develop. The Thermo-hydro-chemical code ToughReact can be used to describe the coupled effects of fluid flow, heat transfer, and fluid-rock chemical interactions that occur in MOR systems. We describe the results of 2-dimensional simulations of steady state flow in fractured diabase with mineral-fluid chemical reactions. Basal heating and specified permeability yield maximum temperature of 400°C. Total fluid flux and high fracture flow velocities are in accord with observations. Fluid chemistry, mineralogical changes and 87 Sr/ 86 Sr ratios can be compared to observations to assess and calibrate models. Simulated high temperature fracture fluids have Mg and SO 4 near zero, elevated Ca and 87 Sr/ 86 Sr of about 0.7040. Total alteration is 10%–50% for simple models of spreading. Anhydrite forms mainly near the base of the upwelling zone and results in substantial local fracture porosity reduction. A calibrated model is used to predict how Sr isotopes and other features of altered oceanic crust would be different in the Cretaceous (95 Ma) early Proterozoic (1,800 Ma) and Archean (3,800 Ma), when seawater may have had high Ca and Sr concentrations, lower pH, higher temperature, and lower Na, Mg, and SO 4 . The simulations are offered as a start on what ultimately may require a longer-term community effort to better understand the role of MOR thermo-hydro-chemical systems in Earth evolution.

58 GEOSCIENCES↗
DOE OSTI2020

Application of Neutron Multiplicity Counting Experiments to Optimal Cross-Section Adjustments

This paper presents the first application of model calibration to neutron multiplicity counting (NMC) experiments for cross-section optimization that is informed by adjoint-based sensitivity analysis (SA) and first-order uncertainty quantification (UQ). We summarize previous work on SA applied to NMC and describe notable modifications and additions. We give the procedure for first-order UQ and Bayesian-inference-based parameter estimation (PE). We then discuss model calibration applied to NMC of a 4.5-kg sphere of weapons-grade, alpha-phase plutonium metal (the BeRP ball) with the nPod neutron multiplicity counter. For the BeRP ball in bare and polyethylene-reflected configurations, we discuss the sensitivity of the first- and second-moment detector responses (i.e., first and second moments of the NMC distribution, respectively) to the cross sections. We describe the sources of uncertainty in the measured and simulated responses. Specifically, uncertainty in the measured responses is due to both random and systematic sources. Uncertainty in the simulated responses is due to the cross-section covariances. We describe in detail the adjustment to the cross sections and cross-section covariances due to the optimization. Due to the contribution of systematic uncertainties to the measured response uncertainties, the adjustment to the cross sections is similar in trend but larger in magnitude compared to that recommended by previous work. We compare the measured responses to responses simulated with nominal and optimized cross sections, demonstrating that the best-estimate cross sections produce simulations of NMC experiments that are more accurate with reduced uncertainty.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Null Space Monte Carlo Evaluation of the Plateau to River Model

The Plateau to River Groundwater Model (P2R Model) is a groundwater flow and contaminant fate and transport (F&T) simulation model used to support remedial activities conducted by CH2M HILL Plateau Remediation Company at the Hanford Site in Washington State. Figure 1-1 illustrates the P2R Model extents, discretization, and boundary conditions. The P2R Model provides a computational framework to simulate the F&T of contaminants in groundwater associated with the 200-PO-1, 200-UP-1, 200-BP-5, and 200-ZP-1 Groundwater Operable Units (OUs) in the Hanford Site Central Plateau. In addition, the model includes adjacent areas and facilities (e.g., the State Approved Land Disposal Site). Intended and anticipated uses of the model include calculating water levels, hydraulic gradients, and groundwater flows throughout the model domain (encompassing the 200 West and 200 East Areas) for use in subsequent F&T calculations for contaminants of concern and developing scale-appropriate, telescopic-mesh refinement models for detailed evaluation of areas within the model domain where required. The overall objective of the modeling effort is to provide a basis for making informed remedial action decisions based on descriptions of current and expected future contaminant concentrations in groundwater at decision points within the OU boundaries. The objective for the model development phase is to create a common modeling platform that can be used for investigations of the Central Plateau groundwater OUs and areas downgradient toward the Columbia River. The P2R Model calibration to historical data observed at the Hanford Site is documented in CP-57037, Model Package Report for the Plateau to River Model Version 8.3. The purpose of this environmental calculation is to describe a null space Monte Carlo (NSMC) evaluation was conducted with the historic calibration of the P2R model. Use of numerical groundwater models is always accompanied with uncertainty in the results produced by a model because models are approximations of reality. Thus, by definition, lack the detail to fully represent observed behavior. Use of numerical techniques, such as a NSMC analysis, can help in identifying and quantifying the potential uncertainties associated with a numerical model such as the P2R Model. Use of the NSMC approach results in 100 groundwater flow models that are variants of the calibrated P2R Model. These variant models can be used to evaluate uncertainty in model predictions made by the calibrated P2R Model for other analyses. A secondary purpose of the environmental calculation is to establish these variant models for use with other applications.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2022

Oil Recovery Prediction for Polymer Flood Field Test of Heavy Oil on Alaska North Slope Via Machine Assisted Reservoir Simulation

The first ever polymer flood field pilot to enhance the recovery of heavy oils on the Alaska North Slope is ongoing. This study constructs and calibrates a reservoir simulation model to predict the oil recovery performance of the pilot through machine-assisted reservoir simulation techniques. To replicate the early water breakthrough observed during waterflooding, transmissibility contrasts are introduced into the simulation model, forcing viscous fingering effects. In the ensuing polymer flood, these transmissibility contrasts are reduced to replicate the restoration of injection conformance during polymer flooding, as indicated by a significant decrease in water cut. Later, transmissibility contrasts are reinstated to replicate a water surge event observed in one of the producing wells during polymer flooding. This event may represent decreased injection conformance from fracture overextension; its anticipated occurrence in the other production well is included in the final forecast. The definition of polymer retention in the simulator incorporates the tailing effect reported in laboratory studies; this tailing effect is useful to the simultaneous history match of producing water cut and produced polymer concentration. The top 24 best-matched simulation models produced at each stage of the history matching process are used to forecast oil recovery. The final forecast clearly demonstrates that polymer flooding significantly increases the heavy oil production for this field pilot compared to waterflooding alone. This exercise displays that a simulation model is only valid for prediction if flow behavior in the reservoir remains consistent with that observed during the history matched period. Critically, this means that a simulation model calibrated for waterflooding may not fully capture the benefits of an enhanced oil recovery process such as polymer flooding. Therefore, caution is recommended in using basic waterflood simulation models to scope potential enhanced oil recovery projects.

Keith, Cody Douglas↗
DOE OSTI2026

Development of a Test-Bed for Testing and Refining EarthEn’s Supercritical CO 2 Based Energy Storage System

EarthEn’s energy storage concept leverages supercritical carbon dioxide (sCO 2 ) as a working fluid and relies on compact, high-performance components operating at elevated pressures and temperatures. To accelerate component development and reduce technical risk prior to larger-scale demonstrations, Oak Ridge National Laboratory (ORNL) developed a 100 kW-scale sCO 2 test-bed under a Cooperative Research and Development Agreement with EarthEn (CRADA NO. NFE-24-10050). The objective of the work was to design and construct a flexible experimental facility capable of reproducing key thermodynamic state points and heat-transfer conditions relevant to EarthEn’s thermal energy storage (TES) cycle, with particular emphasis on enabling development and evaluation of next-generation heat exchangers and TES concepts. The test-bed consists of a closed-loop sCO 2 circulation system housed within an open-topped enclosure. In its as-installed configuration, dense-phase sCO 2 is recirculated through a printed circuit recuperator, an electrically heated section, a throttling device used to simulate turbine expansion, and a water-cooled printed circuit heat exchanger that rejects heat to the building chilled-water system before returning to the pump. The pump is driven by a variable frequency drive, enabling controlled adjustment of flow and operating point. A comprehensive instrumentation suite was integrated to support both safe operation and high-quality data collection. Installed sensors include Coriolis flow meters for sCO 2 flow rate and density, resistance temperature detectors and thermocouples distributed throughout the loop (including the heated section and key heat exchanger ports), and pressure transducers for absolute and differential pressure measurements. The facility was designed to support high-pressure (19 MPa nominal) and high-temperature (575°C nominal) operation with credited overpressure protection provided by a rupture disk. Nominal operating conditions were selected to support 100 kW-class testing while maintaining flexibility for non-heated and heated shakedown, control development, and future integration of advanced TES test sections. In parallel with facility development, a system-level thermal-hydraulic model was created using Modelica-based tools to support component sizing, anticipate performance over targeted test conditions, and establish a framework for future model calibration against experimental data. At the conclusion of the project performance period, the facility was in final assembly, and the pressure boundary was nearly completed. However, several practical challenges associated with high-pressure/high-temperature systems and specialized component procurement impacted schedule and prevented initial pump-driven operation and full commissioning within the available resources. This report documents the as-built design, operating capabilities, and instrumentation, and it summarizes key lessons learned related to heater fabrication and testing, first-of-a-kind assembly factors, specialty flange supply constraints, and fill pump corrective actions. Finally, it outlines a phased plan for future commissioning and experimental campaigns, including control and instrumentation shakedown, heater characterization, model calibration, and testing at state points representative of EarthEn’s TES cycle.

25 ENERGY STORAGE↗
DOE OSTI2022

Goal-oriented a-posteriori estimation of model error as an aid to parameter estimation

In this work, a Bayesian model calibration framework is presented that utilizes goal-oriented a-posterior error estimates in quantities of interest (QoIs) for classes of high-fidelity models characterized by PDEs. It is shown that for a large class of computational models, it is possible to develop a computationally inexpensive procedure for calibrating parameters of high-fidelity models of physical events when the parameters of low-fidelity (surrogate) models are known with acceptable accuracy. The main ingredients in the proposed model calibration scheme are goal-oriented a-posteriori estimates of error in QoIs computed using a so-called lower fidelity model compared to those of an uncalibrated higher fidelity model. The estimates of error in QoIs are used to define likelihood functions in Bayesian inversion analysis. A standard Bayesian approach is employed to compute the posterior distribution of model parameters of high-fidelity models. As applications, parameters in a quasi-linear second-order elliptic boundary-value problem (BVP) are calibrated using a second-order linear elliptic BVP. In a second application, parameters of a tumor growth model involving nonlinear time-dependent PDEs are calibrated using a lower fidelity linear tumor growth model with known parameter values.

A-posterior estimates↗