SEARCH · Search NASA

Results for “two-phase simulation”

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 181 records · Page 10

All sources NASA NTRS DOE OSTI

NASA NTRS1988

Fluid spray simulation with two-fluid nozzles

Two-phase interacting flow inside a two-fluid fuel atomizer was investigated and a correction of aerodynamic and liquid-surface forces with characteristic drop diameter was obtained for liquid-jet breakup in Mach 1 gas flow. Nitrogen gas mass-flux was varied from 6 to 50 g/sq cm sec by using four differently sized two-fluid atomizers with nozzle diameters varyig from 0.32 to 0.56 cm. The correlation was derived by using the acoustic gas velocity, V sub c, as a basic parameter in defining and evaluating the dimensionless product of the Weber (We) and Reynolds (Re) numbers. By using the definition of WeRe, it was found that the ratio of orifice diameter to Sauter mean drop diameter could be correlated with the dimensionless ratio WeRe and the gas to liquid density ratio.

Ingebo, Robert D.↗

NASA NTRS1989

Fuel spray simulation with two-fluid nozzles

Two-phase interacting flow inside a two-fluid fuel atomizer was investigated and a correction of aerodynamic and liquid-surface forces with characteristic drop diameter was obtained for liquid-jet breakup in Mach 1 gas flow. Nitrogen gas mass-flux was varied from 6 to 50 g/sq cm sec by using four differently sized two-fluid atomizers with nozzle diameters varying from 0.32 to 0.56 cm. The correlation was derived by using the acoustic gas velocity, V sub c, as a basic parameter in defining and evaluating the dimensionless product of the Weber (We) and Reynolds (Re) numbers. By using the definition of WeRe, it was found that the ratio of orifice diameter to Sauter mean drop diameter could be correlated with the dimensionless ratio WeRe and the gas to liquid density ratio.

Ingebo, Robert D.↗

NASA NTRS1988

A knowledge-based approach to identification and adaptation in dynamical systems control

Artificial intelligence techniques are applied to the problems of model form and parameter identification of large-scale dynamic systems. The object-oriented knowledge representation is discussed in the context of causal modeling and qualitative reasoning. Structured sets of rules are used for implementing qualitative component simulations, for catching qualitative discrepancies and quantitative bound violations, and for making reconfiguration and control decisions that affect the physical system. These decisions are executed by backward-chaining through a knowledge base of control action tasks. This approach was implemented for two examples: a triple quadrupole mass spectrometer and a two-phase thermal testbed. Results of tests with both of these systems demonstrate that the software replicates some or most of the functionality of a human operator, thereby reducing the need for a human-in-the-loop in the lower levels of control of these complex systems.

Glass, B. J.↗

NASA NTRS1989

A two-phase methodology for technology selection and system design

A two-phase methodology that can be used to guide R&D managers in the evaluation and selection of competing technologies is presented. Deterministic multiattribute utility theory is used in the first phase to rank the technological alternatives; the example presented involves the evaluation of electric and hybrid passenger vehicles. In all, 39 individuals from eight automotive firms were interviewed to assess their risk preferences and attitudes toward the vehicle design. In the second phase, the decision-maker must allocate a fixed amount of resources to different projects for the technology selected, some of which may be undertaken in parallel, to maximize a given measure of performance. When parallel funding is pursued the best outcome is chosen. The problem is formulated as a probabilistic network and solved heuristically using Monte Carlo simulation. Results are presented for two decision-makers and three budget options. In each case, the heuristic finds the optimal allocation of funds.

Bard, Jonathan F.↗

NASA NTRS1989

Ground evaluation of seeding an in-flight wingtip vortex using infrared imaging flow visualization technique

An experimental simulation of an in-flight wingtip vortical flow visualization technique uses infrared imaging to observe strong and concentrated vortices. This experiment is phase 1 of a two-phase infrared evaluation program. The system includes a vortex generator (model 320 Vortec Vortex Tube) which generates the required vortex. The mouth of the unit is mounted close to the free end of a half-inch diameter, sixteen and a half foot long stainless steel tubing (sized after tubing currently installed in the wings of an experimental Beechcraft Sundowner 180 aircraft). Dichloro difluoromethane (Freon-12) is entrained into the generated vortex. A breakdown of the vortices is indicated by the rapid diffusion and the resulting pattern is tracked using the infrared imager and video systems. Flow rates (volume and mass) are estimated at the laboratory and proposed flight conditions. The nominal flight altitude is expected to be 2500 feet.

Akinyanju, Ted↗

NASA NTRS1992

A two-phase restricted equilibrium model for combustion of metalized solid propellants

An Eulerian-Lagrangian two-phase approach was adopted to model the multi-phase reacting internal flow in a solid rocket with a metalized propellant. An Eulerian description was used to analyze the motion of the continuous phase which includes the gas as well as the small (micron-sized) particulates, while a Lagrangian description is used for the analysis of the discrete phase which consists of the larger particulates in the motor chamber. The particulates consist of Al and Al2O3 such that the particulate composition is 100 percent Al at injection from the propellant surface with Al2O3 fraction increasing due to combustion along the particle trajectory. An empirical model is used to compute the combustion rate for agglomerates while the continuous phase chemistry is treated using chemical equilibrium. The computer code was used to simulate the reacting flow in a solid rocket motor with an AP/HTPB/Al propellant. The computed results show the existence of an extended combustion zone in the chamber rather than a thin reaction region. The presence of the extended combustion zone results in the chamber flow field and chemical being far from isothermal (as would be predicted by a surface combustion assumption). The temperature in the chamber increases from about 2600 K at the propellant surface to about 3350 K in the core. Similarly the chemical composition and the density of the propellant gas also show spatially non-uniform distribution in the chamber. The analysis developed under the present effort provides a more sophisticated tool for solid rocket internal flow predictions than is presently available, and can be useful in studying apparent anomalies and improving the simple correlations currently in use. The code can be used in the analysis of combustion efficiency, thermal load in the internal insulation, plume radiation, etc.

Sabnis, J. S.↗

NASA NTRS1993

Computational techniques for flows with finite-rate condensation

A computational method to simulate the inviscid two-dimensional flow of a two-phase fluid was developed. This computational technique treats the gas phase and each of a prescribed number of particle sizes as separate fluids which are allowed to interact with one another. Thus, each particle-size class is allowed to move through the fluid at its own velocity at each point in the flow field. Mass, momentum, and energy are exchanged between each particle class and the gas phase. It is assumed that the particles do not collide with one another, so that there is no inter-particle exchange of momentum and energy. However, the particles are allowed to grow, and therefore, they may change from one size class to another. Appropriate rates of mass, momentum, and energy exchange between the gas and particle phases and between the different particle classes were developed. A numerical method was developed for use with this equation set. Several test cases were computed and show qualitative agreement with previous calculations.

Candler, Graham V.↗

NASA NTRS1995

NASAwide electronic publishing system: Electronic printing and duplicating, stage-2 evaluation report (GSFC)

The NASA Scientific and Technical Information Office was assigned the responsibility to continue with the expansion of the NASAwide networked electronic duplicating effort by including the Goddard Space Flight Center (GSFC) as an additional node to the existing configuration of networked electronic duplicating systems within NASA. The subject of this report is the evaluation of a networked electronic duplicating system which meets the duplicating requirements and expands electronic publishing capabilities without increasing current operating costs. This report continues the evaluation reported in 'NASA Electronic Publishing System - Electronic Printing and Duplicating Evaluation Report' (NASA TM-106242) and 'NASA Electronic Publishing System - Stage 1 Evaluation Report' (NASA TM-106510). This report differs from the previous reports through the inclusion of an external networked desktop editing, archival, and publishing functionality which did not exist with the previous networked electronic duplicating system. Additionally, a two-phase approach to the evaluation was undertaken; the first was a paper study justifying a 90-day, on-site evaluation, and the second phase was to validate, during the 90-day evaluation, the cost benefits and productivity increases that could be achieved in an operational mode. A benchmark of the functionality of the networked electronic publishing system and external networked desktop editing, archival, and publishing system was performed under a simulated daily production environment. This report can be used to guide others in determining the most cost effective duplicating/publishing alternative through the use of cost/benefit analysis and return on investment techniques. A treatise on the use of these techniques can be found by referring to 'NASA Electronic Publishing System -Cost/Benefit Methodology' (NASA TM-106662).

Tuey, Richard C.↗

NASA NTRS1996

Phase Distribution Phenomena for Simulated Microgravity Conditions: Experimental Work

This report summarizes the work accomplished at Rensselaer to study phase distribution phenomenon under simulated microgravity conditions. Our group at Rensselaer has been able to develop sophisticated analytical models to predict phase distribution in two-phase flows under variety of conditions. These models are based on physics and data obtained from carefully controlled experiments that are being conducted here. These experiments also serve to verify the models developed.

Singhal, Maneesh↗

NASA NTRS1996

Phase Distribution Phenomena for Simulated Microgravity Conditions: Experimental Work

This report summarizes the work accomplished at Rensselaer to study phase distribution phenomenon under simulated microgravity conditions. Our group at Rensselaer has been able to develop sophisticated analytical models to predict phase distribution in two-phase flows under a variety of conditions. These models are based on physics and data obtained from carefully controlled experiments that are being conducted here. These experiments also serve to verify the models developed.

Singhal, Maneesh↗

NASA NTRS1996

Multigrid Methods for Fully Implicit Oil Reservoir Simulation

In this paper we consider the simultaneous flow of oil and water in reservoir rock. This displacement process is modeled by two basic equations: the material balance or continuity equations and the equation of motion (Darcy's law). For the numerical solution of this system of nonlinear partial differential equations there are two approaches: the fully implicit or simultaneous solution method and the sequential solution method. In the sequential solution method the system of partial differential equations is manipulated to give an elliptic pressure equation and a hyperbolic (or parabolic) saturation equation. In the IMPES approach the pressure equation is first solved, using values for the saturation from the previous time level. Next the saturations are updated by some explicit time stepping method; this implies that the method is only conditionally stable. For the numerical solution of the linear, elliptic pressure equation multigrid methods have become an accepted technique. On the other hand, the fully implicit method is unconditionally stable, but it has the disadvantage that in every time step a large system of nonlinear algebraic equations has to be solved. The most time-consuming part of any fully implicit reservoir simulator is the solution of this large system of equations. Usually this is done by Newton's method. The resulting systems of linear equations are then either solved by a direct method or by some conjugate gradient type method. In this paper we consider the possibility of applying multigrid methods for the iterative solution of the systems of nonlinear equations. There are two ways of using multigrid for this job: either we use a nonlinear multigrid method or we use a linear multigrid method to deal with the linear systems that arise in Newton's method. So far only a few authors have reported on the use of multigrid methods for fully implicit simulations. Two-level FAS algorithm is presented for the black-oil equations, and linear multigrid for two-phase flow problems with strong heterogeneities and anisotropies is studied. Here we consider both possibilities. Moreover we present a novel way for constructing the coarse grid correction operator in linear multigrid algorithms. This approach has the advantage in that it preserves the sparsity pattern of the fine grid matrix and it can be extended to systems of equations in a straightforward manner. We compare the linear and nonlinear multigrid algorithms by means of a numerical experiment.

Molenaar, J.↗

NASA NTRS2000

Computations of Boiling in Microgravity

The absence (or reduction) of gravity, can lead to major changes in boiling heat transfer. On Earth, convection has a major effect on the heat distribution ahead of an evaporation front, and buoyancy determines the motion of the growing bubbles. In microgravity, convection and buoyancy are absent or greatly reduced and the dynamics of the growing vapor bubbles can change in a fundamental way. In particular, the lack of redistribution of heat can lead to a large superheat and explosive growth of bubbles once they form. While considerable efforts have been devoted to examining boiling experimentally, including the effect of microgravity, theoretical and computational work have been limited. Here, the growth of boiling bubbles is studied by direct numerical simulations where the flow field is fully resolved and the effects of inertia, viscosity, surface deformation, heat conduction and convection, as well as the phase change, are fully accounted for. Boiling involves both fluid flow and heat transfer and thus requires the solution of the Navier-Stokes and the energy equations. The numerical method is based on writing one set of governing transport equations which is valid in both the liquid and vapor phases. This local, single-field formulation incorporates the effect of the interface in the governing equations as source terms acting only at the interface. These sources account for surface tension and latent heat in the equations for conservation of momentum and energy as well as mass transfer across the interface due to phase change. The single-field formulation naturally incorporates the correct mass, momentum and energy balances across the interface. Integration of the conservation equations across the interface directly yields the jump conditions derived in the local instant formulation for two-phase systems. In the numerical implementation, the conservation equations for the whole computational domain (both vapor and liquid) are solved using a stationary grid and the phase boundary is followed by a moving unstructured two-dimensional grid. While two-dimensional simulations have been used for preliminary studies and to examine the resolution requirement, the focus is on fully three-dimensional simulations. The numerical methodology, including the parallelization and grid refinement strategy is discussed, and preliminary results shown. For buoyancy driven flow, the heat transfer is in good agreement with experimental correlations. The changes when gravity is turned off and/or fluid shear is added are discussed, as well as the difference between simulations of a layer freely releasing bubbles versus simulations using only one wavelength initial perturbation. Figure 1 shows the early stages of the formation of a three-dimensional bubble from a thin vapor layer. The boundary conditions are periodic in the x and y direction, the bottom is a hot and the top allows a free outflow. The jagged edge of the surface close to the bottom of the computational domain is due to some of the surface elements being on the other side of the domain and some elements not plotted by our plotting routine. In the second figure, we show the temperature distribution through two perpendicular planes.

Tryggvason, G.↗

NASA NTRS2001

Robust Strategy for Rocket Engine Health Monitoring

Monitoring the health of rocket engine systems is essentially a two-phase process. The acquisition phase involves sensing physical conditions at selected locations, converting physical inputs to electrical signals, conditioning the signals as appropriate to establish scale or filter interference, and recording results in a form that is easy to interpret. The inference phase involves analysis of results from the acquisition phase, comparison of analysis results to established health measures, and assessment of health indications. A variety of analytical tools may be employed in the inference phase of health monitoring. These tools can be separated into three broad categories: statistical, rule based, and model based. Statistical methods can provide excellent comparative measures of engine operating health. They require well-characterized data from an ensemble of "typical" engines, or "golden" data from a specific test assumed to define the operating norm in order to establish reliable comparative measures. Statistical methods are generally suitable for real-time health monitoring because they do not deal with the physical complexities of engine operation. The utility of statistical methods in rocket engine health monitoring is hindered by practical limits on the quantity and quality of available data. This is due to the difficulty and high cost of data acquisition, the limited number of available test engines, and the problem of simulating flight conditions in ground test facilities. In addition, statistical methods incur a penalty for disregarding flow complexity and are therefore limited in their ability to define performance shift causality. Rule based methods infer the health state of the engine system based on comparison of individual measurements or combinations of measurements with defined health norms or rules. This does not mean that rule based methods are necessarily simple. Although binary yes-no health assessment can sometimes be established by relatively simple rules, the causality assignment needed for refined health monitoring often requires an exceptionally complex rule base involving complicated logical maps. Structuring the rule system to be clear and unambiguous can be difficult, and the expert input required to maintain a large logic network and associated rule base can be prohibitive.

Santi, L. Michael↗

NASA NTRS2002

Measurements of Turbulence Attenuation by a Dilute Dispersion of Solid Particles in Homogeneous Isotropic Turbulence

This research addresses turbulent gas flows laden with fine solid particles at sufficiently large mass loading that strong two-way coupling occurs. By two-way coupling we mean that the particle motion is governed largely by the flow, while the particles affect the gas-phase mean flow and the turbulence properties. Our main interest is in understanding how the particles affect the turbulence. Computational techniques have been developed which can accurately predict flows carrying particles that are much smaller than the smallest scales of turbulence. Also, advanced computational techniques and burgeoning computer resources make it feasible to fully resolve very large particles moving through turbulent flows. However, flows with particle diameters of the same order as the Kolmogorov scale of the turbulence are notoriously difficult to predict. Some simple flows show strong turbulence attenuation with reductions in the turbulent kinetic energy by up to a factor of five. On the other hand, some seemingly similar flows show almost no modification. No model has been proposed that allows prediction of when the strong attenuation will occur. Unfortunately, many technological and natural two-phase flows fall into this regime, so there is a strong need for new physical understanding and modeling capability. Our objective is to study the simplest possible turbulent particle-laden flow, namely homogeneous, isotropic turbulence with a uniform dispersion of monodisperse particles. We chose such a simple flow for two reasons. First, the simplicity allows us to probe the interaction in more detail and offers analytical simplicity in interpreting the results. Secondly, this flow can be addressed by numerical simulation, and many research groups are already working on calculating the flow. Our detailed data can help guide some of these efforts. By using microgravity, we can further simplify the flow to the case of no mean velocity for either the turbulence or the particles. In fact the addition of gravity as a variable parameter may help us to better understand the physics of turbulence attenuation. The experiments are conducted in a turbulence chamber capable of producing stationary or decaying isotropic turbulence with nearly zero mean flow and Taylor microscale Reynolds numbers up to nearly 500. The chamber is a 410 mm cubic box with the corners cut off to make it approximately spherical. Synthetic jet turbulence generators are mounted in each of the eight corners of the box. Each generator consists of a loudspeaker forcing a plenum and producing a pulsed jet through a 20 mm diameter orifice. These synthetic jets are directed into ejector tubes pointing towards the chamber center. The ejector tubes increase the jet mass flow and decrease the velocity. The jets then pass through a turbulence grid. Each of the eight loudspeakers is forced with a random phase and frequency. The resulting turbulence is highly Isotropic and matches typical behavior of grid turbulence. Measurements of both phases are acquired using particle image velocimetry (PIV). The gas is seeded with approximately 1 micron diameter seeding particles while the solid phase is typically 150 micron diameter spherical glass particles. A double-pulsed YAG laser and a Kodak ES-1.0 10-bit PIV camera provide the PIV images. Custom software is used to separate the images into individual images containing either gas-phase tracers or large particles. Modern high-resolution PIV algorithms are then used to calculate the velocity field. A large set of image pairs are acquired for each case, then the results are averaged both spatially and over the ensemble of acquired images. The entire apparatus is mounted in two racks which are carried aboard NASA's KC-135 Flying Microgravity Laboratory. The rack containing the turbulence chamber, the laser head, and the camera floats freely in the airplane cabin (constrained by competent NASA personnel) to minimize g-jitter.

Eaton, John↗

NASA NTRS2004

Fluid Physics and Transport Phenomena in a Simulated Reduced Gravity Environment

We describe a ground-based apparatus that allows the cancellation of gravity on a fluid using magnetic forces. The present system was designed for liquid oxygen studies over the range 0.001 - 5 g s. This fluid is an essential component of any flight mission using substantial amounts of liquid propellant, especially manned missions. The apparatus has been used to reduce the hydrostatic compression near the oxygen critical point and to demonstrate inverted phase separation. It could also be used to study pool boiling and two-phase heat transfer in Martian, Lunar or near-zero gravity, as well as phenomena such as Marangoni flow and convective instabilities. These studies would contribute directly to the reliability and optimization of the Moon and Mars flight programs.

Lipa, J.↗

NASA NTRS2004

Thermal performance of a multi-evaporator loop heat pipe with thermal masses and thermal electrical coolers

This paper describes thermal performance of a loop heat pipe (LHP) with two evaporators and two condensers in ambient testing. Each evaporator has an outer diameter of 15mm and a length of 76mm, and has an integral compensation chamber (CC). An aluminum mass of 500 grams is attached to each evaporator to simulate the instrument mass. A thermal electric cooler (TEC) is installed on each CC to provide heating as well as cooling for CC temperature control. A flow regulator is installed in the condenser section to prevent vapor from going back to the evaporators in the event that one of condenser is fully utilized. Ammonia was used ad the working fluid. Tests conducted included start-up, power cycle, heat load sharing, sink temperature cycle, operating temperature control with TECs, and capillary limit tests. Experimental data showed that the loop could start with a heat load of less than 1OW even with added thermal masses. The loop operated stably with even and uneven evaporator heat loads, and even and uneven condenser sink temperatures. The operating temperature could be controlled within +/-0.5K of the set point temperature using either or both TECs, and the required TEC control heater power was less than 2W under most test conditions. Heat load sharing between the two evaporators was also successfully demonstrated. The loop had a heat transport capability of 120W to 140W, and could recover from a dry-out when the heat load was reduced. The 500-gram aluminum mass on each evaporator had a negligible effect on the loop operation. Existing LHPs servicing the orbiting spacecraft have a single evaporator with an outer diameter of about 25mm. Important performance characteristics demonstrated by this LHP included: 1) Operation of an LHP with 15mm diameter evaporators; 2) Robustness and reliability of an LHP with multiple evaporators and multiple condensers under various test conditions; 3) Heat load sharing among LHP evaporators; 4) Effectiveness of TECs in controlling the LHP operating temperature; and 5) Effectiveness of the flow regulator in preventing vapor from going back the evaporators.

loop heat pipes↗

NASA NTRS2006

Scalar Dissipation Modeling for Passive and Active Scalars: a priori Study Using Direct Numerical Simulation

Transitional databases from Direct Numerical Simulation (DNS) of three-dimensional mixing layers for single-phase flows and two-phase flows with evaporation are analyzed and used to examine the typical hypothesis that the scalar dissipation Probability Distribution Function (PDF) may be modeled as a Gaussian. The databases encompass a single-component fuel and four multicomponent fuels, two initial Reynolds numbers (Re), two mass loadings for two-phase flows and two free-stream gas temperatures. Using the DNS calculated moments of the scalar-dissipation PDF, it is shown, consistent with existing experimental information on single-phase flows, that the Gaussian is a modest approximation of the DNS-extracted PDF, particularly poor in the range of the high scalar-dissipation values, which are significant for turbulent reaction rate modeling in non-premixed flows using flamelet models. With the same DNS calculated moments of the scalar-dissipation PDF and making a change of variables, a model of this PDF is proposed in the form of the (beta)-PDF which is shown to approximate much better the DNS-extracted PDF, particularly in the regime of the high scalar-dissipation values. Several types of statistical measures are calculated over the ensemble of the fourteen databases. For each statistical measure, the proposed (beta)-PDF model is shown to be much superior to the Gaussian in approximating the DNS-extracted PDF. Additionally, the agreement between the DNS-extracted PDF and the (beta)-PDF even improves when the comparison is performed for higher initial Re layers, whereas the comparison with the Gaussian is independent of the initial Re values. For two-phase flows, the comparison between the DNS-extracted PDF and the (beta)-PDF also improves with increasing free-stream gas temperature and mass loading. The higher fidelity approximation of the DNS-extracted PDF by the (beta)-PDF with increasing Re, gas temperature and mass loading bodes well for turbulent reaction rate modeling.

scalar dissipation modeling↗

NASA NTRS2007

Finite Element Simulation of Solid Rocket Booster Separation Motors During Motor Firing

One of the toughest challenges facing Solid Rocket Booster (SRB) engineers is to ensure that any design changes made to the Shuttle-Derived Booster Separation Motors (BSM) for future space exploration vehicles is able to withstand the increasingly hostile motor firing environment without cracking its critical component - the graphite throat. This paper presents a critical analysis methodology and techniques for assessing effects of BSM design changes with great accuracy and precision. For current Space Shuttle operation, the motor firing occurs at SRB separation - approximately 125 seconds after Shuttle launch at an altitude of about 28 miles. The motor operation event lasts about two seconds, however, the surface temperature of the graphite throat increases approximately 3400 F in less than one second with a corresponding increase in surface pressure of approximately 2200 pounds per square inch (psi) in less than one-tenth of a second. To capture this process fully and accurately, a two-phase sequentially coupled thermal-mechanical finite element approach was developed. This method allows the time- and location-dependent pressure fields to interact with the spatial-temporal thermal fields throughout the operation. The material properties of graphite throat are orthotropic and temperature-dependent. The analysis involves preload and multiple body contacts.

Yu. Weiping↗