Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fast reactors”

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 1,369 records · Page 76

SCALE 6.2.4 Validation: Nuclear Criticality Safety

The computational bias of criticality safety computer codes must be established through the validation of the codes to critical experiments. A large collection of suitable experiments has been vetted by the International Criticality Safety Benchmark Evaluation Project (ICSBEP) and made available in the International Handbook of Evaluated Criticality Safety Benchmark Experiments (ICSBEP Handbook). More than 600 cases from this handbook have been prepared and reviewed within the Verified, Archived Library of Inputs and Data (VALID), which is maintained by the Reactor and Nuclear Systems Division at Oak Ridge National Laboratory. The performance of the KENO V.a and KENO-VI Monte Carlo codes within the SCALE 6.2.4 code system is assessed using the VALID models of benchmark experiments. A range of nuclear cross section libraries based on Evaluated Nuclear Data File (ENDF)/B-VII.1 in both multigroup (MG) and continuous energy (CE) formats is considered. The critical experiments available to validate the KENO V.a code cover 15 broad categories of systems. These systems use a range of fissile materials, including a range of uranium enrichments, various plutonium isotopic vectors, and some mixed uranium/plutonium oxides. The physical forms of the fissile material also vary and are represented as metal, solutions, or arrays of rods or plates in a water moderator. The neutron energy spectra of the systems also vary and cover fast, intermediate, mixed, and thermal spectra. Over 550 of the total cases use the KENO V.a code for the four nuclear data libraries considered in this report.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

AGC-4 As-Run Uncertainty Analysis

The concept of measuring neutron flux (or fluence rate) with activation foils has existed since the late 1940s and is a well-used methodology. Activation foils are selected for their unique isotopic cross sections that create reaction products with desirable properties: relatively simple radioactive decay scheme, measurable quantities of gamma-rays or betas, and half-lives appropriate for the irradiation and decay durations of the experiment. The activation materials, deemed flux wires, fluence wires, dosimeters, or flux foils, are fabricated to have a mass appropriate for the irradiation as well so the measurement using a radiation detector following the irradiation is sufficient but not overwhelming in terms of total activity. Activation measurement results are expressed in terms of reaction rate, which can then be used to “adjust” the neutron spectrum. The Advanced Graphite Creep (AGC)-4 experiment contained 26 flux wire packages, each containing an iron and a niobium wire inside a vanadium capsule. For all stackup positions, the value of fast to thermal neutrons is higher than expected, but the total neutron fluence is near the predicted value of ~2 x 1022 n/cm2. As compared to the predicted neutron spectrum, there are four takeaways in the results as compared to the predicted ATR irradiation: 1) Additional fast (fission) neutrons present, 2) Spectrum hardening due to temperature, 3) Spectrum hardening due to container attenuation, 4) Total neutron fluence is close to predicted values. Since the neutron energy spectrum is a key parameter when attributing neutron damage in materials, the adjusted energy-dependent spectrum and the total neutron fluence results from the dosimeters can be incorporated for future calculations.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Thermonuclear neutron emission from a sheared-flow stabilized Z-pinch

We report the fusion Z-pinch experiment (FuZE) is a sheared-flow stabilized Z-pinch designed to study the effects of flow stabilization on deuterium plasmas with densities and temperatures high enough to drive nuclear fusion reactions. Results from FuZE show high pinch currents and neutron emission durations thousands of times longer than instability growth times. While these results are consistent with thermonuclear neutron emission, energetically resolved neutron measurements are a stronger constraint on the origin of the fusion production. This stems from the strong anisotropy in energy created in beam-target fusion, compared to the relatively isotropic emission in thermonuclear fusion. In dense Z-pinch plasmas, a potential and undesirable cause of beam-target fusion reactions is the presence of fast-growing, “sausage” instabilities. This work introduces a new method for characterizing beam instabilities by recording individual neutron interactions in plastic scintillator detectors positioned at two different angles around the device chamber. Histograms of the pulse-integral spectra from the two locations are compared using detailed Monte Carlo simulations. These models infer the deuteron beam energy based on differences in the measured neutron spectra at the two angles, thereby discriminating beam-target from thermonuclear production. An analysis of neutron emission profiles from FuZE precludes the presence of deuteron beams with energies greater than 4.65 keV with a statistical uncertainty of 4.15 keV and a systematic uncertainty of 0.53 keV. This analysis demonstrates that axial, beam-target fusion reactions are not the dominant source of neutron emission from FuZE. These data are promising for scaling FuZE up to fusion reactor conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of Multiresolution Capabilities for the Holistic Energy Resource Optimization Network (HERON) tool A progress update

INL researchers work on technoeconomic analyses for integrated energy systems (IES) using the Framework for Optimization of ResourCes and Economics (FORCE). Within FORCE, researchers use the Holistic Energy Resource Optimization Network (HERON) tool to conduct optimization of grid portfolios under uncertain market conditions. These optimizations determine optimal capacities for all IES components and strategies for resource dispatch which maximize some economic metric (e.g., net present value). Resource dispatch occurs on finer timescales (typically hours) and thus are asked to respond to a given time series (e.g. hourly load demand profiles for a grid, or pre-determined electricity prices). Volatile and complex bidding dynamics as well as poorly forecasted weather events within deregulated markets add uncertainty to the time series; FORCE can address this uncertainty by training a reduced order model on historical time series and generate unique synthetic time series which represent individual scenarios or realizations of the market. The IES configuration can be simulated under these different sampled realizations and a stochastic optimization is conducted which optimizes the expected value of the desired economic metric. The training of a synthetic time series generator is limited by the chosen time resolution; dynamics can occur on different time scales. Seasonal demand trends can dominate faster dynamical events (such as power outages from certain sectors or severe weather events) which might not get captured correctly by the trained model. In this report, we investigate different ways of addressing the training and generation of time series on multiple time scales using three main algorithms: wavelet decomposition, dynamic mode decomposition, and generative adversarial networks for time series. We demonstrate a time series analysis that yields information on not just the frequency space but also temporal space: where a fast Fourier transform can provide what frequencies dominate, the new algorithms can provide when the frequencies dominate as well. These analyses can help improve IES optimization by allowing researchers to couple simulations at different timescales when it is most needed - seasonal, day-ahead, and real time optimization - with greater computational efficiency. Future work will include implementation of a subset of the proposed algorithms into the FORCE toolset and application of these analyses into multiple timescale optimization.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Catalytic co-pyrolysis of LDPE and PET with HZSM-5, H-beta, and HY: experiments and kinetic modelling

In this study, the catalytic pyrolysis of low-density polyethylene (LDPE), polyethylene terephthalate (PET) and their mixture (1 : 1 wt/wt) with three zeolite catalysts (HZSM-5, H-beta, HY) was investigated using a thermogravimetric analyzer (TGA) and a Pyroprobe® micro-reactor coupled to a gas chromatograph mass spectrometer (Py-GC/MS). The TGA results demonstrated that during pyrolysis at 10 °C min –1 , on average, zeolite catalysts decreased the maximum decomposition temperature by 149 °C for LDPE while only decreasing by 8 °C for PET. The derivative thermogravimetric (DTG) curve evidenced interactions when the two polymers were catalytically co-pyrolyzed for all the three catalysts. A lumped nth order reaction scheme was able to accurately model both non-catalytic and catalytic pyrolysis and co-pyrolysis by using the least squares fitting approach for determining the kinetic parameters. The kinetic model was able to model well the interaction effects observed during catalytic co-pyrolysis of LDPE and PET with HZSM-5, H-beta, and HY (Fit% Wt% > 96%, Fit% DTG > 94%). Py-GC/MS experiments for the catalytic fast pyrolysis of LDPE revealed that HZSM-5 resulted in the highest selectivity to aromatic hydrocarbons (31.6%) and HY resulted in the highest selectivity to gasoline range C 5 –C 10 paraffins and olefins (40.9%). Catalytic fast pyrolysis of PET showed high selectivity to benzene for all the catalysts (>43%) and that HZSM-5 resulted in the highest selectivity to polyaromatic hydrocarbons (24.7%). Finally, the catalytic fast co-pyrolysis of LDPE and PET revealed interaction effects for all the three catalysts evidenced by a positive synergy% for alkylated benzenes (3–142%) and polyaromatics (105–187%) with a concomitant negative synergy% for benzene (24–36%) and C 5 –C 10 paraffins and olefins (27–53%).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Advanced two-phase subchannel method via non-linear iteration

A fast-running, robust two-phase flow, sub-channel model is presented based on non-linear solution of the steady-state subchannel fluid flow equations. The drift-flux model solves for conservation of liquid and vapor mass, mixture energy, and axial and transverse mixture momentum as part of an efficient planar marching scheme and nonlinear, nested outer and inner iteration. Here, models based on mechanistic subcooled boiling, two-phase turbulent void mixing, and drift are included. Solution verification and mesh convergence studies were performed for modern GE 10 × 10 fuel geometry and are shown to have excellent convergence behavior. Run time performance for a 50 axial mesh model showed 2.2 seconds on a single CPU core to tightly converge all 3D distributions (flow, void, pressure) for the GE 10 × 10 fuel geometry, supporting its efficient use within the Virtual Environment for Reactor Applications boiling water reactor framework.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials: Part II Resveratrol Exemplar

This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.

36 MATERIALS SCIENCE↗

AGR-5/6/7 Irradiation Test Final As-run Report

This document presents the as-run analysis of the Advanced Gas Reactor (AGR)-5/6/7 irradiation experiment. AGR-5/6/7 is the last of a series of experiments conducted in the Advanced Test Reactor (ATR) at Idaho National Laboratory in support of the development and qualification of tri-structural isotropic low-enriched fuel for use in high-temperature gas-cooled reactors. The test train contained five separate capsules that were independently controlled and monitored. Each capsule contained multiple 24.91-mm-long and 12.25-mm-dimeter compacts filled with low-enriched uranium carbide/oxide tri-structural isotropic fuel particles. The objectives of the AGR-5/6/7 experiment were to: • Irradiate reference-design fuel particles to support fuel qualification. • Establish operating margins for the fuel, beyond normal operating conditions. • Provide irradiated-fuel performance data and irradiated-fuel samples for post-irradiation examination and safety testing. The primary objective of the AGR-5/6 test (Capsules 1, 2, 4, and 5) was to verify the successful performance of the reference-design fuel under normal operating conditions. The AGR-7 test (Capsule 3) was designed to explore fuel performance at higher temperatures. Its primary objective was to demonstrate the capability of the fuel to withstand conditions beyond normal operating conditions, in support of plant design and licensing. AGR-5/6/7 will also provide irradiated-fuel performance data based on the fission gas release from particles during irradiation. To achieve the test objectives, the AGR-5/6/7 experiment was irradiated in the northeast flux trap of the ATR with a planned duration of 500 effective full-power days. The northeast flux trap was selected because its larger diameter provided greater flexibility for test-train design compared to the Large B positions used for the AGR-1 and AGR-2 irradiations, significantly enhancing test capabilities for the combined irradiation campaigns. Due to delays in the ATR schedule, the AGR-5/6/7 irradiation was significantly shorter than the originally planned 13-cycle schedule. Irradiation began on February 16, 2018 and ended on July 22, 2020, spanning nine ATR cycles (162B–168A) over two and a half years. Thus, the AGR-5/6/7 fuel compacts were irradiated for a total of approximately 360.9 effective full-power days. Final burnup values, on a per-compact basis, ranged from 5.66 to 15.26% fissions per initial heavy metal atom, while fast fluence values ranged from 1.62 to 5.55 × 1025 n/m2 (E >0.18 MeV). Time-averaged volume-averaged fuel temperatures on a capsule basis at the end of irradiation ranged from 756°C in Capsule 5 to 1313°C in Capsule 3 excluding days with significantly lower temperature during the two short powered axial locator mechanism cycles, 163A and 167A. By the end of irradiation, 48 out of 54 installed thermocouples had failed (the bottom three capsules lost all thermocouples). During the first five cycles (162B – 165A), the fission-gas isotope release-rate-to-birthrate (R/B) ratios were stable in the 10-8–10-6 range, and no in-pile particle failures were observed based on the gross gamma counts. During this time, the high exposed kernel fraction and high fuel particle temperatures in Capsule 1 led to a maximum R/B value of around 2 ? 10-6 for Kr-85m. The fission gas release in all capsules started to increase from the second half of Cycle 166A, when a large number of in-pile particle failures occurred in Capsule 1 and a gas line problem in this capsule caused fission gas leakage at various degrees into the other four capsules. This gas line problem also prevented a fission gas release measurement for Capsule 1 during the last three cycles due to gas flow isolation. By the end of irradiation, it is estimated that approximately 15 particles failed in Capsule 3, which was considered possible because the experiment was designed to operate beyond the high-temperature gas-cooled reactor normal operating temperature range. A few hundred in-pile particle failures were estimated for Capsule 1 by the end of Cycle 166A, but the total number of failures is unknown due to the lack of fission gas release data in the later cycles. Additionally, four potential in-pile failures were identified for Capsule 2 during the last cycle, Cycle 168A. In contrast, no in-pile failures were identified in the top two capsules (4 and 5) based on the absence of the typical spikes in gross gamma counts and low failure estimates using the AGR-3/4 R/B per exposed kernel model.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Software Framework for Federated Science Instruments

There is an unprecedented promise of enhanced capabilities for federations of leadership computing systems and experimental science facilities by leveraging software technologies for fast and efficient operations. These federations seek to unify different science instruments, both computing and experimental, to effectively support science users and operators to execute complex workflows. The FedScI project addresses the software challenges associated with the formation and operation of federated environments by leveraging recent advances in containerization of software and softwarization of hardware. We propose a software framework to streamline the federation usage by science users and it’s provisioning and operations by facility providers. A distinguishing element of our work is the support for improved interaction between experimental devices, such as beam-line instruments, and more traditional high-performance computing resources, including compute, network, storage systems. We present guiding principles for the software framework and highlight portions of a current prototype implementation. We describe our science use case involving neutron imaging beam-lines (SNAP/BL-3, Imaging/CG-1D) at the Spallation Neutron Source and High Flux Isotope Reactor facilities at Oak Ridge National Laboratory. Additionally, we detail plans for a more direct instrument interaction within a federated environment, which could enable more advanced workflows with feedback loops to shorten the time to science.

Naughton III, Thomas↗

Assessment of a Detailed Biomass Pyrolysis Kinetic Scheme in Multiscale Simulations of a Single-Particle Pyrolyzer and a Pilot-Scale Entrained Flow Pyrolyzer

A detailed biomass pyrolysis kinetic scheme was assessed in the multiscale simulations of a single-particle pyrolyzer with slow pyrolysis and a pilot-scale entrained flow pyrolyzer with fast pyrolysis. The detailed kinetic scheme of biomass pyrolysis developed by the CRECK group consists of 32 reactions and 58 species. A multiscale simulation model was developed, where the CRECK kinetics was employed to simulate biomass pyrolysis reactions, a one-dimensional particle model was utilized to simulate the intraparticle transport phenomena, and the particle-in-cell (PIC) model was employed to simulate the hydrodynamics. The multiscale model was first applied to simulate a single-particle pyrolysis experiment. The simulation with nonisothermal particles matched the experimental data better than the simulation with isothermal particles. Then the multiscale model was applied to simulate the pilot-scale entrained flow pyrolyzer. In this case, the simulation with isothermal particles matched the experimental data better than the simulation with nonisothermal particles. The reason for this difference might be that the kinetics itself already partially included the intraparticle transport effect as it was fitted using both TGA data (slow pyrolysis of small size biomass) and fluidized bed data (fast pyrolysis of relatively large size biomass). This study provides some insights into biomass pyrolysis kinetics development and pyrolyzer multiscale simulation for a future study.

09 BIOMASS FUELS↗

AGR-5/6/7 Irradiation Test Final As-run Report

This document presents the as-run analysis of the Advanced Gas Reactor (AGR) 5/6/7 irradiation experiment. AGR-5/6/7 is the last of a series of experiments conducted in the Advanced Test Reactor (ATR) at Idaho National Laboratory (INL) in support of the development and qualification of tri-structural isotropic (TRISO) low enriched fuel for use in high-temperature gas-cooled reactors (HTGR). The test train contained five separate capsules that were independently controlled and monitored. Each capsule contained multiple 24.91-mm-long and 12.25-mm-dimeter compacts filled with low-enriched uranium carbide/oxide tri-structural isotropic fuel particles. The objectives of the AGR 5/6/7 experiment were to: 1. Irradiate reference-design fuel particles to support fuel qualification. 2. Establish operating margins for the fuel, beyond normal operating conditions. 3. Provide irradiated-fuel performance data and irradiated-fuel samples for post irradiation examination and safety testing. The primary objective of the AGR-5/6 test (Capsules 1, 2, 4, and 5) was to verify the successful performance of the reference-design fuel under normal operating conditions. The AGR-7 test (Capsule 3) was designed to explore fuel performance at higher temperatures. Its primary objective is to demonstrate the capability of the fuel to withstand conditions beyond normal operating conditions, in support of plant design and licensing. AGR 5/6/7 will also provide irradiated-fuel performance data based on the fission gas release from particles during irradiation. To achieve the test objectives, the AGR-5/6/7 experiment was irradiated in the northeast (NE) flux trap of the ATR with a planned duration of 500 effective full-power days (EFPDs). The NE flux trap was selected because its larger diameter provided greater flexibility for test-train design compared to the Large B positions used for the AGR-1 and AGR-2 irradiations, significantly enhancing test capabilities for the combined irradiation campaigns. Due to delays in the ATR schedule, the AGR-5/6/7 irradiation was significantly shorter than the originally planned 13-cycle schedule. Irradiation began on February 16, 2018 and ended on July 22, 2020, spanning nine ATR cycles (162B?168A) over two and a half years. Thus, the AGR-5/6/7 fuel compacts were irradiated for a total of approximately 360.9 EFPDs. Final burnup values, on a per-compact basis, ranged from 5.66% to 15.26% fissions per initial heavy metal atom, while fast fluence values ranged from 1.62 to 5.55 × 1025 n/m2 (E >0.18 MeV). Time-averaged volume-averaged fuel temperatures on a capsule basis at the end of irradiation ranged from 756°C in Capsule 5 to 1313°C in Capsule 3 excluding days with significantly lower temperature during the two short powered axial locator mechanism cycles, 163A and 167A. By the end of irradiation, 48 out of 54 installed thermocouples (TC) had failed (the bottom three capsules lost all TCs). During the first five cycles (162B ? 165A), fission-gas isotope release-rate-to-birthrate (R/B) ratios were stable in the 10-8?10-6 range, and no in-pile particle failures were observed based on the gross gamma counts. During this time, the high exposed kernel fraction and high fuel particle temperatures in Capsule 1 led to the maximum R/B value of around 2 ? 10-6 for Kr-85m. The fission gas release in all capsules started to increase from the second half of Cycle 166A, when a large number of in-pile particle failures occurred in Capsule 1 and a gas line problem in this capsule caused fission gas leakage at various degrees into the other four capsules. This gas line problem also prevented fission gas release measurement for Capsule 1 during the last three cycles due to gas flow isolation. By the end of irradiation, it is estimated that approximately 15 particle failures in Capsule 3, which was considered possible at the start of the irradiation, since the experiment was designed to operate beyond the high-temperature gas-cooled reactor (HTGR) normal operating temperature range. A few hundred in-pile particle failures were estimated for Capsule 1 by the end of Cycle 166A, but the total number of failures is unknown due to lack of fission gas release data in the later cycles. Additionally, four potential in-pile failures were identified for Capsule 2 during the last cycle, 168A. In contrast, no in-pile failures were identified in the top two capsules (4 and 5) based on absence of the typical spikes in gross gamma counts and low failure estimates using the AGR-3/4 R/B per exposed kernel model.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Energy-Efficient Microbial Electrochemical Lignin and Alkaline Hydroxide Recovery from DMR Black Liquor

Black liquor is a complex waste stream generated from the fast-expanding biorefinery industry, and the high management cost has been known as a major development barrier. In this work, we demonstrated a rationally designed 4-chamber microbial electrolysis process that was capable of recovering the value-added lignin, NaOH, and H2 products while removing waste organics from the Deacetylation and Mechanical Refining (DMR) black liquor, which could potentially reduce the disposal cost by up to 40%. The lab-scale MEC achieved 52% organic removal from a highly concentrated feedstock (28,960 ± 212 mg/L COD) and recovered 83.2 mM NaOH under a small voltage application (1.2 V). Lignin recovery was observed to be as high as 13.98 +/- 0.28 g/L, and 0.35 L/L/day H2 was produced with a faradaic efficiency of 93.5%. This work provided a good example on how complex liquid waste stream can be valorized via multi-functional microbial electrochemical technology, and how the traditionally unfavored pH gradient in an electrochemical reactor can be manipulated for lignin precipitation and alkaline recovery at the same time. The versatility and good performance in concurrent contaminants removal, lignin precipitation, H2 production and NaOH recovery demonstrate a good potential for microbial electrolysis in black liquor treatment.

BIOMASS FUELS,INORGANIC, ORGANIC, PHYSICAL, AND AN↗

Optimizing Ammonia Separation via Reactive Absorption for Sustainable Ammonia Synthesis

Metal halide salts such as magnesium chloride have been shown to be promising candidates for ammonia storage materials for energy storage and agriculture applications due to their ability to incorporate several moles of ammonia per mole of salt. Ammonia exiting a synthesis reactor can be separated from nitrogen and hydrogen by absorption into magnesium chloride. Such an absorption can be more complete and hotter than separation via ammonia condensation, the current standard in the Haber–Bosch process. Here, we discuss the optimal conditions for the cyclic uptake and release of ammonia from the supported magnesium chloride absorbents. An automated system was designed for measuring the nonequilibrium working capacity of the absorbent, as well as the impact of important operating conditions such as absorption and desorption temperature, pressure, and desorption time. Measurements of absorption and desorption kinetics provide insight into the mechanisms involved. The temperatures and pressures during absorption and desorption were designed to use minimal energy input to maximize the uptake and release of ammonia within a reasonable amount of time. In a laboratory-scale bed, absorption has a small unused bed length, so it is largely independent of mass transfer; it is dominated by how fast ammonia is fed to the bed. Yet, desorption is restricted both by the speed of heating the bed and by diffusion out of the absorbent. These measurements provide guidelines for ammonia separations and cycling sorbent materials on a larger scale.

10 SYNTHETIC FUELS↗

Multiphysics Modeling of Microreactors with NEAMS codes, and Validation Based on KRUSTY Reactivity Insertion

The NEAMS Multiphysics Applications team continues to assess code usability and functionality for microreactor design and safety analyses, while demonstrating that NEAMS tools capture both steady-state and transient behavior across distinct microreactor concepts. In FY2025, the team advanced full-core, high-fidelity, multiphysics models that solve more complex problems and strengthen verification/validation for several microreactor systems: heat-pipe microreactor (HPMR), gas-cooled microreactor (GCMR), and the KRUSTY experiment. These models employ the MOOSE MultiApp/Transfers architecture with Griffin for neutronics, BISON for heat conduction/thermomechanics, Sockeye for heat pipes, SAM/THM for coolant channels and loops, and SWIFT for hydride behavior, with meshes generated via the MOOSE Reactor Module. The graphite models available in the Grizzly code were also investigated for future analyses. For the HPMR, a Na-HPMR variant was constructed to align with recently validated heat-pipe experiments and Sockeye’s LCVF capability, enabling mechanistic heat-pipe transients and startup modeling. The Na-HPMR will serve as the primary model for HPMR investigations in upcoming tasks. The load-following and single heat-pipe failure scenarios (Griffin/BISON/Sockeye), which were previously modeled for the K-HPMR, were replicated for the Na-HPMR, showing strong negative temperature feedback and highly localized thermal effects, respectively, while the startup case captured vapor-front progression and heat-removal activation. Solid mechanics was added to the previously built K-HPMR full-core model in BISON, showing minimal impact on steady-state reactivity yet enabling stress-field predictions that prepare the path for full-core TRISO performance analyses. For the GCMR, automated steady-state and four transient scenarios were executed using Griffin/BISON/SAM/SWIFT. Results confirm robust inherent safety: power collapses promptly in loss-of-cooling events, the inlet-temperature drop settles to a new equilibrium, and a single-channel blockage yields only a ~30 K local fuel-temperature rise with <0.4% power decrease. SWIFT-predicted hydrogen redistribution affects reactivity during both steady-state and transient conditions, underscoring its importance. A Brayton-cycle balance of plant (BOP) model in SAM/THM demonstrated stable startup behavior, and xenon-driven reactivity during load following was analyzed. To improve TRISO-compact temperature fidelity, a fast multiscale Heat Source Decomposition (HSD) treatment was implemented. Against heterogeneous benchmarks, HSD reduces underprediction of kernel temperatures and lowers predicted peak powers in reactivity-insertion transients compared to previous homogenized models. KRUSTY warm-critical validation progressed from FY2024 baselines: the 15Ȼ insertion shows excellent agreement in peak power (~2% high) and temperature trends, and the 30Ȼ case was automated via a feedback controller that maintained power near 3 kW for ~150 s with close agreement to data. The successful modeling of the warm critical tests has laid a strong foundation for simulating more complex nuclear system tests in the years ahead. Throughout FY2025, developer feedback was provided (e.g., MOOSE batch mesh generation, distributed pre-split meshes, Griffin sweeper on displaced meshes), several new models were contributed to the Virtual Test Bed, and an OECD-NEA WPRS multiphysics benchmark based on the HPMR was initiated to enable broader cross-comparison and best-practice development with the nuclear community at large.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Radiation-induced transients on solid and molten Li and K iodides

The study, quantification, and monitoring of iodine in nuclear power generation systems have been paramount due to its environmental and health implications. Thorough speciation and transport studies strengthened the development and deployment of water reactors, which are now required for the new generation of Molten Salt Reactors (MSR), where the fuel and coolant are a complex liquid system of suitable halides. The speciation, radiation-induced chemistry, and transport of iodine (I-127, I-129 and I-131) under extreme conditions of temperature and ionizing radiation are among the challenges for the future implementation of MSR’s. In aqueous systems, the triiodide-iodine-iodide (I3--I2-I-) equilibrium is strongly dependent on the conditions of the medium. In molten iodides and chlorides, the triiodide detection (and its disproportionation) poses challenges due to the extreme conditions of temperature and radiation, but studies have shown a dependence on temperature, absorbed dose and salt composition. Through absorption spectroscopy techniques (diffuse reflectance and high-temperature UV-Vis), as well as resonance (Electron Paramagnetic Resonance) it is possible to identify primary transients induced by ionizing radiation trapped in the solid phase (Solid-state radiolysis), and steady-state main absorption components with a fast-acquisition real-time spectroscopy. The presence of heterogeneous phase (colloids), viscosity changes, and gas evolution are among the observed processes. A specific case of study is speciation in the presence of Ni. A corrosion product with magnetic properties due to unpaired electrons in its electronic structure. Ferromagnetic materials are monitored in solidified salts to identify aggregation and phase separation. Comparative studies of eutectic mixtures of LiI-KI and LiCl-KCl with known addition of dopants will be presented, with a perspective of radiation chemistry and its possible implications in complex systems.

36 - MATERIALS SCIENCE↗

Open‐Source Anaerobic Digestion Modeling Platform, Anaerobic Digestion Model No. 1 Fast (ADM1F)

An open‐source modeling platform, called Anaerobic Digestion Model No. 1 Fast (ADM1F), is introduced to achieve fast and numerically stable simulations of anaerobic digestion processes. ADM1F is compatible with an iPython interface to facilitate model configuration, simulation, data analysis, and visualization. Faster simulations and more stable results are accomplished by implementing an advanced open‐source library of numerical methods called Portable Extensive Toolkit for Scientific Computation (PETSc) to solve the ADM1 system of equations. Leveraging PETSc, ADM1F can consistently complete a steady‐state simulation under 0.2 s, over 99% faster than a benchmark ADM1 model implemented with MATLAB while achieving agreement of model outputs within 1% of those obtained with the benchmark model. For dynamic simulations, however, ADM1F has a computational speed advantage only when the influent characteristics update more frequently than every 4 h. The ability of ADM1F to be useful as a tool to study anaerobic digestion systems is demonstrated through two example implementations of ADM1F: (1) a two‐phase co‐digestion scenario evaluating the impact of the organic loading rate and the substrate composition on reactor performance and stability, and (2) a conventional digester scenario assessing the effectiveness of recovery strategies after disruptions that led to instability. These examples demonstrate how the high simulation speed and the convenience of the iPython interface allow ADM1F to complete complex analyses within minutes, much faster than computational strategies currently reported in the literature.

anaerobic co-digestion↗

Optimization of lithium vapor box divertor evaporator location on NSTX-U using SOLPS-ITER

Commercial fusion reactors will be faced with extremely high divertor target heat fluxes that will require mitigation. Simulations of detachment in an NSTX-U scenario projected to have 92 MW/m 2 unmitigated peak target heat flux are presented, which reaches sub-10 MW/m 2 target heat flux using a highly dissipating lithium vapor box divertor design. The lithium vapor box is a detached divertor design which employs lithium vapor evaporation and condensation to contain lithium below the X-point. Previous SOLPS modeling has indicated a lithium vapor box can reduce the heat flux down to 10 MW/m 2 via simultaneous evaporation from the Private Flux Region (PFR) and the Common Flux Region (CFR) sides of the vapor box. It is found here that PFR evaporation has improved access to the separatrix leading to significantly more efficient power dissipation than CFR evaporation. Simulations of target evaporation with an evaporation distribution that is self-consistent with the temperature of a Capillary Porous System with Fast flowing liquid lithium (CPSF) could reach n Li / n e ~ 0.025-0.030 at the Last Closed Flux Surface (LCFS) depending on the liquid metal flow speeds and lithium sputtering yield, while PFR-side evaporation can reach acceptable heat fluxes with n Li / n e ~ 0.038 at the LCFS. However, PFR evaporator performance can be improved if the target is allowed to be hot enough such that it reflects lithium, reaching n Li / n e ~ 0.028 and reducing required lithium evaporation. Ultimately PFR evaporation and target evaporation are found to have similar ability to produce acceptable heat flux solutions with minimal upstream concentration.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simple Heat Pipe Model

This is a simple model designed to run fast but still maintain the key physics and feedback mechanisms of a heat pipe. First, the capillary pressure is a function of the liquid working fluid volume fraction. Second, the boiling and condensation are based on the saturation temperature that is based on the heat pipe pressure. When the pressure goes up, the saturation temperature goes up and the vapor rains on the wick. When the pressure goes down, the saturation temperature goes down and the liquid in the entire wick boils. This is how the heat pipe adjusts to stay robust under different temperatures and heat fluxes.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗