Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Temperature reactivity coefficient”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
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 OSTI2026

Ab Initio Modeling of Aqueous Methanol Mixtures at DFT-SCAN Level Using Machine Learning Interatomic Potentials

Abstract Methanol–water mixtures find use in many applications, particularly catalytic energy conversion processes. Their importance has motivated numerous computational studies, most of which employed molecular dynamics based on classical force fields. These enable simulations of large systems on long time scales but do not reliably describe reactive dynamics involving bond breaking and bond formation. In contrast, ab initio molecular dynamics (AIMD) based on density functional theory (DFT) is generally more reliable for such applications but has a high computational cost, which discourages systematic studies of alcohol-water mixtures. To remedy this, we trained a machine learning interatomic potential capable of probing the properties of aqueous methanol mixtures at the DFT level using the SCAN functional. Our results show that SCAN qualitatively reproduces multiple key experimental features arising from the amphiphilic nature of methanol, including density, diffusion coefficients, X-ray structure factors, and Kirkwood–Buff integrals. We also find that structural correlations between water molecules are somewhat overestimated, leading to a stronger preferential association than that predicted by experiments. However, increasing the temperature by 30 K mitigates this effect and also recovers the correct mobilities of both methanol and water. These results indicate that SCAN provides an accurate description of methanol–water mixtures, making it a reliable choice for investigating the reactive dynamics in such systems.

Park, Sanghyun J. [Princeton University , , , ,]↗
DOE OSTI2022

Chemical kinetic interactions and sensitivity analyses for 2-ethylhexyl nitrate-doped PRF91 using a reduced mechanism

A numerical and experimental investigation about the chemical kinetic interactions between 2-ethylhexylnitrate (EHN) and PRF91 was performed in this study. Rapid compression machine experiments were conducted to investigate the effect of EHN on the autoignition reactivity of the fuel, and a reduced chemical kinetic mechanism was developed including an EHN sub-model. Experiments showed that the ignition delay decreases as the fuel is doped with EHN, but the effect of the doping level of EHN on the ignition is highly non-linear. Moreover, experiments showed that the EHN effectiveness is lowest during the transition between the low-temperature regime and the negative temperature coefficient (NTC) regime, and it rapidly increases as the temperature increases. Both detailed and (developed) reduced mechanisms were validated against the experimental results, allowing a more in-depth EHN-fuel chemistry study. Additionally, ignition delay sensitivity and EHN effectiveness sensitivity analyses were performed with the reduced mechanism to identify the reactions that control the interaction between EHN and the fuel. As the result, EHN thermal decomposition is only relevant for very low temperatures. The chemistry of EHN-doped fuel is more sensitive to intermediate temperature reactions than that of straight fuel, especially at lower temperatures, due to the effect of EHN on the NTC behavior of the fuel. Finally, the chemistry of EHN-doped fuel is very sensitive to NO 2 -to-NO reactions, especially at high temperatures, because these reactions transform the NO 2 generated by EHN into NO, which is a very effective fuel reactivity enhancer.

2-ethylhexyl nitrate↗
DOE OSTI2022

A multiphysics model of the versatile test reactor based on the MOOSE framework

The traditional modeling approach for sodium fast reactor cores relies on separate physics models, where the fuel performance, thermal–hydraulics, and neutronics calculations required to predict the core physics characteristics for nominal conditions are decoupled by relying on user-imposed boundary conditions. Here, this paper aims at evaluating the impact of multiphysics simulations for predicting the core characteristics of the Versatile Test Reactor, which is being designed as a 300-MWt sodium-cooled fast reactor. The purpose of the Versatile Test Reactor is to accelerate the testing of advanced nuclear materials in the United States. The proposed multiphysics model relies on the Griffin reactor physics code, the SAM thermal–hydraulic system code, the BISON fuel performance code, as well as generic Multiphysics Object-Oriented Simulation Environment capabilities implemented in the open-source tensor mechanics module. For k eff calculations, the introduction of a tight coupling between the neutronics, thermo-mechanical and thermal–hydraulics models induces a change of around 543 pcm in the eigenvalue, compared to the traditional standalone neutronics calculation where approximate temperature profiles are used. The multiphysics model is then employed for quantifying the impact of the thermal conductivity uncertainties on some of the key figures of merit, such as the fuel centerline temperature, assembly powers, and keff for nominal core conditions. As anticipated, uncertainties on fuel thermal conductivity mostly impact the fuel centerline temperature, and to a lesser extend the k eff .

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Electrodeposition of Magnonic V(tetracyanoethylene) 2 Thin Films

Molecule-based magnetic materials have been identified as promising candidates for application in magnonic technologies, owing not only to their solution processability but also because they can exhibit narrow ferromagnetic resonance (FMR) linewidths and low Gilbert damping coefficients—crucial prerequisites for the transmission of coherent magnons over macroscopic distances. In particular, V(TCNE) 2 , a compound with a threedimensional network structure composed of vanadium(II) centers linked by tetracyanoethylene (TCNE •− ) radical anions, displays magnonic properties comparable to yttrium iron garnet, the quintessential magnonic material in the field. However, existing solution and chemical vapor deposition methods for synthesizing V(TCNE) 2 require the use of highly reactive zero-valent molecular vanadium precursors, stymying research on this important material. Herein, we report a facile electrochemical method for the deposition of thin films of V(TCNE) 2 using readily obtainable and stable divalent vanadium precursors and TCNE •− anions generated by electrochemical reduction. Magnetization measurements reveal that the films exhibit ferrimagnetic ordering above room temperature, consistent with V(TCNE) 2 films synthesized via other methods. Moreover, the electrodeposited films exhibit narrow FMR linewidths as low as 17.5 G and a low Gilbert damping coefficient of 1.1 × 10 −3 , values that are on par with some currently integrated metallic magnonic materials. More generally, these results demonstrate that electrodeposition can provide a straightforward means of generating highperformance magnonic materials using readily available molecular precursors.

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

Pulse radiolysis and transient absorption spectra of aqueous solutions of sodium sulfamate

Chemical kinetics for the reactions of sulfamate ions (NH 2 SO 3 − ) with the primary products of water and nitric acid radiolysis were measured in aqueous solutions at ambient temperature. Using time-resolved electron pulse radiolysis techniques with a custom multichannel detection system, we examined the reactivity of NH 2 SO 3 − with the hydroxyl radical ( • OH), hydrogen atom (H • ), and nitrate radical (NO 3 • ). The sulfamate ion was found to react with • OH and H • with second-order rate coefficients of k • OH = (5.60 ± 0.04) × 10 6 M −1 s −1 and k H • = (7.96 ± 0.10) × 10 6 M −1 s −1 , respectively, and with NO 3 • with a rate coefficient of k NO 3 • = (1.67 ± 0.06) × 10 7 M −1 s −1 . The reactions of NH 2 SO 3 − with • OH and H • resulted in the formation of two transient radical species, one with maximum absorbance at 300 nm and a second with maxima at both 300 nm and 600 nm. These spectra are tentatively assigned to • NH 2 SO 3 and • NHSO 3 − , respectively. By measuring the absorbance of these radicals as a function of pH, the radical pK a was determined to be 9.5 ± 0.1. Overall, this work has implications for the longevity and performance of ferrous sulfamate, Fe(NH 2 SO 3 ) 2 , as a plutonium reductant in the reprocessing of used nuclear fuel.

Conrad, Jacy K. [Idaho National Laboratory (INL), ↗
DOE OSTI2020

Solvent selection for biphasic extraction of 5-hydroxymethylfurfural via multiscale modeling and experiments

We introduce a comprehensive conceptual framework for selecting solvents for reactive extraction in biphasic organic-water systems and demonstrate it for the separation of HMF (5-hydroxylmethylfurfural), a platform chemical produced in the acid-catalyzed dehydration of hexoses. We first perform in silico screening of ~2500 solvents, from the ADFCRS-2018 database using the ADF COSMO-RS implementation, and classification, based on the solvent partition coefficient. We then determine experimentally the partition coefficients for HMF, fructose, and products of HMF rehydration (levulinic acid (LA), and formic acid (FA)), the mutual water-organic solvent solubilities, and the separation factors in >50 select solvents spanning multiple homologous series at room temperature and a typical reaction temperature with in situ sampling. We find that COSMO-RS is excellent for screening purposes (typical error in most cases within a factor of ~2). Increased temperatures lead to significant reduction in partitioning, and room temperature measurements are clearly inadequate for solvent selection. Upon down selecting classes of solvents based on separation performance, we perform experimental thermal stability and reaction compatibility studies of a small set of solvents at relevant reactive-extraction temperatures. We discover that many substituted phenols exhibit an order-of-magnitude increase in partitioning compared to conventional solvents due chiefly to hydrogen bond interactions and show the necessary stability but retain a significant fraction of water and LA, factors that need to be considered in technoeconomic analysis. In contrast, anilines, aldehydes, and acids are good to excellent regarding separation but incompatible with this specific reaction medium. Furthermore, this multifaceted framework can be extended to other biomass-derived products and processes.

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

TOUGHREACT V4.12

TOUGHREACT V4.12 is a 3-D reactive transport code that handles to non-isothermal, multiphase/multicomponent fluid flow, heat transport, aqueous and gaseous species advection-diffusion, and equilibrium/kinetic water-gas-rock-biological reactions. TOUGHREACT V4.12 is based on TOUGHREACT V3.3, with several new features including: an option to compute activity coefficients using the Pitzer ion interaction model; implementation of ECO2N V2.0 for temperatures up to 300 C; options to simulate forward and reverse osmosis, including 2D Poiseuille flow; addition of EOS8 and revisions for EOS5 and EWASG with H2(g); addition of gas species decay, various new output formats and time step control options; many bug fixes.

Spycher, Nicolas↗
DOE OSTI2023

Impact of carbamate formation on the surface tension of epoxy‐amine curing systems

Abstract Epoxy‐amine reactive systems, including their surface tension, have been studied due to their use and importance in composites and adhesives in laboratory and industrial applications. Understanding how the system's surface tension changes during curing could help target applications for tailored composites. Using the Wilhelmy plate method, surface tensions for epoxy‐curing agent systems were investigated under different environments. While the epoxy exhibited an expected negative temperature coefficient of −0.096 mN/m/°C, the curing agent exhibited a positive temperature coefficient of 0.073 mN/m/°C, which is due to the curing agent reacting with atmospheric CO 2 to form carbamate through a zwitterion pathway and functional groups accumulating at the surface of the reacting systems. Carbamate formation results in a surface tension below that of the individual constituents in the epoxy‐amine curing agent system. Fourier transform infrared (FTIR) spectroscopy was used to verify in situ formation of carbamate, and by reducing the system's temperature, it was demonstrated that the system surface tension profiles displayed micellar surfactant‐like characteristics. FTIR was used to explore the difference in composition between the surface and the bulk, in particular carbamate formation at the surface, resulting in a lowering of the surface tension of reacting systems.

Shenk, Timothy M.↗
DOE OSTI2024

Kinetics of Radiation-induced Cr( II ) and Cr( III ) Redox Chemistry in Molten LiCl-KCl Eutectic

Chromium (Cr) is a frequent constituent of the metal alloys proposed for molten salt nuclear reactor (MSR) applications, and is typically the least noble metal ion present. Consequently, chromium is preferentially corroded into molten salt solutions. Here, the redox poise and redox cycling of chromium ions in the salt can greatly influence its corrosivity towards structural alloys, ultimately impacting the longevity of MSR systems. Radiation-induced chemistry is expected to play a significant role in determining the chromium oxidation state distribution during MSR operations. In the present research, electron pulse radiolysis techniques were employed to characterize the reactivity of Cr(II) and Cr(III) ions with primary radiolysis products in molten lithium chloride–potassium chloride (LiCl–KCl) eutectic over a temperature range of 400–600 °C. Both chromium oxidation states were found to rapidly react with the primary products of molten chloride salt radiolysis, i.e., the solvated electron (e S - ) and the dichlorine radical anion (Cl 2 ˙ - ). For reactions with the e S - , second-order rate coefficients (k) of k = (4.1 ± 0.2) and (6.1 ± 0.3) × 10 10 M -1 s -1 at 400 °C for Cr(II) and Cr(III), respectively, were determined. Temperature-dependent measurements allowed for the derivation of activation parameters for electron capture by Cr(II) and Cr(III). Both chromium ions also react with Cl 2 ˙ - , k = (7.2 ± 0.3) and (1.4 ± 0.1) × 10 9 M -1 s -1 at 400 °C for Cr(II) and Cr(III), respectively.

37 - INORGANIC, ORGANIC, PHYSICAL AND ANALYTICAL C↗
DOE OSTI2021

Ignition delay times of NH 3 /DME blends at high pressure and low DME fraction: RCM experiments and simulations

Here, autoignition delay times of ammonia/dimethyl ether (NH 3 /DME) mixtures were measured in a rapid compression machine with DME fractions of 0, 2 and 5 and 100% in the fuel. The measurements were performed at equivalence ratios $\varphi$= 0.5, 1.0 and 2.0 and pressures in the range 10–70 bar; depending on the fuel composition, the temperatures after compression varied from 610 K to 1180 K. Admixture of DME is seen to have a dramatic effect on the ignition delay time, effectively shifting the curves of ignition delay vs. temperature to lower temperatures, up to ~250 K compared to pure ammonia. Two-stage ignition is observed at $\varphi$= 1.0 and 2.0 with 2% and 5% DME in the fuel, despite the pressure being higher than that at which pure DME shows two-stage ignition. At $\varphi$= 0.5, a reproducible pre-ignition pressure rise is observed for both DME fractions, which is not observed in the pure fuel components. A novel NH 3 /DME mechanism was developed, including modifications in the NH 3 subset and addition of the NH 2 + CH 3 OCH 3 reaction, with rate coefficients calculated from ab initio theory. Simulations faithfully reproduce the observed pre- ignition pressure rise. While the mechanism also exhibits two-stage ignition for NH 3 /DME mixtures, both qualitative and quantitative improvement is recommended. The overall ignition delay times for ammonia/DME mixtures are predicted well, generally being within 50% of the experimental values, although reduced performance is observed for pure ammonia at $\varphi$= 2.0. Simulating the ignition process, we observe that the DME is oxidized much more rapidly than ammonia. Analysis of the mechanism indicates that this ‘early DME oxidation’ generates reactive species that initiate the oxidation of ammonia, which in turn begins heat release that raises the temperature and accelerates the oxidation process towards ignition. The reaction path analysis shows that the low-temperature chain-branching reactions of DME are important in the early oxidation of the fuel, while the sensitivity analysis indicates that several reactions in the oxidation of DME, including cross reactions between DME and NH 3 species presented here, are critical to ignition, even at fractions of 2% DME in the fuel.

10 SYNTHETIC FUELS↗
DOE OSTI2020

Predicting the effects of thermally-induced density gradients on the hydrodynamic behavior of PBX 9502

High explosives are often exposed to thermal and mechanical stimuli within normal and abnormal environments. These stimuli can introduce density gradients in the explosive charge but the effect of such gradients is not well-understood. The goal of this work is to develop a comprehensive methodology to model the evolution of density gradients under various thermal conditions and then to predict shock initiation and detonation behavior of the gradated sample. PBX 9502 was chosen for these scoping studies because it has a large coefficient of thermal expansion and because the available shock initiation data indicates that initial temperature strongly affects the shock-to-detonation transition (SDT). The multiphysics code COMSOL was first used to model induced density gradients in a slab of PBX 9502 with asymmetric heating. Then, the gradated sample was constructed in the Pagosa hydrocode to simulate steel flyer impact scenarios, with SDT and detonation propagation following a calibrated Scaled Uniform Reactive Front (SURF) burn model. Furthermore, the gradated material exhibited very different initiation and corner turning behavior compared to homogeneous materials, and the orientation of the gradation was also found to have some effect. The methodology used here can estimate thermomechanical response of the material for non-uniform thermal environments and compute the resulting detonation properties of the PBX.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Exploration of Fuel Property Impacts on the Combustion of Late Post Injections Using Binary Blends and High-Reactivity Ether Bioblendstocks

In this study, the impacts of fuel volatility and reactivity on combustion stability and emissions were studied in a light-duty single-cylinder research engine for a three-injection catalyst heating operation strategy with late post-injections. Here, N-heptane and blends of farnesane/2,2,4,4,6,8,8-heptamethylnonane were used to study the impacts of volatility and reactivity. The effect of increased chemical reactivity was also analysed by comparing the baseline #2 diesel operation with a pure blend of mono-ether components (CN > 100) representative of potential high cetane oxygenated bioblendstocks and a 25 vol.% blend of the mono-ether blend and #2 diesel with a cetane number (CN) of 55. At constant reactivity, little to no variation in combustion performance was observed due to differences in volatility, whereas increased reactivity improved combustion stability and efficiency at late injection timings. Fuels with higher reactivity were found to reduce engine-out hydrocarbon and carbon monoxide emissions while also achieving stable combustion at post-injection timings later than those achievable with #2 diesel fuel. The pure ether blend had the latest achievable post-injection timing of +30.5 CAD while still maintaining stable combustion (coefficient of variation of gross-indicated mean effective pressure < 5%). With post-injection timing adjusted to achieve a matched exhaust temperature of 300 °C, the ether-diesel CN 55 blend was observed to have slightly higher thermal efficiency in comparison to the baseline #2 diesel fuel. The results also indicate that cetane number may serve as a good indicator of combustion characteristics at late injection timings used for aftertreatment thermal management operation.

09 BIOMASS FUELS↗
DOE OSTI2023

Exploration of Fuel Property Impacts on the Combustion of Late Post Injections Using Binary Blends and High-Reactivity Ether Bioblendstocks

In this study, the impacts of fuel volatility and reactivity on com- bustion stability and emissions were studied in a light-duty single- cylinder research engine for a three-injection catalyst heating oper- ation strategy with late post-injections. N-heptane and blends of farnesane/2,2,4,4,6,8,8-heptamethylnonane were used to study the im- pacts of volatility and reactivity. The effect of increased chemical reac- tivity was also analysed by comparing the baseline #2 diesel operation with a pure blend of mono-ether components (CN > 100) representa- tive of potential high cetane oxygenated bioblendstocks and a 25 vol.% blend of the mono-ether blend and #2 diesel with a cetane number (CN) of 55. At constant reactivity, little to no variation in combustion performance was observed due to differences in volatility, whereas in- creased reactivity improved combustion stability and efficiency at late injection timings. Fuels with higher reactivity were found to reduce engine-out hydrocarbon and carbon monoxide emissions while also achieving stable combustion at post-injection timings later than those achievable with #2 diesel fuel. The pure ether blend had the latest achievable post-injection timing of +30.5 CAD while still maintaining stable combustion (coefficient of variation of gross-indicated mean ef- fective pressure < 5%). With post-injection timing adjusted to achieve a matched exhaust temperature of 300 °C, the ether-diesel CN 55 blend was observed to have slightly higher thermal efficiency in comparison to the baseline #2 diesel fuel. The results also indicate that cetane num- ber may serve as a good indicator of combustion characteristics at late injection timings used for aftertreatment thermal management opera- tion.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2023

Experimental investigation of power transient flow boiling

To advance the mechanistic understanding of power transients in flow boiling, experimental investigations are performed to capture the thermal response of four different kinds of cladding materials, including FeCrAls, to Fuchs reactivity initiated accident transients and linear ramp power transients. The transition boiling regime is phenomenologically reflected in power-transient flow boiling, which is different from steady-state boiling. In light of this, the critical heat flux that is usually featured with temperature overshooting does not indicate thermal safety margin of cladding materials as conservatively as the maximum heat flux, which is greater than critical heat flux in the power transient boiling curve. This is physically attributed to the thermal energy deposition on the cladding wall. In addition, the thermal energy released from the cladding wall results in different boiling heat transfer coefficient for the decreasing and increasing power stages respectively during the Fuchs power transient. The maximum heat flux difference gap, which appears due to changes in cladding materials and transient time scales, is appreciable under an intermediate heat convection regime. However, this difference gap is gradually weakened by the progressive increasing of mass flux and/or inlet subcooling due to the enhanced dominance of heat convection over heat conduction. Finally, the small difference gap of maximum heat flux is found to be insignificant under the weak heat convection because of power transient induced annular flow instability.

42 ENGINEERING↗
DOE OSTI2024

Final Report on Predictive Analyses of PRD as Function of Anomalies

This report summarizes the work completed in FY-2024 to analyze the power reactivity decrement (PRD) concepts of the ARC-100 core. The PRD has been traditionally defined by the reactivity change from a hot zero power (HZP) to a particular power. Consequently, the PRD accounts for the core reactivity changes due to increase coolant temperature gradient axially and radially across the core, and increased fuel temperature. The coolant temperature gradient leads to sodium and structure density changes, to radial core expansion from assembly flowering and bowing (due to axial and radial temperature gradients within the assemblies), and to control rod driveline thermal expansion. The fuel temperature increase associated with coolant temperature and power increases leads to Doppler effect and axial thermal expansion. In this work, the normal operating Hot Full Power (HFP) state is the only power level of interest, so analyses focus on the PRD calculated from HZP to HFP. The PRD has been used to assess the reactor safety features asymptotically in unprotected accident scenarios, including the loss of heat sink (LOHS), loss of flow (LOF), and transient overpower (TOP) without scram. The PRD concept relies on the “global” reactivity coefficients A, B, and C that are estimated based on “individual” reactivity effects (Doppler, sodium density, etc.). The objectives of this work are: 1) to improve and verify the methodology used to calculate the ABC coefficients used in the PRD, 2) to assess if the PRD can be used to reliably identify abnormal events. This report fulfills the FY-2024 scope of WBS#1.15.8.3 activity, “ANL0120 – Predictive Analyses of PRD as function of deformation”. The PRD concept is described in Sections 2. Additional effort in refining the methodology for core bowing modeling is performed in Sections 3. Then, two verification exercises are proposed in Section 4 to benchmark these coefficients based on direct neutronic-only calculations and on dynamic core transient simulations. Finally, the PRD approach is assessed for the detection of several unexpected events, such as primary flow perturbation or improper fuel loading, in Section 5.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2026

Material‐Driven Neuronal Oscillators and Filters via Active Reactance in CC‐NDR and VC‐NDR Electro‐Thermal Memristors

The continued scaling of artificial intelligence and telecommunications hardware is increasingly constrained by the power, bandwidth, and area limitations of transistor-based circuits. Neuromorphic processor units, analog oscillators, and active inductors and capacitors rely on complex multi-transistor architectures restricting material choices and incurring energy and footprint overhead. Here, we show that active reactance in electro-thermal memristors provides an intrinsic, material driven route to neuronal oscillator dynamics and signal processing. Using a physics-based compact modeling framework, we bridge negative differential resistance (NDR) and bias-tunable reactance, which underlies spiking dynamics in electro-thermal memristors. Memristors with negative temperature coefficients of resistance (TCR) manifest current-controlled (CC-) NDR and act as active inductors, thus generating spiking above a critical circuit capacitance; whereas memristors with positive TCR manifest voltage-controlled (VC-) NDR and active capacitance, leading to spiking above a critical inductance. By creating a compact model for La 0.7 Ca 0.3 MnO 3 as a representative VC-NDR material and comparing it with LaCoO 3 manifesting CC-NDR, we explain the physical origins of their distinct current-voltage characteristics, reactive phase shifts and consequent spiking behaviors. Finally, we demonstrate tunable filtering enabled by the active reactance of electro-thermal memristors, establishing them as a compact hardware platform for neuronal oscillator functionality and integrated filtering beyond conventional CMOS.

active reactance↗
DOE OSTI2024

Uncertainty Quantification and Sensitivity Analysis of Non-Nuclear Advanced Controls Testbed Reactor Mockup

The research presented in this report describes our progress in applying stochastic methods and uncertainty quantification, parametric study, and variance-based sensitivity analysis (also known as Sobol sensitivity analysis) to a full-core model of a nuclear thermal propulsion (NTP) system simulated with Griffin, with the goal of developing a reduced order (surrogate) model which can be rapidly sampled while perturbing multiple input parameters. In this NTP system, reactivity and power feedback affect the rotation of control drums, which are controlled by a hybrid proportional, integral and derivative (PID) controller, actuated by the power demand and reactivity feedback from the numerical model. This model uses reactor kinetic feedback (mean generation time and $\beta$ from a transient Griffin simulation executed with the improved quasi-static method to provide the kinetic parameters) as inputs to functions which control the CD rotation angle. Using a number of stochastic method approaches, we developed a dual purpose training-surrogate model of the NTP system using polynomial regression. The trained model can be rapidly sampled while simultaneously perturbing various input parameters of the model, such as coefficients on the PID control, or temperature (directly affect the neutron cross section). The surrogate model delivers accurate results orders-of-magnitude faster (minutes, not days) than the base model. Once the base model has been trained, distributions of the uncertain parameters can be changed at will to investigate the effects of perturbing multiple inputs and their effect on the output. For example, coefficients used in the PID control system may vary due to some physical interference, or there may be uncertainty in the temperature of the neutron cross sections in various regions of the reactor. A distribution can be placed on these parameters and operational boundaries can be determined. The goal of this work is to support development of an advanced control system to operate CDs in a functioning NTP system.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗