Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “two-phase flows”

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 577 records · Page 32

Computer simulation of dust grain evolution

The latest results are reported from a Monte Carlo code that is being developed at NASA Ames. The goal of this program, is to derive from the observed and presumed properties of the interstellar medium (ISM) the following information: (1) the size spectrum of interstellar dust; (2) the chemical structure of interstellar dust; (3) interstellar abundances; and (4) the lifetime of a dust grain in the ISM. Presently this study is restricted to refractory interstellar material, i.e., the formation and destruction of ices are not included in the program. The program is embedded in an analytic solution for the bulk chemical evolution of a two-phase interstellar medium in which stars are born in molecular clouds, but new nucleosynthesis products and stellar return are entered into a complementary intercloud medium. The well-mixed matter of each interstellar phase is repeatedly cycled stochastically through the complementary phase and back. Refractory dust is created by thermal condensation as stellar matter flows away from sites of nucleosynthesis such as novae and supernovae and/or from the matter returned from evolved intermediate stars. The history of each particle is traced by standard Monte Carlo techniques as it is sputtered and fragmented by supernova shock waves in the intercloud medium. It also accretes an amorphous mantle of gaseous refractory atoms when its local medium joins with the molecular cloud medium. Finally it encounters the possibility of astration (destruction by star formation) within the molecular clouds.

Liffman, K.↗

Acquisition and correlation of cryogenic nitrogen mass flow data through a multiple orifice Joule-Thomson device

Liquid nitrogen mass flow rate, pressure drop, and temperature drop data were obtained for a series of multiple orifice Joule-Thomson devices, known as Visco Jets, over a wide range of flow resistance. The test rig used to acquire the data was designed to minimize heat transfer so that fluid expansion through the Visco Jets would be isenthalpic. The data include a range of fluid inlet pressures from 30 to 60 psia, fluid inlet temperatures from 118 to 164 R, outlet pressures from 2.8 to 55.8 psia, outlet temperatures from 117 to 162 R and flow rate from 0.04 to 4.0 lbm/hr of nitrogen. A flow rate equation supplied by the manufacturer was found to accurately predict single-phase (noncavitating) liquid nitrogen flow through the Visco Jets. For cavitating flow, the manufacturer's equation was found to be inaccurate. Greatly improved results were achieved with a modified version of the single-phase equation. The modification consists of a multiplication factor to the manufacturer's equation equal to one minus the downstream quality on an isenthalpic expansion of the fluid across the Visco Jet. For a range of flow resistances represented by Visco Jet Lohm ratings between 17,600 and 80,000, 100 percent of the single-phase data and 85 percent of the two-phase data fall within + or - 10 percent of predicted values.

Papell, S. Stephen↗

The modelling of heat, mass and solute transport in solidification systems

The aim of this paper is to explore the range of possible one-phase models of binary alloy solidification. Starting from a general two-phase description, based on the two-fluid model, three limiting cases are identified which result in one-phase models of binary systems. Each of these models can be readily implemented in standard single phase flow numerical codes. Differences between predictions from these models are examined. In particular, the effects of the models on the predicted macro-segregation patterns are evaluated.

Voller, V. R.↗

Effect of a crystal-melt interface on Taylor-vortex flow

The linear stability of circular Couette flow between concentric infinite cylinders is considered for the case that the stationary outer cylinder is a crystal-melt interface rather than a rigid surface. A radial temperature difference is maintained across the liquid gap, and equations for heat transport in the crystal and melt phases are included to extend the ordinary formulation of this problem. The stability of this two-phase system depends on the Prandtl number. For small Prandtl number the linear stability of the two-phase system is given by the classical results for a rigid-walled system. For increasing values of the Prandtl number, convective heat transport becomes significant and the system becomes increasingly less stable. Previous results in a narrow-gap approximation are extended to the case of a finite gap, and both axisymmetric and nonaxisymmetric disturbance modes are considered. The two-phase system becomes less stable as the finite gap tends to the narrow-gap limit. The two-phase system is more stable to nonaxisymmetric modes with azimuthal wavenumber n = 1; the stability of these n = 1 modes is sensitive to the latent heat of fusion.

Mcfadden, G. B.↗

Active thermal control systems for lunar and Martian exploration

Several ATCS options including heat pumps, radiator shading devices, and single-phase flow loops were considered. The ATCS chosen for both lunar and Martian habitats consists of a heat pump integral with a nontoxic fluid acquisition and transport loop, and vertically oriented modular reflux-boiler radiators. The heat pump operates only during the lunar day. The lunar and Martian transfer vehicles have an internal single-phase water-acquisition loop and an external two-phase ammonia rejection system with rotating inflatable radiators. The lunar and Martian excursion vehicles incorporate internal single-phase water acquisition, which is connected via heat exchangers to external body-mounted single-phase radiators. A water evaporation system is used for the transfer vehicles during periods of high heating.

Ewert, Michael K.↗

Heavy rain effects on airplane performance

The objective is to determine if the aerodynamic characteristics of an airplane are altered while flying in the rain. Wind-tunnel tests conducted at the NASA Langley Research Center (LaRC) have shown losses in maximum lift, reduction in stall angle, and increases in drag when a wing is placed in a simulated rain spray. For these tests the water spray concentration used represented a very heavy rainfall. A lack of definition of the scaling laws for aerodynamic testing in a two-phase, two-component flow makes interpolation of the wind-tunnel test uncertain. Tests of a large-scale wing are to be conducted at the LaRC. The large-scale wing is mounted on top of the Aircraft Landing Dynamics Facility (ALDF) carriage. This carriage (which is 70-foot long, 30-foot wide, and 30-foot high) is propelled with the wing model attached down a 3000-foot long test track by a water jet at speeds of up to 170 knots. A simulated rain spray system has been installed along 500 feet of the test track and can simulate rain falls from 2 to 40 inches/hour. Operational checks are underway and the initial tests should be completed by the Fall of 1989.

Dunham, R. E., Jr.↗

Two-dimensional simulation of a two-phase, regenerative pumped radiator loop utilizing direct contact heat transfer with phase change

An innovative pumped loop concept for 600 K space power system radiators utilizing direct contact heat transfer, which facilitates repeated startup/shutdown of the power system without complex and time-consuming coolant thawing during power startup, is under development. The heat transfer process with melting/freezing of Li in an NaK flow was studied through two-dimensional time-dependent numerical simulations to characterize and predict the Li/NaK radiator performance during startup (thawing) and shutdown (cold-trapping). Effects of system parameters and the criteria for the plugging domain are presented together with temperature distribution patterns in solid Li and subsequent melting surface profile variations in time.

Rhee, Hyop S.↗

Modeling the use of a binary mixture as a control scheme for two-phase thermal systems

Two-phase thermal loops using mechanical pumps, capillary pumps, or a combination of the two have been chosen as the main heat transfer systems for the space station. For these systems to operate optimally, the flow rate in the loop should be controlled in response to the vapor/liquid ratio leaving the evaporator. By substituting a mixture of two non-azeotropic fluids in place of the single fluid normally used in these systems, it may be possible to monitor the temperature of the exiting vapor and determine the vapor/liquid ratio. The flow rate would then be adjusted to maximize the load capability with minimum energy input. A FLUINT model was developed to study the system dynamics of a hybrid capillary pumped loop using this type of control and was found to be stable under all the test conditions.

Benner, S. M.↗

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

Program Computes Flows Of Fluids And Heat

SINDA'85/FLUINT incorporates lumped-parameter-network and one-dimensional-flow mathematical models. System enables analysis of mutual influences of thermal and flow phenomena. Offers two finite-difference numerical solution techniques: forward-difference explicit approximation and Crank-Nicholson approximation. Enables simulation of nonuniform heating and facilitates mathematical modeling of thin-walled heat exchangers. Ability to model nonequilibrium behavior within two-phase volumes included. Recent changes in program improve modeling of real evaporator pumps and other capillary-assist evaporators. Written in FORTRAN 77.

Cullimore, Brent↗

Comparison of the initial development of shear layers in two-dimensional and axisymmetric ejector configurations

A two-phase experimental investigation designed to study the development of shear layers in axisymmetric and two-dimensional single-nozzle ejectors has been completed. In this study, combinations of similar and dissimilar gases were used as the supersonic primary and subsonic secondary. Test cases included combinations of air/air, argon/air and helium/air as the supersonic primary and subsonic secondary, respectively. Similar flow conditions were studied for each ejector configuration. Mixing of the gases occurred in a constant-area tube, where the inlet pressure was maintained at 34.5 kPa. The cases studied resulted in convective Mach numbers that range between 0.06 and 1.9. The data gathered shows differences between the initial shear-layer development for the two ejector geometries, and also between the different test cases studied for each ejector configuration. The measured growth rates for the axisymmetric ejector are more than twice those measured for the two-dimensional ejector. However, in both cases the results show that compressibility has a reducing effect on the growth rate. Further, in the region immediately after the inlet to the mixing tube, compressibility seems to affect the ejector shear layers in a manner similar to that of two-stream two-dimensional mixing layers.

Dufflocq, M.↗

Advanced Active Thermal Control Systems Architecture Study

The Johnson Space Center (JSC) initiated a dynamic study to determine possible improvements available through advanced technologies (not used on previous or current human vehicles), identify promising development initiatives for advanced active thermal control systems (ATCS's), and help prioritize funding and personnel distribution among many research projects by providing a common basis to compare several diverse technologies. Some technologies included were two-phase thermal control systems, light-weight radiators, phase-change thermal storage, rotary fluid coupler, and heat pumps. JSC designed the study to estimate potential benefits from these various proposed and under-development thermal control technologies for five possible human missions early in the next century. The study compared all the technologies to a baseline mission using mass as a basis. Each baseline mission assumed an internal thermal control system; an external thermal control system; and aluminum, flow-through radiators. Solar vapor compression heat pumps and light-weight radiators showed the greatest promise as general advanced thermal technologies which can be applied across a range of missions. This initial study identified several other promising ATCS technologies which offer mass savings and other savings compared to traditional thermal control systems. Because the study format compares various architectures with a commonly defined baseline, it is versatile and expandable, and is expected to be updated as needed.

Hanford, Anthony J.↗

Premixed Atmosphere and Convection Influences on Flame Inhibition and Combustion (Pacific)

Flame spread over flat solid fuel beds is a useful paradigm for studying the behavior of more complex two-phase nonpremixed flames. For practical applications, two of the most important elements of flame spreading are the effects of (1) the ambient atmosphere (e.g. pressure and composition) and (2) the flow environment on the spread rate and extinction conditions. Concerning (1), studies of flame spread in vitiated air and non-standard atmospheres such as those found in undersea vessels and spacecraft are particularly important for the assessment of fire hazards in these environments as well as determination of the effectiveness of fire suppressants. Concerning (2), the flow environment may vary widely even when no forced flow is present because of buoyancy effects. Consequently, the goal of this work is to employ microgravity (micro g) experiments to extend previous studies of the effects of ambient atmosphere and the flow environment on flame spread through the use of microgravity (micro g) experiments. Because of the considerable differences between upward (concurrent-flow) and downward (opposed-flow) flame spread at 1g (Williams, 1976, Fernandez-Pello, 1984), in this work both upward and downward 1g spread are tested. Two types of changes to the oxidizing atmosphere are considered in this work. One is the addition of sub-flammability-limit concentrations of a gaseous fuel ('partially premixed' atmospheres). This is of interest because in fires in enclosures, combustion may occur under poorly ventilated conditions, so that oxygen is partially depleted from the air and is replaced by combustible gases such as fuel vapors, H2 or CO. Subsequent fire spread over the solid fuel could occur under conditions of varying oxygen and gaseous fuel content. The potential significance of flame spread under vitiated or partially premixed conditions has been noted previously (Beyler, 1984). The second change is the diluent type, which affects the radiative properties of the mixture as well as the Lewis number (Le) of the reactants in the atmosphere, which for oxygen is defined as the thermal diffusivity of the bulk mixture divided by the mass diffusivity of oxygen into the bulk mixture. Understanding the effect of diluent type is desirable because in some undersea and spaceborne habitations, it is desirable to use diluent gases other than nitrogen. Prior experiments have shown that both radiation (Bhattacharjee and Altenkirch, 1993) and Lewis number (Zhang et al, 1992) effects are important in flame spreading problems.

Honda, Linton K.↗

Microstructural Development during Directional Solidification of Peritectic Alloys

A thorough understanding of the microstructures produced through solidification in peritectic systems has yet to be achieved, even though a large number of industrially and scientifically significant materials are in this class. One type of microstructure frequently observed during directional solidification consists of alternating layers of primary solid and peritectic solid oriented perpendicular to the growth direction. This layer formation is usually reported for alloy compositions within the two-phase region of the peritectic isotherm and for temperature gradient and growth rate conditions that result in a planar solid-liquid interface. Layered growth in peritectic alloys has not previously been characterized on a quantitative basis, nor has a mechanism for its formation been verified. The mechanisms that have been proposed for layer formation can be categorized as either extrinsic or intrinsic to the alloy system. The extrinsic mechanisms rely on externally induced perturbations to the system for layer formation, such as temperature oscillations, growth velocity variations, or vibrations. The intrinsic mechanisms approach layer formation as an alternative type of two phase growth that is inherent for certain peritectic systems and solidification conditions. Convective mixing of the liquid is an additional variable which can strongly influence the development and appearance of layers due to the requisite slow growth rate. The first quantitative description of layer formation is a model recently developed by Trivedi based on the intrinsic mechanism of cyclic accumulation and depiction of solute in the liquid ahead of the interface, linked to repeated nucleation events in the absence of convection. The objective of this research is to characterize the layered microstructures developed during ground-based experiments in which external influences have been minimized as much as possible and to compare these results to the current the model. Also, the differences between intrinsic and externally influenced layer formation were explored. The choice of alloy system is critical to a study of the formation of layered microstructures. The ideal system would have a well-characterized phase diagram, equal densities of both elements in the liquid state to minimize compositionally-driven convective flows, a low peritectic temperature to simplify directional solidification and the achievement of a high temperature gradient in the liquid, a broad composition range for the peritectic reaction, and a reasonable hardness at room temperature to facilitate handling and metallographic preparation. The In-Sn system was selected initially due to a very low peritectic temperature and the nearly equal densities of In and Sn in the liquid state. Since the In-rich peritectic reaction had apparently not been utilized previously for solidification research, experiments were conducted to check the phase diagram in the region of interest. The alloys in this system proved to be difficult to handle and prepare in bulk form with the equipment available, so experiments were initiated with the Sn-Cd system. Layered microstructures had been observed previously in Sn-Cd.

Lograsso, Thomas A.↗

Modeling Cyclic Phase Change and Energy Storage in Solar Heat Receivers

Numerical results pertaining to cyclic melting and freezing of an encapsulated phase change material (PCM), integrated into a solar heat receiver, have been reported. The cyclic nature of the present melt/freeze problem is relevant to latent heat thermal energy storage (LHTES) systems used to power solar Brayton engines in microgravity environments. Specifically, a physical and numerical model of the solar heat receiver component of NASA Lewis Research Center's Ground Test Demonstration (GTD) project was developed and results compared with available experimental data. Multi-conjugate effects such as the convective fluid flow of a low-Prandtl-number fluid, coupled with thermal conduction in the phase change material, containment tube and working fluid conduit were accounted for in the model. A single-band thermal radiation model was also included to quantify reradiative energy exchange inside the receiver and losses through the aperture. The eutectic LiF-CaF2 was used as the phase change material (PCM) and a mixture of He/Xe was used as the working fluid coolant. A modified version of the computer code HOTTube was used to generate results for comparisons with GTD data for both the subcooled and two-phase regimes. While qualitative trends were in close agreement for the balanced orbit modes, excellent quantitative agreement was observed for steady-state modes.

Hall, Carsie A., III↗

Parametric Analysis of Cyclic Phase Change and Energy Storage in Solar Heat Receivers

A parametric study on cyclic melting and freezing of an encapsulated phase change material (PCM), integrated into a solar heat receiver, has been performed. The cyclic nature of the present melt/freeze problem is relevant to latent heat thermal energy storage (LHTES) systems used to power solar Brayton engines in microgravity environments. Specifically, a physical and numerical model of the solar heat receiver component of NASA Lewis Research Center's Ground Test Demonstration (GTD) project was developed. Multi-conjugate effects such as the convective fluid flow of a low-Prandtl-number fluid, coupled with thermal conduction in the phase change material, containment tube and working fluid conduit were accounted for in the model. A single-band thermal radiation model was also included to quantify reradiative energy exchange inside the receiver and losses through the aperture. The eutectic LiF-CaF2 was used as the phase change material (PCM) and a mixture of He/Xe was used as the working fluid coolant. A modified version of the computer code HOTTube was used to generate results in the two-phase regime. Results indicate that parametric changes in receiver gas inlet temperature and receiver heat input effects higher sensitivity to changes in receiver gas exit temperatures.

Hall, Carsie A., III↗

Selection and Assessment of Predictions of the Mars Pathfinder Landing Site

The successful landing of the Mars Pathfinder spacecraft on Mars allows the review of the process of selecting the landing site and assessing predictions made for the site based on Viking and Earth-based data. Selection of the landing site for Mars Pathfinder was a two-phase process. The first phase took place from October 1993 to June 1994 and involved: initial identification of engineering constraints, definition of environmental conditions at the site for spacecraft design, and evaluation of the scientific potential of different landing sites. This phase culminated with the first "Mars Pathfinder Landing Site Workshop", held at the Lunar and Planetary Institute in Houston, Texas on April 18-19, 1994, in which suggested approaches and landing sites were solicited from the entire scientific community. A preliminary site was selected by the project for design purposes in June 1994. The second phase took place from July 1994 to March 1996 and involved: developing criteria for evaluating site safety using images and remote sensing data, testing of the spacecraft and landing subsystems (with design improvements) to establish quantitative engineering constraints on landing site characteristics, evaluating all potential landing sites on Mars, and certification of the site by the project. This phase included a second open workshop, "Mars Pathfinder Landing Site Workshop II: Characteristics of the Ares Vallis Region and Field Trips in the Channeled Scabland, Washington" held in Spokane and Moses Lake September 24-30, 1995 and formal acceptance of the site by NASA Headquarters. Engineering constraints on Pathfinder landing sites were developed from the initial design of the spacecraft and the entry, descent and landing scenario. The site must be within 5 degrees of the subsolar latitude at the time of landing (15N for maximum solar power and flexible communications with Earth. It also must be below 0 km elevation to enable enough time for the parachute to bring the lander to the proper terminal velocity for landing. The entire landing ellipse, which is 70 km by 200 km due to navigational, ephemeris and atmospheric uncertainties, must be free of steep slopes, scarps and obvious hazards in Viking orbiter images, have acceptable radar reflectivity, moderate rock abundances and have little or no dust. Scientific considerations of the Mars Pathfinder payload and mission indicate that analyses of "grab bag" samples at the mouths of outflow channels can offer a first order assessment of a variety of rock types on Mars. Highland sites offer the advantage of in situ analysis of ancient rocks on Mars that record crustal differentiation and the nature of the early environment. Dark gray sites offer the potential of analyzing unweathered and unoxidized materials. Following a general assessment of the safety of different sites, a preliminary selection of a "grab bag" site was made. This site, Ares Vallis, is near the mouth of an outflow channel that may contain ancient Noachian terrain, Hesperian ridged plains, and reworked channel materials. All potential landing sites on Mars that met basic safety criteria were analyzed in detail. Sites (100 by 200 km target ellipses) were considered safe if they were below 0 km elevation, were free of obvious hazards (high relief surface features) in high-resolution (< 50 m/pixel) Viking orbiter images and had acceptable reflectivity and roughness at radar wavelengths, high thermal inertia, moderate rock abundance, low red to violet ratio, and low albedo. Only 4 sites on Mars met all the above criteria, which included 1995 opposition 3.5 cm delay-Doppler radar data. Complete data were evaluated for 7 sites and the Viking landing sites for comparison for all the above criteria as well as crater abundance, hill and mesa abundance, slopes over meter to kilometer scales, low altitude winds (from global circulation models and slopes), the size-frequency distribution of large rocks, as well as rover trafficability and science potential. Discussion of potential hazards at Ares Vallis using a variety of data sets (including radar) at a second open workshop, indicated this site cannot be shown to be any more hazardous than the Viking landing sites. Field trips to the Channeled Scabland and the Ephrata Fan, analogs for Ares Vallis and the landing site, respectively, provided valuable insight into possible geologic processes and potential surface characteristics. Three sites met all the data requirements and safety criteria for landing Pathfinder. Ares Vallis was selected by the project because it appeared acceptably safe (although it appeared to have greater rock abundances than other sites, its elevation was likely the best known) and offered the prospect of analyzing a variety of rock types expected to be deposited by catastrophic floods, which would enable addressing first-order scientific questions such as differentiation of the crust, the development of weathering products, and the nature of the early martian environment and its subsequent evolution. The selection was reviewed by an external board at a number of meetings and accepted, and the site was approved by NASA Headquarters. Data gathered by the Pathfinder lander' and rover provides the opportunity to test the predictions made for the site in the selection process based on remote observations from Earth, orbit, and the surface. The discussion below is taken from Golombek et al. to which the reader is referred for a more complete discussion and a complete list of references, which are omitted here for brevity. Many characteristics of the landing site are consistent with its being shaped and deposited by the Ares and Tiu catastrophic floods. The rocky surface is consistent a depositional plain comprising semi-rounded pebbles, cobbles and tabular boulders (some of which appear imbricated and/or inclined in the direction of flow) that appear similar to depositional plains in terrestrial catastrophic floods. The Twin Peaks appear to be streamlined hills in lander images, which is consistent with interpretations of larger hills in Viking orbiter images of the region that suggest the lander is on the flank of a broad, gentle ridge trending northeast from Twin Peaks. This ridge, which is the rise to the north of the lander, is aligned in the downstream direction from the Ares and Tiu Valles floods, and may be a debris tail deposited in the wake of the Twin Peaks. Channels visible throughout the scene may be a result of late stage drainage. As predicted by delay-Doppler radar measurements and tracking results, the average elevation of the center of the site was about the same as Viking Lander I relative to the 6.1 mbar geoid. The Doppler tracking and two-way ranging estimate for the elevation of the spacecraft is only 45 in lower than the Viking I Lander and within 100 in of that expected, which is within the uncertainties of the measurements. After landing, surface pressures and winds (5-10 m/s) were found to be similar to expectations based on Viking data, although temperatures were about 10 K warmer. The temperature profile below 50 km was also roughly 20 K warmer. As a result, predicted densities were 5% higher near the surface and up to 40% lower at 50 km but within the entry, descent and landing design margins. The populations of craters and small hills and the slopes of the hills measured in high-resolution (38 m/pixel) Viking orbiter images and the radar derived slopes of the landing site are all consistent with observations of these properties in the lander images. A rocky surface was expected from Viking Infra-Red Thermal Mapper (IRTM) observations and comparisons with the Viking landing sites. The observed cumulative fraction of area covered by rocks with diameters greater than 3 cm and heights greater than 0.5 in (potentially hazardous to landing) at Ares is similar to that predicted by IRTM observations and models of Viking lander and Earth analog rock size-frequency distributions. The IRTM prediction postulated an effective thermal inertia of 30 (10(exp -3) cgs units - cal/cubic cm/s(exp 0.5)/K) for the rock population, but we obtain a slightly different effective thermal inertia for the actual rock population. The validity of interpretations of radar echoes prior to landing are supported by a simple radar echo model, an estimate of the reflectivity of the soil from its bulk density, and the fraction of area covered by rocks. In the calculations, the soil produces the quasi-specular echo and the rocks produce the diffuse echo. The derived quasispecular cross section is comparable to the cross-sections and reflectivities reported for 3.5-cm wavelength observations. The model yields a diffuse echo that is modestly larger than the polarized diffuse echo reported for 3.5-cm wavelength observations. At 12.5-cm wavelength, similar rock populations at Ares and the Viking I site were expected because the diffuse echoes are comparable, but the large normal reflectivities suggests that bulk densities of the soils at depth are greater than those at the surface. We also obtain a fine-component inertia near 8.4 which agrees with the fine-component inertia of 8.7 (in 10(exp -3) cgs units) estimated from thermal observations from orbit by the IRTM; for this estimate, we used a bulk thermal inertia of 10.4 for the landing site, an effective thermal inertia near 40 (10(exp -3) cgs units) for the rock population, and a graphical representation of Kieffer's model. Color and albedo data for Ares suggested surfaces of materials at Ares Vallis would be relatively dust free or unweathered prior to landing compared with the materials at the Viking landing sites. This suggestion is supported by the abundance of relatively dark-gray rocks at Ares and their relative rarity at the Viking landing sites, where rocks are commonly coated with bright red dust. Finally, the 40 km long Ephrata Fan of the Channeled Scabland in Washington state, which was deposited where c

Golombek, M. P.↗

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