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.

SourceResearch connectionFrom yearThrough year

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

At least 469 records · Page 26

All sources NASA NTRS DOE OSTI

NASA NTRS1991

A Local Condensation Analysis Representing Two-phase Annular Flow in Condenser/radiator Capillary Tubes

NASA's effort for the thermal environmental control of the Space Station Freedom is directed towards the design, analysis, and development of an Active Thermal Control System (ATCS). A two phase, flow through condenser/radiator concept was baselined, as a part of the ATCS, for the radiation of space station thermal load into space. The proposed condenser rejects heat through direct condensation of ATCS working fluid (ammonia) in the small diameter radiator tubes. Analysis of the condensation process and design of condenser tubes are based on the available two phase flow models for the prediction of flow regimes, heat transfer, and pressure drops. The prediction formulas use the existing empirical relationships of friction factor at gas-liquid interface. An attempt is made to study the stability of interfacial waves in two phase annular flow. The formulation is presented of a stability problem in cylindrical coordinates. The contribution of fluid viscosity, surface tension, and transverse radius of curvature to the interfacial surface is included. A solution is obtained for Kelvin-Helmholtz instability problem which can be used to determine the critical and most dangerous wavelengths for interfacial waves.

Karimi, Amir↗

NASA NTRS2006

Multiphysics Thrust Chamber Modeling for Nuclear Thermal Propulsion

The objective of this effort is to develop an efficient and accurate thermo-fluid computational methodology to predict environments for a solid-core, nuclear thermal engine thrust chamber. The computational methodology is based on an unstructured-grid, pressure-based computational fluid dynamics formulation. A two-pronged approach is employed in this effort: A detailed thermo-fluid analysis on a multi-channel flow element for mid-section corrosion investigation; and a global modeling of the thrust chamber to understand the effect of heat transfer on thrust performance. Preliminary results on both aspects are presented.

Wang, Ten-See↗

NASA NTRSDate not supplied

Integrated Computational Fluid Dynamics and Thermal Desktop Thermal Modeling for Assessment of the EMU in Support of ISS EVA 80

Following the reports of water accumulation in the Extravehicular Mobility Unit (EMU) helmet during the ISS US EVA-80, various efforts for mitigation and further understanding of this phenomenon have been undertaken with the goal of prevention and to ensure crew safety in future EVAs. In support of this goal, a combination of Thermal Desktop thermal modeling and Computational Fluid Dynamics (CFD) has been performed in order to characterize the performance of the EMU components, specifically those within the ventilation/cooling water loops. This modeling effort evaluates the dry/wet gas pressure drop and flow distribution through gas and liquid flow paths using CFD, as well as a two-phase flow assessment of condensation production and water separation/removal performance using an integrated Thermal Desktop model of important EMU ventilation loop components. The models assess this overall thermal-fluid performance in comparison with previous existing models, design points, and historical EVA data using their inputs.

Blain Ryan Lancaster↗

NASA NTRSDate not supplied

Integrated Computational Fluid Dynamics and Thermal Desktop Thermal Modeling for Assessment of the EMU in Support of ISS EVA 80

EVA 80↗

NASA NTRS2001

Modeling of Chill Down in Cryogenic Transfer Lines

A numerical model to predict chill down in cryogenic transfer lines has been developed. Three chill down cases using hydrogen as the working fluid are solved: 1) a simplified model amenable to analytical solution, 2) a realistic model of superheated vapor flow, and 3) a realistic model of initially subcooled liquid flow. The first case compares a numerical model with an analytical solution with very good agreement between the two. Additionally, the analytical solution provides a convenient way to look at parametric effects on the chill down. The second and third cases are numerical models which provide temperature histories of the fluid and solid tube wall during chill down as well as several other quantities of interest such as pressure and mass flow rate. Of great interest is the ability to predict accurate values of chill down time (the time required to achieve steady-state cryogenic flow). The models predict that a 26 in. long, 3/16 in. ID aluminum tube has a shorter chill down time (approx. equal to 100 sec) and uses less hydrogen with superheated vapor flow than with initially subcooled liquid flow (greater than 200 sec for chill down).

Cross, Matthew F.↗

NASA NTRS1993

Computational fluid dynamic modelling of cavitation

Models in sheet cavitation in cryogenic fluids are developed for use in Euler and Navier-Stokes codes. The models are based upon earlier potential-flow models but enable the cavity inception point, length, and shape to be determined as part of the computation. In the present paper, numerical solutions are compared with experimental measurements for both pressure distribution and cavity length. Comparisons between models are also presented. The CFD model provides a relatively simple modification to an existing code to enable cavitation performance predictions to be included. The analysis also has the added ability of incorporating thermodynamic effects of cryogenic fluids into the analysis. Extensions of the current two-dimensional steady state analysis to three-dimensions and/or time-dependent flows are, in principle, straightforward although geometrical issues become more complicated. Linearized models, however offer promise of providing effective cavitation modeling in three-dimensions. This analysis presents good potential for improved understanding of many phenomena associated with cavity flows.

Deshpande, Manish↗

NASA NTRSDate not supplied

Modeling Electrolytic O2 Recovery from Metabolic CO2 for Advanced Closed Loop Life Support Systems in Extraterrestrial Human Missions

The International Space Station (ISS) is currently equipped with a complex, heavy, and power consuming system that recovers approximately 50% of O2 from metabolic CO2. Future long duration human missions to the Moon and Mars will necessitate a sustainable and highly efficient metabolic oxygen recovery system capable of yielding a minimum of 75% O2 recovery. A Macrofluidic Electrochemical Reactor (MFECR) technology development effort is currently underway at NASA Marshall Space Flight Center (MSFC) to significantly increase current metabolic O2 recovery efficiency, expand mission sustainability, and reduce complexity of the system. The novel design combines CO2 conversion to O2 along with C2H4 as byproduct and water electrolysis (currently conducted in two separate units) into a single compact unit that runs at standard conditions and is theoretically capable of generating O2 with a theoretical maximum metabolic CO2 conversion of 73% while consuming less than metabolic water. This paper presents a comprehensive multi-physic 3D model developed at MSFC on CO2 conversion to O2 and C2H4 at standard conditions via MFECR. The 3D spatial domain of the model is a replica of the actual MFECR’s 3D drawing generated for the MFECR fabrication and operated to recover O2 from CO2 yielding C2H4 as byproduct. Electrochemical (EC) physics that includes EC multicomponent reaction mechanisms, mass transport, and electrical current density distributions is coupled in the model with all the other physics phenomena involved in the MFECR’s process, such as two-phase flow, free and porous fluid regimes, multicomponent mass transfer, heat transfer, and DC electrical current generation along with Joule heating effect. The EC reaction sections of the MFECR consists of two porous gas diffusion electrodes (GDE) and an electrolyte serpentine channel sandwiched in the middle. The CO2 feeds the cathode serpentine channel and part of the O2 product is fed back to the anode serpentine chamber. An alkaline solution feeds the electrolyte serpentine chamber wetting the GDEs of both, the anode and cathode allowing the OH- ionic transport between them. The EC reactions in the cathode’s GDE yield C2H4 from CO2 and H2 from water while the EC reaction in the anode’s GDE yields O2. The authors will present in this paper the validation of the model using experimental data and the utilization of the validated model in building a reliable simulator that will not only assist the authors on the MFECR design but also the optimization of its operation in the ISS and future spatial human missions.

Jesus A Dominguez↗

NASA NTRSDate not supplied

Conceptual Design of a Test Section for a Cryogenic Transfer Line Experiment

This project presents the development of a comprehensive computer-aided design (CAD) model for the Transfer Line test section, which is an upcoming test module for the Flow Boiling and Condensation Experiment (FBCE) aboard the International Space Station. The FBCE aims to serve as the primary platform to study two-phase flow in microgravity, and the Transfer Line module aims to utilize its existing hardware to study cryogenic chilldown; a critical component of cryogenic propellant transfer in space. The CAD model was created using SolidWorks software and includes the test section’s main components: a heat transfer line, a flow visualization line, insulation, high-speed visible imaging, and infrared imaging. The heat transfer line primarily consists of a stainless steel tube that will be heated via Joule heating. The flow visualization line consists of segmented Indium-Titanium-Oxide coated sapphire tubes, which allow for visualization of flow regimes by high-speed and infrared cameras, similar thermal properties as the stainless steel tube, and Joule heating via the ITO coating resistance. Polyether Ether Ketone was used to insulate and encapsulate the two sections. High-speed visible imaging and infrared imaging were also incorporated to view the propagation of fluid during inverted annular flow and map out the thermal profile of the ITO-coated sapphire tube. This model serves as the guiding concept design for the Transfer Line project and will assist the project team during future safety reviews and concept design reviews. Overall, this work is significant in two-phase flow research and offers a promising path towards an established facility for cryogenic transfer line research in microgravity.

Flow boiling↗

NASA NTRS1996

Animal models and their importance to human physiological responses in microgravity

Two prominent theories to explain the physiological effects of microgravity relate to the cascade of changes associated with the cephalic shifts of fluids and the absence of tissue deformation forces. One-g experiments for humans used bed rest and the head-down tilt (HDT) method, while animal experiments have been conducted using the tail-suspended, head-down, and hindlimbs non-weightbearing model. Because of the success of the HDT approach with rats to simulate the gravitational effects on the musculoskeletal system exhibited by humans, the same model has been used to study the effects of gravity on the cardiopulmonary systems of humans and other vertebrates. Results to date indicate the model is effective in producing comparable changes associated with blood volume, erythropoiesis, cardiac mass, baroreceptor responsiveness, carbohydrate metabolism, post-flight VO2max, and post-flight cardiac output during exercise. Inherent with these results is the potential of the model to be useful in investigating responsible mechanisms. The suspension model has promise in understanding the capillary blood PO2 changes in space as well as the arterial PO2 changes in subjects participating in a HDT experiment. However, whether the model can provide insights on the up-or-down regulation of adrenoreceptors remains to be determined, and many investigators believe the HDT approach should not be followed to study gravitational influences on pulmonary function in either humans or animals. It was concluded that the tail-suspended animal model had sufficient merit to study in-flight and post-flight human physiological responses and mechanisms.

NASA Discipline Cardiopulmonary↗

NASA NTRS2016

Passive Rocket Diffuser Theory: A Re-Examination of Minimum Second Throat Size

Second-throat diffusers serve to isolate rocket engines from the effects of ambient back pressure during testing without using active control systems. Among the most critical design parameters is the relative area of the diffuser throat to that of the nozzle throat. A smaller second throat is generally desirable because it decreases the stagnation-to-ambient pressure ratio the diffuser requires for nominal operation. There is a limit, however. Below a certain size, the second throat can cause pressure buildup within the diffuser and prevent it from reaching the start condition that protects the nozzle from side-load damage. This paper presents a method for improved estimation of the minimum second throat area which enables diffuser start. The new 3-zone model uses traditional quasi-one-dimensional compressible flow theory to approximate the structure of two distinct diffuser flow fields observed in Computational Fluid Dynamics (CFD) simulations and combines them to provide a less-conservative estimate of the second throat size limit. It is unique among second throat sizing methods in that it accounts for all major conical nozzle and second throat diffuser design parameters within its limits of application. The performance of the 3-zone method is compared to the historical normal shock and force balance methods, and verified against a large number of CFD simulations at specific heat ratios of 1.4 and 1.25. Validation is left as future work, and the model is currently intended to function only as a first-order design tool.

Jones, Daniel R.↗

NASA NTRS1981

Optimal time splitting for two- and three-dimensional Navier-Stokes equations with mixed derivatives

A new explicit, time splitting algorithm has been developed for finite difference modelling of the full two and three-dimensional time-dependent, compressible, viscous Navier-Stokes equations of fluid mechanics. The scheme is optimal in the sense that the split operators achieve their maximum allowable time step, i.e., the corresponding Courant number. The algorithm allows a conservation-form formulation. Stability is proven analytically and verified numerically. In proving stability it was shown that all nine matrix coefficients of the Navier-Stokes equations are simultaneously symmetrizable by a similarity transformation. Two such transformations and their resulting symmetric matrix coefficients are presented explicitly.

Abarbanel, S.↗

NASA NTRS1981

An experimental investigation into the feasibility of a thermoelectric heat flux gage

An experiment was conducted to determine the feasibility of using a commerically available thermoelectric device as a heat flux gage. In certain research applications, the thermoelectric heat flux gage can provide a relatively simple means to model a warm fluid--cold wall. The experiment showed that heat flux through the gage could be correlated within two percent to the applied thermoelectric current through the device and the hot and cold side temperature with a simple algebraic equation.

Jones, J. C.↗

NASA NTRS1990

Evidence for the presence of quasi-two-dimensional nearly incompressible fluctuations in the solar wind

Assuming that the slab and isotropic models of solar wind turbulence need modification (largely due to the observed anisotropy of the interplanetary fluctuations and the results of laboratory plasma experiments), this paper proposes a model of the solar wind. The solar wind is seen as a fluid which contains both classical transverse Alfvenic fluctuations and a population of quasi-transverse fluctuations. In quasi-two-dimensional turbulence, the pitch angle scattering by resonant wave-particle interactions is suppressed, and the direction of minimum variance of interplanetary fluctuations is parallel to the mean magnetic field. The assumed incompressibility is consistent with the fact that the density fluctuations are small and anticorrelated, and that the total pressure at small scales is nearly constant.

Matthaeus, William H.↗

NASA NTRS2016

Thermal Ablation Modeling for Silicate Materials

A thermal ablation model for silicates is proposed. The model includes the mass losses through the balance between evaporation and condensation, and through the moving molten layer driven by surface shear force and pressure gradient. This model can be applied in ablation simulations of the meteoroid or glassy Thermal Protection Systems for spacecraft. Time-dependent axi-symmetric computations are performed by coupling the fluid dynamics code, Data-Parallel Line Relaxation program, with the material response code, Two-dimensional Implicit Thermal Ablation simulation program, to predict the mass lost rates and shape change. For model validation, the surface recession of fused amorphous quartz rod is computed, and the recession predictions reasonably agree with available data. The present parametric studies for two groups of meteoroid earth entry conditions indicate that the mass loss through moving molten layer is negligibly small for heat-flux conditions at around 1 MW/cm(exp. 2).

Ablators↗

NASA NTRS2002

Constitutive Relation in Transitional Granular Flows

To study the constitutive behavior of granular materials, the presence of gravity is detrimental. Although empirical relations have been obtained for engineering designs to control granular flows on Earth, it is not known how well these Earthbound relations can be used in another gravity field. Fundamental understanding must be derived to reliably design for granular flows in space exploration. There are two extremes of granular flows of which significant amount of knowledge is available. One deals with a dense and quasi-static situation where the deformation rate nearly vanishes. The other deals with dilute and rapidly fluctuating grain velocities where particle inertia dominates. This project, funded by the NASA Microgravity Fluid Physics Program, aims to study this transitional regime via physical experiments and computer simulations. A conceptual model has been established as described below. There are two natural time scales in a granular flow. One is the travel time between two consecutive collisions and the other is the duration of a collision contact. At a very low shear-rate, the shear-induced particle velocity is low. Hence the travel time between collisions is longer than the contact time between colliding particles. Binary collisions prevail. As the shear-rate increases, the traveling time between collisions reduces and the probability of multiple collisions goes up. These particle groups disperse shortly after and new groups form. When shear-rate is further increased, clusters grow in size due to an increasing chance for free particles to join before groups have the time to disperse. The maximum cluster size may depend on the global concentration and material properties. As the solid concentration approaches zero, the cluster size goes to one particle diameter. The maximum possible cluster size under any condition is the container size, provided that the shear flow is inside a container. The critical shear-rates that dictate the initiation of the multiple contacts, and the size and lifetime of the collision clusters, are functions of the concentration also. A 'regime' theory has been proposed by Babic et al. This theory suggested that both the solid concentration, C, and the non-dimensional shear-rate, B, are important in determining the regimes of the granular constitutive law. Additional information is included in the original extended abstract.

Shen, Hayley H.↗

NASA NTRS1979

A simplified model for the viscous crossflow in a slotted test section

A simplified physical model was constructed which simulates the viscous crossflow in a fluid layer near the slots at a fixed streamwise location in a slotted wind tunnel. For low to moderate Reynolds numbers, numerical solutions of the two dimensional, incompressible Navier-Stokes equations in stream function and vorticity, which govern the model flow, were obtained. Fairly general slot geometry was incorporated by means of the Thompson-Thames-Mastin transformation. An approximate factorization scheme with cyclic acceleration parameters was employed to solve a finite difference analog of the stream function equation. The vorticity equation was numerically solved with a modified version of the classical alternating direction implicit scheme. Although no quantitative assessment of solution accuracy can be made, numerical results for variations in incremental wall pressure around the slat are at least qualitatively similar to some experimental results.

Cooke, C. H.↗

NASA NTRS1996

Analysis of Two-Phase Flow in Damper Seals for Cryogenic Turbopumps

Cryogenic damper seals operating close to the liquid-vapor region (near the critical point or slightly su-cooled) are likely to present two-phase flow conditions. Under single phase flow conditions the mechanical energy conveyed to the fluid increases its temperature and causes a phase change when the fluid temperature reaches the saturation value. A bulk-flow analysis for the prediction of the dynamic force response of damper seals operating under two-phase conditions is presented as: all-liquid, liquid-vapor, and all-vapor, i.e. a 'continuous vaporization' model. The two phase region is considered as a homogeneous saturated mixture in thermodynamic equilibrium. Th flow in each region is described by continuity, momentum and energy transport equations. The interdependency of fluid temperatures and pressure in the two-phase region (saturated mixture) does not allow the use of an energy equation in terms of fluid temperature. Instead, the energy transport is expressed in terms of fluid enthalpy. Temperature in the single phase regions, or mixture composition in the two phase region are determined based on the fluid enthalpy. The flow is also regarded as adiabatic since the large axial velocities typical of the seal application determine small levels of heat conduction to the walls as compared to the heat carried by fluid advection. Static and dynamic force characteristics for the seal are obtained from a perturbation analysis of the governing equations. The solution expressed in terms of zeroth and first order fields provide the static (leakage, torque, velocity, pressure, temperature, and mixture composition fields) and dynamic (rotordynamic force coefficients) seal parameters. Theoretical predictions show good agreement with experimental leakage pressure profiles, available from a Nitrogen at cryogenic temperatures. Force coefficient predictions for two phase flow conditions show significant fluid compressibility effects, particularly for mixtures with low mass content of vapor. Under these conditions, an increase on direct stiffness and reduction of whirl frequency ratio are shown to occur. Prediction of such important effects will motivate experimental studies as well as a more judicious selection of the operating conditions for seals used in cryogenic turbomachinery.

Arauz, Grigory L.↗

NASA NTRS1994

Whirl measurements on leakage flows in turbomachine models

The beneficial effects claimed for whirl control devices demonstrate that the dynamic behavior of rotors is influenced by the fluid whirl in shaft and balance drum seals. The present paper reports results from two series of experiments, the first on the factors affecting the whirl at the seal inlet, and the second on the variation of whirl velocity along the seal. In both cases the LDA measurement technique required the clearance between the fixed and rotating parts of the models to be substantially greater than occurs in real machines, but the results are indicative nevertheless. Experimental and theoretical results are given for the radial distribution of whirl velocity in the gap between impeller shroud and pump casing. Results of tests with modified stator surfaces are also shown. This work leads naturally into the second series of experiments where some preliminary measurements of velocity distribution in the clearance between a fixed stator and a rotating shaft are reported for a range of inlet whirl conditions.

Addlesee, A. J.↗