Variable Dynamic Mode Decomposition for Estimating Time Eigenvalues in Nuclear Systems
Not provided.
Nuclear Science & Technology↗
SEARCH · Engineering Papers
Results for “Reactor Kinetics”
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.
Rapid depletion analysis of flowing-pebble reactor systems at equilibrium using SCALE
Several high-temperature gas-cooled reactor concepts (and more recently, salt-cooled designs such as the fluoride salt-cooled high-temperature reactor) feature core designs employing continuously circulating fuel pebbles. These reactor designs permit both continuous online refueling of fuel elements as well as higher overall achievable discharge burnups. However, rapid calculation of time-dependent fuel isotopic inventories proves challenging for this class of dynamic systems with current analysis tools. While iterative approaches employing coupled neutron transport have been developed to solve this issue, rapid depletion analysis techniques are needed to calculate time-dependent inventories for individual pebbles and batches (and thus the construction of full- core inventory at equilibrium). We propose a depletion analysis strategy for this type of system for cores at equilibrium. Drawing upon previous neutronic analysis of the PBMR-400 equilibrium core, we demonstrate the viability of developing collapsed one-group cross section libraries suitable for performing rapid depletion analyses with SCALE. (authors)
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Stereolithographic geometry model of the IBR-2M experimental facility
The IBR-2M is a fast research reactor that operates in supercritical condition for ∼ 800 μs every 200 ms. Two reflector parts in nickel rotate in opposite directions generating 1.8 GWth peak power when they align with the fuel zone changing the reactor status from deep subcritical to supercritical. The reactor core uses high-enriched plutonium fuel and is cooled by sodium. This reactor has been modeled by MCNP and SERPENT computer programs. The MCNP model uses combinatorial geometry, whereas the SERPENT model employs the Stereolithographic (STL) geometry representation that can be used by 3D printers. The STL geometry was constructed using the CUBIT computer program. The CUBIT program was also used for a three-dimensional visualization of the Monte Carlo models. SERPENT and MCNP models use the same geometry, material specifications, and nuclear data. The latter are based on the ENDF/B-7.0 library. SERPENT and MCNP using same geometry and same material specifications produce similar k{sub eff} values within 120 pcm.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Exponential time differencing scheme for mass transport and depletion in molten salt reactors
This work extends the capability previously shown for addressing the problem of computing depletion and mass transport calculations in molten salt reactors (MSRs) by calculating matrix exponentials. Additional algorithms are implemented to compute the matrix exponential and the action of the matrix exponential on a matrix. These algorithms include two methods based on the Pade approximation, a Taylor series method, and three methods based on Cauchy's integral formula. In addition to the added matrix exponential solvers, a variable-order total variation diminishing scheme is applied to the convective flux approximation to provide enhanced accuracy. Finally, a simplified MSR problem is shown for each of the exponential time differencing solvers along with classical backwards differencing integrators. The results show excellent convergence for exponential time differencing methods. Computation time is a key element for selecting the optimal solver in these problems, and this work shows that Pade and Cauchy-based solvers may provided the fastest and most accurate solutions. (authors)
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Alternative analysis of the MINERVE ZPR oscillation experiments
Delayed neutrons are of fundamental importance in the field of nuclear reactor dynamics and control. However, the precursor yield fraction for a given nuclear reactor are dependent on the properties of the reactor. Thus, in-pile experiments, such as oscillation experiments are conducted in order to measure those values. In this work, an alternative analysis of the piston oscillation experiments that have been conducted in the MINERVE reactor in 2013 is performed. A new method which evolves effective terms that cancel out undesired drifts of the flux during the experiments is presented. The evaluation of the uncertainty on the values of the response function is also presented. Moreover, the effective delayed neutron fraction β{sub eff} is evaluated and is compared to results in previous works. As the analysis has led to an estimation of β{sub eff} with a large uncertainty, it has been deduced that the oscillation experiments that have been conducted in MINERVE are not a reliable method of experimentation to determine the value of β{sub eff} in the reactor and that the noise experiments are better suited for that purpose.
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Expanding Benchmarks for Nuclear Data Validation [Slides]
Abstract not provided.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
BETO 2021 Peer Review - Thermochemical Platform Analysis WBS: 2.1.0.302
The objective of the NREL Thermochemical Platform Analysis (WBS 2.1.0.302) project is to inform and guide R&D priorities for thermal and catalytic conversion processes by providing process design and techno-economic analysis (TEA). This is achieved through close collaboration with researchers and external experts, along with the use of both commercially available modeling tools and the development or use of collaboration-derived domain-specific tools and resources, such as refinery integration, kinetic and reactor models, phase equilibrium models, and pertinent bio-products market studies. This project is directly aligned with DOE-BETO goals, with the enabling of technology advancements and cost reduction for biomass derived biofuels being one of its primary objectives. TEA-guided research facilitated by this project has helped achieve significant modeled cost reductions for the ex situ catalytic fast pyrolysis (CFP) pathway and the indirect liquefaction (IDL) pathway for the conversion of syngas to high-octane gasoline (HOG). Cost reduction through refinery integration, development of valuable co-products, and other options are being identified for future research to help reduce the modeled MFSP to $2.50/GGE by 2030. Additional priorities anticipated in the future, such as the use of renewable electricity for liquid fuels and products, and emphasis on waste utilization are also being explored in conjunction with research on catalytic utilization of syngas and other gases (including CO2). Industry-relevant parameters are given deliberate attention as part of the work done under this project to help answer questions important for future commercialization and address associated risks.
29 ENERGY PLANNING, POLICY, AND ECONOMY↗
Thermal feedback coupling in a transient Monte Carlo high-order/low-order scheme
Transient simulations of nuclear systems face the computational challenge of resolving both space and time during reactivity changes. A common strategy for tackling this issue is to split the neutron flux into shape and amplitude functions. This split can be solved with high- order/low-order methods. While this multi-fidelity approach has traditionally been reserved for deterministic methods, it is also possible to implement in Monte Carlo as an efficient alternative to Dynamic Monte Carlo. This work implements the frequency transform method with thermal feedback in high-order/low-order Monte Carlo by blending static coupling methods such as single-batch Monte Carlo, with a simple thermal-fluids calculation. While previous work focused solely on prescribed transients, the addition of time-dependent thermal-fluids allows transients to be self-propagating. Tests were run in a fluids-initiated transient to showcase the basic functionalities of this methodology. Preliminary results behave as expected, paving the way for studying more sophisticated thermal-fluids coupling methods. (authors)
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Transient MOC with frequency transform and DSA on unstructured mesh
We present an implementation of the transient method of characteristics (MOC) with isotropic time derivatives, accelerated by diffusion synthetic acceleration (DSA). The fully implicit frequency transform method is used to solve the transient problem with analytic precursor integration. The code works on meshes composed of almost any of the commonly used non-curvilinear finite element types, and can handle the deformation of geometry in time-dependent transport calculations. We present results of a continuous Fourier analysis for the transient multigroup DSA problem, and representative benchmarking results are presented for the C5G7-TD benchmark in 2D showing reasonable performance and agreement compared to other codes. (authors)
22 GENERAL STUDIES OF NUCLEAR REACTORS↗
Transient convective delayed neutron precursor of {sup 235}U for the molten salt reactor experiment
This work investigates the delayed neutron precursors from the fission of {sup 235}U in molten salt reactors. The six delayed neutron groups predict the spatial distribution of the advective transport for radioactive mass transfer in the molten-salt reactor experiment (MSRE). The Mole code, using the framework of the Multiphysics Object-Oriented Simulation Environment (MOOSE), was used to compute the delayed neutron concentration of the whole system. This approach can be used for analysis of reactivity in circulating conditions. The effect of delayed neutron concentrations on the static or dynamic behavior of the system in MSRE was analyzed. Using this approach, the pattern and validation of delayed neutron for a range of parameters of interest are reported and discussed herein. (authors)
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Importance of 3-D S{sub N} depletion in non-proliferation using BSOLVE
We present a single Pressurized Water Reactor (PWR) 3-D fuel rod design for depletion analysis using BSOLVE, our newly developed Runge-Kutta-Fehlberg based depletion code. BSOLVE is coupled with the deterministic 3-D S{sub N} particle transport code, PENTRAN, applied here with a 4-neutron energy group comparison to Continuous Energy (C/E) SERPENT2 Monte Carlo results. Differences are expected, as PENTRAN+BSOLVE retains full (multi-group) energy information for reactions, nuclide specific fission contributions, and energy dependent fission yields, using the latest available ENDF-BVIII data, important to retain accurate burned fuel inventories; SERPENT2 collapses burnup reactions to a single energy value. For depletion times up to ∼ 700 days and typical PWR power densities, relative differences between multigroup 3-D S{sub N} with full energy data and Monte Carlo one group burnup for trans-uranium nuclide concentrations and fission products are up to ∼20%. System eigenvalues are consistent, but with differences early and late in the cycle attributed to multigroup vs. C/E Monte Carlo cross sections. This work highlights the importance of low variance transport driven burnup for non-proliferation concerns, since plutonium quality varies significantly along axial lengths, and is more challenging to converge using Monte Carlo; details of depletion steps with spatial/zone dependent plutonium quality are provided. (authors)
21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
Physics-Informed Recurrent Neural Networks to Predict Reactor Operations of the AGN-201 Nuclear Reactor
4 page paper submitted to ANS Student conference. Summary of paper similar to the following abstract: The ability to predict how a reactor will operate, understand when anomalous conditions arise, and ensure a reactor is being operated as expected is crucial for deploying new nuclear facilities. Digital twins serve as a unique solution to recognizing reactor behavior; however, they require data to be useful. For next-generation reactors, this data may not currently be available. To explore how synthetic physics-informed reactor data can be used to predict reactor operations, a recurrent neural network was implemented for the Idaho State University AGN-201 digital twin. The goal of this work is to determine how synthetic data can be used to train a recurrent neural network model for predicting the reactor power of the AGN-201. The recurrent neural network was validated using both synthetic and real operational data. We envision this approach will help bridge the gap between the virtual and physical sides of a digital twin, where reactor physics models based on as-built data can be corrected for actual operating parameters to ensure the virtual model mirrors reality.
98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
High-Temperature Fe-Based Fischer–Tropsch Synthesis: Experimentally Validated Kinetic Models Implemented at Pellet and Reactor Scales
Pellet- and reactor-scale models for Fischer–Tropsch synthesis (FTS) with a Fe–K/silica catalyst were developed to investigate the sensitivity of the hydrocarbon products and carbon dioxide selectivity to process conditions and feed composition at high temperature (350–400 °C), moderate pressure (1–10 bar), and a range of H 2 /CO ratios (3–1). The major objective of this paper is to develop, validate, and evaluate a high-temperature FTS model that is then used to assess the feasibility of process integration with syngas production. Since there is limited kinetic data available, in literature at these conditions, bench-scale reactor tests were conducted to obtain operational data for parameter fitting of kinetic expressions used in the model. This resulting kinetic model demonstrated agreement with the experimental data with an R 2 of 0.97 to the testing data set and, thus, was feasible to apply at pellet and reactor scales. Here, multiple pellet sizes were modeled to detail the role of transport limitations as the sphere’s diameter approached and exceeded 1 mm. Application of the reactor model indicated that hydrocarbon selectivity depended strongly on temperature, whereas the ratio of olefin to paraffin products decreased with increasing temperature, pressure, and H 2 /CO ratio. Product selectivity was not sensitive to the conversion of carbon monoxide. Furthermore, the roles of the pressure and H 2 /CO ratio were closely coupled. At a H 2 /CO ratio of 3, only slight variations in selectivity occurred over a pressure range of 1–20 bar, whereas at a ratio of 1, selectivity could vary by as much as 30% over the same pressure range. At pressures below 5 bar and temperatures above 350 °C, minimal selectivity to heavy hydrocarbons (C 12+ ) is obtained, and selectivity to midrange products (C 5–11 ) rapidly declined as pressure dropped below 5 bar, which indicated that an operational pressure of at least 5 bar is needed to achieve reasonable yields in this temperature range. These results, while tentative, provide guidelines for further experimentation and evaluation of integrated FTS processes.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Accelerated kinetic model for global macro stability studies of high-beta fusion reactors
The field reversed configuration (FRC), such as studied in the C-2W experiment at TAE Technologies, is an attractive candidate for realizing a nuclear fusion reactor. In an FRC, kinetic ion effects play the majority role in macroscopic stability, which allows global stability studies to make use of fluid-kinetic hybrid (also referred to as Ohm's law) models wherein ions are treated kinetically while electrons are treated as a fluid. The development and validation of such a hybrid particle-in-cell algorithm in the Exascale Computing Project code WarpX are reported here. Implementation of this model in the WarpX framework benefits from the numerical efficiency of WarpX as well as its scalability on large HPC systems and portability to different architectures. Performance benchmarks of the new algorithm for large, 3-dimensional, full device simulations from the Perlmutter supercomputer are presented. Results of a series of FRC simulations are discussed in which the impact of two-fluid effects on the tilt-mode growth rate was studied. It was observed that, in agreement with previous Hall-MHD studies, two-fluid effects have a stabilizing impact on the tilt mode.
Physics↗
Design improvements for a recirculating reactor: Enhanced temperature measurement and sample-isolated reactivity in steady-state kinetic studies
Building upon a previous recirculating reactor design [S.A. Tenney, K. Xie, J.R. Monnier, A. Rodriguez, R.P. Galhenage, S. Audrey, D.A. Chen, Rev. Sci. Instrum. 84, 104101 (2013)], we present significant improvements that address key limitations in steady-state kinetic measurements for heterogeneous catalysis. The enhanced reactor design features direct sample heating with a focused IR lamp and temperature measurement in direct contact with the sample, enabling more accurate temperature control and improved kinetic analysis. A critical advancement is the isolation of sample reactivity from reactor wall contributions, ensuring that only the sample contributes to measured reaction rates. This was a limitation in earlier designs where the entire reactor contributed to the observed reactivity. The system incorporates a bypass flow cell for direct comparison with powder catalysts under identical conditions using a standard plug-flow reactor configuration. We demonstrate these capabilities through CO oxidation experiments on Pt(111) single crystals and graphene-passivated Pt(111), highlighting the system's ability to differentiate catalytic activity in model systems and directly compare them with high surface area powder catalysts. This reactor is particularly suited for thin films and low surface area catalysts that are not effectively evaluated in traditional flow reactors, especially for samples with low numbers of active sites or slow reaction rates.
36 MATERIALS SCIENCE↗
Multiscale Modeling of a Direct Nonoxidative Methane Dehydroaromatization Reactor with a Validated Model for Catalyst Deactivation
Due to the recent boom in shale gas production, aromatics production using direct nonoxidative methane dehydroaromatization (DHA) is being investigated extensively. However, due to rapid coke formation, catalysts in the nonoxidative methane DHA reactors get deactivated, which is one of the critical issues for the commercial success of the methane DHA process. In this paper, a model for catalyst deactivation is developed. Rate models for other DHA reactions are developed by considering the decrease in the catalyst activity with time. Due to the very fast coke formation rate on the fresh catalyst, there is coke formation immediately upon the introduction of the feed. Therefore, an algorithm is developed for estimation of the initial state of the reactor and the kinetic parameters by coupling an iterative direct substitution approach with an optimization approach. Transient experimental data from an in-house reactor are first reconciled and then used for developing the kinetic model including the coke formation model. Using the rate model, a dynamic, heterogeneous, multiscale reactor model with embedded heating is developed. Here, the model couples the catalyst pellet level model with a reactor level model. Impacts of temperature, L/D ratio, and scheduling of reactors on variability in conversion and yield with time are studied.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Role of surface diffusion in formation of unique reactivity for graphite oxidation: Time-resolved measurements in a pulsed diffusion reactor
Quantification of oxidation kinetics is essential to develop graphitic materials for diverse applications: from refractories found in gas-cooled nuclear reactors to catalysts needed for chemical manufacturing. In this work, using well-defined highly oriented pyrolytic graphite, low-pressure isotopic transient experiments combined with controlled annealing periods, we resolve the role of surface diffusion and quantify oxidation kinetics with nanomole-precision. We observe an unexpected increase in reactivity following annealing which is explained by the role of surface diffusion increasing the probability for trapping mobile oxygen at more reactive edge sites. Here, the locus of adsorption and spillover to the basal plane is distinct from the trapping location creating a more active oxygen species. Isotopic products reflect the population dynamics of oxygen added at the edge and surface diffusion that relocates basal plane oxygen to more reactive edge sites. Since this process proceeds in parallel with direct oxidation reactions, it is not likely to be observed using steady-state or conventional ‘bulk’ characterization techniques. Our unique time-resolved non-equilibrium measurement in a well-defined transport regime, enables observation of three distinct behaviors: short-term deactivation due to the balance of rates in oxygen supply/product formation, reactivity increases due to surface diffusion and longer-term reactivity increase with oxygen accumulation.
36 MATERIALS SCIENCE↗
On the novel 3-D neutron transport kinetic tRAPID algorithm and its validation
The Real-time Analysis for Particle-transport and In-situ Detection (RAPID) Code System, based on the Multi-stage Response-function Transport (MRT) methodology, allows for real-time simulation of nuclear systems based on 3-D continuous-energy particle transport. RAPID's steady-state (criticality) neutron transport algorithm is based on the Fission Matrix (FM) method, and has been extensively verified and validated against computational benchmarks and experiments. This paper introduces the novel 3-D time-dependent transport algorithm that has been implemented into the code, tRAPID, and its validation using the JSI TRIGA Mark-II reactor. tRAPID accurately and efficiently calculates neutron kinetics parameters (such as β{sub eff}, l{sub eff} , Λ, α{sub Rossi}) and 3-D time-dependent neutron fission source distribution and neutron importances for both prompt and delayed neutrons. tRAPID is used to simulate a rod insertion experiment performed at the JSI TRIGA Mark-II reactor, during which signals from four fission chambers at four different locations in the core were collected. The results demonstrate how tRAPID is capable of calculating detailed and accurate results with only a minimal use computational resources and time.
22 GENERAL STUDIES OF NUCLEAR REACTORS↗