Search NASA⌕ Search

SEARCH · Search NASA

Results for “two fluid model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 325 records · Page 18

Characterization of the Ignition Over-Pressure/Sound Suppression Water in the Space Launch System Mobile Launcher Using Volume of Fluid Modeling

The Space Launch System (SLS) Vehicle consists of a Core Stage with four RS-25 engines and two Solid Rocket Boosters (SRBs). This vehicle is launched from the Launchpad using a Mobile Launcher (ML) which supports the SLS vehicle until its liftoff from the ML under its own power. The combination of the four RS-25 engines and two SRBs generate a significant Ignition Over-Pressure (IOP) and Acoustic Sound environment. One of the mitigations of these environments is the Ignition Over-Pressure/Sound Suppression (IOP/SS) subsystem installed on the ML. This system consists of six water nozzles located parallel to and 24 inches downstream of each SRB nozzle exit plane as well as 16 water nozzles located parallel to and 53 inches downstream of the RS-25 nozzle exit plane. During launch of the SLS vehicle, water is ejected through each water nozzle to reduce the intensity of the transient pressure environment imposed upon the SLS vehicle. While required for the mitigation of the transient pressure environment on the SLS vehicle, the IOP/SS subsystem interacts (possibly adversely) with other systems located on the Launch Pad. One of the other systems that the IOP/SS water is anticipated to interact with is the Hydrogen Burn-Off Igniter System (HBOI). The HBOI system's purpose is to ignite the unburned hydrogen/air mixture that develops in and around the nozzle of the RS-25 engines during engine start. Due to the close proximity of the water system to the HBOI system, the presence of the IOP/SS may degrade the effectiveness of the HBOI system. Another system that the IOP/SS water may interact with adversely is the RS-25 engine nozzles and the SRB nozzles. The adverse interaction anticipated is the wetting, to a significant degree, of the RS-25 nozzles resulting in substantial weight of ice forming and water present to a significant degree upstream of the SRB nozzle exit plane inside the nozzle itself, posing significant additional blockage of the effluent that exits the nozzle upon motor start leading to detrimental effects. The purpose of the CFD simulations were to i) characterize the location of the IOP/SS water after it is ejected from the IOP/SS nozzles, ii) characterize the interaction of the IOP/SS system with the HBOI system and iii) characterize the interaction of the IOP/SS water with the RS-25 nozzles and the SRB nozzles.

West, Jeff↗

Numerical simulations of loops heated to solar flare temperatures. III - Asymmetrical heating

A numerical model is defined for asymmetric full solar flare loop heating and comparisons are made with observational data. The Dynamic Flux Tube Model is used to describe the heating process in terms of one-dimensional, two fluid conservation equations of mass, energy and momentum. An adaptive grid allows for the downward movement of the transition region caused by an advancing conduction front. A loop 20,000 km long is considered, along with a flare heating system and the hydrodynamic evolution of the loop. The model was applied to generating line profiles and spatial X-ray and UV line distributions, which were compared with SMM, P78-1 and Hintori data for Fe, Ca and Mg spectra. Little agreement was obtained, and it is suggested that flares be treated as multi-loop phenomena. Finally, it is concluded that chromospheric evaporation is not an effective mechanism for generating the soft X-ray bursts associated with flares.

Cheng, C.-C.↗

Zero-Boil-Off Tank (ZBOT) Experiment: Ground-Based Validation of Self-Pressurization and Pressure Control Two-Phase CFD Model

Integral to all phases of NASA's projected space and planetary expeditions is affordable and reliable cryogenic fluid storage for use in propellant or life support systems. Cryogen vaporization due to heat leaks into the tank from its surroundings and support structure can cause self-pressurization relieved through venting. This has led to a desire to develop innovative pressure control designs based on mixing of the bulk tank fluid together with some form of active or passive cooling to allow storage of the cryogenic fluid with zero or reduced boil-off. The Zero-Boil-Off Tank (ZBOT) Experiments are a series of small scale tank pressurization and pressure control experiments aboard the International Space Station (ISS) that use a transparent volatile simulant fluid in a transparent sealed tank to delineate various fundamental fluid flow, heat and mass transport, and phase change phenomena that control storage tank pressurization and pressure control in microgravity. The hardware for ZBOT-1 flew to ISS on the OA-7 flight in April 2017 and operations are planned to begin in September 2017, encompassing more than 90 tests. This paper presents preliminary results from ZBOT's ground-based research delineating both pressurization and pressure reduction trends in the sealed test tank. Tank self-pressurization tests are conducted under three modes: VJ heating, strip heating and simultaneous VJ and strip heating in attempt to simulate heat leaks from the environment, the support structure and both. The jet mixing pressure control studies are performed either from an elevated uniform temperature condition or from thermally stratified conditions following a self-pressurization run. Jet flow rates are varied from 2-25 cm/s spanning a range of jet Re number in laminar, transitional, and turbulent regimes and a range of Weber numbers covering no ullage penetration, partial penetration and complete ullage penetration and break-up (only in microgravity). Numerical prediction of a two-phase CFD model are compared to experimental 1g results to both validate the model and also indicate the effect of the residual non-condensable gas on evolution of pressure and temperature distributions in the tank during pressurization and pressure control.

microgravity evaporation condensatio↗

Numerical Simulation of No Vent Chill and Fill of a Large Liquid Hydrogen Tank

On orbit transfer of cryogenic liquids is necessary for many future space missions. The techniques required to transfer cryogens in micro gravity are quite different from ground system applications. In normal gravity, cryogenic tanks are usually filled from the bottom at nearly atmospheric pressure. The vapor generated due to heat transfer from warm tank walls, is allowed to vent from the top of the tank while the tank is being filled. If the normal gravity technique is used on orbit, the uncertainty of liquid and vapor distributions in unsettled condition may result in dumping of large amounts of liquid propellant overboard. The Vented Chill and No Vent Fill (VCNVF) process is a methodology [1] used to reduce propellant loss by holding the vent valve open during initial wall chilldown, preventing over-pressurization of the tank and stalling of the incoming flow, but closing the tank vent valve after chilling the tank wall sufficiently. An accurate simulation of no vent tank chill and fill is necessary to estimate the amount of propellant necessary to chill the tank, determine the maximum sustainable flowrate to fill the tank without stalling, and establish a timeline for the filling process. Numerical modeling of filling a cryogenic tank is complex because the flow is two-phase and the process involves both boiling and condensation with heat transfer between solid and fluid. In recent years, progress has been made to develop a numerical modeling technique [2,3,4] that has been verified by comparing with the test data from a small flight tank using liquid nitrogen [5]. The purpose of the present paper is to apply this technique in a larger tank with liquid hydrogen and compare with the test data of Chato [1]. The K-site Test Facility [1], shown in Figure 1, consists of a Test Tank, spray system, instrumentation, and the vacuum chamber. The test tank selected was ellipsoidal with an 87 inch major diameter and a 1.2 to 1 major to minor axis ratio. The two ends are joined by a short 1.5 inch cylindrical section. The tank is made of 2219 aluminum chemically milled to a nominal thickness of 0.087 inches. The tank weighs 329.25 pounds, and the tank volume is 175 ft 3 . The GFSSP model of the Test Tank is shown in Figure 2. Node 13 is a boundary node that represents the supply tank which is supplying liquid hydrogen at a given pressure and temperature. The total measured flowrate is evenly distributed through Branches 131 through 139. Nodes 1 through 9 represent the test tank and they are connected with metal solid nodes 14 through 22 through fluid-to-solid conductors that allow convective heat transfer between the solid and fluid nodes. The predicted tank pressure, filling rate, wall temperature and fluid temperature will be compared with test data.

Alok Majumdar↗

Observations on Cavitation Damage in a Flowing System

Cavitation damage to specimens of stainless steel, carbon steel, aluminum, and plexiglas, placed in a cavitating venturi using water and mercury as test fluids is mostly in the form of irregularly shaped pits which do not change with additional exposure to the cavitating field within the limited durations utilized. The rate of damage is very high initially, decreases for a relatively short period of time, then increases again up to the maximum test durations of 150 hours with water and 270 hours with mercury. Observation of damage effects by several independent techniques, using a variety of specimen materials, with two different fluids under various fluid dynamic conditions, leads to a suggested correlating model in terms of the cavitation bubble density and energy and specimen material strength.

VENTURI TUBE↗

Simulation of space plasma phenomena

The importance of computer modeling in investigating physical phenomena is discussed, noting that modern computers can follow the time evolution of systems containing many millions of degrees of freedom, all of which are simultaneously interacting with each other. The two types of models used in simulating plasma, that is, particle and fluid, are described. Attention is given to particle models and some recent results obtained concerning the diffuse aurora. The simulations are shown to give clear evidence that a type of nonlinear cyclotron resonance process for the cold particles is responsible for the saturation observed in these calculations. It is also responsible for the cold electron heating and it is found that at saturation the rate of heating of the cold electrons essentially balances the rate at which energy is fed into unstable waves.

Dawson, J. M.↗

A model to determine open or closed cellular convection

The horizontal asymmetry of cellular convection is demonstrated to be detectable by consideration of the vertical asymmetry of the driving force. A two-dimensional numerical model for convection in an internally heated and cooled fluid is presented, based on equations for the conservation of temperature and vorticity. Attention is focused on the steady-state finite-amplitude solutions for fixed Rayleigh number, Prandtl number, and aspect ratio. Temperature in the model corresponds to potential temperature under dry conditions and the equivalent potential temperature under saturated conditions. Slowly varying convection driven by asymmetric boundary fluxes is expressed in terms of steady-state convection driven by an asymmetric internal heat source. It is shown that heating near the ground produces open convection patterns where most of the fluid is descending, while cooling near the top of the flow leads to closed cellular patterns with a preponderance of ascending fluid. Extension of the model to three dimensions is indicated.

Helfand, H. M.↗

Rheological properties, shape oscillations, and coalescence of liquid drops with surfactants

A method was developed to deduce dynamic interfacial properties of liquid drops. The method involves measuring the frequency and damping of free quadrupole oscillations of an acoustically levitated drop. Experimental results from pure liquid-liquid systems agree well with theoretical predictions. Additionally, the effects of surfactants is considered. Extension of these results to a proposed microgravity experiment on the drop physics module (DPM) in USML-1 are discussed. Efforts are also underway to model the time history of the thickness of the fluid layer between two pre-coalescence drops, and to measure the film thickness experimentally. Preliminary results will be reported, along with plans for coalescence experiments proposed for USML-1.

Apfel, R. E.↗

Bubble trajectories in rotating drops

The investigation summarized in this paper describes recent normal gravity experiments involving the behavior of compound drops in rotating flows and in particular, the subsequent migration of the less dense phase (air bubble) toward the rotation axis. The data are compared to two models. The first was developed to predict the trajectory of a fluid particle in an infinitely large drop in the presence of both gravitational and rotational fields at the limit of quasi-steady creeping flow. The second predicts the trajectory of a fluid particle in a compound drop in the presence of a rotational field and at the limit of creeping flow. Gravity has not yet been incorporated into this second model.

Brone, D.↗

On the dynamics of a shock-bubble interaction

We present a detailed numerical study of the interaction of a weak shock wave with an isolated cylindrical gas inhomogenity. Such interactions have been studied experimentally in an attempt to elucidate the mechanisms whereby shock waves propagating through random media enhance mixing. Our study concentrates on the early phases of the interaction process which are dominated by repeated refractions of acoustic fronts at the bubble interface. Specifically, we have reproduced two of the experiments performed by Haas and Sturtevant : M(sub s) = 1.22 planar shock wave, moving through air, impinges on a cylindrical bubble which contains either helium or Refrigerant 22. These flows are modelled using the two-dimensional, compressible Euler equations for a two component fluid (air-helium or air-Refrigerant 22). Although simulations of shock wave phenomena are now fairly commonplace, they are mostly restricted to single component flows. Unfortunately, multi-component extensions of successful single component schemes often suffer from spurious oscillations which are generated at material interfaces. Here we avoid such problems by employing a novel, nonconservative shock-capturing scheme. In addition, we have utilized a sophisticated adaptive mesh refinement algorithm which enables extremely high resolution simulations to be performed relatively cheaply. Thus we have been able to reproduce numerically all the intricate mechanisms that were observed experimentally (e.g., transitions from regular to irregular refraction, cusp formation and shock wave focusing, multi-shock and Mach shock structures, jet formation, etc.), and we can now present an updated description for the dynamics of a shock-bubble interaction.

Quirk, James J.↗

Transport phenomena of reactive fluid flow in heterogeneous combustion processes.

A previously developed computer program was used to model two transient hybrid combustion processes involving tubes of solid Plexiglas. In the first study, representing combustion of a hybrid rocket, the oxidizing gas was oxygen, and calculations were continued sufficiently long to obtain steady-state values. Systematic variations were made in reaction rate constant, mass flow rate, and pressure, alternatively using constant and temperature dependent regression rate models for the fuel surface. Consistent results were obtained, as is evidenced by the values for the mass function of the reaction product and the flame temperature, for which plots are supplied. In the second study, fire initiation in a duct was studied, with an air mixture as the oxidizing gas. It was demonstrated that a satisfactory flame spread mechanism could be reproduced on the computer. In both of the above applications, the general, transient, two-dimensional conservation equations were represented, together with chemical reactions, solid-fuel interface conditions, and heat conduction in the solid fuel.

Hung, W. S. Y.↗

A Wall-Distance Method for Turbulence Modeling

The distance from a grid point to the closest wall surface, wall distance, is a funda- mental quantity in turbulence modeling. Efficiency of wall-distance calculations has become more critical as the size of computational grids has significantly increased in recent years. This paper reports on an initial implementation of a new search-based wall-distance method that is suitable for general unstructured computational fluid dynamics (CFD) grids and tailored for requirements specific for turbulence modeling. The method represents a two-step approach to calculate the wall distance. In the first step, the wall distance is approximated for each grid point as the minimum distance from this point to a vertex of a triangular face at the wall. The point-to-vertex distance calculation is relatively inexpensive but may lead to a significant error in the wall-distance ap- proximation, especially for grid points near the wall. In the second step, for grid points located within a predefined distance ( threshold ) from the wall, the wall distance is computed as the minimum distance to wall faces. As a result, the wall distance is exact for all grid points within the threshold. This two-step approach reduces the computational cost yet achieves high and controllable accuracy in the evaluation of the wall distance. Algorithmic enhancements are presented to improve efficiency of wall-distance computations. Comprehensive assessment of the new method is reported for large-scale unstructured CFD grids generated for the Fifth AIAA CFD High-Lift Prediction Workshop. The performance of the new wall-distance method compares favorably with performance of two established methods implemented in high-performance CFD codes.

Wall Distance↗

Computation of thermo-chemical nonequilibrium Martian atmospheric entry flows

The aerothermodynamics of a Martian atmospheric entry trajectory are studied. A chemical kinetics model for the CO2-N2 mixture found in the Martian atmosphere is implemented in a two-dimensional computational fluid dynamics method. The gas mixture is allowed to be in thermo-chemical nonequilibrium and to be composed of eight species characterized by two temperatures. The CO2 reaction model is validated by comparing computations to experimental results. Results for typical Martian entry conditions are presented and limitations of the current thermophysical model are discussed.

Candler, Graham V.↗

Progress of High Efficiency Centrifugal Compressor Simulations Using TURBO

Three-dimensional, time-accurate, and phase-lagged computational fluid dynamics (CFD) simulations of the High Efficiency Centrifugal Compressor (HECC) stage were generated using the TURBO solver. Changes to the TURBO Parallel Version 4 source code were made in order to properly model the no-slip boundary condition along the spinning hub region for centrifugal impellers. A startup procedure was developed to generate a converged flow field in TURBO. This procedure initialized computations on a coarsened mesh generated by the Turbomachinery Gridding System (TGS) and relied on a method of systematically increasing wheel speed and backpressure. Baseline design-speed TURBO results generally overpredicted total pressure ratio, adiabatic efficiency, and the choking flow rate of the HECC stage as compared with the design-intent CFD results of Code Leo. Including diffuser fillet geometry in the TURBO computation resulted in a 0.6 percent reduction in the choking flow rate and led to a better match with design-intent CFD. Diffuser fillets reduced annulus cross-sectional area but also reduced corner separation, and thus blockage, in the diffuser passage. It was found that the TURBO computations are somewhat insensitive to inlet total pressure changing from the TURBO default inlet pressure of 14.7 pounds per square inch (101.35 kilopascals) down to 11.0 pounds per square inch (75.83 kilopascals), the inlet pressure of the component test. Off-design tip clearance was modeled in TURBO in two computations: one in which the blade tip geometry was trimmed by 12 mils (0.3048 millimeters), and another in which the hub flow path was moved to reflect a 12-mil axial shift in the impeller hub, creating a step at the hub. The one-dimensional results of these two computations indicate non-negligible differences between the two modeling approaches.

turbomachinery↗

Multiphysics Simulation of Active Hypersonic Lip Cooling

This article describes the application of the Multidisciplinary Analysis (MDA) solver, Spectrum, in analyzing a hydrogen-cooled hypersonic cowl leading-edge structure. Spectrum, a multiphysics simulation code based on the finite element method, addresses compressible and incompressible fluid flow, structural, and thermal modeling, as well as the interactions between these disciplines. Fluid-solid-thermal interactions in a hydrogen impingement-cooled leading edge are predicted using Spectrum. Two- and semi-three-dimensional models are considered for a leading edge impingement coolant, concept under either specified external heat flux or aerothermodynamic heating from a Mach 5 external flow interaction. The solution accuracy is demonstrated from mesh refinement analysis. With active cooling, the leading edge surface temperature is drastically reduced from 1807 K of the adiabatic condition to 418 K. The internal coolant temperature profile exhibits a sharp gradient near channel/solid interface. Results from two different cooling channel configurations are also presented to illustrate the different behavior of alternative active cooling schemes.

Melis, Matthew E.↗

Ice-atmosphere interactions in the central Arctic: Remotely-sensed and simulated ice concentration and motion

The response of the Beaufort Sea ice pack to the passage of strong low pressure systems is studied using SAR (synthetic aperture radar), AVHRR (advanced very high resolution radiometer), and SSM/I (special sensor microwave/imager) data combined with model simulations. The SAR, AVHRR, and modeled concentrations concur generally in showing a 1 to 5 percent decrease in ice fraction during the passage of the lows. SSM/I derived concentrations appear to underestimate overall concentration and overestimate the opening of the pack relative to the other data sets. Ice motion derived from SAR and AVHRR and model simulations using two ice rheologies agree reasonably well in direction and magnitude, although the cavitating fluid approximation tends to overestimate ice motion compared to a viscous plastic rheology. Ice motions from SAR, AVHRR, and buoy observations can be merged to yield uniform grids particularly well suited to synoptic scale studies and model validation.

Maslanik, J. A.↗

Multiphase Flow Technology Impacts on Thermal Control Systems for Exploration

The Two-Phase Flow Facility (TPHIFFy) Project focused on bridging the critical knowledge gap by developing and demonstrating critical multiphase fluid products for advanced life support, thermal management and power conversion systems that are required to enable the Vision for Space Exploration. Safety and reliability of future systems will be enhanced by addressing critical microgravity fluid physics issues associated with flow boiling, condensation, phase separation, and system stability. The project included concept development, normal gravity testing, and reduced gravity aircraft flight campaigns, in preparation for the development of a space flight experiment implementation. Data will be utilized to develop predictive models that could be used for system design and operation. A single fluid, two-phase closed thermodynamic loop test bed was designed, assembled and tested. The major components in this test bed include: a boiler, a condenser, a phase separator and a circulating pump. The test loop was instrumented with flow meters, thermocouples, pressure transducers and both high speed and normal speed video cameras. A low boiling point surrogate fluid, FC-72, was selected based on scaling analyses using preliminary designs for operational systems. Preliminary results are presented which include flow regime transitions and some observations regarding system stability.

McQuillen, John↗

Physics Based Model for Online Fault Detection in Autonomous Cryogenic Loading System

We report the progress in the development of the chilldown model for rapid cryogenic loading system developed at KSC. The nontrivial characteristic feature of the analyzed chilldown regime is its active control by dump valves. The two-phase flow model of the chilldown is approximated as one-dimensional homogeneous fluid flow with no slip condition for the interphase velocity. The model is built using commercial SINDAFLUINT software. The results of numerical predictions are in good agreement with the experimental time traces. The obtained results pave the way to the application of the SINDAFLUINT model as a verification tool for the design and algorithm development required for autonomous loading operation.

Cryogenics↗