Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ramping”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Fermilab Main Injector and Recycler Operations in the Megawatt Era

Significant upgrades to Fermilab's accelerator complex have accompanied the development of LBNF (the Long Baseline Neutrino Facility) and DUNE (Deep Underground Neutrino Experiment). These improvements will facilitate 1-MW operation of the NuMI (Neutrinos at the Main Injector) beam for the first time this year through changes to the Recycler slip-stacking procedure and shortening of the Main Injector ramp time. The modifications to the Recycler slip-stacking and efforts to reduce the duration of the Main Injector ramp will be discussed. Additionally, details regarding further shortening of the ramp time and the subsequent impacts on future accelerator operations are presented.

43 PARTICLE ACCELERATORS↗
DOE OSTI2023

EVs@Scale Next-Gen Profiles - EV Profile Capture 2023

As part of the U.S. DOE EVs@Scale consortium Next-Gen Profiles (NGP) project, the profile capture of production electric vehicles undergoing high power charging (HPC) is conducted over a wide range of conditions to explore variance and performance. Charge session parameters are collected from both the electric vehicle (EV) and electric vehicle supply equipment (EVSE) at a rate of 10Hz and entered into a time-series database for analysis. These charge profiles are captured under nominal and off-nominal conditions, exploring the impact of battery state of charge (SOC), battery temperature, vehicle condition, smart charge management (SCM), and EVSE limitations. Nominal conditions are defined to be ideal conditions that should transfer the maximum allowable energy in the minimum possible amount of time. Nominal condition profiles are compared across EVs to characterize state-of-the-art EV charging performance against one another. Off-nominal condition profiles are compared against its nominal condition profile counterpart to highlight the variance across less desirable starting conditions within a single EV. Under nominal conditions, most EVs can achieve the original equipment manufacturer (OEM) rated peak-power and charge times. Peak power across EVs has variance due to the vehicle design, battery topology, charging strategy, etc. These 10-100% nominal preconditioned charge profiles across 13 EVs (11 light-duty (LD), 2 heavy-duty (HD)) were used for analysis in power curve analysis, power distribution, SOC and range comparison, thermal impacts of current draw, battery pack size and energy charged, and ramp rates. Power curve analysis show the uniqueness in power vs time curves across all EVs, breaking down the features of a typical profile and where variation is typically seen. Power distribution results showed that over 50% of charge time is spent below 50kW and only 12.1% is spent above 200kW. SOC and range performance yielded different top performing EVs when exploring goals of performance from SOC and range gained after 10-min (EV2) and 20-min (EV8), and time-to-achieve 80% SOC (EV8) and 200 miles of range (EV1). Current draw from 400-volt EVs had a higher thermal impact to cable/connector temperatures when compared to 800-volt EVs, due to a higher current requirement to achieve similar power levels. There was high variance in C rating, a useful metric when comparing the relationship between peak/average charge session power and relative battery pack size, across EVs under test. Ramp rates during initial power transfer were examined, fastest and slowest speeds ranging from 192.5kW/second to 2.6kW/second respectively. A similar analysis of ramp rates was also conducted for OCPP curtailment testing, where EVs underwent a 2-minute 65A curtailment request before returning to full-power charge. Under off-nominal conditions, most EVs experienced variation from the nominal condition profiles. EVs that underwent the full set of NGP defined testing conditions were compared, analysis of which was categorized by 800-volt and 400-volt EVs. It should be noted that the 800-volt EVs had significantly higher peak power ratings than the 400-volt EVs, and thus were more prone to variance. Initial state of charge, battery temperature, and vehicle condition had a considerable impact on peak power levels and charge time for the 800-volt vehicles under test. EVSE limited tests lowered achievable power down to 200kW, greatly impacted 800-volt charging power, and had little to no effect on 400-volt charging power. Adapter testing was performed on a single 400-volt EV, where power and current limitations were found but overall charge time was not significantly impacted. Adapter and boost converter testing was performed on a single 800-volt EV, where both charge power and charge time were greatly impacted. Charge profiles are unique, and comparing such requires analysis that examines starting conditions, vehicle and battery topologies,

Charging↗
DOE OSTI2025

Fuel Bonding and its Impact on Axial Gas Communication Behavior in Light-Water Reactor Fuel Rods

Axial gas communication concerns the flow along the axial axis of nuclear fuel rods during ramp and loss of coolant accident (LOCA) conditions. During power ramps, the higher linear heat generation rate may cause fuel-to-clad gap closure that may prevent transport of released fission gases to the plenum. Upon reduction in power the gas then can communicate to the plenum. This phenomenon has been experimentally observed by short power dips during ramp experiments completed at the Risø reactor. At higher burnups it is observed that the UO2 fuel and Zircaloy cladding forms a chemical bond. This bond results in complete closure of the gap. When these high burnup rods are subjected to a LOCA, the bond has implications on both the mechanical response (i.e., ballooning) of the cladding and subsequent fuel relocation and axial gas communication. In the LOCA scenario, gas communication is of interest in two different regimes: 1) pre-rupture communication from the plenum towards the lower pressure ballooning area and 2) the post-rupture depressurization of the plenum to the external system pressure. In both regimes the presence of a fuel-to-cladding bond will impact the rate of depressurization. In this work we present a fuel-to-clad bonding model that is coupled to an existing axial gas communication model framework in the BISON fuel performance code. The effect of considering the bond on fuel performance modeling predictions is presented through comparisons to existing experimental data. Experiments considered include several rods from the Halden IFA-650 test series. An evaluation on a full-length rod that explores the combined effect of plenum size and bonding status on axial gas communication behavior is also presented.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Evaluating the Impact of Managed EV Charging for Reliable Operation of Bulk Power Systems with High Non-Dispatchable Generation

The growth of electric vehicles (EVs) and variable-generation (VG) sources introduces new challenges for power-system operations. This study introduces a modeling framework and evaluates five EV charging strategies under projected 2040 grid conditions in the Evergy service territory with high non-dispatchable generation. Using realistic EV behavior and generation models, their impacts on system peak demand, ramp rate, and reserve capacity are evaluated. Results show that only the peak-avoidance strategy effectively reduces system peak demand, while decentralized strategies-particularly cost based dynamic charging-can exacerbate peaks due to synchronized user behavior. However, ramp-rate minimization strategy significantly reduce the stress on dispatchable generation achieving the lowest maximum absolute ramp rate (MARR) (56.81 MW) and lowest reserve requirement (2.39 GW). In contrast, unmanaged and TOU random strategies increase the stress on dispatchable generation sources with increased MARR and reserve requirements. These findings highlight the importance of coordinated, system-aware managed charging strategies to ensure reliable and affordable grid operation in the presence of EVs and VG sources.

14 - SOLAR ENERGY↗
DOE OSTI2023

Merging flux-variance with surface renewal methods in the roughness sublayer and the atmospheric surface layer

In this article, Two micrometeorological methods utilizing high-frequency sampled air temperature were tested against eddy covariance (EC) sensible heat flux (𝐻) measurements at three sites representing agricultural, agro-forestry, and forestry systems. The two methods cover conventional and newly proposed forms of the flux-variance (FV) and surface renewal (SR) schemes of differing complexities. The sites represent measurements in surface, roughness, and roughness to surface transitional layers. Regression analyzes against EC show that the most reliable FV and SR forms estimate 𝐻 with slopes within ±10% from unity and coefficient of determination 𝑅2 > 0.9 across all the three sites. The best performance of both FV and SR was found at the agricultural site with measurements well within the surface layer, while the worst was found for the tall forest with measurements within the roughness sublayer where its thickness needed to be additionally estimated. The main variable driving 𝐻 in FV is the temperature variance, whereas in SR, it is the geometry of ramp-like structures. Since these structures are also responsible for most of the temperature variance, a novel FV-SR approach emerging from combining the methods is proposed and evaluated against EC measurements and conventional FV and SR schemes. The proposed FV-SR approach requiring only a single fast response thermocouple is potentially independent of calibration and ameliorates some of the theoretical objections that arise when combining ramp statistics with similarity arguments. The combination of methods also provides new insights into the contribution of coherent structures to the temperature variance and its dependence on atmospheric stratification. Other potential utility of the new method is to include it in multi-tool assessments of surface energy fluxes, since a convergence or divergence of the results has a high diagnostic value.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

A Fast Flexibility-Driven Generation Portfolio Planning Method for Sustainable Power Systems

We are witnessing an acceleration in the uptake of renewable energy in power systems. Because of the associated variability and uncertainty of renewables, power systems need to have an adequate supply of flexibility to allow for suitable management of short-term operations. So far most of the work in this area has neglected how flexibility needs associated with renewables are fulfilled as part of dispatchable generation capital investments decisions. To address this challenge, we propose an approach to plan the dispatchable generation mix of a power system as needed to counteract variability and uncertainty associated with significant shares of variable renewable generation. The approach exploits the linear time-invariant feature of variable generation variability using historical phase planes of capacity (in MW) and ramp (in MW/min) to bridge the gap between long-term capacity planning and short-term intra-hour flexibility needs. This approach is much more computationally tractable than other proposals, while also being able to capture adequately short-term operational features like ramping and net load variability. Numerical tests are performed on realistic datasets to substantiate the effectiveness of the approach.

bulk power system planning↗
DOE OSTI2024

ANS Winter 2024 Summary: Optimizing the ATF-2Ramp Power Profile

When the Halden Boiling Water Reactor closed down in 2018, a need to restore the capability for in-reactor power ramp testing arose. Such testing is valuable for studying pellet-clad interaction phenomena in nuclear fuels. The data from these studies is of great interest to a number of research programs, including the accident-tolerant fuel (ATF) program at Idaho National Laboratory (INL). In 2022, Woolstenhulme et al. proposed several power ramp testing ideas using facilities at INL, including irradiation in the Transient Reactor Test Facility (better known as TREAT) and the Advanced Test Reactor (ATR) [1]. Worrall et al. [2] and Labossiere-Hickman et al. [3] subsequently performed feasibility studies for the ATR testing options in 2023. This summary further investigates the three-pin trefoil design (Fig. 1) for the proposed ATF-2Ramp Experiment discussed in Labossiere-Hickman et al. [3]. ATF-2Ramp is designed to operate in the center flux trap (CFT) of the ATR during a powered axial locator mechanism (PALM) cycle: a short, variable-powered cycle with an asymmetric power distribution. Previously, it was shown that tailoring the thickness of the hafnium (Hf) neutron shields (“mini-shrouds”) surrounding each pin offered a degree of control sufficient to achieve the programmatic linear heat generation rate (LHGR) targets for ATF-2Ramp during the high-power period of a PALM cycle. New work involves shortening the experiment test train for consistency with the fuel pins in ATF-2D [4] and then shaping the axial power profile of the three test pins.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Effect of Sample Mass, Confinement, and Preheating Time on the Thermal Response of LLM‐105: Experiments and Kinetic Analysis

Various small-scale experiments were performed to provide data for developing a model to predict the thermal response of LLM-105 over a wide range of conditions. The thermal decomposition of LLM-105 was studied as a function of sample mass, confinement of volatile products, and preheating time in both isothermal and ramped heating experiments. The thermal decomposition of LLM-105 is a two-step process, as shown by the two exothermic peaks in the heat flow profiles, which were fitted to two nth-order autocatalytic reaction models with a similar activation energy of ∼289 kJ/mol. The magnitude and shape of these peaks varied with sample mass and confinement. Increasing sample mass enhanced the second exotherm with respect to the first one, while increasing the level of confinement promoted a transition from a sublimation-dominated regime towards thermal decomposition. The effect of LLM-105 particle size on the rate of weight loss was evident for open-pan experiments, where bigger particles sublimed at lower temperatures than smaller particles. Thermal response and solid residue composition of LLM-105 samples were analyzed following preheating for different durations. Longer preheating times caused a shift of the second exotherm to lower temperatures and a decrease in the reaction enthalpy, confirming that LLM-105 decay is a consecutive reaction mechanism, probably autocatalytic. In conclusion, the kinetic model derived from ramped experiments was validated against the measured LLM-105 fraction remaining and the enthalpy remaining of the solid residue as a function of preheating times and showed good agreement.

Chemistry - Chemical explosives↗
DOE OSTI2022

Adversarial classification via distributional robustness with Wasserstein ambiguity

Abstract We study a model for adversarial classification based on distributionally robust chance constraints. We show that under Wasserstein ambiguity, the model aims to minimize the conditional value-at-risk of the distance to misclassification, and we explore links to adversarial classification models proposed earlier and to maximum-margin classifiers. We also provide a reformulation of the distributionally robust model for linear classification, and show it is equivalent to minimizing a regularized ramp loss objective. Numerical experiments show that, despite the nonconvexity of this formulation, standard descent methods appear to converge to the global minimizer for this problem. Inspired by this observation, we show that, for a certain class of distributions, the only stationary point of the regularized ramp loss minimization problem is the global minimizer.

Ho-Nguyen, Nam (ORCID:0000000344647730)↗
DOE OSTI2022

Multiscale simulation of integrated energy system and electricity market interactions

Accelerating the deep decarbonization of the world's electric grids requires the coordination of complex energy systems and infrastructures across timescales from seconds to decades. Here, we present a new multiscale simulation framework that integrates process- and grid-centric modeling paradigms to better design, operate, and control integrated energy systems (IESs), which combine multiple technologies, in wholesale energy markets. Traditionally, IESs are analyzed with a process-centric paradigm such as levelized cost of electricity (LCOE) or annualized net revenue, ignoring important interactions with electricity markets. This framework explicitly models the complex interactions between an IES's bidding, scheduling, and control decisions and the energy market's clearing and settlement processes, while incorporating operational uncertainties. Through two case studies, we show the importance of understanding and quantifying complex resource-grid interactions. In case study 1, we demonstrate that optimized bidding from one resource shifts the profit distribution for all energy systems in the market. This result suggests new and more flexible IES technologies can disrupt the economics of all market participants, possibly leading to accelerated retirements of less flexible resources. Interestingly, the optimized bidding has little impact on grid-level aggregate statistics, such as total generation costs and renewable penetration rate. While aggregate modeling strategies may remain valid under some IES adoption scenarios for analysis focused on regional outcomes, direct comparisons of IES technologies at specific locations without considering these interactions may lead to misleading or incorrect conclusions. In case study 2, we consider the design and flexible operation of IESs that hybridize conventional generators with energy storage. Through a sensitivity analysis, we find that as the size of the storage system increases, the total number of start-ups for coal- and natural gas-based IESs reduced by 25% and 33.6%, and the total thermal generator ramping (i.e., mileage) reduced by 86.5% and 62.5%, respectively. This shows the primary benefit of storage may not be reduced operational costs (which do not change significantly) but fewer start-ups and less ramping, which may greatly simplify the design, operation, and control of carbon capture systems. The new modeling and optimization capabilities from this work enable the coupling of rigorous, dynamic process models with grid-level production cost models to quantitatively identify the nuanced interdependencies across these vast timescales that must be addressed to realize clean, safe, and secure energy production. Moreover, the proposed general multiscale simulation framework is applicable to all IES technologies and can be easily extended to consider other energy carriers (e.g., hydrogen, ammonia) and energy infrastructures (e.g., natural gas pipelines).

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

A bilevel multistage stochastic self-scheduling model with indivisibilities for trading in the continuous intraday electricity market

In this paper, we study the profit maximization problem of a virtual power plant trading in the continuous intraday electricity market. Our virtual power plant model is compatible with renewable, and thermal assets, covering a range of virtual power plants currently participating in energy markets. We model the trading problem as a bilevel multistage stochastic program. The upper level of the problem accounts for the profit maximization of the virtual power plant with explicit modeling of the technical constraints of the operational status of the thermal power plant including minimum start-up and shut-down times, ramp-up and ramp-down rates, and minimum generation level. The upper level also decides which continuous and indivisible (fill-or-kill) orders are submitted to the market. The lower-level problem accounts for the clearing of the continuous intraday market, i.e., matching of buy and sell orders. Because of the presence of fill-or-kill orders, the lower-level problem is mixed-integer, which prevents its direct conversion to a single-level problem using duality. In order to solve this challenging problem, we develop a convex-hull extended formulation for the lower-level problem, apply duality theory to obtain a single-level stochastic equivalent formulation, and employ McCormick envelopes to turn the problem into a multistage stochastic mixed-integer linear problem, which we solve using the stochastic dual dynamic integer programming algorithm. We conduct numerical experiments and analyze the optimal trading behavior of a virtual power plant trading in an ideal continuous market without arbitrage.

Bilevel multistage stochastic programming problem↗
DOE OSTI2025

Plasma shape and position control development for NSTX-U using the GSEvolve plasma simulator

NSTX-U discharges have been reproduced using the GSEvolve plasma simulator in order to improve their control and reduce oscillations during the plasma current (Ip) ramp-up phase, and with the ultimate goal of enabling reliable H-mode access with a double-null shape at high elongation, low internal inductance, and high plasma performance. The control objective is to eliminate undesired oscillations in the plasma shape and vertical position that start when the plasma diverts during the Ip ramp-up phase, and that happened recurrently during the first NSTX-U experimental campaign. Model-based control solutions have been developed, including alternative H-mode access recipes and a Multi-Input Multi-Output (MIMO) decoupling controller. These solutions have demonstrated their capability to reduce or eliminate the aforementioned oscillations in GSEvolve simulations, enabling reproducible H-mode access. The results of this work aim to inform NSTX-U experimental operations on potential control developments to maximize its scientific output, and highlight the powerful capabilities of high-fidelity control-simulation frameworks like GSEvolve.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Thermal mechanical assessment of a SiC-SiC-composite clad fuel pin concept in a light water reactor environment

Accident Tolerant Fuels (ATFs) are designed to increase coping time following an accident scenario while preserving or improving current steady state reactor operational performance. A potential ATF concept is SiC-SiC composite claddings. Fuel performance simulations were conducted on a SiC-SiC based cladding concept utilizing a multilayered approach for improved performance. This cladding concept referred to as the Duplex concept is a duplex structure composed of a monolithic SiC layer placed on the outside of a SiC-SiC composite. A liquid metal is added to fuel-cladding gap for improved heat dissipation from the fuel. The monolithic SiC layer is used to improve the coolant corrosion characteristics and protect the SiC-SiC composite layer from exposure to the coolant. The fuel performance code BISON was used to conduct fuel performance simulations on the cladding concepts. Comparisons are made with a current prototypic fuel rod design (UO 2 fuel enclosed in Zircaloy-4 cladding). Representative steady-state cases were considered for normal power and two cycle power histories. Additionally, a PCI ramp case was simulated to analyze potential anticipated operational occurrences. Transient response during a Loss of Coolant Accident and a Reactivity Initiated Accident were also simulated. This computational study demonstrated that for normal operating conditions, the SiC concept cladding performed as well as the baseline for the standard power cases evaluated. The ramping evaluations indicate potential fracturing of the SiC-SiC composite of the composite cladding compared to the Zircaloy-4 cladding due to the temperature gradient and the subsequent differential thermal conductivity degradation and swelling across the composite thickness. In conclusion, the rod fails early at low enthalpy for RIA but survives a LOCA with minimal material loss due to high temperature steam corrosion.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Enhancement of Uranium Ionization Efficiencies Using Zn-MOF-74 Derived Nanoporous Ion Emitters for Thermal Ionization Mass Spectrometry

The recent introduction of nano-porous ion emitters (nano-PIEs) formed from metal organic frameworks (MOFs) has demonstrated the potential to enhance sensitivity for thermal ionization mass spectrometry (TIMS). Nano-PIEs take advantage of the parent MOF’s chemical and structural tunability to form scaffolds for ion emitters. A study by Barpaga et al. 2023 on MOF-74 as the parent material found that with high volatility metals in their framework uranium sample utilization efficiency (SUE) increases by up to nine-fold (e.g., Zn-MOF-74) compared to that of the analyte on a bare filament (~0.05%). Here, in this study, we investigate the performance of Zn-MOF-74 to maximize uranium efficiencies at the trace level (= 10 -12 g) by altering the parent MOF morphology and chemistry (i.e., MOF crystal sizes and thermal degradation) and optimizing its integration with TIMS (i.e., MOF mass on a filament and ramp conditions). We observed improvement in SUE up to 20 times (≤1.0%) that of a bare filament load when nano-PIEs derived from nanocrystals of Zn-MOF-74 were heated under a specific current ramp condition. This demonstrates that rates of nano-PIE structural collapse during TIMS analysis and the subsequently formed nanomaterials (and their features) can be tuned to control analyte ionization.

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

Dynamic Modeling and Simulation of a Subcritical Coal-Fired Power Plant under Load-Following Conditions

Dynamic models for power plants that capture realistic general process trends and effects of manipulated variables are needed to improve load-following, while minimizing carbon footprint. In this work, a dynamic modeling approach and simulation results for subcritical coal-fired power plant components are presented. These encompass simulation of the dynamics in the fireside, including the effects of fuel, air combustion, and the dynamics of the entire waterside and power generation sections. This model development enables the simulation and analysis of the important short and long-time scale dynamics of components such as heaters, evaporative loop, and power generation units. Furthermore, additional variables in the power generation section are introduced to improve model accuracy, extending the prediction capability of subcritical power plant models and opening new opportunities for research in operator training, optimization, and advanced model-based controller design that are based on these models. The change in process gain for different ramp rates associated with disturbance signals that affect process variables is also explored and a correlation developed. This provides opportunities to study disturbance rejection control implementation and adaptation for scenarios with such variations in ramp rates. The prediction capabilities of selected components are compared to data available in literature, with the obtained root mean squared error ranges that reflect the model performance and quality of predictions.

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

Effect of Radical Initiators on Polypropylene Thermal Deconstruction

A new route for polypropylene (PP) deconstruction through radical pathways based on thermodynamic and kinetic considerations has been investigated. Radical polypropylene (PP) deconstruction, activated through small quantities of initiators, can enable the deconstruction of waste POs into unsaturated products. We found that stirring was detrimental to the β-scission extent because it accelerates radical termination reactions through mixing. The best results were achieved by using a semibatch process that included mechanical mixing during the temperature ramp, static heating of the polymer/initiator mixture at the final temperature, and volatile product removal with N 2 . Dicumyl peroxide and alkylated dicumene initiators produced similar liquid and solid products after a 30 min thermal treatment. Terminal alkenes were formed in the liquid products of reactions at 375 and 400 °C, with about 5% of the protons in the liquid product ascribed to terminal alkenes. Yields of liquid products increased with temperature, reaching 60% w/w at 400 °C; at the same time, yields of solid products decreased with temperature to 16% w/w at 400 °C. Although radical initiators decrease PP molecular weight during the temperature ramp, at 400 °C, initiators only marginally increased the liquid product fraction compared to control experiments. Finally, the terminal double bonds in the liquid product mixture (C8–C36) provide multiple pathways for upgrading to surfactants, plasticizers, lube oils, and other valuable products.

liquids↗
DOE OSTI2025

Machine-Learned Force Field for Molecular Dynamics Simulations of Nonequilibrium Ammonia Synthesis on Iron Catalysts

Ammonia (NH 3 ) is one of the most important industrial chemicals. The conventional NH 3 synthesis method-the Haber–Bosch process-converts atmospheric nitrogen (N 2 ) into NH 3 using H 2 with an iron (Fe) catalyst. However, this process requires high pressures (100–200 atm) and temperatures (700–800 K) near thermal equilibrium. Recently, Fe-based nanocatalysts have been reported to produce promising NH 3 yields under atmospheric pressures and temperature-modulated nonequilibrium conditions. Understanding the mechanism of nonequilibrium catalysis with programmed temperature variation could help to optimize this fully electrified and less energy-intensive process. Although reactive molecular dynamics (RMD) simulations can be a useful tool to model nonequilibrium catalytic processes, they require the development of accurate force fields (i.e., interatomic potentials). Here, we present a machine-learned (ML) force field within the Deep Potential MD (DPMD) framework, trained using periodic density functional theory (DFT) calculations, to model NH 3 synthesis on Fe catalysts with various surface adsorbates such as *N, *H, *N 2 , *H 2 , *NH, *NH 2 , and *NH 3 . Here, we generated the DFT data from static models of elementary reactions on the most stable (110) surface of body-centered cubic Fe, which then were augmented by data from constant number of particles–volume–temperature (NVT) DFT-MD trajectories at various temperatures. Finally, we utilized the fully optimized ML force field to investigate reaction dynamics at an Fe(110) surface at linearly increasing temperatures using NVT-DPMD simulations. Our simulations indicate that pulsed temperature ramping could prove favorable for NH3 synthesis. For example, we conducted ramping under multiple sets of conditions: (i) from 900 to 1200 K over periods of 0.1–0.3 ns for Fe surfaces precovered with N or NH along with H; and (ii) from 300 to 600 K over 0.1–0.3 ns for Fe surfaces precovered with NH 3 . While our simulations so far are limited to short time scales (very rapid heating), these observations shed light on the mechanism of the high NH 3 synthesis rate achieved in a novel temperature-modulated nonequilibrium catalytic reactor using pulsed heating and cooling.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗