Engineering PapersSearch

SEARCH · Engineering Papers

Results for “REACTOR”

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 91 records · Page 5

Analysis of the Gas Core Actinide Transmutation Reactor (GCATR)

Design power plant studies were carried out for two applications of the plasma core reactor: (1) As a breeder reactor, (2) As a reactor able to transmute actinides effectively. In addition to the above applications the reactor produced electrical power with a high efficiency. A reactor subsystem was designed for each of the two applications. For the breeder reactor, neutronics calculations were carried out for a U-233 plasma core with a molten salt breeding blanket. A reactor was designed with a low critical mass (less than a few hundred kilograms U-233) and a breeding ratio of 1.01. The plasma core actinide transmutation reactor was designed to transmute the nuclear waste from conventional LWR's. The spent fuel is reprocessed during which 100% of Np, Am, Cm, and higher actinides are separated from the other components. These actinides are then manufactured as oxides into zirconium clad fuel rods and charged as fuel assemblies in the reflector region of the plasma core actinide transmutation reactor. In the equilibrium cycle, about 7% of the actinides are directly fissioned away, while about 31% are removed by reprocessing.

Clement, J. D.

Coupled Reactor and Engine Nuclear Thermal Propulsion Modeling Methodology

The design and development process of a Nuclear Thermal Propulsion (NTP) system requires extensive multiphysics modeling to couple the neutron physics and thermal feedback effects to determine the reactor’s power shape. Propulsion system performance codes utilize this power shape to determine NTP key performance parameters. While the power shape is heavily dependent on the temperature profile and geometry of the reactor, many analyses either assume a constant power shape, or use neutronics analysis to determine a power shape for a specific reactor configuration. The development of a coupling interface for a propulsion system performance code and a Monte Carlo neutron transport code (OpenMC) allows for the reactor power shape to be calculated in an iteration loop. The interface utilizes a file share system to transfer geometry dimensions, temperatures, and material identifiers to OpenMC, which is used to perform a neutron transport simulation of a design like the government Testing Reference Design reactor. The interface is then able to post-process the results from OpenMC and use the same file share system to share a power shape and other important neutron transport parameters to the system performance code. Initial results show that neglecting the changes to power shape when comparing reactor configurations can yield inaccurate results. Furthermore, utilizing propellants other than hydrogen gas can cause significant changes to the power shape, and thus, the thermal performance of a specific reactor design. This methodology is being expanded to allow for multiple families of NTP reactors to be analyzed, including block moderator, particle bed, and NERVA-derived reactors.

multiphysics coupling

Coupled Reactor and Engine Nuclear Thermal Propulsion Modeling Methodology

The design and development process of a Nuclear Thermal Propulsion (NTP) system requires extensive multiphysics modeling to couple the neutron physics and thermal feedback effects to determine the reactor’s power shape. Propulsion system performance codes utilize this power shape to determine NTP key performance parameters. While the power shape is heavily dependent on the temperature profile and geometry of the reactor, many analyses either assume a constant power shape, or use neutronics analysis to determine a power shape for a specific reactor configuration. The development of a coupling interface for a propulsion system performance code and a Monte Carlo neutron transport code (OpenMC) allows for the reactor power shape to be calculated in an iteration loop. The interface utilizes a file share system to transfer geometry dimensions, temperatures, and material identifiers to OpenMC, which is used to perform a neutron transport simulation of a design like the government Testing Reference Design reactor. The interface is then able to post-process the results from OpenMC and use the same file share system to share a power shape and other important neutron transport parameters to the system performance code. Initial results show that neglecting the changes to power shape when comparing reactor configurations can yield inaccurate results. Furthermore, utilizing propellants other than hydrogen gas can cause significant changes to the power shape, and thus, the thermal performance of a specific reactor design. This methodology is being expanded to allow for multiple families of NTP reactors to be analyzed, including block moderator, particle bed, and NERVA-derived reactors.

multiphysics coupling

Vacuum-assisted carbon molecular sieve membrane reactor for non-oxidative ethane dehydrogenation

Non-oxidative ethane dehydrogenation (EDH) is equilibrium-limited and endothermic. Selective hydrogen removal using a gas-permeable membrane within the EDH reaction zone can overcome the thermodynamic equilibrium, enabling higher ethane conversions. Employing vacuum as the permeation driving force, rather than a sweep gas, enhances the industrial viability of membrane reactors by eliminating additional post-reaction separation units. This study presents a membrane reactor that integrates H 2 -permeable carbon molecular sieve (CMS) hollow fiber membranes embedded in a fixed bed of cobalt in a dealuminated beta zeolite (Co@DeAl-BEA) catalyst, utilizing a vacuum to remove hydrogen efficiently. The CMS membrane exhibits high hydrogen permeance and an excellent H 2 /C 2 H 6 separation factor. The membrane reactor significantly enhanced the ethane conversion under reaction conditions comparable to those reported in the literature. A Langmuir-Hinshelwood kinetic rate expression was developed and incorporated into a one-dimensional steady-state reactor model. The experimentally validated model indicates that increasing the number of hollow fibers improves ethane conversion, although ethane loss to the permeate limits the benefit. The contact area between the catalyst and the membrane limits the reactor performance more than the catalytic throughput. Furthermore, we find that the location of the catalyst packing relative to the hollow fiber membranes influences ethane loss and conversion. Higher reactor pressures and inlet ethane flow rates improve space-time yield at the expense of lower ethane conversion. Increasing reactor temperature or packing length promotes both performance metrics. The EDH membrane reactor demonstrated durability over 200 h of continuous operation, maintaining record-low deactivation rates and high ethylene selectivity. Protocols for catalyst regeneration were developed.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI

Technoeconomic Design Optimization for Fast Reactors. Part I: Workflow Development and Case Study for Small LFR District Energy Application

The nuclear industry is developing small reactor designs that can target a variety of deployment locations and energy products. Smaller nuclear designs have traditionally struggled to handle the steep trade-offs between size and cost that have historically incentivized large reactors. This motivates computational optimization of small reactors to minimize costs and quantify the trade-off between size and cost. In this paper, the cost/size trade-off for a small fast reactor is derived using a multi-objective genetic algorithm optimization, with steady-state, transient, and cost analysis of the fast reactor being performed. Specifically, the method is demonstrated on a small 10- to 120-MW(thermal) U-Pu-Zr–fueled lead-cooled fast reactor with a 10-year core life for district energy applications, which can have a thermal load compatible with this range. The results reinforced that fast reactor cores at the lower end of this power range suffer cost penalties due to critical mass considerations. It was found that high power density cores with strong reactivity swings and many control rods were favored over designing to minimize reactivity swing. Furthermore, this contrasts with some traditional configurations designed using engineering judgment and demonstrates that optimizers can find nontraditional but realistic solutions, along with demonstrating the value of incorporating cost functions into whole-reactor design optimization.

Fast reactor

SCALE 6.3 Validation: Reactor Physics

This study was performed to validate the SCALE/Polaris v6.3.0–PARCS v3.4.2 code procedure with the ENDF/B-VII.1 AMPX 56-group library for light-water reactor analysis by comparing the simulated results with the measured data for critical experiments and operating light-water reactors. Uncertainties of the SCALE/Polaris–PARCS code procedure for light-water reactor physics analysis were evaluated in the validation for key nuclear parameters such as reactivity, control bank work, temperature coefficients, and pin and assembly power peaking factors. In addition, the SCALE/TRITON v6.3.1 procedure with the ENDF/B-VII.1 and VIII.0 252-group and continuous-energy cross sections was validated for non-lightwater reactors including the HTR-10 reactor, the High-Temperature Test Reactor, the Molten Salt Reactor Experiment, and the Experimental Breeder Reactor II.

22 GENERAL STUDIES OF NUCLEAR REACTORS

Improvement and Verification of Online Cross Section Generation Capability of Griffin for TRISO-fueled Reactors

Griffin, a MOOSE-based reactor multiphysics code jointly developed by Idaho National Laboratory and Argonne National Laboratory under the DOE Office of Nuclear Energy’s NEAMS program, has pursued the development of an online multigroup cross section generation capability for a few years to enable high-fidelity, problem-dependent neutronics analyses of advanced thermal reactors. Recent advancements in Griffin’s online multigroup cross section generation capability have significantly improved the accuracy, robustness, and efficiency of self-shielding calculations for both prismatic and pebble-bed TRISO-fueled reactor applications. Key developments include a unified fuel self-shielding method applicable to both TRISO and annular compact/spherical shell fuel zone geometries; an advanced Dancoff Category-based Equivalence Theory using a bell function for non-fuel resonance treatment, achieving more than an order-of-magnitude speedup compared to the Tone method; an on-the-fly multigroup equivalence approach to mitigate group condensation errors; and a streaming correction method for pebble-bed homogenization. A proof-of-concept demonstration of on-the-fly group condensation with consistent P0 transport correction was also achieved. The method reproduced direct fine-group solutions with excellent accuracy (eigenvalue errors within 10 pcm and pin-power differences within 0.5%), but due to performance limitations of the current fixed-source solver, improvements to solver efficiency will be addressed in future work. Verification tests were performed on graphite-moderated TRISO-fueled two-dimensional core benchmark problems representing gas-cooled microreactors, heat pipe-cooled microreactors, gas-cooled pebble-bed reactors, and fluoride salt-cooled high-temperature reactors. Across all cases, Griffin showed excellent agreement with Serpent2 continuous energy Monte Carlo solutions: eigenvalue errors within 200 pcm, pin-power root-mean-square errors within 2%, and control rod and drum worth errors less than 2%. It should be noted that, for the benchmark problem, cross section generation contributed less than 3% of the total simulation times. These results demonstrate that Griffin’s online cross section generation capability delivers accurate and efficient reactor physics solutions across a wide spectrum of TRISO-fueled advanced reactor designs. With further improvements to the fine-group fixed-source solver and planned extensions to depletion, transients, and coupled neutron–gamma transport, Griffin will be well-positioned to become a powerful and comprehensive tool for advanced reactor analysis.

22 GENERAL STUDIES OF NUCLEAR REACTORS

Developing Source Term Database for Advanced Reactors

A source term database is crucial to informing nuclear emergency response measures, enabling emergency responders to assess the potential severity of nuclear and radiological consequences. In recent times, various advanced reactor designs have come into operation, are under construction, or are being designed and developed. This report documents an effort carried out to develop a source term database for advanced reactors. The report covers key design features of these reactors and discusses radioactivity buildup and source term inventories of dose-significant radionuclides in the reactor core. For neutronic and depletion analyses, we used the SCALE code system, a computational suite for reactor physics, depletion, criticality, and sensitivity/uncertainty quantification. We used SCALE/TRITON to perform depletion calculations to predict cycle length and discharge burnup and to generate the ORIGEN reactor library. Subsequently, we used SCALE/ORIGAMI to calculate radioactivity buildup and, thereby, the source term inventories at the targeted discharge burnup, using the ENDF/B-VII.1 nuclear data library. This report covers several advanced reactors, including the KLT-40S, RITM-200N, VOYGR, and eVinci. However, other reactors, such as the RITM-200S and ARC-100, have yet to be investigated and will be explored in future efforts.

22 GENERAL STUDIES OF NUCLEAR REACTORS

Modeling a Packed Bed Reactor Utilizing the Sabatier Process

A numerical model is being developed using Python which characterizes the conversion and temperature profiles of a packed bed reactor (PBR) that utilizes the Sabatier process; the reaction produces methane and water from carbon dioxide and hydrogen. While the specific kinetics of the Sabatier reaction on the RuAl2O3 catalyst pellets are unknown, an empirical reaction rate equation1 is used for the overall reaction. As this reaction is highly exothermic, proper thermal control is of the utmost importance to ensure maximum conversion and to avoid reactor runaway. It is therefore necessary to determine what wall temperature profile will ensure safe and efficient operation of the reactor. This wall temperature will be maintained by active thermal controls on the outer surface of the reactor. Two cylindrical PBRs are currently being tested experimentally and will be used for validation of the Python model. They are similar in design except one of them is larger and incorporates a preheat loop by feeding the reactant gas through a pipe along the center of the catalyst bed. The further complexity of adding a preheat pipe to the model to mimic the larger reactor is yet to be implemented and validated; preliminary validation is done using the smaller PBR with no reactant preheating. When mapping experimental values of the wall temperature from the smaller PBR into the Python model, a good approximation of the total conversion and temperature profile has been achieved. A separate CFD model incorporates more complex three-dimensional effects by including the solid catalyst pellets within the domain. The goal is to improve the Python model to the point where the results of other reactor geometry can be reasonably predicted relatively quickly when compared to the much more computationally expensive CFD approach. Once a reactor size is narrowed down using the Python approach, CFD will be used to generate a more thorough prediction of the reactors performance.

Reactor

Coupled Reactor Multiphysics and Mass Scalability Assessment for Crewed Megawatt-Class NEP System Architectures

Nuclear Electric Propulsion (NEP) is an in-space propulsion technology capable of enabling opposition and conjunction class crewed Mars missions. NEP subsystems include the reactor for heat generation, a power conversion system (PCS), power management and distribution (PMAD), electric propulsion subsystem (EPS), and a primary heat rejection system. Specific mass, or αe (kg / kWe), is a key performance parameter (KPP) of the propulsion system which is directly scalable with the performance and mass estimates for individual components. To inform technology maturation planning, full system and component level parametric modeling is ongoing to explore the design trade space and illustrate the effect of subsystem design choices on the system KPPs. In this study, scaling of high-assay, low-enriched uranium (HALEU) reactor designs is assessed through coupled reactor physics and thermal hydraulics analyses. Scaling analyses evaluate the impact of system performance parameters (power level, interface temperatures) on mass for direct gas cooled, pumped liquid metal, and passively cooled heat pipe reactor concepts. Each concept requires specific geometries and working fluids to reach the performance goals of PCS interface conditions (temperature, pressure, flow rate) and system mass. The reactor assembly includes the active core (fuel, moderator, cladding, working fluid), axial and radial neutron reflectors, control drums, structural support / pressure vessel, and external radiation shielding. Each of these components are parametrically sized based on performance parameters for a megawatt-class power cycle. Results of this scaling analysis increase NEP propulsion system modeling fidelity and ultimately aim to support technology down-selection along with related technology development planning. The reactor and shield αe are a function of several PCS design choices, and reactor scaling with these parameters must be considered to enable an informed decision on reactor geometry and working fluid combination.

Nuclear Electric Propulsion

Coupled Reactor Multiphysics and Mass Scalability Assessment for Crewed Megawatt-Class NEP System Architectures

Nuclear Electric Propulsion (NEP) is an in-space propulsion technology capable of enabling opposition and conjunction class crewed Mars missions. NEP subsystems include the reactor for heat generation, a power conversion system (PCS), power management and distribution, electric propulsion system, and heat rejection system. Specific mass, or α (kg/kWe), is a key performance parameter (KPP) of the propulsion system which is directly scalable with the performance and mass predictions for each individual component. To inform technology maturation planning activities, full system and component level parametric modeling is ongoing to explore the design trade space and illustrate the effect of subsystem design choices on the system KPPs. In this study, scaling of high-assay, low-enriched uranium reactor designs is assessed through coupled reactor physics and thermal hydraulics analyses. Scaling analyses evaluate the impact of system performance parameters (power level, interface temperatures) on mass for direct gas cooled, pumped liquid metal, and passively-cooled heat pipe reactor concepts. Each concept requires specific geometries, fluids, and power conversion interface conditions (temperature, pressure, flow rate) to meet desired performance and mass. The reactor assembly includes the active core (fuel, moderator, cladding, working fluid), axial and radial neutron reflectors, control drums, structural support / pressure vessel, and external radiation shielding. Each of these components are parametrically sized based on performance parameters for a megawatt-class power cycle. Results of this scaling analysis increase NEP propulsion system modeling fidelity and ultimately aim to support concept down-selection along with related technology development planning. The reactor and shield α are a function of several PCS and heat rejection system design choices, and reactor scaling with these parameters must be considered to enable an informed decision on an optimal reactor geometry and working fluid combination.

Nuclear Electric Propulsion

CFD Validation of a Methanation Reactor With Cooling Jacket Used in an ISRU Propellant Production System

A digital model of an OxEon Energy methanation reactor (with a cooling jacket) planned for use in a NASA designed In-Situ Resource Utilization (ISRU) propellant production system was developed by the UTEP Aerospace Center using STAR-CCM+ software, a Siemens Computational Fluid Dynamics (CFD) package. The reactor digital model created includes two primary components: (i) a cooling jacket that maintains the required reactor chamber surface temperature and (ii) a reactor core chamber where the exothermic methanation reactions occur. The StarCCM+ digital model has two objectives: (i) maintain the surface temperature of the reactor below a set value so that the extra heat accumulation does not melt or crack the chamber, and (ii) produce the desired amount of methane (CH4) according to test cases completed by OxEon Energy for NASA JSC. The computational digital model was validated with the test data enabling the cooling jacket scale to be predicted based on the reactor's operating temperature. The StarCCM+ simulation of the digital model shows a good agreement with the cooling jacket performance and the CH4 production. One test case for the cooling jacket maintained the reactor wall temperature around 310 °C was validated and the reactor produced 185 g/hr CH 4 .

ISRU

Multi-cycle reload analysis of a long cycle gas-cooled fast modular reactor

There is currently significant interest in deploying HALEU-fueled fast reactors, including the General Atomics (GA) Fast Modular Reactor (FMR). Such reactors can achieve very long fuel cycles, but with multi-batch loading will take decades to reach equilibrium. This motivates design and analysis of both the initial core and multi-cycle reload, which is typically performed using fast-running, deterministic fast reactor codes such as the Argonne Reactor Computation (ARC) codes. In this paper, multicycle reload of the GA FMR is analyzed using the ARC codes. The GA FMR utilizes 19.75 % enriched fuel in a 16 year cycle with a three-batch strategy, with twice-burned fuel placed on the core periphery. The GA FMR has a softened neutron spectrum due to reflecting elements in the core, so the neutronic solution is first benchmarked against the OpenMC Monte Carlo code. Discrepancy on k eff is 400–600 pcm, likely due to the softened neutron spectrum, heterogeneous fuel assembly design and central reflector. However, the rms discrepancy on the assembly power distribution is only 0.6 %, despite the presence of the central reflector. A reload strategy is devised for the first three cycles of such a reactor, ultimately spanning the first 45–48 years of its operation. The fresh core uses 19.75 %, 19.25 % and 16.75 % enriched fuel in place of fresh, once-burned and twice-burned and is then subsequently refueled with only 19.75 % enriched fuel. The cycle length is varied over 3 cycles of operation to balance fuel utilization and reactor availability, specifically with use of an extended 18-year Cycle 1, followed by a shortened 11-year Cycle 2. Cycle 3 is close to the target 16-year length. Finally, placing twice burned assemblies next to the GA FMR central reflector can reduce power peaking by 3 %, at the expense of slightly reducing the cycle length.

22 GENERAL STUDIES OF NUCLEAR REACTORS

Technoeconomic Design Optimization for Fast Reactors. Part II: Impact of Technoeconomic Constraints on Optimal Design

There is a current drive toward optimizing reactors, particularly small/micro reactors to minimize cost and maximize performance. Previous work has investigated the development of technoeconomic workflows for the design optimization of pool-type fast reactors that aim to deploy into district energy grids. Initial scoping studies verified that the workflow was capable of capturing design trends throughout a variety of design configurations and problem formulations while remaining sufficiently flexible. In this paper, this methodology is applied to understand how cost functions and technoeconomic constraints can drive optimal reactor design. Specifically, the UPu10Zr-fueled fast reactor model from Part I is adapted to include changes in the fissile content limits, control rod worth limits, control rod drive cost, and assumed fuel form. In the case of constraint relaxation at fixed power (fissile content and control rod worth limits), cost sensitivities of 5% to 10% were uncovered. Multi-objective optimization at varying reactor power levels with individualized control rod drives for each assembly (as opposed to one operational and one safety drive) increased cost by $\$10$ to $\$25$ million and substantially altered the optimal core geometry, favoring geometries with substantially fewer control rod placements relative to baseline optimization. Finally, a multi-objective optimization was performed at varying power levels with the fuel form overhauled to metallic, high-assay low-enriched uranium–based U10Zr with more refined fuel cost models. In the case of uranium fueling, the costs increased by at least $50 million relative to the baseline case. Furthermore, economic fuel zoning and lower reactivity swing cores were recovered. Each case serves to demonstrate the value of applying technoeconomic workflows to initial reactor design scoping studies to better understand the trade-off for a proposed concept between different design options.

Argonne Reactor Computation (ARC) codes

Engineering design of a kW-scale continuous reactor-heat exchanger for high temperature discharge of particle-based thermochemical energy storage

This study investigates the theoretical design parameters and thermal performance of a kW-scale continuous oxidation reactor for high temperature (~1000 °C) thermochemical energy storage (TCES) applications. The concept comprises a counter-current particle-based system that includes a reaction zone with a heat exchanger to extract the heat produced from the oxidation reaction. Both above and below the hot reactive volume are sensible heat recuperation zones to enable the feed and removal of particles and oxidizing gas near ambient temperature during steady state operation. Two operation types for the reaction zone are studied, a fluidized bed reactor (FBR) and a moving bed reactor (MBR). The results of the parametric analysis suggest that the MBR requires a smaller volume per kW of heat produced, achieving power densities in excess of 2500 kW/m3 compared to ~ 900 kW/m 3 in the FBR. Additionally, the MBR achieves between 0.71 and 0.99 oxidation conversions compared to between 0.23 and 0.38 conversions in the FBR with the same volumes and flowrates. However, the FBR has the potential to maintain a uniform reactor temperature which can produce heat transfer fluid (HTF) outlet temperatures as high as the reactor temperature, i.e., ~1000 °C, whereas the MBR produces variable reactor temperatures that can create overheating zones and low HTF outlet temperatures (< 800 °C) depending on the operating conditions selected. Future work should aim at understanding the coupled fluid dynamics, heat and mass transfer, and thermochemical reaction for any given combination of reactor volume and contacting patterns. Here, these studies should be complemented by experimental work on particle-gas TCES reactors.

25 ENERGY STORAGE

Renewable low-Z wall for fusion reactors with built-in tritium recovery (Final Technical Report)

This project pursued development of a novel renewable plasma-facing wall technology for fusion reactors. The technology is based on a slurry which can be easily delivered by delivery tubes to the reactor wall. The slurry dries at the hot reactor wall into pebble rods which are extruded out into the hot plasma, where the pebbles break off and fall along the reactor wall and can be recovered by gravity and re-used. The falling pebbles carry away heat and tritium and also protect the wall against large scale erosion or redeposition of material. The research focused on carbon-based pebble rods and demonstrated that pebble rods could be produced from slurry with tolerable levels of outgassing on a reactor-relevant timescale (< 5 minutes). Steady-state handling of reactor relevant (up to 50 MW/m 2 ) normal-incidence heat loads was demonstrated. Pebble release velocities were found to be sufficiently small (< 1 m/s) to allow recovery below the vacuum chamber. Tunability of the pebble rod breaking rate was demonstrated by changing the fill fraction of the interpebble matrix which binds the pebbles together. This work could benefit the public by helping move forward the design of commercially viable fusion energy reactors. Designing a first wall for magnetic fusion reactors which can handle the huge heat loads present and also avoid buildup of tritium-containing deposits is extremely challenging and requires novel approaches like the one being investigated here.

36 MATERIALS SCIENCE

Applicability of 100kWe-class of space reactor power systems to NASA manned space station missions

An assessment is made of a manned space station operating with sufficiently high power demands to require a multihundred kilowatt range electrical power system. The nuclear reactor is a competitor for supplying this power level. Load levels were selected at 150kWe and 300kWe. Interactions among the reactor electrical power system, the manned space station, the space transportation system, and the mission were evaluated. The reactor shield and the conversion equipment were assumed to be in different positions with respect to the station; on board, tethered, and on a free flyer platform. Mission analyses showed that the free flyer concept resulted in unacceptable costs and technical problems. The tethered reactor providing power to an electrolyzer for regenerative fuel cells on the space station, results in a minimum weight shield and can be designed to release the reactor power section so that it moves to a high altitude orbit where the decay period is at least 300 years. Placing the reactor on the station, on a structural boom is an attractive design, but heavier than the long tethered reactor design because of the shield weight for manned activity near the reactor.

Silverman, S. W.

A Computational Fluid Dynamic and Heat Transfer Model for Gaseous Core and Gas Cooled Space Power and Propulsion Reactors

A computational model based on the axisymmetric, thin-layer Navier-Stokes equations is developed to predict the convective, radiation and conductive heat transfer in high temperature space nuclear reactors. An implicit-explicit, finite volume, MacCormack method in conjunction with the Gauss-Seidel line iteration procedure is utilized to solve the thermal and fluid governing equations. Simulation of coolant and propellant flows in these reactors involves the subsonic and supersonic flows of hydrogen, helium and uranium tetrafluoride under variable boundary conditions. An enthalpy-rebalancing scheme is developed and implemented to enhance and accelerate the rate of convergence when a wall heat flux boundary condition is used. The model also incorporated the Baldwin and Lomax two-layer algebraic turbulence scheme for the calculation of the turbulent kinetic energy and eddy diffusivity of energy. The Rosseland diffusion approximation is used to simulate the radiative energy transfer in the optically thick environment of gas core reactors. The computational model is benchmarked with experimental data on flow separation angle and drag force acting on a suspended sphere in a cylindrical tube. The heat transfer is validated by comparing the computed results with the standard heat transfer correlations predictions. The model is used to simulate flow and heat transfer under a variety of design conditions. The effect of internal heat generation on the heat transfer in the gas core reactors is examined for a variety of power densities, 100 W/cc, 500 W/cc and 1000 W/cc. The maximum temperature, corresponding with the heat generation rates, are 2150 K, 2750 K and 3550 K, respectively. This analysis shows that the maximum temperature is strongly dependent on the value of heat generation rate. It also indicates that a heat generation rate higher than 1000 W/cc is necessary to maintain the gas temperature at about 3500 K, which is typical design temperature required to achieve high efficiency in the gas core reactors. The model is also used to predict the convective and radiation heat fluxes for the gas core reactors. The maximum value of heat flux occurs at the exit of the reactor core. Radiation heat flux increases with higher wall temperature. This behavior is due to the fact that the radiative heat flux is strongly dependent on wall temperature. This study also found that at temperature close to 3500 K the radiative heat flux is comparable with the convective heat flux in a uranium fluoride failed gas core reactor.

Anghaie, S.