Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CONTROL ROD WORTHS”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Transition core modeling for extended enrichment, accident tolerant fuel using PARCS/Polaris

Current plans and efforts from reactor operators and vendors to include extended enrichment (EE) and accident-tolerant fuel (ATF) in current reactor fleets motivate the study of these changes in reactor physics analysis. This work uses the US Nuclear Regulatory Commission core simulator PARCS to do the core calculation and SCALE Polaris lattice physics code to generate the required homogenized, few-group constants. The lattice model used is based on the GE-14 10x10 assemblies with UO{sub 2} fuels. Both nominal core and transition core are studied in this work, and both cores use GE-14 10x10 assemblies with UO{sub 2} fuel. The nominal core uses regular UO{sub 2} fuel with a maximum enrichment of 5 wt % and ZIRC-2 cladding. The ATF transition core uses FeCrAl cladding and regular UO{sub 2} fuel, while the EE-ATF transition core uses FeCrAl cladding with UO{sub 2} fuel with 8 max wt % enrichment. The accuracy presented in the colorset models verified the capability of the PARCS/Polaris procedures for the transition core analysis. For the whole core calculation, the ATF and EE-ATF transition core models were made in addition to the nominal core model. The core parameters to study are the core power distribution and power peaking factor, doppler temperature coefficients, and control rod worth at cold zero power and hot full power. Comparing the core parameters of the transition cores with the nominal core in PARCS, the results suggest that there is no unexpected outcome for the implementation of ATF and EE fuels. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

Reactor physics benchmark experiments at the JSI TRIGA MARK II reactor - Current status and future outlook

Full text of publication follows. With the development of new high-fidelity computational methods, improvement of nuclear data, and multiphysics modelling, there is an increased need for benchmark experiments to experimentally validate the models, methods and input data. Many of the nuclear facilities designed to perform reactor physics benchmark experiments have been shut down. Therefore, research reactors offer a great opportunity for benchmark experiments, if they are well designed and performed with great care and accuracy. In this presentation we provide an overview of the past and ongoing activities related to benchmark experiments at the Jozef Stefan Institute TRIGA Mark II research reactor. The following experiments have been performed: criticality with fresh fuel, {sup 197}Au(n,γ) and {sup 27}Al(n,α) reaction rates in irradiation channels, absolute and relative {sup 197}Au(n,γ), {sup 235}U(n,f) and {sup 238}U(n,f) reaction rates in the core, burnup, kinetic parameters, control rod worth, isothermal reactivity coefficient, self-shielding, slow and fast (pulse) transients, nuclear heating, delayed and prompt gamma ray production, temperature profiles for multi-physics. Since the existing fleet of research reactors is ageing very rapidly and new experiments are needed, new research reactors should be designed and built to meet the needs of future advanced reactors, education and training, and other technologies in the coming years. We will review planned activities at the JSI TRIGA reactors and plans for the new research reactor in Slovenia. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Whole core analysis of OECD/NEA TVA-WB1 benchmark exercise 1 using continuous-energy Monte Carlo code MVP with JENDL-4.0u and ENDF/B-VII.0

The TVA WB1 benchmark problem was precisely modeled in accordance with the geometry and material specifications, and benchmark exercise 1 was solved with the MVP code and the nuclear data libraries based on JENDL-4.0u (JENDL-4.0 updated file) and ENDF/B-VII.0. The criticality, the control rod bank worth, the isothermal temperature coefficient (ITC), the integral rod worth of bank D, and differential soluble boron worth (DBW) were obtained. The criticality and the control rod bank worth were compared with the measured data and published calculation results with other codes, Serpent and KENO. The criticality and the control rod bank worth calculated by MVP and the libraries agree with the measured and the published calculation results within the criteria specified in the ANSI/ANS standard. From these results, the calculation model for the MVP and the libraries was verified. The other parameters such as ITC, the integral rod worth, and DBW were also compared between JENDL-4.0u and ENDF/B-VII.0. The comparison of the results shows JENDL-4.0u and ENDF/B-VII.0 give comparable results within the standard deviation, the difference between libraries is small. (author)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

High-performance and high-fidelity Monte Carlo solutions to the BEAVRS benchmark

The BEAVRS (Benchmark for Evaluation and Validation of Reactor Simulation) benchmark is solved by PRAGMA, the GPU-based continuous energy Monte Carlo code. The resulting solutions are comprised of the detailed simulation results of two cycles, each of which consists of the zero power physics test (ZPPT) and the core depletion calculations. The ZPPT consists of characteristic parameters, such as critical boron concentration (CBC), control rod bank worth, isothermal temperature coefficients, and assembly-wise detector signal, which are compared with measured data provided by the benchmark administration. The core depletion calculations were performed for both the hot full power and the load follow modes, and the comparison was made with the measured or deduced CBCs and assembly-wise detector signals. In the load follow calculations, the operating power history was approximately applied to simulate the real operation as closely as possible. PRAGMA performed the various calculations with a tremendous number of histories ranging up to hundreds of millions per cycle, exploiting GPUs' massively parallel performance. The load follow run time was shorter than 16 hours on a single rack of computing nodes mounded with 24 gaming GPUs with a remarkable agreement with the measurements for most comparisons. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

MCS solutions for the TVA Watts Bar Unit 1 multi-physics depletion benchmark

The high-fidelity solution for exercises 1-3 of the recently developed TVA Watts Bar Unit 1 multi-physics and multi-cycle depletion benchmark using MCS was coupled with two different thermal-hydraulics (TH) codes: TH1D and COBRA-TF (CTF), and their solutions were compared. The solutions from MCS/CTF and MCS/TH1D for typical Pressurized Water Reactor (PWR) problems are quite close. And to evaluate MCS solutions for the benchmark, their results were also compared against measured data as well as publicly available high-fidelity solutions. The MCS criticality solutions at Hot Zero Power (HZP) deviate within 50 pcm against the reference solutions. While for the critical boron concentration (CBC) search calculations, the difference on the MCS calculated CBC is about 7 ppm, and the root-mean-squared (RMS) errors for the assembly power distribution and outlet coolant temperature are less than 0.5% and 1 Celsius degree respectively. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

CEFR simulation using diffusion code system RAST-F

In this study, the RAST-F code system, which is based on the two-step approach that couples a multi-group cross-section generation Monte-Carlo (MC) code and a multi-group nodal diffusion solver, was used for the neutronic simulation of the CEFR start-up experiments. The numerical solution of the RAST-F system was verified against the full core MC solution MCS at all control rods fully inserted and withdrawn states. The RAST-F solution of the selected experimental simulations was compared against the measurement data. A good agreement between RAST-F and MCS solutions was observed with less than 120 pcm discrepancies and 1.2% root-mean-square error in terms of k{sub eff} and power distribution, respectively. Meanwhile, the RAST-F result agreed well with the experimental data within two-sigma of experimental uncertainty. The good agreement of these results indicates that RAST-F can be used to neutronic steady-state simulations for small core-size SFR, which was challenged to deterministic code systems. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2023

A Predictive Transient Model of the TREAT-SIRIUS Experiments

To support experiments modeling and transient analysis of the NASA-sponsored SIRIUS experiments in the Transient Reactor Test Facility (TREAT), an innovative predictive transient model has been developed to simulate various transient experiments. The developed model is based on a new approach that utilizes steady-state Monte Carlo calculations along with a surrogate model based on polynomial regression to determine the core reactivity. The reactivity is supplied to an exact point kinetics model that determines the reactor power, and it is coupled to an adiabatic feedback model to determine the fuel average temperature. In this work, the newly developed model is introduced along with validation test results considering SIRIUS-1 experiments at different power levels. The initial test results show a very good agreement with experimental data at different power levels, which indicates that the model can be used to predict the power of other transient tests.

42 ENGINEERING↗
DOE OSTI2022

Simulation of CEFR neutronic start-up tests with FENNECS

This paper presents simulation results of selected Neutronic Start-up Tests of the China Experimental Fast Reactor (CEFR) obtained by the neutronics code FENNECS that have been performed within the frame of the IAEA Coordinated Research Program I31032. The Finite Element Neutronics code FENNECS is developed at GRS and solves the few-group steady-state and transient diffusion equation using a Galerkin-based finite element approach. Its main purpose is the safety assessment of Small Modular Reactors and Micro Reactors with complex geometry (e.g., rotating control drums) which gain increased interest internationally. Serpent has been applied to create reference models of the CEFR and cross-section libraries in 10 energy groups for FENNECS. Using these libraries, the following experiments have been simulated with FENNECS: net criticality, control rod integral and differential worth, void reactivity effect, subassembly exchange reactivity effects and reaction rate distribution. The obtained satisfactory agreements with the measurements represent a valuable contribution to the validation of FENNECS. (author)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Modeling of the TRIGA IPR-R1 research reactor with the Serpent2/RINNOVO Nodal core analysis package

The Serpent2/RINNOVO nodal core analysis code system, a dedicated tool for modeling research reactors, has the capability to accurately predict important core physics parameters involved in the safety of reactor operation, such as various reactivity coefficients, control rod and bank worths including the shutdown margin, power distributions as well as local neutron flux predictions at various core locations of high importance, e.g., at irradiation rigs. In this work, a hexagonal model of the unrodded initial core of the Brazilian IPR-R1 Mark I type TRIGA nuclear reactor has been created using the Serpent2/RINNOVO code system. The choice of employing a hexagonal core geometry representation has mainly been made to facilitate subsequent fuel shuffling operations and core follow calculations of this reactor. Numerical results in terms of the core eigenvalue and the assembly power distribution have then been compared against corresponding full core Serpent2 results to prove feasibility of using the Serpent2/RINNOVO code system for modeling small and highly heterogeneous TRIGA reactors. Overall, RINNOVO predicts the power distribution very accurately but the eigenvalue error still remains quite large. In light of being the very first evaluation of a TRIGA reactor with the RINNOVO nodal core simulator, these results are also considered to be very preliminary. Improved accuracy is expected by incorporating a proper methodology in the current code system for computing discontinuity factors for multi-assembly configurations. Furthermore, improved results are expected by increasing the number of energy groups used by RINNOVO in these core calculations.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Improvement and Verification of Online Cross Section Generation Capability of Griffin for TRISO-fueled Reactors

Griffin, a MOOSE-based reactor multiphysics code jointly developed by Idaho National Laboratory and Argonne National Laboratory under the DOE Office of Nuclear Energy’s NEAMS program, has pursued the development of an online multigroup cross section generation capability for a few years to enable high-fidelity, problem-dependent neutronics analyses of advanced thermal reactors. Recent advancements in Griffin’s online multigroup cross section generation capability have significantly improved the accuracy, robustness, and efficiency of self-shielding calculations for both prismatic and pebble-bed TRISO-fueled reactor applications. Key developments include a unified fuel self-shielding method applicable to both TRISO and annular compact/spherical shell fuel zone geometries; an advanced Dancoff Category-based Equivalence Theory using a bell function for non-fuel resonance treatment, achieving more than an order-of-magnitude speedup compared to the Tone method; an on-the-fly multigroup equivalence approach to mitigate group condensation errors; and a streaming correction method for pebble-bed homogenization. A proof-of-concept demonstration of on-the-fly group condensation with consistent P0 transport correction was also achieved. The method reproduced direct fine-group solutions with excellent accuracy (eigenvalue errors within 10 pcm and pin-power differences within 0.5%), but due to performance limitations of the current fixed-source solver, improvements to solver efficiency will be addressed in future work. Verification tests were performed on graphite-moderated TRISO-fueled two-dimensional core benchmark problems representing gas-cooled microreactors, heat pipe-cooled microreactors, gas-cooled pebble-bed reactors, and fluoride salt-cooled high-temperature reactors. Across all cases, Griffin showed excellent agreement with Serpent2 continuous energy Monte Carlo solutions: eigenvalue errors within 200 pcm, pin-power root-mean-square errors within 2%, and control rod and drum worth errors less than 2%. It should be noted that, for the benchmark problem, cross section generation contributed less than 3% of the total simulation times. These results demonstrate that Griffin’s online cross section generation capability delivers accurate and efficient reactor physics solutions across a wide spectrum of TRISO-fueled advanced reactor designs. With further improvements to the fine-group fixed-source solver and planned extensions to depletion, transients, and coupled neutron–gamma transport, Griffin will be well-positioned to become a powerful and comprehensive tool for advanced reactor analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

SCALE Modeling of the Fast Spectrum Heat Pipe Reactor

As part of the severe accident analysis collaboration with Sandia National Laboratories (SNL) and the Nuclear Regulatory Commission (NRC), SCALE models were developed for a fast-spectrum heat pipe reactor. These models were based on the Idaho National Laboratory (INL) Design A concept, which is an alternative design to the Los Alamos National Laboratory (LANL) Special Purpose Reactor (SPR), also known as the Megapower reactor. The model contains 1,134 heat pipes, surrounded by hexagonal fuel elements, with a potassium working fluid; the fuel is UO 2 with 19.75 wt% 235 U enrichment. The model contains axial beryllium oxide (BeO) reflectors above and below the active fuel region along with a radial alumina reflector containing 12 B 4 C control drums. The center of the core is left unfueled to make room for two shutdown control rods, one annular and one solid. The active region of the core was discretized into twenty axial and five radial zones to analyze spatial variations in power and burnup. Infinite lattice unit cell sensitivity studies were used to perform verification between the SCALE and INL models. The eigenvalue results agreed well with the reported results to within roughly 50 percent mille (pcm). Full-core model verification was performed by analyzing system eigenvalues with differing configurations of control drum and shutdown rod positions. These full core results all had eigenvalue differences less than 310 pcm. Control drum and shutdown rod worths were also compared, with differences of 3.2% or less. Using the verified model, the isotopic inventory and decay heat, as well as temperature feedback coefficients, were calculated and provided to SNL as input to the MELCOR severe accident code to analyze potential releases from this class of reactor. The results of the MELCOR analysis are provided in a different report.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Innovative control mechanism for research and test reactors using mandrel-shaped control rods

Research and test reactors have historically played a pivotal role in supporting the initial development of nuclear reactors. They continue to provide essential data for enhancing fuel designs and material knowledge. However, with many such reactors aging and the growing demand for data to bolster advanced reactor development, it is more necessary to research potential design attributes of the next generation of research and test reactors. For test reactors dedicated to fuel and material testing, the design of control mechanisms significantly influences the stabilization of neutron flux levels in irradiation positions while sustaining criticality. This study presents an innovative control mechanism for potential research and test reactor designs. It employs small absorber rods that move in opposite axial directions to maintain axial symmetry of power and neutron flux during burnup cycles. These rods maximize reactivity worth while also offering flexibility to flatten the radial power distribution. An axial translation of the control mechanisms’ absorbers, as compared to the rotational movement of absorbers in control cylinders, also provides a benefit to available excess reactivity and cycle length. Additionally, this work utilizes a simplified core model of the Advanced Test Reactor to assess the performance of this control mechanism. Compared to the current control system based on rotating control cylinders, the new control mechanism has the potential to enhance, or at least maintain, neutronic performance parameters in this reactor design.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Benchmark Exercise for the Control Rod Swelling Evaluation

The VTR core has six reactivity control assemblies and three safety assemblies. The control assemblies or primary control rods are adjusted during the normal operation to balance the core reactivity and to control the reactor power. A typical control assembly radial layout is presented in Figure 1. The figure shows the swelled absorber (B 4 C) rod. Initially, helium gas fills the gap between the pin and the cladding before irradiation swelling takes place. For VTR, HT9 steel was selected as the cladding and duct material. The main neutron absorbing material used in the VTR is B 4 C. When residing in the core, the neutronics, thermophysical, and mechanical properties of the materials used in a control assembly will degrade due to accumulated neutron damage. Material degradation limits how long a control assembly can reside in the core. Many phenomena affect the control assembly lifetime, such as the loss of reactivity worth due to B 4 C depletion, the mechanical interaction of the absorber rod and the cladding due to B 4 C swelling, the helium gas buildup in the pin due to B-10 capture, etc. B 4 C swelling, which causes closure of the gap between the absorber rod and the cladding, is usually considered as the main limiting factor from past experience. An initial study was conducted at PNNL to evaluate the irradiation behavior of a VTR control assembly. The evaluation was performed using the CNRD2 code that was initially developed for the FFTF. The study also included an assessment of the VTR control assembly and focused on a 61-pin control assembly design, which is different from that used (37-pin design) in the core design study. The study conducted by PNNL was reviewed independently by ANL. A Python script referred to as the Control Assembly Evaluation Script (CAES) was developed for the independent review and additional assessment of 37-pin control assembly design. The script has focused on the assessment of the absorber rod swelling for its importance in determining the control assembly lifetime. CAES uses geometry, neutronics, materials data as input to predict the swelling of the absorber rod during its residence in the reactor core. The results from CAES showed some non-negligible differences against the PNNL results. Some of the differences can be attributed to the different interpretation of the control rod assembly dimensions. To resolve this issue, a benchmark exercise was proposed. The benchmark specification was developed by PNNL. The benchmark exercise was performed independently at PNNL and ANL using different codes/scripts (CRND2 and CAES). This memo documents the results calculated using the different codes. However, this report is limited to presenting the results obtained. Further investigation of the cause of the observed difference will be performed as part of future activities, pending continuation of the VTR program.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Chicago Pile-1 paved the way for nuclear science and a lab in Los Alamos First self-sustaining nuclear chain reaction was nearly 80 years ago

On a bitter-cold winter day, 43 scientists gathered at an abandoned squash court at the University of Chicago where they would ultimately enable a secret lab in Los Alamos to change the world just years later. It was December 2, 1942. The group, led by Italian physicist and Nobel laureate Enrico Fermi, stacked graphite bricks, piling 57 layers that totaled more than 770,000 pounds. Later named Chicago Pile-1, their goal was to create the world’s first self-sustaining, controlled nuclear chain reaction. Inside the approximately 20-feet-tall pile were smaller blocks of uranium and control rods that, when removed, would cause the reaction to go critical – meaning create a nuclear chain reaction. It was roughly $1 million worth of materials, equivalent to nearly $16 million today, and a concept that a nuclear chain reaction would allow the weaponization of the atom. “Its success would be the crucial proof needed to know it would be possible to create an atomic bomb,” said LANL Historian Roger Meade (C-NR). “This was the precursor to the Lab we have today, nearly 80 years later.”

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

High-Fidelity Measurements for Flattop-HEU Benchmark Reevaluation

Flattop was first built in the 1950’s at Los Alamos National Laboratory. Flattop-HEU is composed of a sphere of highly enriched uranium (HEU) surrounded by a thick spherical natural uranium (NU) reflector. The reflector is composed of three parts: a stationary hemisphere and two movable quarter spheres. For fine control of the reactivity of the system, there are three control rods of natural uranium located in voids in the stationary hemisphere. The final components that make Flattop a useful critical assembly are the glory hole and mass adjustment pieces. These pieces can be loaded in various configurations into the glory hole and the core pedestal to control the known worth of the system. The glory hole and mass adjustment pieces are mostly small pieces of HEU with some mass adjustment pieces fabricated from NU. This allows for the irradiation of samples to a specified level. To better document the system, Flattop was evaluated and included in the International Criticality Safety Benchmark Evaluation Project (ICSBEP) handbook. The original benchmark evaluation of Flattop-HEU was written in 1999 based on an experiment completed in the 1960’s. This original evaluation was written to provide a single diameter that defined critical mass; however, as computational capabilities have increased, the focus for benchmark evaluations has shifted to include detailed modelswith all physical dimensions. Thus, as Flattop is a lynchpin in critical experiment work, the benchmark is being reevaluated at current standards. This summary discusses some of the largest known uncertainties from the evaluation and the high-fidelity measurements taken to reduce these uncertainties.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2022

Investigation of the OECD/NEA PWR MOX/UO{sub 2} core transient benchmark using a coupled whole core pin-by-pin route in WIMS

This paper investigates the new CAMELOT coupling route developed in WIMS, which uses prepared tabulated cross-sections combined with the flux solver MERLIN module and the integrated sub-channel thermal hydraulics solver ARTHUR module to solve for the coupled core state. This coupling route is applied to the OECD/NEA pressurized water reactor (PWR) core transient benchmark exercise, involving 3-D neutronics and thermal-hydraulics analysis of a full-sized PWR core featuring MOX fuel at various core states, as well as simulating a control rod ejection transient. Analysis of the results for each part of the benchmark demonstrate that the high-resolution pin-by-pin analysis used in WIMS closely matches the other participants in the benchmark, with differences of between 20-200 pcm for k-effective, 1-2.5% RMS difference for assembly-averaged powers, 2-8% difference in rod worth calculations, < 200 ppm for calculations of critical boron concentration, zero difference in delayed neutron fraction calculation and a very close match for power response to the rod ejection transient. This is the first application of the WIMS-CAMELOT approach for the transient analysis of a full-sized reactor and as such further improvements such as reducing computational cost through parallelization, optimisation and memory reduction, are currently in active development. This paper demonstrates the accuracy and flexibility of the CAMELOT coupling route, providing a straightforward and consistent user image that can be easily used for a variety of modelling problems involving coupled neutronics and thermal-hydraulics. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2024

Determining Reactor Operating Conditions to Prevent Xenon Walkaway Accident

The xenon walk-away is a postulated uncontrolled criticality that could occur if operators are forced to evacuate the control room in certain circumstances. This poster shows how previously held assumptions are not entirely valid and that less restrictions should be placed on operators and engineers than initially assumed.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗