Search NASA⌕ Search

SEARCH · Search NASA

Results for “SubChannel”

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

Interface Tracking Investigation of Geometric Effects on the Bubbly Flow in PWR Subchannels

We report that absorbing heat from the fuel rod surface, water as coolant can undergo subcooled boiling within a pressurized water reactor (PWR) fuel rod bundle. Because of the buoyancy effect, the vapor bubbles generated will then rise along and interact with the subchannel geometries. Reliable prediction of bubble behavior is of immense importance to ensure safe and stable reactor operation. However, given a complex engineering system like a nuclear reactor, it is very challenging (if not impossible) to conduct high-resolution measurements to study bubbly flows under reactor operation conditions. The lack of a fundamental two-phase-flow database is hindering the development of accurate two-phase-flow models required in more advanced reactor designs. In response to this challenge, first-principles–based numerical simulations are emerging as an attractive alternative to produce a complementary data source along with experiments. Leveraged by the unprecedented computing power offered by state-of-the-art supercomputers, direct numerical simulation (DNS), coupled with interface tracking methods, is becoming a practical tool to investigate some of the most challenging engineering flow problems. In the presented research, turbulent bubbly flow is simulated via DNS in single PWR subchannel geometries with auxiliary structures (e.g., supporting spacer grid and mixing vanes). The geometric effects these structures exert on the bubbly flow are studied with both a conventional time-averaging approach and a novel dynamic bubble tracking method. Lastly, the new insights obtained will help inform better two-phase models that can contribute to safer and more efficient nuclear reactor systems.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Assessment of the CTF subchannel code for modeling a large-break loss-of-coolant accident reflood transient

With increased industry interest in extending reactor operating cycles, the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has been investigating the behavior of high-burnup fuel during design basis accidents such as the large-break loss-of-coolant accident (LBLOCA) with consideration for risk of fuel fragmentation, relocation, and dispersal (FFRD). As part of that activity, the NEAMS subchannel thermal/ hydraulics (T/H) code, CTF, is being used for modeling of LBLOCA and to determine the impact of subchannel resolution on results. Although CTF includes a wide range of models for LBLOCA conditions, the code has not been used for this application while maintained at Oak Ridge National Laboratory (ORNL) until now. Therefore, here, in this work, a preliminary assessment of several of these models was performed using openly available reflood experimental data from the Flooding Experiments in Blocked Arrays (FEBA) tests. One coarse mesh and one fine mesh model were set up in CTF for high and low flooding rate tests performed in the unblocked FEBA facility. A coarse TRACE model was set up to be as consistent as possible with the coarse CTF model to allow for code-to-code benchmarking. The assessment shows a tendency of the codes to over-predict peak cladding temperature (PCT) near the top of the bundle and to quench early. Advanced spacer grid models were shown to improve upper bundle predictions in CTF. The resolved CTF model over-predicted PCT by a larger degree in the center channels in the low-flooding rate test, and it is believed that the radiative heat transfer model, which was not used in this study, may be needed to correct this over-prediction. Finally, this work demonstrates the importance of the droplet model in determining quench time and vapor temperature and PCT prediction, which necessitates a more in-depth validation of these models in the future.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Comparison of experimental and simulation results on interior subchannels of a 61-pin wire-wrapped hexagonal fuel bundle

As part of a joint U.S. Department of Energy project, Framatome, Argonne National Laboratory, Terrapower, and Texas A&M University have collaborated to produce experimental data and computational results to characterize the flow behavior of a 61-pin wire-wrapped hexagonal fuel bundle. Here, comparisons are made between Texas A&M’s experimental particle image velocimetry velocity field measurements and Argonne National Laboratory’s large-eddy simulation results in an interior subchannel for a Reynolds number of 19,000. It can be concluded that the shape and magnitude of the mean vertical velocity component is in excellent agreement for the selected horizontal line plots. Maximum relative errors are less than 10% until the subchannel walls are approached. Finally, the shape and, to a lesser extent, magnitude of the mean horizontal velocity component is in satisfactory agreement between the particle image velocimetry and large-eddy simulation results.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Implementation of an Orificing Optimization Algorithm in the DASSH Subchannel Analysis Code

The Ducted Assembly Steady-State Heat transfer code (DASSH) performs full-core subchannel thermal hydraulics calculations in liquid metal fast reactors. One of the applications of subchannel codes is to optimize coolant flow orificing. As a design activity, the primary task is to determine the best way to divide assemblies into groups and distribute coolant flow rates among them. This report documents an algorithm implemented in DASSH to automatically optimize coolant orificing. Over the course of multiple iterations, DASSH determines the orifice grouping and flow distribution that minimizes peak coolant, clad, or fuel temperatures across all timesteps for a user-specified number of assembly groups. The total coolant flow rate in the reactor is constrained to achieve the specified core-average outlet temperature. The flow rate to each orifice group may also be constrained by the allowable pressure drop. The distribution of coolant flow among groups is accelerated using a predictor-corrector algorithm based on interpolated results from single-assembly parametric calculations. The assembly orificing grouping is initially predicted based on assembly power but can be refined if results demonstrate that an assembly would fit better in another group. The algorithm is demonstrated with two case studies. The first is a simple model for a reactor core consisting of just fuel assemblies; the pin power distributions are specified to create a situation where the initial assembly grouping prediction is suboptimal. This example is used to describe the initial grouping, demonstrate convergence over multiple iterations, and highlight the impact of regrouping. Then, the algorithm is applied to minimize peak clad and fuel temperatures in an example sodium-cooled fast reactor, the Versatile Test Reactor. The multicycle optimization confirms prior calculations for the reference core design. The example highlights how optimizing for different peak temperatures affects the results and demonstrates the use of the pressure drop constraint to limit the maximum flow rate.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Comparison Between Pin-by-Pin Subchannel and System Level Thermal Hydraulic Results for High Burnup Loss-of-Coolant Applications

This milestone report summarizes recent work to investigate higher fidelity modeling and simulation practices for large-break loss-of-coolant accident (LBLOCA) analysis in high-burnup pressurized water reactor (PWR) cores. Because of current industry interest in extending fuel cycle lengths ranging from 18 to 24-months, the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program has been investigating the susceptibility of high-burnup core designs to fuel fragmentation, release, and dispersal (FFRD) during accident conditions such as LBLOCA. This work has prioritized developing and demonstrating a methodology for calculating core-wide susceptibility to FFRD and addressing uncertainties identified in the analysis and the US Nuclear Regulatory Commission (NRC) Research Information Letter on FFRD. Part of this investigation seeks to quantify differences between system-level thermal hydraulic behavior and higher fidelity subchannel modeling methods to identify potential safety concerns or opportunities to minimize FFRD susceptibility. To this end, the NEAMS subchannel code, CTF, is being used, along with the NRC system analysis code TRACE, for analysis of LBLOCA in a core containing high-burnup fuel. This project includes two thrusts: (1) improving on the existing TRACE model for a 4-loop PWR for LBLOCA so that the model is higher fidelity and consistent with current USNRC best practices, and (2) using CTF to perform pin-resolved modeling of the core region of the reflood phase of a LBLOCA in a high-burnup PWR.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Oklo Sponsored Testing using PELICAN: System, Subchannel, and High-Fidelity Software Validation (Final CRADA Report)

Argonne National Laboratory (the Contractor), located in Lemont IL, and Oklo, Inc., (the Participant), headquartered in Santa Clara, CA, propose to enter into a Cooperative Research and Development Agreement (CRADA) to perform a gap analysis of thermal hydraulic data, perform prototypical fuel assembly pressure drop and cavitation model validation, generate the experimental data as well as the corresponding uncertainties for this matrix and, finally, develop the validation models with the Argonne system level code SAS4A/SASSYS-1, the Argonne subchannel analysis code DASSH, the Argonne high fidelity code Nek5000, and/or the Idaho National Laboratory code Pronghorn Subchannel. The work outlined below will significantly improve the experimental and validation database currently available for liquid metal fast reactors, thus making it a viable part of a comprehensive reactor design and licensing suite to be used by the participant.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Extension of Clad Damage Propagation Model for Fission Gas Dispersal and Two-Phase Flow Effects in MOOSE SubChannel Module

This report presents an extension of the Clad Damage Propagation (CDAP) model implemented in the MOOSE SubChannel Module (SCM) to capture post-failure fission-gas dispersal and two-phase flow effects in sodium-cooled fast reactor assemblies. The extended model tracks discharged gas axially and radially, computes channel-averaged flow quality and void fraction using a Lockhart–Martinelli framework, evaluates two-phase frictional pressure-drop multipliers, determines inlet mass-flow degradation under fixed core pressures, and applies an intensified-void-based heat-transfer degradation to affected fuel pins. Radial plume expansion is parameterized using mineral-oil jet experiments mapped to sodium conditions via Reynolds–Weber similarity. Implementation details are documented, along with the new methods and user inputs needed to control plume mapping and two-phase behavior. Demonstration simulations for 19- and 37-pin bundles show that breach size and inlet velocity strongly influence propagation potential: small breaches (≤0.5 mm) produce limited degradation while larger breaches (~1 mm) can drive oscillatory temperature spikes and enhanced failure propagation, especially at higher velocities. These results demonstrate that the extended CDAP model provides a more complete framework for quantifying cladding damage propagation and evaluating propagation potential in transient scenarios. The approach remains computationally efficient, consistent with subchannel-level analysis, yet incorporates sufficient physics to bridge localized post-failure effects with bundle- and assembly-scale degradation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Development of an advanced multi-scale subchannel model for CRUD risk assessments

AREVA NP has developed a process, based upon standard licensing methods and specialized applications, for the assessment of the risk of crud induced corrosion and crud induced axial off-set. This streamlined process allows for the evaluation of the crud risk for full core subchannel by subchannel models, by applying local conditions on each fuel rod/subchannel, with axial nodes spanning several centimeters. AREVA NP has applied these tools successfully to units from various NSSS vendors. To date, there have been no plant performance issues that have been outside the performance as predicted by the AREVA NP risk assessment tools. The study presented in this paper utilized the core simulator and thermal-hydraulic packages of this tool set to evaluate the effects of changing various fuel cycle design (e.g. radial zoning patterns) and fuel assembly design parameters (e.g. fuel rod diameter, spacer grid placement) on thermal-hydraulic key decision variables (KDVs) of interest. The study may be valuable as a first step towards optimizing the fuel cycle design, optimizing the fuel lattice design, or evaluating potential operating changes such as a power up-rate or modification of the cycle length for a given PWR plant. The various geometric inputs to the AREVA NP risk assessment process were tested by individually adjusting a single parameter and observing the effect. Multi-parameter studies are also included to assess the impact of optimized fuel cycle designs. The sensitivity of the model, in terms of thermal-hydraulic results, is presented. Provided a relationship between the thermal-hydraulic results to both crud deposition thickness and lithium uptake in the zirconium cladding oxide layer, approximate impacts to the risk of crud induced corrosion and crud induced axial off-set can be estimated. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Time-Resolved PIV/PTV Measurements on Interior Subchannels of a Wire-Wrapped 61-pin Hexagonal Fuel Bundle

Liquid metal fast reactor fuel bundles using sodium as a coolant typically utilize a tightly packed triangular lattice of fuel pins enclosed in a hexagonal duct. Often, the fuel pins are helically wrapped with a wire spacer. The primary function of the wire spacers is to maintain a gap between adjacent fuel pins. They also enhance sub-channel mixing to increase convective heat transfer and mitigate vortex-induced vibration. Texas A and M University has conducted isothermal flow experiments in a wire-wrapped 61-pin hexagonal fuel bundle to support the research on advanced nuclear fuel development sponsored by the Department of Energy (DOE). These tests produced high spatial and temporal resolution measurements of the velocity field and pressure at different locations in the fuel bundle. The experimental facility employs matched-index-of-refraction (MIR) techniques and laser diagnostic velocity measurement techniques such as particle image velocimetry (PIV), particle tracking velocimetry (PTV), and laser Doppler velocimetry (LDV). The result of this experimental activity will be an experimental database of pressure and velocity field measurements that will be suitable for validating computational fluid dynamics (CFD) methods, such as Reynolds-averaged Navier-Stokes (RANS) and highly resolved large eddy simulation (LES). PIV results on the set of edge subchannels near the hexagonal duct have been presented in another paper. This paper presents PIV results on two sets of interior subchannels, defined as Interior- 1 and Center-2, at a Reynolds number (Re) of 19,000 in the turbulent regime.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Steady-State Subchannel Thermalhydraulic Assessment of Thorium-Based Fuel Concepts for Use in Pressure Tube Heavy Water Reactors

Thorium is an abundant alternative fertile nuclear fuel that could potentially be implemented in pressure tube heavy water reactors (PT-HWR) in the near future. One of the key stages in developing the thorium-based fuels for PT-HWRs is the reactor core concept, which requires the collaboration of different disciplines, such as neutron physics, thermal-hydraulics, materials, chemistry, safety, and others. Currently there are numerous core concepts under investigation at the Canadian Nuclear Laboratories that involve several types of thorium-based fuels, and that could be implemented in a 700-MWe-class PT-HWR. Four fuel bundle concepts are considered, including i) the reference (or nominal) 37-element bundle; ii) a 37-element modified bundle, with the center element using a different fuel material; iii) a 35-element bundle design, and iv) an 18-element internally cooled annular fuel (ICAF-18) fuel bundle. This study presents preliminary results of the steady-state thermalhydraulic subchannel analyses of the hot channels of different PT-HWR core concepts using thorium-augmented or thorium-based fuels. A follow up study will be initiated to improve the assessment of the proposed cores and associated fuel bundle concepts for thorium-augmented and thorium-based fuels. The computational model will be revisited to include; - Expanding the analysis to cover different power profiles at different burnup stages of the bundle. - Expand the analysis to cover 28.5 kg/s. - Reduce the eccentricity of the 35-element fuel bundle. - Assess the effect of the phenomena not yet included in the subchannel model (such as the effect of cold wall, and gamma heat). - Expand the analysis to include the effect of spacers and junctions. - Compare bundle concepts using other performance indicators (such as burnup, breeding ratio), as optimizing the fuel bundle geometry is a multidisciplinary problem.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

A study on the impact of using a subchannel resolution for modeling of large break loss of coolant accidents

The nuclear industry is investigating the feasibility of transitioning from 18- to 24-month fuel cycles because of the positive impact it would have on the operational costs for the current fleet of light-water reactors. A challenge to making this change is the increased risk of fuel fragmentation, relocation, and dispersal (FFRD) due to the known potential for ceramic fuel to pulverize into fine particles at the higher discharge burnups. Previous work has been performed by the Nuclear Energy Advanced Modeling and Simulation program to assess FFRD risk in high-burnup cores using the BISON fuel performance code and a coarse mesh thermal hydraulics (T/H) solution for a loss-of-coolant accident (LOCA) using the TRACE system T/H code. Because of the importance of the T/H solution for FFRD assessment, this study seeks to investigate the impact of using higher-fidelity subchannel techniques for modeling of the LOCA transient. CTF was used to model a subregion of a high-burnup core that was depleted by the Virtual Environment for Reactor Applications (VERA) multiphysics core simulator. Both coarse-mesh and pin-resolved models were created in CTF, and a consistent coarse-mesh TRACE model was also developed to allow for benchmarking the code results. Further, a large-break loss-of-coolant accident (LBLOCA) reflood transient was simulated using these three models, and results were compared. Results showed some consistent differences between the CTF and TRACE coarse models, including a higher peak cladding temperature (PCT) prediction in CTF and later quenching in CTF; however, the transient clad temperature behavior was similar, and these differences are likely due to post-critical heat flux heat transfer modeling differences and minimum film boiling temperature model differences. The pin-resolved results indicate that the PCT in the lumped model is often under-predicted by as much as 70 °C and that PCT occurs at a different location than the high-power pin in the assembly. The lumped model predicts a difference of 10 °C or less between the average and hot pins in the assembly, whereas the pin-resolved model predicts a range of over 100 °C. These results indicate that higher-fidelity T/H results may have an impact on predicted core behavior during LOCA, which may be important to consider when assessing FFRD risk.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Advanced two-phase subchannel method via non-linear iteration

A fast-running, robust two-phase flow, sub-channel model is presented based on non-linear solution of the steady-state subchannel fluid flow equations. The drift-flux model solves for conservation of liquid and vapor mass, mixture energy, and axial and transverse mixture momentum as part of an efficient planar marching scheme and nonlinear, nested outer and inner iteration. Here, models based on mechanistic subcooled boiling, two-phase turbulent void mixing, and drift are included. Solution verification and mesh convergence studies were performed for modern GE 10 × 10 fuel geometry and are shown to have excellent convergence behavior. Run time performance for a 50 axial mesh model showed 2.2 seconds on a single CPU core to tightly converge all 3D distributions (flow, void, pressure) for the GE 10 × 10 fuel geometry, supporting its efficient use within the Virtual Environment for Reactor Applications boiling water reactor framework.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Implementation of a Spacer Grid Rod Thermal-Hydraulic Reconstruction (ROTHCON) Capability into the Thermal-Hydraulic Subchannel Code CTF

The Consortium for Advanced Simulation of Light Water Reactors is developing a core simulator capability known as the Virtual Environment for Reactor Applications (VERA) to address nuclear industry challenge problems such as crud-induced power shift (CIPS). The CTF thermal-hydraulic (T/H) subchannel code provides thermal feedback in the coupled neutronics, T/H, crud chemistry simulation that VERA performs. It has been discovered that the coarse meshing approach used by CTF (in which fuel rods are discretized into four azimuthal segments) can be a source of error in predicting crud growth and boron distribution in VERA CIPS calculations. Spacer grid effects lead to complex rod-to-fluid heat transfer behavior that, when not resolved, can lead to error in the prediction of crud growth and boron deposition. A higher-fidelity computational fluid dynamics approach can be used instead of CTF, but this leads to excessive simulation times. This paper presents an approach for using high-fidelity computational fluid dynamics data to create shape functions that are used in CTF to reconstruct rod surface heat transfer behavior as a function of spacer grid geometry. The approach is demonstrated for a 5 × 5 rod bundle facility with five mixing vane grids under a range of operating conditions encountered in nominal pressurized water reactor conditions. Furthermore, it is demonstrated that the grid heat transfer maps are successful at introducing a higher-fidelity heat transfer modeling capability into CTF.

CFD-informed↗

Multiscale Thermal-hydraulic analysis of the MARVEL micro-reactor using coupled MOOSE Subchannel (SCM) and SAM

MARVEL is a natural-convection-cooled sodium-potassium microreactor that is anticipated to generate 85 kilowatts of thermal energy. It will operate within Idaho National Laboratory’s Transient Reactor Test Facility and is being developed by the DOE Microreactor Program. MARVEL will be used to test microreactor applications, evaluate systems for remote monitoring, and develop autonomous control technologies. A thermal-hydraulic computational model of this facility is a valuable tool to study important transients and calculate the safety limits of the micro-reactor design. For this purpose, the authors have chosen to use a multiscale coupled simulation: SCM for modeling the reactor core and SAM for the reactor’s primary cooling system. SCM is MOOSE physics module for subchannel analysis, which was designed to model single-phase flows through liquid-metal cooled, wire-wrapped fuel pin sub-assemblies, ordered in a triangular lattice. The SCM code was modified to be able to model MARVEL’s unique geometry. SAM is a systems analysis module based on the MOOSE framework. It aims to provide fast-running, whole-plant transient analyses capability with improved-fidelity for various advanced reactor types. The coupling between the two SCM and SAM for MARVEL modeling is done implementing a domain over-lapping approach. The resulting coupled simulation can model transients such as reactor startup/shutdown and provide an intermediate fidelity picture of the temperature field and other variables, in the core. Results for the steady-state simulations are presented in the article as well as flow blockage transient.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

L3.PHI.CTF.P10.02-rev2 Coupling of Subchannel T/H (CTF) and CRUD Chemistry (MAMBA1D)

The purpose of this milestone is to create a preliminary capability for modeling light water reactor (LWR) thermal-hydraulic (T/H) and CRUD growth using the CTF subchannel code and the subgrid version of the MAMBA CRUD chemistry code, MAMBA1D. In part, this is a follow-on to Milestone L3.PHI.VCS.P9.01, which is documented in Report CASL-U-2014-0188-000, titled "Development of CTF Capability for Modeling Reactor Operating Cycles with Crud Growth". As the title suggests, the previous milestone set up a framework for modeling reactor operation cycles with CTF. The framework also facilitated coupling to a CRUD chemistry capability for modeling CRUD growth throughout the reactor operating cycle. To demonstrate the capability, a simple CRUD \surrogate" tool was developed and coupled to CTF; however, it was noted that CRUD growth predictions by the surrogate were not considered realistic. This milestone builds on L3.PHI.VCS.P9.01 by replacing this simple surrogate tool with the more advanced MAMBA1D CRUD chemistry code. Completing this task involves addressing unresolved tasks from Milestone L3.PHI.VCS.P9.01, setting up an interface to MAMBA1D, and extracting new T/H information from CTF that was not previously required in the simple surrogate tool. Speci c challenges encountered during this milestone include (1) treatment of the CRUD erosion model, which requires local turbulent kinetic energy (TKE) (a value that CTF does not calculate) and (2) treatment of the MAMBA1D CRUD chimney boiling model in the CTF rod heat transfer solution. To demonstrate this new T/H, CRUD modeling capability, two sets of simulations were performed: (1) an 18 month cycle simulation of a quarter symmetry model of Watts Bar and (2) a simulation of Assemblies G69 and G70 from Seabrook Cycle 5. The Watts Bar simulation is merely a demonstration of the capability. The simulation of the Seabrook cycle, which had experienced CRUD-related fuel rod failures, had actual CRUD-scrape data to compare with results. As results show, the initial CTF/MAMBA1D-predicted CRUD thicknesses were about half of their expected values, so further investigation will be required for this simulation.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

CTF (Subchannel) Calculations and Validation L3:VVI.H2L.P15.01

The goal of the Verification and Validation Implementation (VVI) High to Low (Hi2Lo) process is utilizing a validated model in a high resolution code to generate synthetic data for improvement of the same model in a lower resolution code. This process is useful in circumstances where experimental data does not exist or it is not sufficient in quantity or resolution. Data from the high-fidelity code is treated as calibration data (with appropriate uncertainties and error bounds) which can be used to train parameters that affect solution accuracy in the lower-fidelity code model, thereby reducing uncertainty. This milestone presents a demonstration of the Hi2Lo process derived in the VVI focus area. The majority of the work performed herein describes the steps of the low-fidelity code used in the process with references to the work detailed in the companion high-fidelity code milestone (Reference 1). The CASL low-fidelity code used to perform this work was Cobra Thermal Fluid (CTF) and the high-fidelity code was STAR-CCM+ (STAR). The master branch version of CTF (pulled May 5, 2017 – Reference 2) was utilized for all CTF analyses performed as part of this milestone. The statistical and VVUQ components of the Hi2Lo framework were performed using Dakota version 6.6 (release date May 15, 2017 – Reference 3). Experimental data from Westinghouse Electric Company (WEC – Reference 4) was used throughout the demonstrated process to compare with the high-fidelity STAR results. A CTF parameter called Beta was chosen as the calibration parameter for this work. By default, Beta is defined as a constant mixing coefficient in CTF and is essentially a tuning parameter for mixing between subchannels. Since CTF does not have turbulence models like STAR, Beta is the parameter that performs the most similar function to the turbulence models in STAR. The purpose of the work performed in this milestone is to tune Beta to an optimal value that brings the CTF results closer to those measured in the WEC experiments.

97 MATHEMATICS AND COMPUTING↗

Code and Solution Verification Assessment of the CTF Thermal Hydraulic Subchannel Code

CTF is a thermal-hydraulics subchannel code jointly developed by Oak Ridge National Laboratory and North Carolina State University. Over the past seven years, the Consortium for Advanced Simulation of Light Water Reactors (CASL) has made a significant investment in developing CTF so it can be used to model light water reactors, including nominal operating conditions, departure from nucleate boiling analysis, and transients ranging from loss of flow to reactivity insertion accidents. In addition to implementing new modeling capabilities and developing the user input and output interface, extensive work has been performed to improve the code’s quality assurance program, resulting in a development process that conforms with NQA-1 requirements. The CASL program follows the Predictive Capability Maturity Model (PCMM) approach for assessing code quality, which emphasizes performing code verification(ensuring the code converges to the correct answer) and solution verification (ensuring the code converges for the intended application). Code and solution verification are used to identify uncertainty errors introduced by numerical approximations in the code and are important for demonstrating that the model is coded without error, which is an important aspect of the Best Estimate plus Uncertainty method. This paper presents a comprehensive overview of the code and solution verification testing that has been performed on CTF. A top-down approach is taken in which the intended CTF applications are presented, followed by the code features required for their modeling. These features are then linked to the applicable code and solution verification tests that demonstrate proper functioning. Past testing efforts are summarized, and new tests are added to help close gaps in the presented test matrix. Rather than performing “one-off” exercises, these tests are added to the automated CTF regression test suite to ensure continual code quality.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Development and Integration of a Stochastic Clad Damage Propagation Model into PRONGHORN-SC Subchannel Analysis Code

The failure of fuel pins in nuclear reactors is intrinsically stochastic. Typically, a combination of variation in manufacturing that affects the material characteristics and the fuel assembly dimensions, variation in operating conditions, such as local power, coolant flow rate, and irradiation induced changes in material properties lead to a large uncertainty in failure margin of the fuel pins. Failure, therefore, may occur in exceptional pins with adverse combinations of these variations. Upon a metal fuel pin (U-Pu-Zr/HT9) failure, depressurization of the fuel pin takes place by release of fission gas, liquid sodium bond, and potentially solid fuel particles or molten/eutectic fuel droplets through the hole in cladding. The effect of a fission gas jet on neighbor fuel pins and possible propagation of a clad damage during normal operation was studied experimentally in 1970s and it was found that the post-failure fission gas jet insulates the jet impingement area of the target fuel pin surface and could increase the target pin’s surface temperature by as much as 100 – 200 K during the failed pin depressurization. It was concluded that the effect should not lead to fuel pin failure propagation during normal operation. In accident scenarios of sodium and lead fast reactors such as Unprotected Loss-Of-Flow (ULOF) or Unprotected Transient Over Power (UTOP), the fuel pins can be subjected to higher clad temperatures and fuel pin pressures or fuel clad mechanical/chemical interaction where thermal creep margin becomes significantly lower compared to the normal operation conditions. Therefore, possible stochastic failure and the post-failure fission gas/fuel jet impingement could be critical in order to predict fuel pin failure propagation. Pin depressurization due to fission gas release may degrade the heat transfer by formation of a gas blanket on a neighboring pin surface, which is a local phenomenon, and by causing coolant flow deceleration and starvation, which could affect a surrounding region as well. Furthermore, the potential presence of solid fuel particles or molten fuel at the time of clad failure could boost post-failure jet induced degradation even further. The present study models the U-Pu-Zr/HT9 metal fuel pin failure and stochastic clad damage propagation by biased sampling based on a Cumulative Damage Fraction (CDF) type clad failure criterion and the normal distribution of fuel failure probability density as a function of logarithm of Cumulative Damage Fraction. In addition, the effect of post-failure fission gas jet on heat transfer degradation is modeled for the target pins. This model is called stochastic Clad Damage Propagation (CDAP). The CDAP model is now fully integrated into developmental version of PRONGHORN-SC subchannel analysis code, allowing for modeling local failures and its propagation potential. Section 2 describes the components of the CDAP models. Section 3 describes the model implementation to PRONGHORN-SC and input specifications. Section 4 describes the CDAP model validation coupled to PRONGHORN-SC.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗