Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “NATURAL CONVECTION”

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 397 records · Page 22

CFD Extraction of Heat Transfer Coefficient in Cryogenic Propellant Tanks

Current reduced-order thermal model for cryogenic propellant tanks is based on correlations built for flat plates collected in the 1950's. The use of these correlations suffers from inaccurate geometry representation; inaccurate gravity orientation; ambiguous length scale; and lack of detailed validation. This study uses first-principles based CFD methodology to compute heat transfer from the tank wall to the cryogenic fluids and extracts and correlates the equivalent heat transfer coefficient to support reduced-order thermal model. The CFD tool was first validated against available experimental data and commonly used correlations for natural convection along a vertically heated wall. Good agreements between the present prediction and experimental data have been found for flows in laminar as well turbulent regimes. The convective heat transfer between the tank wall and cryogenic propellant, and that between the tank wall and ullage gas were then simulated. The results showed that the commonly used heat transfer correlations for either vertical or horizontal plate over-predict heat transfer rate for the cryogenic tank, in some cases by as much as one order of magnitude. A characteristic length scale has been defined that can correlate all heat transfer coefficients for different fill levels into a single curve. This curve can be used for the reduced-order heat transfer model analysis.

Yang, H. Q.↗

Numerical and Experimental Study of Transport Phenomena in Directional Solidification of Succinonitrile

A numerical and experimental study of the growth of succinonitrile (SCN) using a horizontal Bridginan furnace and transparent glass ampoule was conducted. Two experiments were considered: one in which the temperature profile was fixed relative to the ampoule (no-growth case); and a second in which the thermal profile was translated at a constant rate (steady growth case). Measured temperature profiles on the outer surface of the ampoule were used as thermal boundary conditions for the modelling. The apparent heat capacity formulation combined with the variable viscositymeth~ was used to model the phase change in SeN. Both 2-D and 3-D models were studied and numerical solutions obtained using the commercial finite element code, FIDAP1. Comparison of the numerical results to experimental data showed excellent agreement. The complex 3-D shallow-cavity flow in the melt, differences between 2-D and 3-D models, effects of natural convection on the thermal gradient and shape of the solid/liquid interface, and the sensitivity of simulations to specific assumptions, are also discussed.

de Groh, Henry C., III↗

Numerical Investigation of Microgravity Tank Pressure Rise Due to Boiling

The ability to control self-pressurization in cryogenic storage tanks is essential for NASAs long-term space exploration missions. Predictions of the tank pressure rise in Space are needed in order to inform the microgravity design and optimization process. Due to the fact that natural convection is very weak in microgravity, heat leaks into the tank can create superheated regions in the liquid. The superheated regions can instigate microgravity boiling, giving rise to pressure spikes during self-pressurization. In this work, a CFD model is developed to predict the magnitude and duration of the microgravity pressure spikes. The model uses the Schrage equation to calculate the mass transfer, with a different accommodation coefficient for evaporation at the interface, condensation at the interface, and boiling in the bulk liquid. The implicit VOF model was used to account for the moving interface, with bounded second order time discretization. Validation of the models predictions was carried out using microgravity data from the Tank Pressure Control Experiment, which flew aboard the Space Shuttle Mission STS-52. Although this experiment was meant to study pressurization and pressure control, it underwent boiling during several tests. The pressure rise predicted by the CFD model compared well with the experimental data. The ZBOT microgravity experiment is scheduled to fly on February 2016 aboard the ISS. The CFD model was also used to perform simulations for setting parametric limits for the Zero-Boil-Off Tank (ZBOT) Experiments Test Matrix in an attempt to avoid boiling in the majority of the test runs that are aimed to study pressure increase rates during self-pressurization. *Supported in part by NASA ISS Physical Sciences Research Program, NASA HQ, USA

pressure rise↗

Multi-Node Modeling of Cryogenic Tank Pressurization System using Generalized Fluid System Simulation Program

This paper presents a multi-node model of autogenous pressurization of cryogenic propellant in a flight tank using the Generalized Fluid System Simulation Program (GFSSP), a general purpose flow network code developed at NASA/Marshall Space Flight Center. Tests were conducted to measure the pressure and temperatures at the various axial locations of the stratified ullage at 75% and 45% fill level. Liquid nitrogen was pressurized by gaseous nitrogen from a supply tank while the drain valve from the tank remained closed during the pressurization process. The ullage was discretized into 25 uniformly distributed nodes: 5 in the radial direction and 5 in the axial direction assuming the flow to be axisymmetric. Heat and mass transfer between the liquid and vapor has been modeled at the liquid vapor interface. Heat transfer between wall and vapor at the ullage has been accounted for by assuming heat transfer occurs by natural convection. The model also accounts for heat leak to the tank through the insulation and metal wall by heat conduction. The predicted pressures and temperatures are compared with the measured data.

Bandyopadhyay, Alak↗

Two-Stage Autoignition Dynamics of N-Dodecane Droplets Under Normal Gravity at High Pressures

Experimental observations of two-stage, autoignition characteristics of fiber-supported, n-dodecane fuel droplets under normal gravity are reported for a wide range of ambient temperatures, and pressures varying from 1 atm to 25 atm, which encompasses the fuel critical pressure (~ 18 atm). A high-speed camera is used to obtain backlit images of the autoignition process. During two-stage ignition, a cool-flame front forms below the droplet in the natural convective plume of the fuel vapor originating from the vaporizing droplet, and propagates upward, followed by the formation of a hot-flame kernel that expands rapidly and surrounds the droplet, eventually leading to the classical diffusion flame. Experimentally measured cool, and hot flame formation locations and propagation speeds are presented for varying initial temperatures and pressures. Measured first, and second stage ignition delay times are also reported. The influences of negative temperature coefficient (NTC) regime chemistry on flame starting location, propagation speeds, and ignition delay times are described. The first-stage cool flame characteristics are found to be controlled by both the chemistry, and physical aspects such as droplet evaporation and fuel-oxidizer mixing, while the second-stage hot flame ignition is primarily controlled by chemical kinetics. Further work is currently underway to develop models to quantitatively predict the observed behavior.

Rose, Evan↗

Marshall Space Flight Center Faculty Fellowship Program

This paper presents a multi-node model of autogenous pressurization of liquid nitrogen in a flight tank using the Generalized Fluid System Simulation Program (GFSSP), a general purpose flow network code developed at NASA/Marshall Space Flight Center. This study considers two different models: (a) pressurization of the tank by nitrogen gas from a supply tank and (b) pressurization of the tank by using the boil-off propellant through an IVF (Integrated Vehicle Fluids) loop. The tank pressurization model considers two different liquid fill-levels of the tank (a) 75% and (b) 45%. Heat and mass transfer between the liquid and vapor has been modeled at the liquid vapor interface. Heat transfer between wall and vapor at the ullage has been accounted for by assuming heat transfer occurs by natural convection. The model also accounts for heat leak to the tank through the insulation and metal wall by heat conduction. The predicted pressures and temperatures are compared with the measured data and good agreement has been observed. In the second case, the tank is pressurized by using the boil-off fluid that recirculates through an IVF loop. Preliminary results from this model are also reported.

N.F. Six↗

Numerical Modeling of No Vent Filling of a Cryogenic Tank

This paper presents a multi-node finite volume model of the chilling and filling of a cryogenic tank using the Generalized Fluid System Simulation Program (GFSSP), a general purpose flow network code. Vented Chill and No Vent Fill (VCNVF) tests were conducted at Marshall Space Flight Center where a flight-like tank was filled with liquid nitrogen from a supply tank. In the VCNVF tests, the tank was partially chilled while the vent valve was open. After the partial chilling, the vent valve was closed and the tank was filled without any venting. An integrated numerical model of the test set up was developed. The model included the transfer line from the supply tank, the target tank with the spray nozzle and solid walls, and the discharge line with the vent valve. The tank was discretized into multiple fluid nodes and branches to represent the ullage and liquid nitrogen and multiple solid nodes to represent the tank wall and structure. The heat transfer between solid to fluid was calculated from pool boiling correlations which include film, transition, and nucleate boiling, as well as natural convection during pre- and post-boiling. The model also accounts for the condensation of vapor in the tank when it comes in contact with the liquid spray. The predicted pressure, resident mass, wall and ullage temperature in the tank were compared with the test data.

Alok Majumdar↗

Numerical Modeling of No Vent Filling of a Cryogenic Tank

This paper presents a multi-node finite volume model of the chilling and filling of a cryogenic tank using the Generalized Fluid System Simulation Program (GFSSP), a general purpose flow network code. Vented Chill and No Vent Fill (VCNVF) tests were conducted at Marshall Space Flight Center where a flight-like tank was filled with liquid nitrogen from a supply tank. In the VCNVF tests, the tank was partially chilled while the vent valve was open. After the partial chilling, the vent valve was closed and the tank was filled without any venting. An integrated numerical model of the test set up was developed. The model included the transfer line from the supply tank, the target tank with the spray nozzle and solid walls, and the discharge line with the vent valve. The tank was discretized into multiple fluid nodes and branches to represent the ullage and liquid nitrogen and multiple solid nodes to represent the tank wall and structure. The heat transfer between solid to fluid was calculated from pool boiling correlations which include film, transition, and nucleate boiling, as well as natural convection during pre- and post-boiling. The model also accounts for the condensation of vapor in the tank when it comes in contact with the liquid spray. The predicted pressure, resident mass, wall and ullage temperature in the tank were compared with the test data.

Alok Majumdar↗

Assessment of Governing Heat and Mass Transfer Coefficients for Cryogenic No-Vent Top-Off Modeling

No-vent fill / no-vent top-off (NVF/NVTO) is a method to enable the trans-fer of cryogenic propellants without the need to vent. Transferring propellants inmicrogravity environments is necessary to support the long-duration space missionsplanned for the coming decades. A key developmental milestone to developing thiscapability fully is the simulation and experimentation of NVF/NVTO in 1-g settings.This work seeks to advance the former by assessing the effectiveness of different heattransfer correlations used in a 1-g non-equilibrium model. Heat transfer correlationsfor natural convection and two-phase boiling heat transfer are compared against ex-perimental results for 34 different tests. Two distinct experiments, comprised of 21tests, are presented in detail here. These tests cover a variety of different injectionmethods, receiver tank geometries, receiver tank materials, cryogenic propellants, andinitial tank states. Overall for the cases tested in this work, the model was able topredict the pressure response within the receiver tank to a mean absolute percentageerror of 25.9%. The temperature response error rate was 29.5% and 24.9% for wall-gasand two-phase nodes, respectively. Compared to the worst case correlation set tested,the pressure prediction error rate represents a 31% improvement. For the cases pre-sented in this work, the mean pressure prediction error rate was 15.5%. Future workcould evaluate the effectiveness of this model against experiments conducted withdifferent propellants or different initial conditions. In its current form, however, themodel can still be used to help reduce design and testing time for 1-g experiments,enabling the quicker iteration that is necessary to meet the technological demands offuture space missions.

No-Vent Top Off↗

Numerical Modeling of No Vent Filling of a Cryogenic Tank with Thermo-dynamic Vent System Assisted Injector

This paper presents a multi-node finite volume model of No Vent Filling (NVF) of a cryogenic tank with Thermo-dynamic Vent System (TVS) assisted injector, using the Generalized Fluid System Simulation Program (GFSSP), a general purpose flow network code. NVF tests were conducted in a CRYOgenic Orbital Testbed (CRYOTE) tank which was filled by liquid nitrogen using a TVS assisted injector with vent valve closed during the entire filling process. In TVS assisted injector, the liquid flow splits into two streams: one stream is routed to a Joule-Thomson (J-T) orifice where the flow immediately flashes from liquid to vapor or two phase mixture because the downstream to J-T leg is maintained at vacuum level.; the other stream is injected into the tank after being cooled by cold vapor of J-T leg in a heat exchanger. The flow through J-T leg is also used to cool the outer metal matrix of the injector which in turn cools the vapor in the tank ullage. The cooling of vapor in the ullage by cold injector surface reduces ullage pressure that allows liquid to enter and fill the tank. An integrated numerical model of the test set up was developed. The model included the dual lines of TVS assisted injector, tank and tank wall. The tank was discretized into multiple fluid nodes and branches to represent the ullage and liquid nitrogen and multiple solid nodes to represent the tank wall and structure. The heat transfer between solid to fluid was calculated from pool boiling correlations which include film, transition, and nucleate boiling, as well as natural convection during pre- and post-boiling. The model also accounts for the condensation of vapor at ullage-injector interface and when it comes in contact with the liquid spray. The predicted pressure, resident mass, wall temperature in the tank were compared with the test data.

Alok Majumdar↗

Validation of Gas-Liquid Sharp Interface Model in Loci-Stream for Propellant Tank Self-pressurization under Normal Gravity

Deep space missions require advances in cryogenic fluid management (CFM) for long term storage of propellants. One of the important phenomena is self-pressurization due to heat leakages into the tank. Managing self-pressurization is one of the key technologies for deep space exploration and long-term space missions. The complex interactions involving natural convection, thermal gradients, turbulence, and phase change near the gas-liquid interface cannot be modeled using reduced order or nodal analysis models, and 3-D CFD analyses are necessary to fully characterize the dynamics. CFD analyses pose their own difficulties. The requisite CFD tool to tackle this problem need to be modular with the ability to incorporate various physics models, efficient, and computationally scalable for simulating flight size tanks. In this paper, we discuss modeling approach for self-pressurization in situations where the liquid interface is static and validate NASA MSFC’s Loci-Stream CFD solver for predicting self-pressurization in a test tank in normal gravity. This tank was designed to capture quantifiable and accurate data for understanding various two-phase fluid phenomena and to validate modeling tools. Specifically, we demonstrate the capability of the Loci-Stream solver with a two-phase sharp-interface treatment to predict self-pressurization of a tank with an unperturbed gas-liquid interface. This validation highlights the reliability of our modeling approach and our solver to serve as a design and analysis tool for NASA’s CFM application needs.

cryogenic fluid management↗

Validation of Gas-Liquid Sharp Interface Model in Loci-Stream for Propellant Tank Self-pressurization under Normal Gravity

Deep space missions require advances in cryogenic fluid management (CFM) for long term storage of propellants. One of the important phenomena is self-pressurization due to heat leakages into the tank. Managing self-pressurization is one of the key technologies for deep space exploration and long-term space missions. The complex interactions involving natural convection, thermal gradients, turbulence, and phase change near the gas-liquid interface cannot be modeled using reduced order or nodal analysis models, and 3-D CFD analyses are necessary to fully characterize the dynamics. CFD analyses pose their own difficulties. The requisite CFD tool to tackle this problem need to be modular with the ability to incorporate various physics models, efficient, and computationally scalable for simulating flight size tanks. In this paper, we discuss modeling approach for self-pressurization in situations where the liquid interface is static and validate NASA MSFC's Loci-Stream CFD solver for predicting self-pressurization in a test tank in normal gravity. This tank was designed to capture quantifiable and accurate data for understanding various two-phase fluid phenomena and to validate modeling tools. Specifically, we demonstrate the capability of the Loci-Stream solver with a two-phase sharp-interface treatment to predict self-pressurization of a tank with an unperturbed gas-liquid interface. This validation highlights the reliability of our modeling approach and our solver to serve as a design and analysis tool for NASA's CFM application needs.

cryogenic fluid management↗

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 cooling processes in a large-scale liquid hydrogen storage tank under normal gravity. The self-pressurization model utilizes the kinetics-based Schrage equation with the Volume-Of-Fluid (VOF) model to account for the evaporative and condensing interfacial mass fluxes. It uses a laminar approach to model natural convection and compares the pressures and temperatures predicted by the CFD model with experimental data from the Multipurpose Hydrogen Test Bed (MHTB) self-pressurization experiment. Additionally, the predicted interfacial mass flux is presented. A CFD model for simulating pressure control using a spray bar is also presented. An Eulerian-Lagrangian approach models the interactions between the discrete droplets and the continuous ullage phase. The spray model is coupled with the VOF model by tracking particles in the ullage, removing particles from the ullage when they reach the interface, and then adding their contributions to the liquid. The T-sat model is used to calculate the droplet-ullage heat and mass transfer. In this model, droplets warm up to the saturation temperature corresponding to the ullage vapor pressure thereafter evaporating while remaining at the saturation temperature. The tank pressure, vapor, and liquid temperature evolutions are validated against data from the MHTB spray bar mixing experiment.

Self-Pressurization↗

Nodal Modeling of Tank Pressurization and Draining using a Multi-Node-Ullage Approach

The purpose of the pressurization system in liquid rocket propulsion is to control the pressure in the gas space of the propellant tank (known as the ullage space) and the propellant mass flowrate to the engine. A mathematical model is required to predict the amount of pressurant necessary to ensure that pressure and temperature levels inside the tank remain within acceptable limits and that the propellant pressure leaving the tank satisfies the net positive suction pressure (NPSP) requirement of the pump feeding the engine. Nodal codes typically model tank pressurization and draining using a single node to represent the ullage and a single node to represent the propellant. As the tank drains, the ullage node grows and the propellant node shrinks. The heat transfer between ullage to wall and ullage to propellant is governed by natural convection. Designers of liquid propulsion systems often use empirical correlations to estimate the “Collapse Factor” which represents the ratio of pressurant required with heat transfer and the amount of pressurant required without heat transfer. A single node ullage model of tank pressurization was developed using GFSSP to compute the collapse factor reasonably well and later was used to model tank pressurization during test firing of the FASTRAC rocket engine. The predicted tank pressure compared well with the test data. In the early 1970’s, pressurization and drain tests with liquid methane were performed at NASA Lewis Research Center in a vacuum chamber. A 5 ft diameter spherical aluminum tank was tested to drain from 95% to 5% full using gaseous helium, hydrogen, nitrogen and methane as pressurant. Tests were conducted with different pressurant inlet temperatures and drain times. Measured data include pressurant requirement, amount of pressurant condensed, and ullage and wall temperatures at various heights in the ullage space at the end of draining. A single node GFSSP model was developed to simulate helium pressurization of the methane tank. Predicted helium consumption was 8-23% less than measured. The average error of the six test cases was 16%. A single ullage node with multiple solid node model was developed using Thermal Desktop. Predicted helium consumption compares with the test data within 2%. This paper describes the development of a GFSSP multi-node ullage model of the test configuration and compares the predicted pressurant consumption for both helium and autogenous pressurization using gaseous methane with experimental as well as TD predictions.

Nodal Model↗

Simulations of Large-Scale Zero Boiloff, Densification, and Solidification of Hydrogen

An Integrated Refrigeration and Storage (IRAS) experimental system called the Ground Operations and Demonstration Unit for Liquid Hydrogen (GODU-LH2) demonstrated the ability store cryogens in a zero boiloff (ZBO) process, and to densify liquid hydrogen by reducing the temperature and pressure down to the triple point resulting in solidification. The incompressible Navier-Stokes was utilized to simulate ZBO, densification, and solidification of liquid hydrogen (LH2) in the IRAS tank. The simulations were performed using a commercially available Computational Fluid Dynamics (CFD) pressure-based mass and momentum flow model and an enthalpy-porosity energy model. Results demonstrated the simulation’s ability to predict time-dependent flow and temperature fields and solid-liquid phase locations for hydrogen during ZBO, densification, and solidification. The simulations showed good agreement with experimental data, with errors within acceptable ranges for temperature and pressure predictions, provided detailed insights into natural convection and solid-liquid phase change dynamics. These findings are critical for the design of future cryogenic fluid management systems.

Computational Fluid Dynamics↗

Thermal Insulation for a High Temperature Molten Salt Storage Tank in a CSP Plant: Preprint

A ternary molten chloride has been suggested as a high temperature (550-720 °C) sensible heat storage medium for a two-tank system. An effective thermal insulation is proposed to provide an opportunity for the fabrication of both the hot and cold storage tanks from lower cost metals, e.g. A240-347H, by maintaining the tank wall below 500 °C. To achieve this, the hot storage tank mainly comprises of inner shell (11.4 cm of a refractory brick as liner (KX-99 ®), 30 cm of insulation brick) and outer metal shell, where 5 cm of frozen layer of PCM800 (NaCl), and 5 cm of ceramic wool fill the annular space between the shells. Another layer of ceramic wool was introduced at the outer side of the metal shell to keep the temperature lower than 60 °C. ANSYS Fluent was used for the numerical study of the cool down behavior of the system during idle period when tank is quarter-filled. The heat transfer modes were considered as radiation and natural convection heat transfer from the hot molten salt at 720 °C, conduction through the solid layers and convection heat transfer from the outer shell to the environment at 15 °C. The results showed that the proposed insulation system including a frozen layer of PCM800 maintains the temperature at the metal wall lower than 500 °C, and lower than 60 °C at the outer side of the tank. Apart from the insulation property of the frozen layer of PCM800, it absorbs any leak of hot molten salt, acting as a barrier between salt and the metal shell. The thermal performance analysis provide an insight into the benefits of the multilayer insulation.

41 EE - Solar Energy Technologies Office (EE-4S)↗

A Study on Computational Cost Reduction of Simulations of Phase-Change Material (PCM) Embedded Heat Exchangers

Thermal storage can be implemented using Phase-Change Materials (PCM), which absorb significant latent heat with a relatively small temperature change. PCM phase-change processes are transient and are driven by thermal diffusion and natural convection – the latter, especially for melting process. Modeling and simulation of PCM heat exchangers (HX’s) is typically computationally intensive due to the relatively complex time-dependent physics. Most of the PCM modeling work in the literature uses high-order modeling tools such as Computational Fluid Dynamics (CFD) and Lattice-Boltzmann Method (LBM). For design purposes, the existing PCM modeling approaches are not practical, limiting researchers in their ability to investigate new ideas and different PCM’s with faster turnarounds. This paper presents a study investigating the reduction of computational cost of PCM embedded HX’s CFD models by evaluating the feasibility of spatial reduction without losing accuracy. The analysis consists of comparing full and partial domain under full melting conditions. The subject of this study is a single straight tube with circular transverse fins in the vertical orientation, using PCM’s with 35oC nominal melting temperature. Different tube and fin dimensions are investigated. Results indicate that the reduced domain reproduces - in half the run time - the same behavior as the full domain since the buoyancy effects are localized and patterned. The outputs from the partial domain simulation were used to build a non-general correlation for the PCM heat transfer characteristics and demonstrated how it can be implemented in a Finite Control Volume Reduced Order Model (ROM). The ROM can accurately reproduce the CFD simulations at 4 to 5 orders of magnitude faster.

25 ENERGY STORAGE↗

Automated Parameterized CFD Simulations of Phase-Change Material Embedded Heat Exchangers

Thermal storage can be implemented using Phase-Change Materials (PCM), which absorb significant latent heat with a relatively small temperature change. PCM phase-change processes are transient and are driven by thermal diffusion and natural convection – the latter, especially for melting process. Modeling and simulation of PCM heat exchangers is typically computationally intensive due to the relative complex time-dependent underlying physics. This paper presents a study investigating a single straight tube with circular transverse fins in the vertical orientation, using PCM’s with 37oC nominal melting temperature, leveraging automated parametric CFD simulations. The computational domain is reduced to 2D axisymmetric for a segment of the HX comprising a region in between two fins with constant wall temperature as boundary condition. This approach allows for faster evaluations of the transient thermal characteristics for a given design space of interest. Correlations and reduced order models (ROM’s) are developed. Results indicate that the ROM’s reproduce the CFD simulations with high accuracy for designs from the training data set, however, designs from cross-validation and test sets do not always have good agreement. The limited number of designs in the training set is a possible cause for the lack of robustness of the ROM’s and thus additional training designs will be added in future work. The ROM’s solve the simulations on average 500,000 times faster than CFD.

25 ENERGY STORAGE↗