Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “calibrating distribution system 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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Better calibration of cloud parameterizations and subgrid effects increases the fidelity of the E3SM Atmosphere Model version 1

Abstract. Realistic simulation of the Earth's mean-state climate remains a major challenge, and yet it is crucial for predicting the climate system in transition. Deficiencies in models' process representations, propagation of errors from one process to another, and associated compensating errors can often confound the interpretation and improvement of model simulations. These errors and biases can also lead to unrealistic climate projections and incorrect attribution of the physical mechanisms governing past and future climate change. Here we show that a significantly improved global atmospheric simulation can be achieved by focusing on the realism of process assumptions in cloud calibration and subgrid effects using the Energy Exascale Earth System Model (E3SM) Atmosphere Model version 1 (EAMv1). The calibration of clouds and subgrid effects informed by our understanding of physical mechanisms leads to significant improvements in clouds and precipitation climatology, reducing common and long-standing biases across cloud regimes in the model. The improved cloud fidelity in turn reduces biases in other aspects of the system. Furthermore, even though the recalibration does not change the global mean aerosol and total anthropogenic effective radiative forcings (ERFs), the sensitivity of clouds, precipitation, and surface temperature to aerosol perturbations is significantly reduced. This suggests that it is possible to achieve improvements to the historical evolution of surface temperature over EAMv1 and that precise knowledge of global mean ERFs is not enough to constrain historical or future climate change. Cloud feedbacks are also significantly reduced in the recalibrated model, suggesting that there would be a lower climate sensitivity when it is run as part of the fully coupled E3SM. This study also compares results from incremental changes to cloud microphysics, turbulent mixing, deep convection, and subgrid effects to understand how assumptions in the representation of these processes affect different aspects of the simulated atmosphere as well as its response to forcings. We conclude that the spectral composition and geographical distribution of the ERFs and cloud feedback, as well as the fidelity of the simulated base climate state, are important for constraining the climate in the past and future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Study of the Cosmic Rays and Interstellar Medium in Local H I Clouds Using Fermi -LAT Gamma-Ray Observations

An accurate estimate of the interstellar gas density distribution is crucial to understanding the interstellar medium (ISM) and Galactic cosmic rays (CRs). To comprehend the ISM and CRs in a local environment, a study of the diffuse $γ$-ray emission in a midlatitude region of the third quadrant was performed. Specifically, the $γ$-ray data in the 0.1–25.6 GeV energy range of the Fermi Large Area Telescope (LAT) and other interstellar gas tracers such as the HI4PI survey data and the Planck dust thermal emission model were used, and the northern and southern regions were analyzed separately. The variation of the dust emission D em with the total neutral gas column density N H was studied in high dust temperature areas, and the N H /D em ratio was calibrated using $γ$-ray data under the assumption of a uniform CR intensity in the studied regions. The measured integrated $γ$-ray emissivities above 100 MeV are (1.58 ± 0.04) x 10 -26 photons s -1 sr -1 H-atom-1 and (1.59 ± 0.02) x 10-26 photons s-1 sr-1 H-atom -1 in the northern and southern regions, respectively, supporting the existence of a uniform CR intensity in the vicinity of the solar system. While most of the gas can be interpreted to be H I with a spin temperature of T S = 125 K or higher, an area dominated by optically thick H I with T S ~ 40 K was identified.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Decision support for United States—Canada energy integration is impaired by fragmentary environmental and electricity system modeling capacity

The renewable energy transition is leading to increased electricity trade between the United States and Canada, with Canadian hydropower providing firm lower-carbon power and buffering variability of wind and solar generation in the U.S. However, long-term power purchase agreements and transborder transmission projects are controversial, with two of four proposed transmission lines between Quebec, Canada and the northeast U.S. cancelled since 2018. Here, we argue that controversies are exacerbated by a lack of open-source data and tools to understand tradeoffs of new hydropower generation and transmission infrastructure in comparison to alternatives. This gap includes impacts that incremental transmission and generation projects have on the economics of the entire system, for example, how new transmission projects affect exports to existing markets or incentivize new generation. We identify priority areas for data synthesis and model development, such as integrating linked hydropower and hydrologic interactions in energy system models and openly releasing (by utilities) or back-calculating (by researchers) hydropower generation and operational parameters. Publicly available environmental (e.g. streamflow, precipitation) and techno-economic (e.g. costs, reservoir size,) data can be used to parameterize freely usable and extensible models. Existing models have been calibrated with operational data from Canadian utilities that are not publicly available, limiting the range of scientific and commercial questions these tools have been used to answer and the range of parties that have been involved. Studies conducted using highly resolved, national-scale public data exist in other countries, notably, the United States, and demonstrate how greater transparency and extensibility can drive industry action. Improved data availability in Canada could facilitate approaches that (1) increase participation in decarbonization planning by a broader range of actors; (2) allow independent characterizations of environmental, health, and economic outcomes of interest to the public; and (3) identify decarbonization pathways consistent with community values.

13 HYDRO ENERGY↗
DOE OSTI2023

Model certainty in cellular network-driven processes with missing data

Mathematical models are often used to explore network-driven cellular processes from a systems perspective. However, a dearth of quantitative data suitable for model calibration leads to models with parameter unidentifiability and questionable predictive power. Here we introduce a combined Bayesian and Machine Learning Measurement Model approach to explore how quantitative and non-quantitative data constrain models of apoptosis execution within a missing data context. We find model prediction accuracy and certainty strongly depend on rigorous data-driven formulations of the measurement, and the size and make-up of the datasets. For instance, two orders of magnitude more ordinal (e.g., immunoblot) data are necessary to achieve accuracy comparable to quantitative (e.g., fluorescence) data for calibration of an apoptosis execution model. Notably, ordinal and nominal (e.g., cell fate observations) non-quantitative data synergize to reduce model uncertainty and improve accuracy. Finally, we demonstrate the potential of a data-driven Measurement Model approach to identify model features that could lead to informative experimental measurements and improve model predictive power.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Enabling machine learning-ready HPC ensembles with Merlin

With the growing complexity of computational and experimental facilities, many scientific researchers are turning to machine learning (ML) techniques to analyze large scale ensemble data. With complexities such as multi-component workflows, heterogeneous machine architectures, parallel file systems, and batch scheduling, care must be taken to facilitate this analysis in a high performance computing (HPC) environment. Here, we present Merlin, a workflow framework to enable large ML-friendly ensembles of scientific HPC simulations. By augmenting traditional HPC with distributed compute technologies, Merlin aims to lower the barrier for scientific subject matter experts to incorporate ML into their analysis. As a producer–consumer workflow model, Merlin enables multi-machine, cross-batch job, dynamically allocated yet persistent workflows capable of utilizing surge-compute resources. Key features of Merlin are a flexible HPC-centric interface, low per-task overhead, multi-tiered fault recovery, and a hierarchical sampling algorithm that allows for $\mathscr{O}$(N) task execution and $\mathscr{O}$(N ln N) task queuing to ensembles of millions of tasks. In addition to Merlin’s design, we test the algorithm’s performance in an HPC center and demonstrate the ability to enqueue 40 million simulations in 100 s, with a 30 millisecond per-task overhead that is independent of ensemble size. Finally, we describe some example applications that Merlin has enabled on leadership-class HPC resources, such as the ML-augmented optimization of nuclear fusion experiments and the calibration of infectious disease models to study the progression of and possible mitigation strategies for COVID-19.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Addressing Experimental Challenges in Probing Dark Energy with Accuracy and Precision with the Rubin Observatory Legacy Survey of Space & Time (LSST)

With support from this award, the Stanford group under the leadership of PI Patricia Burchat focused on pixel-level algorithms that are key to the investigation of dark energy with the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). The work was carried out through strong engagement in the DOE-supported LSST Dark Energy Science Collaboration (DESC), for which the PI served as Deputy Spokesperson (07/2019 - 07/2021) and Technical Coordinator (07/2021 - 07/2023). Members of the Burchat group are Full DESC members and are among scientists who have been accepted to join the Rubin Observatory System Integration, Test, and Commissioning (SIT-Com) effort. The specific intellectual focus of the group was the optimization of gravitational lensing as a probe of the evolution of the distribution of dark matter in space and time to learn about the nature of dark energy. The Stanford group pursued an integrated research program to understand in detail, and then remediate or calibrate, systematic effects that could limit the sensitivity of the LSST to cosmic shear as a dark energy probe -- particularly those related to our understanding of the atmospheric point spread function (PSF), challenges associated with blended galaxy images, and impacts of residual camera artifacts. Key to this work was developing, enhancing, and validating the fidelity of simulations that play a key role in the development of our understanding of subtle systematic biases and then using these simulations to develop techniques to characterize and mitigate biases, and associated algorithms for implementation in the LSST or DESC science pipelines. The group’s work included contributions to the open-source GalSim simulation package, which is widely used in the community; simulations and studies of spatial correlations of atmospheric PSFs; studies of temporal correlations of the atmospheric PSF; impacts of chromatic PSFs on cosmic shear calibration; continued studies of anisotropic Gaussian process interpolation of PSF models across the focal plane; and measuring and mitigating impacts of residual camera artifacts on cosmic shear. A primary impact of this work is that it increases the scientific discovery potential of the DOE and NSF investments in the Rubin Observatory Legacy Survey of Space \& Time, by addressing challenges in interpreting astronomical images through more accurate descriptions of the blurring of images due to the atmosphere. Specifically, this work impacts our ability to understand the fundamental nature of dark energy.

Burchat, Patricia Rose↗
DOE OSTI2020

Sensitivity of meteorological-forcing resolution on hydrologic variables

Projecting the spatiotemporal changes in water resources under a no-analog future climate requires physically based integrated hydrologic models which simulate the transfer of water and energy across the earth's surface. These models show promise in the context of unprecedented climate extremes given their reliance on the underlying physics of the system as opposed to empirical relationships. However, these techniques are plagued by several sources of uncertainty, including the inaccuracy of input datasets such as meteorological forcing. These datasets, usually derived from climate models or satellite-based products, are typically only resolved on the order of tens to hundreds of kilometers, while hydrologic variables of interest (e.g., discharge and groundwater levels) require a resolution at much smaller scales. In this work, a high-resolution hydrologic model is forced with various resolutions of meteorological forcing (0.5 to 40.5 km) generated by a dynamical downscaling analysis from the regional climate model Weather Research and Forecasting (WRF). The Cosumnes watershed, which spans the Sierra Nevada and Central Valley interface of California (USA), exhibits semi-natural flow conditions due to its rare undammed river basin and is used here as a test bed to illustrate potential impacts of various resolutions of meteorological forcing on snow accumulation and snowmelt, surface runoff, infiltration, evapotranspiration, and groundwater levels. Results show that the errors in spatial distribution patterns impact land surface processes and can be delayed in time. Localized biases in groundwater levels can be as large as 5–10 m and 3 m in surface water. Most hydrologic variables reveal that biases are seasonally and spatially dependent, which can have serious implications for model calibration and ultimately water management decisions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Hybrid Solar System (Final Scientific/Technical Report)

GTI Energy (GTI) teamed with the University of California at Merced (UCM) to scaleup the hybrid solar system (HSS) technology for demonstrating its performance at the US Gypsum (USG) plant in Plaster City, California. The technology integrates two-stage concentrating solar collector with matching particle thermal transport and storage (TSS) system to deliver cost-effective, and on-demand distributed high temperature industrial process heat up to 600°C with solar thermal, in this case to a gypsum kettle, to reduce its fuel use and carbon footprint. Current solar technologies, which reach these temperatures, are not distributable (towers) or cost-effective (dish). The research team developed a conceptual system design for host site retrofit, including preliminary heat balance, process flow diagram, particle to process heat exchanger and equipment placements at the site. Subsequently, parallel efforts were carried out at UCM to design, build and test a 12 m long commercial scale prototype concentrating thermal-only collector system and at GTI to design, build and test a matching 650°C capable particle TTS system. The nominal 50 kWth collector consists of a parabolic trough and three 4 m long two-stage receivers in series. Prior to on-sun testing, a 4 m long receiver was fabricated and successfully tested at 650 °C in a laboratory setting for 100 hrs of continuous operation showing less than 15% radiation loss. A 7 m wide x 17 m long parabolic trough was then installed at UCM for on-sun testing of the 12 m long receiver, and concurrently several 4 m long receivers were built. The optics of the parabolic trough were calibrated, and on-sun test were carried out on 12 m long receivers. During tests, the intense solar radiation (53x) caused the absorber tubes in the receivers to bend, reducing the overall optical efficiency. To address the bending issue, a self-consistent algorithm that includes ray tracing, thermal and deformation models was developed to perform thermal stress analysis on absorbers for parabolic solar collectors. Results obtained with this algorithm showed a dramatic rise in deformation as absorber tube length increases. A combined efficiency parameter that includes the occluded area for the mounts was developed to obtain an optimized tube length obtained. Based on the results, a length of 2.7 m for the absorber + 0.2 m for the coupler was chosen to minimize any bending and optimize optical efficiency while maintaining ease of mounting. The associated particle TTS system was designed, built and successfully tested at GTI. It includes storage, receiving and lock hoppers and piping that simulates the transfer of captured solar energy to an actual industrial furnace. Tests over 77 charge-discharge cycles demonstrated <2% particle degradation, with no problematic particle accumulations and no flow interruptions. The piping pressure drop was about 5 psi. The team also worked with Stanley Consultants (Stanley) to prepare conceptual and preliminary engineering packages to facilitate follow-on development and commercialization efforts. These include process and instrumentation diagram’s (P&ID’s), general arrangements, electrical one-line, project definitions document, equipment data sheets, schedule, and construction cost estimate for 2 MWth system. Updated HSS technology commercialization and customer engagement plans and detailed costs and evaluated market trade-offs and manufacturing.

03 NATURAL GAS↗
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↗
DOE OSTI2024

Simulation-based inference for parameter estimation of complex watershed simulators

High-resolution, spatially distributed process-based (PB) simulators are widely employed in the study of complex catchment processes and their responses to a changing climate. However, calibrating these PB simulators using observed data remains a significant challenge due to several persistent issues, including the following: (1) intractability stemming from the computational demands and complex responses of simulators, which renders infeasible calculation of the conditional probability of parameters and data, and (2) uncertainty stemming from the choice of simplified representations of complex natural hydrologic processes. Here, we demonstrate how simulation-based inference (SBI) can help address both of these challenges with respect to parameter estimation. SBI uses a learned mapping between the parameter space and observed data to estimate parameters for the generation of calibrated simulations. To demonstrate the potential of SBI in hydrologic modeling, we conduct a set of synthetic experiments to infer two common physical parameters – Manning's coefficient and hydraulic conductivity – using a representation of a snowmelt-dominated catchment in Colorado, USA. We introduce novel deep-learning (DL) components to the SBI approach, including an “emulator” as a surrogate for the PB simulator to rapidly explore parameter responses. We also employ a density-based neural network to represent the joint probability of parameters and data without strong assumptions about its functional form. While addressing intractability, we also show that, if the simulator does not represent the system under study well enough, SBI can yield unreliable parameter estimates. Approaches to adopting the SBI framework for cases in which multiple simulator(s) may be adequate are introduced using a performance-weighting approach. The synthetic experiments presented here test the performance of SBI, using the relationship between the surrogate and PB simulators as a proxy for the real case.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

X-ray Inspection Model Validation with Physical Dosimetry

A non-functional printed circuit board assembly was developed using ULTIboard and a simple computational transport model was created for use with the Monte Carlo software MCNP. EBT 3 film radiography was used to compare with Monte Carlo simulation to validate the board representation. The computational transport model was then revised to include connection pins, solder balls, highly attenuating internal structures, and copper trace distribution. Results from the revised model were compared to the EBT 3 films to observe improvements to dose profiles. It was found that this method was useful in verifying the placement of components, as the dose profiles were observed to follow the same trends. The experiment was then repeated using XRQA 2 films to achieve the same level of contrast with 1% the dose of EBT 3 films. It was found that high contrast may be achieved using these films to identify major issues with the model geometry, at a cost of dose profile accuracy. A second validation method was applied to the model using 37 CaF 2 thermoluminescent dosimeters (TLDs). TLD measurements were compared with the simplified and complex transport models to identify the features that have the greatest impact on simulation accuracy. The TLD calibration to CaF 2 was found to be accurate within 5.6%, while calibration to dose in Si was found to be accurate within 4.7%. Finally, it was observed that the accurate representation of solder balls and proper modeling of highly attenuating internal structures had the greatest impact on simulation accuracy.

61 RADIATION PROTECTION AND DOSIMETRY↗
DOE OSTI2023

Cummins R-SOFC System Development

The overall purpose of this project was to reduce the Reversible-Solid Oxide Fuel Cell (R-SOFC) system cost by developing two technologies, an improved cell design and the incorporation of an ejector in the fuel recycle loop instead of a blower. A Simulink model of the baseline SOFC system was developed and calibrated with experimental test data. The R-SOFC system model was built by integrating GT Suite developed models of the steam generation components into the baseline Simulink SOFC system model. The ability to run the stack in SOEC operating mode was also added to the model. The system model was used to explore the ability of the R-SOFC system to meet operational constraints on Steam/Carbon ratio and H2 concentration on the fuel side electrode. A CFD ejector model was developed and used to explore a range of ejector design parameters, leading to the final ejector design that was prototyped for testing. A prototype steam ejector was first tested in a laboratory environment using room temperature air. The steam ejector was subsequently tested using the full hot recycle loop with all relevant heat exchangers and steam generation components. The test conditions utilized temperatures, pressures, and flow rates expected in an R-SOFC application. Throughout the experimental testing work, ejector performance test data was used to improve and then validate the CFD ejector model. A CFD cell model was developed and used to optimize thermal gradients, voltage, and cost of a new cell substrate design. A CFD comparison of co-flow and cross-flow cell designs informed the decision to use a co-flow design for the new cell substrate. Multiple rounds of CFD simulation were used to improve the cell design to minimize the variation in air and fuel distribution across cell channels and to minimize the variation in air and fuel distribution across different cells in the stack. A few prototypes of the new cell substrate design were produced and validated in a laboratory environment by thermally spraying and verifying that they met established manufacturing specifications. The cell manufacturing process was adjusted in order to bring these metrics within acceptable tolerances. Cummins’ internal calculations show that the new cell design reduces cost ~50% compared to the baseline cell, while the ejector + superheater/boiler concept reduces cost of the recycle loop by ~40%. The impact of these cost reductions on the cost of producing H2 will depend on the specific system where they are applied. Therefore, a Techno-Economic Analysis was completed using system cost as a variable, and showing how the NREL Current and Future system costs translate into H2 production cost.

Henrichsen, Lars↗
DOE OSTI2020

Wireless Sensor Modalities at a Nuclear Plant Site to Collect Vibration Data

One of the major contributors to the total operating costs of domestic nuclear fleet of reactors today is the operation and maintenance (O&M) costs. These include labor-intense preventive maintenance programs involving manually-performed inspection, calibration, testing, and maintenance of plant assets at periodic frequency and time-based replacement of assets at periodic frequency, irrespective of their conditions. This has resulted in a labor-centric business model to achieve high capacity factors. To build an optimal maintenance program, it’s time to transition from this labor-centric business model to a technology-centric business model. Fortunately, there are technologies (advanced sensor, data analytics, and risk assessment methodologies) that will support this transition. The technology-centric business model will result in significant plant life extension and reduction of time-based maintenance activities. This will drive down O&M costs as labor is a rising cost and technology is a declining cost. This approach will lay the foundation for real-time condition assessment of plant assets, allowing condition-based maintenance to enhance plant safety, reliability, and economics of operation. The goal of this project is to address challenges in the area of digital monitoring, i.e., the application of advanced sensor technologies (particularly wireless sensor technologies) and science-based data analytic capabilities to advance online monitoring and predictive maintenance in nuclear plants to improve plant performance (efficiency gain and economic competitiveness). To achieve the project goal, in partnership with Exelon Generating Company (Exelon), researchers from Idaho National Laboratory (INL) and Oak Ridge National Laboratory (ORNL) are performing research and development (R&D) to demonstrate application of wireless sensors using the distributed antenna system and advanced data analytics to achieve predictive maintenance. In the report, wireless vibration sensors, vibration data and its indicator are described. The wireless vibration sensors presented in this report support three types of wireless communication, namely, Wi-Fi, cellular, and 900 MHz. These wireless vibration sensors are considered by partner plant site for installation on plant asset to enable online vibration monitoring to replace periodic measurements. These vibration data along with other plant process data will be utilized to develop diagnostic and prognostic models.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

AI-assisted detector design for the EIC (AID(2)E)

Artificial Intelligence is poised to transform the design of complex, large-scale detectors like ePIC at the future Electron Ion Collider. Featuring a central detector with additional detecting systems in the far forward and far backward regions, the ePIC experiment incorporates numerous design parameters and objectives, including performance, physics reach, and cost, constrained by mechanical and geometric limits. This project aims to develop a scalable, distributed AI-assisted detector design for the EIC (AID(2)E), employing state-of-the-art multiobjective optimization to tackle complex designs. Supported by the ePIC software stack and using G EANT 4 simulations, our approach benefits from transparent parameterization and advanced AI features. The workflow leverages the PanDA and iDDS systems, used in major experiments such as ATLAS at CERN LHC, the Rubin Observatory, and sPHENIX at RHIC, to manage the compute intensive demands of ePIC detector simulations. Tailored enhancements to the PanDA system focus on usability, scalability, automation, and monitoring. Ultimately, this project aims to establish a robust design capability, apply a distributed AI-assisted workflow to the ePIC detector, and extend its applications to the design of the second detector (Detector-2) in the EIC, as well as to calibration and alignment tasks. Additionally, we are developing advanced data science tools to efficiently navigate the complex, multidimensional trade-offs identified through this optimization process.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Probabilistic projections of the Amery Ice Shelf catchment, Antarctica, under conditions of high ice-shelf basal melt

Abstract. Antarctica's Lambert Glacier drains about one-sixth of the ice from the East Antarctic Ice Sheet and is considered stable due to the strong buttressing provided by the Amery Ice Shelf. While previous projections of the sea-level contribution from this sector of the ice sheet have predicted significant mass loss only with near-complete removal of the ice shelf, the ocean warming necessary for this was deemed unlikely. Recent climate projections through 2300 indicate that sufficient ocean warming is a distinct possibility after 2100. This work explores the impact of parametric uncertainty on projections of the response of the Lambert–Amery system (hereafter “the Amery sector”) to abrupt ocean warming through Bayesian calibration of a perturbed-parameter ice-sheet model ensemble. We address the computational cost of uncertainty quantification for ice-sheet model projections via statistical emulation, which employs surrogate models for fast and inexpensive parameter space exploration while retaining critical features of the high-fidelity simulations. To this end, we build Gaussian process (GP) emulators from simulations of the Amery sector at a medium resolution (4–20 km mesh) using the Model for Prediction Across Scales (MPAS)-Albany Land Ice (MALI) model. We consider six input parameters that control basal friction, ice stiffness, calving, and ice-shelf basal melting. From these, we generate 200 perturbed input parameter initializations using space filling Sobol sampling. For our end-to-end probabilistic modeling workflow, we first train emulators on the simulation ensemble and then calibrate the input parameters using observations of the mass balance, grounding line movement, and calving front movement with priors assigned via expert knowledge. Next, we use MALI to project a subset of simulations to 2300 using ocean and atmosphere forcings from a climate model for both low- and high-greenhouse-gas-emission scenarios. From these simulation outputs, we build multivariate emulators by combining GP regression with principal component dimension reduction to emulate multivariate sea-level contribution time series data from the MALI simulations. We then use these emulators to propagate uncertainty from model input parameters to predictions of glacier mass loss through 2300, demonstrating that the calibrated posterior distributions have both greater mass loss and reduced variance compared to the uncalibrated prior distributions. Parametric uncertainty is large enough through about 2130 that the two projections under different emission scenarios are indistinguishable from one another. However, after rapid ocean warming in the first half of the 22nd century, the projections become statistically distinct within decades. Overall, this study demonstrates an efficient Bayesian calibration and uncertainty propagation workflow for ice-sheet model projections and identifies the potential for large sea-level rise contributions from the Amery sector of the Antarctic Ice Sheet after 2100 under high-greenhouse-gas-emission scenarios.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

An Integrated Multiscale Experimental-Numerical Analysis on Reconsolidation of Salt-Clay Mixture for Disposal of Heat-Generating Waste (Final NEUP Technical Report)

The overall purpose of this research is to improve understanding of THMC coupling effect on the reconsolidation of granular (or crushed) salt-clay mixture used for seal systems of shafts and drifts in salt repositories. This proposed work is partially motivated by the recent work on the Waste Isolation Pilot Plant (WIPP) that shows the promising sealing capability of clay-salt mixture compared to crushed salt. In particular, primary emphasis is to develop a fully integrated multiscale experiment-numerical study to determine and explain what leads to the superior sealing ability of the clay-salt mixture. These research activities are designed to seek further understanding of (1) why clay additives may enhance the fluid trapping and (2) whether this flow barrier effect may prevail under different combinations of temperature, confining pressure, deviatoric stress and other foreseeable environmental factors. If successful, this enhanced flow trapping ability of the seal provides significant improvement to the seal and repository performance and therefore make the repository safer in the long-term. The experiment component includes microstructural investigation and macroscopic tests on a reconsolidated salt-clay mixture. In the microstructural study, the goal is to (1) characterize microscopic distributions of distinct phases (e.g., clay, salt crystal boundaries, trapped brine, and pore) to examine the connectivity of the pore network inside the salt-clay mixture with different amounts of clay additive and moisture content and (2) analyze multiscale imaging data to reconstruct the polycrystalline microstructures for numerical simulations. Meanwhile, macroscopic tests are performed to analyze how clay alters the failure/creep mechanisms in the salt-clay mixture. Microscopic and macroscopic experimental observations will both be used to calibrate and validate a multiscale model that explicitly simulates the capillary and multiphase flow in the connected pores and the deformation due to the presence of intra-crystalline brine at the pore scale via a new polyhedral discrete element–lattice Boltzmann method (DEM-LBM) coupling model. The pore-scale simulations are homogenized via an upscaling procedure that converts pore-scale information (e.g. force exerted on grain boundary, sliding, pressure-solution) to continuum measures (e.g. Cauchy stress, Darcy’s flow) at each integration point in the macroscopic multiphase TMHC model. This multiscale scheme will allow coupling be- tween high-fidelity simulations of brine-salt-clay interaction and the macroscopic TMHC model. The multiscale model helps the understanding of how the trapped brine inclusion affects the pressure-solution mechanism with the presence of clay and moisture. This work brings new insight into the sealing capacity of salt-clay mixture under elevated temperature over a long period of time - a key to evaluating the potential of salt-clay mixture usage for salt repositories.

42 ENGINEERING↗
DOE OSTI2022

The WRF-Solar Ensemble Prediction System: Development, Test, and Validation

Providing reliable probabilistic solar radiation information is needed to improve management of the uncertainty and variability of solar generation. Thus, guidance on how to develop skillful and accurate ensemble forecasts is essential and it will ultimately contribute to integration of high amounts of solar energy on the grid. A team from the National Renewable Energy Laboratory and the National Center for Atmospheric Research had been collaborating to develop the WRF-Solar ensemble prediction system (WRF-Solar EPS) in the past three years to produce probabilistic solar irradiance forecasts and better predict solar energy by quantifying forecast uncertainty. The WRF-Solar EPS basically generates ensemble members for solar irradiance based on stochastic perturbations to provide intraday and day-ahead probabilistic forecasts. This study will present main research steps in developing the WRF-Solar EPS including: (a) tangent linear analysis for identifying key input variables of six WRF-Solar modules significantly related to predicting of cloud and solar irradiance, (b) combining stochastic perturbation technique with the WRF-Solar model, and (c) ensemble calibration method to decrease error and uncertainty of ensemble-based solar forecasts. The capability of WRF-Solar EPS is now updated to the most recent version of standard WRF model. This presentation will summarize comprehensive results from the evaluation of forecasts against the National Solar Radiation Data Base as well as ground-measured observations. Moreover, we will introduce the user's guide for WRF-Solar EPS (e.g., parameters to configure stochastic perturbations) and future extension of this research.

day-ahead forecast↗
DOE OSTI2023

Correlated Trajectory Uncertainty for Adaptive Sequential Decision Making

One of the great challenges with decision making tasks on real world systems is the fact that data is sparse and acquiring additional data is expensive. In these cases, it is often crucial to make a model of the environment to assist in making decisions. At the same time, limited data means that learned models are erroneous, making it just as important to equip the model with good predictive uncertainties. In the context of learning sequential decision making policies, these uncertainties can prove useful for informing which data to collect for the greatest improvement in policy performance \citep{mehta2021experimental, mehta2022exploration} or informing the policy about unsure regions of state and action space to avoid during test time \citep{yu2020mopo}. Additionally, assuming that realistic samples of the environment can be drawn, an adaptable policy can be trained that attempts to make optimal decisions for any given possible instance of the environment \citep{ghosh2022offline, chen2021offline}. In this work, we examine the so-called ``probabilistic neural network'' (PNN) model that is ubiquitous in model-based reinforcement learning (MBRL) works. We argue that while PNN models may have good marginal uncertainties, they form a distribution of non-smooth transition functions. Not only are these samples unrealistic and may hamper adaptability, but we also assert that this leads to poor uncertainty estimates when predicting multiple step trajectory estimates. To address this issue, we propose a simple sampling method that can be implemented on top of pre-existing models.We evaluate our sampling technique on a number of environments, including a realistic nuclear fusion task, and find that, not only do smooth transition function samples produce more calibrated uncertainties, but they also lead to better downstream performance for an adaptive policy.

Offline Reinforcement Learning↗