Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Fuel assembly bow”

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.

Modeling fuel assembly distortion in light water reactors using BISON

Individual fuel rods in a light water reactor fuel assembly (FA) often have different power profiles that result in a spatial temperature gradient in the FA. This temperature gradient causes some rods to elongate more than others which could in turn distort the entire FA. Such FA distortion has been known to cause operational difficulties with refueling and with the insertion of control rods [1]. Numerical modeling of FA distortion due to differential elongation of fuel rods requires: (i) modeling the thermal elongation as well as the creep and irradiation induced elongation of the fuel rods, and (ii) modeling the bending response of the FA due to the differential elongation of the fuel rods. Past studies on this topic have been focused on the mechanical deformation of the FA with simplistic modeling of fuel rod elongation. In this paper, a simulation approach is presented that directly couples 2D axisymmetric fuel performance models of the individual fuel rods in a FA with structural beam models of the fuel rods in a FA using the Bison code [2]. This allows all of the mechanisms internal to the fuel rod leading to cladding elongation to be accounted for in the FA distortion model. To demonstrate this capability, the bowing of a simplified and idealized 6x6 FA (Fig. 1a), modeled using beam elements, is analyzed under different spatially varying power profiles. For this study, the guide tubes are rigidly connected to the spacer grid. Also, the fuel rods are rigidly connected to the spacer grid in all directions except the vertical direction (Y) along which slipping can occur between the fuel rods and spacer grids, depending on the friction coefficient between them. Coupling the thermo-mechanical response of individual fuel rods to the global mechanical behavior of the FA can shed light on the mechanisms that cause distortion of the FA and potential solutions to address this issue.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Simulation of the Fast Reactor Fuel Assembly Duct-Bowing Reactivity Effect Using Monte Carlo Neutron Transport and Finite Element Analysis

This paper discusses a new method of simulating the fuel assembly duct-bowing reactivity coefficient for EBR-II run 138B. Quantification of the fuel assembly duct-bowing reactivity effect in liquid metal–cooled fast reactors has been a persistent problem since they were first designed and operated. Simulation of the duct-bowing reactivity effect is difficult because the level of detail required to simulate the effect has exceeded most modeling capabilities. The new method outlined in this paper utilizes the finite element analysis code ANSYS to analyze the thermal and structural components. Here, the displacement of the fuel assembly duct due to thermal expansion and mechanical interaction was calculated by ANSYS using recorded EBR-II run 138B temperature and power boundary value data. The displacement values were incorporated into to a Monte Carlo model of EBR-II run 138B and keff was calculated. Multiple Monte Carlo calculations were performed with duct displacement values corresponding to different reactor temperatures. Using the calculated keff values associated with the different duct displacement results allowed calculation of the duct-bowing reactivity coefficient. The duct-bowing reactivity coefficient was calculated to be –14.5 × 10 –4 $/°C/ ± 4.4%.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Assessment of MOOSE-Based Tools for Calculating Radial Core Expansion

Radial core expansion in liquid-metal cooled fast reactor systems is a well-known phenomenon that produces strong reactivity feedback effects. An inherently safe reactor design takes advantage of negative reactivity feedback in accident scenarios by utilizing a core restraint system which produces a bowed shape that allows for radial expansion of the fuel regions. Detailed modeling and simulation of core radial expansion itself as well as subsequent reactivity feedback is a challenging task involving contact of many fuel assembly elements and physics feedback from neutronics, thermal hydraulics, and thermal mechanical response. A variety of physics codes have been developed to model aspects of radial core expansion but in general invoke geometrical or physics approximations. No code system currently exists which tightly and robustly couples these physics with enough detail to fully resolve the complex core radial expansion reactivity feedback effects. The future availability of such a code system is of vital importance to fully understanding the reactivity feedback effects that occur due to radial expansion, and consequently to optimizing the design of the core restraint system. A high-fidelity code will also be used to benchmark existing lower fidelity, faster running models to understand their benefits, limitations, and range of applications. A code development path using MOOSE-based tools is proposed in order to leverage the detailed geometry capabilities and natural tight coupling and robustness of MOOSE-based applications for modeling this complex phenomena. While simulation of the full phenomenon involves several physics, an assessment has been initiated on the capabilities and readiness of the currently available Tensor Mechanics module within MOOSE for calculation of the structural mechanical responses which occur within the reactor core. This report focuses on modeling the force-deformation response which mimics the physics of a duct contact deformation, as well as differential thermal expansion which produces a thermal bowed shaped for the fuel assemblies in a core. Simple examples were initially performed such as simple supported beam bending under load. The complexity of examples was progressively increased to better mimic the duct behavior by including a differential thermal example and inclusion of hexagonal cross-sections in the geometry. Further assessment of the structural mechanical response simulation capability is still required for modeling duct contact interactions and irradiation creep and swelling. Companion thermal hydraulic and neutronics assessments will also be required; these activities are planned for future years. Finally, integration of the multiple physics components through MOOSE is required to predict the core radial expansion and subsequent feedback effects.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Long-Length SiC/SiC Tube Irradiation Experiment at LWR Temperatures Inserted into HFIR

This report describes the pre-irradiation characterization and successful assembly of a High Flux Isotope Reactor (HFIR) irradiation experiment designed to assess radiation-induced lateral bowing of silicon carbide fiber–reinforced, silicon carbide matrix composite (SiC/SiC) components under a radial fast neutron flux gradient at representative light-water reactor (LWR) temperatures of approximately 300 °C. The over-arching goal of this work is to validate thermomechanical models of SiC/SiC composite fuel assembly distortion used to predict the extent of lateral bowing and dose at which bowing reaches a maximum value. Ultimately, these models are used to ensure that fuel assembly distortion does not block coolant channels, interfere with control rod/blade movements, and/or impact local reactivity. The experiment contains six cladding tube specimens that are approximately 558.5 mm long with representative pressurized water reactor (PWR) outer and inner diameters of 9.5 mm and 7.5 mm, respectively. The specimens were thoroughly characterized prior to irradiation using traditional dimensional inspection and surface profilometry so that these measurements can later be compared with similar measurements that will be made post-irradiation to determine radiation-induced deformations. Furthermore, fine engraving markers were inscribed along the outer surfaces of the specimen and mapped using analytic photography and a three-dimensional stage. This technique allowed for accurate measurements of the marker spacings, which can be compared with similar measurements that will be made post-irradiation to provide local radiation-induced strain mapping. A light curtain optical micrometer was also employed to measure the pre-irradiation bowing of each specimen. The experiment was successfully assembled and inserted into position VXF-19 for HFIR cycle 516, which began May 12, 2026. The experiment will be removed after cycle 519, which is currently scheduled to end in December 2026.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Final Report on Predictive Analyses of PRD as Function of Anomalies

This report summarizes the work completed in FY-2024 to analyze the power reactivity decrement (PRD) concepts of the ARC-100 core. The PRD has been traditionally defined by the reactivity change from a hot zero power (HZP) to a particular power. Consequently, the PRD accounts for the core reactivity changes due to increase coolant temperature gradient axially and radially across the core, and increased fuel temperature. The coolant temperature gradient leads to sodium and structure density changes, to radial core expansion from assembly flowering and bowing (due to axial and radial temperature gradients within the assemblies), and to control rod driveline thermal expansion. The fuel temperature increase associated with coolant temperature and power increases leads to Doppler effect and axial thermal expansion. In this work, the normal operating Hot Full Power (HFP) state is the only power level of interest, so analyses focus on the PRD calculated from HZP to HFP. The PRD has been used to assess the reactor safety features asymptotically in unprotected accident scenarios, including the loss of heat sink (LOHS), loss of flow (LOF), and transient overpower (TOP) without scram. The PRD concept relies on the “global” reactivity coefficients A, B, and C that are estimated based on “individual” reactivity effects (Doppler, sodium density, etc.). The objectives of this work are: 1) to improve and verify the methodology used to calculate the ABC coefficients used in the PRD, 2) to assess if the PRD can be used to reliably identify abnormal events. This report fulfills the FY-2024 scope of WBS#1.15.8.3 activity, “ANL0120 – Predictive Analyses of PRD as function of deformation”. The PRD concept is described in Sections 2. Additional effort in refining the methodology for core bowing modeling is performed in Sections 3. Then, two verification exercises are proposed in Section 4 to benchmark these coefficients based on direct neutronic-only calculations and on dynamic core transient simulations. Finally, the PRD approach is assessed for the detection of several unexpected events, such as primary flow perturbation or improper fuel loading, in Section 5.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Assembly Bowing Reactivity Calculation Methodology Applied to Lead Fast Reactor

Ducted assemblies bow during operation due to power and temperature gradients which can be influenced by operating flow rates. For fast spectrum cores using ducted assemblies, the bowing behavior follows that of the duct and because there are gaps between the ducts, the bowing can result in compaction or expansion of the active core. This local displacement can have a positive or negative impact and knowing the net effect during transients is important for system reactivity control. Keeping the net bowing reactivity worth low is possible with attentive placement of load pads above the active core and selecting load pad gap thicknesses that result in a desired behavior at standard operating conditions. This paper considers a Lead Fast Reactor (LFR) concept fueled by HALEU UO2 developed by Westinghouse Electric Company (WEC) and applies a workflow of Argonne-developed codes to estimate the core bowing reactivity worth. Using orifice flow rates grouped by assembly type, the net reactivity impact due to core assembly bowing for the LFR was found to be small and in line with other liquid metal fast reactors: +29/+32/+35 pcm, or about +$0.049/+$0.055/+$0.059, for BOEC/MOEC/EOEC, respectively.

bowing reactivity↗

Impact of Control Blade Insertion on the Deformation Behavior of SiC-SiC Channel Boxes in Boiling Water Reactors

This report describes the analysis of distortion of a Silicon carbide fiber-reinforced silicon carbide matrix (SiC-SiC) composite channel box under in-reactor conditions of a boiling water reactor (BWR). The BWR core has significant gradients in the fast neutron flux across the channel box due to the presence of water rods within the fuel assemblies, and these gradients increase further with the insertion of control blades. As a result of temperature and neutron flux dependent irradiation induced swelling of SiC, the SiC-SiC composite channel box can undergo significant distortion. In this work, we evaluate the SiC-SiC channel box distortion for three different control blade positions. This analysis is based on the neutron flux and temperature distributions in the BWR core calculated using the neutronics code MPACT and thermalhydraulics code CTF. This calculation is coupled through temperature feedback. Subsequently, we have performed structural analysis based on the calculated neutron flux and temperature distributions, to determine the deformation and stress development in the channel box. The structural analysis was performed using the fuel performance modeling code BISON and the commercial finite element analysis software Abaqus. The results indicate that the channel box will undergo time dependent lateral bowing for all the control blade positions in the assembly. The bowing behavior is dominated by the swelling of SiC-SiC material under non-uniform neutron flux, and changes with variation in the control blade position. The lateral bending will cause a temporary interference between the channel box and the control blade. The developed stresses exceed the proportional limit stress of the material, which may cause matrix microcracking in the channel box. However, the stresses remain below the ultimate tensile strength of the material, and therefore, development of a full, through thickness crack in the channel box is not expected.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

HFIR SiC Bowing Test Ready to Insert

This report describes the successful assembly of a High Flux Isotope Reactor (HFIR) irradiation experiment designed to assess radiation-induced lateral bowing of silicon carbide fiber–reinforced, silicon carbide matrix composite (SiC/SiC) components under a radial fast neutron flux gradient. Excessive bowing of a SiC/SiC channel box in a boiling water reactor could potentially interfere with control blade movements. Similar concerns exist for SiC/SiC fuel cladding in light water reactors. The experiment described herein will provide experimental validation of the structural response of a miniature SiC/SiC channel box and tube specimens with pressurized water reactor diameters during irradiation. The significant radial fast neutron flux gradients that exist in the permanent reflector of HFIR were characterized using detailed three-dimensional neutronic calculations. The three-dimensional displacement damage dose rate profile and the resulting volumetric swelling in SiC were used as inputs to structural analyses that predicted the deformation and stresses in the channel box specimen. The specimens were thoroughly characterized prior to irradiation using traditional dimensional inspection and surface profilometry so that these measurements can later be compared with similar measurements that will be made post-irradiation to determine radiation-induced deformations. Furthermore, fine engraving markers were inscribed along all outer surfaces of the specimen and mapped using a digital microscope and a three-dimensional stage. This technique allowed for accurate measurements of the marker spacings, which can be compared with similar measurements that will be made post-irradiation to provide local radiation-induced strain mapping. The experiment was successfully assembled and is scheduled for insertion during HFIR cycle 492, which is currently scheduled to run from May 25, 2021 to June 18, 2021.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Impact of control blade insertion on the deformation behavior of SiC-SiC channel boxes in BWRs

This article describes the analysis of distortion of a silicon carbide fiber-reinforced, silicon carbide matrix (SiC-SiC) composite channel box under in-reactor conditions of a boiling water reactor (BWR). The BWR core has significant gradients in the fast neutron flux across the channel box due to the presence of water rods within the fuel assemblies, and these gradients increase further with the insertion of control blades. As a result of the temperature and neutron flux dependent irradiation-induced swelling of SiC, the SiC-SiC composite channel box can undergo distortion. In this work, we evaluate the SiC-SiC channel box distortion for three different control blade positions. This analysis is based on the neutron flux and temperature distributions in the BWR core calculated using the neutronics code MPACT and thermal-hydraulics code CTF. This calculation is coupled through temperature feedback. Subsequently, we have performed structural analysis based on the calculated neutron flux and temperature distributions to determine the deformation and stress development in the channel box. The structural analysis was performed using the fuel performance modeling code BISON and the commercial finite element analysis software Abaqus. The results indicate that large gradients in fast neutron flux (up to 35–40% across a single axial level) will develop across the channel box. Due to these gradients, the channel box will undergo time-dependent bending for all the control blade positions in the assembly. The time-dependent bowing behavior is dominated by the transient swelling of SiC-SiC material under non-uniform neutron flux, and changes with variation in the control blade position. The bending will cause temporary interference between the channel box and control blade, and the interference is expected to be most severe for the fully inserted control blade position. The developed stresses due to differential swelling in the channel box exceed the proportional limit stress of the material, which may cause matrix microcracking in the channel box. However, the stresses remain below the tensile strength of the material, and therefore, development of a full, through-thickness crack in the channel box is not expected. Finally, further work is recommended to explore and evaluate the mitigation strategies.

42 ENGINEERING↗

FY26 Progress on Demonstration of a Multiphysics Steady State Capability for Modeling Core Radial Expansion in SFRs

Under the U.S. Department of Energy Office of Nuclear Energy Advanced Modeling and Simulation (NEAMS) Program, an integrated multiphysics approach is being developed to model the core bowing phenomena important to liquid metal-cooled fast reactors. Core bowing is an important passive safety mechanism in liquid metal-cooled fast reactors and involves multiphysics effects including radiation transport, fluid flow, heat transfer, and mechanical response to temperature and flux gradients. This report summarizes recent progress on developing a multiphysics, MOOSE-based workflow to predict core bowing and associated reactivity feedback. Significant new capabilities in the reactor physics code Griffin - sodium backfill and pin power reconstruction for deformed geometries - were applied in this effort. This year’s work included verification, code comparisons, sensitivity studies, and coupled demonstrations that advance the state of MOOSE-based core bowing workflow. Griffin’s sodium backfill capability was verified by demonstrating that its automated treatment of geometry expansion and material-density updates reproduces manual calculations exactly, confirming solid mass conservation and proper coolant backfilling in expanded geometries. Reconstructed pin powers were compared for Griffin’s ductheterogeneous and ring-heterogeneous treatments in single-, seven-, and nineteen-assembly cases, with best agreement observed in lower-leakage configurations and the duct-heterogeneous approach offering substantially lower computational cost. Thermal-hydraulic sensitivity sensitivities showed that MOOSE SCM, SAM, and CFD are expected to produce similar deformation predictions despite variances in their temperature predictions, and that explicit treatment of inter-assembly flow becomes increasingly important as gap flow rate increases. Finally, coupled demonstrations on small multi-assembly configurations using Griffin, MOOSE Solid Mechanics, MOOSE SCM, and Heat Conduction produced physically consistent reactivity feedback from thermal expansion and bowing. The coupled demonstrations simulated grid plate expansion as well as resultant core bowing at full power conditions. Simplifications were made in current workflow, namely the assumption of instantaneous full power conditions following hot zero power, and pre-expanding the Griffin geometry axially due to lack of an axial fuel pin expansion model and temperature feedback to Griffin.

Wozniak, Nicholas↗