Search NASA⌕ Search

SEARCH · Search NASA

Results for “CFD simulation”

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 325 records · Page 18

Performance Improvement Through Indexing of Turbine Airfoils: Numerical Simulation - Part 2

An experimental/analytical study has been conducted to determine the performance improvements achievable by circumferentially indexing succeeding rows of turbine stator airfoils. A series of tests was conducted to experimentally investigate stator wake clocking effects on the performance of the space shuttle main engine (SSME) alternate turbopump development (ATD) fuel turbine test article (TTA). The results from this study indicate that significant increases in stage efficiency can be attained through application of this airfoil clocking concept. Details of the experiment and its results are documented in part 1 of this paper. In order to gain insight into the mechanisms of the performance improvement, extensive computational fluid dynamics (CFD) simulations were executed. The subject of the present paper is the initial results from the CFD investigation of the configurations and conditions detailed in part 1 of the paper. To characterize the aerodynamic environments in the experimental test series, two-dimensional (2D), time accurate, multistage, viscous analyses were performed at the TTA midspan. Computational analyses for five different circumferential positions of the first stage stator have been completed. Details of the computational procedure and the results are presented. The analytical results verify the experimentally demonstrated performance improvement and are compared with data whenever possible. Predictions of time-averaged turbine efficiencies as well as gas conditions throughout the flow field are presented. An initial understanding of the turbine performance improvement mechanism based on the results from this investigation is described.

Griffin, Lisa W.↗

Numerical Modeling of a Two-Stage Ocean Current Turbine

The Equinox Ocean Turbines (EQOT) current energy converter has a unique design with power generation in two small-diameter turbines attached to the tips of a large-diameter passive rotor. This configuration offers some key advantages for capturing ocean currents. With no centrally placed generator, almost no reaction torque is required at the nacelle of the main large-diameter rotor, and the small-diameter tip turbine generators operate at a higher speed and lower torque. The physics that determine the performance and loads on the turbine are also unique. The interactions of the flow field between the two stages and the general architecture of the system cannot be captured with traditional mid-fidelity modeling tools. For design iterations and large sets of load cases, it is important to have mid-fidelity models that can capture the important phenomenon with enough accuracy to identify global trends. This work uses a limited set of high-fidelity computational fluid dynamics (CFD) simulations to help inform the selection of and construction of a custom mid-fidelity model. Mid-fidelity modeling approaches were verified by comparing key turbine performance quantities to those found with the CFD model. Hydrodynamic interactions of the two-stage rotor were identified through high-fidelity CFD modeling. This highlighted the impact of the main rotor tip vortex and wake on the secondary rotor apparent inflow. This results in a relative flow rotation and sharp deficit, that change the optimal secondary rotor rotation speed and adds unsteadiness to the blade loading respectively. Multiple mid-fidelity approaches were evaluated for their ability to capture these effects. A simple approximation of the combined-stage performance based on single-stage BEM provides a reasonable rough prediction, especially near the peak TSR values, with some larger discrepancy at higher TSRs. Predicting the combined-stage performance based on single-stage CFD data improves this prediction across the TSR range. Although the combined-stage modeling in OLAF was not successful in this stage of the project, it showed promise as a mid-fidelity method, assuming the parameters can be tuned to account for the significant differences in time and length scales between the main and secondary rotors. This may be addressed through code changes in future work. A significant finding from the OLAF work was the agreement between the vortex core radius values found independently via a parameter space search and via CFD. The technique of using single-stage secondary rotor BEM, with a custom inflow taken from single-stage main rotor CFD or OLAF, provides an efficient method to capture one-way coupled flow interactions. This method provided generally good predictions of the impact of the flow rotation on the secondary rotor but struggled to accurately predict the peaks of the unsteady load progression. Future work could include some superposition of a tuned main rotor trailing edge viscous wake into the custom inflow to better predict this interaction.

16 TIDAL AND WAVE POWER↗

Aerothermodynamic Analyses of Towed Ballutes

A ballute (balloon-parachute) is an inflatable, aerodynamic drag device for application to planetary entry vehicles. Two challenging aspects of aerothermal simulation of towed ballutes are considered. The first challenge, simulation of a complete system including inflatable tethers and a trailing toroidal ballute, is addressed using the unstructured-grid, Navier-Stokes solver FUN3D. Auxiliary simulations of a semi-infinite cylinder using the rarefied flow, Direct Simulation Monte Carlo solver, DSV2, provide additional insight into limiting behavior of the aerothermal environment around tethers directly exposed to the free stream. Simulations reveal pressures higher than stagnation and corresponding large heating rates on the tether as it emerges from the spacecraft base flow and passes through the spacecraft bow shock. The footprint of the tether shock on the toroidal ballute is also subject to heating amplification. Design options to accommodate or reduce these environments are discussed. The second challenge addresses time-accurate simulation to detect the onset of unsteady flow interactions as a function of geometry and Reynolds number. Video of unsteady interactions measured in the Langley Aerothermodynamic Laboratory 20-Inch Mach 6 Air Tunnel and CFD simulations using the structured grid, Navier-Stokes solver LAURA are compared for flow over a rigid spacecraft-sting-toroid system. The experimental data provides qualitative information on the amplitude and onset of unsteady motion which is captured in the numerical simulations. The presence of severe unsteady fluid - structure interactions is undesirable and numerical simulation must be able to predict the onset of such motion.

Gnoffo, Peter A.↗

Attitude Reconstruction of Free-Flight CFD Generated Trajectories Using Non-Linear Pitch Damping Coefficient Curves

Attitude history reconstruction of Free-flight CFD generated trajectories with non-linear pitch damping coefficient curves is investigated. Free-flight CFD simulations of the capsule shape used for the Genesis sample return mission and the upcoming Dragonfly mission to Titan are conducted for 1-, 2-, and 3-degree-of-freedom cases. Two different data reduction methodologies are employed to derive a pitch damping curve as a function of instantaneous angle of attack. These curves are then used to reconstruct the attitude history of the body which is compared to the raw simulation results. While both data reduction methods produce pitch damping curves that can generally reconstruct the motion seen in the Free-flight simulations, it is found that optimization of the pitch damping curve using an inverse estimation process yields superior and more generalizable results. Further refinement of this technique could allow pitch damping curves derived using CFD to serve as a capability complementary to existing techniques for dynamic stability characterization.

entry↗

Attitude Reconstruction of Free-Flight CFD Generated Trajectories Using Non-Linear Pitch Damping Coefficient Curves

Attitude history reconstruction of Free-flight CFD generated trajectories with non-linear pitch damping coefficient curves is investigated. Free-flight CFD simulations of the capsule shape used for the Genesis sample return mission and the upcoming Dragonfly mission to Titan are conducted for 1-, 2-, and 3-degree-of-freedom cases. Two different data reduction methodologies are employed to derive a pitch damping curve as a function of instantaneous angle of attack. These curves are then used to reconstruct the attitude history of the body which is compared to the raw simulation results. While both data reduction methods produce pitch damping curves that can generally reconstruct the motion seen in the Free-flight simulations, it is found that optimization of the pitch damping curve using an inverse estimation process yields superior and more generalizable results. Further refinement of this technique could allow pitch damping curves derived using CFD to serve as a capability complementary to existing techniques for dynamic stability characterization.

entry↗

Development of a Digital Twin for Hydrogen Dispersion and Safety Assessment in an Electrolyzer Based Hydrogen Production Facility

Digital twin models are virtual representations of physical systems that use real-time data to simulate and optimize performance. This study presents the development and initial implementation of a digital twin (DT) for the electrolyzer-based hydrogen production facility at NREL's Advanced Research on Integrated Energy Systems (ARIES), focused on enhancing safety and optimizing sensor placement through physics-based simulations and metadata integration. The DT incorporates detailed facility-specific information, including component layout, leak locations, and controlled release parameters, to model hydrogen dispersion under varying environmental conditions. Using steady-state computational fluid dynamics (CFD) simulations informed by real meteorological data, such as wind speed, direction, and vertical wind profiles, the DT enables visualization of hydrogen plume behavior and spatial concentration distributions. Comparative analysis between high and low wind speed scenarios illustrates the significant influence of wind dynamics on plume shape and extent, with horizontal momentum dominating dispersion at higher speeds, while buoyancy effects become more prominent under low wind conditions. These simulations generate a rich dataset embedded within the DT, allowing users to assess potential leak outcomes and identify optimal sensor locations based on concentration thresholds. The model supports scenario-based analysis to guide safety strategies and equipment deployment for open-area hydrogen infrastructure. The digital twin thus serves as a dynamic platform for virtual prototyping, providing predictive insight into hydrogen behavior and enhancing risk-informed decision-making. This initial phase establishes a validated foundation for future integration of transient, uncontrolled leak scenarios and real-time sensor feedback, positioning the DT as a critical tool for safety design, operational planning, and adaptive monitoring in hydrogen systems. Overall, the approach demonstrates the value of combining environmental data with digital simulations to inform safer and more efficient deployment of hydrogen technologies.

08 HYDROGEN↗

Two-Pendulum Model of Propellant Slosh in Europa Clipper PMD Tank

The objective of this fluids analysis is to model propellant slosh for the Europa Clipper mission using a two-pendulum model, such that controls engineers can predict slosh behavior during the mission. Propellant slosh causes shifts in center of mass and exerts forces and torques on the spacecraft which, if not adequately controlled, can lead to mission failure. The two-pendulum model provides a computationally simple model that can be used to predict slosh for the Europa Clipper tank geometry. The Europa Clipper tank is cylindrical with a domed top and bottom and includes a propellant management device (PMD). Due to the lack of experimental data in low gravity environments, computational fluid dynamics (CFD) simulation results were used as 'real' slosh behavior for two propellants at three fill fractions. Key pendulum parameters were derived that allow the pendulum model's center of mass, forces, and moments to closely match the CFD data. The parameter trends were examined as a function of tank fill fraction and compared with solutions to analytic equations that describe the frequency of slosh in tanks with simple geometries. The trends were monotonic as expected, and parameters resembled analytical predictions; any differences could be explained by the specific differences in the geometry of the tank. This paper summarizes the new method developed at Goddard Space Flight Center (GSFC) for deriving pendulum parameters for two-pendulum equivalent sloshing models. It presents the results of this method and discusses the validity of the results. This analysis is at a completed stage and will be applied in the immediate future to the evolving tank geometry as Europa Clipper moves past its preliminary design review (PDR) phase.

Propellant Slosh↗

Optimizing Performance and Accuracy with LIFQWIK

For the 2020 Summer Internship period, optimization and analysis was performed of a new iteration of the LIFQWIK PLIF model, used to evaluate the accuracy of CFD simulations of fuel injector performance. This presentation details the accuracy and performance optimization steps taken and presents the results of that optimization.

CFD↗

Boundary-Layer Flow Simulations Over Ablating Woven Thermal Protection System Material

Spallation is the mechanical removal of small chunks of material gets removed typically due to high shear conditions of the flow field. This reduces the ability of the thermal protection system (TPS) material to protect the spacecraft as well as cause turbulence in the flow causing higher heating rates. In this work, we focus on the material removal through ablation and high shear flow within the boundary layer region of woven TPS material. Woven TPS (WTPS) material is the latest class of material developed by NASA, to be used within the next generation of space flights. They are complex interlocked weaves designed to create a rigid structure that is highly resistant to heat and can be easily designed and tailored for a wide variety of entry environments. Due to material removal resulting from chemical degradation, the structural integrity of TPS material is affected. Spallation occurs when this structurally compromised material is exposed to the high shear flow conditions within the boundary layer. In order to understand the spallation mechanism within WTPS material, we first perform the material removal simulations which occur primarily through oxidation to obtain the microstructure at various stages of degradation. These simulations are performed using the Porous Microstructure Analysis (PuMA) software developed at NASA Ames. The micro-structure geometry used within these simulations were generated artificially to be similar to the 3D weave architecture of MSR-EEV (Mars Sample Return - Earth Entry Vehicle). The various eroded TPS micro-structures are then subjected to the boundary layer flow conditions to obtain critical surface quantities which contribute to the structural failure mechanism such as heat flux, pressure, and shear stress. The direct simulation Monte Carlo (DSMC) methodology is used to perform these simulations in order to accurately capture the strong gradients within the high-temperature boundary layer flow over the intricate geometry of WTPS material. The boundary layer profile is directly taken from the Computational Fluid Dynamics (CFD) simulation and provided as boundary conditions to the DSMC inlet and outlet. Further, the variation of these properties as the microstructure undergoes changes due to oxidation is also investigated. Finally, these quantities are used as input in PuMA to understand the material expansion/compression and strain within the woven TPS geometry and help in developing a comprehensive spallation and structure failure model.

microstructure↗

Boundary-Layer Flow Simulations Over Ablating Woven Thermal Protection System Material

Spallation is the mechanical removal of small chunks of material gets removed typically due to high shear conditions of the flow field. This reduces the ability of the thermal protection system (TPS) material to protect the spacecraft as well as cause turbulence in the flow causing higher heating rates. In this work, we focus on the material removal through ablation and high shear flow within the boundary layer region of woven TPS material. Woven TPS (WTPS) material is the latest class of material developed by NASA, to be used within the next generation of space flights. They are complex interlocked weaves designed to create a rigid structure that is highly resistant to heat and can be easily designed and tailored for a wide variety of entry environments. Due to material removal resulting from chemical degradation, the structural integrity of TPS material is affected. Spallation occurs when this structurally compromised material is exposed to the high shear flow conditions within the boundary layer. In order to understand the spallation mechanism within WTPS material, we first perform the material removal simulations which occur primarily through oxidation to obtain the microstructure at various stages of degradation. These simulations are performed using the Porous Microstructure Analysis (PuMA) software developed at NASA Ames. The micro-structure geometry used within these simulations were generated artificially to be similar to the 3D weave architecture of MSR-EEV (Mars Sample Return - Earth Entry Vehicle). The various eroded TPS micro-structures are then subjected to the boundary layer flow conditions to obtain critical surface quantities which contribute to the structural failure mechanism such as heat flux, pressure, and shear stress. The direct simulation Monte Carlo (DSMC) methodology is used to perform these simulations in order to accurately capture the strong gradients within the high-temperature boundary layer flow over the intricate geometry of WTPS material. The boundary layer profile is directly taken from the Computational Fluid Dynamics (CFD) simulation and provided as boundary conditions to the DSMC inlet and outlet. Further, the variation of these properties as the microstructure undergoes changes due to oxidation is also investigated. Finally, these quantities are used as input in PuMA to understand the material expansion/compression and strain within the woven TPS geometry and help in developing a comprehensive spallation and structure failure model.

microstructure↗

Development of a Digital Twin for Hydrogen Dispersion and Safety Assessment in an Electrolyzer-Based Hydrogen Production Facility: Preprint

Digital twin models are virtual representations of physical systems that use real-time data to simulate and optimize performance. This study presents the development and initial implementation of a digital twin (DT) for the electrolyzer-based hydrogen production facility at the National Renewable Energy Laboratory (NREL)'s Advanced Research on Integrated Energy Systems (ARIES), focused on enhancing safety and optimizing sensor placement through physics-based simulations and metadata integration. The DT incorporates detailed facility-specific information, including component layout, leak locations, and controlled release parameters, to model hydrogen dispersion under varying environmental conditions. Using steady-state computational fluid dynamics (CFD) simulations informed by real meteorological data, such as wind speed, direction, and vertical wind profiles, the DT enables visualization of hydrogen plume behavior and spatial concentration distributions. Comparative analysis between high and low wind speed scenarios illustrates the significant influence of wind dynamics on plume shape and extent, with horizontal momentum dominating dispersion at higher speeds, while buoyancy effects become more prominent under low wind conditions. These simulations generate a rich dataset embedded within the DT, allowing users to assess potential leak outcomes and identify optimal sensor locations based on concentration thresholds. The model supports scenario-based analysis to guide safety strategies and equipment deployment for open-area hydrogen infrastructure. The digital twin thus serves as a dynamic platform for virtual prototyping, providing predictive insight into hydrogen behavior and enhancing risk-informed decision-making. This initial phase establishes a validated foundation for future integration of transient, uncontrolled leak scenarios and real-time sensor feedback, positioning the DT as a critical tool for safety design, operational planning, and adaptive monitoring in hydrogen systems. Overall, the approach demonstrates the value of combining environmental data with digital simulations to inform safer and more efficient deployment of hydrogen technologies.

08 HYDROGEN↗

Pressurization of a Flightweight, Liquid Hydrogen Tank: Evaporation & Condensation at a Liquid/Vapor Interface

This paper presents an analysis and simulation of evaporation and condensation at a motionless liquid/vapor interface. A 1-D model equation, emphasizing heat and mass transfer at the interface, is solved in two ways, and incorporated into a subgrid interface model within a CFD simulation. Simulation predictions are compared with experimental data from the CPST Engineering Design Unit tank, a cryogenic fluid management test tank in 1-g. The numerical challenge here is the physics of the liquid/vapor interface; pressurizing the ullage heats it by several degrees, and sets up an interfacial temperature gradient that transfers heat to the liquid phase-the rate limiting step of condensation is heat conducted through the liquid and vapor. This physics occurs in thin thermal layers O(1 mm) on either side of the interface which is resolved by the subgrid interface model. An accommodation coefficient of 1.0 is used in the simulations which is consistent with theory and measurements. This model is predictive of evaporation/condensation rates, that is, there is no parameter tuning.

Stewart, Mark E. M.↗

A Computational Fluid Dynamics Study of Swirling Flow Reduction by using Anti-vortex Baffle

An anti-vortex baffle is a liquid propellant management device placed adjacent to an outlet of the propellant tank. Its purpose is to substantially reduce or eliminate the formation of free surface dip and vortex, as well as prevent vapor ingestion into the outlet, as the liquid drains out through the flight. To design an effective anti-vortex baffle, Computational Fluid Dynamic (CFD) simulations were undertaken for the NASA Ares I vehicle LOX tank subjected to the simulated flight loads with and without the anti-vortex baffle. The Six Degree-Of-Freedom (6- DOF) dynamics experienced by the Crew Launch Vehicle (CLV) during ascent were modeled by modifying the momentum equations in a CFD code to accommodate the extra body forces from the maneuvering in a non-inertial frame. The present analysis found that due to large moments, the CLV maneuvering has significant impact on the vortical flow generation inside the tank. Roll maneuvering and side loading due to pitch and yaw are shown to induce swirling flow. The vortical flow due to roll is symmetrical with respect to the tank centerline, while those induced by pitch and yaw maneuverings showed two vortices side by side. The study found that without the anti-vortex baffle, the swirling flow caused surface dip during the late stage of drainage and hence early vapor ingestion. The flow can also be non-uniform in the drainage pipe as the secondary swirling flow velocity component can be as high as 10% of the draining velocity. An analysis of the vortex dynamics shows that the swirling flow in the drainage pipe during the Upper Stage burn is mainly the result of residual vortices inside the tank due to conservation of angular momentum. The study demonstrated that the swirling flow in the drainage pipe can be effectively suppressed by employing the anti-vortex baffle.

Yang, H. Q.↗

A Computational Fluid Dynamics Study of Swirling Flow Reduction by Using Anti-Vortex Baffle

An anti-vortex baffle is a liquid propellant management device placed adjacent to an outlet of the propellant tank. Its purpose is to substantially reduce or eliminate the formation of free surface dip and vortex, as well as prevent vapor ingestion into the outlet, as the liquid drains out through the flight. To design an effective anti-vortex baffle, Computational Fluid Dynamic (CFD) simulations were undertaken for the NASA Ares I vehicle LOX tank subjected to the simulated flight loads with and without the anti-vortex baffle. The Six Degree-Of-Freedom (6-DOF) dynamics experienced by the Crew Launch Vehicle (CLV) during ascent were modeled by modifying the momentum equations in a CFD code to accommodate the extra body forces from the maneuvering in a non-inertial frame. The present analysis found that due to large moments, the CLV maneuvering has a significant impact on the vortical flow generation inside the tank. Roll maneuvering and side loading due to pitch and yaw are shown to induce swirling flow. The vortical flow due to roll is symmetrical with respect to the tank centerline, while those induced by pitch and yaw maneuverings showed two vortices side by side. The study found that without the anti-vortex baffle, the swirling flow caused surface dip during the late stage of drainage and hence early vapor ingestion. The flow can also be non-uniform in the drainage pipe as the secondary swirling flow velocity component can be as high as 10% of the draining velocity. An analysis of the vortex dynamics shows that the swirling flow in the drainage pipe during the Upper Stage burn is mainly the result of residual vortices inside the tank due to the conservation of angular momentum. The study demonstrated that the swirling flow in the drainage pipe can be effectively suppressed by employing the anti-vortex baffle.

Yang, H. Q.↗

A Computational Fluid Dynamics Study of Swirling Flow Reduction by Using Anti-Vortex Baffle

An anti-vortex baffle is a liquid propellant management device placed adjacent to an outlet of the propellant tank. Its purpose is to substantially reduce or eliminate the formation of free surface dip and vortex, as well as prevent vapor ingestion into the outlet, as the liquid drains out through the flight. To design an effective anti-vortex baffle, Computational Fluid Dynamic (CFD) simulations were undertaken for the NASA Ares I vehicle LOX tank subjected to the simulated flight loads with and without the anti-vortex baffle. The Six Degree-Of-Freedom (6-DOF) dynamics experienced by the Crew Launch Vehicle (CLV) during ascent were modeled by modifying the momentum equations in a CFD code to accommodate the extra body forces from the maneuvering in a non-inertial frame. The present analysis found that due to large moments, the CLV maneuvering has a significant impact on the vortical flow generation inside the tank. Roll maneuvering and side loading due to pitch and yaw are shown to induce swirling flow. The vortical flow due to roll is symmetrical with respect to the tank centerline, while those induced by pitch and yaw maneuverings showed two vortices side by side. The study found that without the anti-vortex baffle, the swirling flow caused surface dip during the late stage of drainage and hence early vapor ingestion. The flow can also be non-uniform in the drainage pipe as the secondary swirling flow velocity component can be as high as 10% of the draining velocity. An analysis of the vortex dynamics shows that the swirling flow in the drainage pipe during the Upper Stage burn is mainly the result of residual vortices inside the tank due to the conservation of angular momentum. The study demonstrated that the swirling flow in the drainage pipe can be effectively suppressed by employing the anti-vortex baffle.

Yang, H. Q.↗

Semi-Supervised Learning of Lift Optimization of Multi-Element Three-Segment Variable Camber Airfoil

This chapter describes a new intelligent platform for learning optimal designs of morphing wings based on Variable Camber Continuous Trailing Edge Flaps (VCCTEF) in conjunction with a leading edge flap called the Variable Camber Krueger (VCK). The new platform consists of a Computational Fluid Dynamics (CFD) methodology coupled with a semi-supervised learning methodology. The CFD component of the intelligent platform comprises of a full Navier-Stokes solution capability (NASA OVERFLOW solver with Spalart-Allmaras turbulence model) that computes flow over a tri-element inboard NASA Generic Transport Model (GTM) wing section. Various VCCTEF/VCK settings and configurations were considered to explore optimal design for high-lift flight during take-off and landing. To determine globally optimal design of such a system, an extremely large set of CFD simulations is needed. This is not feasible to achieve in practice. To alleviate this problem, a recourse was taken to a semi-supervised learning (SSL) methodology, which is based on manifold regularization techniques. A reasonable space of CFD solutions was populated and then the SSL methodology was used to fit this manifold in its entirety, including the gaps in the manifold where there were no CFD solutions available. The SSL methodology in conjunction with an elastodynamic solver (FiDDLE) was demonstrated in an earlier study involving structural health monitoring. These CFD-SSL methodologies define the new intelligent platform that forms the basis for our search for optimal design of wings. Although the present platform can be used in various other design and operational problems in engineering, this chapter focuses on the high-lift study of the VCK-VCCTEF system. Top few candidate design configurations were identified by solving the CFD problem in a small subset of the design space. The SSL component was trained on the design space, and was then used in a predictive mode to populate a selected set of test points outside of the given design space. The new design test space thus populated was evaluated by using the CFD component by determining the error between the SSL predictions and the true (CFD) solutions, which was found to be small. This demonstrates the proposed CFD-SSL methodologies for isolating the best design of the VCK-VCCTEF system, and it holds promise for quantitatively identifying best designs of flight systems, in general.

Intelligent Systems↗

COMPUTATION FLUID DYNAMICS ANALYSIS FOR GENERIC SMALL MODULAR REACTOR CONTAINMENT SEPARATE EFFECTS TEST

It is desirable for fourth-generation Small Modular Reactors to be passively cooled in standard and accident operations. Passive Containment Cooling Systems can reject heat from the containment structure, without using pumps or blowers. The targeted design containment structure is a large, domed, stainless steel, cylindrical vessel. In a postulated Design Basis Accident, steam will flash inside containment. Steam condensation occurs on the inner containment wall and transfers heat through the steel containment into a large body of water known as the annular reservoir (AR) surrounding the vessel serving as the ultimate heat sink. Natural circulation drives the flow in the AR and heat will be released to the environment by evaporation of water. Unique containment geometry requires a separate effects test (SET) facility for the verification and validation of the computer code and evaluation model development and assessment for reactor licensing efforts. In this study, STAR-CCM+, a computational fluid dynamics (CFD) code was used to inform the decision-making process on the design of the SET. The CFD simulation modeled, a two-phase turbulent flow with fluid film development and heat transfer for different containment geometries. The Reactor Excursion and Leak Analysis Program will also be used in a code-to-code verification against the CFD results.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

Computational Fluid Dynamics Simulations of Glass Vitrification Refractory Coupon Tests

The Waste Treatment and Immobilization Plant (WTP) at the Hanford site is nearing the start of the Direct-Feed Low-Activity Waste (DFLAW) operations. DFLAW is destined to convert a pretreated low activity waste portion of the 56 million gallons of tank waste into a stable solid glass. In the subsequent decade completion of the high-level waste (HLW) facility is anticipated. Sustained operational missions of both LAW and HLW melter facilities are expected over multiple decades. In high-temperature glass melters, the refractory lining corrodes over time, which could potentially be an issue for longer term operations, this refractory corrosion is higher at the level of the glass-air interface due to surface tension driven flow. The glass viscosity, melt pool temperature, and glass chemical composition can impact the rate at which the refractory corrodes. This rate is important to quantify for the various waste glasses to be produced at the WTP since the integrity of the refractory should not be a limiting factor affecting the lifetime of the melter. To this end, a series of glasses representative of the first batches of waste glass produced by the WTP will be melted in small-scale crucibles with Monofrax® K-3 coupons inserted. The corrosion of the K-3 will be measured in the melt and at the meltline (or neckline). A model for the corrosion rate will be constructed and implemented into a previously developed framework for a computational fluid dynamics (CFD) model of the full-scale WTP. To assist with experimental design and validate the implementation of the model in the full-scale melter, CFD simulations of the small-scale crucible tests were performed. The bubbling that occurs in the small-scale crucible is initially validated here with a model that uses silicone oil at room temperature. The viscosity of the oil ranges from 1 to 100 Pa•s, which corresponds to operating glass pool temperatures near 1150 °C down to idling temperatures near 950 °C. The simulation results show good agreement with the bubble sizes that form during experiments. CFD modeling of the crucible setup was used to determine bubbling characteristics to match the range of near-wall velocities expected in the full-scale WTP. This study presents the initial CFD modeling results, corrosion testing plan, and some preliminary corrosion samples with an outline for the next steps for the development of the corrosion model.

Abboud, Alexander W. [Idaho National Lab]↗