Engineering Papers⌕ Search

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.

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

At least 217 records · Page 12

Cardiovascular and Body Fluid Adjustments During Bed Rest and Space Flight

Although a few scientific bed rest (BR) studies were conducted soon after World War II, advent of the space program provided impetus for utilizing prolonged (days-months) BR, which employed the horizontal or 6 degree head-down tilt (HDT) body positions, to simulate responses of healthy people to microgravity. Shorter (hours) HDT protocols were used to study initial mechanisms of the acclimation-deconditioning (reduction of physical fitness) syndromes. Of the major physiological factors modified during BR, reduced force on bones, ligaments, and muscles, and greatly reduced hydrostatic pressure within the cardiovascular system, the latter: which involves shifts of blood from the lower extremities into the upper body, increase in central venous pressure, and diuresis, appears to be the initial stimulus for acclimation. Increase in central venous pressure occurs in subjects during weightless parabolic flight, but not in astronauts early during orbital flight. But significant reduction in total body water (hypohydration) and plasma volume (hypovolemia) occurs in subjects during both BR and microgravity. Response of interstitial fluid volume is not as clear, It has been reported to increase during BR, and it may have increased in Skylab II and IV astronauts. Reduction of total body water, and greater proportional reduction of extracellular volume, indicates increased cellular volume which may contribute to inflight cephalic edema. Cerebral pressure abates after a few days of HDT, but not during flight. accompanied by normal (eugravity) blood constituent concentrations suggesting some degree of acclimation had occurred. But during reentry, with moderately increased +Gz (head-to-foot) acceleration and gravitational force, the microgravity "euhydration" becomes functional progressive dehydration contributing to the general reentry syndrome (GRS) which, upon landing the Shuttle, can and often results in gastrointestinal distress, disorientation, vertigo, fatigue, and fainting. Various pre-reentry hyperhydration procedures have been utilized to counteract the GRS. Thus, the somewhat decreased central venous pressure and lack of diuresis early in spaceflight suggests mechanisms other than the Gauer-Henry reflex are more important for maintaining fluid volume homeostasis in astronauts. Inflight hypohydration and hypovolemia are more important for maintaining fluid volume homeostasis in astronauts.

Greenleaf, John E.↗

Fluid and protein shifts after postural changes in humans

The effects of head-up, one-step, slow-graded, and rapid-graded tilt positions on blood density (BD), plasma density (PD), and erythrocyte density (ED) in male and female subjects are evaluated using a mass density detection method. The preparation of the samples and procedures for density and hematocrit measurements are described. The formula for calculating the density and volume of fluid shifted between intravascular and extravascular compartments is derived. The data reveal that ED does not change with varying tilt positions; however, PD and BD increase linearly with increasing angles of tilt, thereby resembling the time course of hemodilution and hemoconcentration.

Hinghofer-Szalkay, H.↗

Microburst modelling and scaling

A microburst can be modeled by releasing a volume of fluid that is slightly heavier than the ambient fluid, allowing it to fall onto a horizontal surface. Vorticity develops on the sides of this parcel as it descends and causes it to roll up into a turbulent vortex ring which impinges on the ground. Such a model exhibits many of the features of naturally occurring microbursts which are a hazard to aviation. In this paper this model is achieved experimentally by releasing a volume of salt water into fresh water from a cylindrical dispenser. When care is taken with the release the spreading rate of the surface outflow is measurable and quite repeatable despite the fact that the flow is turbulent. An elementary numerical approximation to this model, based on inviscid vortex dynamics, has also been developed. A scaling law is proposed which allows experiments with different fluid densities to be compared with each other and with the numerical results. More importantly the scaling law makes it possible to compare the model results with real microbursts.

Lundgren, T. S.↗

A Volume Flux Approach to Cryolava Dome Emplacement on Europa

We previously modeled a subset of domes on Europa with morphologies consistent with emplacement by viscous extrusions of cryolava. These models assumed instantaneous emplacement of a fixed volume of fluid onto the surface, followed by relaxation to form domes. However, this approach only allowed for the investigation of late-stage eruptive processes far from the vent and provided little insight into how cryolavas arrived at the surface. Consideration of dome emplacement as cryolavas erupt at the surface is therefore pertinent. A volume flux approach, in which lava erupts from the vent at a constant rate, was successfully applied to the formation of steep-sided volcanic domes on Venus. These domes are believed to have formed in the same manner as candi-date cryolava domes on Europa. In order to gain a more complete understanding of the potential for the emplacement of Europa domes via extrusive volcanism, we have applied this new volume flux approach to the formation of putative cryovolcanic domes on Europa. Assuming as in that europan cryolavas are briny, aqueous solutions which may or may not contain some ice crystal fraction, we present the results of this modeling and explore theories for the formation of low-albedo moats that surround some domes.

Cryolava Dome Emplacement↗

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↗

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↗

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↗

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↗

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↗

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↗

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.↗

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↗

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↗

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↗

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↗

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↗

Characterization of fluid physics effects on cardiovascular response to microgravity (G-572)

The recognition and understanding of cardiovascular adaptation to spaceflight has experienced substantial advancement in the last several years. In-flight echocardiographic measurements of astronaut cardiac function on the Space Shuttle have documented a 15 percent reduction in both left ventricular volume index and stroke volume with a compensatory increase in heart rate to maintain cardiac output. To date, the reduced cardiac size and stroke volume have been presumed to be the consequence of the reduction in circulating fluid volume following diuresis and other physiological processes to reduce blood volume within a few days after orbital insertion. However, no specific mechanism for the reduced stroke volume has been elucidated. The following investigation proposes the use of a hydraulic model of the cardiovascular system to examine the possibility that the observed reduction in stroke volume may, in part, be related to fluid physics effects on heart function. The automated model is being prepared to fly as a GAS payload. The experimental apparatus consists of a pneumatically actuated, elliptical artificial ventricle connected to a closed-loop, hydraulic circuit with compliance and resistance elements to create physiologic pressure and flow conditions. The ventricle is instrumented with high-fidelity, acceleration-insensitive, catheter-tip pressure transducers (Millar Instruments) in the apex and base to determine the instantaneous ventricular pressures and (delta)P(sub LV) across the left ventricle (LVP(sub apex)-LVP(sub base). The ventricle is also instrumented with a flow probe and pressure transducers immediately upstream of the inflow valve and downstream of the outflow valve. The experiment will be microprocessor controlled with analog signals stored on the FM data tape recorder. By varying the circulating fluid volume, ventricular function can be determined for varying preload pressures with fixed afterload pressure. Pilot experiments on board the NASA KC-135 aircraft have demonstrated proof-of-concept and provided early support for the proposed hypothesis. A review of the pilot experiments and developmental progress on the GAS version of this experiment will be presented.

Pantalos, George M.↗