Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “reactor cooling”

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 721 records · Page 40

Progress on the US-Japan Source Term Benchmark Analysis Collaboration

There has been a resurgence in global interest in advanced reactor technology, with a variety of innovative concepts being proposed and developed. Sodium-cooled Fast Reactors (SFRs) have garnered considerable attention given their beneficial characteristics and extensive historical development programs and operating experience. Central to the licensing of non-light water reactor (non-LWR) technology is the characterization of the safety case. At its basics, reactor safety focuses on the prevention and mitigation of the release of radioactive material to the environment.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Development and application of two-step uncertainty propagation and sensitivity analysis methodology for fast reactor safety analysis

Uncertainty quantification (UQ) in nuclear reactors for transients is directly linked with safety assessment through the cross-sections uncertainties, provided as a covariance matrix, which are propagated through the reactor system to output of interest pertaining to reactor safety, such as peak temperatures in fuel/clad/coolant. Using a two-step approach, uncertainties are first quantified and propagated from basic input variables (such as reaction cross-sections) to intermediate quantities (such as reactivity feedback coefficients) through lattice level calculations. Uncertainties of intermediate quantities (from the first step) are then propagated through the system transient calculations, in the second step, to obtain uncertainties on reactor safety output parameters of interest. The scope of this work consists of Uncertainty Quantification & Propagation of nuclear data uncertainties that are highly correlated through unprotected transient overpower and unprotected loss of flow to assess their impact on core safety parameters. This two-step approach in the presence of covariance renders the sensitivity analysis very challenging. In fact, usually the sensitivity analysis is restricted to each step, which limits its application since the sensitivities between the system output quantities and the basic input variables are difficult to obtain. Here, in this work, we address this issue by proposing a simple, general methodology to combine the sensitivity indices obtained in each step by assuming the model behavior being linear. For the first step Generalized Perturbation theory based indices are used while in the second step the recently studied Johnson indices. The uncertainty quantification and sensitivity methodologies discussed here are demonstrated on a generic LFR design which is based on the 500 MWth demonstration Lead-cooled fast reactor (DLFR) using oxide fuel, developed by Westinghouse Electric Company (WEC).

42 - ENGINEERING↗

Modeling of HTTF test PG-26 using RELAP5-3D and SAM

The High Temperature Test Facility (HTTF) at Oregon State University (OSU) is a scaled integral effects experiment designed to investigate transient behavior in high-temperature gas-cooled nuclear reactors (HTGR) with prismatic fuel and reflector blocks [1]. Several tests have been completed at the HTTF including depressurized conduction cooldown (DCC) and pressurized conduction cooldown (PCC) transients. This summary reports on the analysis of test PG-26 using the INL system code REALP5-3D [2] as well as the ANL system code SAM [3]. Test PG-26 is a progression of the Double Ended Inlet-Outlet Crossover Duct Break transient that is referred to as a DCC [4]. Core initial conditions (i.e., before the DCC started) have been met using low power (<100 kW) and two of ten available electric heater banks. The DCC transient was initiated during the 50th hour of the test. The break valves were opened, and hot helium from the core and cold helium from the reactor cavity simulation tank (RCST) started mixing. The gases flowed in a countercurrent fashion, where the top half of the hot duct contained hot helium that flowed in one direction and cold helium that flowed in the other direction in the bottom half of the duct. After the pressure and density reached equilibrium, the event entered a diffusion mode. The onset of a reverse natural circulation was not observed during the DCC period of the test.

97 MATHEMATICS AND COMPUTING↗

Tellurium Corrosion of Type 304/304L Stainless Steel, Iron, Chromium, and Nickel in High-Temperature Liquid Sodium

Investigating tellurium (Te) corrosion on structural materials is crucial for sodium-cooled fast reactors (SFRs) due to radionuclide presence and knowledge gaps. In this study, Type 304/304L stainless steel (SS304), chromium (Cr), iron (Fe), and nickel (Ni) samples were immersed in low-oxygen environments with Te in liquid sodium at 773 K for 30 days. At 10 ppm oxygen, SS304 showed multiple oxide layers, including a compact NaCrO 2 interlayer and porous Na-Fe-Ni-O outer layers. Tellurium penetrated through the porous layers but was hindered by the NaCrO 2 interlayer. At 0.01 ppm oxygen, Cr had no oxide layer, while Fe and Ni had unstable ones. Tellurium-induced pitting was deeper in Fe and Ni compared to Cr. Oxygen levels and Cr composition are critical factors affecting stable oxide compound layer formation and mitigating Te-induced pitting.

36 MATERIALS SCIENCE↗

An investigation of the phase behaviors for quaternary U-Nb/Mo-Ti-Zr metallic fuel alloys

We report quaternary fuel alloys containing U, Nb/Mo, Ti, and Zr are proposed as fuel candidates for sodium-cooled fast reactors (SFRs). In this work, two Nb-bearing alloys, i.e., U-NT5Z (U-2.5Nb-2.5Ti-5.0Zr in wt%) and U-NT7Z (U-1.5Nb-1.5Ti-7.0Zr in wt%), and two Mo-bearing alloys, i.e., U-MT5Z (U-2.5Mo-2.5Ti-5.0Zr in wt%) and U-MT7Z (U-1.5Mo-1.5Ti-7.0Zr in wt%) were characterized and compared. The characterization techniques were differential scanning calorimetry (DSC), X-ray powder diffraction (XRD), and scanning electron microscopy/energy dispersive X-ray spectroscopy (SEM/EDS). DSC was performed to obtain the transition behaviors, and XRD and SEM/EDS were applied for phase identification. The results were combined to obtain the solid-state phase transitions between 500 °C and 850 °C for the alloys. It is found that the Nb-bearing alloys comprise similar phase transition behaviors as the Mo-bearing alloys. The phase transitions in U-NT5Z are a →γ at 608 °C and U 2 Ti → γ at 627 °C, which are ~40 °C higher than that of U-MT5Z. The phase transition in U-NT7Z is a → γ at 645 °C, and is 23 °C higher than that of U-MT7Z.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Reactor Containment Passive Safety Analysis: Steam Condensation in Presence of Non-condensable Gas Scaled Experiment and Modeling

This study presents steam condensation scaled experiments and semi-empirical models in presence of nitrogen (N)—a noncondensable gas (NCG), simulating air in the reactor containment—to support water-cooled small modular reactors (SMRs) passive containment cooling system (PCCS) design and analysis. Previous experimental studies on PCCS are focused on fixed and smaller tube (mostly 2-in.) geometries and specific test condition variations, bringing challenges with geometric scaling and mismatching with SMR prototypic design. To address these challenges, this study presents steam condensation test dataset obtained from three scaled test sections of 1-, 2-, and 4-in.-diameter steam condensers with an annular/jacket cooling of 2-, 3-, and 6 in.-diameter tubes, respectively. Test data were collected for steam ranges from 58 to 63 kg/hr., and NCG flow of 4.4 to 13.3 kg/hr. Annular cooling water flow was varied to obtain required testing conditions of saturated steam inlet and fully condensed outlet. Axial temperature test data of bulk cooling water, steam and condensate were collected by thermocouples for three test sections and various steam-NCG mixing/testing conditions. A standard data reduction method was adopted—utilizing iterative and nodalized mass and heat transfer calculation—to estimate axial local heat fluxes, heat transfer coefficients (HTCs), condensation rates, film thickness, and Nusselt number. Based on the obtained dataset semi-empirical model results—a ratio of experimental and Nusselt’s theoretical HTC are presented. Such results and findings are supportive of developing scaled-up testing facility, to enable model validations and accelerate next generation of reactors development and deployment

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Presentation: Reactor Containment Passive Safety Analysis: Steam Condensation in Presence of Non-condensable Gas Scaled Experiment and Modeling

This study presents steam condensation scaled experiments and semi-empirical models in presence of nitrogen--a noncondensable gas (NCG), simulating air in the reactor containment--to support water-cooled small modular reactors (SMRs) passive containment cooling system (PCCS) design and analysis. Previous experimental studies on PCCS are focused on fixed and smaller tube (mostly 2-in.) geometries and specific test condition variations, bringing challenges with geometric scaling and mismatching with SMR prototypic design. To address these challenges, this study presents steam condensation test dataset obtained from three scaled test sections of 1-, 2-, and 4-in.-diameter steam condensers with an annular/jacket cooling of 2-, 3-, and 6 in.-diameter tubes, respectively. Test data were collected for steam ranges from 58 to 63 kg/hr., and NCG flow of 4.4 to 13.3 kg/hr. Annular cooling water flow was varied to obtain required testing conditions of saturated steam inlet and fully condensed outlet. Axial temperature test data of bulk cooling water, steam and condensate were collected by thermocouples for three test sections and various steam-NCG mixing/testing conditions. A standard data reduction method was adopted--utilizing iterative and nodalized mass and heat transfer calculation to estimate axial local heat fluxes, heat transfer coefficients (HTCs), condensation rates, film thickness, and Nusselt number. Based on the obtained dataset semi-empirical model results--a ratio of experimental and Nusselt's theoretical HTC are presented. Such results and findings are supportive of developing scaled-up testing facility, to enable model validations and accelerate next generation of reactors development and deployment.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Propagating neutronic uncertainties for FFTF LOFWOS Test #13

The safety evaluation conducted for licensing a Sodium-cooled Fast Reactor (SFR) may require modeling transients with best-estimate calculation tools that must first be validated against real-world measurements. To provide the community with a valuable benchmarking opportunity for validating SFR analysis tools and methods, the International Atomic Energy Agency (IAEA) initiated a coordinated research project (CRP) in 2018 for the analysis of the Fast Flux Test Facility (FFTF) Loss of Flow Without Scram (LOFWOS) Test #13. The impact of nuclear data uncertainties on neutronics parameters was previously investigated based on the COMMARA-2.0 covariance matrix. Since the transient simulation results are very sensitive to certain reactivity coefficients, it was decided to employ rigorous uncertainty propagation methods to quantify the impact of nuclear data uncertainties on the best-estimate predication of FFTF LOFWOS Test #13. The DAKOTA code is used to propagate neutronic uncertainties through SAS4A/SASSYS-1 transient simulations, while taking into account spatial and reaction-wise correlations within these uncertainties. This study shows that the remaining discrepancies observed between the Argonne National Laboratory (ANL) best-estimate results and the experimental measurements can be partly explained by the uncertainty associated with Gas Expansion Module (GEM) worth, which contributes the majority of the overall nuclear data uncertainty on the output from the FFTF LOFWOS Test #13 transient simulation. This study also confirmed the importance of including spatial and reaction-wise correlations of nuclear data uncertainties on feedback coefficients in the uncertainty propagation to avoid under-estimating their impact during the transient simulations.

42 ENGINEERING↗

Comparative post-irradiation examination of high burnup U-19Pu-10Zr: Assessing steady-state irradiation behavior against historical and modeled fuel performance

Here, the development of next-generation sodium-cooled fast reactors necessitates comprehensive research on metallic fuels to maximize economic performance while ensuring safe operation. In this study, we investigated the steady-state irradiation behavior of two high burnup U-19Pu-10Zr fuel pins, DP-36 and DP-40, in preparation for planned safety testing. Post-irradiation examination (PIE) was performed to quantify fuel column elongation, regions of low-density at the top of the fuel column, pin deformation, fission product distribution, fractional fission gas release, microstructural evolution, and fuel constituent redistribution. Benchmarking against existing PIE data from U-19Pu-10Zr fuel pins irradiated in EBR-II revealed consistent patterns in fuel column elongation and cladding diametral strain. However, both pins exhibited longer low-density structures, and destructive examination of DP-36 revealed more complex constituent redistribution patterns compared to previously reported data for ternary fuel pins. The steady-state irradiation of both pins was also modeled using BISON. Comparisons of PIE results with modeled predictions showed overall agreement in fractional fission gas release but consistent overestimation of axial and radial swelling due to gaseous and solid swelling models. These findings underscore the critical importance of pre-test characterization on test and sibling pins to accurately capture steady-state fuel behavior ahead of transient testing, thus establishing a baseline for post-test comparison. Additionally, these analyses identified key data gaps that warrant further investigation to improve the understanding and prediction of fuel swelling, thereby enhancing the synergy between modeling and experimental efforts in supporting accident testing.

BISON↗

Carburization and decarburization behavior of Grade 91 ferritic-martensitic steel in liquid sodium environments

Here, this paper presents a study of carbon transfer and its effect on microstructure and tensile properties of Grade 91 (G91) ferritic-martensitic steel exposed to sodium at 550–650 °C. Sodium exposure tests were conducted in Argonne's forced convection sodium loops up to exposure times of ~40,000 h. Thermal aging study of G91 steel was conducted in parallel to isolate the thermal aging effect from the sodium effect. It was found that sodium exposures at 650 °C dissolved M 23 C 6 carbides, eliminated the martensite subgrain structure resulting in excessive grain growth and reduced the tensile strength by >50%, while sodium exposures at 550 and 600 °C had an insignificant effect on its microstructure and tensile properties. These effects were attributed to the carburization/decarburization process of G91 steel in sodium environments. Carbon concentrations in sodium were determined by a foil equilibration method. The estimated carbon concentration was in the range of 0.8–1.2 ppm in the SMT-1 loop and 0.3–0.7 ppm in the SMT-2 loop. Thermodynamic analysis of the carburization – decarburization process was conducted for G91 steels exposed in sodium environments. The carbon activity-concentration relationship for G91 was evaluated by considering four phases in G91, i.e. bcc ferrite, M 23 C 6 , NbC and VC carbides. It was found that the carburization-decarburization process in G91 steel was dictated by M 23 C 6 carbides at high carbon activities, while NbC and VC carbides dominated the process at low carbon activities. The calculated carburization-decarburization boundary showed that G91 would undergo decarburization at 650 °C and carburization at 550 °C in the sodium loop environments, which was consistent with our experimental observations. This experimental and theoretical analysis provided a basis for predicting the effect of carbon transfer on the integrity of reactor components in sodium environments and for the design of new alloys used in sodium-cooled fast reactors.

36 MATERIALS SCIENCE↗

Phase 1 Specifications and preliminary sensitivity analyses of the OECD-NRC liquid metal fast reactor core thermal-hydraulics benchmark

The developments in wire-wrapped fuel bundle thermal-hydraulics research have become increasingly diverse and this makes comparisons across studies difficult. This paper describes an ongoing benchmark that will address this issue by directly comparing computational results with experimental data from two unique wire-wrapped thermal-hydraulic facilities. Preliminary subchannel sensitivity / uncertainty analyses show largely expected behavior for the presented cases, where rod diameter and friction factor show the greatest sensitivities to geometrical perturbations. Guidelines for computational fluid dynamics evaluations have been prepared and highlight the need for consolidated modeling strategies. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Coupled Monte Carlo and thermal-hydraulics modeling of a prismatic gas reactor fuel assembly using Cardinal

Cardinal is a MOOSE application that couples OpenMC Monte Carlo transport and NekRS computational fluid dynamics to the MOOSE framework, closing the neutronics and thermal-fluid gaps in conducting tightly-coupled, high-resolution multiscale and multiphysics analyses. By leveraging MOOSE's interfaces for wrapping external codes, Cardinal overcomes many challenges encountered in earlier multiphysics coupling works, such as file-based I/O or overly-restrictive geometry mapping requirements. In this work, we leverage a subset of the multiphysics interfaces in Cardinal to perform coupling of OpenMC neutron transport, MOOSE heat conduction, and THM thermal-fluids for steady-state modeling of a prismatic gas reactor fuel assembly.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Assessment of the Griffin Reactor Multiphysics Application Using the Empire Micro Reactor Design Concept

In late 2019, INL and ANL agreed to jointly develop the reactor physics code named Griffin based on the integration of the two code suites, MAMMOTH/Rattlesnake (INL) and MC2 - 3/PROTEUS (ANL). Griffin is being developed based on the MOOSE framework and MOOSE quality assurance procedures. This decision was made to be able to allow DOE-NE to efficiently invest funding to this area and to provide effective and timely support for existing and potential users; the latter includes industry and government organizations who are developing various types of advanced reactors in the near and long term. Since MAMMOTH/Rattlesnake has been developed based on the MOOSE framework, the INL/ANL Griffin development team agreed to build Griffin beginning with a merger of MAMMOTH and Rattlesnake into a single code and moving forward by implementing capabilities from the PROTEUS suite into Griffin. Moving forward, both ANL and INL efforts are equally invested in the Griffin project, with management support, to provide an advanced reactor multiphysics tool to assist in reactor design, optimization, and safety analysis. Much work remains in moving Griffin forward to migrate PROTEUS capabilities and to optimize performance to meet user needs. The main objective of this work is to assess the current status of Griffin capabilities in terms of performance and accuracy, to determine priorities for PROTEUS migration, and to identify capabilities and features to improve for supporting the code integration effort. For this assessment, the Empire micro reactor problem that was developed in the ARPA-E MEITNER program was selected as an advance reactor concept of interest to the technical community. The Empire reactor problem was expanded from its original incomplete specification to be a small heat-pipe-cooled micro reactor core with ~113 cm radius and 70 cm in height, composed of 18 fuel assemblies, 12 control drums, and beryllium radial and axial reflectors. In the current model, using 5 cm axial reflectors specified in the original Empire assembly model, more than 10% of neutrons leak axially and through the empty center safety hole, as well as through heat pipe channels in fuel assembly elements that extend through the top reflector region. Several calculation models of the core were defined for systematic assessment, including 2-D and 3-D fuel assemblies and whole cores with cylindrical boundaries. Cross sections were generated using Serpent 2, and meshes were produced using the Argonne mesh tool or the INL neutronics meshing tools combined with CUBIT. Cross sections and meshes were converted to the ISOXML and Exodus formats, respectively, so that Griffin and PROTEUS could use consistent data for solving the reactor problems. With the prepared cross sections and meshes, PROTEUS was run first to ensure that all input data were correctly generated and input options in terms of angle, mesh, and energy group were accurately determined. Comparisons against Serpent 2 solutions were made in terms of eigenvalue and pin power. The same calculations and comparisons were then conducted using Griffin. For the fuel assembly and whole core problems, the PROTEUS eigenvalues agreed well with reference Serpent 2 solutions within 100 and 30 pcm, respectively, and pin power differences relative to Serpent 2 were overall less than 2.2% and RMS 0.8% for the whole core models. This indicated that all input data were properly prepared. Using the same data, Griffin was run selecting the SAAF-CFEM SN solver with Legendre-Gaussian quadrature and NDA and DSA for acceleration. It was found that the SAAF-CFEM solver of Griffin required finer meshes to achieve eigenvalue and pin power solutions in good agreement with Serpent 2, consequently requiring more memory requirement and longer computation time. On the other hand, the SPH-Diffusion 2-D core calculations performed using Griffin were able to recover the exact eigenvalue from the reference Serpent 2 solutions, resulting in a pin-power distribution with an RMS of 0.6% and maximum absolute difference of less than 1.4%. The runtimes for SPH-Diffusion for the 2-D core were less than 3 minutes on 40 cores. During this evolution of this evaluation, many updates were made in Griffin by the Griffin development team of INL (focusing on software updates) and ANL (reviewing and supporting software updates) to complete this assessment. Observations from the code assessment are presented in the conclusion section of this report, followed by a discussion of recommendations for future work.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Development and benchmarking of transient nodal code SIMULATE5-K neutron kinetics solver

SIMULATE5-K is Studsvik's next generation best estimate transient code. The time dependent diffusion equation is solved with a nodal method consistent with that implemented in the licensed core design code SIMULATE5. Arbitrary number of neutron and delayed neutron precursor groups can be used. For the solution of the spatial problem, the coupling coefficients used to relate the node leakages are found by first converting the time dependent diffusion equation to a static diffusion equation with the use of flux and delayed neutron precursor dynamic frequencies. Once the static-like equations are obtained, the multi-group analytical nodal model is used to obtain the coupling coefficients, expressing the node leakage in terms of adjacent node average fluxes. The coupling coefficients are then inserted into the time dependent nodal balance equation. For the time integration, the time dependent neutron balance equation is solved with the frequency transformation method. The treatment of the temporal dependence yields a fixed source problem which can be solved utilizing the existing fixed-source methodology. The primary purpose of this paper is to describe the neutron kinetics methodology implemented in SIMULATE5-K. The accuracy of the method is demonstrated for a series of well-known, neutronic-only benchmark problems. (author)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Performance evaluation of a Terry GS-2 steam impulse turbine with air-water mixtures

Terry steam turbines are widely used in various industries because of their robust design. Within the nuclear power generation industry, they are used in the Reactor Core Isolation Cooling System to remove decay heat during reactor isolation events. During the Fukushima Daiichi nuclear power station disaster in Japan in 2011, the Reactor Core Isolation Cooling System and associated Terry turbine operated for over 70 hours in Unit 2; this runtime is well beyond the expected operating duration. Theories suggest the turbine was subjected to a two-phase inlet flow, which could degrade the turbine performance. In this work, an experimental test rig was constructed to test a full-scale Terry model GS-2 steam turbine under two-phase air/water flows. Steady-state efficiency and torque performance maps of the turbine were developed over a range of turbine inlet pressures (1.38–4.83 bar or 20–70 psia), air mass fractions (0.05–1.0) and rotational speeds up to 4000 RPM. Furthermore, turbine performance followed expected trends with torque varying linearly and efficiency varying quadratically with rotational speed. In addition, high-speed images of the two-phase flow entering the turbine were also analyzed to understand how changes in inlet pressure and air mass fraction affect the flow regime and homogenization. The present tests with air–water two-phase mixtures are an important step towards providing an understanding of the full-scale Terry turbine’s behavior and performance curves under two-phase conditions. The results of this work will be combined with air/water and steam/water data gathered using a small-scale Terry ZS-1 steam turbine in order to understand the scaling relationship between large and small size Terry turbines and fluid pairs. The combined data set will enable further development of analytical models over a wide range of conditions and may be used to provide technical justification for expanded use of the Terry turbines in nuclear power plant safety systems and other systems.

42 ENGINEERING↗

Segmentation and Classification of Fission as Pores in Reactor Irradiated Annular U–10Zr Metallic Fuel Using Machine Learning Models

Metallic fuels, particularly U—10Zr, are promising candidates for next-generation sodium-cooled fast reactors. Irradiation of nuclear fuels in reactors can lead to the formation of solid and gas fission product which subsequently forms microstructural pores, deteriorating fuel performance. Due to the massive amount of pores and complex phases formed, a quantitative description of fission gas pores is not yet available, preventing the development of microstructure-informed fuel performance modeling for fuel qualification. This paper applied a pre-trained deep learning model to ~10,260 high magnification scanning electron microscopy images. This method increased the accuracy of fission gas pore segmentation and allows statistical features to be extracted which cannot be achieved manually. A pre-trained decision tree model worked on the segemenation results and further classified the pores into different categories to produce a correlation between the pores, movement of lanthanides, and temperature gradient during irradiation. Finally, this paper emphasizes the potentials of machine learning models to accelerate fuel research, development, and qualification for advanced reactors.

36 MATERIALS SCIENCE↗