Analytical continuation for calculating multigroup-multizone eigenvalue problems using Lie series.
Critical volume of cylindrical reactors calculated using n-group diffusion theory
SEARCH · Search NASA
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.
Critical volume of cylindrical reactors calculated using n-group diffusion theory
The prediction of boiling incipience is a critical issue for a reactor. Vapor lock causes operating instability and undesirable temperature rise, especially for small channels. This study uses the Eulerian multiphase flow boiling model to predict the onset of nucleate boiling in a millimeter scale rectangular channel. Two wall boiling models, namely the RPI (Rensselaer Polytechnic Institute) and the non-equilibrium sub cooled model are applied. Computational fluid dynamics is employed for this study to investigate the impact of operating conditions like liquid flow rate, operating pressure, heat flux, velocity profile, turbulence model, and inlet liquid temperature. Here, the effect of fundamental boiling parameters on the surface wall temperature and vapor volume fraction is also studied. The parameters under consideration are the bubble departure diameter, bubble departure frequency, nucleate site density, quenching time period, and interface heat transfer coefficients. A series of CFD calculations is conducted by varying the considered variables systematically in a wide range of flow conditions covering laminar, transition, and turbulent flows. The sub-component heat fluxes like the liquid convective, the vapor convective, the evaporative, and the quenching terms are monitored to infer the boiling dynamics in the transition region from the single-phase flow to the nucleate boiling zone. It is found that the CFD approach to detecting the boiling incipience point agrees reasonably with available experimental data. However, its limitations, like inaccuracy for lower flow rates and early transition to boiling flow, are also noticed due to the nature of the considered RPI model.
Sandia National Laboratories (SNL) personnel operate a low power research reactor (the Annular Core Research Reactor, or ACRR), and a zero-power critical experiment assembly referred to as CX. In accordance with 10 CFR 830, Subpart B, Appendix A, the acceptable methodology for developing a Documented Safety Analysis (DSA) for DOE nuclear reactors is the Nuclear Regulatory Commission’s (NRC’s) Regulatory Guide 1.70 (RG 1.70). RG 1.70 does not address certain areas required by 10 CFR 830 and expected by DOE (e.g., full facility hazard analysis).Thus, the current DSAs for SNL’s reactor nuclear facilities are based on RG 1.70, but also of necessity supplemented by DOE-STD-3009-94 methods. SNL personnel, in consultation with the National Nuclear Security Administration (NNSA) Sandia Field Office (SFO), have concluded that an alternate methodology is preferred to RG 1.70. The details of the proposal, and the reasons motivating its development, are discussed in the order described below. The proposed alternate methodology will be applicable to the ACRR and the CX (i.e., it will be applicable to nuclear facilities in which a reactor and/or a critical assembly will be operated).
The first edition of the current American National Standards Institute (ANSI)/American Nuclear Society (ANS) standard ANSI/ANS-8.1-2014, “Nuclear Criticality Safety in Operations with Fissionable Materials Outside Reactors,” was published in 1964 as ASA N6.1-1964, “Safety Standard for Operations with Fissionable Materials Outside Reactors.” In 1969, that standard was revised as ANSI N16.1-1969, “Nuclear Criticality Safety in Operations with Fissionable Materials Outside Reactors.” ANSI N16.1-1969 includes a variety of subcritical limits (SCLs) for uniform aqueous solutions and metals containing fissile nuclides of 233 U, 235 U, and 239 Pu. Furthermore, SCLs are also included for uranium-water lattices. In the 1983 version of ANSI/ANS-8.1 (a revision of ANSI N16.1-1975), the suite of SCLs in the standard was expanded to include 235 U enrichment limits for homogeneous uranium-water mixtures and dry/damp oxides, uniform aqueous solutions of low-enriched uranium, and uniform aqueous mixtures of Pu(NO 3 ) 4 containing 240Pu in addition to the SCLs included in ANSI N16.1-1969. The SCLs have changed little in subsequent revisions (ANSI/ANS-8.1-1998 and ANSI/ANS-8.1-2014). The ANSI/ANS-8.1-2014 standard is currently being revised to include new SCLs (uranium metal and compounds with enrichments up to 20 wt% 235 U) and possible updates to the current SCLs already in the standard, although these new/updated SCLs will not be available to the nuclear criticality safety (NCS) community for a number of years. The original bases for these SCLs were documented in papers in journals such as Nuclear Science and Engineering and in ANS meeting transactions; however, these bases are ambiguous enough that sites and regulators in the United States have not widely endorsed them for safety purposes. Overall, the purpose of this paper is to present the results of a comparison study for the SCLs in the ANSI/ANS-8.1-2014 standard using modern codes (SCALE and MCNP) and cross sections (ENDF/B-VIII.0) to provide some assurance about their quality (bias and bias uncertainty) for use in NCS applications and for consideration by the ANS-8.1 Working Group as a reference for future revisions.
Heat transfer between cladding and coolant during transient scenarios remains a critical area of uncertainty in understanding nuclear reactor safety. To advance the understanding of transient and accident scenarios involving critical heat flux (CHF), an in-pile experiment for the Transient Reactor Test (TREAT) facility at Idaho National Laboratory was developed. The experiment named CHF-Static Environment Rodlet Transient Test Apparatus (CHF-SERTTA), consists of a hollow borated stainless steel heater rod submerged in a static water pool which is heated via the (n, a) reaction in boron-10. This paper presents a novel inverse heat transfer method to determine CHF by using the optimization and uncertainty software DAKOTA to calibrate a RELAP5-3D model of CHF-SERTTA to temperature measurements obtained from a thermocouple welded to the surface of the rod.
Critical experiments, sometimes referred to zero-power reactors, are crucial tools in developing and validating nuclear physics predictions and thus an indispensable capability to ensure criticality safety during all parts of the nuclear fuel cycle and in reducing uncertainties in reactor physics predictions toward optimizing nuclear energy production. The System Physics Advanced Reactor Facility (SPARC) project was recently initiated to enable large-scale criticality experiments using a horizontal split table machine well suited toward solid core materials systems (fuel, moderator, reflector). The facility selected for this mission was a former pool-type research reactor building and thus also well suited toward a second critical experiment capability able to house full-size light-water reactor (LWR) fuel bundles. A conceptual design study was undertaken to review past water tank critical experiments used for LWR physics experiments and to develop an early engineering design for a new critical assembly tank (CAT). The work described here shows that a relatively simple CAT concept can be constructed and deployed in the SPARC facility to meet the urgent demands for new critical experiments on advanced LWR fuel bundles designs. The SPARC facility layout is conducive to the receipt and upending of LWR fuel bundles using existing containers and equipment from the LWR industry. The facility’s overhead crane can then be used to handle fuel bundles and place them in a vertical storage rack or in the CAT for critical experiments, both of which fit within the building’s “open basement” alongside other equipment planned for SPARC and the horizontal split table. A slightly lower area in the basement can serve as a large drain tank so that fail-safe valves drain the CAT reactor tank for safe shutdown. Neutronic configurations were determined where a 3 × 3 array of fuel bundles can be surrounded by full-length “loose rods” to adjust reactivity so that critical is achieved when the bundles are fully submerged. Viable configurations were determined for both pressurized- and boiling-water-reactor-type fuel bundles. This design concept was used to develop an early planning basis for establishing the CAT capability alongside the otherwise planned SPARC project in order to help streamline the process. Recent presidential executive orders have highlighted the need to achieve power uprates in LWR plants and the CAT capability will be a crucial element of these initiatives. Based on the work described herein, it is recommended that an earnest and timely project begin in order to establish this urgently needed capability.
Deep generative modeling provides a powerful pathway to overcome data scarcity in energy-related applications where experimental data are often limited. By learning the underlying probability distribution of the training dataset, deep generative models, such as the diffusion model, can generate high-fidelity synthetic samples that statistically resemble the training data. Such synthetic data generation can significantly enrich the size and diversity of the available training data, and more importantly, improve the robustness of downstream machine learning models in predictive tasks. The objective of this paper is to investigate the effectiveness of diffusion models for overcoming data scarcity in nuclear energy applications. By leveraging a public dataset on critical heat flux which covers a wide range of commercial nuclear reactor operational conditions, we developed a diffusion model that can generate an arbitrary amount of synthetic samples. Since a vanilla diffusion model can only generate samples randomly, we also developed a conditional diffusion model capable of generating targeted critical heat flux data under user-specified thermal-hydraulic conditions. The performance of the diffusion model was evaluated based on its ability to capture empirical feature distributions and pair-wise correlations, as well as to maintain physical consistency. The results showed that both the diffusion model and conditional diffusion model can successfully generate realistic and physics-consistent critical heat flux data. Furthermore, uncertainty quantification results demonstrate that the conditional diffusion model is highly effective in augmenting critical heat flux data while maintaining acceptable levels of uncertainty.
Given the existential climate crisis faced by mankind and the world, the lifetime and sustainability of nuclear reactors as a carbon-free source of renewable energy depend on the susceptibility of their structural components to environmental degradation. In particular, critical components for light water reactors (LWRs) evolve over decades in service, losing ductility and toughness due to thermal and irradiation aging. Techniques to monitor their health cannot be easily applied in the field due to their destructive, expensive, or immobile nature. Thus, non-destructive evaluation (NDE) methods are sought to monitor and evaluate the health of major LWR components such as core barrels, steam generator tubes, or primary coolant pipes and are often required by policy, such as NRC policy #10-CFR-50.65. In this work we demonstrate the use of gigahertz, non-contact ultrasonics to gauge the state of cast austenitic stainless steels (CASS), used in some of the largest components in LWR primary systems. We do so by linking changes in their surface acoustic wave (SAW) characteristics using transient grating spectroscopy (TGS) to transmission electron microscopy (TEM)-verified evidence of spinodal decomposition and G-phase precipitation. In this thesis, thermal aging is shown to induce SAW peak splitting in spinodally decomposed CASS alloys, correlated strongly with lowered toughness and decreased ductility. Furthermore, statistical testing on the number of SAW peak splits observed show that the second SAW peak significantly appears more frequently and is significantly different in frequency in comparison to counts and frequencies measured in unaged specimens. The ability of this technique to non-destructively detect microstructural degradation at a distance in a predictive manner in the case of CASS motivates extending gigahertz ultrasonics to detect other LWR material degradation modes as an in-vessel inspection technique, such as reactor pressure vessel (RPV) embrittlement. This allows for the greater use of NDE techniques for confident monitoring of LWR structural material health to 80 years and beyond, saving costs by minimizing structural replacements until needed and maximizing energy production by preventing early decommission until necessary.
Mass accountancy in a molten salt reactor (MSR) system is the ability to track and quantify the radioactive species throughout all portions of the reactor, not just the core. This is critical for the design, analysis, and regulation of these novel reactors. This report discusses a chemistry-based approach for understanding the chemical state, physical characteristics, and time-dependent location of a species throughout an MSR system. This report does not represent a completed capability; rather, it indicates that specific activities are needed (i.e., chemistry coupling and benchmark creation) and require sustained support to enable MSR development and deployment.
The Transient Reactor Test Facility (TREAT) at Idaho National Laboratory (INL) serves a vital role in nuclear fuel safety research, enabling transient experiments that simulate reactivity excursions and accident scenarios. Central to these operations is the transient control rod drive system (TCRDS), which drives rapid motion of the transient control rods such that TREAT can simulate rapid power changes typical of reactor accidents. The reliability and performance of this system are critical for protecting both fuel specimens and reactor infrastructure. This study presents the initial phase of a two-year investigation into the dynamics and reliability of the TREAT hydraulic TCRDS. Conducted in collaboration with INL, the research employs a combined computational and experimental approach to analyze the system's response time, pressure transients, and potential failure modes. Emphasis is placed on understanding how fluid characteristics influence the TCRDS’s ability to achieve both rapid power changes and mechanical stability. The TRDS and the skid that powers it will be analyzed throughout this investigation. Computational modeling using computational fluid dynamics (CFD) will simulate the hydraulic response under varying conditions. In parallel, experimental testing planned at INL will validate these models and capture key performance metrics. This paper outlines the system design, analytical framework, and modeling strategies that form the foundation for later testing. Ultimately, this work aims to support improvements to the TCRDS’s design and reliability, contributing to the broader goal of enhancing nuclear fuel safety and sustaining TREAT’s mission as a premier nuclear fuel test facility.
Under the extreme environments of light water reactors, the cladding materials demonstrate evolutionary transformations in terms of material compositions, surface morphologies, and thermal-physical properties. This study investigates how the temperature overshooting of Critical Heat Flux (CHF) affects the re-occurrence of CHF on the post-CHF surfaces. Experimental results reveal that the CHF difference between as-received and post-CHF surfaces can be significant, which depends on materials. CHF is enhanced on post-CHF Stainless Steel (SS 316) specimens while it deteriorates on post-CHF Inconel-600 specimens. In comparison with the as-received surfaces, the post-CHF surface of SS 316 is more hydrophobic, while Inconel 600 is more hydrophilic. The scanning electron microscope micrographs of post-CHF surfaces show that surface morphologies that are structured by oxide crystals are formed on SS 316 samples but not on Inconel 600 samples. The physical rationale behind the CHF difference could be boiling heat transfer conjugated with thermal-physical proprieties of investigated surface materials and surface morphologies features. In conclusion, the results imply that it is indispensable to bookkeep the evolutionary changes of cladding surface materials, especially for cladding candidates of Accident Tolerant Fuel since operation histories can leave footprints on cladding surfaces, thus resulting in CHF re-occurrence differences.
Accurately predicting the criticality of an experiment before interacting with the experimental components is very important for criticality safety. Radiation transport software can be utilized to calculate the effective neutron multiplication factor of a nuclear system. Because of the integral nature of the effective neutron multiplication factor, the value calculated contains various sources of nuclear-data induced uncertainty. The sensitivity analysis and data assimilation technique presented in this paper exhibit one possible method of identifying and reducing the effective neutron multiplication factor nuclear-data induced uncertainty. The results presented in this work show that it is possible to use relative sensitivity coefficients of the prompt neutron decay constant and the effective delayed neutron fraction to 239 Pu nuclear data to reduce nuclear-data induced uncertainties in the effective neutron multiplication factor. This work has been utilized by members of the Los Alamos National Laboratory project EUCLID (Experiments Underpinned by Computational Learning for Improvements in Nuclear Data) for optimally designing a new experiment, which will be used to reduce compensating errors in 239 Pu nuclear data.
Microbial electrolysis cells (MECs) are appealing for recovering the chemical energy contained in domestic and industrial liquid wastes as hydrogen gas. Despite several years of research in the field, there is still a lack of critical analysis of how the reactor architecture dictates the electrochemical performance of the cell. In this study, internal resistance and onset voltage from the electrode potential slope analysis (EPS) were used in combination with current density, hydrogen production rate, reactor packing density, electrode spacing, membrane type and composition from 23 different studies to identify the reactor design parameters that primarily govern electrochemical performance of MECs. Using anion exchange membranes resulted in smaller internal resistances (AEM R int = 41± 40 mΩ m 2 ) and larger current density (18 ± 14 A m −2 ) compared to single chamber reactors (SC R int = 68 ± 58 mΩ m 2 ; 22 ± 16 A m −2 ) or MECs with cation exchange membranes (CEM R int = 376 ± 280 mΩ m 2 ; 3.0 ± 2.1 A m −2 ). Higher electrochemical performance for AEM- and SC-MECs translated in larger hydrogen gas production rates (0.122 mL H 2 C −1 for AEM vs 0.117 mL H 2 C −1 for SC), but only when inhibitors against hydrogen scavengers were added in single chamber systems (0.080 mL H 2 C −1 for SC without inhibitors). Following membrane type and composition, maintaining a small electrode spacing was the most critical parameter to improve MEC performance, indicating that the low conductivity of the media primarily limit performance by increasing ohmic resistance. Here, reactor volume and electrode surface area negatively correlated with internal resistance and current density, indicating that better performance of scaled-up reactors can likely be obtained by stacking multiple smaller units rather than just increasing reactor size. Although challenges remain in the implementation of MECs for hydrogen production from liquid wastes, advances in electrochemical engineering of the reactors can facilitate scale up and performance prediction at scale.
A pseudo real-time interactive reactor simulator for the Lady Godiva reactor has been developed using the GNU Octave software package. The simulator tracks the time-dependent power and fuel temperature behavior of the Lady Godiva reactor during sub-prompt critical transients. Users interact with the simulator by controlling reactivity input. This input occurs through users repositioning a control rod as desired during simulator operation. Various reactor parameters are displayed graphically in near “real time,” facilitating user assessment and feedback input. The reactor model consists of a lumped-fuel core, a single control rod, and a simple heat transfer mechanism that allows heat to be rejected from the core. The goal of the simulator is to demonstrate in real-time the system response to reactivity perturbations, the effect of reactivity feedback, and the effect of temperature-dependent equilibrium conditions. The dynamic behavior of the reactor is primarily driven by the relative characteristic time-constants between the delayed neutrons and the overall heat transfer parameter. The simulator will allow users to gain experience performing standard reactor operations, such as establishing steady-state at different fuel temperatures and making power level changes
The Advanced Test Reactor (ATR) creates a unique, high neutron flux materials testing environment that subsequently causes neutron embrittlement damage to its beryllium reflector. The reflector is eventually required to be replaced during a reactor Outage. Concurrently, a piece-for-piece core replacement is conducted to maintain reactor functionality called a core internals changeout (CIC). Once changeout is complete, ATR Reactor Engineering utilizes previous CIC data, notably 1977-1994 CIC data, to conduct predictive criticality analysis on the new, unirradiated core. However, CIC data has been inconsistently reported over the years making it difficult, and in some cases impossible, to utilize. These inconsistencies arise in various forms of not using correct units, inconsistent diagram usage, and not portraying the same kind of data similarly across the reports. To work around this inconsistency issue, Reactor Engineering is implementing the use of several new neutronic modeling codes that utilize modern computation speeds and methods to model the new core in different criticality environments. Further validation is then conducted by nuclear testing and bringing the reactor critical and conducting various irradiation tests. Upon completion and certification that the reactor passed all testing requirements, ATR can once again re-start normal operations.
The validation of the Ants nodal neutronics code for VVER applications is started by modelling the zero power physics tests for the initial core of the Khmelnitsky 2 nuclear power plant as described in the X2 benchmark using the Serpent-Ants two step neutronics calculation chain. The Ants prediction compare favorably against an earlier continuous energy Monte-Carlo reference solution as well as the measured data, except for discrepancies in the SCRAM worth between predicted and measured values. Similar discrepancies have been previously reported for VVER-1000 reactors, and the Ants predicted SCRAM worths match well with the Serpent results. Additionally, the assembly power distribution and axial power distribution predicted by Ants for the critical HZP state of the reactor are also compared to a reference Serpent solution with very good results. Lastly, an investigation into the effects of the few-group structure used in the nodal calculations and the use or lack of leakage correction for group constants on the results shows that the best accuracy is reached with seven or eight energy groups without leakage correction, although a reasonable accuracy can also be obtained with two or three energy groups using fundamental mode leakage correction. (authors)
An experimental investigation was conducted using the United Technologies Research Center (UTRC) 80 kW and 1.2 MW RF induction heater systems to aid in developing the technology necessary for designing a self-critical fissioning uranium plasma core reactor (PCR). A nonfissioning, steady-state RF-heated argon plasma seeded with pure uranium hexafluoride (UF6) was used. An overall objective was to achieve maximum confinement of uranium vapor within the plasma while simultaneously minimizing the uranium compound wall deposition. Exploratory tests were conducted using the 80 kW RF induction heater with the test chamber at approximately atmospheric pressure and discharge power levels on the order of 10 kW. Four different test chamber flow configurations were tested to permit selection of the configuration offering the best confinement characteristics for subsequent tests at higher pressure and power in the 1.2 MW RF induction heater facility.
Extensive efforts have been carried out at ANL for the verification and validation of the Argonne Reactor Codes (ARC) software package currently used for the design of Versatile Test Reactor (VTR). The ARC software package consists of steady state neutronics and thermal hydraulics modeling capabilities which are being used by the VTR program to develop most of the VTR reactor design details which will be part of the licensing application. It is anticipated that this software will continue to be used for the design work and for initial operations although additional software may be introduced at a later time. The validation work was focused primarily on obtaining validation data consistent with VTR and usable for the ARC software. Because no critical facilities or operating fast spectrum reactors are available to do experiments for the VTR, the next best option is to identify historical experimental data that can be used as validation data. Early on in VTR, the ZPPR-15 set of experiments was identified as good validation data because of 1) the availability and quality of the data, 2) existing staff that are already familiar with the experimental machine and measurements, 3) most of the ZPPR-15 loadings of interest have already been processed into ARC models, and 4) a full uncertainty quantification has already been done for several loadings of ZPPR-15. The FFTF startup and operations data was identified as the most consistent reactor type that has validation data usable for VTR. Finally, the EBR-II fuel depletion measurements were identified as the best available validation data for VTR. It is important to note that both the FFTF and EBR-II reactors typically come with higher uncertainties than the ZPPR. In the frame of the discussed verification and validation efforts, the present document discusses the analysis of selected FFTF measurements included in the benchmark specifications of the International Reactor Physics Experiment (IRPhE) handbook. The FFTF reactor core configurations from the benchmark specification are presented in Section 2. The analysis is performed with the use of the ARC code suite available at ANL for fast reactor studies and is discussed in Section 3. The reactor parameters from the benchmark include criticality, neutron spectra, effective delayed neutron spectra, control rod worth, isothermal temperature coefficient and low energy gamma-ray spectra. The calculated values and the comparison with the experimental data is discussed in Sections 4 to 9 for each considered reactor parameter. Finally, conclusions are presented in Section 10.