Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “CFD”

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 1,585 records · Page 88

Improving Fidelity of Launch Vehicle Liftoff Acoustic Simulations

Launch vehicles experience high acoustic loads during ignition and liftoff affected by the interaction of rocket plume generated acoustic waves with launch pad structures. Application of highly parallelized Computational Fluid Dynamics (CFD) analysis tools optimized for application on the NAS computer systems such as the Loci/CHEM program now enable simulation of time-accurate, turbulent, multi-species plume formation and interaction with launch pad geometry and capture the generation of acoustic noise at the source regions in the plume shear layers and impingement regions. These CFD solvers are robust in capturing the acoustic fluctuations, but they are too dissipative to accurately resolve the propagation of the acoustic waves throughout the launch environment domain along the vehicle. A hybrid Computational Fluid Dynamics and Computational Aero-Acoustics (CFD/CAA) modeling framework has been developed to improve such liftoff acoustic environment predictions. The framework combines the existing highly-scalable NASA production CFD code, Loci/CHEM, with a high-order accurate discontinuous Galerkin (DG) solver, Loci/THRUST, developed in the same computational framework. Loci/THRUST employs a low dissipation, high-order, unstructured DG method to accurately propagate acoustic waves away from the source regions across large distances. The DG solver is currently capable of solving up to 4th order solutions for non-linear, conservative acoustic field propagation. Higher order boundary conditions are implemented to accurately model the reflection and refraction of acoustic waves on launch pad components. The DG solver accepts generalized unstructured meshes, enabling efficient application of common mesh generation tools for CHEM and THRUST simulations. The DG solution is coupled with the CFD solution at interface boundaries placed near the CFD acoustic source regions. Both simulations are executed simultaneously with coordinated boundary condition data exchange.

Liever, Peter↗

Boundary Condition Study for the Juncture Flow Experiment in the NASA Langley 14x22-Foot Subsonic Wind Tunnel

Because future wind tunnel tests associated with the NASA Juncture Flow project are being designed for the purpose of CFD validation, considerable effort is going into the characterization of the wind tunnel boundary conditions, particularly at inflow. This is important not only because wind tunnel flowfield nonuniformities can play a role in integrated testing uncertainties, but also because the better the boundary conditions are known, the better CFD can accurately represent the experiment. This paper describes recent investigative wind tunnel tests involving two methods to measure and characterize the oncoming flow in the NASA Langley 14- by 22-Foot Subsonic Tunnel. The features of each method, as well as some of their pros and cons, are highlighted. Boundary conditions and modeling tactics currently used by CFD for empty-tunnel simulations are also described, and some results using three different CFD codes are shown. Preliminary CFD parametric studies associated with the Juncture Flow model are summarized, to determine sensitivities of the flow near the wing-body juncture region of the model to a variety of modeling decisions.

Rumsey, C. L.↗

TPSAS-NF1676L-12494-DND

Outline: NASA Missions and their requirements for Computational Fluid Dynamics (CFD) methodology, CFD state-of-the-art employed in NASA applications, examples of applications challenging NASA’s CFD capability, desired characteristics of future CFD solvers and the role of high-order methods in their development, and a vision for future production CFD methods.

David M Schuster↗

Insights and Lessons Learned from the NASA Juncture Flow Experiment

The NASA Juncture Flow experiment involved both CFD and wind tunnel measurements in its quest to provide CFD validation data for separated flow in a wing-fuselage corner. The experience has produced not only a wealth of valuable validation data and a new version of a turbulence model, it also yielded many lessons learned. This paper conveys those insights, particularly with respect to the qualities we believe to be essential in a CFD validation experiment. These include wind tunnel characterization and use of CFD as an assessment tool during the validation process. With a considerable number of validation tests already run both by the NASA team as well as by independent groups, a brief assessment is made of CFD’s current ability to predict the corner flow separation.

separation↗

Insights and Lessons Learned from the NASA Juncture Flow Experiment

The NASA Juncture Flow experiment involved both CFD and wind tunnel measurements in its quest to provide CFD validation data for separated flow in a wing-fuselage corner. The experience has produced not only a wealth of valuable validation data and a new version of a turbulence model, it also yielded many lessons learned. This paper conveys those insights, particularly with respect to the qualities we believe to be essential in a CFD validation experiment. These include wind tunnel characterization and use of CFD as an assessment tool during the validation process. With a considerable number of validation tests already run both by the NASA team as well as by independent groups, a brief assessment is made of CFD’s current ability to predict the corner flow separation.

separation↗

Fully Coupled Aeroelastic Stability Analysis of Adaptive Shape Memory Alloy Structural Technologies for Airframe Noise Reduction

This final report documents work performed by ATA Engineering, Inc., (ATA) to develop computational models and analyze the coupled fluid-structure response of two types of noise treatments applied to the leading-edge-slat component of a high-lift system typical of modern transport aircraft. The first treatment is a slat-gap filler (SGF), which closes the gap between the suction surfaces of a deployed slat and an aircraft main wing, and the second treatment is a slat-cove filler (SCF), which replaces the recirculating flow on the slat cove with a surface that promotes flow attachment. The representative airframe chosen for this work was NASA’s High-Lift Common Research Model (CRM-HL) in a baseline high-lift configuration. Superelastic shape memory alloys (SMAs) have been identified as enabling materials for these structural treatments. Since the technology elements rely upon having a highly reconfigurable structure, designs were assessed for their static aeroelastic deflection as well as their dynamic aeroelastic stability using coupled computational fluid dynamics (CFD) and nonlinear computational structural dynamics (NL-CSD) tools. Specifically, fluid-structure interaction (FSI) problems were solved computationally using the CFD solver Loci/CHEM and the NL-CSD solver Abaqus. As a part of the overall project, a similar capability was implemented using the CFD solver FUN3D coupled to Abaqus, although that work is documented in a separate report and that FSI framework was not used to analyze any of NASA’s SGF and SCF configurations. The technical approach consisted of solving for the flow field around the entire vehicle using a global CFD model, followed by extraction of relevant local subdomain data for CFD and NL-CSD co-simulations. The SGF design was analyzed using both 2D and 3D co-simulations to predict quasi-static aeroelastic deformations and to assess dynamic aeroelastic stability, whereas the SCF was analyzed in 2D only. SGF static aeroelastic response predictions focused on characterizing the deformed shape, with maximum displacements predicted to be on the order of magnitude of the technology element panel thickness. SGF dynamic aeroelastic response predictions used Partial Floquet analysis of the temporal evolution of selected nodal displacements to quantify the sign and magnitude of aeroelastic damping. Results suggest that the CRM-HL operating conditions would result in a dynamically stable response. The simulated dynamic pressure was also increased up to a factor of about four, and resulting responses suggest that predicted dynamic stability would be achieved with some margin.

Fluid Structure Interaction↗

Predicting Orion Launch Abort Acoustics

Computational fluid dynamics (CFD) scale-resolving simulations were performed to help assess the Orion launch abort system’s vibro-acoustic environment in collaboration with the Orion Loads and Dynamics team at the Johnson Space Center (JSC). The primary objectives are to establish best practices and to validate CFD predictions for surface fluctuating pressure spectra with ground, wind tunnel, and flight test data. Excellent agreement is obtained between CFD predictions by the Cartesian adaptive mesh refinement (AMR) module of the Launch, Ascent, and Vehicle Aerodynamics (LAVA) solver and measurements from a ground test representative of the start of a pad abort scenario, and with measurements from a wind tunnel test case corresponding to a transonic, high total angle of attack ascent abort scenario. A simplified (static and at rest) simulation of a pad abort flight test shows good agreement with the flight data, whereas a more costly simulation where the vehicle is accelerating relative to the ground from rest to Mach 0.4 shows consistent over-prediction of the spectra. Simulation results for a supersonic ascent abort flight test are in reasonable agreement with the measurements, but show larger excursions in regions where non-linear interactions between the plume turbulence, shocks, and turbulent boundary layer are strongest – likely due to CFD modeling assumptions, the effect of which are discussed. Overall, the CFD results demonstrate that scale-resolving simulation methodologies like the one presented can help characterize full-scale launch abort vehicle near-field noise and can complement wind tunnel, ground and flight tests by providing trends and reducing uncertainty for scenarios that cannot or are not planned to be tested.

Orion↗

Validation of Cryogenic Propellant Tank Self-Pressurization

NASA’s near-future long-term space missions necessitate advancements in cryogenic fluid management (CFM), which includes safe and reliable long-term propellant storage. Consequently, NASA STMD (Space Technology Mission Directorate) has established the CFM Portfolio Project to improve CFM technologies for upcoming missions. As a part of CFM Modeling Portfolio, NASA Marshall Space Flight Center’s (MSFC) Fluid Dynamics branch within the Propulsion Systems Department is tasked with assessing and improving computational tools used to support flight projects such as Human Lander System and Commercial Lunar Payload Services. One of the challenging modeling problems is that of self-pressurization of propellant tank due to heat leakage over long time periods. Reduced order and nodal tools find it extremely difficult to accurately predict self-pressurization under transient conditions or where complex flow patterns or thermal gradients exist, and application of 3-D CFD (computational fluid dynamics) simulations is necessary to characterize these problems. Until recently, CFD simulations for these long-term processes (order of hours or days) have been too impractical to conduct due to prohibitive wall time and computational resource requirements. The requisite CFD tool need to be efficient, computationally scalable, modular with ability to incorporate various physics models, and robust enough to not accumulate conservation errors over several hours of simulated time. NASA MSFC’s Loci-Stream CFD tool along with the VOF module is a great candidate to fit this mold. In this paper, we validate Loci-Stream for predicting self-pressurization of a flight scale propellant tank so it can serve as a reliable design and analysis tool for NASA’s CFM application needs. Liquid hydrogen tank pressurization tests carried out at the K-site testing facility provide a reliable data set for this purpose. These tests were simulated using Loci-Stream solver with VOF module as well as a hybrid approach which uses a lumped model for the ullage gas domain and CFD simulation of the liquid propellant. Both are shown to have very good predictive capabilities over multiple K-site experiments.

Hong Q. Yang↗

Validation of Cryogenic Propellant Tank Self-Pressurization

NASA's near-future long-term space missions necessitate advancements in cryogenic fluid management (CFM), which includes safe and reliable long-term propellant storage. Consequently, NASA STMD (Space Technology Mission Directorate) has established the CFM Portfolio Project to improve CFM technologies for upcoming missions. As a part of CFM Modeling Portfolio, NASA Marshall Space Flight Center's (MSFC) Fluid Dynamics branch within the Propulsion Systems Department is tasked with assessing and improving computational tools used to support flight projects such as Human Lander System and Commercial Lunar Payload Services. One of the challenging modeling problems is that of self-pressurization of propellant tank due to heat leakage over long time periods. Reduced order and nodal tools find it extremely difficult to accurately predict self-pressurization under transient conditions or where complex flow patterns or thermal gradients exist, and application of 3-D CFD (computational fluid dynamics) simulations is necessary to characterize these problems. Until recently, CFD simulations for these long-term processes (order of hours or days) have been too impractical to conduct due to prohibitive wall time and computational resource requirements. The requisite CFD tool need to be efficient, computationally scalable, modular with ability to incorporate various physics models, and robust enough to not accumulate conservation errors over several hours of simulated time. NASA MSFC's Loci-Stream CFD tool along with the VOF module is a great candidate to fit this mold. In this paper, we validate Loci-Stream for predicting self-pressurization of a flight scale propellant tank so it can serve as a reliable design and analysis tool for NASA's CFM application needs. Liquid hydrogen tank pressurization tests carried out at the K-site testing facility provide a reliable data set for this purpose. These tests were simulated using Loci-Stream solver with VOF module as well as a hybrid approach which uses a lumped model for the ullage gas domain and CFD simulation of the liquid propellant. Both are shown to have very good predictive capabilities over multiple K-site experiments.

cryogenic fluid management↗

Recommended DSMC Collision Model Parameters for Planetary Entry

Hybrid simulation frameworks combining Computational Fluid Dynamics (CFD) and Direct Simulation Monte Carlo (DSMC) are frequently employed to efficiently perform high-fidelity simulations of environments consisting of both continuum and rarified flow. DSMC is a stochastic, particle-based method which solves the fundamental Boltzmann equation and is therefore necessary for high-Knudsen flow where continuum-based assumptions governing CFD break down. However, the DSMC methodology is generally computationally inefficient to model the continuum regime. In a CFD/DSMC hybrid approach, obtaining an accurate, high-fidelity solution hinges on the consistent treatment of transport properties and the thermo-chemical models employed within the two solvers. In principle, in regions where CFD and DSMC are both employed, the same gas mixture under the same conditions should have the same properties, regardless of simulation type. Observed differences should be due to non-equilibrium processes, rather than differences in physical models. The goal of this work is to provide a comprehensive DSMC transport database encompassing the vast majority of known gas species encountered during Earth or other planetary atmospheric entry. This goal is largely possible due to recently performed ab initio quantum chemistry calculations. Combined with other high-fidelity data, the planned database will consist of collision integral data for over 200 neutral and ionized species and over 20000 binary collisions. From these collision integrals, Nelder-Mead optimization is used to compute collision-specific Variable Soft Sphere (VSS) collision model parameters, fitted from 300 K to 20000 K. Initial comparisons of transport properties of relevant equilibrium gas mixtures show great agreement between CFD and DSMC-derived results. The completed database can be readily applied to model binary collisions of any gas mixture containing the included species over the specified temperature range, making it a valuable tool for future planetary probe modeling efforts. An example is shown below in Fig. 1. Equilibrium mixture transport properties for a 35-species mixture composed originally of 10% air and 90% pyrolysis species of a carbon-phenolic ablator material [4] are computed using both fitted VSS Parameters and the original CFD collision integral values. Deviations in computed properties between the two approaches are less than 5% for the entire temperature range.

M. R. Gosma↗

The ZBOT-NC Experiment - Effects of Non-Condensable Gases on Propellant Tank Pressurization and Pressure Control

Introduction Integral to all phases of NASA’s projected planetary expeditions is affordable and reliable cryogenic fluid storage for use in propellant or life support systems. It is greatly advantageous to develop innovative vent-less pressure control designs based on cooling/mixing of the bulk tank fluid to allow storage of the cryogenic fluid with zero or reduced boil-off. The presence of noncondensable gases, can interfere with the condensation at the interface impacting tank pressure control during subcooled jet mixing especially in microgravity. The Zero-Boil-Off Tank (ZBOT) Experiments are a series of small-scale experiments aboard the International Space Station (ISS) that use a transparent volatile simulant fluid in a transparent sealed tank to delineate various fundamental fluid flow, heat and mass transport, and phase change phenomena associated with storage tank pressurization and pressure control in microgravity. The ZBOT-1 experiment was performed on the ISS 2017-2018 timeframe and collected data to validate a state-of-the-art CFD model for tank pressurization and pressure control for a pure system. The ZBOT-NC Experiment is the second experiment in the series to be performed on the International Space Station (ISS) in 2025. Its goal is to investigate the effects of noncondensable gases on interfacial evaporation and condensation during self -pressurization and jet-mixing pressure control in microgravity for a two-component system. Materials & Methods In this work, we will describe the detailed features of the ZBOT-NC experimental hardware and Diagnostics that includes nonintrusive Quantum Dot Thermometry (QDT) for whole field temperature measurement. All microgravity pressurization and pressure control tests will be performed for both the pure and the two-components systems with Xenon and Neon as the two noncondensable gases spanning small and large molecular weights and sizes. The two-phase CFD model that is developed as part of the project will be also presented and discussed. Results Ground-based pressurization and jet mixing experiments and CFD simulation results are compared to each other to validate both the fidelity of the CFD model predictions and the accuracy of the QDT measurement. Model simulations for noncondensable gas effects will also be compared against large Cryogenic LH2-GHe experiments to indicate the noncondensable gas effects on tank pressure control during jet mixing in 1G. Finally, CFD results will be presented to predict the effects of the noncondensable gas during subcooled jet mixing during the ZBOT-NC in advance of the microgravity experiment in 2015.

Evaporation Condensation↗

Computational Analysis of X-57 Rate Derivatives and Comparison to Wind Tunnel Data

Rate derivates of the wind tunnel model of X-57 were evaluated using hysteresis from computational fluid dynamics simulation and compared with that computed using the wind tunnel experimental data. To accomplish the task, hysteresis of the pitch, roll, and yaw oscillations were produced numerically using two different CFD solvers: STARCCM+ and Launch Ascent and Vehicle Analysis. The rate derivatives were computed using same parameter identification technique on CFD-based and wind tunnel-based hysteresis data. Results show CFD-based pitch rate derivatives have large error, up to 73%, while the lateral and direction derivates compared well with errors within 5% with respect to the that computed using the wind tunnel data. This may be due to lack of wind tunnel wall and sting model in the CFD simulation as well as lack of wall correction in the wind tunnel data which has larger impact on pitch than lateral and directional oscillations. Geometrical features such as the shape, slope of extremities, and enclosed area were compared between CFD-based and wind tunnel-based hysteresis. Results show that area enclosed by the hysteresis had largest influence on the accuracy of the computed rate derivatives given that shape of the hysteresis were roughly equivalent.

Seung Y. Yoo↗

The Zero-Boil-Off Tank Experiment: the Effects of Non-Condensable Gases on Pressurization and Pressure Control of Propellant Tanks in Microgravity

Introduction Integral to all phases of NASA’s projected planetary expeditions is affordable and reliable cryogenic fluid storage for use in propellant or life support systems. It is greatly advantageous to develop innovative vent-less pressure control designs based on cooling/mixing of the bulk tank fluid to allow storage of the cryogenic fluid with zero or reduced boil-off. The presence of noncondensable gases, can interfere with the condensation at the interface impacting tank pressure control during subcooled jet mixing especially in microgravity. The Zero-Boil-Off Tank (ZBOT) Experiments are a series of small-scale experiments aboard the International Space Station (ISS) that use a transparent volatile simulant fluid in a transparent sealed tank to delineate various fundamental fluid flow, heat and mass transport, and phase change phenomena associated with storage tank pressurization and pressure control in microgravity. The ZBOT-1 experiment was performed on the ISS 2017-2018 timeframe and collected data to validate a state-of-the-art CFD model for tank pressurization and pressure control for a pure system. The ZBOT-NC Experiment is the second experiment in the series to be performed on the International Space Station (ISS) in 2025. Its goal is to investigate the effects of noncondensable gases on interfacial evaporation and condensation during self -pressurization and jet-mixing pressure control in microgravity for a two-component system. Materials & Methods In this work, we will describe the detailed features of the ZBOT-NC experimental hardware and Diagnostics that includes nonintrusive Quantum Dot Thermometry (QDT) for whole field temperature measurement. All microgravity pressurization and pressure control tests will be performed for both the pure and the two-components systems with Xenon and Neon as the two noncondensable gases spanning small and large molecular weights and sizes. The two-phase CFD model that is developed as part of the project will be also presented and discussed. Results Ground-based pressurization and jet mixing experiments and CFD simulation results are compared to each other to validate both the fidelity of the CFD model predictions and the accuracy of the QDT measurement. Model simulations for noncondensable gas effects will also be compared against large Cryogenic LH2-GHe experiments to indicate the noncondensable gas effects on tank pressure control during jet mixing in 1G. Finally, CFD results will be presented to predict the effects of the noncondensable gas during subcooled jet mixing during the ZBOT-NC in advance of the microgravity experiment in 2015.

Evaporation Condensation↗

Modeling Self-Pressurization and Spray Bar Pressure Control of A Cryogenic Storage Tank in Normal Gravity

This paper presents computational fluid dynamics (CFD) models for simulating self-pressurization and spray-bar pressure control processes in a large-scale liquid hydrogen storage tank under normal gravity conditions. For self-pressurization, the model employs the kinetics-based Schrage equation alongside the volume-of-fluid (VOF) method to account for interfacial mass transfer. The CFD predictions of pressure and temperature are compared with experimental data from the Multipurpose Hydrogen Test Bed (MHTB) experiment, and the predicted interfacial mass transfer rates are also presented. A CFD model simulating pressure control using a spray bar has also been developed. An Eulerian-Lagrangian approach models the interactions between discrete droplets and the continuous ullage (vapor) phase. The spray model is coupled with the VOF method by tracking droplets in the ullage and removing them when they reach the liquid interface. The T-sat model calculates droplet-ullage heat and mass transfer, where droplets warm up to the saturation temperature corresponding to the ullage vapor pressure before evaporating while remaining at the saturation temperature. The evolution of tank pressure, vapor temperature, and liquid temperature predicted by the CFD model is validated against data from the MHTB spray-bar mixing experiment. Overall, the CFD models agree with experimental data, demonstrating their capability to simulate self-pressurization and pressure control processes in large-scale cryogenic storage tanks. These models can be valuable tools for designing and optimizing cryogenic fluid management systems in future applications.

Computational Fluid Dynamics↗

Modeling Self-Pressurization and Spray Bar Pressure Control of A Cryogenic Storage Tank in Normal Gravity

This paper presents computational fluid dynamics (CFD) models for simulating self-pressurization and spray-bar pressure control processes in a large-scale liquid hydrogen storage tank under normal gravity conditions. For self-pressurization, the model employs the kinetics-based Schrage equation alongside the volume-of-fluid (VOF) method to account for interfacial mass transfer. The CFD predictions of pressure and temperature are compared with experimental data from the Multipurpose Hydrogen Test Bed (MHTB) experiment, and the predicted interfacial mass transfer rates are also presented. A CFD model simulating pressure control using a spray bar has also been developed. An Eulerian-Lagrangian approach models the interactions between discrete droplets and the continuous ullage (vapor) phase. The spray model is coupled with the VOF method by tracking droplets in the ullage and removing them when they reach the liquid interface. The T-sat model calculates droplet-ullage heat and mass transfer, wherein droplets warm up to the saturation temperature corresponding to the ullage vapor pressure before evaporating while remaining at the saturation temperature. The evolution of tank pressure, vapor temperature, and liquid temperature predicted by the CFD model is validated against data from the MHTB spray-bar mixing experiment. Overall, the CFD models agree with experimental data, demonstrating their capability to simulate self-pressurization and pressure control processes in large-scale cryogenic storage tanks. These models can be valuable tools for designing and optimizing cryogenic fluid management systems in future applications.

Self-Pressurization↗

Development and Validation of a Slosh-Based Ullage Collapse Model

Ullage collapse modeling is an essential tool for the success of spacecraft powered by cryogenic propellants. The degree of ullage collapse can have a major impact on the design of cryogenic propellant tanks as well as the mission’s concept of operations. Additional commodities (e.g., helium gas for pressure control) and their quantities are dependent on the amount of heat transfer and phase change of propellants expected during flight. Currently, high-fidelity computational fluid dynamics (CFD) with heat transfer and phase change is the most accurate method for predicting ullage collapse. CFD simulations of this kind, however, can become very costly for problems requiring large domains, long simulation times, or significant liquid motion. Consequently, there is a need for a reduced order model to provide fast responses which can be used in engineering design decisions. The proposed numerical model is a nodal approach, augmented by outputs from CFD analyses that are generated rapidly and reliably. These outputs include transient values for liquid-vapor interface area, interface velocity, and liquid position using a volume of fluid (VoF) method. There are several benefits of combining a nodal approach with slosh predictions from CFD simulations: leveraging of the Loci/STREAM-VoF solver that has been validated in high-g and low-g environments, rapid modeling of VoF with isothermal propellant simplifications and relatively coarse meshes, and improved model predictions from the nodal tool. The nodal approach, informed by the slosh predictions from CFD, assumes uniform conditions in the liquid and gas domains. Modifications to account for temperature stratification in the liquid and gas may be implemented as a future improvement. Thermodynamic changes within the tank are solved at user-defined timesteps. First, impacts due to mass entering or leaving the control volume are captured. Second, heat transfer impacts are accounted for by tracking the volume of the gas which has been contacted by liquid in that timestep and rebalancing the internal energy. Finally, phase change is approximated using the area and velocity of the interface. The tank conditions are re-computed by conserving mass and energy. Calibration of the model to several test flights and experiments are required to baseline the predictive power and uncertainty. The final product will enable designers and engineers to create quick-turnaround predictions for ullage collapse and inform design feasibility.

Computational Fluid Dynamics↗

Reduced Order Model to Predict Dispersion of Flammable Refrigerant into a Space

As the HVAC&R industry mobilizes to deploy more low-GWP refrigerants, relevant standards are being continually reviewed and updated. Those include the general safety standards ISO 5149 and ASHRAE 15, and the equipment standards IEC and UL. The standards systematically set the allowable maximum amount of refrigerant that should be used in different equipment types and different applications. To do so, they rely on predictions of how a leaked refrigerant mass will disperse into a space. Dispersion characteristics, such as total flammable volume and its residence time, determine the risk associated with the presence of the flammable refrigerant. The standards have included provisions for the use of flammable refrigerants for approximately two decades. They relied on limited analytical analyses and test cases in their development. Dispersion of a refrigerant into a space is complex. Computational fluid dynamics (CFD) are the most accurate in predicting a given problem. However, CFD is computationally expensive and requires specialized expertise and resources and is not suitable for use by standards development working group as prediction tool. This paper presents the development of a reduced order model (ROM) that predicts the key dispersion characteristics relevant to the dispersion of a leaked refrigerant into a space for any combination of input variables. The inputs are the refrigerant release height, the total released refrigerant mass and its release flow rate, the refrigerant molecular weight, the ventilation flow rate, the floor area and height of the space, recirculation air flow rate, and the tightness of the space. The outputs are histograms of volume fraction of the room in prescribed concentration bins and the total mass of the refrigerant in each bin normalized by the total refrigerant charge at 13 prescribed simulation time stamps between 1 and 900 seconds. The ROM is constructed from a set of CFD simulations with carefully chosen combinations of input parameters. The selection if done using a multidimensional sparse grid which is a generalization of the classical tensor approach but offers additional flexibility and thus can be more carefully tuned towards a specific model. The tuning is done to improve the accuracy, measured in the difference between the output values of the ROM and the CFD model, while minimizing the computational cost, measured in number of CFD simulations which is orders of magnitude more expensive than the processing the training data.

Edwards, Dean↗

Implementation of a High-Fidelity Interface Resolving Method in Nek5000

The development and utilization of computational fluid dynamics (CFD) models for large, high-temperature electric melters in the Waste Treatment and Immobilization Plant (WTP) in eastern Washington State have proven to be valuable for various purposes. These models allow for a better understanding of the physio-chemical processes occurring within the melter vessels and can contribute to improving operational efficiency, throughput, and addressing operational issues related to vitrification. The CFD models employed for these melter vessels incorporate multiphase fluid flow and heat transfer simulations in different regions, including the plenum, cold cap, and molten glass regions. As the tank waste and glass formers are introduced into the melter, a reacting batch layer known as the cold cap forms on top of the molten glass. To enhance the melt rate, forced convection bubblers located at the bottom of the melters generate convection currents that help homogenize the molten glass and provide heat to the cold cap. As the bubbles rise through the highly viscous glass, they adopt a spherical-cap shape [1]. Meanwhile, the conversion of the batch to glass generates significant amounts of gases (such as water vapor, carbon dioxide, sulfur dioxide, and NOx) due to thermal decomposition [2]. These gases become trapped between the cold cap and molten glass, forming a foam layer [3]. For modeling multiphase flow in CFD and heat transfer simulations of waste glass melters at different scales, efforts are underway to augment the capabilities of the Nek5000 [4] and NekRS [5] open-source codes [6]. Nek5000/NekRS is a scalable and efficient spectral element code that has been successfully applied to a wide range of fluid dynamics problems. By leveraging the Nek5000/NekRS software, it becomes possible to model the melter systems more affordably and with lower computational requirements compared to currently utilized commercial CFD software. The specific objectives of this ongoing effort include: 1. Implementation of a level set method in Nek5000/NekRS: The level set method is a numerical technique commonly employed in CFD simulations to track and represent the interface between different phases or materials accurately. By incorporating this method into Nek5000/NekRS, the ability to simulate multiphase flows in waste glass melters at a high level of fidelity can be achieved. 2. Demonstration of capability for air bubbling through molten glass: As part of the development process, a specific case of air bubbling through molten glass will be simulated using the augmented Nek5000/NekRS code. This demonstration aims to showcase the ability of the software to accurately capture and analyze the complex phenomena involved in the multiphase flow within waste glass melters. By achieving these objectives, the improved Nek5000/NekRS code will offer a powerful computational tool for simulating and analyzing waste glass melter systems, enabling better understanding, optimization, and troubleshooting of these vitrification processes. The ability to accurately model and simulate multiphase flows has broad relevance across many industries and scientific domains, and the improved functionality can contribute to advancements in various fields beyond waste glass melter simulations.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗