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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Fluid Flow Investigations within a 37 Element CANDU Fuel Bundle Supported by Magnetic Resonance Velocimetry and Computational Fluid Dynamics

The current work presents experimental and computational investigations of fluid flow through a 37 element CANDU nuclear fuel bundle. Experiments based on Magnetic Resonance Velocimetry (MRV) permit three-dimensional, three-component fluid velocity measurements to be made within the bundle with sub-millimeter resolution that are non-intrusive, do not require tracer particles or optical access of the flow field. Computational fluid dynamic (CFD) simulations of the foregoing experiments were performed with the hydra-th code using implicit large eddy simulation, which were in good agreement with experimental measurements of the fluid velocity. Greater understanding has been gained in the evolution of geometry-induced inter-subchannel mixing, the local effects of obstructed debris on the local flow field, and various turbulent effects, such as recirculation, swirl and separation. These capabilities are not available with conventional experimental techniques or thermal-hydraulic codes. Finally, the overall goal of this work is to continue developing experimental and computational capabilities for further investigations that reliably support nuclear reactor performance and safety.

42 ENGINEERING↗
DOE OSTI2020

Development and Implementation of a CTF Code Verification Suite.

CTF is a thermal hydraulic subchannel code developed to predict light water reactor (LWR) core behavior. It is a version of Coolant Boiling in Rod Arrays (COBRA) developed by Oak Ridge National Laboratory (ORNL) and North Carolina State University (NCSU) and used in the Consortium for the Advanced Simulation of LWRs (CASL). Here, the existing CTF code verification matrix is expanded, which ensures that the code is a faithful representation of the underlying mathematical model. The suite of code verification tests are mapped to the underlying conservation equations of CTF and significant gaps are addressed. As such, five new problems are incorporated: isokinetic advection, conduction, pressure drop, convection, and pipe boiling. Convergence behavior and numerical errors are quantified for each of the tests and all tests converge at the correct rate to their corresponding analytic solution. A new verification utility that generalizes the code verification process is used to incorporate these problems into the CTF automated test suite.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

CFD analysis of a solid pin-fueled small modular fluoride salt-cooled reactor

A computational fluid dynamics (CFD) study is performed for the pre-conceptual solid pin-fueled small modular fluoride salt-cooled reactor developed by the Oak Ridge National Laboratory (ORNL). The solid and fluid regions in a 1/12th section of a fuel assembly are modeled. The solid region includes fuel and non-fuel pins, with heat generation in the active region of the fuel pins. Two different power profiles, uniform and center-peaked, are considered. The FLiBe coolant flows from the bottom to the top of the core, parallel to the bank of fuel and non-fuel pins. The k-ω shear-stress transport (SST) model is chosen as the baseline turbulence model. The effects of grid refinement, inlet turbulence specification, and turbulence models on the temperature and pressure drop predictions are studied. Turbulence model sensitivity is investigated by comparing the results from the k-ω SST model with other two-equation models (k-ω baseline or BSL and k-ε realizable) as well as with anisotropic Reynolds stress transport models (linear pressure-strain and stress-BSL). Here the results show that the choice of turbulence model has a significant impact on the pin temperature and bundle pressure drop predictions. The results from the baseline turbulence model show good agreement between the fuel pin temperatures and bundle pressure drop values predicted using a subchannel model previously developed by the authors of the present study.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

CTF: A modernized, production-level, thermal hydraulic solver for the solution of industry-relevant challenge problems in pressurized water reactors

CTF is a thermal hydraulic (T/H) subchannel tool that has been extensively developed over the past ten years as part of the Consortium for Advanced Simulation of Light Water Reactors (CASL) program. The code was selected early in the CASL program for support of high-impact challenge problems that were found to be relevant to the nuclear industry and its currently operating fleet of pressurized water reactors (PWRs), including issues such as departure from nucleate boiling (DNB), crud-induced power shifts (CIPSs), and reactivity-insertion accidents (RIAs). By incorporating CTF into the multiphysics Virtual Environment for Reactor Application (VERA) core simulator software developed by CASL, CTF has become the primary means of providing fluid and fuel thermal feedback, as well as T/H figure-of-merits (FOMs) in large-scale reactor simulations. With the goal of solving industry challenge problems, CASL placed great emphasis on developing high-quality, high-performance, validated software tools that offer higher fidelity than what is currently possible with current industry methods. In support of this effort, CTF was developed from a research tool into an nuclear quality assurance (NQA-1)–compliant, production-level software tool that is capable of addressing the stated challenge problems and goals of CASL. This work presents a review of the major technological achievements that were realized in developing CTF over the past decade of the CASL program and presents an overview of the code solution approach and closure models.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Leveraging MARVEL and SPHERE to Demonstrate NEAMS Thermal Hydraulics Codes

The Single Primary Heat Extraction and Removal Emulator (SPHERE) and Microreactor Applications Research Validation and Evaluation (MARVEL) were utilized to demonstrate thermal hydraulics codes of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. A recent long-duration test of a sodium heat pipe conducted at the SPHERE facility was utilized for the validation of the heat pipe code Sockeye. The first 10 hours of operation featured a startup of the heat pipe from room temperature; this portion was modeled with Sockeye and compared to the experimental data. Measurements included several thermocouple measurements across the length of the heat pipe surface. Two of Sockeye’s heat pipe models were applied to experiment, as well as several models of the cooling at the condenser end of the heat pipe. Sockeye results gave a decent match to the experimental data, with notable discrepancy in the condenser section, where experimental data indicated a significant inactive length of the heat pipe. The MARVEL microreactor was modeled through a coupled simulation of System Analysis Module (SAM) and SubChannel Module (SCM). This multiscale coupled simulation was performed as an exercise to demonstrate the compatibility and applicability of these codes for thermal-hydraulic analysis of microreactors. The coupling between the SCM and SAM codes for multiscale modeling of MARVEL, was achieved through a domain overlapping approach. In this example a transient SAM simulation provided boundary conditions to SCM, and SCM provided pressure drop to SAM.

97 - MATHEMATICS AND COMPUTING↗
DOE OSTI2016

Implementation of the LAX-Wendroff Method in Cobra-TF for Solving Two-Phase Flow Transport Equations

COBRA-TF (Coolant Boiling in Rod Arrays Two Fluid), or CTF, is a subchannel code used to conduct the reactor core thermal hydraulic (T/H) solution in both standalone and coupled multi-physics applications. CTF applies the first-order upwind spatial discretization scheme for solving two-phase flow conservation equations. In this work, the second-order Lax-Wendroff (L-W) scheme has been implemented in CTF to solve the two-phase flow transport equations to improve numerical accuracy in both temporal and spatial discretization. To avoid the oscillation issue, a non-linear flux limiter VA (Van Albada) is employed for the convective terms in the transport equations. Assessments have been carried out to evaluate the performance and stability of the implemented second-order L-W scheme. It has been found that the L-W scheme performs better than the upwind scheme for the single-phase and two-phase flow problems in terms of numerical accuracy and computational efficiency.

CTF↗
DOE OSTI2021

Validation of Oregon State University High Temperature Test Facility Experiments Using Pronghorn

The OSU High Temperature Test Facility is a quarter-scale diameter, 1/64 scale volume test facility meant to replicate thermophysical phenomena in the prototypical General Atomics Modular High Temperature Gas Reactor. Tests pertaining to conduction cooldown events were performed from 2016-2019, providing a large database by which computational methods that are applicable to different length scales can be validated. One of these codes is Pronghorn, which is a coarse-meshed, porous-based subchannel thermal hydraulics code based on the MOOSE application, with the intention of better capturing the physics of both conduction and convection heat transfer within the OSU HTTF core. The goal of this summer project is to develop the framework by which Pronghorn can perform validation exercises of the HTTF core for benchmarking, by generating a mesh appropriate to the geometry of the HTTF core, developing input decks that accurately capture the initial and boundary conditions, materials, and relevant equations to the physics seen in the HTTF core, and using a postprocessor to compare simulation results to various experimental data. While validation of codes is a multi-year project, a mesh has been generated and tested in Pronghorn that meets mass conservation and basic heat transfer principles. The next step is to accurate depict the fluid inlet and outlet boundary conditions, which will be performed using computational fluid dynamics software.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

SCM overview & the EBR-II shutdown heat removal tests validation.

Pronghorn is an engineering-scale, coarse-mesh, thermal-hydraulics tool for supporting reactor-core simulations of advanced nuclear reactors. Most of the current efforts in Pronghorn have been devoted in developing porous finite-volume capabilities and adapting closure correlations for coarse-mesh thermal-hydraulics modeling. However, for liquid-metal reactors (LMRs) with wire-wrapped fuel pin assemblies, a pin-level thermal-hydraulic resolution is required for most safety case studies (pin rupture, channel blockage, etc.). For this purpose, a new Subchannel application is developed in MOOSE, which affords the required flow field resolution, while still preserving an engineering-scale approach. This new solver can be natively coupled to Pronghorn and other MOOSE objects to enable full-core, multi-physics, multi-scale engineering studies. This presentation presents the main features of the SCM code and demonstrates a validation case based on the EBR-II SHRT tests.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗
DOE OSTI2019

Numerical exploration of three relativistic particles in a finite volume including two-particle resonances and bound states

In this work, we use an extension of the quantization condition, given in ref. [1], to numerically explore the finite-volume spectrum of three relativistic particles, in the case that two-particle subsets are either resonant or bound. The original form of the relativistic three-particle quantization condition was derived under a technical assumption on the two- particle K matrix that required the absence of two-particle bound states or narrow two- particle resonances. Here we describe how this restriction can be lifted in a simple way using the freedom in the definition of the K-matrix-like quantity that enters the quantization condition. With this in hand, we extend previous numerical studies of the quantization condition to explore the finite-volume signature for a variety of two- and three-particle interactions. We determine the spectrum for parameters such that the system contains both dimers (two-particle bound states) and one or more trimers (in which all three particles are bound), and also for cases where the two-particle subchannel is resonant. We also show how the quantization condition provides a tool for determining infinite-volume dimer- particle scattering amplitudes for energies below the dimer breakup. We illustrate this for a series of examples, including one that parallels physical deuteron-nucleon scattering. All calculations presented here are restricted to the case of three identical scalar particles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

A simulation of I 2 S-LWR selected transients

The Integral Inherently Safe (I 2 S-LWR) reactor was designed to improve the safety performance of PWR types reactors, and several design basis accidents were eliminated using innovative approaches to the I 2 S-LWR plant layout. However, some accidents such as station blackout (SBO), main steam line break (MSLB) type accident – the secondary hot leg break (SHLB) – were deemed statistically relevant, and the work performed here was to simulate these events using the system thermal–hydraulics code RELAP5 coupled to the 3D spatial kinetics code PARCS. The safety response of the reactor to a station blackout is predicted to be acceptable at both beginning of cycle (BOC) and end of cycle EOC core conditions with either the UO 2 or U 3 Si 2 fuel options considered here. In the analysis of the SHLB performed here, an unprotected (no SCRAM) condition was assumed to highlight the I 2 S-LWR capabilities. Results indicate that with the UO 2 core, the fuel temperatures rise notably, but no fuel melting was predicted and departure from nucleate boiling ratio (DNBR) limits were not exceeded. Similarly, the U 3 Si 2 fuel performed well, but with lower fuel temperatures and decreased DNBR. Following the break of the steam line in the I 2 S-LWR, shutdown is achieved after about two minutes via negative reactivity from loss of primary flow. In all cases considered, the I 2 S-LWR appears to accomplish the goals of increased safety performance. Further design improvements and other fuel options may marginally improve on the observed safety performance. Additionally, higher fidelity analysis with the next generation subchannel, Computational Fluid Dynamics (CFD), or full core transport may be useful in refining the prediction of the peak fuel and cladding temperatures, as well limiting DNBR during the accidents.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2021

Comprehensive framework for data-driven model form discovery of the closure laws in thermal-hydraulics codes

The two-phase two-fluid model is a basis of many thermal-hydraulics codes used in design, licensing, and safety considerations of nuclear power plants. Thermal-hydraulics codes rely on the closure laws to close the system of conservation equations and describe the interactions between phases. These laws, derived from years of experimental investigations, are semi-empirical correlations that lack generality and have a limited range of applicability. Increase of computational power, availability of new experiments, and development of high-fidelity simulations has increased the number of validation data. The discrepancies between the code predictions and the validation data are a great source of knowledge. Missing physics that are not included in the model but are important for the considered phenomena can be discovered by propagating the information from the experimental results through the model. Furthermore, physics-discovered data-driven model form (P3DM) methodology integrates available integral effect tests and separate effects tests to determine the necessary corrections to the model form of the closure laws. In contrast to existing calibration techniques, the methodology modifies the functional form of the closure laws. Based on the functional form of the correction, the missing physics that were not included in the original model can be discovered. The methodology provides the alternative to the machine learning approach, in which the model is discovered in the form of the intractable black-box relation. In this work, the methodology was applied to the CTF subchannel code to improve the prediction of the two-phase flow phenomena.

42 ENGINEERING↗
DOE OSTI2016

Modeling and simulation challenges pursued by the Consortium for Advanced Simulation of Light Water Reactors (CASL)

The Consortium for the Advanced Simulation of Light Water Reactors (CASL), the first Energy Innovation Hub of the Department of Energy, was established in 2010 with the goal of providing modeling and simulation (M&S) capabilities that support and accelerate the improvement of nuclear energy's economic competitiveness and the reduction of spent nuclear fuel volume per unit energy, and all while assuring nuclear safety. To accomplish this requires advances in M&S capabilities in radiation transport, thermal-hydraulics, fuel performance and corrosion chemistry. To focus CASL's R&D, industry challenge problems have been defined, which equate with long standing issues of the nuclear power industry that M&S can assist in addressing. To date CASL has developed a multi-physics “core simulator” based upon pin-resolved radiation transport and subchannel (within fuel assembly) thermal-hydraulics, capitalizing on the capabilities of high performance computing. CASL's fuel performance M&S capability can also be optionally integrated into the core simulator, yielding a coupled multi-physics capability with untapped predictive potential. Material models have been developed to enhance predictive capabilities of fuel clad creep and growth, along with deeper understanding of zirconium alloy clad oxidation and hydrogen pickup. Understanding of corrosion chemistry (e.g., CRUD formation) has evolved at all scales: micro, meso and macro. CFD R&D has focused on improvement in closure models for subcooled boiling and bubbly flow, and the formulation of robust numerical solution algorithms. For multiphysics integration, several iterative acceleration methods have been assessed, illuminating areas where further research is needed. Finally, uncertainty quantification and data assimilation techniques, based upon sampling approaches, have been made more feasible for practicing nuclear engineers via R&D on dimensional reduction and biased sampling. Industry adoption of CASL's evolving M&S capabilities, which is in progress, will assist in addressing long-standing and future operational and safety challenges of the nuclear industry. - Highlights: • Complexity of physics based modeling of light water reactor cores being addressed. • Capability developed to help address problems that have challenged the nuclear power industry. • Simulation capabilities that take advantage of high performance computing developed.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2016

Best estimate plus uncertainty analysis of departure from nucleate boiling limiting case with CASL core simulator VERA-CS in response to PWR main steam line break event

VERA-CS (Virtual Environment for Reactor Applications, Core Simulator) is a coupled neutron transport and thermal-hydraulics subchannel code under development by the Consortium for Advanced Simulation of Light Water Reactors (CASL). VERA-CS was used to simulate a typical pressurized water reactor (PWR) full core response with 17x17 fuel assemblies for a main steam line break (MSLB) accident scenario with the most reactive rod cluster control assembly stuck out of the core. The accident scenario was initiated at the hot zero power (HZP) at the end of the first fuel cycle with return to power state points that were determined by a system analysis code and the most limiting state point was chosen for core analysis. The best estimate plus uncertainty (BEPU) analysis method was applied using Wilks’ nonparametric statistical approach. In this way, 59 full core simulations were performed to provide the minimum departure from nucleate boiling ratio (MDNBR) at the 95/95 (95% probability with 95% confidence level) tolerance limit. The results show that this typical PWR core remains within MDNBR safety limits for the MSLB accident.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

Application of regression, variance, and density based global sensitivity methods to integrated VERA-CS and BISON simulations [Application of regression, variance, and density based global sensitivity methods to integrated VERA-CS and BISON simulations]

Uncertainty quantification (UQ) and sensitivity analyses (SA) are performed for coupled simulations between VERA-CS, a coupled pin resolved neutron transport and subchannel thermal hydraulics code, and the fuel performance code BISON. The interface between VERA-CS and BISON is performed in a multiphysics environment known as the LOCA Toolkit for U.S. light water reactors (LOTUS) currently under development at Idaho National Laboratory (INL). A focus is placed on using a variety of SA measures, including two regression based (Pearson and Spearman), one variance based (Sobol indices), and three moment independent measures (Delta moment independent measures with L 1 , L 2 , and L ∞ norms). The problem under inspection is a single assembly depletion case for three fuel cycles. The figures of merit are the minimum departure from nucleate boiling ratio (MDNBR), maximum fuel centerline temperature (MFCT), and gap conductance at peak power (GCPP). SA results show MDNBR to be linear with consistent rankings throughout the fuel cycles. MFCT is linear, but with a change in rankings at the switch from open gap to closed gap models. GCPP is nonlinear at intermediate states that coincide with the onset of contact between fuel and cladding. Furthermore, these nonlinear states allow for the showcasing of higher order SA measures over first order methods.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2018

A new fuel modeling capability, CTFFuel, with a case study on the fuel thermal conductivity degradation

Here, a new fuel modeling capability, CTFFuel, is developed from the subchannel code, CTF. This code is a standalone interface to the CTF fuel rod models, allowing for fuel rod simulations to be run independently from the fluid. This paper provides an overview of the code with a case study on the thermal conductivity degradation of LWR fuels to demonstrate its capabilities. The modeling of fuel thermal conductivity degradation in the code is improved through the addition of new modeling options to account for the irradiation effects via globally defined parameters. After the initial implementation, a variety of order-of-accuracy tests and code comparisons are performed to test software quality. A controlled analysis is allowed by CTFFuel to verify the numerical scheme of CTF’s conduction solution and to benchmark its fuel temperature predictions against FRAPCON-4.0’s. Overall, the software quality and verification procedure ensures that the new model is coded correctly, that it properly interacts with the rest of the code.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2019

CTF-R: A novel residual-based thermal hydraulic solver

The traditional scientific process has been revolutionized by the advent of computational modeling, but the nuclear industry generally uses "legacy codes," which were developed early in the evolution of computers. One example of a legacy code, the thermal hydraulic subchannel code CTF, is modernized in this work through the development of a novel residual-based version, CTF-R. Unlike its predecessor, CTF-R is not limited by the strict computational limitations of the early 1980's, and can therefore be designed such that it is inherently flexible and easy to use. A case study is examined to demonstrate how the flexibility of the code can be used to improve simulation results. In this example, the coupling between the solid and liquid fields is examined. Traditionally, this coupling is modeled explicitly, which imposes numerical stability limits on the time step size. These limits are derived and it is shown that they are removed when the coupling is made implicit. Further, the development of CTF-R will enable future improvements in next generation reactor modeling, numerical methods, and coupling to other codes. Through the further development of CTF-R and other residual-based codes, state-of-the-art simulation is possible.

42 ENGINEERING↗
DOE OSTI2019

Comparison of experimental and simulated critical heat flux tests with various cladding alloys: Sensitivity of iron-chromium-aluminum (FeCrAl) to heat transfer coefficients and material properties

In this paper we analyze the differences between previous transient critical heat flux (CHF) experiments using iron-chromium-aluminum (FeCrAl), Inconel 600, and stainless steel 316 (SS316) alloy test sections and best-estimate modeling results from widely-used nuclear engineering systems and subchannel analysis tools. FeCrAl is an Accident Tolerant Fuel (ATF) candidate cladding material. The thermal hydraulic performance and safety characteristics of FeCrAl are being evaluated to determine viability as a cladding material in Light Water Reactors (LWRs). In this study, the results of the CHF experiments conducted at atmospheric pressure and fixed inlet coolant temperature and mass flux are compared to models built in the fifth version of the Reactor Excursion and Leak Analysis Program (RELAP5-3D) and CTF, the modernized version of COBRA-TF developed by the Consortium for Advanced Simulation of LWRs (CASL). Results from RELAP5-3D and CTF showed differences from the experiments and from each other in predicting CHF. In the Inconel 600 case, both computational tools overpredicted CHF, which led to an underprediction in the tube outer surface temperature. In the SS316 and FeCrAl cases, CHF was underpredicted by the codes, leading to an overprediction of the tube outer surface temperature. To understand the discrepancies in CHF and post-CHF predictions, studies were performed using RELAP5-3D and RAVEN to determine the sensitivity of CHF and peak test section temperature, an analog to peak cladding temperature (PCT), to heat transfer coefficients, a CHF multiplier, and uncertainties in the thermal conductivity and volumetric heat capacity. We found that CHF depends most strongly on the CHF multiplier and thermophysical properties. A combination of these factors that produced the best match to the experiment based on CHF, PCT, and the total energy deposited into the tube was determined. The best match parameters were able to provide best-estimate predictions of the CHF and integral heat flux, but were still conservative when predicting the PCT. The best match set of parameters developed in this paper are intended only as a demonstration of an approach that could be applied in the future with a larger set of experiments to produce more accurate models of CHF and post-CHF behavior.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

One-way coupled simulation of FIV in a 7-pin wire-wrapped fuel pin bundle

A normal fuel geometry found in sodium cooled reactor designs is the hexagonal lattice of wire-wrapped fuel pins. The wire-wrappers, which replace the more common spacer grids, act as a spacer between the fuel pins, while simultaneously increasing the fluid mixing between subchannels. The induced transverse flow, while good for mixing, increases the force load experienced by the fuel pins, increasing the possibility of problematic flow induced vibrations. The primary source of vibrations is turbulent buffeting, so very small in amplitude, but capable of rupturing the fuel cladding by fretting. The vast majority of research on fretting has been for bare rods with spacer grids. There are several significant ways that wire-wrapped fuel pins differ from bare rods that may have significant impact on the vibrations present. The increased transverse flow and asymmetrical flow pattern may result in larger and imbalanced force loads on the pins. The helical bracing of each pin may result in unique vibration patterns, and gaps may form between pins affecting how they are braced. Simulation allows for the observation of the forces present, how the flow field affects these forces, and how the structure responds and can be repeated for multiple cases for relatively low cost. Here, simulation of the flow induced vibrations for a 7-pin wire-wrapped fuel pin bundle is carried out by coupling Nek5000, a spectral element LES/DNS solver with DIABLO, a finite element CSM solver. The small amplitude vibrations allow for the assumption that the structure motion has negligible effect on the fluid flow, permitting the much less computationally expensive one-direction coupling. In order to gather significant force history data for the 2.63 helical pitch structure, a scheme was developed to artificially create force history from the power spectrum density of the one-pitch CFD simulation forces. The vibrations could then be simulated in the structure multiple times varying the contact and gaps between various pins. This provides for the search of worst-case scenarios and informs about the types of fretting that might be present. The simulation geometry is based on an experiment for comparison. However, the viscosity of the simulation was increased by a factor of four to make the simulation computationally feasible. While validation is not possible, comparisons are still valuable.

42 ENGINEERING↗