Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “building energy code”

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

Demonstration of μCHP in Light Commercial Hot Water Applications

Internal combustion engine driven combined heat and power (CHP) systems produce power while the waste heat is recovered and used for another purpose. For µCHP systems (less than 50kWe output) overall efficiencies can reach over 90% for products currently on the market throughout the world. On a fuel basis this offers an advantage over using grid-based power and another fuel for heat. Additionally, the cost of natural gas on an energy unit basis is significantly less expensive than electric power. This difference can lead to substantial savings for the end user of a CHP system. In some parts of the world, particularly Germany and Japan, µCHP has become a common part of the heat and distributed generation equipment mix. Market analysis conducted in 2012 showed that in North America there is potential of a 5,300 units per year market for µCHP, resulting in an annual energy savings of approximately 13.5 TBtu and more than $84 million in annual customer savings. However, µCHP still has not seen much market penetration in North America. This study helps to identify some of the reasons why, address roadblocks that still need to be overcome and provide experience to established best practices. Nine commercial sites were initially identified (fitness club, surgical center, senior community, resort and spa, commercial laundry, two multi-family housing, restaurant and manufacturing facility). Systems from two different manufacturers were installed at five sites total, but due to unforeseen circumstances, only three of those sites were finally commissioned: test data were collected by Oak Ridge National Laboratory. While not all systems were installed as intended, the project resulted in the following key findings: Micro combined heat and power is an unfamiliar technology for building owners and plumbing and electrical contractors that serve the buildings most suited for µCHP. Manufacturers will need to be aggressive in providing training and providing broad support to successfully bring the technology to market.; The proper integration of µCHP into other building systems is the most critical part of a successful installation. Manufacturers will need to be able to provide the expertise to train, advise and support the mechanical contractors during the design and installation of the heat integration system.; Regulations governing µCHP varies across the country. There are a few UL standards (UL 2200 and UL 1741) that apply nationally. Emissions are governed by EPA except in California where the California Air Resources Board (CARB) has established their own standard. Requirements that dictate installation and interconnection are dependent on local building codes and local utility processes and requirements. Some out of date building codes do not offer clear guidance for these requirements, leading to confusion. This makes developing training and determining broadly applicable best practices difficult.; Installation costs of µCHP are a key factor in making an installation successful. It is critical that these costs and all national, state and local regulations be understood and accounted for by all parties involved in the installation before a project is undertaken.; A robust customer service organization is important to meet the needs of the customers. These systems require annual maintenance and service and local capability is needed to keep costs inline.

03 NATURAL GAS↗
DOE OSTI2024

Metric DBSCAN

SAND2025-11725O Metric DBSCAN is an implementation of the popular DBSCAN clustering algorithm that works in general metric spaces. DBSCAN is a clustering algorithm, a fundamental building block in machine learning. It takes a set of objects and, given some notion of distance, identifies coherent groups of objects. With Metric DBSCAN, users can provide an arbitrary function to compute distance. Nearly all existing implementations of DBSCAN restrict distance to one of a few formulations. Metric DBScan accomplishes this cleanly and efficiently. The Python source code is on Github. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Dalbey, Keith↗
DOE OSTI2020

Spitfire

Spitfire is a Python and C++ code for solving canonical flame structure problems to build tabulated chemistry models for the simulation of large-scale fires. Spitfire is a modular software with performance critical code written in C++, abstract numerical methods provided in Python, and chemistry tabulation techniques made convenient for high-level use by non-expert users. It facilitates both high-level usage and novel research in tabulated chemistry approaches. SAND2019-11791 M Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE-NA0003525.

Armstrong, Elizabeth↗
DOE OSTI2021

REopt Lite User Manual

REopt Lite evaluates the economic viability of grid-connected solar photovoltaics, wind, combined heat and power, and electric and thermal storage at commercial and small industrial sites. It allows building owners to identify the system sizes and dispatch strategies that minimize the site’s life cycle cost of energy. REopt Lite also estimates the amount of time on-site generation and storage can sustain the site's critical load during a grid outage and allows the user the choice of optimizing for energy resilience. It is primarily used to inform project development decisions and to support research on the factors that drive project feasibility for market development and policy analysis. It is available through a web interface, application programing interface, and open-source code. This user manual provides an overview of the model, including its capabilities and typical applications; inputs and outputs; economic calculations; technology descriptions; and model parameters, variables, and equations. The model is highly flexible and is continually evolving to meet the needs of each analysis. Therefore, this report is not an exhaustive description of all capabilities, but rather a summary of the core components of the model.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2022

Los Alamos National Laboratory Wetland Assessment for the Replacement Water Lines from Technical Area 48 to Technical Area 55 Project

The National Nuclear Security Administration (NNSA), a semi-autonomous agency within the U.S. Department of Energy (DOE), is proposing to take action at Los Alamos National Laboratory (LANL) near a wetland at Technical Area (TA) 48 and under a wetland at TA-55. The project proposes to install two 12-inch water lines along an alignment from TA-48 to Pecos Road on the east side of TA-55 (Figure 1). The proposed water lines will replace aging cast iron pipes that provide water for fire suppression to buildings within TA-55, and potable water to other areas in the Pajarito Corridor. NNSA has prepared this wetland assessment in accordance with 10 Code of Federal Regulations (CFR) Part 1022 Compliance with Wetland and Wetland Environmental Review Requirements (10 CFR Part 1022) (CFR 2003) which was promulgated to implement DOE requirements under Executive Order 11988 Wetland Management (EO 1977). A wetland is defined in 10 CFR 1022 as “the lowlands adjoining inland and coastal waters and relatively flat areas and flood prone areas of offshore islands,” and a base wetland as “the 100-year wetland, that is, a wetland with a 1.0 percent chance of flooding in any given year (CFR 2003).” This wetland assessment evaluates potential impacts to wetland values and functions from implementation of the proposed action, identifies alternatives to the Proposed Action, and allows for meaningful public comment.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRSDate not supplied

Modeling Radiation Sources for Experimental Validation

To build on previously completed experimental work, I have been working to implement a model of our lab’s Americium-Beryllium neutron source and Ludlum detection instruments computationally to provide secondary verification of our experimental results. This has been done using the radiation transport code MCNP6 (Monte Carlo N-Particle). I have developed files to calculate the mass ratio of various composite materials that we have exposed to the radiation source, as well as the components of the detection instruments and converted all of these into geometric coordinates that define the experimental setup. Identifying the proper energy distribution and detector efficiencies as well as the correct interpretation of the detector physics have proven to be the most difficult and time-consuming challenges of the work. It is currently too early for results based on the most updated and accurate models. With this experience, the next step is to develop another source model of a proton beam to simulate exposures on other composite targets. This will be developed based on existing experimental datasets, but not used to validate them. Each of these projects will further define the radiation shielding needs for humans in space radiation environments that astronauts will someday be exposed to beyond Low Earth Orbit.

Juliana Simon↗
DOE OSTI2022

Sensitivity Coefficients Calculated for the Prompt Neutron Decay Constant at or Near Delayed Critical [Abstract]

Experimenters at Los Alamos National Laboratory (LANL) measure the prompt neutron decay constant for many experiments at the National Criticality Experiments Research Center (NCERC) to infer reactivity and the effective neutron multiplication factor. These quantities are very important for nuclear criticality safety and validating nuclear data. Uncertainty in measures of criticality of an experimental configuration can be determined prior to physically performing the experiment by applying first order perturbation theory to Monte Carlo codes, such as MCNP®. The first order perturbation theory produces first derivatives of some nuclear parameter to nuclear data (e.g., cross section data). This first derivative is commonly referred to as a sensitivity coefficient. Currently, the MCNP® software has the capability of computing effective neutron multiplication factor sensitivity coefficients to cross section data. This work builds off of this MCNP® capability and the first order perturbation theory to provide a method of calculating sensitivity coefficients for the prompt neutron decay constant at or near delayed critical to cross section data. The prompt neutron decay constant sensitivity coefficient calculated in this work does not depend on any modification of the MCNP® source code. Prompt neutron decay constant sensitivity coefficient calculations can be used to infer reactivity and effective neutron multiplication factor sensitivity coefficient values as well. By investigating the trends of prompt neutron decay constant sensitivity coefficients for nuclide-reaction pairs across energy spectra, experiments can be designed to maximize or minimize the uncertainty in the prompt neutron decay constant in a particular energy region, which can lead to further optimization studies. Prompt neutron decay constant sensitivity coefficients will be calculated for the Jezebel benchmark. Subsequently, these sensitivity coefficients will be used in a data assimilation process to determine if there is or are optimal experiments that can be performed to provide insight into adjustments of uncertain/inaccurate cross section data. Specifically, the effect of the prompt neutron decay constant sensitivity coefficients on the nuclear data-induced uncertainty in the effective neutron multiplication factor will be examined

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2015

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL Multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. This work builds upon previous efforts to verify the use of the acoustic velocity potential equation (AVPE) laid out by Campos. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, del squared psi - (lambda/c) squared psi - M x [M x del((del)(psi))] - 2((lambda)(M)/c + M x del(M) x (del)(psi) - 2(lambda)(psi)[M x del(1/c)] = 0. with M as the Mach vector, c as the speed of sound, and λ as the complex eigenvalue. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, axial velocity, exit velocity). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The one way coupled analysis is perform four times utilizing geometries determined through traditional SRM modeling procedures. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, S. R.↗
DOE OSTI2026

Heat Pump Retrofits for Central Plant Hydronic Heating Systems: A Software Toolkit for Screening and Design

Retrofitting existing central plants with high-efficiency heat pump technologies can play a crucial role in achieving long-term planning goals. Modern heat pump technologies are able to use waste heat recovery to meet a building's heating demand, but there is a lack of accessible tools designed for non-HVAC experts, such as building owners, to quickly and easily conduct what-if analysis, e.g., estimating retrofit costs and payback period for their partial or full equipment replacement. This paper introduces an open-source software toolkit designed to facilitate the initial screening and decision-making of heat pump retrofits in existing central plants using a building's yearly load profile from metered or utility bill data. The toolkit evaluates the technical and economic viability of replacing traditional central plant equipment with various options including water-to-water or air-to-water heat pumps, which can provide efficient and lower-cost heating and cooling. It allows users to compare current central plant configurations with retrofit scenarios, assessing energy consumption, life-cycle costs, and environmental impact. The toolkit offers (1) a web-based tool designed for user-friendly access by a broad audience and (2) Python-based source code for researchers and engineers conducting parametric studies and design parameter optimization. The toolkit compares a typical central plant configuration to a configuration that uses a heat pump to supply hydronic heating and cooling. The output metrics include energy consumption and output of each equipment, life-cycle cost analyses and metrics, and environmental impact of the system.

Excell, L↗
NASA NTRS1993

Collisions of small Spacewatch asteroids with the Earth

The discovery, with the Spacewatch Telescope, of 8 Earth-approaching objects smaller than 100 m in diameter has been reported. I have calculated the probability and velocity of collision with Earth for each of these objects. Using a code we have successfully employed to model the Tunguska and Revelstoke atmospheric explosions, I have simulated the entry of these objects into Earth's atmosphere, assuming iron, stony, carbonaceous, and cometary compositions. The smallest of these objects, with tens to hundreds of kilotons of kinetic energy, pose a substantial threat at the surface only if they are iron objects. An object is taken to 'pose a substantial threat' if it either craters the ground, or explodes in the atmosphere with sufficient energy at low enough altitude to create an overpressure at the surface capable of felling trees or destroying buildings. Larger objects, with energies greater than about a megaton, devastate the surface regardless of whether they are of iron, stony, or carbonaceous composition. Iron objects crater the ground, whereas stony and carbonaceous objects explode in the atmosphere low enough to fell trees and damage buildings over thousands of square kilometers. Spacewatch objects, if presumed to be of carbonaceous composition, are as dangerous as their stony counterparts, as the former objects' lower yield strengths (and hence, higher explosion altitudes) are roughly compensated by their larger masses (as derived from their lower albedos for a given observed magnitude). Although comets are intrinsically less dangerous than asteroids, the 90 m diameter Spacewatch objects would devastate hundreds of square kilometers at the surface, even if cometary.

Chyba, C. F.↗
NASA NTRS1991

A study of physical mechanisms for filament eruption and coronal mass ejection via numerical simulation

It is well-known from both observation and theory that photospheric shear motion has played a key role in physical processes of the energy build-up and release for the solar flare. In order to further our understanding of the occurrence of solar flares it is necessary to investigate the triggering mechanism. One popular scenario for the onset of eruptive solar flares is that in response to photospheric shear motions the magnetic field evolves slowly through a series of magnetohydrodynamic-equilibria until a threshold is reached where magnetohydrodynamic (MHD) non-equilibrium sets in. Thus, a magnetic eruption occurs, causing the solar flare. To substantiate this claim we have employed our newly developed three-dimensional time-dependent MHD code with gravity to simulate the evolution of the coronal field. We use plasma beta = 0.1 to closely approximate the condition in the actual corona. Some preliminary results are presented.

Wu, S. T.↗
DOE OSTI2025

Gerischer Electrochemistry Today

Semiconductor photoelectrochemistry is a dynamic and interdisciplinary field at the forefront of research in solar fuels, energy conversion, and catalysis. Here, this Perspective captures the collective insights from the second Gerischer Electrochemistry Today Symposium, held at Colorado State University in Fort Collins, CO, in August 2024, which convened leading researchers, early-career scientists, and industry partners to define the critical next steps for the field. Through interactive sessions, technical talks, panel discussions, and training initiatives─including a Semiconductor Electrochemistry Bootcamp─the symposium emphasized three pillars of advancement: (i) facilitating the exchange of new ideas in semiconductor electrochemistry and charge separation; (ii) fostering the development of future researchers, research topics, and participation in the semiconductor workforce; and (iii) building community. This Energy Focus distills key themes from the meeting and identifies major knowledge gaps in the following areas: mechanisms of charge separation and recombination, role of defects and disorder, dynamic and operando characterization methods, interfacial chemistry and surface passivation, theoretical and modeling limitations, and standardization and benchmarking. The inclusive and collaborative structure of the symposium enabled the generation of this comprehensive report that will serve as a roadmap for fundamental and applied research in the rapidly evolving field of semiconductor electrochemistry over the next decade.

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

Status of New Models Hosted on the Virtual Test Bed (VTB) in 2024

The National Reactor Innovation Center (NRIC) mission is to support deployment of novel reactor concepts. This is achieved by providing physical and virtual spaces for building and testing various components, systems, and complete pilot plants. The Virtual Test Bed (VTB) represents the virtual counterpart to the physical test bed. It is in development in collaboration with the Department of Energy’s (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. The mission of the VTB is to accelerate the deployment and licensing of advanced reactors by leveraging state-of-the-art modeling and simulation (M&S) tools developed by the DOE NEAMS program. This is accomplished by three primary means: (1) openly hosting simulations that showcase analysis capabilities, (2) continuously testing the models hosted against code updates to avoid deprecation, and (3) filling key M&S gaps that are relevant for the physical NRIC test beds. The VTB repository consists of two sub-entities: 1. A documentation website detailing the models (https://mooseframework.inl.gov/virtual_test_bed). 2. A GitHub repository that hosts the corresponding files (https://github.com/idaholab/virtual_test_bed). Previous documentation on the models hosted in the VTB can be found in [1,2,3,4]. These references also include additional background information on the various NEAMS codes showcased in the VTB (which is omitted here for brevity). This paper primarily provides a status update of the most recent additions to the repository.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Core performance predictions in projected SPARC first-campaign plasmas with nonlinear CGYRO

This work characterizes the core transport physics of SPARC early-campaign plasmas using the PORTALS-CGYRO framework. Empirical modeling of SPARC plasmas with L-mode confinement indicates an ample window of breakeven (Q > 1) without the need of H-mode operation. Extensive modeling of multi-channel (electron energy, ion energy, and electron particle) flux-matched conditions with the nonlinear CGYRO code for turbulent transport coupled to the macroscopic plasma evolution using PORTALS reveals that the maximum fusion performance to be attained will be highly dependent on the near-edge pressure. Stiff core transport conditions are found, particularly when fusion gain approaches unity, and predicted density peaking is found to be in line with empirical databases of particle source-free H-modes. Impurity optimization is identified as a potential avenue to increase fusion performance while enabling core-edge integration. Extensive validation of the quasilinear TGLF model builds confidence in reduced-model predictions. The implications of projecting L-mode performance to high-performance and burning-plasma devices is discussed, together with the importance of predicting edge conditions.

Rodriguez-Fernandez, P. (ORCID:0000000273611131)↗
DOE OSTI2024

Infrastructure improvements in the National Reactor Innovation Center Virtual Test Bed

The National Reactor Innovation Center’s (NRIC) mission is to support deployment of novel reactor concepts. This is achieved by providing physical and virtual spaces for building and testing various components, systems, and complete pilot plants. The Virtual Test Bed (VTB) represents the virtual counterpart to the physical test bed. The VTB is being developed in collaboration with the Department of Energy’s (DOE) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. The mission of the VTB is to accelerate the deployment and licensing of advanced reactors by leveraging state-of-the-art modeling and simulation (M&S) tools developed by the DOE NEAMS program. This is accomplished via three primary means: (1) openly hosting simulations that showcase analysis capabilities, (2) continuously testing the models hosted against code updates to avoid deprecation, and (3) filling key M&S gaps that are relevant for the physical NRIC test beds. The VTB repository consists of two sub-entities: • A documentation website detailing the models: https://mooseframework.inl.gov/virtual_test_bed • A GitHub repository that hosts the corresponding files: https://github.com/idaholab/virtual_test_bed. This paper presents the infrastructure added to the VTB in the last two fiscal years. Additional information about the VTB can be found in various publications [1, 2, 3, 4, 5, 6, 7, 8].

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Case Study: Leveraging GenAI to Build AI-based Surrogates and Regressors for Modeling Radio Frequency Heating in Fusion Energy Science

This work presents a detailed case study on using Generative AI (GenAI) to develop AI surrogates for simulation models in fusion energy research. The scope includes the methodology, implementation, and results of using GenAI to assist in model development and optimization, comparing these results with previous manually developed models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

MOOSE Reactor Module Meshing Enhancements to Support Reactor Physics Analysis

The U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) program develops an integrated suite of advanced reactor physics tools built upon the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework. Each code generally requires an input finite element mesh on which the physics solution is calculated, reported, and transferred to other physics codes. The meshing process is often burdensome for the complex geometries present in reactors due to lack of easy-to-use, open-source meshing tools. To address the bottleneck associated with meshing complex geometries found in nuclear reactors, NEAMS initiated the development of the MOOSE Reactor Module starting in FY21. The Reactor Module builds off the existing MOOSE Mesh System to include targeted meshing capabilities such as the ability to generate hexagonal pin cells, assemblies with ducts, rotating control drums, cores, peripheral zones around a core, as well as the automatic labeling (“reporting IDs”) of pin, assembly, and planar zones to simplify post-processing of results. As a Physics Module in MOOSE, the Reactor Module is open-source, available with any MOOSE installation, directly compatible with MOOSE-based tools, and can be invoked from MOOSEbased applications to generate meshes. Functionality from the Reactor Module has been applied to several advanced reactor concepts to demonstrate user workflow improvements and accuracy. The primary objective of the Reactor Module is to improve useability of MOOSEbased tools by streamlining mesh generation and output inspection processes. During FY22, the functionality of the Reactor Module (and accompanying Mesh System) has been expanded based on user needs. First, the Reactor Geometry Mesh Builder capability developed primarily in FY21 has been refactored and merged to the public MOOSE repository. This capability wraps underlying Reactor Module mesh generators into a “Pin – Assembly – Core” workflow appropriate for conventional Cartesian and hexagonal geometries, and notably assigns material IDs during mesh generation stage and generates only the minimal number of blocks needed in order to reduce computational burden. Biasing and boundary layer options have been added to the base mesh generators as required by thermal hydraulics solvers. The reporting ID functionality has been expanded to differentiate ring-wise and azimuthal sectors within a pin for use with depletion algorithms, and VectorPostProcessor and Reporter objects are now available to integrate solution variables across zones based on ID combinations. Functionality to trim hexagonal meshes along the center or periphery has been developed so users may leverage symmetry and reflective boundary conditions to reduce the mesh size. A flexible and powerful tool to fill the space between two sidesets has been introduced to the framework and can be used for transition layers such as stitching two assemblies together with different numbers of pins, or for complex geometries which do not follow conventional Cartesian/hexagonal patterns. Finally, additional verification problems were performed with NEAMS physics tools in complement with existing NEAMS work.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

GPR_calculator: An on-the-fly surrogate model to accelerate massive nudged elastic band calculations

We present GPR_calculator, a package based on Python and C++ programming languages to build an on-the-fly surrogate model using Gaussian Process Regression (GPR) to approximate computationally expensive electronic structure calculations. The key idea is to dynamically train a GPR model during the simulation that can accurately predict energies and forces with uncertainty quantification. When the uncertainty is high, the costly electronic structure calculation is performed to obtain the ground truth data, which is then used to update the GPR model. To illustrate the effectiveness of GPR_calculator, we demonstrate its application in Nudged Elastic Band (NEB) simulations of surface diffusion and reactions, achieving 3-10 times acceleration compared to pure ab initio calculations. The source code is available at https://github.com/MaterSim/GPR_calculator.

Gaussian process regression↗