Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simplified analysis”

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

Analytical and Numerical Analysis of Lithium Plating Onset in Single and Bilayer Graphite Electrodes during Fast Charging

Lithium plating is commonly observed in anodes charged at fast rates, and can lead to capacity loss and battery safety issues. The increased risk of plating has been attributed to transport limitations, and architectured electrodes may reduce plating risk. However, while theoretical studies have shown that reaction non-uniformity arises due to interplay of transport limitations, anode open circuit voltage behavior and reaction kinetics, its effect on lithium plating has not been studied. We use analytic and numerical simulations to predict onset of plating in graphite anode half-cells at high C-rates and demonstrate how anodes with layered porosities can delay plating. Simplified analytical models identify trends for plating onset and predictions are calibrated against numerical models. A calibrated numerical model of graphite demonstrates qualitative agreement with analytical model predictions. This reaction inhomogeneity mechanism occurs in the absence of lithium ion depletion, indicating that these mechanisms may contribute to capacity loss independently or simultaneously. A bilayer model of graphite exhibits delayed plating onset, and an optimization procedure is presented. This theoretical work presents quantitative and mechanistic insight on how reaction inhomogenity affects lithium metal plating onset and can be used as a guide to engineer anodes resistant to lithium plating.

Mijailovic, Aleksandar S.↗
DOE OSTI2024

Towards a data-driven model of hadronization using normalizing flows

We introduce a model of hadronization based on invertible neural networks that faithfully reproduces a simplified version of the Lund string model for meson hadronization. Additionally, we introduce a new training method for normalizing flows, termed MAGIC, that improves the agreement between simulated and experimental distributions of high-level (macroscopic) observables by adjusting single-emission (microscopic) dynamics. Our results constitute an important step toward realizing a machine-learning based model of hadronization that utilizes experimental data during training. Finally, we demonstrate how a Bayesian extension to this normalizing-flow architecture can be used to provide analysis of statistical and modeling uncertainties on the generated observable distributions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

MHD flow in liquid metal blankets: Major design issues, MHD guidelines and numerical analysis

The design of breeding blankets represents the major challenge for fusion reactor engineering because of performance requirements and severe operating conditions in terms of heat load and neutron flux. Liquid metal alloys such as lead-lithium, due to their lithium content, can be used to breed tritium, one of the plasma fuel components, and owing to their high thermal conductivity, they may serve as coolants. On the other hand, there are technical issues related to the fact that the liquid metals are electrically conducting and interact with the plasma-confining magnetic field. Induced electric currents and generated electromagnetic forces affect velocity and pressure distribution in the blankets. Magnetohydrodynamic (MHD) flows for fusion applications have been often investigated in simplified geometries, such as pipes, ducts, bends, with focus on their fundamental features. These analyses are essential, since results remain valid as background for the development of blanket designs, even when a concept is dismissed. However, the conceptual study of fusion blankets requires to take into account the global multiple effects, that arise when the full system is considered. Progress made in fusion-related MHD research results from combined numerical and experimental activities. In this paper we review and summarize features of 2D and 3D MHD flows that are typical in liquid metal blankets, together with available correlations for MHD pressure losses. This knowledge can provide simple design MHD guidelines that support a preliminary estimate of MHD effects in a blanket concept, in terms of pressure drop and flow distribution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Analytic results for scalar-mediated Higgs boson production in association with two jets

Abstract We present compact analytic formulae for all one-loop amplitudes representing the production of a Higgs boson in association with two jets, mediated by a colour triplet scalar particle. Many of the integral coefficients present for scalar mediators are identical to the case when a massive fermion circulates in the loop, reflecting a close relationship between the two theories. The calculation is used to study Higgs boson production in association with two jets in a simplified supersymmetry (SUSY) scenario in which the dominant additional contributions arise from loops of top squarks. The results presented here facilitate an indirect search for top squarks in this channel, by a precision measurement of the corresponding cross section. However, we find that the potential for improved discrimination between the Standard Model and SUSY cases suggested by the pattern of results in the one- and two-jet samples is unlikely to be realized due to the loss in statistical power compared to an inclusive analysis.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

RELAP-7 Application and Enhancement for FLEX Strategies and ATF Behavior under Extended Loss of AC Power Conditions

This report summarizes the results of a three-year research project sponsored by the U.S. Department of Energy (DOE) Nuclear Energy University Program (NEUP) to enhance and apply the RELAP-7 code by adding and improving several important components (e.g., a mechanistic Reactor Core Isolation Cooling (RCIC) system model) for thermal hydraulic studies of LWRs under ELAP conditions and evaluating the time available for transition to portable FLEX equipment. The project team included University of Massachusetts–Lowell, The Ohio State University, Texas A&M University, Idaho National Laboratory and Oak Ridge National Laboratory. In the Fukushima accident, it was found that the RCIC system played a crucial role in delaying core meltdown by almost three days in Fukushima Daiichi Unit 2, because of self-regulated operation of the steam driven RCIC turbine-pump injection system. Steam flow in the convergent-divergent nozzles of the RCIC Terry turbine is two-phase non-equilibrium transonic flow with homogenous nucleation condensation. To more accurately predict the dynamic process and behavior of the transonic compressible steam flow, a one-dimensional transient two-phase analytical model is presented. A simplified four-fluid model was employed in the present work with the consideration of four separate fluid fields: vapor, liquid film, entrained droplets and condensed droplets. The mass, momentum and energy interactions between the fluids were considered and modeled. An extended seven-equation non-equilibrium critical flow model was developed to obtain the critical pressure and velocities of each phase at the nozzle throat. To predict the wetness in the divergent section, a mechanistic nucleation condensation model was integrated in the nozzle analysis model, considering the generation and consequent growth of droplets. The governing differential equations on a staggered grid were discretized using the second-order Lax-Wendroff scheme with a flux limiter, and the Semi-Implicit Method for Pressure-Linked Equation (SIMPLE) algorithm was employed to solve the discrete linear system. To demonstrate the predictability and reliability of the physical models and the numerical method proposed in the present work, three representative nozzles were modeled and simulated. The results show good agreement with the available experimental data, even for condensation shock. Then, the 1D nozzle model was employed to obtain nozzle flow tables of the Terry turbine nozzle for different working pressures which can cover the operation pressure range of the RCIC system. A mechanistic RCIC turbine-pump system model was developed and implemented in the system code TRACE to simulate dynamic responses of the RCIC system under Beyond Design Basis Accident (BDBA) conditions. The turbine-pump governing equations are based on the control volume approach of the angular momentum balance. The physics based mechanistic RCIC model was developed using the TRACE control system components (i.e., signal variables, control blocks, and tables), and incorporated into a TRACE boiling water reactor (BWR) model. The TRACE model in this report has a detailed nodalization of the reactor pressure vessel (RPV), and all of the major flow paths and system components, including the safety relief valves (SRVs) and the containment suppression pool and drywell. Based on the nozzle flow tables generated from the 1D nozzle model developed, the turbine drive torque can be calculated from table lookup. Since the detailed specifications of the RCIC pump are unavailable, the homologous curves for a Bingham pump were used in the current pump component. A station black-out (SBO) accident test problem was selected to demonstrate the TRACE RCIC model. The short-term SBO simulations were performed for two cladding materials: Zircaloy and FeCrAl, to demonstrate the effect of the accident tolerant fuel cladding on fuel heat-up under BDBA accident conditions. The wetwell plays a vital safety role in SBO and other BWR accident scenarios in that it can reduce containment pressure and supply additional core make-up water. The suppression pool temperature distribution has a very large impact on both RPV and containment pressure. Thus, another novel contribution of the project comes mainly from an improved, systems-level wetwell model which can capture buoyancy-induced thermal stratification effects due to steam injection and condensation. A two-zone stratified wetwell model has been implemented in RELAP-7 and some results from that model are presented. This wetwell model is capable of simulating thermal stratification due to a low steam mass injection rate. With a low mass flow rate, the model assumes that all the steam condenses within the pipe and the resulting plume can be approximated with a purely buoyant, heat-source driven model. The wetwell model developed with these assumptions is adequate to simulate slow transients such as extended SBO transients.

42 ENGINEERING↗
DOE OSTI2021

A field–particle correlation analysis of a perpendicular magnetized collisionless shock

Using the field–particle correlation technique, we examine the particle energization in a three-dimensional (one spatial dimension and two velocity dimensions; 1D-2V) continuum Vlasov–Maxwell simulation of a perpendicular magnetized collisionless shock. The combination of the field–particle correlation technique with the high-fidelity representation of the particle distribution function provided by a direct discretization of the Vlasov equation allows us to ascertain the details of the exchange of energy between the electromagnetic fields and the particles in phase space. We identify the velocity-space signatures of shock-drift acceleration of the ions and adiabatic heating of the electrons arising from the perpendicular collisionless shock by constructing a simplified model with the minimum ingredients necessary to produce the observed energization signatures in the self-consistent Vlasov–Maxwell simulation. We are thus able to completely characterize the energy transfer in the perpendicular collisionless shock considered here and provide predictions for the application of the field–particle correlation technique to spacecraft measurements of collisionless shocks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Localized Evaluation for Constructing Discrete Vector Fields

Topological abstractions offer a method to summarize the behavior of vector fields, but computing them robustly can be challenging due to numerical precision issues. One alternative is to represent the vector field using a discrete approach, which constructs a collection of pairs of simplices in the input mesh that satisfies criteria introduced by Forman's discrete Morse theory. While numerous approaches exist to compute pairs in the restricted case of the gradient of a scalar field, state-of-the-art algorithms for the general case of vector fields require expensive optimization procedures. This paper introduces a fast, novel approach for pairing simplices of two-dimensional, triangulated vector fields that do not vary in time. The key insight of our approach is that we can employ a local evaluation, inspired by the approach used to construct a discrete gradient field, where every simplex in a mesh is considered by no more than one of its vertices. Specifically, we observe that for any edge in the input mesh, we can uniquely assign an outward direction of flow. We can further expand this consistent notion of outward flow at each vertex, which corresponds to the concept of a downhill flow in the case of scalar fields. Working with outward flow enables a linear-time algorithm that processes the (outward) neighborhoods of each vertex one-by-one, similar to the approach used for scalar fields. Here, we couple our approach to constructing discrete vector fields with a method to extract, simplify, and visualize topological features. Empirical results on analytic and simulation data demonstrate drastic improvements in running time, produce features similar to the current state-of-the-art, and show the application of simplification to large, complex flows.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Characteristics of Vertical Ground Motions and Their Effect on the Seismic Response of Bridges in the Near-Field: A State-of-the-Art Review

Despite the evidence from past earthquakes and several numerical investigations demonstrating the detrimental impact of vertical ground motions (VGMs) on the integrity of bridge structures, incorporating their effects into seismic assessment and design procedures has traditionally been given limited consideration. Current codes utilize rather simplistic approaches to account for the concurrent effects of vertical and horizontal motions in structural performance evaluations, potentially leading to unconservative estimates of structural demands. This paper reviews the main features of VGMs and their effect on the seismic response of bridges. The methods and empirical models available to estimate vertical motions for design purposes are discussed, and research gaps and related research needs are identified. Finally, the emerging role of physics-based ground-motion simulations, as well as their limitations, in supporting future research and informing the development of simplified design procedures is examined. The main areas of interest for future research are identified in the need to carry out systematic sensitivity studies to gain insight into the main earthquake parameters that influence key VGMs features, understand the influence of soil nonlinearities on VGMs amplitude and frequency content, inform the development of empirical models that cover a range of site conditions and source-to-site distances where current models are poorly constrained, generate arrays of motions to update coherency models to properly inform the analysis of distributed infrastructure, investigate the impulsive character of VGMs, and assess the approximations made in estimating VGMs with 1D site response analyses. Furthermore, specific focus is laid on large-magnitude earthquakes in the near-field.

Asynchronous ground motion↗
DOE OSTI2022

Bay Area Regional Energy: Network Integrated Commercial Retrofits (BRICR) Project. Final Report

The BRICR project applied large-scale building energy modeling concepts with the aim of reducing the cost of energy efficiency targeting, design, and project development, and measurement of energy savings for energy efficiency programs implemented by local governments that serve small and medium commercial buildings (SMB). The project leveraged the services and resources of existing local government energy programs serving disadvantaged and hard-to-reach SMB customers. In contrast to programs run by utilities, local government programs generally do not have direct access to energy billing records for an entire class of customers in a geographic area, which prior research demonstrated useful for large-scale building energy model baseline development and calibration. , However, local governments are rich in public records that offer important clues about physical attributes and uses that, along with behavior, determine energy use. Relying only on public records, BRICR demonstrated development of credible baseline energy models for 3,792 office, retail, and hotel buildings. Publicly disclosed annual energy use data from a local energy benchmarking program and anonymized data from the Building Performance Database, the nation’s largest dataset about energy-related characteristics of buildings, were utilized to validate and calibrate energy models via an innovative method comparing distributions of energy intensity by fuel type for portfolios of buildings of similar size, vintage, and use. Portfolio calibration does not provide certainty that an energy model fits an individual building; the method is useful when billing data is not accessible – a common situation for researchers, energy service providers and ESCOs, local governments, and any party other than a utility. A software component was developed, the BRICR gem, which automates simulation when relevant data is added or edited by the user to a file saved in the standardized BuildingSync XML schema for energy audit data. The component was demonstrated as a simplified means to generate a mass of energy models corresponding to public records containing basic attributes such as building scale, location, use, year built, and aspect ratio in combination with building energy code prototype data corresponding to use and vintage. The component was also demonstrated as a simplified means to automate energy simulation when attributes are revised; the intention was to enable iterative improvement of the baseline model and energy savings estimates for common energy conservation measures as users revise relevant attributes based on their observations. In the context of institutional change and uncertainty for the participating local government energy programs, 13 whole building retrofits were completed. Impacts were measured by applying the CalTRACK2.0 methods to standardize measurement of normalized metered energy consumption. The GRIDMeter methods of stratified sampling and individual load shape analysis were applied to adjust for impacts of the effect of COVID-19 on retrofitted buildings in the context of all local buildings of similar size and use. Excluding impacts of the pandemic, retrofitted buildings demonstrated between 1.6% and 25.1% reduction in energy use. The project contributed use cases and feedback that helped inform evolution of the software tools and data formats that were combined for the first time in the BRICR project.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2025

2025 Advances in NekRS: Supporting improved performance for nuclear applications

This report presents several 2025 advancements in NekRS, a high-fidelity spectral element CFD code developed at Argonne National Laboratory to support the NEAMS thermal-hydraulics program. The forthcoming v25 release consolidates several of these advances, adding new features for portability across heterogeneous GPU architectures, real-time in situ visualization, improved turbulence modeling, and conjugate heat transfer coupling. Over the past year, NekRS has demonstrated strong scalability and performance on DOE’s leading exascale platforms, including Aurora and Frontier, confirming its readiness for some of the largest and most complex simulations attempted to date. These achievements provide a powerful new platform for high-fidelity data generation, which in turn supports the development and validation of advanced closure models critical for reactor safety and design. Significant algorithmic innovations have also been introduced. A new global runtime h-refinement capability simplifies workflows by reducing mesh preparation burdens and enabling coarse-to-fine restarts. Building on this, a novel multigrid strategy was implemented to accelerate pressure and transport solves at scale, addressing long-standing bottlenecks in exascale CFD. Together, these developments improve both the efficiency and accessibility of high-fidelity simulations for reactor-relevant problems. Collectively, these enhancements represent a major step forward in simulation technology, positioning NekRS as a cornerstone of NEAMS efforts to enable accurate, efficient, and scalable high-fidelity analysis of advanced nuclear systems.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Transient kinetic insights into selective propene oxidation over industrial bismuth molybdate catalysts

Selective oxidation of propene to acrolein over industrial multicomponent bismuth molybdate (BMO) catalysts significantly depends on reaction conditions that include operating parameters and catalyst state. Here, this work investigates selective propene oxidation in the intrinsic kinetic paradigm of Temporal Analysis of Products (TAP) reactor by systematically varying catalyst redox state, temperature, and oxygen-to-propene feed ratio. A 93 % propene conversion with an acrolein yield of 80 % is achieved at elevated temperatures (450 °C) on oxidized catalysts under oxygen-rich conditions (O 2 :C 3 H 6 = 10). However, these conditions diminish acrolein-to-CO 2 selectivity due to enhanced total oxidation to CO 2 . In contrast, a reduced catalyst state, moderate temperature (350 °C), and lower oxygen feed (O 2 :C 3 H 6 = 1) nearly doubles the acrolein to-CO 2 selectivity, albeit at a lower acrolein yield (33 %). Transient kinetic studies together with a kinetic model that simplify the major product formation pathways in lumped non-elementary forms reveal that the availability of surface oxygen species plays a pivotal role in governing reaction pathways. Additionally, density functional theory (DFT) calculations on pure BMO catalysts inform the role of surface redox states on propene and oxygen activation barriers. Lattice oxygen at acrolein-selective sites drives both acrolein and CO 2 formation, while adsorbed oxygen at activation sites favors unselective CO 2 generation. This work establishes a critical relationship between transient product selectivity and surface oxygen availability, which is strongly influenced by catalyst redox state, feed ratio, and reaction temperature. These insights underscore the importance of dynamic reactor operation strategies and offer a foundation for designing next-generation propene oxidation processes with tunable acrolein selectivity.

42 - ENGINEERING↗
DOE OSTI2021

Mixed Tin-Titanium Oxides by Atomic Layer Deposition on Planar Substrates: Physical and Electronic Structure

Ternary atomic layer deposition of Sn-doped TiO x (Sn:TiO x ) and Ti-doped SnO x (Ti:SnO x ) was performed using tetrakisdimethylamidotitanium (IV) and tin(IV) with water at 150 °C on (1 0 0) Si and quartz substrates. The physical structure of the films was investigated using scanning electron, atomic force, and Kelvin probe force microscopies as well as spectroscopic ellipsometry, Raman spectroscopy, and X-ray diffraction. All as-deposited films are amorphous and possessed flat and conformal topographies with low surface roughness (<1 nm). Annealed Sn:TiO x and Ti:SnO x crystallize into the tetragonal anatase TiO 2 and cassiterite SnO 2 structures, respectively, without large variations to the film topography. The chemical composition of the as-deposited and annealed films was investigated through X-ray photoelectron spectroscopy. While dopants are homogenously mixed, the ternary ALD to form Sn:TiO x and Ti:SnO x surprisingly does not follow the rule of mixtures, despite similarities in both the growth rates of binary SnO x and TiO x and the precursors used. This finding is also supported by spectroscopic ellipsometry. Here, we postulate the deviation from the rule of mixtures is due to non-ideal surface chemistry of the binary ALD of SnO x . The optical bandgap of all films was quantified by the Tauc analysis and the valence band position was investigated by valence band X-ray photoelectron spectroscopy. Intra-bandgap states above the valence band maximum were observed for as-deposited Ti:SnO x and Sn:TiO x and ascribed to 5 s states from Sn 2+ ions. The presence of these intra-bandgap states induces enhanced visible absorption and bathochromic shifts to the apparent bandgap for all as-deposited films. Annealing the films removes the intra-bandgap states and increases the optical bandgap. Simplified band diagrams created using the data available demonstrated that the conduction and valence band positions did not vary between the positions of pure TiO 2 or SnO 2 with monotonic increases to the dopant incorporation.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A quantum hamiltonian simulation benchmark

Hamiltonian simulation is one of the most important problems in quantum computation, and quantum singular value transformation (QSVT) is an efficient way to simulate a general class of Hamiltonians. However, the QSVT circuit typically involves multiple ancilla qubits and multi-qubit control gates. In order to simulate a certain class of n-qubit random Hamiltonians, we propose a drastically simplified quantum circuit that we refer to as the minimal QSVT circuit, which uses only one ancilla qubit and no multi-qubit controlled gates. We formulate a simple metric called the quantum unitary evolution score (QUES), which is a scalable quantum benchmark and can be verified without any need for classical computation. Under the globally depolarized noise model, we demonstrate that QUES is directly related to the circuit fidelity, and the potential classical hardness of an associated quantum circuit sampling problem. Under the same assumption, theoretical analysis suggests there exists an ‘optimal’ simulation time t opt ≈ 4.81, at which even a noisy quantum device may be sufficient to demonstrate the potential classical hardness.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Characterizing electrical demand and load diversity of low-power water and space heating appliances in US homes

Home renovation and remodeling projects can involve costly and time consuming electrical infrastructure upgrades at the household level. From the grid perspective they also lead to costly replacement of local infrastructure, such as transformers, and can add stress to the grid at peak times. The emergence of innovative, power-efficient household appliances offers a way to minimize these problems. These appliances are designed for lower power consumption, simplifying installation through standard plug-in connections, reducing the need for new electric circuits/panels/service, and minimizing the peak power demand for the home. Key examples include low-power heat pump water heaters (HPWHs) and cold climate window heat pumps that operate on standard 120V outlets. To assess the real-world impact of these solutions, we compiled and analyzed power metering data from several US field studies. This data provides insights into the effects on peak power demand of selecting lower-power appliances. Our analysis focuses on several key metrics, including the maximum power demand of individual appliances, their operational runtime, continuous operation and load diversity. While individual low-power 120V space and water heating appliances offer significant peak demand reductions compared to 240V heat pump or resistance alternatives, their longer runtimes might increase the likelihood of operation during whole-dwelling peak events, albeit at lower power levels.

Less, Brennan↗
DOE OSTI2022

Validation of wind resource and energy production simulations for small wind turbines in the United States

Abstract. Due to financial and temporal limitations, the small wind community relies upon simplified wind speed models and energy production simulation tools to assess site suitability and produce energy generation expectations. While efficient and user-friendly, these models and tools are subject to errors that have been insufficiently quantified at small wind turbine heights. This study leverages observations from meteorological towers and sodars across the United States to validate wind speed estimates from the Wind Integration National Dataset (WIND) Toolkit, the European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis v5 (ERA5), and the Modern-Era Retrospective Analysis for Research and Applications, version 2 (MERRA-2), revealing average biases within ±0.5 m s−1 at small wind hub heights. Observations from small wind turbines across the United States provide references for validating energy production estimates from the System Advisor Model (SAM), Wind Report, MyWindTurbine.com, and Global Wind Atlas 3 (GWA3), which are seen to overestimate actual annual capacity factors by 2.5, 4.2, 11.5, and 7.3 percentage points, respectively. In addition to quantifying the error metrics, this paper identifies sources of model and tool discrepancies, noting that interannual fluctuation in the wind resource, wind speed class, and loss assumptions produces more variability in estimates than different horizontal and vertical interpolation techniques. The results of this study provide small wind installers and owners with information about these challenges to consider when making performance estimates and thus possible adjustments accordingly. Looking to the future, recognizing these error metrics and sources of discrepancies provides model and tool researchers and developers with opportunities for product improvement that could positively impact small wind customer confidence and the ability to finance small wind projects.

17 WIND ENERGY↗
DOE OSTI2025

Cost analysis of alternative large-scale high-temperature solid oxide electrolysis hydrogen production facilities

We extend our past cost analysis of gigawatt-scale solid oxide electrolysis (SOE) facilities that produce high purity hydrogen gas from water by estimating construction and operating costs for three new alternative design cases: (1) offsite feed steam generation; (2) near-atmospheric pressure (NAP) stack; and (3) onsite electric boiler feed steam generation. Pressure effects on hydrogen electrode-(cathode-)supported SOE cell (SOEC) stack performance are estimated for the same assumed cell and stack construction and used to determine facility-wide stack capital costs for achieving a fixed H2 production at different pressures. Modular balance of plant (BOP) process equipment capital costs are estimated for each new alternative design case using our past equipment sizing, design, and cost data and scaling relationships. Furthermore, we update BOP equipment sizing and design for the NAP case using Aspen®. Vendor quotes for electric boilers are used to estimate costs for the electric boiler design case. Factory and onsite assembly and installation costs for SOEC stacks and BOP equipment are calculated using our past simplified first-principles approach. First-of-a-kind (FOAK) and N th -of-a-kind (NOAK) production maturity cost estimates are included for all results. The case with NAP stacks offers the lowest facility total capital cost (TCC, ~23% lower than base) while use of small electric boilers requires the highest TCC (~3% higher than base). H 2 production prices decrease from the base of ~$\$2.17$ /kgH 2 to ~$\$1.92$/kgH 2 for 1 GW e DC SIP facilities utilizing NAP stacks supplied by offsites steam situated in large modules and blocks for $\$0.030$/kWh e and $0.009/kWh t prices for electricity and thermal energy, respectively. We report all costs in 2021 US dollars.

Balance of plant (BOP) process equipment↗
DOE OSTI2021

Search for dark matter produced in association with a leptonically decaying Z boson in proton-proton collisions at $\sqrt{s} =$ 13 TeV

A search for dark matter particles is performed using events with a Z boson candidate and large missing transverse momentum. The analysis is based on proton–proton collision data at a center-of-mass energy of 13$\,\text {Te}\text {V}$, collected by the CMS experiment at the LHC in 2016–2018, corresponding to an integrated luminosity of 137$\,\text {fb}^{-1}$. The search uses the decay channels ${\mathrm{Z}} \rightarrow {\mathrm{e}} {\mathrm{e}} $ and ${\mathrm{Z}} \rightarrow {{\upmu }{}{}} {{\upmu }{}{}} $. No significant excess of events is observed over the background expected from the standard model. Limits are set on dark matter particle production in the context of simplified models with vector, axial-vector, scalar, and pseudoscalar mediators, as well as on a two-Higgs-doublet model with an additional pseudoscalar mediator. In addition, limits are provided for spin-dependent and spin-independent scattering cross sections and are compared to those from direct-detection experiments. The results are also interpreted in the context of models of invisible Higgs boson decays, unparticles, and large extra dimensions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Application of FARM to an IES scenario within the FORCE ecosystem

The FARM (Feasible Actuator Range Modifier) software module is a component of the RAVEN-based FORCE framework for analysis of Integrated Energy Systems (IES). FARM supports the HERON software module in the evaluation of the optimal dispatch by evaluating feasible set-points for the different IES unit components. Set-points are required to satisfy limits on both production variables (i.e., the variables to be optimized such as the electrical power, the hydrogen production rate, etc.) and process variables tied to the service life of equipment (e.g., steam flowrate, vessel pressure, turbine firing temperature, etc.). This problem is addressed by adopting a two-stage approach. First, the HERON power dispatcher determines set-points that meet the constraints on the former variables (e.g., power levels and power ramp rate limits). These constraints are called explicit constraints. Then, FARM adjusts these set-points to ensure the respect of the limits on the latter variables given the knowledge of the system physics acquired through machine learning algorithms. These constraints are called implicit constraints. From this standpoint, FARM constitutes a bridge between the HERON power dispatcher that adopts a simplified description of the IES unit (low-resolution physics) and the HYBRID high-fidelity models (high-resolution physics). The original version of the FARM software module (FARM-Alpha) was released by Argonne National Laboratory in January 2021. In the latest version of the code released in July 2022 (FARM-Delta), the Reference Governor (RG) algorithm was upgraded to a Multi-Input Multi-Output version from its original Single-Input Single-Output form. The RG algorithm acts to enforce constraints. With this improvement an IES unit is now treated as a single dynamic system from the standpoint of control. The crosstalk among components in an IES unit is now fully considered thereby ensuring a true optimization is obtained for those units that have multiple set-points. In this report, the capabilities of FARM-Delta operating within the FORCE ecosystem are demonstrated for an IES test case. The specific configuration of IES unit for this case was selected by the IES team with consultation from the Advanced Reactor IES Expert Group. A full TEA analysis that invoked HERON, HYBRID, FARM, and RAVEN was performed and serves to demonstrate how the latest modification to FARM algorithms (i.e., state variable selection, state-space matrices derivation, set-point verification) can shape setpoints that might otherwise compromise the health of equipment through accelerated wear and tear. In this specific test case, it was demonstrated that these algorithms ensure a more efficient utilization of steam resources to be shared by two different subsystems, namely Balance of Plant (BOP) and High-Temperature Steam Electrolysis (HTSE). Finally, some code improvements that can further enhance the user-friendliness are suggested.

97 MATHEMATICS AND COMPUTING↗