Search NASA⌕ Search

SEARCH · Search NASA

Results for “lower plenum”

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 19 records

Investigation of thermal hydraulic behavior of the High Temperature Test Facility's lower plenum via large eddy simulation

A high-fidelity computational fluid dynamics (CFD) analysis was performed using the Large Eddy Simulation (LES) model for the lower plenum of the High–Temperature Test Facility (HTTF), a ¼ scale test facility of the modular high temperature gas-cooled reactor (MHTGR) managed by Oregon State University. In most next–generation nuclear reactors, thermal stress due to thermal striping is one of the risks to be curiously considered. This is also true for HTGRs, especially since the exhaust helium gas temperature is high. In order to evaluate these risks and performance, organizations in the United States led by the OECD NEA are conducting a thermal hydraulic code benchmark for HTGR, and the test facility used for this benchmark is HTTF. HTTF can perform experiments in both normal and accident situations and provide high-quality experimental data. However, it is difficult to provide sufficient data for benchmarking through experiments, and there is a problem with the reliability of CFD analysis results based on Reynolds–averaged Navier–Stokes to analyze thermal hydraulic behavior without verification. To solve this problem, high-fidelity 3-D CFD analysis was performed using the LES model for HTTF. It was also verified that the LES model can properly simulate this jet mixing phenomenon via a unit cell test that provides experimental information. As a result of CFD analysis, the lower the dependency of the sub-grid scale model, the closer to the actual analysis result. In the case of unit cell test CFD analysis and HTTF CFD analysis, the volume-averaged sub-grid scale model dependency was calculated to be 13.0% and 9.16%, respectively. As a result of HTTF analysis, quantitative data of the fluid inside the HTTF lower plenum was provided in this paper. As a result of qualitative analysis, the temperature was highest at the center of the lower plenum, while the temperature fluctuation was highest near the edge of the lower plenum wall. The power spectral density of temperature was analyzed via fast Fourier transform (FFT) for specific points on the center and side of the lower plenum. FFT results did not reveal specific frequency-dominant temperature fluctuations in the center part. It was confirmed that the temperature power spectral density (PSD) at the top increased from the center to the wake. The vortex was visualized using the well-known scalar Q-criterion, and as a result, the closer to the outlet duct, the greater the influence of the mainstream, so that the inflow jet vortex was dissipated and mixed at the top of the lower plenum. Additionally, FFT analysis was performed on the support structure near the corner of the lower plenum with large temperature fluctuations, and as a result, it was confirmed that the temperature fluctuation of the flow did not have a significant effect near the corner wall. In addition, the vortices generated from the lower plenum to the outlet duct were identified in this paper. It is considered that the quantitative and qualitative results presented in this paper will serve as reference data for the benchmark.

97 MATHEMATICS AND COMPUTING↗

Design of an experimental facility with a unit cell test section for studies of the lower plenum in prismatic high temperature gas reactors

A proposed design for a Generation IV reactor is the very-high temperature reactor (VHTR), which employs helium as the primary coolant. Although many advantages exist with this design, a number of challenges remain. One such area of concern is the turbulent mixing of non-isothermal flows in the lower plenum, where an array of coolant jets collectively mix and traverse around a series of structural support posts before exiting the core. Neighboring jets in the lower plenum could experience temperature differences as high as 300–400 K, causing potential hot streaking in the support posts. Additional problems could occur further downstream if the level of non-uniformity in the temperatures is too high (e.g, inducing thermal stresses on gas turbine blades). Therefore, it is important that adequate tools exist to predict the turbulent mixing within the VHTR lower plenum. This task remains difficult to achieve due to the complex geometry and flow physics that dictate the mixing. Various modeling approaches can only be assessed using data acquired from a carefully designed experimental facility. This paper presents details of such a facility, where high-fidelity measurements can be acquired for a scaled portion of the lower plenum defined herein as a ‘‘unit cell”, which consists of a hexagonal array of jets whose flow is directed into a test section with support posts positioned between neighboring jets, and subjected to a crossflow. Through adjusting velocity and temperature conditions at the jet inlets, different locations within the lower plenum can be experimentally simulated. System response quantities (SRQs) include both temperature and velocity measurements, captured via multiple thermocouples embedded in each post, and particle image velocimetry (PIV), respectively. These SRQs provide assessment metrics for both the momentum and energy models present in current computational solvers, and therefore the experimental facility is a valuable source for validation data. Another benefit is a more comprehensive understanding of the thermal-hydraulic loading conditions and flow physics present in the VHTR lower plenum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

CFD Simulations of Lower Plenum Mixing

Review of model development and validation performed in the Advanced Reactor Technologies (ART) program for thermal mixing at the outlet of High Temperature Gas Reactors (HTGRs). Understanding the mixing that occurs in the lower plenum in an HTGR is necessary to facilitate design improvements and to perform reactor safety analysis. Numerical models are one possible approach to gain a better understanding of mixing in the lower plenum. Given the complexity of the geometry and the intense mixing present, it is important to perform validation of numerical models. Three models have been developed during FY2025: a porous media with Pronghorn, a Reynolds Averaged Navier Stokes (RANS) with STAR-CCM+, and a Large Eddy Simulation (LES) with NekRS. The reference facility is a scaled-down version of the lower plenum of the High Temperature Gas-Cooled Reactor - Pebble-bed Module (HTR-PM) demonstration reactor. Preliminary results of the porous media and the RANS shows general good agreement against experimental benchmark data. Future work will leverage high-fidelity results obtained through LES to guide model selection and improvements to the lower-fidelity models, with particular attention to the Pronghorn porous media.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗

Computational study of full-scale VHTR lower plenum for turbulent mixing assessment

Next-generation nuclear reactors are poised to efficiently provide reliable power at enhanced safety levels for many years to come. Among the fundamental designs for these reactors is the very-hightemperature reactor (VHTR), which employs helium as the primary coolant and has the potential of reaching elevated temperatures such that chemical processing (e.g., electrolysis for hydrogen production) is viable alongside traditional electricity generation. This and other next-generation reactors also represent a strategic bridge for ultimately transitioning from fossil fuel-dominated energy portfolios to one where renewable options abound. For the VHTR, the ultimate goal of constructing a new plant will be preceded by additional fundamental research aimed at generating trusted models for behavior prediction under normal and accident scenarios. This paper presents a detailed computational fluid dynamics simulation of the lower plenum, where hot coolant from the core mixes together in a turbulent fashion before traveling to power conversion equipment. Because the flow is expected to enter the lower plenum across a wide range of temperatures and velocities, concerns exist when the mixing is incomplete. The potential hot spots on cylindrical support pedestals and the uniformity of the main outlet are the primary metrics of interest in this study, along with flow velocities and temperatures at different locations within the lower plenum. Three locations are probed in detail and suggest a large range exists for approach velocities and temperatures. A large degree of stratification is also seen on the surfaces of the support pedestals, suggesting a facility should accommodate testing for such behavior. The characterizations presented provide the valuable data needed in order to design appropriate experimental testbeds, where scaled modeling can be carried out in a manner meaningful in predicting the full-scale behavior.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Use of scaled BWR lower plenum boron mixing tests to qualify the boron transport model used in TRACG

In 2001 GEH applied best estimate methods combined with a statistical methodology to determine upper bound limits for key licensing parameters for anticipated operation occurrence (AOO) transient and anticipated transients without scram (ATWS) overpressure analyses for operating Boiling Water Reactors (BWRs). The methodology was subsequently extended for ESBWR AOO, ATWS, loss of coolant, and stability analyses. GEH is extending the methodology to long-term ATWS analyses for the operating BWRs. A long-term ATWS scenario uses injection of borated water to achieve reactor shutdown. Predicting the mixing and transport of boron is important for calculating the impact on the key licensing parameters. For the many operating BWRs where the denser boron solution is injected into the lower plenum, stratification may occur, delaying boron transport to the core region. CFD modeling can be used to model the stratification and mixing of the boron solution, but such calculations are extremely computer intensive and not cost effective; therefore, a more-empirical approach supported by a theoretical scaling of the dominant phenomena and backed by test data and benchmark calculations is used. The paper presents the TRACG lower plenum boron transport model qualification effort. The scaling basis used to implement the TRACG boron transport model for BWR applications is discussed. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

OSU High Temperature Test Facility Test Acceptance Report, PG-28 Low Power (<350kW) Lower Plenum Mixing Test

Test Acceptance Report for Test PG-28 at the Oregon State University High Temperature Test Facility. The purpose of this test is to investigate the impact on Lower Plenum thermal profiles and mixing due to changes in the system mass flow rate in an integral test facility scaled to General Atomics’ Modular High Temperature Gas Reactor (MHTGR) design.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

OSU High Temperature Test Facility Test Acceptance Report, PG-30 Low Power (<350kW) Lower Plenum Mixing, Constant Temperature Test

Test Acceptance Report for Test PG-30 at the Oregon State University High Temperature Test Facility. The purpose of this test is to investigate the impact on Lower Plenum thermal profiles and mixing due to changes in the system mass flow rate in an integral test facility scaled to General Atomics’ Modular High Temperature Gas Reactor (MHTGR) design.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Coupled computational fluid dynamics and heat transfer analysis of the VHTR lower plenum.

The very high temperature reactor (VHTR) concept is being developed by the US Department of Energy (DOE) and other groups around the world for the future generation of electricity at high thermal efficiency (> 48%) and co-generation of hydrogen and process heat. This Generation-IV reactor would operate at elevated exit temperatures of 1,000-1,273 K, and the fueled core would be cooled by forced convection helium gas. For the prismatic-core VHTR, which is the focus of this analysis, the velocity of the hot helium flow exiting the core into the lower plenum (LP) could be 35-70 m/s. The impingement of the resulting gas jets onto the adiabatic plate at the bottom of the LP could develop hot spots and thermal stratification and inadequate mixing of the gas exiting the vessel to the turbo-machinery for energy conversion. The complex flow field in the LP is further complicated by the presence of large cylindrical graphite posts that support the massive core and inner and outer graphite reflectors. Because there are approximately 276 channels in the VHTR core from which helium exits into the LP and a total of 155 support posts, the flow field in the LP includes cross flow, multiple jet flow interaction, flow stagnation zones, vortex interaction, vortex shedding, entrainment, large variation in Reynolds number (Re), recirculation, and mixing enhancement and suppression regions. For such a complex flow field, experimental results at operating conditions are not currently available. Instead, the objective of this paper is to numerically simulate the flow field in the LP of a prismatic core VHTR using the Sandia National Laboratories Fuego, which is a 3D, massively parallel generalized computational fluid dynamics (CFD) code with numerous turbulence and buoyancy models and simulation capabilities for complex gas flow fields, with and without thermal effects. The code predictions for simpler flow fields of single and swirling gas jets, with and without a cross flow, are validated using reported experimental data and theory. The key processes in the LP are identified using phenomena identification and ranking table (PIRT). It may be argued that a CFD code that accurately simulates simplified, single-effect flow fields with increasing complexity is likely to adequately model the complex flow field in the VHTR LP, subject to a future experimental validation. The PIRT process and spatial and temporal discretizations implemented in the present analysis using Fuego established confidence in the validation and verification (V and V) calculations and in the conclusions reached based on the simulation results. The performed calculations included the helicoid vortex swirl model, the dynamic Smagorinsky large eddy simulation (LES) turbulence model, participating media radiation (PMR), and 1D conjugate heat transfer (CHT). The full-scale, half-symmetry LP mesh used in the LP simulation included unstructured hexahedral elements and accounted for the graphite posts, the helium jets, the exterior walls, and the bottom plate with an adiabatic outer surface. Results indicated significant enhancements in heat transfer, flow mixing, and entrainment in the VHTR LP when using swirling inserts at the exit of the helium flow channels into the LP. The impact of using various swirl angles on the flow mixing and heat transfer in the LP is qualified, including the formation of the central recirculation zone (CRZ), and the effect of LP height. Results also showed that in addition to the enhanced mixing, the swirling inserts result in negligible additional pressure losses and are likely to eliminate the formation of hot spots.

08 HYDROGEN↗

An intelligent scanning algorithm for lost fuel rod recovery in the MYRRHA lower plenum

This paper presents a system for visualizing the interior structure of a vessel filled with a liquid metal. Due to the opaque and conductive nature of the medium, regular visualization systems based on cameras, electromagnetic or laser based techniques can't be used. The operation, strengths and weaknesses of an existing set-up are investigated by creating a scenario to find and extract a rectangular object inside the vessel in the shortest possible time. The main focus of this paper is not on designing a new system but to investigate the bottlenecks that are faced and need to be kept in mind in future designs. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Swirling jets for the mitigation of hot spots and thermal stratification in the VHTR lower plenum

The Department of Energy (DOE) has a special interest in the development of Generation IV next generation nuclear plants (NGNPs) in the form of very high (gas cooled) temperature reactors (VHTRs) for the production of electricity, hydrogen, and process heat. As of 2010, the prismatic and pebble bed reactors are the top VHTR candidates [ORNL, 2010; Wikipedia, 2010a; World Nuclear News, 2010; Kinsey, 2010]. However, for reasons elucidated in Section 2, this research will concentrate on the prismatic core VHTR.

Rodriguez, Salvador B.↗

Experimental and CFD Studies of Coolant Flow Mixing within Scaled Models of the Upper and Lower Plenums of NGNP Gas-Cooled Reactors

A 1/16th scaled VHTR experimental model was constructed and the preliminary test was performed in this study. To produce benchmark data for CFD validation in the future, the facility was first run at partial operation with five pipes being heated. PIV was performed to extract the vector velocity field for three adjacent naturally convective jets at statistically steady state. A small recirculation zone was found between the pipes, and the jets entered the merging zone at 3 cm from the pipe outlet but diverged as the flow approached the top of the test geometry. Turbulence analysis shows the turbulence intensity peaked at 41-45% as the jets mixed. A sensitivity analysis confirmed that 1000 frames were sufficient to measure statistically steady state. The results were then validated by extracting the flow rate from the PIV jet velocity profile, and comparing it with an analytic flow rate and ultrasonic flowmeter; all flow rates lie within the uncertainty of the other two methods for Tests 1 and 2. This test facility can be used for further analysis of naturally convective mixing, and eventually produce benchmark data for CFD validation for the VHTR during a PCC or DCC accident scenario. Next, a PTV study of 3000 images (1500 image pairs) were used to quantify the velocity field in the upper plenum. A sensitivity analysis confirmed that 1500 frames were sufficient to precisely estimate the flow. Subsequently, three (3, 9, and 15 cm) Y-lines from the pipe output were extracted to consider the output differences between 50 to 1500 frames. The average velocity field and standard deviation error that accrued in the three different tests were calculated to assess repeatability. The error was varied, from 1 to 14%, depending on Y-elevation. The error decreased as the flow moved farther from the output pipe. In addition, turbulent intensity was calculated and found to be high near the output. Reynolds stresses and turbulent intensity were used to validate the data by comparing it with benchmark data. The experimental data gave the same pattern as the benchmark data. A turbulent single buoyant jet study was performed for the case of LOFC in the upper plenum of scaled VHTR. Time-averaged profiles show that 3,000 frames of images were sufficient for the study up to second-order statistics. Self-similarity is an important feature of jets since the behavior of jets is independent of Reynolds number and a sole function of geometry. Self-similarity profiles were well observed in the axial velocity and velocity magnitude profile regardless of z/D where the radial velocity did not show any similarity pattern. The normal components of Reynolds stresses have self-similarity within the expected range. The study shows that large vortices were observed close to the dome wall, indicating that the geometry of the VHTR has a significant impact on its safety and performance. Near the dome surface, large vortices were shown to inhibit the flows, resulting in reduced axial jet velocity. The vortices that develop subsequently reduce the Reynolds stresses that develop and the impact on the integrity of the VHTR upper plenum surface. Multiple jets study, including two, three and five jets, were investigated.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Standard Problems for CFD Validation for NGNP - Status Report

The U.S. Department of Energy (DOE) is conducting research and development to support the resurgence of nuclear power in the United States for both electrical power generation and production of process heat required for industrial processes such as the manufacture of hydrogen for use as a fuel in automobiles. The project is called the Next Generation Nuclear Plant (NGNP) Project, which is based on a Generation IV reactor concept called the very high temperature reactor (VHTR). The VHTR will be of the prismatic or pebble bed type; the former is considered herein. The VHTR will use helium as the coolant at temperatures ranging from 250°C to perhaps 1000°C. While computational fluid dynamics (CFD) has not previously been used for the safety analysis of nuclear reactors in the United States, it is being considered for existing and future reactors. It is fully recognized that CFD simulation codes will have to be validated for flow physics reasonably close to actual fluid dynamic conditions expected in normal operational and accident situations. The “Standard Problem” is an experimental data set that represents an important physical phenomenon or phenomena, whose selection is based on a phenomena identification and ranking table (PIRT) for the reactor in question. It will be necessary to build a database that contains a number of standard problems for use to validate CFD and systems analysis codes for the many physical problems that will need to be analyzed. The first two standard problems that have been developed for CFD validation consider flow in the lower plenum of the VHTR and bypass flow in the prismatic core. Both involve scaled models built from quartz and designed to be installed in the INL’s matched index of refraction (MIR) test facility. The MIR facility employs mineral oil as the working fluid at a constant temperature. At this temperature, the index of refraction of the mineral oil is the same as that of the quartz. This provides an advantage to the optics used for data gathering. Particle image velocimetry (PIV) is used to take the data. The first standard problem represents several flow physics expected to be present in the lower plenum of the prismatic VHTR. In the lower plenum, heated helium coolant in the form of jets issues downward into the plenum and is then forced to turn ninety degrees and flow toward the exit duct. The lower plenum is filled with cylindrical graphite posts that hold up the core. Figure S-1 provides a plan view of the lower plenum. The red circles represent support posts holding up columns of heated blocks. Grey circles represent support posts under columns of reflector blocks. Helium enters the lower plenum at the junctions of the hexagonal blocks.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Experimentally Validated Numerical Models of Non-Isothermal Turbulent Mixing in High Temperature Reactors

In the Very High Temperature Reactor (VHTR), the non-uniform heat generation and flow maldistribution of coolant in the core generates conditions leading into the lower plenum with many jets at different temperatures and velocities. The hope is that the combined flow becomes fully mixed in a short distance so that the fluctuating temperatures in the mixing region do not cause unwanted thermal stresses in the supporting structures. These cyclic thermal stresses can cause degradation of the wall microstructure including crack initiation and growth, ultimately leading to the fatigue failure in the lower plenum, and can be accelerated by the high temperature environment in the VHTR. The focus of this project was to provide experimental data for future validation of computational efforts aimed at capturing the turbulent mixing in the lower plenum. Computational fluid dynamics (CFD) is also used as part of this effort to assess the requirements to model fundamental flows with application to the lower plenum as well at the applied flow conditions that mimic the experimental facilities. These efforts will lay the groundwork for future incorporation of stress analysis as well as fatigue and failure studies for multiple VHTR components. An additional emphasis is placed on reduced order models (ROM) using canonical flows, which will help to balance computational expense and accuracy and enable future work in uncertainty quantification. The flow in the lower plenum is extremely complex with hundreds of cooling jets entering and then being forced to change directions by 90° while navigating through an array of cylindrical support posts. Generating experimental data for such a complicated flow is difficult, not to mention that the computational requirements to generate an adequate temporal and spatial sampling is extremely expensive and time consuming. Therefore, the ultimate focus of this research is a “unit cell” which consists of six jets in the presence of a crossflow, and seven cylindrical support posts. It is scaled to achieve similar ranges of key dimensionless numbers. This serves to mimic flow conditions in different regions of the lower plenum by adjusting velocities, temperatures, and orientations, and represents a compromise between the complexity inside the full-scale application and the ideal conditions of a more fundamental study. In addition to the unit cell experiments and simulations, the complex flows in the VHTR lower plenum are broken down into a number of simplified (in some cases even fundamental) flows, two of which are turbulent jets and flow past cylinders. Each of these are studied, computationally and experimentally for isothermal and non-isothermal conditions. Flows where trusted experimental data is available in the literature (e.g., turbulent free jet), additional experiments are not needed. The result of these efforts is a database of experimental data and multiple computational studies, which provide an outlook to the requirements, needed to capture turbulent mixing flow physics The ROM results provide assessments for appropriate resource-light models that balance accuracy with reduced computational burdens.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Computational Fluid Dynamic Analysis of Core Bypass Flow Phenomena in a Prismatic VHTR

The core bypass flow in a prismatic very high temperature gas-cooled reactor (VHTR) is one of the important design considerations which impacts considerably on the integrity of reactor core internals including operating fuels. The interstitial gaps are an inherent presence in the reactor core because of tolerances in manufacturing the blocks and the inexact nature of their installation. Furthermore, the geometry of the graphite blocks changes over the lifetime of the reactor because of thermal expansion and irradiation damage. The occurrence of hot spots in the core and lower plenum and hot streaking in the lower plenum (regions of very hot gas flow) will be affected by the bypass flow. In the present study, three-dimensional computational fluid dynamic (CFD) calculations of a typical prismatic VHTR are conducted to understand better the bypass flow phenomenon and establish the evaluation method in the reactor core using commercial CFD code FLUENT. Parametric calculations changing several factors in a on-twelfth sector of a fuel column are performed. The simulations show the impact of each factor on bypass flow and the flow and temperature distributions in the prismatic core. The factors inlcude inter-column gap-width, turbulence model, axial heat generation profile and geometry change from irradiation-induced shrinkage in the graphite block region. It is shown that bypass flow provides a significant cooling effect on the prismatic block and that the maximum fuel and coolant channel outlet temperatures increase with an increase in gap-width, especially when a peak radial factor is applied to the total heat generation rate. Also, the presence of bypass flow causes a large lateral temperature gradient in the block that may have repurcussions on the structural integrity of the block and on the neutronics. These results indicate that bypass flow has a significant effect on hot spots in the core and on the temperature of jets flowing from the core into the lower plenum.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Material Interactions in Severe Accidents – Benchmarking the MELCOR V2.2 Eutectics Model for a BWR-3 MARK-I Station Blackout: Part I – Single Case Analysis

Here in this analysis, the two material interaction models available in the MELCOR code are benchmarked for a severe accident at a BWR under representative Fukushima Daiichi boundary conditions. This part of the benchmark investigates the impact of each material interaction model on accident progression through a detailed single case analysis. It is found that the eutectics model simulation exhibits more rapid accident progression for the duration of the accident. The slower accident progression exhibited by the interactive materials model simulation, however, allows for a greater degree of core material oxidation and hydrogen generation to occur, as well as elevated core temperatures during the ex-vessel accident phase. The eutectics model simulation exhibits more significant degradation of core components during the late in-vessel accident phase – more debris forms and relocates to the lower plenum before lower head failure. The larger debris bed observed in the eutectics model simulation also reaches higher temperatures, presenting a more significant thermal challenge to the lower head until its failure. At the end of the simulated accident scenario, however, core damage is comparable between both simulations due to significant core degradation that occurs during the ex-vessel phase in the interactive materials model simulation. A key difference between the two models’ performance is the maximum temperatures that can be reached in the core and therefore the maximum ΔT between any two components. When implementing the interactive materials model, users have the option to modify the liquefaction temperature of the ZrO 2 -interactive and UO 2 -interactive materials as a way to mimic early fuel rod failure due to material interactions. Through modification of the liquefaction of high melting point materials with significant mass, users may inadvertently limit maximum core temperatures for fuel, cladding, and debris components.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Experimental Investigation of Forced Convection and Natural Circulation Cooling of a VHTR Core under Normal Operation and Accident Scenarios

This NEUP project #15-8205 entitled “Experimental Investigation of Forced Convection and Natural Circulation Cooling of a VHTR Core under Normal Operation and Accident Scenarios”, has investigated important phenomena in the safety of Very High Temperature Reactors (VHTRs) using the experimental facilities at City College of New York (CCNY) and Kansas State University (KSU), in collaboration with the Idaho State University (ISU) and Idaho National Laboratory (INL). The main objectives of this project were to investigate the forced convection and bypass flow phenomena, natural circulation flow and heat transfer, and graphite oxidation due to air ingress which could occur in a Very High Temperature Reactor (VHTR) with a prismatic core. High pressure/high temperature experiments have been conducted to obtain data that could be used for validation of VHTR design and safety analysis codes. The focus of these experiments was on the generation of benchmark data for design and off-design heat transfer for forced and natural circulation in a VHTR core. This work has built on the previous NEUP project #11-3218 which investigated forced convection and flow laminarization phenomena in a prismatic core of a VHTR and natural circulation of a pure gas under Pressurized Conduction Cooldown (PCC) and Depressurized Conduction Cooldown (DCC) conditions. Under each Task, extensive literature reviews on natural circulation, bypass flow, air ingress and graphite oxidation phenomena were conducted. In Task 1, scaling analyses were conducted by Idaho State University followed by design and construction of experimental facilities needed to conduct both bypass flow and natural circulation experiments at CCNY. Natural circulation experiments were conducted in Task 2 to investigate the air ingress phenomena, and bypass flow phenomena in Task 3, both at CCNY. In Task 4, radiation heat transfer and graphite oxidation tests were performed at Kansas State University. The experimental data were then analyzed to obtain quantitative results on natural circulation of pure gases and gas mixtures under DCC conditions, bypass flow rates, radiation heat transfer and graphite oxidation. In addition, CFD and Direct Numerical Simulations were conducted by INL and CCNY to investigate the flow laminarization phenomenon in strongly heated gas flows. The scaling analyses have shown that the key nondimensional numbers and parameters all have acceptable distortion levels principally because the CCNY facility’s test sections have representative cooling channel dimensions and are using either helium or a scalable gas such as nitrogen at representative state conditions. The natural circulation experiments involving helium and nitrogen gas mixtures at CCNY have elucidated the effect of nitrogen (simulating air) ingress into a lower plenum on the onset of natural circulation, mass flow rates and heat transfer coefficients under different power input levels in the riser section. Also, the transport of nitrogen from the lower plenum to upper plenum due to natural circulation was quantified for the first time. Nusselt number correlations for the heat transfer coefficient under natural circulation as a function of Reynolds and Prandtl numbers was also developed. The experiments at KSU revealed that prior to the onset of natural circulationin a VHTR core. This work has built on the previous NEUP project #11-3218 which investigated forced convection and flow laminarization phenomena in a prismatic core of a VHTR and natural circulation of a pure gas under Pressurized Conduction Cooldown (PCC) and Depressurized Conduction Cooldown (DCC) conditions. Under each Task, extensive literature reviews on natural circulation, bypass flow, air ingress and graphite oxidation phenomena were conducted. In Task 1, scaling analyses were conducted by Idaho State University followed by design and construction of experimental facilities needed to conduct both bypass flow and natural circulation experiments at CCNY. Natural circulation experiments were conducted in Task 2 to investigate the air ingress phenomena, and bypass flow phenomena in Task 3, both at CCNY. In Task 4, radiation heat transfer and graphite oxidation tests were performed at Kansas State University. The experimental data were then analyzed to obtain quantitative results on natural circulation of pure gases and gas mixtures under DCC conditions, bypass flow rates, radiation heat transfer and graphite oxidation. In addition, CFD and Direct Numerical Simulations were conducted by INL and CCNY to investigate the flow laminarization phenomenon in strongly heated gas flows. The scaling analyses have shown that the key nondimensional numbers and parameters all have acceptable distortion levels principally because the CCNY facility’s test sections have representative cooling channel dimensions and are using either helium or a scalable gas such as nitrogen at representative state conditions.The natural circulation experiments involving helium and nitrogen gas mixtures at CCNY have elucidated the effect of nitrogen (simulating air) ingress into a lower plenum on the onset of natural circulation, mass flow rates and heat transfer coefficients under different power input levels in the riser section. Also, the transport of nitrogen from the lower plenum to upper plenum due to natural circulation was quantified for the first time. Nusselt number correlations for the heat transfer coefficient under natural circulation as a function of Reynolds and Prandtl numbers was also developed. The experiments at KSU revealed that prior to the onset of natural circulation(ONC), molecular diffusion plays a significant role in air-ingress, however, at higher temperatures convection currents may be influencing air ingress and in-turn ONC times. Bypass flow experiments faced difficulties in the measurement of gas velocity in the main channel due to the failure of a high temperature Hot Wire Anemometer probe, so an energy balance method was instead employed to estimate the flow distribution between the main and bypass flow channels. In the graphite oxidation experiments conducted at KSU, oxidation was found to increase the thermal diffusivity of nuclear grade graphite. This change is believed to be driven by oxidation decreasing graphite sample density rather than raising the thermal conductivity. Oxidation was also shown to greatly alter graphite emissivity. If large amounts of impurities have built up within the graphite over time and have significantly increased its effective ash content, then air ingress is likely to severely inhibit the radiative heat transfer abilities of the graphite due to the formation of low emissivity, ash layers on its surface. However, if the purity of the graphite at the time of the air ingress is still close to that of its virgin specifications, the oxidation induced roughening of the graphite’s outer surface would have the effect of slightly raising its emissivity. The flow laminarization phenomenon was investigated by INL by examining the results of Direct Numerical Simulations of previous experiments, and by CCNY conducting its own DNS simulations using DOE’s Nek5000 code

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗