SEARCH · Engineering Papers

Results for “surface complexation modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Tensorized Interior Radiative Heat Transfer for a Scalable and Calibrated Building Energy Simulator

Building energy simulation is a critical tool for developing and testing advanced control strategies, such as Reinforcement Learning (RL), to provide demand flexibility and affordable energy costs. The recently introduced Smart Buildings Control Suite (sbsim) provides a lightweight, scalable, and data-calibrated simulation environment based on a 2D finite-difference model. However, the initial model primarily focused on conductive and convective heat transfer, neglecting the significant impact of long-wave radiative heat exchange between interior surfaces. This paper presents a significant extension to the sbsim framework by incorporating a physically-grounded model for interior radiative heat transfer. Our primary contribution is the development and integration of a fully tensorized radiative heat transfer module, which preserves the computational efficiency and scalability of the original simulator. This was achieved by developing a pipeline for view factor calculation, including an algorithm to identify directly seeing surfaces within complex floor plans, and formulating the net radiation equations for efficient execution on modern hardware accelerators. We validate the numerical accuracy of our tensorized implementation by comparing its results against a traditional iterative approach, demonstrating identical outcomes. This enhancement increases the physical fidelity of sbsim, enabling more accurate training of RL agents for building energy optimization.

Ham, Sang woo↗

DOE OSTI2023

Improving snow albedo modeling in the E3SM land model (version 2.0) and assessing its impacts on snow and surface fluxes over the Tibetan Plateau

With the highest albedo of the land surface, snow plays a vital role in Earth's surface energy budget and water cycle. Snow albedo is primarily controlled by snow grain properties (e.g., size and shape) and light-absorbing particles (LAPs) such as black carbon (BC) and dust. The mixing state of LAPs in snow also has impacts on LAP-induced snow albedo reduction and surface radiative forcing (RF). However, most land surface models assume that snow grain shape is spherical and LAPs are externally mixed with the snow grains. This study improves the snow radiative transfer model in the Energy Exascale Earth System Model version 2.0 (E3SM v2.0) Land Model (ELM v2.0) by considering non-spherical snow grain shapes (i.e., spheroid, hexagonal plate, and Koch snowflake) and internal mixing of dust–snow, and it systematically evaluates the impacts on the surface energy budget and water cycle over the Tibetan Plateau (TP). A series of ELM simulations with different treatments of snow grain shape, mixing state of BC–snow and dust–snow, and sub-grid topographic effects (TOP) on solar radiation are performed. Compared with two remote sensing snow products derived from the Moderate Resolution Imaging Spectroradiometer, the control ELM simulation (ELM_Control) with the default configurations of spherical snow grain shape, internal mixing of BC–snow, external mixing of dust–snow, and without TOP as well as the ELM simulation with new model features (ELM_New) can both capture the overall snow distribution reasonably. Additionally, ELM_New overall shows smaller biases in snow cover fraction than ELM_Control in spring when snowmelt is important for water management. The estimated LAP-induced RF in ELM_New ranges from 0 to 19.3 W m -2 with the area-weighted average value of 1.5 W m -2 that is comparable to the reported values in existing studies. The Koch snowflake shape, among other non-spherical shapes, shows the largest difference from the spherical shape in spring when snow processes related to the surface energy budget and water cycle have high importance. The impacts of the mixing state of LAP in snow are smaller than the shape effects and depend on snow grain shape. Compared to external mixing, internal mixing of LAP–snow can lead to larger snow albedo reduction and snowmelt, which further affect the surface energy budget and water cycle. The individual contributions of non-spherical snow shape, mixing state of LAP–snow, and local topography impacts on the snow and surface fluxes have different signs and magnitudes, and their combined effects may be negative or positive due to complex and nonlinear interactions among the factors. Overall, the changes in net solar radiation in spring due to individual and combined effects range from -28.6 to 16.9 W m -2 and -29.7 to 12.2 W m -2 , respectively. This study advances understanding of the role of snow grain shape and mixing state of LAP–snow in land surface processes and offers guidance for improving snow simulations and RF estimates in Earth system models under climate change.

58 GEOSCIENCES↗

DOE OSTI2022

Investigation of Coupled Processes in Fractures and the Bordering Matrix via a Micro‐Continuum Reactive Transport Model

Abstract In multi‐mineral fractured rocks, the altered porous layer on the fracture surface resulting from preferential dissolution of the fast‐reacting minerals can have profound impacts on subsequent chemical‐physical alteration of the fractures. This study adopts the micro‐continuum approach to provide further understanding of reactive transport processes in the altered layer (AL), and mass exchanges with the bordering matrix and fracture. The modeling framework couples the Darcy‐Brinkman‐Stokes (DBS) solver in COMSOL Multiphysics and the geochemical modeling capability of CrunchFlow. Three‐dimensional steady state simulations with systematically varied chemical‐physical parameters of the AL were performed to examine the impacts of individual factors and processes. Our simulation results confirm previous observations that dissolution of the fast‐reacting mineral (i.e., calcite) is largely controlled by diffusion across the AL. We also show that dissolution of the slow‐reacting mineral (i.e., dolomite), which controls AL development and fracture enlargement, increases with surface area and has a complex dependence on different local rate‐limiting processes. In particular, advection can result in evident spatial variations in the local dissolution rates of dolomite, although it does not affect the bulk chemistry significantly. The difference in the spatial patterns between simulations with and without advection in the AL is more noticeable in the locations with smaller apertures, with up to 20% difference in local reaction rates. Therefore, it is important to include a full depiction of advection, diffusion, and reactions for accurately capturing local dynamics that control long‐term fracture evolution.

58 GEOSCIENCES↗

DOE OSTI2021

Uranium (VI) Adsorbate Structures on Portlandite [Ca(OH)2] Type Surfaces Determined by Computational Modelling and X-ray Absorption Spectroscopy

Portlandite [Ca(OH) 2 ] is a potentially dominant solid phase in the high pH fluids expected within the cementitious engineered barriers of Geological Disposal Facilities (GDF). This study combined X-ray Absorption Spectroscopy with computational modelling in order to provide atomic-scale data which improves our understanding of how a critically important radionuclide (U) will be adsorbed onto this phase under conditions relevant to a GDF environment. Such data are fundamental for predicting radionuclide mass transfer. Surface coordination chemistry and speciation of uranium with portlandite [Ca(OH) 2 ] under alkaline groundwater conditions (ca. pH 12) were determined by both in situ and ex situ grazing incidence extended X-ray absorption fine structure analysis (EXAFS) and by computational modelling at the atomic level. Free energies of sorption of aqueous uranyl hydroxides, [UO 2 (OH) n ] 2–n (n = 0–5) with the (001), (100) and (203) or (101) surfaces of portlandite are predicted from the potential of mean force using classical molecular umbrella sampling simulation methods and the structural interactions are further explored using fully periodic density functional theory computations. Although uranyl is predicted to only weakly adsorb to the (001) and (100) clean surfaces, there should be significantly stronger interactions with the (203/101) surface or at hydroxyl vacancies, both prevalent under groundwater conditions. The uranyl surface complex is typically found to include four equatorially coordinated hydroxyl ligands, forming an inner-sphere sorbate by direct interaction of a uranyl oxygen with surface calcium ions in both the (001) and (203/101) cases. In contrast, on the (100) surface, uranyl is sorbed with its axis more parallel to the surface plane. The EXAFS data are largely consistent with a surface structural layer or film similar to calcium uranate, but also show distinct uranyl characteristics, with the uranyl ion exhibiting the classic dioxygenyl oxygens at 1.8 Å and between four and five equatorial oxygen atoms at distances between 2.28 and 2.35 Å from the central U absorber. These experimental data are wholly consistent with the adsorbate configuration predicted by the computational models. These findings suggest that, under the strongly alkaline conditions of a cementitious backfill engineered barrier, there would be significant uptake of uranyl by portlandite to inhibit the mobility of U(VI) from the near field of a geological disposal facility.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Nonlinear ionization dynamics of hot dense plasma observed in a laser-plasma amplifier

Understanding the behaviour of matter under conditions of extreme temperature, pressure, density and electromagnetic fields has profound effects on our understanding of cosmologic objects and the formation of the universe. Lacking direct access to such objects, our interpretation of observed data mainly relies on theoretical models. However, such models, which need to encompass nuclear physics, atomic physics and plasma physics over a huge dynamic range in the dimensions of energy and time, can only provide reliable information if we can benchmark them to experiments under well-defined laboratory conditions. Due to the plethora of effects occurring in this kind of highly excited matter, characterizing isolated dynamics or obtaining direct insight remains challenging. High-density plasmas are turbulent and opaque for radiation below the plasma frequency and allow only near-surface insight into ionization processes with visible wavelengths. Here, the output of a high-harmonic seeded laser-plasma amplifier using eight-fold ionized krypton as the gain medium operating at a 32.8 nm wavelength is ptychographically imaged. A complex-valued wavefront is observed in the extreme ultraviolet (XUV) beam with high resolution. Ab initio spatio-temporal Maxwell–Bloch simulations show excellent agreement with the experimental observations, revealing overionization of krypton in the plasma channel due to nonlinear laser-plasma interactions, successfully validating this four-dimensional multiscale model. This constitutes the first experimental observation of the laser ion abundance reshaping a laser-plasma amplifier. The presented approach shows the possibility of directly modelling light-plasma interactions in extreme conditions, such as those present during the early times of the universe, with direct experimental verification.

36 MATERIALS SCIENCE↗

DOE OSTI2023

Dynamic In‐Plane Heterogeneous and Inverted Response of Graphite to Fast Charging and Discharging Conditions in Lithium‐Ion Pouch Cells

Solutions for improving fast charging of lithium‐ion batteries have largely focused on alleviating through‐plane lithiation gradients while little is understood about in‐plane heterogeneities and how to resolve them. Herein, high‐speed synchrotron X‐ray diffraction (XRD) resolves graphite lithiation spatially and temporally during 6 C charging and 2 C discharging. At every point during operation, considerable differences in the state of lithiation across the pouch cell are present. Some regions are more responsive to operation than others, reaching full lithiation early during charge and full delithiation during discharge. Other regions within the cell never fully delithiate during discharge, despite a prolonged voltage hold at 2.8 V. Using time‐resolved XRD data, the calculated local current density (mA cm −2 ) at the graphite surface shows an unexpected occurrence of local inverted current densities where regions of graphite are observed to delithiate during charging and lithiate during discharging. A pseudo‐3D model is developed for the graphite electrode with spatially varying microstructural tortuosity to show how microstructural heterogeneity could influence spatial charge dynamics. The model could not predict the complex in‐plane charge behavior observed within the cell. Consequently physics‐based charging protocols based on homogeneous electrode assumptions may underestimate the local variations in charge dynamics and occurrence of lithium plating.

25 ENERGY STORAGE↗

DOE OSTI2025

Online thermal profile prediction for large format additive manufacturing: A hybrid CNN-LSTM based approach

Large format additive manufacturing (LFAM) is an advanced 3D printing technique that efficiently fabricates large-scale components through a layer-by-layer extrusion and deposition process. Accurate surface layer temperature monitoring is essential to prevent manufacturing failures and ensure final product quality. Traditional physics-based offline approaches for simulating thermal behavior are often inefficient and complex, posing challenges on real-time, in-situ monitoring. Here, to address this, we propose a data-driven hybrid CNN-LSTM model to predict sequential thermal images of arbitrary length using real-time infrared thermal imaging. In this approach, a Convolutional Neural Networks (CNN) is trained offline to capture spatial features, reduce dimensional complexity, and enhance time efficiency, while a stacked Long Short-Term Memory (LSTM) is applied online to capture temporal information for improved prediction of future thermal behavior in subsequent printing layers. Model performance is evaluated using MSE, SSIM, and PSNR metrics and is benchmarked against stacked LSTM and convolutional LSTM models, demonstrating superior accuracy and applicability. Additionally, to mitigate noise from moving extruders and gantry backgrounds in thermal images, a fine-tuned semantic segmentation model is implemented offline to extract printing geometry, enabling precise temperature tracking along the tool path for further thermal analysis. The frameworks developed in this study significantly advance temperature monitoring, thermal analysis, and in-situ manufacturing control for LFAM, bridging the gap between theoretical modeling and practical application.

Geometry extraction↗

DOE OSTI2024

Active and Transfer Learning of High-Dimensional Neural Network Potentials for Transition Metals

Classical molecular dynamics (MD) simulations represent a very popular and powerful tool for materials modeling and design. The predictive power of MD hinges on the ability of the interatomic potential to capture the underlying physics and chemistry. There have been decades of seminal work on developing interatomic potentials, albeit with a focus predominantly on capturing the properties of bulk materials. Such physics-based models, while extensively deployed for predicting the dynamics and properties of nanoscale systems over the past two decades, tend to perform poorly in predicting nanoscale potential energy surfaces (PESs) when compared to high-fidelity first-principles calculations. These limitations stem from the lack of flexibility in such models, which rely on a predefined functional form. Machine learning (ML) models and approaches have emerged as a viable alternative to capture the diverse size-dependent cluster geometries, nanoscale dynamics, and the complex nanoscale PESs, without sacrificing the bulk properties. Here, in this study, we introduce an ML workflow that combines transfer and active learning strategies to develop high-dimensional neural networks (NNs) for capturing the cluster and bulk properties for several different transition metals with applications in catalysis, microelectronics, and energy storage, to name a few. Our NN first learns the bulk PES from the high-quality physics-based models in literature and subsequently augments this learning via retraining with a higher-fidelity first-principles training data set to concurrently capture both the nanoscale and bulk PES. Our workflow departs from status-quo in its ability to learn from a sparsely sampled data set that nonetheless covers a diverse range of cluster configurations from near-equilibrium to highly nonequilibrium as well as learning strategies that iteratively improve the fingerprinting depending on model fidelity. All the developed models are rigorously tested against an extensive first-principles data set of energies and forces of cluster configurations as well as several properties of bulk configurations for 10 different transition metals. Our approach is material agnostic and provides a methodology to transfer and build upon the learnings from decades of seminal work in molecular simulations on to a new generation of ML-trained potentials to accelerate materials discovery and design.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Identifying Key Drivers of Wildfires in the Contiguous US Using Machine Learning and Game Theory Interpretation

Abstract Understanding the complex interrelationships between wildfire and its environmental and anthropogenic controls is crucial for wildfire modeling and management. Although machine learning (ML) models have yielded significant improvements in wildfire predictions, their limited interpretability has been an obstacle for their use in advancing understanding of wildfires. This study builds an ML model incorporating predictors of local meteorology, land‐surface characteristics, and socioeconomic variables to predict monthly burned area at grid cells of 0.25° × 0.25° resolution over the contiguous United States. Besides these predictors, we construct and include predictors representing the large‐scale circulation patterns conducive to wildfires, which largely improves the temporal correlations in several regions by 14%–44%. The Shapley additive explanation is introduced to quantify the contributions of the predictors to burned area. Results show a key role of longitude and latitude in delineating fire regimes with different temporal patterns of burned area. The model captures the physical relationship between burned area and vapor pressure deficit, relative humidity (RH), and energy release component (ERC), in agreement with the prior findings. Aggregating the contribution of predictor variables of all the grids by region, analyses show that ERC is the major contributor accounting for 14%–27% to large burned areas in the western US. In contrast, there is no leading factor contributing to large burned areas in the eastern US, although large‐scale circulation patterns featuring less active upper‐level ridge‐trough and low RH two months earlier in winter contribute relatively more to large burned areas in spring in the southeastern US.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Fast All-sky Radiation Model for Solar applications (FARMS) [SWR-16-18]

The Fast All-sky Radiation Model for Solar applications (FARMS) is used to compute cloudy irradiance. Radiative transfer (RT) models simulating broadband solar radiation have been widely used by atmospheric scientists to model solar resources for various energy applications such as operational forecasting. Due to the complexity of solving the RT equation, the computation under cloudy conditions can be extremely time consuming though many approximations (e.g. two-stream approach and delta-M truncation scheme) have been utilized. Thus, a more efficient RT model is crucial for model developers as a new option for approximating solar radiation at the land surface with minimal loss of accuracy. We have developed a fast all-sky radiation model for solar applications (FARMS) using the simplified clear-sky RT model, REST2, and simulated cloud transmittances and reflectances from the Rapid Radiation Transfer Model (RRTM) with a sixteen-stream Discrete Ordinates Radiative Transfer (DISORT). Simulated lookup tables (LUTs) of cloud transmittances and reflectances were created by varying cloud optical thicknesses, cloud particle sizes, and solar zenith angles. Equations with optimized parameters were fitted to the cloud transmittances and reflectances to develop the model. Using this model the all-sky solar irradiance at the land surface can be computed rapidly by combining REST2 with the cloud transmittances and reflectances. This new RT model is more than 1000 times faster than those currently utilized in solar resource assessment and forecasting since it does not explicitly solve the RT equation for each individual cloud condition. Our results indicate the accuracy of the fast radiative transfer model is comparable to or better than two-stream approximation in term of computing cloud transmittance and solar radiation.

Xie, Yu↗

DOE OSTI2021

Binding of Campylobacter jejuni FliW Adjacent to the CsrA RNA-Binding Pockets Modulates CsrA Regulatory Activity

Campylobacter jejuni CsrA is an mRNA-binding, post-transcriptional regulator that controls many metabolic- and virulence-related characteristics of this important pathogen. In contrast to E. coli CsrA, whose activity is modulated by binding to small non-coding RNAs (sRNAs), C. jejuni CsrA activity is controlled by binding to the CsrA antagonist FliW. In this study, we identified the FliW binding site on CsrA. Deletion of the C-terminus of C. jejuni CsrA, which is extended relative to sRNA-binding CsrA proteins, abrogated FliW binding. Bacterial two-hybrid experiments were used to assess the interaction of FliW with wild-type CsrA and mutants thereof, in which every amino acid was individually mutated. Two CsrA mutations (V51A and N55A) resulted in a significant decrease in FliW binding. The V51A and N55A mutants also showed a decrease in CsrA-FliW complex formation, as assessed by size-exclusion chromatography and surface plasmon resonance. These residues were highly conserved in bacterial species containing CsrA orthologs whose activities are predicted to be regulated by FliW. The location of FliW binding was immediately adjacent to the two RNA-binding sites of the CsrA homodimer, suggesting the model that FliW binding to CsrA modulates its ability to bind to its mRNA targets either by steric hindrance, electrostatic repulsion, or by altering the overall structure of the RNA-binding sites.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2023

Fast All-sky Radiation Model for Solar applications (FARMS) 2023 [SWR-16-18]

Xie, Yu↗

DOE OSTI2022

Middle to Late Holocene Sea Surface Temperature and Productivity Changes in the Northeast Pacific

Variations of the sea surface temperature (SST) and primary productivity in the northeast Pacific have far‐reaching implications. In addition to influencing the regional and global temperature and hydroclimate, these conditions also control marine ecosystems and their services, which subsequently impact regional economies. Yet, our understanding of the variability and controls of northeast Pacific SST and productivity on timescales exceeding observational records remains limited. Here, we use marine sediment records from seven locations, spanning 25.2°N–59.6°N, in the northeast Pacific to characterize the millennial‐scale variability of SST and productivity from 9,000 to 1,000 years BP. We explore the dynamics of their spatiotemporal evolution and compare these data with transient climate model outputs to identify potential drivers. Through a heat budget analysis and optimal fingerprinting analysis, we characterize the spatial pattern of forcings. We find that SST varied spatially in the northeast Pacific, with higher latitudes exhibiting greater magnitude changes than lower latitudes, which differs from previous work suggesting regional synchronicity and coherence during the Holocene. Our analysis did not find evidence for coherent variability of primary producer community nor carbon export, highlighting the difficulty of identifying the complex interactions between environmental conditions, producers, and carbon export. Model‐proxy disagreement demonstrates the need for higher resolution model frameworks, but shows nonetheless that observed variability in the proxy records can be explained by a combination of greenhouse gas and orbital forcing. Here we suggest that the complex SST variations and marine ecosystem responses to forced changes are important factors that can drive disagreements in model projections.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Enabling the broader adoption of fusion simulation on complex geometry

This project addressed a key barrier to advanced fusion and nuclear simulation: the difficulty of performing high-fidelity Monte Carlo neutronics directly on complex, real-world CAD geometry. Traditional workflows require engineers to rebuild CAD models as simplified constructive solid geometry, a time-consuming and error-prone process that limits design iteration and broader adoption of simulation tools. The goal of this Phase I SBIR was to make CAD-based neutronics practical, accessible, and robust for industrial and research users. During the project, Coreform significantly enhanced the Direct Accelerated Geometry Monte Carlo (DAGMC) workflow and fully integrated it into Coreform Cubit as a first-class capability. Major achievements include optimized material assignment and surface meshing workflows, substantial performance improvements to geometry imprinting and preparation, native export of DAGMC models, and new visualization tools to support OpenMC source definition and lost-particle debugging. Coreform also expanded Cubit’s capabilities as a full OpenMC preprocessor, including the ability to convert OpenMC constructive solid geometry models back into CAD for visualization, multiphysics coupling, and debugging. In collaboration with Argonne National Laboratory, the project delivered comprehensive new DAGMC documentation and training materials, transforming DAGMC from a research-oriented tool into a production-ready workflow. Results were disseminated through tutorials, conference training, and multiple well-attended webinars demonstrating integrated CAD-based neutronics and multiphysics workflows. Overall, this project demonstrated that high-fidelity Monte Carlo simulations can be performed directly on complex CAD geometry, reducing setup time, improving usability, and enabling faster, more informed design decisions for fusion and nuclear energy systems.

42 ENGINEERING↗

DOE OSTI2024

Engineering Active Sites of Metal/Metal Oxide Catalysts by Oxide Ligand Overlayers

Abstract Engineering sites of supported metal catalysts is essential to enhancing activity and selectivity. Such enhancement is typically achieved by particle size modification, surface alloying, or attaching molecular ligands. Yet, control strategies for complex, multifunctional molecules and catalysts, where selectivity is crucial, are lacking. Here, we demonstrate that submonolayer WO x with tunable coverage preferentially decorates well‐coordinated Pt terrace sites as a stable ligand. By combining experimental kinetics with probe molecules, in situ spectroscopies, and first‐principles modeling, we show that the WO x coverage on Pt modifies the metal‐to‐acid site balance while retaining the acid strength intact and results in optimal reactivity for metal‐acid catalyzed reactions at a specific metal, size, and support‐dependent WO x coverage. The oxide can also alter the reactant adsorption mode, reversing selectivity and pathways from terrace‐ to step‐dominated, as evidenced in furfural decarbonylation and hydrogenation. The insights open avenues for improving metal/metal oxide catalysts beyond the specific system.

Zhou, Jiahua↗

DOE OSTI2020

The influence of wintertime SST variability in the Western North Pacific on ENSO diversity

Abstract Sea surface temperature anomalies (SSTa) in the Western North Pacific (WNP) have been linked to the development of El Niño Southern Oscillation (ENSO) events a full year in advance. However, the contribution of the WNP precursor to the temporal evolution and spatial complexity of ENSO remains unclear. Using the preindustrial experiment of the Community Earth System Model as the control climate, a partially coupled experiment is conducted in which WNP SSTa are restored to the model climatology. By comparing the perturbed experiment to the control, we are able to clearly characterize ENSO’s response to WNP SST variability. We find that the WNP is predominantly linked to eastern Pacific (EP) ENSO events. Without SST variability in the WNP, central Pacific (CP) ENSO events are more likely to develop. This variation in ENSO flavor is controlled by how the WNP projects onto wind stress anomalies in the western equatorial Pacific, which in turn impacts the discharge and recharge of ocean heat during the ENSO cycle. Specifically, the removal of SSTa in the WNP weakens the buildup of ocean heat in the western equatorial Pacific, which then hinders the development of EP-type events.

Meteorology & Atmospheric Sciences↗

DOE OSTI2021

Comparison of Observations and Predictions of Daytime Planetary-Boundary-Layer Heights and Surface Meteorological Variables in the Columbia River Gorge and Basin During the Second Wind Forecast Improvement Project

The second Wind Forecast Improvement Project (WFIP2) is an 18-month field campaign in the Pacific Northwest U.S.A., whose goal is to improve the accuracy of numerical-weather-prediction forecasts in complex terrain. The WFIP2 campaign involved the deployment of a large suite of in situ and remote sensing instrumentation, including eight 915-MHz wind-profiling radars, and surface meteorological stations. The evolution and annual variability of the daytime convective planetary-boundary-layer (PBL) height is investigated using the wind-profiling radars. Three models with different horizontal grid spacing are evaluated: the Rapid Refresh, the High-Resolution Rapid Refresh, and its nested version. The results are used to assess errors in the prediction of PBL height within the experimental and control versions of the models, with the experimental versions including changes and additions to the model parametrizations developed during the field campaign, and the control version using the parametrizations present in the National Oceanic and Atmospheric Administration/National Centers for Environmental Prediction operational version of the models at the start of the project. Results show that the high-resolution models outperform the low-resolution versions, the experimental versions perform better compared with the control versions, model PBL height estimations are more accurate on cloud-free days, and model estimates of the PBL height growth rate are more accurate than model estimates of the rate of decay. Finally, using surface sensors, we assess surface meteorological variables, finding improved surface irradiance and, to a lesser extent, improved 2-m temperature in the experimental version of the model.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Bayesian learning of orthogonal embeddings for multi-fidelity Gaussian Processes

Uncertainty propagation in complex engineering systems often poses significant computational challenges related to modeling and quantifying probability distributions of model outputs, as those emerge as the result of various sources of uncertainty that are inherent in the system under investigation. Gaussian Processes regression (GPs) is a robust meta-modeling technique that allows for fast model prediction and exploration of response surfaces. Multi-fidelity variations of GPs further leverage information from cheap and low fidelity model simulations in order to improve their predictive performance on the high fidelity model. In order to cope with the high volume of data required to train GPs in high dimensional design spaces, a common practice is to introduce latent design variables that are typically projections of the original input space to a lower dimensional subspace, and therefore substitute the problem of learning the initial high dimensional mapping, with that of training a GP on a low dimensional space. Here in this paper, we present a Bayesian approach to identify optimal transformations that map the input points to low dimensional latent variables. The \projection" mapping consists of an orthonormal matrix that is considered a priori unknown and needs to be inferred jointly with the GP parameters, conditioned on the available training data. The proposed Bayesian inference scheme relies on a two-step iterative algorithm that samples from the marginal posteriors of the GP parameters and the projection matrix respectively, both using Markov Chain Monte Carlo (MCMC) sampling. In order to take into account the orthogonality constraints imposed on the orthonormal projection matrix, a Geodesic Monte Carlo sampling algorithm is employed, that is suitable for exploiting probability measures on manifolds. We extend the proposed framework to multi-fidelity models using GPs including the scenarios of training multiple outputs together. We validate our framework on three synthetic problems with a known lower-dimensional subspace. The benefits of our proposed framework, are illustrated on the computationally challenging aerodynamic optimization of a last-stage blade for an industrial gas turbine, where we study the effect of an 85-dimensional shape parameterization of a three-dimensional airfoil on two output quantities of interest, specifically on the aerodynamic efficiency and the degree of reaction

42 ENGINEERING↗