Search NASA⌕ Search

SEARCH · Search NASA

Results for “Equilibrium programming”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Nematic Colloidal Micro‐Robots as Physically Intelligent Systems

Abstract Physically intelligent micro‐robotic systems exploit information embedded in micro‐robots, their colloidal cargo, and their milieu to interact, assemble, and form functional structures. Nonlinear anisotropic fluids such as nematic liquid crystals (NLCs) provide untapped opportunities to embed interactions via their topological defects, complex elastic responses, and ability to dramatically restructure in dynamic settings. Here a four‐armed ferromagnetic micro‐robot is designed and fabricated to embed and dynamically reconfigure information in the nematic director field, generating a suite of physical interactions for cargo manipulation. The micro‐robot shape and surface chemistry are designed to generate a nemato‐elastic energy landscape in the domain that defines multiple modes of emergent, bottom‐up interactions with passive colloids. Micro‐robot rotation expands the ability to sculpt interactions; the energy landscape around a rotating micro‐robot is dynamically reconfigured by complex far‐from‐equilibrium dynamics of the micro‐robot's companion topological defect. These defect dynamics allow transient information to be programmed into the domain and exploited. Robust micro‐robotic manipulation strategies are demonstrated that exploit these diverse modes of nemato‐elastic interaction to achieve cargo docking, transport, release, and assembly of complex reconfigurable structures at multi‐stable sites. Such structures are of great interest to future developments of LC‐based advanced optical device and micro‐manufacturing in anisotropic environments.

Yao, Tianyi↗
DOE OSTI2024

Bubble Mass Transport Measurement in the Large-Scale Test Loops at Oak Ridge National Laboratory

Molten salts are complex fluids that incorporate multi-phase behavior depending on the chemistry and the physical properties of the entrained components. These components include the carrier salt, the actinide fuel, and fission, activation, and corrosion products, the concentrations of which depends on the burnup history of the salt. Radionuclide transport from the salt into the cover gas / off-gas system depends on volatility as predicted by thermochemistry, but data from the Molten Salt Reactor Experiment (MSRE) conducted in the late 1960s suggest that bubble formation and transport are also important. Anomalously high fractions of noble metals were found in the off-gas system and were attributed to transportation with parent salt aerosols and their association with rising noble gas bubbles. The prediction of such phenomena requires coupled neutronic, thermal hydraulic, and chemical equilibrium calculations, the framework of which is being developed within the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program supported by the US Department of Energy (DOE) Office of Nuclear Energy (NE). However, separate effects tests and models of experiments using tools such as SAM and Thermochimica can guide model developers through important processes. Conversely, model development guides the choice of experiments and systems to provide data that are relevant for validation. This report describes several experiments that tracked gas transport in molten salts, ranging from small-scale systems to large-scale loops. Gas transportation in a molten salt, LiCl-KCl, has been studied using the shadowgraph technique. Sensors such as residual gas analysis, Raman spectroscopy, and laser-induced breakdown spectroscopy have been tested for off-gas measurements. These data were used to interpret how the gases move through the upstream salt / cover gas and the interface between them. Differential pressure measurements were able to detect individual gas bubbles as they popped at the liquid–gas interface. Salt aerosols were collected on a cascade impactor. Their formation was also observed directly via the shadowgraph method, and most of these aerosols were launched ballistically into the cover gas. Fine mists could also be observed. Convection currents through the salt were visualized and can be used to calculate the thermophysical properties of the salt itself. The apparatuses described in this report and in a previous work (McFarlane et al. 2023) have been commissioned and are available for use in making further measurements of salt/surrogate fission product behavior. Novel approaches using neutron imaging are planned for the study of fluoride salts, which cannot be contained in quartz, so shadowgraph visualization is not available. Bubble transport in convective flow and in a slow-moving salt column are planned. The mobile laser-induced breakdown spectroscopy (LIBS) system is available for several applications, including iodine capture in a molten hydroxide scrubber, H2 transport though molten salts to complement Raman analysis, and online tracking of salt aerosol generation, transport, and deposition. This report summarizes the findings from FY24 and the plans for FY25. The work completes milestone M2AT-24OR0702013 of the DOE-NE Advanced Reactor Technology, Molten Salt Reactor Campaign, DOE-NE-5.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

SorpPropLib : An Open-Source Database for Sorption Equilibrium Properties

Hundreds of sorption working pair equilibrium measurements have been reported in the literature, and many experimental results have been fitted to equilibrium equations. However, an open database does not exist for the equilibrium equations. Therefore, this work compiled 438 published equilibrium equations of sorption working pairs into an open-source database called SorpPropLib. SorpPropLib is accessible from six programming environments commonly used in sorption cycle simulation (C++, Excel®, Python®, Matlab®, Modelica® and Lab-VIEW™). Plotting capabilities allow fast cycle visualization in each programming environment. Thereby, the open-source database SorpPropLib aims to enable and encourage easy data sharing to advance sorption research.

Engelpracht, Mirko↗
DOE OSTI2022

A parametric approach to identify synergistic domains of process intensification for reactive separation

Process intensification aims to combine multiple tasks within multi-functional units to drastically improve economic, energy or sustainability metrics of a chemical process. Limited work exists to systematically identify the synergistic domains where intensification outperforms its nonintensified counterpart. In this work, we computationally derive the synergistic domains of a reactive separation system. Specifically, we first postulate general models for both intensified and nonintensified systems. We use these models to generate data to train a ReLU-type artifical neural network (ANN). Further, the trained ReLU-NN model is formulated as a multi-parametric mixed-integer linear program (mp-MILP), and the critical regions of this mp-MILP define the synergistic feasible domains of intensification. We have derived these synergistic domains of vapor–liquid equilibrium (VLE)-based reactive separation for several industrial applications. These synergistic domains enable quick screening of properties that favor intensification.

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

Preparation and Characterization of Materials for Low- to Intermediate-Temperature CO2 Adsorption

Global carbon dioxide emissions are rising and the use of fossil fuels in several sectors are the leading causes. As global population and economies continue to grow significantly, the most practical method of lowering such emissions is to capture CO2. Although other technologies are more developed, adsorption is very promising and has attracted much attention. To ensure this technology’s success, it is essential to have suitable CO2 adsorbent materials. In this work, several new hydrotalcites (HTs) with different initial concentrations of ion precursors were prepared for the first time by the co-precipitation method—it was possible to verify that the ion concentrations influence the characteristics of the materials. The prepared HTs were characterized by thermogravimetric analysis (TG), X-Ray diffraction (XRD), surface area measurements and temperature-programmed desorption of CO2 (TPD-CO2) to relate their CO2 capture capacity to their physicochemical properties; the CO2 adsorption equilibrium isotherms were determined at 35 and 300 °C for the prepared samples, as well as for some commercial materials: magnesium oxide, calcium oxide, aluminium oxide and Zeolite 13X. After determining which materials present the best CO2 adsorption capacity, these were submitted to adsorption-desorption cycles to study their stability. The main objective of the work was to prepare and study different CO2 adsorbents for processes that are carried out at low and intermediate temperatures. From the experimental results, it was possible to conclude that the Zeolite 13X showed the best capacity at 35 °C, 3.38 mmol·g−1 (@ pCO2 = 1 bar), and a prepared calcined HT (c-HT2) was the best at 300 °C, 0.97 mmol·g−1 (@ pCO2 = 1 bar). Moreover, it seems there is an optimum initial concentration of the ions’ solutions for the tested HTs, which depends on the final application—c-HT1 showed a better capacity at 35 °C and c-HT2 at 300 °C. From the adsorption-desorption cycles—performed at 35 and 300 °C with the best materials using a magnetic suspension microbalance at 1 bar of CO2 partial pressure —, a working cyclic capacity of 2.69 mmol∙g−1 was achieved by the Zeolite at 35 °C; in turn, c-HT2 showed a working cyclic capacity of 0.79 mmol∙g−1 at 300 °C.

Figueiredo, Anabela↗
DOE OSTI2024

coh3

CoH3 (CoH ver.3) is an optical model, exciton pre-equilibrium, and Hauser-Feshbach statistical model code, which calculates nuclear reaction cross sections for medium to heavy targets in the keV to MeV energy region. This program is written in standard C++, divided into approximately 200 source and header files. CoH solves the Schroedinger equation for optical potentials defined in the code, and calculates differential elastic scattering, reaction, and total cross sections, for neutron, proton, deuteron, triton, 3He, and alpha-particle. Deformed optical potentials are solved with the coupled-channels method, in which the ground state rotational band members, or vibrational phonon states are coupled. The optical model gives particle transmission coefficients that are fed into the statistical model calculations. CoH includes the pre-equilibrium model (exciton model), the direct/semidirect capture model, and the multi-stage Hauser-Feshbach statistical decay with width fluctuation correction based on the Gaussian orthogonal ensemble. For weakly coupled levels, the DWBA (distorted wave Born approximation) method is used to calculate the direct inelastic scattering process to the excited states.

Kawano, Toshihiko↗
DOE OSTI2021

CHARMM-GUI Nanomaterial Modeler for Modeling and Simulation of Nanomaterial Systems

Molecular modeling and simulation are invaluable tools for nanoscience that predict mechanical, physicochemical, and thermodynamic properties of nanomaterials and provide molecular-level insight into underlying mechanisms. However, building nanomaterial-containing systems remains challenging due to the lack of reliable and integrated cyberinfrastructures. Here we present Nanomaterial Modeler in CHARMM-GUI, a web-based cyberinfrastructure that provides an automated process to generate various nanomaterial models, associated topologies, and configuration files to perform state-of- the-art molecular dynamics simulations using most simulation packages. The nanomaterial models are based on the interface force field, one of the most reliable force fields (FFs). The transferability of nanomaterial models among the simulation programs was assessed by single-point energy calculations, which yielded 0.01% relative absolute energy differences for various surface models and equilibrium nanoparticle shapes. Three widely used Lennard-Jones (LJ) cutoff methods are employed to evaluate the compatibility of nanomaterial models with respect to conventional biomolecular FFs: simple truncation at r = 12 Å (12 cutoff), force-based switching over 10 to 12 Å (10-12 fsw), and LJ particle mesh Ewald with no cutoff (LJPME). The FF parameters with these LJ cutoff methods are extensively validated by reproducing structural, interfacial, and mechanical properties. As a result, we find that the computed density and surface energies are in good agreement with reported experimental results, although the simulation results increase in the following order: 10-12 fsw <12 cutoff < LJPME. Nanomaterials in which LJ interactions are a major component show relatively higher deviations (up to 4% in density and 8% in surface energy differences) compared with the experiment. Nanomaterial Modeler's capability is also demonstrated by generating complex systems of nanomaterial-biomolecule and nanomaterial-polymer interfaces with a combination of existing CHARMM-GUI modules. We hope that Nanomaterial Modeler can be used to carry out innovative nanomaterial modeling and simulations to acquire insight into the structure, dynamics, and underlying mechanisms of complex nanomaterial-containing systems.

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

Analytical simulation of meander morphology from equilibrium to long-term evolution: Impacts of channel geometry and vegetation-induced coarsening

This study introduces an innovative approach to modeling meandering river morphology, integrating and investigating the effects of geometric characteristics and vegetation-induced channel coarsening. The developed comprehensive framework combines several advanced techniques: Genetic Programming for refining the scour factor of transverse bed slope, a Leaf Area Index (LAI)-enhanced analytical model for quantifying vegetative flow resistance, and an upstream-weighted moving average method for efficient approximation of the convolution integral in meander migration calculations. The model is validated against both an idealized Kinoshita meander and a natural bend of the Tumen River (China) in equilibrium, demonstrating its robustness across diverse scales and conditions. The model's ability to simulate the long-term evolution, including cutoff events, provides valuable insight for river management strategies. The current findings demonstrate that channel geometry, particularly width-to-depth ratio, plays a dominant role in meander evolution, with wider channels prone to more complex and rapid morphological changes. Vegetation effects are most pronounced in channels with moderate width-to-depth ratios, where they can significantly influence migration rates and bed topography. A combination of channel widening and deepening, coupled with strategic vegetation management, can effectively enhance navigability while maintaining channel stability in the studied Tumen River reach. Sensitivity analyses highlight the complex interplay between hydraulic conditions, sediment characteristics, and vegetation in shaping river morphology. This research advances understanding of the multifaceted nature of meandering river systems and offers practical tools for informed decision-making in river engineering and environmental management, particularly in the context of climate change and increasing anthropogenic pressures on fluvial ecosystems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Purification and Concentration of 99 Mo from a UREX raffinate

The medical isotope 99m Tc ( t ½ = 6.0 h) accounts for over 80% of isotopes used in diagnostic nuclear medicine today and is normally derived from its transient equilibrium parent 99 Mo ( t ½ = 66 h). To reduce U.S. dependence on 99 Mo derived from HEU by foreign suppliers, the NNSA/M3 program—under the American Medical Isotope Production Act—has been tasked with facilitating the work of domestic 99 Mo suppliers that do not utilize HEU. Superconducting electron linear accelerators with high-Z converter targets can generate bremsstrahlung photons and neutron fluxes that can induce photonuclear reactions and LEU fission. After sufficient production intervals, targets can be rotated out and processed while another batch is irradiated. The process flow chemistry is shown in Figure 1. The irradiated U 3 O 8 targets are retrieved and dissolved in HNO 3 ; the volatile fission products are expelled and captured during this step. The HNO 3 liquor bearing the uranyl, fission and activation products is injected into the UREX suite of liquid-liquid extraction banks. The uranyl is partitioned in typical fashion using tri-$\textit{n}$-butyl phosphate (TBP) in a hydrocarbon diluent. The raffinate—comprising mainly 99 Mo and fission products—is then injected into the MoLLE (Molybdenum Liquid-Liquid Extraction) flowsheet, where the Mo is selectively extracted by an organophosphorous acid extractant such as di-(2-ethylhexyl phosphoric acid) (HDEHP). Trace amounts of Nb, Zr, Np, I, and Te are co-extracted. The Mo is stripped using acetohydroxamic acid (AHA) and fed onto an anion exchange column. Following a series of hydroxide, HCl, and oxalic acid wash steps, the Mo can be recovered in NaOH/NaCl. A manuscript describing this process chemistry in more detail was recently published.

07 ISOTOPE AND RADIATION SOURCES↗
DOE OSTI2021

SISGR: Chemomechanics of Far-From-Equilibrium Interfaces (COFFEI)

Portable, reliable, and deployable devices for energy storage and conversion require fundamental changes in design of solid-state composites comprising ceramics and metals. These materials comprise the electrodes and electrolytes of next-generation solid-oxide fuel cells and all solid-state batteries, forming solid-state functional composites. The advent of solid-state batteries – which replace liquid electrolytes with solid electrolytes capable of lithium ion transport for reliable energy storage in portable batteries – and the increased demand space for all solid-state fuel cells capable of oxygen reduction at intermediate temperatures remain important challenges for improved material stability and decreased system cost. However, little is understood about three fundamental facets of materials that enable such solid-state energy applications. First, how do such materials deform, fracture, or delaminate under operando conditions? Second, how does such mechanical deformation limit or facilitate electronic and ionic transport within and across such material interfaces? Third, how we can predictably design interface-rich composites to engineer both structural and electrochemical stability? This COFFEI Group comprised expertise from Materials Science & Engineering and Nuclear Science & Engineering to integrate unique in situ experiments, simulations, and fabricated interfaces that address these fundamental questions in solid-state interfaces of nanoscale composites that will guide solid-state electrochemistry, transport kinetics, and mechanical deformation for nonstoichiometric materials that enable such applications. In particular, we built on COFFEI’s understanding of chemomechanical coupling among defect concentrations, ionic transport, electron transport, and stored elastic energy that is particularly acute in the far-from-equilibrium conditions typical of energy device applications. By tailoring our focus to solid-state interfaces, we addressed these important issues by (a) developing and applying advanced in situ and ex situ characterization tools to characterize model materials and interfaces synthesized with molecular-level control, under both laboratory-controlled and extreme environments representative of energy device operation; and (b) employing computational modeling and simulation frameworks to predict transport mechanisms, reactivity and stability of these model materials and interfaces under significant chemical strains typical of energy device operation. Recent progress provided insights to additional materials systems and electrochemomechanical fatigue and fracture that were not fully envisioned when the program was initiated. Specifically, in the final three years of COFFEI we pursued two integrated thrusts, with complementary focus. Thrust I focused on failure-resistant electrochemomechanical composites, while Thrust II focused on strain-modulated conductivity and reactivity across interfaces. In contrast to our initial COFFEI focus, these thrusts concentrated wholly on solid-state material interfacial interactions and included greater integration of multiscale visualization including in situ electron microscopy of strained structures/interactions and mesoscale simulations. Successful development of functionally superior and long-lived battery and fuel cell systems and stress adaptable oxides requires a deeper, fundamental understanding of the coupling among the historically important subfields of solid-state electrochemistry, transport kinetics, and mechanical deformation for nonstoichiometric metal oxide electrodes. In this program, the understanding and the application of chemomechanical coupling of defect concentrations, ionic transport, electro-catalytic activity and stored elastic energy, particularly acute in the far-from-equilibrium conditions typical of energy device applications, are being refined and implications for device operation clarified, including for miniaturized solid-state batteries and fuel cells.

36 MATERIALS SCIENCE↗
DOE OSTI2024

How numerical simulations helped to achieve breakeven on the NIF

The inertial confinement fusion program relies upon detailed simulations with inertial confinement fusion (ICF) codes to design targets and to interpret the experimental results. These simulations treat as much physics from essential principles as is practical, including laser deposition, cross beam energy transfer, x-ray production and transport, nonlocal thermal equilibrium kinetics, thermal transport, hydrodynamic instabilities, thermonuclear burn, and transport of reaction products. Improvements in radiation hydrodynamic code capabilities and vast increases in computing power have enabled more realistic, accurate 3D simulations that treat all known asymmetry sources. We describe how numerical simulations helped to guide the program, assess the impediments to breakeven, and optimize every aspect of target design. A preshot simulation of the first National Ignition Facility experiment that surpassed breakeven predicted an increased yield that matches the experimental result, within the preshot predicted uncertainty, with a target gain of 1.5. We will cover the key developments in Lawrence Livermore National Laboratory ICF codes that enabled these simulations and give specific examples of how they helped to guide the program.

Marinak, M. M. (ORCID:0009000331127696)↗
DOE OSTI2020

Data-based and secure switched cyber–physical systems

In this work, we develop a completely model-free moving target defense framework for the detection and mitigation of sensor and/or actuator attacks in cyber–physical systems with dynamics that evolve in discrete-time. We incorporate an intrusion detection mechanism based on an approximate dynamic programming technique that learns the policies for optimal regulation and optimal tracking while simultaneously defending against actuator and sensor attacks in a model-free fashion. Switching rules are leveraged to force proactive and reactive defense mechanisms as well as, guarantee the stability of the equilibrium point. Finally, as a case study, we apply the proposed moving target defense framework to a DC–DC converter that is used in electric vehicles.

42 ENGINEERING↗
DOE OSTI2020

Theory Manual for the Fuel Cycle Analysis Code REBUS

REBUS is a system of programs designed for the fuel cycle analysis of fast reactors and accelerator-driven systems. The first version of REBUS was developed in 1970 for calculating the equilibrium operating conditions of liquid metal fast breeder reactors. During the past five decades, various extensions and refinements have been made in the simulation of material motions in fuel cycle facilities and the solution of the neutronics and isotopic transmutation equations that govern the reactor behavior. The latest version complies with the standard code practices and interface dataset specifications of the committee on computer code coordination (CCCC). This report provides a detailed description of the physical processes and their mathematical models and solution methods of the current production version of REBUS. The models of the facilities and physical processes involved in fast reactor fuel cycle are first discussed. Then the mathematical formulation of these physical models is presented, followed by detailed discussions of solution methods of the resulting set of equations. Finally the output edit quantities are discussed along with their definitions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Master equation study of hydrogen abstraction from HCHO by OH via a chemically activated intermediate

The abstraction reaction of hydrogen from formaldehyde by OH radical plays an important role in formaldehyde oxidation. The reaction involves a bimolecular association to form a chemically activated hydrogen-bonded reaction complex followed by a unimolecular reaction of the complex to generate the products. The reaction rate is usually considered to be pressure-independent by assuming equilibrium between the reactants and the complex. However, our nonequilibrium calculations based on the chemically significant eigenmode of the master equation, carried out with our recently developed TUMME program, indicate that the reaction complex makes the rate constant dependent on pressure at low temperatures (T < 200 K). The calculations include anharmonicity, variational effects, and multi-dimensional tunneling. We find that the reaction rate constant reaches a low-pressure limit at pressures below 10 Torr over the whole investigated temperature range (20–1800 K), which explains why the available low-temperature experiments, which are for pressures below 2 Torr, did not observe the pressure dependence. A new extension of the TUMME master-equation program is used to explore the time evolutions of the concentrations of the OH radical and the complex under pseudo-first-order conditions. The time-dependent evolution of the concentrations of the complex at a low temperature provide direct evidence for the stabilization of the reaction complex at high pressures, and it shows the negligible role of the stabilized reaction complex at low pressures. Finally, the picture that emerges is qualitatively consistent with our previous study of the reaction of methanol with OH in that the tunneling in the unimolecular step from the complex to the products affects the phenomenological reaction rate constants differently at high and low pressures and leads to a significant pressure effect.

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

The DOE E3SM Model Version 2: Overview of the Physical Model and Initial Model Evaluation

Abstract This work documents version two of the Department of Energy's Energy Exascale Earth System Model (E3SM). E3SMv2 is a significant evolution from its predecessor E3SMv1, resulting in a model that is nearly twice as fast and with a simulated climate that is improved in many metrics. We describe the physical climate model in its lower horizontal resolution configuration consisting of 110 km atmosphere, 165 km land, 0.5° river routing model, and an ocean and sea ice with mesh spacing varying between 60 km in the mid‐latitudes and 30 km at the equator and poles. The model performance is evaluated with Coupled Model Intercomparison Project Phase 6 Diagnosis, Evaluation, and Characterization of Klima simulations augmented with historical simulations as well as simulations to evaluate impacts of different forcing agents. The simulated climate has many realistic features of the climate system, with notable improvements in clouds and precipitation compared to E3SMv1. E3SMv1 suffered from an excessively high equilibrium climate sensitivity (ECS) of 5.3 K. In E3SMv2, ECS is reduced to 4.0 K which is now within the plausible range based on a recent World Climate Research Program assessment. However, a number of important biases remain including a weak Atlantic Meridional Overturning Circulation, deficiencies in the characteristics and spectral distribution of tropical atmospheric variability, and a significant underestimation of the observed warming in the second half of the historical period. An analysis of single‐forcing simulations indicates that correcting the historical temperature bias would require a substantial reduction in the magnitude of the aerosol‐related forcing.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Development of High-Fidelity and Efficient Modeling Capabilities for Enabling Co-Optimization of Fuels and Multi-Mode Engines

The objective of this closely coordinated research project is to develop improved physical sub-models and innovative numerical algorithms to enable accurate and efficient predictions of multi-mode combustion to support EERE’s Co-Optima program. To this end, a research pro-gram is proposed that seeks • to develop accurate submodels for predicting multi-mode combustion regimes, wall-heat transfer, non-equilibrium plasma ignition, and combustion-mode transition, • to develop dynamic load-rebalancing algorithms and efficient semi-implicit time-integration schemes for enabling exascale computing, • to validate the computational submodels against experiments and direct numerical simulation (DNS) data in conjunction with sensitivity analysis, and • to support the direct transition of the resulting sub-models into other engine-simulation tools through the development of self-contained modules. All sub-models will be integrated into the open-source exascale engine-simulation platform Nek5000.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
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↗
DOE OSTI2022

Efficient Region of Attraction Characterization for Control and Stabilization of Load Tap Changer Dynamics

In this article, we study the monitoring and control of long-term voltage stability considering load tap changer (LTC) dynamics. We show that under generic conditions, the LTC dynamics admit a unique stable equilibrium. For the stable equilibrium, we characterize an explicit inner approximation of its largest region of attraction (ROA). Compared to existing results, the computational complexity of the ROA characterization is drastically reduced. We propose a quadratically constrained linear program formulation for the ROA characterization problem. In addition, we formulate a second-order cone program for online voltage stability monitoring and control exploiting the proposed ROA characterization. Finally, we demonstrate the efficacy of the proposed formulations on the ROA characterization and stability monitoring and control using a standard IEEE test system.

24 POWER TRANSMISSION AND DISTRIBUTION↗