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 127 records · Page 7

Improved departure from nucleate boiling prediction in rod bundles using a physics-informed machine learning-aided framework

The critical heat flux (CHF) corresponding to the departure from nucleate boiling (DNB) crisis is a regulatory limit for the licensing of pressurized water reactors (PWRs) worldwide. Despite the abundance of predictive tools available to the reactor thermal-hydraulics community, the path for an accurate CHF model remains elusive. This work approaches the prediction of DNB through a physics-informed machine learning-aided framework (PIMLAF) with the objective of achieving superior predictive capabilities for a rod bundle. In view of the limitations in existing macro-scale physics-driven tools, an improved mechanistic model is first proposed, leveraging key concepts in the liquid sublayer dryout and bubble crowding mechanisms. Furthermore, the proposed mechanistic model is able to predict DNB in different heater geometries for a broad range of flow conditions without the need for recalibration. This model is then incorporated as the physics-informed component of the hybrid framework PIMLAF, which takes advantage of established understanding in the field (i.e., domain knowledge [DK]) and uses machine learning (ML) to capture undiscovered information from the mismatch between the actual and DK-predicted output. Two bundle-related case studies using the PWR subchannel and bundle tests (PSBT) database are carried out to illustrate the PIMLAF’s improved performance over traditional approaches for both interpolation and extrapolation purposes. In light of the PIMLAF’s promising potential to reduce prediction error, reactor vendors are encouraged to leverage their in-house experimental efforts and apply the hybrid framework to potentially achieve margin reductions in the minimum DNB ratio (MDNBR) for the designs of interest.

42 ENGINEERING↗

Demonstration of RANS models with wall functions in the spectral element code Nek5000

The spectral element based computational fluid dynamics (CFD) code Nek5000 has been traditionally used for high-fidelity applications, such as direct numerical simulation (DNS) and large eddy simulation (LES). These techniques require very fine numerical resolution to accurately capture turbulent fluctuations which can be prohibitively expensive for users without access to leadership class computing facilities. For broader application and adoption, significant effort has been invested to develop Reynolds-averaged Navier–Stokes (RANS) capabilities in Nek5000. Here, this work presents details of the implementation and demonstration of the standard wall functions for the κ–τ model in Nek5000. Results using the wall-modeled approach are compared to a wall-resolved approach for cases with negligible pressure gradient, viz., channel flow, pipe flow and flow in a reactor subchannel. Results show reasonably good agreement between the two approaches for friction factor and Nusselt number. Some expected differences are identified near the wall. These cases demonstrate the potential for significant computational savings by using much coarser meshes for the wall-modeled approach, with only minor differences between the predicted result. Additionally, several Reynolds numbers up to 1,000,000 are demonstrated for pipe flow and predicted friction factors and Nusselt numbers compared well to available correlations, with the worst below 10%. As the Reynolds number is increased, better agreement is observed between the correlations and the wall-modeled approach. In addition, flow in a molten salt fast reactor (MSFR) core is considered which features an adverse pressure gradient and flow separation. It showcases the inability of standard wall functions to accurately predict flows with adverse pressure gradients. The results, however, match reasonably well in trend in regions of the flow where the boundary layer is attached. Ongoing research is dedicated to include a pressure gradient correction to wall functions to improve the accuracy of flows with separation or reattachment and adverse or favorable pressure gradients.

42 ENGINEERING↗

Evaluation of RANS vs. LES simulation of fluid flow through 3 × 3 rod bundle with a simple spacer grid as a precursor to coupled fluid–structure interaction simulations

The research literature on Computational Fluid Dynamics (CFD) of coolant flow through rod bundles with spacer-grids and mixing vanes is replete, ranging from high fidelity Large Eddy Simulation (LES)/Direct Numerical Simulation (DNS) simulations to Reynolds-Averaged Navier–Stokes (RANS) modeled studies. The mixing of flow between subchannels and the pressure drop through the bundle are fundamental quantities useful for comparing and evaluating CFD methods. Less commonly observed and compared are the forces exerted onto the structure by the fluid. The present study seeks to evaluate the use of RANS simulations for predicting the structural response to fluid flow. Wall resolved RANS simulations are benchmarked against LES simulations of fluid flow at a Reynolds number of 15,000 through a 3 × 3 fuel rod bundle with a simple spacer grid. Velocity line-plots are compared showing good agreement between RANS and LES results, ascertaining that the former is capable of capturing the essential time-averaged velocity profile. Additionally, the distribution of forces on the spacer grid and fuel rods are collected as a function of time and space. The RANS methods are evaluated using the frequency and magnitude of the fluctuating forces on various portions of the structure as compared to LES. In conclusion, the power spectral density evaluation of the models reveal underprediction of force amplitude on the rod walls by RANS and also discrepancy in the prediction of high frequency spectra, especially in the immediate vicinity of spacer-grid structure, which may be attributed to the lack of random turbulence fluctuation or insufficient modeling of small-scale eddies in RANS simulation.

FIV↗

Large Eddy Simulation of Low-Reynolds-Number Turbulent Flow of Low-Prandtl-Number Fluid in a Tight Lattice Bundle for Assessment of Reynolds-Averaged Navier-Stokes Turbulence Model

The MARVEL (Microreactor Applications Research Validation and Evaluation) microreactor utilizes natural circulation as core cooling mechanism and liquid metal as a primary coolant. Moreover, the reactor core has a pitch-to-diameter ratio of 1.056, which is considered a tight lattice configuration. Numerous studies have widely reported that Reynolds-Averaged Navier-Stokes (RANS) turbulence models inaccurately predict heat transfer in liquid metals and fail to capture flow pulsations that can occur within tight lattices, leading to further inaccuracies in simulation results. Therefore, evaluating the accuracy of RANS turbulence models in the thermal-hydraulic analysis of the MARVEL microreactor core is crucial for assessing reactor safety. In this study, a Large Eddy Simulation (LES) of the MARVEL microreactor core subchannel was conducted and compared with a RANS simulation to evaluate the accuracies and conservatism of the RANS model. The flow pulsation in a tight lattice predicted by LES enhanced the heat transfer, whereas the RANS model underpredicted it. Consequently, the RANS model predicted the peak cladding temperature higher than the LES model, but the discrepancy between the two approaches was not significant due to the good thermal characteristics of the liquid metal. It can be concluded that the steady-state RANS model is effective for the thermal analysis of liquid-metal-cooled MARVEL microreactor core and can provide conservative predictions from a safety analysis perspective.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

Uncertainty Quantification Study of CTF for the OECD/NEA LWR Uncertainty Analysis in Modeling Benchmark

This work describes the results of a quantitative uncertainty analysis of the thermal-hydraulic subchannel code for nuclear engineering applications, Coolant Boiling in Rod Arrays-Three Field (COBRA-TF). CTF is used, which is a version of COBRA-TF developed in cooperation between the Consortium for Advanced Simulation of Light Water Reactors and North Carolina State University. Four steady-state cases from Phase II Exercise 3 of the Organisation for Economic Co-operation and Development/Nuclear Energy Agency Light Water Reactor Uncertainty Analysis in Modeling (UAM) Benchmark are analyzed using the statistical analysis tool, Design Analysis Kit for Optimization and Terascale Applications (Dakota). The input parameters include boundary condition, geometry, and modeling uncertainties, which are selected using a sensitivity study and then defined based on expert judgment. Here, a forward uncertainty quantification method with Latin hypercube sampling (LHS) is used, where the sample size is based on available computational resources.

42 ENGINEERING↗

A Code-Agnostic Driver Application for Coupled Neutronics and Thermal-Hydraulic Simulations

While the literature has numerous examples of Monte Carlo and computational fluid dynamics (CFD) coupling, most are hard-wired codes intended primarily for research rather than as standalone, general-purpose applications. In this work, we describe an open source application, ENRICO, that enables coupled neutronic and thermal-hydraulic simulations between multiple codes that can be chosen at runtime (as opposed to a coupling between two specific codes). The application has been designed such that the control flow logic, domain mapping, nonlinear fixed-point iteration, solution transfers, and convergence checks are all agnostic to the underlying physics solvers used. Special emphasis has also been placed on enabling efficient execution on distributed-memory computing environments. The transfer of solution fields between solvers is performed in memory rather than through filesystem I/O. Additionally, solvers can be configured to run on overlapping or disjoint sets of processes. To date, coupling with the OpenMC and Shift Monte Carlo codes, the Nek5000 CFD code, and a simplified heat diffusion and subchannel solver has been implemented in ENRICO. We present results for coupled simulations of a single light-water reactor fuel assembly based on the NuScale reactor using various combinations of the physics solvers. For this problem, the coupled simulations are shown to converge in about four Picard iterations. A comparison of the heat source and temperature distributions computed by ENRICO using OpenMC coupled with Nek5000 and Shift coupled with Nek5000 illustrates remarkable agreement between the codes.

42 ENGINEERING↗

Coupled-channel approach to isotensor π π π scattering from lattice QCD

The quest to understand three-body dynamics from first-principle QCD includes the study of nonresonant and resonant systems. The isospin I = 2 system is of particular interest having no three-body resonance but featuring a resonance in a subchannel, while also being a coupled-channel problem. In this study, we calculate the finite-volume spectrum from lattice QCD at two different pion masses, map the amplitude to the infinite volume through a generalized Finite-Volume Unitarity three-body quantization condition, investigate the limit of a narrow ρ , and compare with an effective Lagrangian prediction at leading order. Chiral extrapolations between different pion masses are performed.

Feng, Yuchuan [The George Washington University] (↗

Bound states in the B -matrix formalism for the three-body scattering

In this work, we consider a model of relativistic three-body scattering with a bound state in the two-body subchannel. We show that the naiïve K-matrix type parametrization, here referred to as the B-matrix, has nonphysical singularities near the physical region. We show how to eliminate such singularities by using dispersion relations and also show how to reproduce unitarity relations by taking into account all relevant open channels.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Analytic continuation of the relativistic three-particle scattering amplitudes

We investigate the relativistic scattering of three identical scalar bosons interacting via pair-wise interactions. Extending techniques from the nonrelativistic three-body scattering theory, we provide a detailed and general prescription for solving and analytically continuing integral equations describing the three-body reactions. We use these techniques to study a system with zero angular momenta described by a single scattering length leading to a bound state in a two-body subchannel. We obtain bound-state-particle and three-particle amplitudes in the previously unexplored kinematical regime; in particular, for real energies below elastic thresholds and complex energies in the physical and unphysical Riemann sheets. We extract positions of three-particle bound-states that agree with previous finite-volume studies, providing further evidence for the consistency of the relativistic finite-volume three-body quantization conditions. We also determine previously unobserved virtual bound states in this theory. Lastly, we find numerical evidence of the breakdown of the two-body finite-volume formalism in the vicinity of the left-hand cuts and argue for the generalization of the existing formalism.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

High-Fidelity Numerical Investigation on Elucidating Sodium Heat Transfer Characteristics for 37-Pin Wire-Wrapped Fuel Bundle in the PLANDTL Facility

This study involved a Reynolds-averaged Navier-Stokes- (RANS-) based computational fluid dynamics (CFD) analysis of the 37-pin wire-wrapped fuel bundle of the PNC Plant dynamics test loop (PLANDTL) facility. Previously, mainly the hydrodynamic phenomena of the wire-wrapped fuel bundle were analyzed, but the present study additionally included heat transfer analysis through conjugate heat transfer. The main purpose of the study was to benchmark the experimental data of the PLANDTL 37-pin wire-wrapped fuel bundle to investigate the heat transfer phenomena. In addition, the aim was to verify the accuracy of the RANS-based CFD analysis method using the STAR-CCM+ simulation software in comparison with the experimental data. The grid used for verification was an innovative grid system consisting of hexahedra using Fortran-based code. The development of the RANS-based CFD methodology included grid sensitivity analysis, turbulence model sensitivity analysis, and turbulent Prandtl number sensitivity analysis. Information on the temperature, mass flow rate, and area of the CFD results for each subchannel was provided for the top of the heated section and is expected to serve as a reference for future studies aiming to perform the validation and verification of a PLANDTL facility. In addition, the dependence of the peak temperature on the azimuth angle of each pin was analyzed.

97 MATHEMATICS AND COMPUTING↗

SubKit

Generating a CTF (thermal-hydraulics subchannel code) model can be time consuming and error prone because of the rigid formatting of the input file and the amount of input requirements. Because of the structured input format, making small changes to a model can require large scale changes to the input file. Furthermore, a user typically has to perform calculations to develop CTF input parameters from actual model geometry. Changing model geometry requires re-calculating input parameters and copying them from wherever the calculation was performed, be it a script, spreadsheet, or piece of paper. This process complicates the quality assurance process for that model and makes changing the model more difficult. CtfDeckBuilder is an object-oriented Python interface to CTF that allows the user to build a CTF model in Python. The user no longer needs to interface with the CTF input file when using the CtfDeckBuilder. The interface reduces the information required of the user by offering defaults for most modeling options provided by CTF. The user can perform calculations directly in Python when building the model so that calculations stay with the model itself.

Salko, Jr., RobertK↗

DUCTED ASSEMBLY STEADY STATE HEAT TRANSFER SOFTWARE

The Ducted Assembly Steady State Heat transfer software (DASSH) performs full core, steady state, subchannel temperature calculations for hexagonal, ducted assemblies. DASSH is targeted for use in the early stages of the reactor design process when assembly characteristics may not be defined or are undergoing considerable changes. DASSH provides a relatively rapid assessment of the temperature and flow distribution to inform reactor design updates and analysis.

ATZ, MILOSIVO↗

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↗

Development of a high-fidelity multi-physics coupling between MCNP6.2 and CTF4.0 for VVER applications

Ensuring system safety in the design, licensing, and operation phases is a priority in the nuclear industry. Performing extensive, full-scale reactor safety experiments is often prohibitive due to the large associated costs. Computational simulations offer an alternative safety analysis method, typically with significant cost reductions. Recent high-level developments in technology and increased availability of computational resources have allowed the development of high-fidelity, high-resolution multi-physics coupled codes. Such developments may be used to generate reference models for deterministic core calculations. Under this framework, the high-fidelity continuous energy Monte Carlo-based neutron transport code MCNP6.2 was coupled externally with the state-of-the-art thermal-hydraulics subchannel code, CTF4.0 for VVER (Water-Water Energetic Reactor) applications. A VVER-1000 fuel assembly model was used to demonstrate the capability of the coupled code. The converged coupled solution is compared to initial results, consisting of the first MCNP evaluation and first CTF evaluation after initialization. The VVER-1000 type assembly results are compared to other evaluations of the same assembly model. The findings indicate good agreement with expectations and reference cases, where available. The initial results of the coupled MCNP6.2/CTF4.0 calculations at the assembly level for steady-state calculations are presented in this study, which may support future work toward high-fidelity coupled full-core modeling capabilities. The multi-physics model may be further improved, enhanced, and expanded for both cycle depletion and transient applications. Such a multi-physics system will also be applicable to the VVER-1200 and other triangular lattice designs and support their deployment and operation safely and economically. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

CTF Void Drift Validation Study

This milestone report is a summary of work performed in support of expansion of the validation and verification (V&V) matrix for the thermal-hydraulic subchannel code, CTF. The focus of this study is on validating the void drift modeling capabilities of CTF and verifying the supporting models that impact the void drift phenomenon. CTF uses a simple turbulent-diffusion approximation to model lateral cross-flow due to turbulent mixing and void drift. The void drift component of the model is based on the Lahey and Moody model. The models are a function of two-phase mass, momentum, and energy distribution in the system; therefore, it is necessary to correctly model the ow distribution in rod bundle geometry as a first step to correctly calculating the void distribution due to void drift.

97 MATHEMATICS AND COMPUTING↗

STAR-CCM+ (CFD) Calculations and Validation L3:VVI.H2L.P15.02

This milestone presents a demonstration of the High-to-Low (Hi2Lo) process in the VVI focus area. Validation and additional calculations with the commercial computational fluid dynamics code, STAR-CCM+, were performed using a 5x5 fuel assembly with non-mixing geometry and spacer grids. This geometry was based on the benchmark experiment provided by Westinghouse. Results from the simulations were compared to existing experimental data and to the subchannel thermal-hydraulics code COBRA-TF (CTF). An uncertainty quantification (UQ) process was developed for the STAR-CCM+ model and results of the STAR UQ were communicated to CTF. Results from STAR-CCM+ simulations were used as experimental design points in CTF to calibrate the mixing parameter β and compared to results obtained using experimental data points. This demonstrated that CTF’s β parameter can be calibrated to match existing experimental data more closely. The Hi2Lo process for the STAR-CCM+/CTF code coupling was documented in this milestone and closely linked L3:VVI.H2LP15.01 milestone report.

97 MATHEMATICS AND COMPUTING↗

Report on UQ Assessments to support SESAME wire-wrapped bundle experiment

This work assesses the influence of assumptions made when generating a mesh of a wire-wrapped geometry. The contact region between a wire and its adjacent pin is commonly modeled by either embedding the wire to the adjacent pin or trimming the wire so that a gap separates the wire from its adjacent pin. These models are referred to as close-gap and open-gap approaches herein and are applied to two geometries. The first geometry consists of a single pin wire-wrapped subchannel. A polyhedral mesh and a hexahedral mesh are generated. The second and third geometry are a 7-pin and a 19-pinwire-wrapped bundles meshed with polyhedral elements only. Pressure drops are obtained with the STAR-CCM+computational fluid dynamic package. Sensitivity analyses of the mesh density, the mesh type, and the turbulent models are performed. Numerical results show that the best match to the experimental data and to the Cheng-Todreas correlation is obtained with the combination of a hexahedral mesh, the shear stress transport (SST) turbulent model, and the open-gap approach. In the case of the 7-pin geometry, the best results are obtained with the open-gap approach and the SST turbulent model. The 19-pin geometry yields contradictory results to the 7-pin geometry results, and thus will require further investigations.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

CASL VMA Milestone Report FY16 (L3:VMA.VUQ.P13.08): Westinghouse Mixing with STAR-CCM+

STAR-CCM+ (STAR) is a high-resolution computational fluid dynamics (CFD) code developed by CD-adapco. STAR includes validated physics models and a full suite of turbulence models including ones from the k-ε and k-ω families. STAR is currently being developed to be able to do two phase flows, but the current focus of the software is single phase flow. STAR can use imported meshes or use the built in meshing software to create computation domains for CFD. Since the solvers generally require a fine mesh for good computational results, the meshes used with STAR tend to number in the millions of cells, with that number growing with simulation and geometry complexity. The time required to model the flow of a full 5x5 Mixing Vane Grid Assembly (5x5MVG) in the current STAR configuration is on the order of hours, and can be very computationally expensive. COBRA-TF (CTF) is a low-resolution subchannel code that can be trained using high fidelity data from STAR. CTF does not have turbulence models and instead uses a turbulent mixing coefficient β. With a properly calibrated β, CTF can be used a low-computational cost alternative to expensive full CFD calculations performed with STAR. During the Hi2Lo work with CTF and STAR, STAR-CCM+ will be used to calibrate β and to provide high-resolution results that can be used in the place of and in addition to experimental results to reduce the uncertainty in the CTF results.

97 MATHEMATICS AND COMPUTING↗