Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ensemble simulations”

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

A generalized approximate control variate framework for multifidelity uncertainty quantification

We describe and analyze a variance reduction approach for Monte Carlo (MC) sampling that accelerates the estimation of statistics of computationally expensive simulation models using an ensemble of models with lower cost. These lower cost models — which are typically lower fidelity with unknown statistics — are used to reduce the variance in statistical estimators relative to a MC estimator with equivalent cost. We derive the conditions under which our proposed approximate control variate framework recovers existing multifidelity variance reduction schemes as special cases. We demonstrate that existing recursive/nested strategies are suboptimal because they use the additional low-fidelity models only to efficiently estimate the unknown mean of the first low-fidelity model. As a result, they cannot achieve variance reduction beyond that of a control variate estimator that uses a single low-fidelity model with known mean. However, there often exists about an order-of-magnitude gap between the maximum achievable variance reduction using all low-fidelity models and that achieved by a single low-fidelity model with known mean. We show that our proposed approach can exploit this gap to achieve greater variance reduction by using non-recursive sampling schemes. The proposed strategy reduces the total cost of accurately estimating statistics, especially in cases where only low-fidelity simulation models are accessible for additional evaluations. Several analytic examples and an example with a hyperbolic PDE describing elastic wave propagation in heterogeneous media are used to illustrate the main features of the methodology.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Reducing uncertainty of polar to midlatitude linkages using DOE’s E3SM in a coordinated model-experiment setting

This project brought DOE’s climate modeling effort with the Energy Exascale Earth System Model into the Polar Amplification Model Intercomparison Project (PAMIP), which is part of the sixth and latest Coupled Model Intercomparison Project, CMIP6. PAMIP examines the causes and consequences of polar amplification, when external forcing results in a larger temperature increase in high latitudes than the global average, in a coordinated set of model experiments with a common modeling protocol. Our teams from UC Irvine and the University of Toronto have designed, carried out, analyzed, and disseminated PAMIP output from the Energy Exascale Earth System Model (E3SM) and the Community Earth System Model (CESM). PAMIP’s ongoing stream of significant new results have advanced progress in the community’s understanding and led to new outstanding research questions that have motivated further work. PAMIP has led to improved consensus on the atmospheric response to sea ice loss. The important finding is that for a similar sea ice anomaly forcing, the simulated atmospheric response in the troposphere is remarkably consistent among the 16 models’ runs analyzed. The zonal-mean tropospheric response consists of a very robust equatorward shift of the westerly flow in mid-latitudes. However, while the multi-model mean response is robust, it has a weak amplitude relative to internal variability. We identified a weakness in the models (including E3SM) in terms of their response to sea-ice forcing that is related to eddy forcing (or nonlinear dynamical effects) at mid- to high latitudes. In fact, E3SM is an outlier in terms of the models participating, and in that sense, it turned out to be a vital participant model. We found that reductions in energy transport due sea-ice loss and involving dry air only are compensated by increases in moist energy transport from warmer sea surface temperature in midlatitudes. This leads to a large spread in energy transport into the Arctic and is a potential source of spread in Arctic amplification. We identified an important role that climate modes, including tropical modes of variability (El Nino and the Southern Oscillation (ENSO); the Quasi-biennial Oscillation (QBO)) play in the response to sea-ice anomalies, including in ocean coupled experiments. Similarly, we identified and quantified the contribution of sea-ice thickness to the atmospheric response compared to the response to sea-ice extent only. We found that it is important to run large ensembles and even with an ensemble size of 100 simulations the response is largely influenced by internal variability. We demonstrated convincingly that Ural blocking, not sea-ice loss, provides the weakening of the stratospheric polar vortex in fall/early winter and a negative phase of the North Atlantic Oscillation that can last for up to two months. However, sea-ice anomalies can influence the background flow so that the response to Ural blocking is more persistent under low sea-ice conditions in the Barents/Kara Sea than high sea ice Atmospheric model hierarchies that progressively add individual processes have a long history in providing dynamical insight for modeling the atmosphere. Similarly, coupled model hierarchies that progressively add individual ocean processes can provide insights into the workings of the coupled climate system, however such hierarchies have not been available except for a non-dynamical slab ocean model. Because of the missing processes, surface flux corrections must be added to produce a target climate. In this project, we managed to overcome this problem and develop a globally coupled ocean model hierarchy in CESM that can turn on and off the processes of mixed-layer entrainment and Ekman flow. We used the hierarchy to study the impact of Arctic sea-ice loss on the climate system. We find that the effect of mixed-layer entrainment on ocean heat uptake influences the atmospheric circulation by shifting the latitudinal positions of the mid-latitude westerly jet and the Intertropical Convergence Zone (ITCZ). In quadrupled CO 2 experiments, we studied how air-sea coupling affects the response of tropical rainfall under global warming. In order to identify the importance of individual ocean processes, we used the hierarchy of ocean models to separate the effects of seasonal mixed-layer entrainment, wind-driven Ekman flows, and frictional flows. We showed that including Ekman and frictional flows allows our simulation to produce the Pacific Ocean's enhanced equatorial warming pattern and equatorward ITCZ contraction noted in previous climate simulations. We also showed that the frictional flow, which has yet to receive much attention, is as important as the Ekman flow in generating equatorial heat convergence.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Regional climate change: consensus, discrepancies, and ways forward

Climate change has emerged across many regions. Some observed regional climate changes, such as amplified Arctic warming and land-sea warming contrasts have been predicted by climate models. However, many other observed regional changes, such as changes in tropical sea surface temperature and monsoon rainfall are not well simulated by climate model ensembles even when taking into account natural internal variability and structural uncertainties in the response of models to anthropogenic radiative forcing. This suggests climate model predictions may not fully reflect what our future will look like. The discrepancies between models and observations are not well understood due to several real and apparent puzzles and limitations such as the “signal-to-noise paradox” and real-world record-shattering extremes falling outside of the possible range predicted by models. Addressing these discrepancies, puzzles and limitations is essential, because understanding and reliably predicting regional climate change is necessary in order to communicate effectively about the underlying drivers of change, provide reliable information to stakeholders, enable societies to adapt, and increase resilience and reduce vulnerability. The challenges of achieving this are greater in the Global South, especially because of the lack of observational data over long time periods and a lack of scientific focus on Global South climate change. To address discrepancies between observations and models, it is important to prioritize resources for understanding regional climate predictions and analyzing where and why models and observations disagree via testing hypotheses of drivers of biases using observations and models. Gaps in understanding can be discovered and filled by exploiting new tools, such as artificial intelligence/machine learning, high-resolution models, new modeling experiments in the model hierarchy, better quantification of forcing, and new observations. Conscious efforts are needed toward creating opportunities that allow regional experts, particularly those from the Global South, to take the lead in regional climate research. This includes co-learning in technical aspects of analyzing simulations and in the physics and dynamics of regional climate change. Finally, improved methods of regional climate communication are needed, which account for the underlying uncertainties, in order to provide reliable and actionable information to stakeholders and the media.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Sensitivity Analysis of Drivers Water Shortage in the Los Angeles Region During Drought

The code and detailed step-by-step instructions for generating the model output data, processing results, and analysis and plotting are provided at https://github.com/IMMM-SFA/Ferencz_et_al_2026_ER_Water. The PyArtes model is a python adaptation of the Artes model. PyArtes uses many of the same input data and optimization model architecture as Artes. Documentation for the PyArtes model is provided in the Supplement to the paper. The primary data product are simulated monthly water shortages for indoor and outdoor demand under a large ensemble of drought scenarios (>13,000). The droughts are hypothetical and are not based on historical time series data of supply sources - though historical data did help inform ranges explored for supply parameters. Demands are informed by recent 2017-2021 water supply data. Demands used for the model can be accessed at https://github.com/IMMM-SFA/Ferencz_et_al_2026_ER_Water. Simulations resolve demand for over 90 water providers in the study region. The results report 36 months of water shortage data for each indoor and outdoor demand node. The study also developed a multilayer perceptron (MLP) neural network trained on a subset of the simulated shortage ensemble to emulate worst annual water shortage for a given set of parameter multipliers -- provided the parameter values fall within the ranges sampled in the ensemble. Emulated water shortages for synthetic ensembles are in the MLP-generated shortages folder. The MLP model was used to generate larger ensembles to support Sobol analysis that would have been extremely computationally expensive to simulate. Datasets provided in this repository*: Simulated shortages. These results are used for the analysis for Figures 5, 8, and 9 in the paper, and also to train the MLP emulator. .zip file containing outputs for the 13,312 scenario ensemble. Separate .csv files for indoor and outdoor shortage for each scenario. Rows = demand ids (~100), Columns = months (36) Units = acre-feet/month of shortage (shortage = monthly demand - supply). 1 acft = 1233.48 m^3 .csv files of aggregated shortages derived from the 13,312 ensemble Rows = scenarios (13,312), Columns = demand ids (~100) Units = acre-feet/year (either worst annual shortage or total shortage over the 3-year drought) .csv file of the parameter multipliers scenarios for the ensemble .csv file of the parameter ranges and baseline values the multipliers were applied to MLP-generated shortages. These results are used for Figures 4, 6, and 7 in the paper. mwd higher folder: scenario ensembles, emulated worst year total shortages (acft), and Sobol results Emulated shortages. Rows = scenarios, columns = demand ids, units acft Sobol results. Rows = demand ids, columns Sobol (S1, ST, or 95% confidence interval) value for each parameter mwd lower folder: scenario ensembles, emulated worst year total shortages (acft), and Sobol results same organization as mwd higher MLP performance: performance metrics (R^2, RMSE, BIAS, MAPE) for the testing subset (20% or 2,662 scenarios) and simulated vs emulated worst year shortage (acre-feet/year) for every demand node, MWD wholesale regions, and the entire study region (LAC). Supporting data for figures. Figure plotting scripts in the associated GitHub repo. These files support analysis and visualization. Geospatial Data used for plotting simulated water shortages and Sobol results. Dictionary of full names for demand nodes in the model and estimates of water supply by source type informed by Artes input files and California Urban Water Management Planning data: https://water.ca.gov/Programs/Water-Use-And-Efficiency/Urban-Water-Use-Efficiency/Urban-Water-Management-Plans *Readme files provided for each folder.

drought↗
DOE OSTI2021

ATHENA: A unique radiation environment platform at the National Ignition Facility

This paper describes the ATHENA platform, an energy tuning assembly, which was developed to spectrally shape the National Ignition Facility (NIF) deuterium–tritium fusion neutron source to a thermonuclear (fusion) plus prompt fission neutron spectrum with a capability to act as a short-pulse neutron source. This unique, otherwise inaccessible radiation environment complements existing experimental facilities and capabilities. Here, the flexible ATHENA irradiation positions were modeled using an ensemble of Monte Carlo simulations with stochastic sampling of the nuclear cross-sections to characterize the radiation environments and uncertainty for the platform. Validation of the internal neutron spectrum produced from fielding ATHENA at NIF occurred through neutron flux unfolding with 20 measured activation products.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
DOE OSTI2022

Future changes in extreme precipitation over the San Francisco Bay Area: Dependence on atmospheric river and extratropical cyclone events

Extreme precipitation poses a major challenge for local governments, including the City and County of San Francisco, California, as flooding can damage and destroy infrastructure and property. As the climate continues to warm, reliable future precipitation projections are needed to provide the best possible information to decision makers. However, future changes in the magnitude of extreme precipitation are uncertain, as current state-of-the-art global climate models are typically run at relatively coarse horizontal resolutions that require the use of convective parameterization and have difficulty simulating observed extreme rainfall rates. Here, we performed ensembles of convection-permitting regional climate model simulations to investigate how five historically impactful extreme precipitation events over the San Francisco Bay Area could change if similar events occurred in future climates. We found that changes in storm-total precipitation depend strongly on storm type. Precipitation associated with an atmospheric river (AR) accompanied by an extratropical cyclone (ETC) is projected to increase at a rate exceeding (by up to 1.5 times) the theoretical Clausius Clapeyron scaling of 6–7% per °C warming. On the other hand, future precipitation changes are weak or negative for events characterized by an AR only, despite increases in precipitable water and integrated vapor transport that are similar to those of the co-occurring AR and ETC events. The differences in the sign of future precipitation change between AR-only events and co-occurring AR and ETC events is instead linked with changes in mid-tropospheric vertical velocity. Given that the majority of observed ARs are associated with an ETC, this research has important implications for future precipitation impacts over the Bay Area, as it indicates that storm-total precipitation associated with the most common type of storm event may increase by up to 26–37% in 2100 relative to historical.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Structure–Function Relationships in Sequence-Controlled Copolymers for Rare Earth Element Chelation

The ability to tune material function through primary sequence is a defining feature of biological macromolecules, allowing precise control over structure and target interactions in complex aqueous environments. However, translating sequence–structure–function relationships to synthetic macromolecules is challenging due to their dispersity in sequence, conformation, and composition. Here, we report systematic studies of amphiphilic polymer chelators designed to probe how composition and patterning influence binding affinity and selectivity for rare earth elements (REEs), a series of technologically relevant metals with challenging separation profiles. A library of copolymers varying hydrophobic monomer composition and patterning was synthesized via reversible addition–fragmentation chain transfer (RAFT) polymerization, spanning statistical, gradient, and block architectures. REE binding was quantified using a high-throughput colorimetric assay, and reconstruction of polymer ensembles using kinetic stochastic simulations enabled quantitative comparisons of sequence heterogeneity, linking local monomer colocalization to emergent REE binding. Further, we investigated the role of different hydrophobic comonomers in tuning metal coordination, with binding trends linked to structural features that influence binding site desolvation. Complementary dynamic light scattering (DLS) and small-angle X-ray scattering (SAXS) measurements showed that both polymer and monomer architecture modulate metal-induced conformational changes, and that multichain assembly behavior emerges beyond critical hydrophobic thresholds. Sequence control also altered REE selectivity, with nonmonotonic differences observed across compositionally identical polymers with different sequence architectures. Together, these findings establish design principles that connect polymer sequence and structure to binding performance, guiding the design of macromolecular chelators with enhanced affinity and selectivity for applications in separations, sensing, and catalysis.

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

Future Changes in Midwest Extreme Precipitation Depend on Storm Type

Midwestern U.S. extreme precipitation is associated with multiple storm types including mesoscale convective systems (MCSs) and/or training thunderstorms, tropical cyclone (TC) remnants, and winter storms. Anthropogenic warming is expected to increase climatological precipitation globally, however, there may be little correspondence with regional storm-based changes. Furthermore, uncertainty remains in precipitation-temperature scaling due to use of convective parameterization in most global models. In this study, we investigated historically impactful extreme precipitation events from multiple types of Midwest storms using the Weather Research and Forecasting model at convection-permitting resolution. We simulated five-member ensembles of historical hindcasts and experiments representing the storms in the future using the pseudo-global warming method. We found that future precipitation changes depend on storm type, with increases near Clausius-Clapeyron (CC) for winter storms, no consensus for MCSs and/or training thunderstorms, and sub-CC increases for TC remnants. This research highlights the importance of considering storm type in future extreme precipitation projections.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Summer Convective Precipitation Changes Over the Great Lakes Region Under a Warming Scenario

To understand future summer precipitation changes over the Great Lakes Region (GLR), we performed an ensemble of regional climate simulations through the Pseudo-Global Warming (PGW) approach. We found that different types of convective precipitation respond differently to the PGW signal. Isolated deep convection (IDC), usually concentrated in the southern domain, shows an increase in precipitation to the north of the GLR. Mesoscale convective systems (MCSs), usually concentrated upwind of the GLR, shift to the downwind region with increased precipitation. Thermodynamic variables such as convective available potential energy (CAPE) and convective inhibition energy (CIN) are found to increase across almost the entire studied domain, creating a potential environment more favorable for stronger convection systems and less favorable for weaker ones. Meanwhile, changes in the lifting condensation level (LCL) and level of free convection (LFC) show a strong correlation with variations in convective precipitation, highlighting the significance of these thermodynamic factors in controlling precipitation over the domain. Our results indicate that the decrease in LCL and LCF in areas with increased convective precipitation is mainly due to increased atmospheric moisture. In response to the prescribed warming perturbation, MCSs occur more frequently downwind, while localized IDCs exhibit more intense rain rates, longer durations, and larger rainfall area.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Connecting energetics to dynamics in particle growth by oriented attachment using real-time observations

The interplay between crystal and solvent structure, interparticle forces and ensemble particle response dynamics governs the process of crystallization by oriented attachment (OA), yet a quantitative understanding is lacking. Using as a model system ZnO, we combine in situ TEM observations of single particle and ensemble assembly dynamics with simulations of interparticle forces and responses to relate experimentally derived interparticle potentials to the underlying interactions. We show that OA is driven by forces and torques due to a combination of electrostatic ion-solvent correlations and dipolar interactions that act at separations well beyond 5 nm. Importantly, coalignment is achieved before particles reach separations at which strong attractions drive the final jump to contact. The observed barrier to attachment is negligible, while dissipative factors in the quasi-2D confinement of the TEM fluid cell lead to abnormal diffusivities with timescales for rotation much less than for translation, thus enabling OA to dominate.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Local-scale Arctic tundra heterogeneity affects regional-scale carbon dynamics

In northern Alaska nearly 65% of the terrestrial surface is composed of polygonal ground, where geomorphic tundra landforms disproportionately influence carbon and nutrient cycling over fine spatial scales. Process-based biogeochemical models used for local to Pan-Arctic projections of ecological responses to climate change typically operate at coarse-scales (1km 2 –0.5°) at which fine-scale (<1km 2 ) tundra heterogeneity is often aggregated to the dominant land cover unit. Here, we evaluate the importance of tundra heterogeneity for representing soil carbon dynamics at fine to coarse spatial scales. We leveraged the legacy of data collected near Utqiagvik, Alaska between 1973 and 2016 for model initiation, parameterization, and validation. Simulation uncertainty increased with a reduced representation of tundra heterogeneity and coarsening of spatial scale. Hierarchical cluster analysis of an ensemble of 21 st -century simulations reveals that a minimum of two tundra landforms (dry and wet) and a maximum of 4km 2 spatial scale is necessary for minimizing uncertainties (<10%) in regional to Pan-Arctic modeling applications.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Projected runoff declines from plant physiological effects on precipitation

Abstract The impact of plants on runoff under high atmospheric CO 2 is a major uncertainty for future water resources. Theory and Earth system models (ESMs) suggest that stricter plant stomatal regulation under high CO 2 will reduce transpiration, potentially boosting runoff. Yet, across a 12-member ensemble of idealized ESM simulations that isolate plant responses to CO 2 , we show that lower transpiration robustly enhances runoff over only 5% of modelled global land area. Precipitation changes are five times more important than transpiration changes in driving runoff responses and are a significant signal of CO 2 physiological forcing over 31–57% of land areas across models. Crucially, ESMs largely disagree on where physiologically forced precipitation changes occur but agree that plant responses in most locations are as likely to reduce runoff as increase it. These results imply that large model uncertainties in precipitation responses, rather than transpiration responses, explain why ESMs disagree on plant physiologically driven runoff changes.

Lesk, Corey S. (ORCID:0000000308684040)↗
DOE OSTI2022

Calculating adsorption isotherms using the two-phase thermodynamic method and molecular dynamics simulations

We describe the calculation of adsorption isotherms from molecular dynamics simulations based on the two-phase thermodynamic (2PT) model. The 2PT model developed for bulk fluid phases treats the gas-like components as hard spheres (HSs), which correctly recovers the limiting behaviors of unconfined fluids. We showed that this treatment, however, does not always lead to the correct zero-loading behavior in strongly confining systems. For methane adsorption into zeolite MFI, the HS reference state underestimates entropy by up to 20% at low loadings and leads to an order-of-magnitude increase in the adsorption onset pressure. To fix these issues, we propose the use of ideal adsorbed gas (IAG) as the gas reference model, the properties of which can be computed using the Widom insertion method on an empty adsorbent. We further describe three routes to compute adsorption isotherms from the Helmholtz free energy at different loadings. Comparing against established Monte Carlo (MC) methods, we found that the adsorption isotherms obtained using the IAG reference state agrees to within 40%, which corresponds to deviations of <5% in adsorption free energy. The isotherms calculated using the HS reference state underestimate the adsorption uptake at low to medium loadings in strongly confining systems, but its accuracy improves at higher loadings and as the pore size increases relative to the sorbate diameter. The methods described here provide an alternative approach for computing adsorption isotherms when MC simulations in an open ensemble are undesirable and enable a direct comparison of computed adsorption thermodynamics with experiments.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Melting point of iron at high pressure: An assessment of uncertainties and effect of electronic temperature

An accurate calculation of the melting point of iron at various pressures in the Earth's core is important for understanding the core structure, geodynamo, and the Earth's history. Previous studies have assessed the melt line of iron at these extreme conditions using various experimental measurement techniques as well as both ab initio and classic molecular dynamics simulations. However, experimental measurements have uncertainties up to several hundred Kelvin, and inconsistencies remain among simulation results. Here in this work, we propose an iterative framework that couples density functional theory (DFT) calculations and molecular dynamics simulations performed using an ensemble of interatomic potentials to assess the effect of electronic temperature on the melting point. We systematically validate the potentials by comparing lattice constants and phonon dispersion curves at 0 K and enthalpy differences between liquid and HCP, FCC, BCC phases of iron close to the melt line at 300 GPa with DFT. Our results show that HCP iron melts at 6144 K (at 300 GPa), BCC phase is thermodynamically unstable, and FCC is metastable at this temperature. The melting points of FCC and BCC phases at 300 GPa are 5858 and 5647 K, respectively.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Variable domain mutational analysis to probe the molecular mechanisms of high viscosity of an IgG 1 antibody

Subcutaneous injection is the preferred route of administration for many antibody therapeutics for reasons that include its speed and convenience. However, the small volume limit (typically ≤2 mL) for subcutaneous delivery often necessitates antibody formulations at high concentrations (commonly ≥100 mg/mL), which may lead to physicochemical problems. For example, antibodies with large hydrophobic or charged patches can be prone to self-interaction giving rise to high viscosity. Here, we combined X-ray crystallography with computational modeling to predict regions of an anti-glucagon receptor (GCGR) IgG 1 antibody prone to self-interaction. An extensive mutational analysis was undertaken of the complementarity-determining region residues residing in hydrophobic surface patches predicted by spatial aggregation propensity, in conjunction with residue-level solvent accessibility, averaged over conformational ensembles from molecular dynamics simulations. Dynamic light scattering (DLS) was used as a medium throughput screen for self-interaction of ~200 anti-GCGR IgG 1 variants. A negative correlation was found between the viscosity determined at high concentration (180 mg/mL) and the DLS interaction parameter measured at low concentration (2–10 mg/mL). Additionally, anti-GCGR variants were readily identified with reduced viscosity and antigen-binding affinity within a few fold of the parent antibody, with no identified impact on overall developability. The methods described here may be useful in the optimization of other antibodies to facilitate their therapeutic administration at high concentration.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Scalable Risk Assessment of Rare Events in Power Systems With Uncertain Wind Generation and Loads

Risk assessment of rare events has become increasingly important in power system planning and operation with the increasing integration of renewable energy and the presence of system uncertainties. However, quantifying the risk posed by rare events via the traditional method, i.e., Monte Carlo sampling (MCS), incurs substantial computational expense stemming from the vast ensemble of power flow simulations. To accelerate the assessment, this paper proposes a Deep Neural Network (DNN)-kernelized vector-valued Gaussian Process (VVGP) approach with excellent computational efficiency while maintaining high accuracy. Consequently, serving as a surrogate model for the power flow solver, the DNN-kernelized VVGP enables significantly faster but accurate risk assessment compared to the power flow solver. The developed surrogate model evaluates low-order N - k events that contain more than 90% instances by adeptly capturing the topological features while the high-order N - k events are assessed via a power flow solver, thereby striking a balance between computational efficiency and uncertainty quantification accuracy. Moreover, the model incorporates a Support Vector Machine (SVM) classifier to resample concerning low-probability tail events to counteract the biases potentially introduced during the DNN-kernelized VVGP evaluations. Simulations conducted on the modified IEEE 24-bus, 118-bus, and European 1354-bus systems demonstrate that the proposed method maintains the accuracy benchmark set by MCS while significantly reducing computational demands in large-scale power systems as compared to other state-of-the-art methods.

17 WIND ENERGY↗
DOE OSTI2021

AI-driven multiscale simulations illuminate mechanisms of SARS-CoV-2 spike dynamics

We develop a generalizable AI-driven workflow that leverages heterogeneous HPC resources to explore the time-dependent dynamics of molecular systems. We use this workflow to investigate the mechanisms of infectivity of the SARS-CoV-2 spike protein, the main viral infection machinery. Our workflow enables more efficient investigation of spike dynamics in a variety of complex environments, including within a complete SARS-CoV-2 viral envelope simulation, which contains 305 million atoms and shows strong scaling on ORNL Summit using NAMD. We present several novel scientific discoveries, including the elucidation of the spike’s full glycan shield, the role of spike glycans in modulating the infectivity of the virus, and the characterization of the flexible interactions between the spike and the human ACE2 receptor. We also demonstrate how AI can accelerate conformational sampling across different systems and pave the way for the future application of such methods to additional studies in SARS-CoV-2 and other molecular systems.

59 BASIC BIOLOGICAL SCIENCES↗