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 55 records · Page 3

Intracore Natural Circulation Study in the High Temperature Test Facility

The development of the Modular High-Temperature Gas-Cooled Reactor is a significant milestone in advanced nuclear reactor technology. One of the concerns for the reactor’s safe operation is the effects of a loss-of-flow accident (LOFA) where the coolant circulators are tripped, and forced coolant flow through the core is lost. Depending on the steam generator placement, loop or intracore natural circulation develops to help transfer heat from the core to the reactor cavity, cooling system. This paper investigates the fundamental physical phenomena associated with intracore coolant natural circulation flow in a one-sixth Computational Fluid Dynamics (CFD) model of the Oregon State University High Temperature Test Facility (OSU HTTF) following a loss-of-flow accident transient. This study employs conjugate heat transfer and steady-state flow along with an SST k-ω turbulence model to characterize the phenomenon of core channel-to-channel natural convection. Previous studies have revealed the importance of complex flow distribution in the inlet and outlet plenums with the potential to generate hot coolant jets. For this reason, complete upper and lower plenum volumes are included in the analyzed computational domain. CFD results also include parametric studies performed for a mesh sensitivity analysis, generated using the STAR-CCM+ software. The resulting channel axial velocities and flow directions support the test facility scaling analysis and similarity group distortions calculation.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Preliminary Investigation of Rapid Depressurization Phenomena Following a Sudden DLOFC in a VHTR

Air ingress has been identified as a potential threat for Very High Temperature gas-cooled Reactors (VHTR). Reactor components constructed of graphite will, at high temperatures, produce exothermic reactions in the presence of oxygen. The danger lies in the possibility of fuel element damage and core structural failure. Previous investigations of air ingress mechanisms have focused on thermal and molecular diffusion, density-driven stratified flow due to hydrodynamic instability, and natural convection. Here, we investigate the possibility of a rapid flow reversal of helium coolant due to a Taylor (rarefaction) wave expansion after a hypothetical sudden Depressurized Loss of Forced Cooling (DLOFC) scenario in a VHTR. Conceivably, flow reversal of the helium coolant could entrain significant quantities of air into the reactor vessel. Our analysis starts with a one-dimensional shock tube simulation to simply illustrate the development of a Taylor wave with resulting reentrant flow. Then, a simulation is performed of an idealized two-dimensional axisymmetric representation of the lower plenum of General Atomics GT-MHR subjected to a hypothetical catastrophic break of the hot duct. Analysis shows the potential for significant and rapid air ingress into the reactor vessel in the case of a large break in the cooling system.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

AIR INGRESS ANALYSIS: PART 2 – COMPUTATIONAL FLUID DYNAMIC MODELS

The Idaho National Laboratory (INL), under the auspices of the U.S. Department of Energy, is performing research and development that focuses on key phenomena important during potential scenarios that may occur in very high temperature reactors (VHTRs). Phenomena Identification and Ranking Studies to date have ranked an air ingress event, following on the heels of a VHTR depressurization, as important with regard to core safety. Consequently, the development of advanced air ingress-related models and verification and validation data are a very high priority. Following a loss of coolant and system depressurization incident, air will enter the core of the High Temperature Gas Cooled Reactor through the break, possibly causing oxidation of the in-the core and reflector graphite structure. Simple core and plant models indicate that, under certain circumstances, the oxidation may proceed at an elevated rate with additional heat generated from the oxidation reaction itself. Under postulated conditions of fluid flow and temperature, excessive degradation of the lower plenum graphite can lead to a loss of structural support. Excessive oxidation of core graphite can also lead to the release of fission products into the confinement, which could be detrimental to a reactor safety. Computational fluid dynamic model developed in this study will improve our understanding of this phenomenon. This paper presents two-dimensional and three-dimensional CFD results for the quantitative assessment of the air ingress phenomena. A portion of results of the density-driven stratified flow in the inlet pipe will be compared with results of the experimental results.

99 GENERAL AND MISCELLANEOUS↗

Natural Circulation Patterns in the VHTR Air-Ingress Accident and Related Issues

Natural circulation patterns in the VHTR during a hypothetical air-ingress accident have been investigated using computational fluid dynamic (CFD) methods in order to compare results from the previous 1-D model which was developed using GAMMA code for the air-ingress analyses. The GT-MHR 600 MWt reactor was selected to be the reference design and modeled by a half symmetric 3-D geometry using FLUENT 6.3, a commercial CFD code. CFD simulations were carried out as the steady-state calculation, and the boundary conditions were either assumed or provided from the 1-D GAMMA code results. Totally, 12 different cases have been reviewed, and many notable results have been obtained through in this work. According to the simulations, natural circulation patterns in the reactor were quite different from the previous 1-D assumptions. A large re-circulation flow with thermal stratification phenomena was clearly observed in the hot-leg and the lower plenum in the 3-D model. This re-circulation flow provided about an order faster air-ingress speed (0.46 m/s in superficial velocity) than previously predicted by 1-D modeling (0.02~0.03 m/s). It indicates that the 1-D air-ingress modeling may significantly distort the air-ingress scenario and consequences. In addition, complicated natural circulation patterns are eventually expected to result in very complex graphite oxidations and corrosion behaviors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Mixed Convection in the VHTR in the Event of a LOFA

The US Department of Energy, Office of Nuclear Energy (DOE-NE) is supporting the development of a very high temperature reactor (VHTR) concept as the primary focus of it next generation nuclear power plant (NGNP) program. The VHTR is cooled by forcing helium downwards through the core into the lower plenum and out the hot duct. In the event that the coolant circulators are lost, the driving pressure drop across the core will reduce to zero and there will be the opportunity for natural circulation to occur. During the time that the circulators are powering down, the heat transfer in the core from the graphite blocks to the helium coolant will transform from turbulent forced convection to mixed convection, where buoyancy effects become important, to free or natural convection, where buoyancy is dominant. Analysis of the nature of the forced, mixed and free convection is best done using computational fluid dynamic (CFD) software that can provide fine details of the flow and heat transfer. However, CFD analysis involves approximations in the results because of the finite nature of the spatial and temporal discretizations required, the inexact nature of the turbulence models that are used and the finite precision of the computers employed. Therefore, it is necessary to validate the CFD computations. Validation is accomplished by comparing results from specific CFD computations to experimental data that have been taken specifically for the purpose of validation and that are related to the physical phenomena in question. The present report examines the flow and heat transfer parameters (dimensionless numbers) that characterize the flow and reports ranges for their values based on specific CFD studies performed for the VHTR.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Parametric Model of Debris Penetration Through Sump Strainers with Concurrent Filtration and Shedding

Two primary debris fates are of key interest to the risk-informed resolution analysis: (1) accumulation on the primary sump strainers with corresponding head loss, and (2) partial penetration of debris through the strainers with subsequent accumulation on fuel channel screens in the lower plenum. Accumulation on the core can impede flow, potentially causing local excessive heating and potentially limiting mixing to the point that boron precipitation can occur. Two opportunities exist for diverting debris away from the core: (1) a fraction of the ECCS recirculation flow feeds containment spray systems directly from the common header, (2) a fraction of debris-laden flow can bypass the core and proceed directly to the break without accumulating debris on the fuel assembly screen. Both of these pathways extend the time necessary to build debris levels of concern on either the strainers or the fuel assemblies.

42 ENGINEERING↗

Modular Accident Analysis Program (MAAP) - MELCOR Crosswalk: Phase II Analyzing a Partially Recovered Accident Scenario

Following the conclusion of the first phase of the crosswalk analysis, one of the key unanswered questions was whether or not the deviations found would persist during a partially recovered accident scenario, similar to the one that occurred in TMI - 2. In particular this analysis aims to compare the impact of core degradation morphology on quenching models inherent within the two codes and the coolability of debris during partially recovered accidents. A primary motivation for this study is the development of insights into how uncertainties in core damage progression models impact the ability to assess the potential for recovery of a degraded core. These quench and core recovery models are of the most interest when there is a significant amount of core damage, but intact and degraded fuel still remain in the cor e region or the lower plenum. Accordingly this analysis presents a spectrum of partially recovered accident scenarios by varying both water injection timing and rate to highlight the impact of core degradation phenomena on recovered accident scenarios. This analysis uses the newly released MELCOR 2.2 rev. 966 5 and MAAP5, Version 5.04. These code versions, which incorporate a significant number of modifications that have been driven by analyses and forensic evidence obtained from the Fukushima - Daiichi reactor site.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Validation Data Acquisition in HTTF during PCC Events

A validation experimental campaign was conducted in an Integral Effect Test (IET) facility of High Temperature Gas Reactors (HTGR), the High-Temperature Test Facility (HTTF) at Oregon State University (OSU). The HTTF simulates Depressurized and Pressurized Conduction Cooldowns (DCC and PCC). This campaign required the development of a new laser spectroscopic diagnostic to measure velocity in the challenging conditions of the HTTF: low speed (~1 m/s) gas flows at HTGR prototypical temperature and 1/10th pressure. This was a collaborative effort between co-PIs at The George Washington University (GW), Oregon State University (OSU), and NASA Langley Research Center. The main accomplishments of this project include the record for dynamic range for velocimetry, highest accuracy obtained with this technique, successful deployment in an IET leading to new validation matrix for CFD. These are detailed below and in manuscript appended to this executive summary. For this project, we introduced a new class of laser spectroscopic diagnostics to Thermal-Hydraulics to measure velocity of gases; furthermore, the diagnostic was demonstrated in-situ in an IET during DCC events. In such scenarios, particles used in mainstream techniques, like Particle Image Velocimetry (PIV) are not appropriate as they settle down too rapidly and also contaminate the experimental facility. Molecular tracers stay mixed within the working gas and can seed the flow in a technique called Molecular Tagging Velocimetry (MTV). In MTV a molecular tracer is photo-dissociated by a first (write) laser into radicals or molecules. The pattern created by the write laser is then interrogated with planar laser-induced fluorescence (PLIF), the read pulse(s), which are recorded with a camera. The pattern is probed and matched at two times (interval or probe time, dt), resulting in a time-of-flight velocimetry technique. This project demonstrated a new application range for MTV in gases. MTV has been extensively used for high-speed gas dynamics: molecules do not suffer from inertia associated with particles in regions of high accelerations (such as shock waves). For high-speed flows, the probe time is on the order of 0.1-10 µs. In time-of-flight techniques, the velocity is reconstructed from the recorded displacement, dx, between time interval dt: U = dx/dt. To first order, the resolution of the pattern matching is typically a (fixed) fraction of a pixel; therefore the resolution of the velocity measurement increases with dx and, for a given flow velocity, dt. For the low-speed flows of interest here, this requires long probe time, and therefore finding tracers that are stable for a long time. After an extensive literature survey we identified 6 tracers that would be suitable. We also worked with Spectra Physics to design a custom laser system that is flexible enough to probe these various tracers. The final instrument is made of 4 lasers that are integrated into a compact and movable cart to facilitate the deployment of the instrument off-site. One unique aspect of the laser system is the tunable dye laser (needed to do PLIF of the radicals). By being pumped by two Nd:YAG lasers, this laser is able to provide two pulses with any probe time (dt) for wavelengths between 200 and 900 nm with a linewidth of 3 pm. Accommodating the dual-pulse capability required customizing the pump lasers and the dye laser, making this system unique. Additionally, a custom UV-intensified CCD camera capable of frame-straddling and exposure time as short as 10 ns was acquired. More than two-third of the instrument cost (~$200k) was covered by a cost sharing from the PI’s startup funds. Of the 6 potential tracers identified, two were tested extensively (H2O and N2O) at high-temperature and elevated pressure at GW prior to conducting the tests at OSU. For both tracers, an Ar-F laser (193 nm, 10 ns, 10 mJ/pulse) photo-dissociates trace amounts (0.1 – 5% molar fraction) of the gases into OH and NO. OH is probed at 281.905 nm and NO at 226.186 nm (4-30mJ, 10 ns pulse). NO is very stable in inert environment and we could reconstruct velocities precisely with probe delay times as long at 40 ms, compared to 2 ms for OH, resulting in increased precision. Therefore, for the tests in the HTTF, N2O was selected as the final tracer and its optimal concentration was determined to be 0.5% molar. In the GW laboratory, we obtained a velocity measurement precision of 4 mm/s at conditions mimicking the HTTF tests. We also demonstrated the benefit of molecular tracers over particles and the flexibility of our instrument by continuously measuring velocity from 1000 to 0.1 m/s in a blow down test. This represents a dynamic range on velocity of 80 dB, which is unprecedented for velocimetry! Once the diagnostic was optimized, two two-week campaigns were conducted in the HTTF. This experimental campaign required numerous practical solutions to successfully and rapidly deploy this very sensitive table-top diagnostic in-situ. This was accomplished thanks to very extensive planning from all parties involved in the project. Due to constraints on availability of the HTTF, only adiabatic DCC tests were conducted. Specifically, the velocity was measured at the exit of the hot leg into the Reactor Cavity Simulation Tank (RCST). We were able to record the gas velocity associated with the air ingress for long times (up to one hour). Depending on the initial conditions (density of the gases in the Reactor Pressure Vessel and the RCST) the lock exchange lasted 30-60 s past the valve openings. Large vortices were identified within the shear layer at the interface of the cross flow. These vortices are an indication of shear-driven instability between the counter streams of different density. After the lock-exchange, a small residual flow was still present and its strength decayed exponentially. It finds its source in a buoyancy difference between the gases in both vessels due to the mixing and diffusion of the gases in the cross-over duct and the RPV lower plenum. We did not expect to capture this part of the transient, which would have been less significant (if not negligible) in the presence of a heated core; the precision of our measurement in the HTTF was 6 mm/s, which exceeded our expectations for such a complex test! Using commercial software (Fluent) we simulated the experiment with Computational Fluid Dynamic (CFD) using Reynolds Averaged Navier-Stokes (RANS) turbulence models. We were not able to accurately reproduce the shear layer instability, mixing, and diffusion captured in our experimental data, even with large 3D meshes (>10 M points). We attributed this to the inadequacy of the RANS models for the flow of interest: our experimental Reynolds number was too low and was outside the applicability range of the RANS models. Instead, we should have employed CFD with less modeling, such as Large Eddy (LES) or even Direct Numerical (DNS) simulations. Unfortunately, we did not have the computational resources to do so. Finally, it should be noted that this experimental technique holds great promises for measuring velocity in thermalhydraulics. MTV has an application domain (highlighted in the narrative) that is complimentary to more established (and simpler) techniques, such as PIV. It has attracted the attention of several thermalhydraulics researchers and sponsors to measure velocity in harsh environments. The spectroscopic nature of this technique also enables to expand it to measure temperature and concentration simultaneously to velocity. In lieu of a narrative, 4 detailed journal publications are appended to this executive summary to describe the various development and accomplishments during this project. They include 2 published manuscripts (Experiments in Fluids and Measurement Science and Technology), and two invited contributions following the NURETH-17 conference (Nuclear Science and Technology). Our current Nuclear Engineering and Design draft is not included at this time.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Fluid Stratification Separate Effects Analysis, Testing and Benchmarking

High Temperature Gas Reactors (HTGRs) are positioned to disrupt local and global markets via their unique ability to produce carbon-free process heat, high efficiency power generation, and passively safe operational features. However, significant impediments still exist to delay deployment of this particular technology, including a lack of experimental data, verified code application, and lack of consensus with regards to severe accident progression. In particular, air ingress accidents represent a particular challenge to designers and engineers, as they represent low probability, but highly complex, accident scenarios. Including phenomena such as molecular diffusion, free convection, and complex heat and mass transfer paths, experimental and traceable data is essential to maturing the state of the industry. Therefore, this work presents an experimental investigation of the transition to natural convection in HTGR applications using the Stratified Flow Separate Effects Test Facility, housed at Oregon State University. In particular, this work will present data that challenges the assumption that molecular diffusion is a significant factor in this severe accident in the reference facility of the General Atomic 600 MW th Gas Turbine-Modular Helium Reactor (GT-MHR). Rather, experimentally produced data shows a statistically suggestive retardant effect of cross duct orientation, indicating that ingress mechanics may be fundamentally altered via facility geometry. Experimentally, this is achieved using a simplified cross duct that may be positioned in one of two ways so as to provide either horizontal or vertical access to the lower plenum area. Onset of natural convection (ONC) is measured using an oxygen sensor probe, immersed in the helium working fluid, so as to provide direct indication of air presence in the upper plenum. This document will present the doctoral dissertation produced by this work, along with its analysis, experimental data, observations, and conclusions, in addition to select documents that will highlight critical changes during the development process. Those documents will include the quality assurance plan, all subsequent change requests, and the engineering transfer from Harris Thermal in order to establish the most complete quality record possible.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Finite Volume Discretization of the Euler Equations in Pronghorn

Modeling flow and heat transfer in high temperature gas reactors (HTGR) requires the ability to model a wide range of flow speeds from slow (natural convection), to intermediate (forced-flow conditions), to supersonic regimes (depressurization) for a wide range of geometries including the pebble bed, upper and lower plenum, and risers. In previous work, Pronghorn has effectively modeled low-to-medium speed flows in scenarios such as the one described in the two-dimensional PBMR-400 benchmark, using its finite-element-based streamline-upwind Petrov-Galerkin (SUPG) stabilized implementation of the Euler equations. However, limitations of this method become apparent when dealing with more complicated geometries (e.g. imposing slip boundary conditions at nodes belonging to two different boundaries) and when gas speeds are fast enough for shocks and supersonic flow to occur. For these problems, the finite-element-based solver lacks robustness and is plagued by slow iterative convergence or even divergence. In order to address these challenges, the Pronghorn code at INL has been updated with new, modified versions of its original equations. The new Pronghorn models are built on the finite volume method with a Harten-Lax-van Leer-Contact (HLLC) Riemann solver based numerical flux method, which (1) allows imposing slip boundary conditions much more robustly and (2) performs well for a wide range of flow speeds. The finite-volume-based flow solver will form the basis for a robust coarse-mesh thermal-hydraulics capability in Pronghorn.

97 MATHEMATICS AND COMPUTING↗

HTTF Benchmarking Activities in FY 2023

This report describes the work conducted in FY23 within the High Temperature Gas-Cooled Reactor (HTGR) Applications Drivers activity under the Multiphysics Applications technical area in the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program within DOE-NE. The focus of this past year’s activity was on the OECD/NEA thermal hydraulic benchmark for high temperature gas-cooled reactors using HTTF data (HTGR T/H). The activities focus on code-to-code comparisons for two benchmark problems: (1) Depressurized Conduction Cooldown using SAM, and (2) Lower Plenum Mixing using Nek5000/RS.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

High-Temperature Gas-Cooled Pebble-Bed Reactors Running In And Transient Modeling Capabilities Demonstration

This study presents a comprehensive benchmarking and verification effort of several thermal-hydraulic and multiphysics capabilities for high-temperature gas-cooled reactor (HTGR) applications. The first part of this effort focuses on the running-in verification of Griffin's multiphysics capabilities, specifically for simulating the evolution of Pebble Bed reactor cores from startup to equilibrium. In the absence of validation data, code-to-code comparisons are conducted with Kugelpy, showing good agreement for key quantities like maximum power density and fresh core k-eigenvalue predictions. However, discrepancies in equilibrium core predictions suggest potential issues with cross sections, underscoring the need for further refinement and evaluation. The HTTF system analysis code benchmark involves RELAP5-3D, SAM, and GAMMA+ to assess their predictive capabilities for HTTF behavior under both normal operation and pressurized conduction cooldown (PCC) transient conditions. While there is good agreement in predicting major parameters such as coolant temperature, solid temperature, and flow distribution, discrepancies in transient behavior highlight differences in modeling approaches, nodalizations, and heat transfer models. The HTTF lower plenum CFD benchmark employs nekRS to simulate flow mixing phenomena, successfully capturing relevant flow physics and demonstrating mesh independence in complex geometries. Preliminary results suggest a relatively uniform temperature field but significant unsteadiness in the flow, requiring time-averaging analyses. The GPBR200 system analysis code benchmark uses SAM's core channel and porous media models, incorporating an RCCS loop for decay heat removal. During steady-state and transient conditions, including protected de-pressurized and pressurized loss of forced cooling (DLOFC and PLOFC), both models show good agreement in predicting temperature profiles and key parameters. Notably, while the core channel model underpredicts convective heat transfer effects, both models maintain temperatures well below the TRISO fuel safety limit. These benchmarking efforts collectively enhance the predictive capabilities of the tools used in HTGR design and safety analysis, guiding developments to improve their accuracy and applicability.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

High-Temperature Gas-Cooled Reactors Multiphysics Simulation Demonstration and Code Validation

This study presents a comprehensive benchmarking and verification effort of several thermal-hydraulic and multiphysics capabilities for high-temperature gas-cooled reactor applications. The first part of this effort focuses on the running-in verification of Griffin’s multiphysics capabilities, specifically for simulating the evolution of pebble-bed reactor cores from startup to equilibrium. Since Fiscal Year 2024, improvements and enhancements have been implemented in Griffin, including simplifying the process to specify streamlines and developing the online cross-section generation capability. In the absence of validation data, code-to-code comparisons are conducted with kugelpy, showing good agreement for integral quantities like k-eff predictions and predictions for maximum power density. However, accuracy issues are noted for more detailed quantities like the spatial distribution of fission rate densities which will require further work to address. The second part of this report presents an improved System Analysis Module (SAM) core channel model where the effects of cross flow are considered during the pressurized loss of forced cooling transient, resulting in an improved agreement of the predicted pebble temperature with respect to the predictions from the SAM 2D porous media model. Additionally, the wall channeling effect due to variable porosity at the near wall region of the core is also investigated. Furthermore, to demonstrate Griffin’s online cross-section generation capability, a Multiphysics simulation is performed by coupling Griffin to the SAM core channel model. In the third part of the report, as a part of the Organisation for Economic Co-operation and Development/Nuclear Energy Agency (OECD/NEA) thermal-hydraulic code validation benchmark activity for a high-temperature gas-cooled reactor, the High Temperature Test Facility (HTTF) is investigated first using the NekRS computational fluid dynamics (CFD) code to study the flow mixing phenomenon in the lower plenum of the facility. Then, code-to-code and code-to-data comparisons are performed for Test PG27, which is a pressurized conduction cooldown (PCC) test, using five different codes by six organizations from five countries. The different simulations show good agreements in terms of the general trend but there are differences in some results such as the peak temperatures of different regions and heat removal rate.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Challenge Problem 1: Preliminary Results of the Direct Numerical Simulation of Transient Flows

This report presents the first direct numerical simulations (DNS) of transient mixed convection in an idealized downcomer-like channel (Challenge Problem 1, Phase II). Using the GPU-accelerated NekRS solver, we modeled a sudden decay in driving pressure, mimicking loss-of-flow events, and tracked the resulting evolution of Reynolds number, boundary-layer structure, turbulence statistics, and heat-transfer metrics. Key findings include the systematic thickening and eventual asymmetry of velocity and thermal boundary layers under buoyant deceleration; minimal “memory” lag in Reynolds shear stress and TKE profiles when sampled at matching Re, yet clear shifts of peak locations toward the cooled wall; overshoots in transient eddy-viscosity and eddy-diffusivity (and corresponding sub-unity turbulent Prandtl numbers) on the cooled side; and a pronounced transient Nusselt-number enhancement driven by wall-temperature inertia and residual eddy mixing. These effects combined to offer a temporary cooling margin above steady-state predictions during reactor LOF transients. Future work will extend this work to a more complex “Case II” geometry (90° turn + lower plenum) and generate multi-Re/Pr datasets for data-driven turbulence closures.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Air Ingress Analysis: Part 1 - Theoretical Approach

Idaho National Laboratory (INL), under the auspices of the U.S. Department of Energy (DOE), is performing research and development that focuses on key phenomena important during potential scenarios that may occur in very high temperature gas-cooled reactors (VHTRs). Phenomena identification and ranking studies to date have ranked an air ingress event, following on the heels of a VHTR depressurization, as important with regard to core safety. Consequently, the development of advanced air-ingress-related models and verification and validation data are a very high priority. Following a loss of coolant and system depressurization incident, air will enter the core of the VHTR through the break, possibly causing oxidation of the graphite core and reflector graphite structure. Simple core and plant models indicate that, under certain circumstances, the oxidation may proceed at an elevated rate with additional heat generated from the oxidation reaction itself. Under postulated conditions of fluid flow and temperature, excessive degradation of lower plenum graphite caused by graphite oxidation can lead to a loss of mechanical strength. Excessive oxidation of core graphite can also lead to a release of fission products into the confinement, which could be detrimental to reactor safety. Analytical models developed in this study will improve our understanding of this phenomenon. This paper presents two sets of analytical models for the qualitative assessment of the air ingress phenomena. The results from the analytical models are compared with results of the computational fluid dynamic models (CFD) in the subsequent paper. The analytical models agree well with those CFD results.

99 GENERAL AND MISCELLANEOUS↗

Natural Circulation Patterns in the VHTR Air-Ingress Accident and Related Issues

A natural circulation pattern in a Very High Gas-Cooled Reactor during a hypothetical air-ingress accident has been investigated using computational fluid dynamic (CFD) methods in order to compare with the previous 1-D flow path model for the air-ingress analyses. The GT-MHR 600 MWt reactor was selected to be the reference design and modeled by a half symmetric 3-D geometry using FLUENT 6.3, a commercial CFD code. The simulation was carried out as steady-state calculations, and the boundary conditions were either assumed or provided from the 1-D GAMMA code results. Totally, 12 different cases have been estimated, and many notable findings and results have been obtained in this study. According to the simulations, the natural circulation pattern in the reactor was quite different from the previous 1-D assumptions. A large re-circulation flow with thermal stratification phenomena was clearly observed in the hot-leg and the lower plenum in the 3-D model. This re-circulation flow provided approximately an order faster air-ingress speed (0.46 m/s in superficial velocity) than previously predicted values by 1-D modeling (0.02~0.03 m/s). It indicates that the 1-D air-ingress modeling may significantly distort the air-ingress scenario and consequences. In addition, the complicated natural circulation pattern is eventually expected to lead to very complex graphite oxidations and corrosion patterns.

99 GENERAL AND MISCELLANEOUS↗

Evaluation of in-vessel corium retention through external reactor vessel cooling for integral reactor

In-vessel corium retention through external reactor vessel cooling (IVR-ERVC) for a small integral reactor has been evaluated to determine the thermal margin for the prevention of a reactor vessel failure. A thermal load analysis from the corium pool to the outer reactor vessel wall in the lower plenum of the reactor vessel has been performed to determine the heat flux distribution. The critical heat flux (CHF) on the outer reactor vessel wall has been determined to fix the maximum heat removal rate through the external coolant between the outer reactor vessel and the insulation of the reactor vessel. Finally, the thermal margin has been evaluated by comparison of the thermal load with the maximum heat removal rate of the CHF on the outer reactor vessel wall. The maximum heat flux from the corium pool to the outer reactor vessel is estimated at approximately 0.25 MW/m{sup 2} in the metallic layer because of the focusing effect. The CHF of the outer reactor vessel is approximately 1.1 MW/m{sup 2} because of a two phase natural circulation mass flow. Since the thermal margin for the IVR-ERVC is sufficient, the reactor vessel integrity is maintained during a severe accident of a small integral reactor. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Development of a coolant channel helium and nitrogen gas ratio sensor for a high temperature gas reactor

To measure the changing gas composition of the coolant during a postulated High Temperature Gas Reactor (HTGR) accident, an instrument is needed. This instrument must be compact enough to measure the ratio of the coolant versus the break gas in an individual coolant channel. This instrument must minimally impact the fluid flow and provide for non-direct signal routing to allow minimal disturbance to adjacent channels. The instrument must have a flexible geometry to allow for the measurement of larger volumes such as in the upper or lower plenum of a HTGR. The instrument must be capable of accurately functioning through the full operating temperature and pressure of a HTGR. This instrument is not commercially available, but a literature survey has shown that building off of the present work on Capacitance Sensors and Cross-Capacitors will provide a basis for the development of the desired instrument. One difficulty in developing and instrument to operate at HTGR temperatures is acquiring an electrical conductor that will not melt at 1600 deg. C. This requirement limits the material selection to high temperature ceramics, graphite, and exotic metals. An additional concern for the instrument is properly accounting for the thermal expansion of both the sensing components and the gas being measured. This work covers the basic instrument overview with a thorough discussion of the associated uncertainty in making these measurements. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗