SEARCH · Engineering Papers

Results for “Volume of fluid”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15

All sources NASA NTRS DOE OSTI

NASA NTRSDate not supplied

Validation of Cryogenic Propellant Tank Filling using Computational Fluid Dynamics Simulation

Validation of cryogenic propellant tank filling was performed using the Computational Fluid Dynamics (CFD) solver Loci/STREAM-VoF. The validation effort helped identify modeling methodologies that enable NASA to best support its partners in both launch pad and on-orbit filling operations. Data from liquid hydrogen ground tests filled via jet injection were used for the validation effort which include a sensitivity to initial tank wall temperature. Initially hot walls are expected to yield rapid evaporation and possibly boiling. A Volume of Fluid (VoF) methodology was used to capture the gas-liquid interface. Rapid breakup of the liquid jet was observed in simulation results as liquid evaporated and expanded. Inflowing liquid transitioned to a contiguous jet as tank temperatures decreased, gas pressure increased, and saturation conditions at the incoming liquid temperature were approached. The final phase of filling was distinguished by rapid gas pressure rise due to a higher rate of gas volume compression than condensation at the liquid surface. Key physics of propellant tank filling were captured in the computational predictions, and opportunities for added simulation robustness and efficiency in future modeling efforts were identified.

CFM↗

NASA NTRSDate not supplied

Validation of CFD Model Prediction of Flow Boiling Regime Transitions during LN2 Line Chilldown

Propellant storage and transfer during future long-duration missions will involve fuel depot operations in which a donor depot tank is used to fill a receiver spacecraft tank with a cryogenic propellant. Prior to the on-Orbit filling operation both the receiver tank and the transfer line must be cooled. The line chilldown process involves transition between boiling regimes in microgravity that will be quite different from their 1g ground-based counterparts. Since the cryogenic propellant itself will be used to perform the chilldown process, the time constants to cool the wall and the amount of fuel that will be used become important design considerations. In this light, the focus of the present work is to understand the flow boiling characteristics of a cryogenic fluid, namely, liquid nitrogen, during the chill-down of a transfer line using CFD modelling and simulations. The cryogenic chill-down process involves different flow boiling regimes: film boiling, transition boiling and nucleate boiling. The prediction of transition between these regimes in a CFD framework is new and challenging. The present work addresses this challenge by employing a volume-of-fluid (VOF) based methodology with Lee phase change model to predict the film boiling regime of the chill-down process in ANSYS Fluent®. The transition and nucleate boiling regimes are predicted by incorporating a sub-grid model developed by Craft Tech that accounts for bubble nucleation, growth, shedding frequency, and departure diameter. The sub-grid model is implemented into Fluent via a user-defined function for wall-fluid heat flux calculations. The sub-grid model is similar in formulation to the well-known Rensselaer Polytechnic Institute (RPI) boiling model. The model constants are tested for different operating conditions and designated values are reported. The CFD model is validated against published experimental data for liquid nitrogen chill-down of a heated stainless-steel pipe in 1g. Predicted results showing good agreement of wall temperature, rewetting temperature, and transition between film and nucleate boiling with the experimental measurements are presented and discussed.

Boiling regimes↗

NASA NTRSDate not supplied

Predictions of Line Chilldown Boiling Regime Transitions by a Coupled CFD-Sub-Grid Boiling Model Validated against 1G LN2 Experiments

Propellant storage and transfer during future long-duration missions will involve fuel depot operations in which a donor depot tank is used to fill a receiver spacecraft tank with a cryogenic propellant. Prior to the on-Orbit filling operation both the receiver tank and the transfer line must be cooled. The line chilldown process involves transition between boiling regimes in microgravity that will be quite different from their 1g ground-based counterparts. Since the cryogenic propellant itself will be used to perform the chilldown process, the time constants to cool the wall and the amount of fuel that will be used become important design considerations. In this light, the focus of the present work is to capture the flow boiling characteristics of a cryogenic fluid, namely, liquid nitrogen, during the chill-down of a transfer line using CFD modelling and simulations. The cryogenic chill-down process involves different flow boiling regimes: film boiling, transition boiling and nucleate boiling. The prediction of transition between these regimes in a CFD framework is new and challenging. The present work addresses this challenge by employing a volume-of-fluid (VOF) based methodology with Lee phase change model to predict the film boiling regime of the chill-down process in ANSYS Fluent®. The transition and nucleate boiling regimes are predicted by incorporating a sub-grid model that accounts for bubble nucleation, growth, shedding frequency, and departure diameter. The sub-grid model is implemented into Fluent via a user-defined function for wall-fluid heat flux calculations. The sub-grid model is similar in formulation to the well-known Rensselaer Polytechnic Institute (RPI) boiling model. The CFD model is validated against published experimental data for liquid nitrogen chill-down of a heated stainless-steel pipe in 1g. Predicted results shows good agreement of wall temperature, rewetting temperature, and transition between film and nucleate boiling with the experimental measurements for several different LN2 flowrates in the vertical orientation. The strengths and weaknesses of this modeling approach are presented and discussed.

Cryogenic Fluid Management↗

DOE OSTI2024

In vivo evaluation of an adaptive resuscitation controller using whole blood and crystalloid infusates for hemorrhagic shock

Introduction Hemorrhage remains the leading cause of preventable death on the battlefield. The most effective means to increase survivability is early hemorrhage control and fluid resuscitation. Unfortunately, fluid resuscitation requires constant adjustments to ensure casualty is properly managed, which is often not feasible in the pre-hospital setting. In this study, we showed how an adaptive closed-loop controller for hemorrhage resuscitation can be used to automate hemodynamic management using a swine hemorrhagic shock injury model. Methods The adaptive resuscitation controller (ARC) was previously developed to track pressure–volume responsiveness in real time and adjust its infusion rate to reach the target mean arterial pressure (MAP). Swine while maintained under a surgical plane of anesthesia and analgesia underwent a splenectomy, followed by two hemorrhage and resuscitation events. For the first resuscitation event, hemorrhage was induced to reduce the MAP to 35 mmHg until arterial lactate reached 4 mmol/L. The ARC system then infused whole blood (WB) to reach the target MAP and maintained the subject using crystalloids for 120 min. For the second resuscitation event, the subjects were hemorrhaged again but resuscitated using only crystalloid infusion to reach the target MAP and 120-min maintenance. Results The ARC was effective at WB resuscitation, reaching the target MAP in 2.0 ± 1.0 min. The median performance error was 1.1% ± 4.6%, and target overshoot was 14.4% ± 7.0% of the target MAP. The ARC maintained all animals throughout the 120 min maintenance period. For the second crystalloid-based resuscitation, ARC required a longer time to reach the target MAP, at an average rise time of 4.3 ± 4.0 min. However, target overshoot was reduced to 8.4% ± 7.3% of the target MAP. Much higher flow rates were required to maintain the target MAP during the second resuscitation event than during the first resuscitation event. Discussion The ARC was able to rapidly reach and maintain the target MAP effectively. However, this sometimes required large volumes of fluid as the ARC’s only goal was to reach the target MAP. Further clinical insight is needed regarding the preferred aggression level to achieve the target MAP. In conclusion, the ARC was successful in its programmed objective of reaching and maintaining the target MAP for extended periods of time in vivo , a critical next step toward improving hemorrhage treatment in the pre-hospital environment.

Snider, Eric J.↗

DOE OSTI2022

Ephemeral antibubbles: Spatiotemporal evolution from direct numerical simulations

Antibubbles, which consist of a shell of a low-density fluid inside a high-density fluid, have several promising applications. We show, via extensive direct numerical simulations (DNSs), in both two and three dimensions, that the spatiotemporal evolution of antibubbles can be described naturally by the coupled Cahn-Hilliard-Navier-Stokes (CHNS) equations for a binary fluid. Our DNSs capture elegantly the gravity-induced thinning and breakup of an antibubble via the time evolution of the Cahn-Hilliard scalar-order-parameter field $\phi$, which varies continuously across interfaces, so we do not have to enforce complicated boundary conditions at the moving antibubble interfaces. To ensure that our results are robust, we supplement our CHNS simulations with sharp-interface volume-of-fluid DNSs. We track the thickness of the antibubble and calculate the dependence of the lifetime of an antibubble on several parameters; we show that our DNS results agree with various experimental results; in particular, the velocity with which the arms of the antibubble retract after breakup scales as σ 1/2 , where σ is the surface tension.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRSDate not supplied

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↗

DOE OSTI2024

High- and Mid-Fidelity Modeling Comparison for a Floating Marine Turbine System

There is a lack of suitable numerical tools, particularly open-source tools, that can be used for designing and optimizing marine turbine systems. The National Renewable Energy Laboratory has added features to their widely used mid-fidelity wind turbine modeling code, OpenFAST, to enable modeling of axial-flow marine turbines. This necessitated the addition of several physical effects relevant to marine turbines that are neglected for wind turbines. These include buoyancy, added mass and inertial loads, wave-current superposition, and changes to the coordinate systems. This updated version of OpenFAST allows for the modeling of both fixed and floating marine turbines at a speed comparable to real time. While efficient for large sets of load cases and design studies, mid-fidelity codes make simplifying assumptions that may impact their accuracy. High-fidelity computational fluid dynamics (CFD) simulations can capture more flow effects with fewer assumptions and provide detailed body pressure mapping and flow-field information. It is important to compare predictions between mid-fidelity and high-fidelity codes, both to verify the models and to understand the limitations. A floating marine turbine system was modeled both with OpenFAST and with the commercial CFD code STAR-CCM+. The CFD model used a three-dimensional unsteady Reynolds-averaged Navier-Stokes solver for a volume-of-fluid numerical wave and current tank. The blade-resolved simulations used the sliding-interface technique for the spinning rotor and an overset grid to accommodate the rigid-body motion of the floating system. The mooring system was modeled with a custom coupling of the CFD solver with the open-source code MoorDyn. This improves upon the existing quasi-static catenary solver in STAR-CCM+, which lacks seabed contact or line-to-line connections. Simulation results for a floating marine turbine are compared between OpenFAST and CFD, highlighting the capabilities of the mid-fidelity code and identifying the areas where a high-fidelity approach is needed.

CFD↗

DOE OSTI2024

High- and Mid-Fidelity Modeling Comparison for a Floating Marine Turbine System: Preprint

There is a lack of suitable numerical tools, particularly open-source tools, that can be used for designing and optimizing marine turbine systems. The National Renewable Energy Laboratory (NREL) has added features to their widely used mid-fidelity wind turbine modeling code, OpenFAST, to enable modeling of axial flow marine turbines. This necessitated the addition of several physical effects relevant to marine turbines that can be neglected for wind turbines. These include buoyant loads, added mass and inertia loads, wave-current superposition, and changes to the coordinate systems. This updated version of OpenFAST allows for the modeling of both fixed and floating marine turbine systems at a speed comparable to real time. While efficient for long simulations, large sets of load cases, and design studies, mid-fidelity codes cannot capture all of the potentially important physical phenomenon impacting marine turbine systems. High-fidelity computational fluid dynamics (CFD) simulations can capture more flow effects with fewer assumptions and provide detailed body pressure mapping and flow-field information. It is important to compare predictions between mid-fidelity and high-fidelity codes, both to verify the models and to understand the limitations. A floating marine turbine system designed by NREL was modeled both with OpenFAST and with the commercial CFD code, STARCCM+. The CFD model used a 3-D unsteady Reynolds-averaged Navier-Stokes solver for a volume-of-fluid numerical wave and current tank. The blade-resolved simulations used the sliding-interface technique for the spinning rotor and an overset grid to accommodate the rigid-body motion of the floating system. The mooring system was modelled with a custom coupling of the CFD solver with the open-source code, MoorDyn. This improves upon the existing quasi-static catenary solver in STARCCM+, which lacks seabed contact or line-to-line connections. Spatial and temporal convergence studies were conducted. The simulation results for a combined current and wave condition are compared between OpenFAST and CFD, highlighting the capabilities of the mid-fidelity code and identifying the areas where a high-fidelity approach is needed.

CFD↗

NASA NTRS2021

Experimental and Computational Investigation of Flow Boiling in Microgravity

This study explores use of Computational Fluid Dynamics (CFD) to predict near-saturated flow boiling of FC-72 in microgravity. The computational method employs transient analysis to predict detailed interfacial behavior and heat transfer characteristics along a rectangular channel heated along two opposite walls. Predicted results are validated against experimental temperature measurements and high-speed video images captured during a series of parabolic aircraft maneuvers for three sets of operating conditions which include variations of both mass velocity and wall heat flux. The computational method is based on the multi-phase volume of fluid (VOF) model, which is combined with appropriate phase change and turbulence models, and accounts for both shear-lift force on bubbles and conjugate heat transfer along the heating walls. A key advantage of the CFD method is ability to capture details that are very difficult to measure experimentally, including detailed spatial variations of bubble shape, void fraction, mixture fluid temperature, liquid velocity, and vapor velocity, results for which are presented for each of the three test cases. Different flow regimes predicted along the heated length exhibit a number of dominant mechanisms including bubble nucleation, bubble growth, coalescence, vapor blankets, interfacial waviness, and residual liquid sub-layer, all of which agree well with experiment. Vapor velocity is shown to increase appreciably along the heated length because of increased void fraction, while liquid velocity experiences large fluctuations. Non-equilibrium effects are accentuated with increasing mass velocity, contributing minor deviations of fluid temperature from simulations compared to those predicted by the analytical method. Predicted wall temperature is fairly uniform in the middle of the heated length but increases in the entrance region, due to sensible heat transfer in the subcooled liquid, and decreases toward the exit, mostly because of flow acceleration resulting from increased void fraction.

two-phase cooling; flow boiling; CFD; microgravity↗

NASA NTRSDate not supplied

Validation of Cryogenic Propellant Tank Self-Pressurization by Leveraging Reduced Order Modeling within Computational Fluid Dynamics Simulation

Validation of cryogenic propellant tank self-pressurization was performed using a hybrid Computational Fluid Dynamics (CFD) and reduced order modeling methodology. Data from a liquid hydrogen ground test conducted at the K-site facility at the National Aeronautics and Space Administration (NASA) Glenn Research Center was used for the validation effort. Liquid phase dynamics were explicitly resolved with a CFD tool. Vapor phase dynamics were modeled as a point mass that communicated heat from the tank wall to the liquid phase via a boundary condition used at the gas-liquid interface. The method proved to be more accurate, robust, and efficient than explicit resolution of the dynamics using a standard Volume of Fluid (VOF) methodology. The subject pressurization process was found to be heavily dependent upon both the relatively high liquid temperature gradient near the gas-liquid interface and the natural convection flow path. Modeling the gas-liquid interface as an immovable surface eliminated temperature gradient destroying gas-liquid interface velocities observed in VOF simulations, and correspondingly enabled more rapid simulation since interface advection was not allowed. The single phase computational domain also facilitated the ability to demonstrate spatial resolution convergence of natural convection cells within the liquid which significantly impacted the tank pressurization rate. This work was used to demonstrate the critical physics for tank self-pressurization and numerical methodologies that may be used to best resolve those physics. The findings informed development and operation of production level CFD tools used in the Fluid Dynamics Branch at NASA Marshall Space Flight Center.

J. M. Brodnick↗

NASA NTRS2013

A Non-Cut Cell Immersed Boundary Method for Use in Icing Simulations

This paper describes a computational fluid dynamic method used for modelling changes in aircraft geometry due to icing. While an aircraft undergoes icing, the accumulated ice results in a geometric alteration of the aerodynamic surfaces. In computational simulations for icing, it is necessary that the corresponding geometric change is taken into consideration. The method used, herein, for the representation of the geometric change due to icing is a non-cut cell Immersed Boundary Method (IBM). Computational cells that are in a body fitted grid of a clean aerodynamic geometry that are inside a predicted ice formation are identified. An IBM is then used to change these cells from being active computational cells to having properties of viscous solid bodies. This method has been implemented in the NASA developed node centered, finite volume computational fluid dynamics code, FUN3D. The presented capability is tested for two-dimensional airfoils including a clean airfoil, an iced airfoil, and an airfoil in harmonic pitching motion about its quarter chord. For these simulations velocity contours, pressure distributions, coefficients of lift, coefficients of drag, and coefficients of pitching moment about the airfoil's quarter chord are computed and used for comparison against experimental results, a higher order panel method code with viscous effects, XFOIL, and the results from FUN3D's original solution process. The results of the IBM simulations show that the accuracy of the IBM compares satisfactorily with the experimental results, XFOIL results, and the results from FUN3D's original solution process.

Sarofeen, Christian M.↗

DOE OSTI2024

Liquid Piston with Spray Cooling Near-Isothermal Compressor

The goal of this project was to prototype and characterize the performance of a liquid-piston spray-cooled gas compressor. The working principle of the compressor enables optimized high-efficiency operation over a very wide range of operating conditions, unlike conventional compressors that are optimized for a narrow range of operating conditions. The compressor technology is suitable for many applications, such as gas pipeline transport, gas storage, and commercial and residential heat pumps. Both physical testing and computational fluid dynamics (CFD) modeling of the processes using the Oak Ridge National Laboratory high-performance computing center were completed. The experimental and CFD studies focused on a near-isothermal liquid-piston compressor (LPC) that uses propylene glycol to compress CO 2 . The first prototype demonstrated isothermal operation during several sequentially executed cycles of CO 2 compression and raised the temperature of the compressed CO 2 by only 2 K, compared with approximately 6 K when the gas was compressed non-isothermally. Isothermal operation was demonstrated at CO 2 flow rates of up to 2 L/min. The second prototype was designed with two compression chambers to allow continuous flow of high-pressure CO 2 . However, the design of the valve train to direct flow between the compression chambers was not sufficient to allow demonstration of CO 2 compression. Numerical simulations of the LPC in which the compression chamber was filled with propylene glycol injected from the bottom inlet were performed using large eddy simulation (LES) with the wall-adapting local eddy-viscosity subgrid-scale model coupled with the multiphase volume of fluid (VOF) model to simulate the transient interface between gas and liquid and to capture the heat and mass transfers within the compression chamber. In this effort, the effects of boundary conditions applied to the LES-VOF calculations (i.e., no wall, an adiabatic wall, and a wall with a heat flux subscribed) on the overall pressure and temperature of the CO 2 gas as well as the transient evolution of flow and heat transfer within the compression chamber were investigated and are discussed in this report. The LES calculation with no wall showed no dynamical flow patterns, and the volume-averaged temperature of CO 2 increased from 305 to 392.7 K, whereas LES calculations with a constant wall temperature or a wall heat flux had similar increases of CO 2 temperatures. The results of the LES simulation using a wall heat flux showed different stages in the compression process and revealed dynamical formation and interaction of CO 2 gas layers and circulation flow patterns within the chamber that contributed to the overall heat transfer between the solid wall, gas, and liquid surface in the compressor. Though an industrial partnership for commercializing the compressor was not secured, the technology was attractive for an industrial partner to use in two research proposals in response to US Department of Energy funding opportunity announcements.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

NASA NTRS1924

The Inertial Coefficients of an Airship in a Frictionless Fluid

This report deals with the investigation of the apparent inertia of an airship hull. The exact solution of the aerodynamical problem has been studied for hulls of various shapes and special attention has been given to the case of an ellipsoidal hull. In order that the results for this last case may be readily adapted to other cases, they are expressed in terms of the area and perimeter of the largest cross section perpendicular to the direction motion by means of a formula involving a coefficient K which varies only slowly when the shape of the hull is changed, being 0.637 for a circular or elliptic disk, 0.5 for a sphere, and about 0.25 for a spheroid of fineness ratio 7. For rough purposes it is sufficient to employ the coefficients, originally found for ellipsoids, for hulls otherwise shaped. When more exact values of the inertia are needed, estimates may be based on a study of the way in which K varies with different characteristics and for such a study the new coefficient possesses some advantage over one which is defined with reference to the volume of fluid displaced. The case of rotation of an airship hull has been investigated also and a coefficient has been defined with the same advantages as the corresponding coefficient for rectilinear motion.

Bateman, H↗

NASA NTRSDate not supplied

Surface Instability of Liquid Propellants in Microgravity During Pulsed Settling Operations

Pulsing reaction control system (RCS) thrusters, vent valves, or other propulsion devices can preserve propellant resources in space, but this operation also effectively introduces a vibration to the vehicle. When the vibration is perpendicular to the liquid propellant surface, Faraday waves may be generated at the liquid-vapor interface. These Faraday instabilities can perturb or break up the liquid surface of cryogenic tanks, leading to inefficiencies in thermal management or even ullage collapse. Drawing from theory and experiments, an engineering model defining the allowable design regions for pulsed settling in microgravity was assembled and verified with computational fluid dynamics (CFD) simulations. A traditional settling metric, the Bond number, was also overlaid in the model to indicate which duty cycles were insufficient to overcome surface tension and aggregate propellant. Mission planners and engineers can consult the tool to rapidly evaluate the stability of a liquid-vapor interface given the pulse frequency and the excitation acceleration. Expressions developed for Faraday waves induced by a sinusoidal forcing input at standard gravity were found to provide excellent predictive capabilities for pulsed, or rectangular, waveforms in the absence of a consistent gravitational acceleration. This study extends the usage of these equations to an alternative forcing function and microgravity environments for the purpose of estimating natural frequencies, surface mode shapes, surface wave amplitudes, and the onset of droplet ejection. CFD simulations with the Loci/STREAM-VoF (Volume of Fluid) solver were initially validated against experimental results in standard gravity. Discrete points on the design map were then investigated with CFD and confirmed that the engineering model reliably indicates surface stability and most Faraday wave characteristics without requiring higher-fidelity tools. The engineering model is highly extensible and can be adapted for various propellant fill fractions, fluid properties, and tank sizes.

Faraday Waves↗

NASA NTRSDate not supplied

Surface Instability of Liquid Propellants in Microgravity During Pulsed Settling Operations

Faraday Waves↗

DOE OSTI2025

Advanced Computational Modeling of High-Level Waste Vitrification at the Hanford Site

The U.S. Department of Energy (DOE) has selected vitrification for stabilizing legacy tank waste at the Hanford site, where radioactive waste from plutonium production was historically stored in underground tanks. This waste will be separated into low-activity waste (LAW) and high-level waste (HLW) fractions and processed at the Waste Treatment and Immobilization Plant (WTP). At WTP, glass melters are used for the vitrification of radioactive tank waste, transforming it into a stable borosilicate glass form for safe long-term storage. The melter vessel is constructed from highly durable and heat-resistant materials, where the vitrification process occurs. The main regions that are modeled are the melt pool, plenum, cold cap, riser/discharge chamber, and surrounding structure with insulation layers. Forced convection induced by air bubblers at the base of the melter ensure uniform temperature distribution and provide heat to the cold cap layer. The cold cap is a region of reacting batch feed that floats on top of the molten glass and is where the batch-to-glass reactions occur. Joule heating provided by electrodes mounted along the vertical walls of the melter and immersed directly in the glass, generates the necessary heat for the net endothermic conversion processes that occur in the cold cap. The high temperatures, radioactivity, and opaque nature of the glass prevent direct observation inside the melters. Therefore, computational models are essential for providing insight into factors that affect melter throughput. Thermocouples in the plenum provide operators with plenum temperature measurements. Operational adjustments include bubbling rate, voltage supplied to the electrodes, feed adjustments, and glass removal rate. Different computational fluid dynamics (CFD) models have been developed, each serving a specific purpose. There are CFD models of different scale melters, as well as models that capture the two-phase flow interfaces of rising bubbles in the molten glass or models with a simplified molten glass region so that the surrounding structure and plenum can be feasibly incorporated. Pilot-scale melter models have been developed to serve as validation of the methods employed in the simulation of the full-scale WTP melters. Models incorporating resolved bubbling are used to develop momentum source terms to implement into a single phase, multi-region, steady-state flow model that is being validated by measured process parameters such as glass production rate, voltage, input power, plenum temperatures, etc. The resolved bubbling model uses the multiphase volume of fluid approach to model the system with a high-resolution interface capturing scheme to maintain sharp interfaces between the molten glass and the air phase. The suite of CFD models is continually being improved to incorporate more realistic physics and achieve faster turnaround time. For example, an incremental controller is implemented to automatically adjust electrode voltage within the simulation to a molten glass set point temperature of 1150°C. Newer models feature improved meshes to ensure conformal meshes between regions and eliminate unnecessary mesh refinement in areas that are not of interest (such as boundary layers in offgas ports). Instead of explicitly modeling the structural, refractory, and insulation layers of the melter, a thermal resistance approach is used with published correlations used for boundary conditions. The development of robust and efficient CFD models will be instrumental in enabling the WTP to successfully fulfill its mission of safely stabilizing legacy nuclear waste.

12 - MGMT OF RADIOACTIVE AND NON-RADIOACTIVE WASTE↗

NASA NTRS2014

Self-Pressurization and Spray Cooling Simulations of the Multipurpose Hydrogen Test Bed (MHTB) Ground-Based Experiment

This paper presents a CFD (computational fluid dynamics) model for simulating the self-pressurization of a large scale liquid hydrogen storage tank. In this model, the kinetics-based Schrage equation is used to account for the evaporative and condensing interfacial mass flows. Laminar and turbulent approaches to modeling natural convection in the tank and heat and mass transfer at the interface are compared. The flow, temperature, and interfacial mass fluxes predicted by these two approaches during tank self-pressurization are compared against each other. The ullage pressure and vapor temperature evolutions are also compared against experimental data obtained from the MHTB (Multipuprpose Hydrogen Test Bed) self-pressurization experiment. A CFD model for cooling cryogenic storage tanks by spraying cold liquid in the ullage is also presented. The Euler- Lagrange approach is utilized for tracking the spray droplets and for modeling interaction between the droplets and the continuous phase (ullage). The spray model is coupled with the VOF (volume of fluid) model by performing particle tracking in the ullage, removing particles from the ullage when they reach the interface, and then adding their contributions to the liquid. Droplet ullage heat and mass transfer are modeled. The flow, temperature, and interfacial mass flux predicted by the model are presented. The ullage pressure is compared with experimental data obtained from the MHTB spray bar mixing experiment. The results of the models with only droplet/ullage heat transfer and with heat and mass transfer between the droplets and ullage are compared.

cryogenic fluid storage↗

NASA NTRS1992

A numerical model for atomization-spray coupling in liquid rocket thrust chambers

The physical process of atomization is an important consideration in the stable operation of liquid rocket engines. Many spray combustion computational fluid dynamics (CFD) codes do not include an atomization sub-model but assume arbitrary drop size distributions, drop initial locations, and velocities. A method of coupling an atomization model with the spray model in a REFLEQS CFD code is presented. This method is based on a jet-embedding technique in which the equations governing the liquid jet core are solved separately using the surrounding gas phase conditions. The droplet initial conditions are calculated using a stability analysis appropriate for the atomization regime of liquid jet break-up. This novel coupling model is used to analyze the SSME fuel preburner single injector flow. Results of the diffusion flame characteristics in a single injection element will be presented. The effect of relative velocity, mixture ratio and droplet initial conditions will be shown. The predictions of present atomization model is compared with that of the widely used CICM correlation. The results are also compared with the predictions of volume-of-fluid method.

M. G. Giridharan↗