Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Modeling and Simulation”

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

Developing Partnership between San Jose State University and DOE Lawrence Livermore National Laboratory to Enhance Climate Research Equity and Inclusion

One of the key objectives of this project was to develop partnership between San Jose State University (SJSU) and the Lawrence Livermore National Laboratory (LLNL), a US Department of Energy (DOE) funded national laboratory. Both institutions are closely located within the San Francisco Bay Area in California and their researchers share overlapping research interests and expertise related to Earth system sciences. By leveraging their connections with LLNL, faculty and students from SJSU gained exposure to state-of-the-art observations and simulations related to Earth system sciences, including but not limited to the usage of the facility data provided by the DOE Atmospheric Radiation Measurement (ARM) program and data analysis techniques for interpreting and analyzing the DOE Energy Exascale Earth System Model (E3SM) simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Assessment of CTF for Steady-state and Transient Post-CHF Conditions in Support of Time-at-Temperature Modeling Applications

The US nuclear industry is exploring options to improve operational economics and uprate the current fleet of light-water reactors by investigating transitioning to cladding performance–based safety criteria as opposed to the current limit, which requires complete avoidance of critical heat flux (CHF)/dryout. Past experience has shown that not all events leading to a dryout are severe enough to cause fuel performance degradation. Allowing temporary dryout of the fuel—that is, using a time-at-temperature (TaT) strategy—could allow for economic improvements via large power uprates and enhanced operational flexibility for current plants without compromising fuel integrity. To support this effort, the US Department of Energy is executing a comprehensive program that includes generating cladding material data under TaT conditions, developing new mechanistic models, and demonstrating modeling and simulation capabilities for transients of interest. This paper presents work performed to assess the CTF thermal-hydraulics subchannel code. CTF is a package used in the VERA core simulator, which will ultimately be used for TaT analysis. CTF will provide the thermal-hydraulic boundary conditions that will be needed for fuel performance analysis in the BISON code. Quantifying both the accuracy and uncertainty of post-CHF models will therefore be necessary. This paper outlines the strategy for the assessment of TaT and presents the results of using the steady-state and transient dryout experiments of the Boiling Fine-mesh Bundle Tests for CTF validation. The results show that the current model tends to overpredict steady-state critical power. This behavior translates to the transient tests, in which CTF is unable to capture transient dryout behavior. Some discussion of sensitivity analysis work being performed is provided to indicate which models must be further analyzed to properly model transient dryout and its uncertainty.

Salko Jr, Robert [ORNL] (ORCID:0000000253566679)↗
DOE OSTI2026

Effect of Surrogate Fission Products on Thermal Conductivity and Viscosity of NaCl-UCl 3

Thermophysical properties of salt systems relevant to molten salt reactors (MSRs) are experimentally measured in support of the US Department of Energy Office of Nuclear Energy (DOE-NE) MSR campaign within the Advanced Reactor Technology (ART) Program. These property measurements also support the development of the thermophysical arm of the Molten Salt Database (MSD) within DOE-NE’s Nu clear Energy Advanced Modeling and Simulation (NEAMS) Program. Several US Department of Energy (DOE) laboratories, including Oak Ridge National Laboratory (ORNL), have been capitalizing on modern methodologies and sample characterization techniques to conduct thermophysical property measurements to support these programs and to provide MSR developers and modelers with more recent and higher-quality data.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2026

Simulation of electron Bernstein waves using FullWave with a 2D non-local hot plasma model

Hot plasma wave simulation capability is expanded in the FullWave code by updating the hybrid iterative solver in the code with a semi-implicit time stepping method. The new approach is used to simulate Electron Bernstein Wave (EBW) heating in over-dense spherical tokamak plasmas. The code’s hybrid iterative solver circumvents the prohibitive memory cost of direct methods by combining a time evolution of Maxwell’s equations with frequency-domain relaxation, while the conductivity kernel, calculated via 3D particle tracking, captures the essential non-local wave–particle interactions. One-dimensional EBW simulations verify the algorithm’s accuracy by demonstrating mode conversion from X-mode wave to EBW at the upper hybrid resonance and a strong cyclotron damping near the plasma core. Two-dimensional simulation reproduces the predicted short EBW wavelength and quantitatively matches the hot-plasma dispersion relation. This study demonstrates the fidelity of the hybrid solver for the electron cyclotron frequency range.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials: Part II Resveratrol Exemplar

This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

An Evaluation and Qualification of U.S.-Based Research Reactors for Irradiation Capabilities Supporting Advanced Nuclear Systems

Irradiation experiments are a prerequisite for evaluating nuclear reactor system designs, analyzing the performance of these systems, and obtaining licenses. Likewise, irradiation facilities are necessary for producing the radioisotopes used in industrial and medical applications. Recent developments in modeling and simulation capabilities and advancements in computational resources have further enabled the design of irradiation experiments for evaluating radiation-induced phenomena and determining nuclear fuel, material, and system design and safety criteria pertaining to both normal and accident scenarios. These computational tools and models require comprehensive experimental datasets acquired under prototypic radiation conditions—for exploring material and system performance under the uniquely harsh environments found in nuclear reactors—to enable verification and validation for qualification and licensing purposes. However, qualification of irradiation experimental facilities, primarily research and test reactors (RTRs), necessitates that their performance be evaluated based on the irradiation environment (e.g. flux, power, testing capabilities) using an appropriate scoring matrix. Although many university campus RTRs are available for research and development (R&D) activities and initiatives, this study focuses on evaluating and qualifying the irradiation facilities (mostly RTRs) within the United States that are suitable for advanced nuclear fuel, material, and system irradiation experiments aimed at establishing operational-performance limits and informing component and fuel designs so as to improve operational efficiencies and mitigate proliferation vulnerabilities, as well as for radioisotope production aimed at multipurpose applications. As a result, the findings of the present study support the acceleration of nuclear fuel and material qualifications, thus hastening new and advanced nuclear energy system demonstrations and radioisotope production efforts by using extended R&D.

irradiation experiment↗
DOE OSTI2025

EVSE Characterization: V2G EVSE Comparison

As part of the U.S. Department of Energy EVs@Scale consortium Next-Generation Profiles project, results and analysis from the characterization of high-power conductive and wireless charging infrastructure are presented. This characterization is conducted over a wide range of direct current (DC) current and DC voltage operation for nominal test conditions and off-nominal test conditions. Test plans and procedures were developed to define the test configurations and requirements, measurement parameters, and test procedures used throughout testing. Results from a 2024 study conducted on electric vehicle supply equipment (EVSE) characterization by the Idaho National Laboratory (INL) include two bi-directional vehicle-to-grid (V2G) capable EVSEs. These EVSE are referred to as V2G-EVSE9 and V2G-EVSE10. Laboratory testing is conducted at nominal test conditions to characterize the power transfer capabilities, efficiency, power factor, and other power quality metrics of the two DC EVSEs capable of V2G bi-directional power transfer. Results from testing show the performance is consistent for V2G-EVSE9 and V2G-EVSE10 when comparing charging to discharging performance, except for V2G-EVSE10 for power transfer when operating above 70% of the rated DC current. V2G-EVSE10 efficiency is >98% while charging and <91% while discharging at the same operating conditions, near maximum-rated current, at 300VDC. In contrast, V2G-EVSE9 results are consistent for charging and discharging. This EVSE is nearly 96% efficient while charging or discharging when operating over 50% of rated AC power. V2G EVSE performance is also characterized during off-nominal AC grid conditions involving AC voltage deviation (426 VAC to 518 VAC), AC frequency deviation of +2% (58.8 Hz to 61.2 Hz), and AC voltage harmonics injection. Many test conditions have little-to-no impact on performance characteristics of the two EVSEs; however, there are a few notable findings with significant power transfer capability impacts. AC voltage harmonics injection resulted in negative impacts on power quality attributes for both EVSEs, but with no impact on power transfer capability. Off-nominal AC voltage and frequency conditions resulted in unstable or lack of power transfer capability for both EVSEs. V2G-EVSE9 is unable to transfer power when AC voltage is >300V L-N. V2G-EVSE10 is unable to transfer power when AC frequency deviation exceeds +0.8%. V2G energy management system transient response and latency are quantified during laboratory testing. V2G-EVSE9 and V2G-EVSE10 utilize cloud-based V2G energy management systems that command the power transfer level between the EVSE and EV. The latency and response characteristics of the entire systems (web-based user interface, V2G energy management system, cellular communications, and EVSE response) are quantified through laboratory testing for V2G-EVSE9 and V2G-EVSE10. V2G-EVSE9 latency ranges from 0.8 to 1.8 seconds, whereas V2G-EVSE10 latency ranges from 3.4 to 8.8 seconds. The ramp rate to a change in power transfer request also differs between the two EVSEs. V2G-EVSE10 ramp rate ranges from 50% to -250% of rated AC power per second, whereas V2G-EVSE9 rate ranges from 95% to -95% of rated AC power per second. At the highest rate of change in power transfer, V2G-EVSE10 can change from full charge power to full discharge power in less than one second. The V2G EVSE characterization presented in this report provides valuable insights and results for use by numerous entities. This includes modeling and simulation organizations, decision makers, fleet planning, industry stakeholders, and many others involved with the development and deployment of electrified transportation technologies. Additional high-power DC chargers, bidirectional chargers, and inductive power transfer EVSE characterization results are anticipated from additional EVSE brands and models, which will be detailed in future publications in support of the U.S. Department of Energy EVs@Scale consortium Next-Gen Profiles project.

25 ENERGY STORAGE↗
DOE OSTI2024

An Active Learning-Based Streaming Pipeline for Reduced Data Training of Structure Finding Models in Neutron Diffractometry

Structure determination workloads in neutron diffractometry are computationally expensive and routinely require several hours to many days to determine the structure of a material from its neutron diffraction patterns. The potential for machine learning models trained on simulated neutron scattering patterns to significantly speed up these tasks have been reported recently. However, the amount of simulated data needed to train these models grows exponentially with the number of structural parameters to be predicted and poses a significant computational challenge. To overcome this challenge, we introduce a novel batch-mode active learning (AL) policy that uses uncertainty sampling to simulate training data drawn from a probability distribution that prefers labelled examples about which the model is least certain. We confirm its efficacy in training the same models with ∼ 75% less training data while improving the accuracy. We then discuss the design of an efficient stream-based training workflow that uses this AL policy and present a performance study on two heterogeneous platforms to demonstrate that, compared with a conventional training workflow, the streaming workflow delivers ∼ 20% shorter training time without any loss of accuracy.

Wang, Tianle [Brookhaven National Laboratory (BNL)↗
DOE OSTI2025

High-Fidelity Modeling of a Type-5 Wind Turbine Gearbox

Type-5 wind turbines are characterized by their use of a hydraulic torque converter and permanent synchronous generator. This combination promotes steady and grid-ready energy without the use of a power converter. Thus, researchers were prompted to study the potential impact on grid reliability, stability, and resilience using a Real Time Digital Simulator (RTDS) model of the type-5 turbine, including a high-fidelity model of its gear box.

13 - HYDRO ENERGY↗
DOE OSTI2025

TRUST Sensors in Environments: Accelerometers (SE-A) Report Release FY25

The objective of the Delivery Environments (DE) Testbeds to Reduce Uncertainty in Simulations and Tests (TRUST) project is to quantify and help increase confidence in specific areas of computa tional and experimental capabilities that are applicable to the development, on-target assessment, and qualification of current and future delivery environments. The TRUST project consists of five "low-complexity", single-feature testbeds used to conduct experiments with accompanying models and simulations. Each of the five single feature testbeds aims to isolate an engineering feature or behavior of interest, then work towards improving the fundamental engineering understanding of that feature. Additionally, the single feature testbeds are used to identify capability development needs that can help reduce model and experimental uncertainty.

42 ENGINEERING↗
DOE OSTI2026

Compound flood simulations indicate rising public exposure to sewage-contaminated waters in Waikīkī, Hawai‘i

Waikīkī, Hawai‘i serves as a valuable case study for flood vulnerability due to its low elevation, present-day exposure to multiple flood sources, and role as the state’s economic engine. This study applies the high-resolution, open-source WRF-Hydro-CUFA (Coastal Urban Flood Application) model to simulate compound flood interactions among pluvial, fluvial, coastal, storm drain-driven, and subsurface processes, and to evaluate how flood characteristics in Waikīkī evolve under projected sea-level rise (SLR). Simulated water levels in an estuarine drainage canal, serving as the primary outlet for gravity-flow drainage from Waikīkī, show strong agreement with observations (r ≥ 0.935, R 2 ≥ 0.874), accurately reproducing flooding observed during recent 5-year and 50-year storm events. Simulations incorporating SLR demonstrate a shift in dominant flood drivers as sea levels rise. In addition, results identify critical thresholds at which elevated estuarine tailwater conditions cause widespread drainage backflow, eventually occurring in the absence of rainfall. These findings suggest that tidally driven flooding will progressively dominate over rainfall-driven events, increasing public exposure to highly contaminated canal and groundwater sources. The results underscore the urgent need to integrate water quality issues and drainage failure into coastal flood management strategies for low-lying urban environments, particularly those reliant on polluted estuarine waterways for stormwater management.

Compound flooding↗
DOE OSTI2025

First-Principles Dissociation Pathways of BCl 3 on the Si(100)-2 × 1 Surface

BCl 3 is a promising acceptor precursor for atomic-precision δ-doping of silicon, as it has been observed to rapidly dissociate into boron doped into the silicon surface and surface chlorine, which can be removed upon annealing. The chemical pathway and the resulting kinetics, through which BCl 3 adsorbs and dissociates on silicon, however, have only been partially explained. Here, in this work, we use density functional theory to expand the dissociation reactions of BCl 3 to include reactions that take place across multiple silicon dimer rows and reactions which end in a bare B atom either at the surface, substituted for a surface silicon, or in a subsurface position. We further simulate the resulting scanning tunneling microscopy images for each of these BCl x dissociation fragments, demonstrating that they often display distinct features that may allow for relatively confident experimental identification. Finally, we input the full dissociation pathway for BCl 3 into a kinetic Monte Carlo model, which simulates realistic reaction pathways as a function of environmental conditions, such as the pressure and temperature of dosing. We find that BCl 2 is broadly dominant at low temperatures, while high temperatures and ample space on the silicon surface for dissociation encourage the formation of bridging BCl fragments and B substitutions on the surface. This work provides the chemical mechanisms for understanding atomic-precision doping of Si with B, enabling a number of relevant quantum applications, such as bipolar nanoelectronics, acceptor-based qubits, and superconducting Si.

Campbell, Quinn T. [Sandia National Laboratories (↗
DOE OSTI2025

End-of-Century Changes in Orographic Precipitation with the Intermediate Complexity Atmospheric Research Model over the Western United States

Abstract Downscaled precipitation projections were created using the Intermediate Complexity Atmospheric Research (ICAR) model over the western United States to increase the physical realism in orographic precipitation changes. End-of-century simulations from eight models in phase 5 of the Coupled Model Intercomparison Project (CMIP5) were downscaled with ICAR and compared to the widely utilized statistically downscaled dataset, localized constructed analogs (LOCAs), to understand where and why projections of cool-season (September–May) precipitation differed. ICAR and LOCA precipitation projections were similar, but their sign differed in hydrologically relevant regions likely due to ICAR’s simulation of microphysics and mesoscale dynamics with high-resolution topography (6 km). In the Pacific Northwest, cool-season precipitation projections from ICAR showed an increase on the windward side of the Cascades and no significant change within the lee. This difference between the windward and leeward side was attributed to reduced zonal wind speeds, allowing more time for microphysical processes within ICAR. This contrast is enhanced by rain’s faster fall speed compared to snow, limiting transport into the lee. Meanwhile, LOCA projected an increase in precipitation across the Cascades. In the Upper Colorado River basin, LOCA projected an increase in precipitation in high elevation regions (>3000 m), but ICAR projected no significant change or a decrease in precipitation. High elevation differences were most evident in the spring and fall and were also attributed to a snow-to-rain transition and dynamical processes that impacted orographic enhancement within ICAR. Idealized, controlled studies are needed to better isolate individual processes, but these results underscore the importance of including microphysics and mesoscale dynamics within regional-scale precipitation projections. Significance Statement A set of global climate model simulations was downscaled using an atmospheric model that contains key physical equations, referred to as Intermediate Complexity Atmospheric Research (ICAR). ICAR was used to examine projected changes in end-of-century cool-season precipitation over mountains in the western United States. Precipitation projections from ICAR were similar to projections that used statistical relationships to downscale climate projections. However, projections differed between ICAR and statistically downscaled datasets in whether they increased, decreased, or stayed the same in specific, hydrologically relevant regions such as the eastern Cascades and high elevation areas of the Upper Colorado River basin. These differences were attributed to the simulation of physical processes in ICAR. The results highlight the importance of kilometer-scale atmospheric processes in regional climate projections.

Currier, William Ryan [NOAA/Physical Sciences Labo↗
DOE OSTI2025

Comparative Evaluation of the Ability of the MYNN‐EDMF PBL Scheme in WRF Model to Reproduce Near Surface Wind Speed Over Different Topographical Types

Abstract This study systematically evaluates the performance of the Mellor‐Yamada‐Nakanishi‐Niino‐Eddy‐Diffusion‐Mass‐Flux planetary boundary layer (PBL) scheme within the Weather Research and Forecasting (WRF) model in simulating near‐surface wind speeds across various topographies in New York State (NYS). Simulated wind speeds are compared with in‐situ measurements from 22 surface sites, grouped into six topographic categories: continental plain (CT), lakeside (LS), river valley (RV), Long Island (LI), Block Island (BI), and offshore ocean (OO). A quantitative evaluation based on Relative Euclidean Distance shows that wind speeds at the OO site are the most accurately reproduced, followed by those at LI sites, while the model performs less accurately for the remaining topographic groups. Wind speeds over CT sites tend to be overestimated by approximately 1 m/s, although their diurnal variability (DV) is well captured. In contrast, the model underestimates wind DV at LS, RV, LI, and BI sites, with the largest biases occurring at LI and BI, resulting in underestimated daytime wind speed and/or overestimated nighttime wind speed. The OO winds exhibit minimal diurnal variation, accurately captured by our WRF model. The surface wind diurnal variation is closely linked to PBL development. Among the indicators of PBL development, surface potential temperature biases most strongly correlate with wind speed biases. Our WRF model faces challenges in capturing the distinctions between winds influenced by local circulations and those over continental plains, and the significantly stronger winds at OO compared to BI. Potential causes for these biases are discussed, offering pathways for improving surface wind simulations in future.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A Parametric, Data-Driven, Non-Intrusive Reduced-Order Model Framework for Crystal Plasticity Simulations of Voids

The influence of the internal structure at micrometer length scales on the deformation of polycrystalline materials can be effectively captured using crystal plasticity finite element methods (CPFEM). However, the complexity and nonlinearity of the deformation equations CPFEM solves demand significant computational power and resources to achieve accurate predictions, limiting its broader application. To address this challenge, we have identified a reduced-order representation of the complex data in order to establish a computationally efficient reduced-order models (ROM) and drastically reduce the computational expense of CPFEM. Specifically, in this work, we developed a parametric, data-driven, and non-intrusive ROM framework for CPFEM using proper orthogonal decomposition (POD) and sparse variational Gaussian process (SVGP) regression for single-crystal microstructures under tensile loading conditions. The developed protocol enables one to compress field into a latent/low-dimensional space described by principal component analysis (PCA) via the singular value decomposition (SVD) algorithm. As a result, the high-dimensional data are reduced to a significantly smaller amount of dimensions with POD bases and POD coefficients. Furthermore, we deployed an ensemble of SVGPs—extended from the classical Gaussian process (GP) regression for scalability and handling big data—in a massively parallel manner to train and predict latent POD coefficients using known POD bases from a set of previously obtained simulations results. Lastly, using the predicted POD coefficients, we reconstructed the full-field results and showed reasonable agreement compared with the true values obtained from running CPFEM. The developed framework is validated with a set of CPFEM simulations of a single embedded void in single-crystal aluminum alloy. While the framework is broadly applicable, this work specifically focuses on single-crystal microstructures, a single load case (e.g., tensile), and a specific void geometry (spherical).

Anisotropy↗
DOE OSTI2025

Measurement of three-dimensional inclusive muon-neutrino charged-current cross sections on argon with the MicroBooNE detector

We report the measurement of the triple-differential cross section d 3 σ/dE vis d cos(θ μ )dP μ for inclusive muon-neutrino charged-current scattering on argon. This measurement utilizes data from 6.4 x 10 20 protons on target of exposure collected using the MicroBooNE liquid argon time projection chamber located along the Fermilab Booster Neutrino Beam with a mean neutrino energy of approximately 0.8 GeV. The mapping from reconstructed kinematics to truth quantities is validated within uncertainties by comparing the distribution of reconstructed hadronic energy in data to that of the model prediction in different muon scattering angle bins after applying a conditional constraint from the muon momentum distribution in data. The success of this validation provides confidence that the energy transfer in the MicroBooNE detector is well-modeled within simulation uncertainties, enabling a reliable unfolding to a triple-differential cross section defined at the nominal neutrino flux over muon momentum, muon scattering angle, and visible neutrino energy. This validation not only supports accurate cross-section extraction, but also establishes a critical foundation for tuning interaction models used in future neutrino oscillation measurements. The unfolded measurement covers an extensive phase space, providing a wealth of information useful for future liquid argon time projection chamber experiments measuring neutrino oscillations. Comparisons against a number of commonly used model predictions are included and their performance in different parts of the available phase-space is discussed.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2025

An Optimization-Based Coupling of Reduced Order Models with an Efficient Reduced Adjoint Basis Generation Approach

Optimization-based coupling (OBC) is an attractive alternative to traditional Lagrange multiplier approaches in multiple modeling and simulation contexts. However, application of OBC to time-dependent problems has been hindered by the computational cost of finding the stationary points of the associated Lagrangian, which requires primal and adjoint solves. This issue can be mitigated by using OBC in conjunction with computationally efficient reduced order models (ROMs). To demonstrate the potential of this combination, in this paper, we develop an optimization-based ROM-ROM coupling for a transient advection-diffusion transmission problem. We pursue the “optimize-then-reduce” path toward solving the minimization problem at each time step and solve reduced space adjoint system of equations, where the main challenge in this formulation is the generation of adjoint snapshots and reduced bases for the adjoint systems required by the optimizer. One of the main contributions of the paper is a new technique for an efficient adjoint snapshot collection for gradient-based optimizers in the context of optimization-based ROM-ROM couplings. In conclusion, we present numerical studies demonstrating the accuracy of the approach along with comparison between various approaches for selecting a reduced order basis for the adjoint systems, including decay of snapshot energy, average iteration counts, and timings.

coupled problems↗
DOE OSTI2025

Modeling the impact of structure and coverage on the reactivity of realistic heterogeneous catalysts

Adsorbates often cover the surfaces of catalysts densely as they carry out reactions, dynamically altering their structure and reactivity. Understanding adsorbate-induced phenomena and harnessing them in our broader quest for improved catalysts is a substantial challenge that is only beginning to be addressed. Here, in this work, we chart a path toward a deeper understanding of such phenomena by focusing on emerging in silico modeling methodologies, which will increasingly incorporate machine learning techniques. We first examine how adsorption on catalyst surfaces can lead to local and even global structural changes spanning entire nanoparticles, and how this affects their reactivity. We then evaluate current efforts and the remaining challenges in developing robust and predictive simulations for modeling such behavior. Last, we provide our perspectives in four critical areas—integration of artificial intelligence, building robust catalysis informatics infrastructure, synergism with experimental characterization, and adaptive modeling frameworks—that we believe can help surmount the remaining challenges in rationally designing catalysts in light of these complex phenomena.

catalytic mechanisms↗