Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reactive Flow”

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 523 records · Page 29
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Quartz Dissolution Effects on Flow Channelization and Transport Behavior in Three‐Dimensional Fracture Networks

We perform a set of reactive transport simulations in three-dimensional fracture networks to characterize the impact of geochemical reactions on flow channelization. Flow channelization, a frequently observed phenomenon in porous and fractured subsurface rock formations, results from the spatially variable hydraulic resistance offered by a geological structure. In addition to geo-structural features such as network connectivity, geometry, and hydraulic resistance, geochemical reactions, for example, dissolution and precipitation, can dynamically inhibit or enhance flow channelization. These geochemical processes can change the fracture permeability leading to increased flow channelization, which are localized connected regions of high volumetric flow rates that are seemingly ubiquitous in the subsurface. In our simulations, fractures partially filled with quartz are gradually dissolved until quasi-steady state conditions are obtained. We compare the flow field's initial unreacted and final dissolved states in terms of flow and transport observations. We observe that the dissolved fracture networks provide less resistance to flow and exhibit increased flow channelization when compared to their unreacted counterparts. However, there is substantial variability in the magnitude of these changes which implies that the channelization strongly depends on the network structure. In turn, we identify the interplay between the particular network structure and the impact of geochemical dissolution on flow channelization. The presented results indicate that geological systems that have been weathering or reactive for longer times in older landscapes are likely to have increased flow channelization compared to their equivalent but younger counterparts, which implies a time dependence on flow channelization in fractured media.

Hyman, Jeffrey D.↗
DOE OSTI2018

Reactivity of heterogeneous surfaces: Modeling precursor–surface interaction using absorbing Markov chains

Heterogeneous reactive surfaces involving a spatial distribution of reactive sites are ubiquitous in thin film growth, from defect nucleation to step flow growth. In this work, the interaction of gas phase species with heterogeneous reactive surfaces was explored using a model combining absorbing Markov chains and local balance equations. Through this model, single particle statistics can be extracted of the adsorption-diffusion-desorption/reaction process even beyond the transport-limited regime, where the surface coverage of adsorbates, and hence adsorbate-adsorbate interaction, is not negligible, and local rates can be connected with macroscopic quantities, such as the sticking probability and desorption rates. In the transport-limited regime, a similarity law was identified for the sticking probability for both step flow growth and nucleation on randomly distributed reactive sites in which the sticking probability is a function of the product of the local diffusion and reaction probabilities. Lastly, the model can be applied to extract local reactivity maps providing the probability that a particle adsorbing on a given site reacts with the surface. This allows one to explore the development of capture zones and the impact of site-blocking and desorption on the local sticking probabilities.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Organizational Principles of Hyporheic Exchange Flow and Biogeochemical Cycling in River Networks Across Scales

Hyporheic zones increase freshwater ecosystem resilience to hydrological extremes and global environmental change. However, current conceptualizations of hyporheic exchange, residence time distributions, and the associated biogeochemical cycling in streambed sediments do not always accurately explain the hydrological and biogeochemical complexity observed in streams and rivers. Specifically, existing conceptual models insufficiently represent the coupled transport and reactivity along groundwater and surface water flow paths, the role of autochthonous organic matter in streambed biogeochemical functioning, and the feedbacks between surface-subsurface ecological processes, both within and across spatial and temporal scales. While simplified approaches to these issues are justifiable and necessary for transferability, the exclusion of important hyporheic processes from our conceptualizations can lead to erroneous conclusions and inadequate understanding and management of interconnected surface water and groundwater environments. This is particularly true at the landscape scale, where the organizational principles of spatio-temporal dynamics of hyporheic exchange flow (HEF) and biogeochemical processes remain largely uncharacterized. This article seeks to identify the most important drivers and controls of HEF and biogeochemical cycling based on a comprehensive synthesis of findings from a wide range of river systems. We use these observations to test current paradigms and conceptual models, discussing the interactions of local-to-regional hydrological, geomorphological, and ecological controls of hyporheic zone functioning. This improved conceptualization of the landscape organizational principles of drivers of HEF and biogeochemical processes from reach to catchment scales will inform future river research directions and watershed management strategies.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Characteristics of flow through randomly packed impermeable and permeable particles using pore resolved simulations

Pore resolved simulations are performed to study the mean flow characteristics in porous media formed by random distribution of impermeable and permeable particles. Here, permeability of the medium and reactive surface areas are analyzed for a range of solid fractions, packing patterns and flow rates. The results are compared with models presented in the literature and with simulations of flow through ordered packings. Permeability of a randomly packed porous medium is found to depend on the heterogeneous distribution of porosity which results from particle agglomeration and is typically ignored by simpler models. The reactive surface area is dependent on the flow field at pore scale in addition to geometry. Moreover, the wake formed downstream of the particles at moderately high flow rates also affects the reactive surface area. Similar analysis is also performed by assuming particles to be permeable to account for dual scale porosity relevant for various engineering applications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1991

Surface reaction and pore diffusion in flow-tube reactors

The interaction of gas diffusion with surface reaction in porous solids is discussed and applied specifically to heterogeneous rate measurements in flow-tube reactors. External diffusion to the outer surface of a reactive solid, internal diffusion within the pores, surface reaction, and laminar flow are considered. A procedure is developed to correct observed surface rate constants for the interaction of these processes. Measured surface areas and bulk densities are used to construct a semiempirical model for porous diffusion in vapor-formed HNO3-H2O ices which are used to simulate polar stratospheric cloud surfaces. The model is tested experimentally by varying the thickness of these ices from about 15 to 120 microns. The results are consistent with the model predictions and show that the HNO3-H2O ices used are highly porous, and the internal surface must be considered in calculating kinetic parameters from observed loss rates. The best fit of the data yields a tortuosity factor of 3.3 +/-1.1 for the ice substrates.

Keyser, Leon F.↗
DOE OSTI2017

Tutorial: Reactive high power impulse magnetron sputtering (R-HiPIMS)

High Power Impulse Magnetron Sputtering (HiPIMS) is a coating technology that combines magnetron sputtering with pulsed power concepts. Furthermore, by applying power in pulses of high amplitude and a relatively low duty cycle, large fractions of sputtered atoms and near-target gases are ionized. In contrast to conventional magnetron sputtering, HiPIMS is characterized by self-sputtering or repeated gas recycling for high and low sputter yield materials, respectively, and both for most intermediate materials. The dense plasma in front of the target has the dual function of sustaining the discharge and providing plasma-assistance to film growth, affecting the microstructure of growing films. Many technologically interesting thin films are compound films, which are composed of one or more metals and a reactive gas, most often oxygen or nitrogen. When reactive gas is added, non-trivial consequences arise for the system because the target may become “poisoned,” i.e., a compound layer forms on the target surface affecting the sputtering yield and the yield of secondary electron emission and thereby all other parameters. It is emphasized that the target state depends not only on the reactive gas' partial pressure (balanced via gas flow and pumping) but also on the ion flux to the target, which can be controlled by pulse parameters. This is a critical technological opportunity for reactive HiPIMS (R-HiPIMS). The scope of this tutorial is focused on plasma processes and mechanisms of operation and only briefly touches upon film properties. It introduces R-HiPIMS in a systematic, step-by-step approach by covering sputtering, magnetron sputtering, reactive magnetron sputtering, pulsed reactive magnetron sputtering, HiPIMS, and finally R-HiPIMS. The tutorial is concluded by considering variations of R-HiPIMS known as modulated pulsed power magnetron sputtering and deep-oscillation magnetron sputtering and combinations of R-HiPIMS with superimposed dc magnetron sputtering.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Impact of lubricant additives on the physicochemical properties and activity of three-way catalysts

As alternative lubricant anti-wear additives are sought to reduce friction and improve overall fuel economy, it is important that these additives are also compatible with current emissions control catalysts. In the present work, an oil-miscible phosphorous-containing ionic liquid (IL), trihexyltetradecylphosphonium bis(2-ethylhexyl) phosphate ([P 66614 ][DEHP]), is evaluated for its impact on three-way catalysts (TWC) and benchmarked against the industry standard zinc-dialkyl-dithio-phosphate (ZDDP). The TWCs are aged in different scenarios: neat gasoline (no-additive, or NA), gasoline+ZDDP, and gasoline+IL. The aged samples, along with the as received TWC, are characterized through various analytical techniques including catalyst reactivity evaluation in a bench-flow reactor. The temperatures of 50% conversion (T50) for the ZDDP-aged TWCs increased by 30, 24, and 25 °C for NO, CO, and C3H6, respectively, compared to the no-additive case. Although the IL-aged TWC also increased in T50 for CO and C3H6, it was notably less than ZDDP, 7 and 9 °C, respectively. Additionally, the IL-aged samples had higher water-gas-shift reactivity and oxygen storage capacity than the ZDDP-aged TWC. Characterization of the aged samples indicated the predominant presence of CePO 4 in the ZDDP-aged TWC aged by ZDDP, while its formation was retarded in the case of IL where higher levels of AlPO 4 is observed. Furthermore, the results in this work indicate that the phosphonium-phosphate IL potentially has less adverse impact on TWC than ZDDP.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2019

Impact of Lubricant Oil Additives on the Performance of Pd-Based Three-Way Catalysts

As alternative lubricant anti-wear additives are sought to reduce friction and improve overall fuel economy, it is important that these additives are also compatible with current emissions control catalysts. In the present work, a second-generation oil-miscible phosphorous-containing ionic liquid (IL), is evaluated for its impact on Pd-based three-way catalyst (TWC) reactivity and benchmarked against the industry standard zinc dialkyl-dithio-phosphate (ZDDP). The TWCs are exposed to the lubricant additives in an engine bench under four different scenarios: base case with no additive (NA), ZDDP-only, IL-only, and IL + ZDDP. The engine-aged TWC samples, along with the as-received TWC, are characterized through various analytical techniques including catalyst reactivity evaluation in a bench-flow reactor. The temperature of 50% conversion (T 50 ) for the ZDDP-aged TWC increases by 11, 21, and 36 °C for CO, C 3 H 6 , and C 3 H 8 , respectively, as compared with the no-additive case. Similarly, the T50 for IL-only and IL + ZDDP-aged TWCs also increase as compared with the no-additive case. Even though the water-gas-shift (WGS) reactivity is similar for all engine-aged samples, the IL-aged TWC had higher oxygen storage capacity than the ZDDP-aged TWC. EPMA analysis reveals penetration of phosphorus deep into the washcoat for all engine-aged TWCs. Results from XRD indicate the presence of CePO 4 and AlPO 4 on the washcoat of IL, ZDDP, and IL + ZDDP-aged TWC samples but not on the fresh and NA TWC samples. Additionally, ICP-OES results show a large amount of phosphorus in the washcoat of IL, ZDDP and IL + ZDDP-aged TWC samples and a lesser amount in NA TWC samples.

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

A High-Order Discontinuous Galerkin Spectral Element Method for Compressible Reacting Flows

High-order methods have recently been shown to be an effective tool for high-fidelity flow computations like direct numerical simulations and large eddy simulations due to their strong balance between accuracy and computational cost. In this work, a high-order discontinuous Galerkin spectral element method (DGSEM) is developed to solve the chemically reactive Euler equations encountered in high-speed combustion. To handle the disparate length and time scales associated with these equations, we develop a novel method which combines the spectral accuracy of the SEM with the flexibility of DG approach. Thus, the framework is well suited to capture turbulence in smooth regions of the flow, while maintaining numerical stability in the presence of shocks. The numerical method is implemented within the spectral element solver Nek5000. Validation cases are conducted for both non-reactive and reactive discontinuous flows to demonstrate the solver capability. In particular, canonical one-dimensional and two-dimensional detonation simulations are performed and the high-order numerical results are validated against available literature data.

computational fluid dynamics (CFD)↗
DOE OSTI2013

Realistic multisite lattice-gas modeling and KMC simulation of catalytic surface reactions: Kinetics and multiscale spatial behavior for CO-oxidation on metal (100) surfaces

A realistic molecular-level description of catalytic reactions on single-crystal metal surfaces can be provided by stochastic multisite lattice-gas (msLG) models. This approach has general applicability, although in this report, we will focus on the example of CO-oxidation on the unreconstructed fcc metal (100) or M(100) surfaces of common catalyst metals M = Pd, Rh, Pt and Ir (i.e., avoiding regimes where Pt and Ir reconstruct). These models can capture the thermodynamics and kinetics of adsorbed layers for the individual reactants species, such as CO/M(100) and O/M(100), as well as the interaction and reaction between different reactant species in mixed adlayers, such as (CO + O)/M(100). The msLG models allow population of any of hollow, bridge, and top sites. This enables a more flexible and realistic description of adsorption and adlayer ordering, as well as of reaction configurations and configuration-dependent barriers. Adspecies adsorption and interaction energies, as well as barriers for various processes, constitute key model input. The choice of these energies is guided by experimental observations, as well as by extensive Density Functional Theory analysis. Model behavior is assessed via Kinetic Monte Carlo (KMC) simulation. We also address the simulation challenges and theoretical ramifications associated with very rapid diffusion and local equilibration of reactant adspecies such as CO. These msLG models are applied to describe adsorption, ordering, and temperature programmed desorption (TPD) for individual CO/M(100) and O/M(100) reactant adlayers. In addition, they are also applied to predict mixed (CO + O)/M(100) adlayer structure on the nanoscale, the complete bifurcation diagram for reactive steady-states under continuous flow conditions, temperature programmed reaction (TPR) spectra, and titration reactions for the CO-oxidation reaction. Extensive and reasonably successful comparison of model predictions is made with experimental data. Furthermore, we discuss the possible transition from traditional mean-field-type bistability and reaction kinetics for lower-pressure to multistability and enhanced fluctuation effects for moderate- or higher-pressure. Behavior in the latter regime reflects a stronger influence of adspecies interactions and also lower diffusivity in the higher-coverage mixed adlayer. We also analyze mesoscale spatiotemporal behavior including the propagation of reaction diffusion fronts between bistable reactive and inactive states, and associated nucleation-mediated transitions between these states. This behavior is controlled by complex surface mass transport processes, specifically chemical diffusion in mixed reactant adlayers for which we provide a precise theoretical formulation. The msLG models together with an appropriate treatment of chemical diffusivity enable equation-free heterogeneous coupled lattice-gas (HCLG) simulations of spatiotemporal behavior. In addition, msLG + HCLG modeling can describe coverage variations across polycrystalline catalysts surfaces, pressure variations across catalyst surfaces in microreactors, and could be incorporated into a multiphysics framework to describe mass and heat transfer limitations for high-pressure catalysis. (C) 2013 Elsevier Ltd. All rights reserved.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Experimental Validation of an Explicit Power-Flow Primary Control in Microgrids

The existing approaches to control electrical grids combine frequency and voltage controls at different time-scales. When applied in microgrids with stochastic distributed generation, grid quality of service problems may occur, such as under- or overvoltages as well as congestion of lines and transformers. The COMMELEC framework proposes to solve this compelling issue by performing explicit control of power flows with two novel strategies: 1) a common abstract model is used by resources to advertise their state in real time to a grid agent; and 2) subsystems can be aggregated into virtual devices that hide their internal complexity in order to ensure scalability. While the framework has already been published in the literature, in this paper, we present the first experimental validation of a practicable explicit power-flow primary control applied in a real-scale test-bed microgrid. We demonstrate how an explicit power-flow control solves the active and reactive power sharing problem in real time, easily allowing the microgrid to be dispatchable in real time (i.e., it is able to participate in energy markets) and capable of providing frequency support, while always maintaining quality of service.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2021

Lyapunov stability of smart inverters using linearized distflow approximation

Fast-acting smart inverters that utilize preset operating conditions to determine real and reactive power injection/consumption can create voltage instabilities (over-voltage, voltage oscillations and more) in an electrical distribution network if set-points are not properly configured. In this work, linear distribution power flow equations and droop-based Volt–Var and Volt–Watt control curves are used to analytically derive a stability criterion using Lyapunov analysis that includes the network operating condition. The methodology is generally applicable for control curves that can be represented as Lipschitz functions. The derived Lipschitz constants account for smart inverter hardware limitations for reactive power generation. A local policy is derived from the stability criterion that allows inverters to adapt their control curves by monitoring only local voltage, thus avoiding centralized control or information sharing with other inverters. The criterion is independent of the internal time-delays of smart inverters. Simulation results for inverters with and without the proposed stabilization technique demonstrate how smart inverters can mitigate voltage oscillations locally and mitigate real and reactive power flow disturbances at the substation under multiple scenarios. The study concludes with illustrations of how the control policy can dampen oscillations caused by solar intermittency and cyberattacks.

42 ENGINEERING↗
DOE OSTI2019

Effect of Oxygen Flow Rate on Structural, Electrical and Optical Properties of Zinc Aluminum Oxide Thin Films Deposited by DC Magnetron Sputtering

Zinc Aluminum Oxide thin films were deposited on glass substrates by reactive DC magnetron sputtering method by varying oxygen flow rates from 1 to 4 sccm. Glancing angle X-ray diffraction patterns ofzinc aluminum oxide thin films exhibits (0 0 2) peak with c-plane preferentially oriented parallel to the substrate. The surface morphology and elemental analysis of the films was observed by field emission scanning electron microscopy attached with energy dispersive X-ray analysis spectroscopy. An average optical transmittance of 83–90% is obtained for the films deposited at various oxygen flow rates. The optical band gap of the films increases from 3.41 to 3.53 eV with the increase of oxygen flow rate due to Burstein–Moss effect. The optical dispersion parameters such as dispersion energy (E{sub d}), oscillator energy (E{sub o}) and static refractive index (n{sub o}) were determined using the Wemple-DiDomenico (W-D) single oscillator model. The nonlinear optical parameters such as optical susceptibility (χ{sup (1)}), third order nonlinear optical susceptibility (χ{sup (3)}) and nonlinear refractive index (n{sub 2}) were also determined.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

A User Guide to PARET/ANL (V.7.6) (Rev.1)

PARET was originally created in 1969 at what is now Idaho National Laboratory (INL), to analyze reactivity insertion events in research and test reactor cores cooled by light or heavy water, with fuel composed of either plates or pins. The use of PARET is also appropriate for fuel assemblies with curved fuel plates when their radii of curvatures are large with respect to the fuel plate thickness. The PARET/ANL version of the code has been developed at Argonne National Laboratory (ANL) under the sponsorship of the U.S. Department of Energy/NNSA, and has been used by the Reactor Conversion Program to determine the expected transient behavior of a large number of reactors. PARET/ANL models the various fueled regions of a reactor core as channels. Each of these channels consists of a single flat fuel plate/pin (including cladding and, optionally, a gap) with water coolant on each side. In slab geometry the coolant channels for a given fuel plate are of identical dimensions (mirror symmetry), but they can be of different thickness in each channel. There can be many channels, but each channel is independent and coupled only through reactivity feedback effects to the whole core. The time-dependent differential equations that represent the system are replaced by an equivalent set of finite-difference equations in space and time, which are integrated numerically. PARET/ANL uses fundamentally the same numerical scheme as RELAP5 for the time-integration of the point-kinetics equations. The one-dimensional thermal-hydraulic model includes temperature-dependent thermal properties of the solid materials, such as heat capacity and thermal conductivity, as well as the transient heat production and heat transfer from the fuel meat to the coolant. Temperature- and pressure-dependent thermal properties of the coolant such as enthalpy, density, thermal conductivity, and viscosity are also used in determining parameters such as friction factors and heat transfer coefficients. The code first determines the steady-state solution for the initial state. Then the solution of the transient is obtained by integration in time and space. Multiple heat transfer, DNB and flow instability correlations are available. The code was originally developed to model reactors cooled by an open loop, which was adequate for rapid transients in pool-type cores. An external loop model appropriate for Miniature Neutron Source Reactors (MNSR’s) was also added to PARET/ANL to model natural circulation within the vessel, heat transfer from the vessel to pool and heat loss by evaporation from the pool. PARET/ANL also contains models for decay heat after shutdown, control rod reactivity versus time or position, time-dependent pump flow, and loss-of-flow event with flow reversal as well as logic for trips on period, power, and flow. Feedback reactivity effects from coolant density changes and temperature changes are represented by tables. Feedback reactivity from fuel heat-up (Doppler Effect) is represented by a four-term polynomial in powers of fuel temperature. Photo-neutrons produced in beryllium or in heavy water may be included in the point-kinetics equations by using additional delayed neutron groups.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2019

A User Guide to PARET/ANL

PARET was originally created in 1969 at what is now Idaho National Laboratory (INL), to analyze reactivity insertion events in research and test reactor cores cooled by light or heavy water, with fuel composed of either plates or pins. The use of PARET is also appropriate for fuel assemblies with curved fuel plates when their radii of curvatures are large with respect to the fuel plate thickness. The PARET/ANL version of the code has been developed at Argonne National Laboratory (ANL) under the sponsorship of the U.S. Department of Energy/NNSA since the inception of the Reactor Conversion Program. Since Reduced Enrichment for Research and Test Reactors (RERTR) began in 1978, PARET/ANL has been benchmarked to experimental data including SPERT testing, and used to determine the expected transient behavior of a large number of reactors regardless of enrichment both inside and outside the Reactor Conversion Program. This document provides the pertinent information for the use of PARET/ANL Version 7.6. PARET/ANL models the various fueled regions of a reactor core as channels. Each of these channels consists of a single flat fuel plate/pin (including cladding and, optionally, a gap) with water coolant on each side. In slab geometry the coolant channels for a given fuel plate are of identical dimensions (mirror symmetry), but they can be of different thickness in each channel. There can be many channels, but each channel is independent and coupled only through reactivity feedback effects to the whole core. The time-dependent differential equations that represent the system are replaced by an equivalent set of finite-difference equations in space and time, which are integrated numerically. PARET/ANL uses fundamentally the same numerical scheme as RELAP5 for the time-integration of the point-kinetics equations. The one-dimensional thermal-hydraulic model includes temperature-dependent thermal properties of the solid materials, such as heat capacity and thermal conductivity, as well as the transient heat production and heat transfer from the fuel meat to the coolant. Temperature- and pressure-dependent thermal properties of the coolant such as enthalpy, density, thermal conductivity, and viscosity are also used in determining parameters such as friction factors and heat transfer coefficients. The code first determines the steady-state solution for the initial state. Then the solution of the transient is obtained by integration in time and space. Multiple heat transfer, DNB and flow instability correlations are available. The code was originally developed to model reactors cooled by an open loop, which was adequate for rapid transients in pool-type cores. An external loop model appropriate for Miniature Neutron Source Reactors (MNSR’s) was also added to PARET/ANL to model natural circulation within the vessel, heat transfer from the vessel to pool and heat loss by evaporation from the pool. PARET/ANL also contains models for decay heat after shutdown, control rod reactivity versus time or position, time-dependent pump flow, and loss-of-flow event with flow reversal as well as logic for trips on period, power, and flow. Feedback reactivity effects from coolant density changes and temperature changes are represented by tables. Feedback reactivity from fuel heat-up (Doppler Effect) is represented by a four-term polynomial in powers of fuel temperature. Photo-neutrons produced in beryllium or in heavy water may be included in the point-kinetics equations by using additional delayed neutron groups.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

A User Guide to PARET/ANL

PARET was originally created in 1969 at what is now Idaho National Laboratory (INL), to analyze reactivity insertion events in research and test reactor cores cooled by light or heavy water, with fuel composed of either plates or pins. The use of PARET is also appropriate for fuel assemblies with curved fuel plates when their radii of curvatures are large with respect to the fuel plate thickness. The PARET/ANL version of the code has been developed at Argonne National Laboratory (ANL) under the sponsorship of the U.S. Department of Energy/NNSA since the inception of the Reactor Conversion Program. Since Reduced Enrichment for Research and Test Reactors (RERTR) began in 1978, PARET/ANL has been benchmarked to experimental data including SPERT testing, and used to determine the expected transient behavior of a large number of reactors regardless of enrichment both inside and outside the Reactor Conversion Program. This document provides the pertinent information for the use of PARET/ANL Version 7.6. PARET/ANL models the various fueled regions of a reactor core as channels. Each of these channels consists of a single flat fuel plate/pin (including cladding and, optionally, a gap) with water coolant on each side. In slab geometry the coolant channels for a given fuel plate are of identical dimensions (mirror symmetry), but they can be of different thickness in each channel. There can be many channels, but each channel is independent and coupled only through reactivity feedback effects to the whole core. The time-dependent differential equations that represent the system are replaced by an equivalent set of finite-difference equations in space and time, which are integrated numerically. PARET/ANL uses fundamentally the same numerical scheme as RELAP5 for the time-integration of the point-kinetics equations. The one-dimensional thermal-hydraulic model includes temperature-dependent thermal properties of the solid materials, such as heat capacity and thermal conductivity, as well as the transient heat production and heat transfer from the fuel meat to the coolant. Temperature- and pressure-dependent thermal properties of the coolant such as enthalpy, density, thermal conductivity, and viscosity are also used in determining parameters such as friction factors and heat transfer coefficients. The code first determines the steady-state solution for the initial state. Then the solution of the transient is obtained by integration in time and space. Multiple heat transfer, DNB and flow instability correlations are available. The code was originally developed to model reactors cooled by an open loop, which was adequate for rapid transients in pool-type cores. An external loop model appropriate for Miniature Neutron Source Reactors (MNSR’s) was also added to PARET/ANL to model natural circulation within the vessel, heat transfer from the vessel to pool and heat loss by evaporation from the pool. PARET/ANL also contains models for decay heat after shutdown, control rod reactivity versus time or position, time-dependent pump flow, and loss-of-flow event with flow reversal as well as logic for trips on period, power, and flow. Feedback reactivity effects from coolant density changes and temperature changes are represented by tables. Feedback reactivity from fuel heat-up (Doppler Effect) is represented by a four-term polynomial in powers of fuel temperature. Photo-neutrons produced in beryllium or in heavy water may be included in the point-kinetics equations by using additional delayed neutron groups.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2023

Constraints on OPF Surrogates for Learning Stable Local Volt/Var Controllers

We consider the problem of learning local Volt/Var controllers in distribution grids (DGs). Our approach starts from learning separable surrogates that take both local voltages and reactive powers as arguments and predict the reactive power setpoints that approximate optimal power flow (OPF) solutions. We propose an incremental control algorithm and identify two different sets of slope conditions on the local surrogates such that the network is collectively steered toward desired configurations asymptotically. Our results reveal the trade-offs between each set of conditions, with coupled voltage-power slope constraints allowing an arbitrary shape of surrogate functions but risking limitations on exploiting generation capabilities, and reactive power slope constraints taking full advantage of generation capabilities but constraining the shape of surrogate functions. AC power flow simulations on the IEEE 37-bus feeder illustrate their guaranteed stability properties and respective advantages in two DG scenarios.

asymptotic stability↗