ALBERT

Description: We have investigated the interaction of Io, Jupiter's innermost Galilean satellite, with the Io plasma torus. The interaction of Io with the plasma surrounding it has been a subject of interest for almost 30 years, dating from the discovery by Bigg (1964) that radio emissions from the Jovian magnetosphere are controlled by Io's position. Since that time, both ground-based and spacecraft observations have shown that Io is a unique satellite that influences the Jovian magnetosphere in important ways. In particular, material from Io is a major source of plasma for the magnetosphere, and the energy that this plasma harnesses from Jupiter's co-rotating magnetic field is an important power source for the magnetosphere. It is apparent that the local interaction of the torus plasma with Io plays a key role in the formation, composition, and energetics of the Io torus; the interaction is also highly nonlinear. We have modeled this interaction using time-dependent three-dimensional magnetohydrodynamic (MHD) simulations. During this past year, we have used NASA support to develop a new MHD code to study the interaction. As part of the Galileo spacecraft's recent successful insertion into orbit around Jupiter, the spacecraft passed within 900 km of Io's surface. Our calculations have focused on using Galileo particles and fields data to examine a question that was not resolved by the Voyager observations: Does Io have an intrinsic magnetic field? In this progress summary, we describe our efforts on this problem to date.

Description: The NPARC Alliance is a partnership between the NASA Lewis Research Center (LeRC) and the USAF Arnold Engineering Development Center (AEDC) dedicated to the establishment of a national CFD capability, centered on the NPARC Navier-Stokes computer program. The three main tasks of the Alliance are user support, code development, and validation. The present paper is a status report on the validation effort. It describes the validation approach being taken by the Alliance. Representative results are presented for laminar and turbulent flat plate boundary layers, a supersonic axisymmetric jet, and a glancing shock/turbulent boundary layer interaction. Cases scheduled to be run in the future are also listed. The archive of validation cases is described, including information on how to access it via the Internet.

Description: We study the onset of a pure Marangoni convection in a liquid layer with two deformable interfaces in the no-gravity environment. Both oscillatory and stationary instabilities are considered for a wide range of parameters. It is shown that only stationary instability is possible when surface tension at the colder interface is lower than that at the hotter one. Oscillatory instability tends to disappear and to be replaced by the stationary instability with increase of the Prandtl number and decrease of surface tension at the colder interface.

Description: Estimates of the mass of dust suspended in the Martian atmosphere are derived from global and regional 9-micrometer opacity maps produced from Viking Infrared Thermal Mapper data. During the peak of the 1977b storm, a total dust mass of approximately 4.3 x 10(exp 14) g was suspended, equivalent to 4.3 x 10(exp -4) g/sq cm, or a layer 1.4 micrometers thick. During a local dust storm near Solis Planum at L(sub s) 227 deg, approximately 1.3 x 10(exp 13) g of dust were lofted, equal to about a 6-micrometer layer in that vicinity.

Description: We model an infrared outburst on Io as being due to a large, erupting lava flow which increased its area at a rate of 1.5 x 10(exp 5)/sq m and cooled from 1225 to 555 K over the 2.583-hr period of observation. The inferred effusion rate of 3 x 10(exp 5) cu m/sec for this eruption is very high, but is not unprece- dented on the Earth and is similar to the high eruption rates suggested for early lunar volcanism. Eruptions occur approxi- mately 6% of the time on Io. These eruptions provide ample resurfacing to explain Io's lack of impact craters. We suggest that the large total radiometric heat flow, 10(exp 14) W, and the size and temperature distribution of the thermal anomalies (McEwen et al. 1992; Veeder et al. 1994) can be accounted for by a series of silicate lava flows in various stages of cooling. We propose that the whole suite of Io's currently observed thermal anomalies was produced by multiple, high-eruptive-rate silicate flows within the past century.

Description: A consistent solution of the radiative transfer equation characterizing photon transport in a semi-infinite medium of refractive index greater than or equal to one is obtained following the method of Sobolev. Fresnel specular reflection, Snell's law and isotropic scattering are assumed. An algorithm is developed and its accuracy is demonstrated. A numerical Laplace transform inversion leads to an efficient evaluation for the interior flux and source function distributions.

Description: We investigate the orbital dynamics of small dust particles generated via the continuous micrometeoroid bombardment of the Martian moons. In addition to Mar's oblateness, we also consider the radiation pressure perturbation that is complicated by the planet's eccentric orbit and tilted rotational axis. Considering the production rates and the lifetimes of dust grains, we show that particles from Deimos with radii of about 15 micrometers are expected to dominate the population of a permanently present and tilted dust torus. This torus has an estimated peak number density of approximately equals 5 x 10(exp -12)/cu cm and an optical depth of approximately equals 4 x 10(exp -8).

Description: The valley network channels on the heavily cratered ancient surface of Mars suggest the presence of liquid water approximately 3.8 Gyr ago. However, the implied warm climate is difficult to explain in the context of the standard solar model, even allowing for the maximum CO2 greenhouse heating. In this paper we investigate the astronomical and planetary implications of a nonstandard solar model in which the zero-age, main-sequence Sun had a mass of 1.05 +/- 0.02 Solar Mass. The excess mass was subsequently lost in a solar wind during the first 1.2(-0.2, +0.4) Gyr of the Sun's main sequence phase. The implied mass-loss rate of 4(+3, -2) x 10(exp -11) M/yr, or about 10(exp 3) x that of the current Sun, may be detectable in several nearby young solar type stars.

Description: Hellas basin on Mars has been the site of volcanism, tectonism, and modification by fluvial, mass-wasting, and eolian processes over its more than 4-b.y. existence. Our detailed geologic mapping and related studies have resulted in the following new interpretations. The asymmetric distribution of highland massifs and other structures that define the uplifted basin rim suggest a formation of the basin by the impact of a low-angle bolide having a trajectory heading S60E. During the Late Noachian, the basin was infilled, perhaps by lava flows, that were sufficiently thick (greater than 1 km) to produce wrinkle ridges on the fill material and extensional faulting along the west rim of the basin. At about the same time, deposits buried northern Malea Planum, which are interpreted to be pyroclastic flows from Amphitrites and Peneus Paterae on the basis of their degraded morphology, topology, and the application of a previous model for pyroclastic volcanism on Mars. Peneus forms a distinctive caldera structure that indicates eruption of massive volumes of magma, whereas Amphitrites is a less distinct circular feature surrounded by a broad, low, dissected shield that suggests generally smaller volume eruptions. During the Early Hesperian, an approximately 1-to 2km-thick sequence of primarily fined-grained, eolian material was deposited on the floor of Hellas basin. Subsequently, the deposit was deeply eroded, except where armored by crater ejecta, and it retreated as much as 200-300 km along its western margin, leaving behind pedestal craters and knobby outliers of the deposit. Local debris flows within the deposit attest to concentrations of groundwater, perhaps in part brought in by outflow floods along the east rim of the basin. These floods may have deposited approximately 100-200m of sediment, subduing wrinkle ridges in the eastern part of the basin floor. During the Late Hesperian and Amazonian, eolian mantles were emplaced on the basin rim and floor and surrounding highlands. Their subsequent erosion resulted in pitted and etched plains and crater fill, irregular mesas, and pedestal craters. Local evidence occurs for the possible former presence of ground ice or ice sheets approximately 100 km across; however, we disagree with a hypothesis that suggest that the entire south rim and much of the floor of Hellas have been glaciated. Orientations of dune fields and yardangs in lower parts of Hellas basin follow directions of the strongest winds predicted by a recently published general circulation model (GCM). Transient frost and dust splotches in the region are, by contrast, related to the GCM prediction for the season in which the images they appear in were taken.

Description: Visible and near-IR refectivity, Moessbauer, and X ray diffraction data were obtained on powders of impact melt rock from the Manicouagan Impact Crater located in Quebec, Canada. The iron mineralogy is dominated by pyroxene for the least oxidized samples and by hematite for the most oxidized samples. Phyllosilicate (smectite) contents up to approximately 15 wt % were found in some heavily oxidized samples. Nanophase hematite and/or paramagnetic ferric iron is observed in all samples. No hydrous ferric oxides (e.g., goethite, lepidocrocite, and ferrihydrite) were detected, which implies the alteration occurred above 250 C. Oxidative alteration is thought to have occurred predominantly during late-stage crystallization and subsolidus cooling of the impact melt by invasion of oxidizing vapors and/or solutions while the impact melt rocks were still hot. The near-IR band minimum correlated with the extent of aleration Fe(3+)/Fe(sub tot) and ranged from approximately 1000 nm (high-Ca pyroxene) to approximately 850 nm (bulk, well-crystalline hematite) for least and most oxidized samples, respectively. Intermediate band positions (900-920 nm) are attributed to low-Ca pyroxene and/or a composite band from hematite-pyroxene assemblages. Manicouagan data are consistent with previous assignments of hematite and pyroxene to the approximately 850 and approximately 1000nm bands observed in Martian reflectivity spectra. Manicouagan data also show that possible assignments for intermediate band positions (900-920 nm) in Martian spectra are pyroxene and/or hematite-pyroxene assemblages. By analogy with impact melt sheets and in agreement with observables for Mars, oxidative alteration of Martian impact melt sheets above 250 C and subsequent erosion could produce rocks and soils with variable proportions of hematite (both bulk and nanophase), pyroxene, and phyllosilicates as iron-bearing mineralogies. If this process is dominant, these phases on Mars were formed rapidly at relativly high temperatures on a sporadic basis throughout the history of the planet. The Manicouagan samples also show that this mineralogical diversity can be accomplished at constant chemical composition, which is also indicated for Mars from the analyses of soil at the two Viking landing sites.

Description: Results are reported of the Surface Tension Driven Convection Experiment (STDCE) aboard USML-1 Spacelab. Steady and transient thermocapillary flows were investigated in a 10 cm dia. circular container filled with 10 Cs silicone oil. The velocity and temperature fields were studied in detail under various conditions. It is shown in this paper how the Marangoni number affects the velocity field. A numerical analysis of the flows was also conducted and its results were compared to the experimental data. Good agreement is shown.

Description: Based on the conservation of chemical elements in chemical reactions, a rule is proved that the number of boundary conditions given by densities and/or nonzero velocities should not be less than the number of chemical elements in the system, and the boundary conditions for species given by densities and velocities should include all elements in the system. Applications of this rule to Mars are considered. It is shown that the problem of the CO2-H2O chemistry in the lower and middle atmosphere of Mars, say, in the range of 0-80 km does not have a unique solution, if only CO2 and H2O densities are given at the lower boundary, and the remaining boundary conditions are fluxes. Two examples of models of this type are discussed. Two models of the photochemistry of the Martian atmosphere, with and without nitrogen chemistry, are considered. The oxygen nonthermal escape ratio of 1.2 x 10(exp 8)/cu cm/s is given at 240 km and is balanced with the total hydrogen escape rate within an uncertainty of 1% for both models. Both models fit the measured O2 and CO mixing ratios, the O3 abundance, and the O2 1.27-micrometer dayglow almost within the uncertainties of the measured values, though the model without nitrogen chemistry fits better. The importance of nitrogen chemistry in the lower and middle atmosphere of Mars depends on a fine balance between production of NO and N in the upper atmosphere which is not known within the required accuracy.

Description: Goethite-bearing samples with values of Mn(s) (Mn/(Mn+Fe) mole fraction) up to 0.206 were synthesized by precipitation from alkaline solution. Samples with Mn(s) less than or equal 0.061 were single-phase Mn-goethites: samples with higher Mn(s) values contained another Mn-bearing phase (probably jacobsite). Mn-hematites were prepared by dehydroxylation of corresponding Mn-goethites at 500 C. Orthorhombic a and b unit cell dimensions of Mn-goethites changed in a linear manner with Mn(s), but not at rates predicted by the Vegrad law. Hexagonal unit cell dimensions of Mn-hematites did not vary with Mn(s). Moessbauer parameters isomer shift (IS), quadrupole splitting (QS), and hyperfine field (B(sub hf)) were measured at 293 and 15 K. For all single-phase Mn-goethites and Mn-hematites (Mn(s) less than or equal 0.061), magnetic splitting was observed at both temperatures. At 293 K, small but systematic decreases in B(sub hf) were observed with increasing Mn substitution; IS and QS were not dependent on Mn(s). Mn substitution strongly lowered the Morin transition temperature of hematite. At 15 K, the Morin transition was not present for Mn(s) greater than 0.020(4). The saturation magnetization of Mn-goethites and Mn-hematites (Mn(s) less than or equal 0.061) was the expected zero (within error) for antiferromagnetic goethite and for hematites obtained from dehydroxylation of goethites. Mn-geothites with Mn(s) greater than 0.061 were magnetic because of the presence of strongly magnetic jacobsite. For reflectivity spectra, bands resulting from MN(3+) were centered near 454 and 596 nm for Mn-goethites and near 545 and 700 nm for Mn-hematites. There is evidence for a approximately 700 nm band in spectral data for Martian bright regions, but association of it with Mn(3+) is not a unique interpretation. Comparison of laboratory and Martian spectral data implies that Mn(s) less than 0.032 for the Mn(3+) content of Martian hematites.

Description: Ferric-iron-bearing materials play an important role in the interpretation of visible to near-IR Mars spectra, and they may play a similarly important role in the analysis of new mid-IR spacecraft spectral observations to be obtained over the next decade. We review exisiting data on mid-IR transmission spectra of ferric oxides/oxyhydroxides and present new transmission spectra for ferric-bearing materials spanning a wide range of mineralogy and crystallinity. These materials include 11 samples of well-crystallized ferric oxides (hematite, maghemite, and magnetite) and ferric oxyhydroxides (goethite, lepidocrocite). We also report the first transmission spectra for purely nanophase ferric oxide samples that have been shown to exhibit spectral similarities to Mars in the visible to near-IR and we compare these data to previous and new transmission spectra of terrestial palagonites. Most of these samples show numerous, diagnostic absorption features in the mid-IR due to Fe(3+) - 0(2-) vibrational transitions, structural and/or bound OH, and/or silicates. These data indicate that high spatial resolution, moderate spectral resolution mid-IR ground-based and spacecraft observations of Mars may be able to detect and uniquely discriminate among different ferric-iron-bearing phases on the Martian surface or in the airborne dust.

Description: Eleven tests were carried out in DLR's high enthalpy tunnel in Goettingen (HEG), Germany, for reservoir conditions ranging from 10 to 23 MJ/kg and a free stream Mach number of approximately 10. A blunted cone model with a cylindrical afterbody (sting) was investigated. To obtain information on the influence of defined parameters on the body back flow, especially of the reacting gas, the heat transfer rate along the body contour was measured with fast response surface thermocouples on the forebody and sensitive thin film heat transfer gauges on the base and sting of the model. Flow visualization with a holographic interferometry system was provided. Tests are described, and a preliminary interpretation of the observed flow effects are given.

Description: Get Away Special (GAS) payload G-093, also called VORTEX (VOrtex Ring Transit EXperiment), is an investigation of the propagation of a vortex ring through a liquid-gas interface in microgravity. This process results in the formation of one or more liquid droplets similar to earth based liquid atomization systems. In the absence of gravity, surface tension effects dominate the drop formation process. The Shuttle's microgravity environment allows the study of the same fluid atomization processes as using a larger drop size than is possible on Earth. This enables detailed experimental studies of the complex flow processes encountered in liquid atomization systems. With VORTEX, deformations in both the vortex ring and the fluid surface will be measured closely for the first time in a parameters range that accurately resembles liquid atomization. The experimental apparatus will record images of the interactions for analysis after the payload has been returned to earth. The current design of the VORTEX payload consists of a fluid test cell with a vortex ring generator, digital imaging system, laser illumination system, computer based controller, batteries for payload power, and an array of housekeeping and payload monitoring sensors. It is a self-contained experiment and will be flown on board the Space Shuttle in a 5 cubic feet GAS canister. The VORTEX Project is entirely run by students at the University of Michigan but is overseen by a faculty advisor acting as the payload customer and the contact person with NASA. This paper summarizes both the technical and programmatic aspects of the VORTEX Project.

Description: This paper focuses on the flight results of the Cryogenic Two-Phase Flight Experiment (CRYOTP), which was a Hitchhiker based experiment that flew on the space shuttle Columbia in March of 1994 (STS-62). CRYOTP tested two new technologies for advanced cryogenic thermal control; the Space Heat Pipe (SHP), which was a constant conductance cryogenic heat pipe, and the Brilliant Eyes Thermal Storage Unit (BETSU), which was a cryogenic phase-change thermal storage device. These two devices were tested independently during the mission. Analysis of the flight data indicated that the SHP was unable to start in either of two attempts, for reasons related to the fluid charge, parasitic heat leaks, and cryocooler capacity. The BETSU test article was successfully operated with more than 250 hours of on-orbit testing including several cooldown cycles and 56 freeze/thaw cycles. Some degradation was observed with the five tactical cryocoolers used as thermal sinks, and one of the cryocoolers failed completely after 331 hours of operation. Post-flight analysis indicated that this problem was most likely due to failure of an electrical controller internal to the unit.

Description: The Capillary Pumped Loop Flight Experiment (CAPL) employs a passive two-phase thermal control system that uses the latent heat of vaporization of ammonia to transfer heat over long distances. CAPL was designed as a prototype of the Earth Observing System (EOS) instrument thermal control systems. The purpose of the mission was to provide validation of the system performance in micro-gravity, prior to implementation on EOS. CAPL was flown on STS-60 in February, 1994, with some unexpected results related to gravitational effects on two-phase systems. Flight test results and post flight investigations will be addressed, along with a brief description of the experiment design.

Description: This paper presents an overview of the Cryogenic Test Bed (CTB) experiments including experiment results, integration techniques used, and lessons learned during integration, test and flight phases of the Cryogenic Heat Pipe Flight Experiment (STS-53) and the Cryogenic Two Phase Flight Experiment (OAST-2, STS-62). We will also discuss the Cryogenic Flexible Diode Heat Pipe (CRYOFD) experiment which will fly in the 1996/97 time frame and the fourth flight of the CTB which will fly in the 1997/98 time frame. The two missions tested two oxygen axially grooved heat pipes, a nitrogen fibrous wick heat pipe and a 2-methylpentane phase change material thermal storage unit. Techniques were found for solving problems with vibration from the cryo-collers transmitted through the compressors and the cold heads, and mounting the heat pipe without introducing parasitic heat leaks. A thermally conductive interface material was selected that would meet the requirements and perform over the temperature range of 55 to 300 K. Problems are discussed with the bi-metallic thermostats used for heater circuit protection and the S-Glass suspension straps originally used to secure the BETSU PCM in the CRYOTP mission. Flight results will be compared to 1-g test results and differences will be discussed.

Description: The Voyager 1 spacecraft flew by Jupiter on March 5, 1979. Spacecraft navigation was performed with radio tracking data from NASA's Deep Space Network. In the years since then, there has been a great deal of progress in the definition of celestial reference frames and in determining the orbit and orientation of the Earth. Using these improvements, the radio metric range and Doppler data acquired from the Voyager 1 spacecraft near its encounter with Jupiter have been reanalyzed to determine the plane-of-sky position of Jupiter with much greater accuracy than was possible at the time of the encounter. The position of Jupiter at the time of encounter has been determined with an accuracy of 40 nrad in right ascension and 140 nrad in declination with respect to the celestial reference frame defined by the International Earth Rotation Service. This position estimate has been done to improve the ephemeris of Jupiter prior to the upcoming encounter of the Galileo spacecraft with Jupiter.

Description: With the recent focus on the needs of design and applications CFD, research groups have begun to address the traditional bottlenecks of grid generation and surface modeling. Now, a host of emerging technologies promise to shortcut or dramatically simplify the simulation process. This paper discusses the current status of these emerging technologies. It will argue that some tools are already available which can have positive impact on portions of the design cycle. However, in most cases, these tools need to be integrated into specific engineering systems and process cycles to be used effectively. The rapidly maturing status of unstructured and Cartesian approaches for inviscid simulations makes suggests the possibility of highly automated Euler-boundary layer simulations with application to loads estimation and even preliminary design. Similarly, technology is available to link block structured mesh generation algorithms with topology libraries to avoid tedious re-meshing of topologically similar configurations. Work in algorithmic based auto-blocking suggests that domain decomposition and point placement operations in multi-block mesh generation may be properly posed as problems in Computational Geometry, and following this approach may lead to robust algorithmic processes for automatic mesh generation.

Description: The last decade has witnessed a vigorous and sustained research effort on unstructured methods for computational fluid dynamics. Unstructured mesh generators and flow solvers have evolved to the point where they are now in use for design purposes throughout the aerospace industry. In this paper we survey the various mesh types, structured as well as unstructured, and examine their relative strengths and weaknesses. We argue that unstructured methodology does offer the best prospect for the next generation of computational fluid dynamics algorithms.

Description: A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations in two dimensions is developed and tested. Grids about geometrically complicated bodies are generated automatically, by recursive subdivision of a single Cartesian cell encompassing the entire flow domain. Where the resulting cells intersect bodies, N-sided 'cut' cells are created using polygon-clipping algorithms. The grid is stored in a binary-tree data structure which provides a natural means of obtaining cell-to-cell connectivity and of carrying out solution-adaptive mesh refinement. The Euler and Navier-Stokes equations are solved on the resulting grids using a finite-volume formulation. The convective terms are upwinded: A gradient-limited, linear reconstruction of the primitive variables is performed, providing input states to an approximate Riemann solver for computing the fluxes between neighboring cells. The more robust of a series of viscous flux functions is used to provide the viscous fluxes at the cell interfaces. Adaptively-refined solutions of the Navier-Stokes equations using the Cartesian, cell-based approach are obtained and compared to theory, experiment and other accepted computational results for a series of low and moderate Reynolds number flows.

Description: Grid related issues of the Chimera overset grid method are discussed in the context of a method of solution and analysis of unsteady three-dimensional viscous flows. The state of maturity of the various pieces of support software required to use the approach is considered. Current limitations of the approach are identified.

Description: This paper presents 'A view from the trenches' on CFD grid generation from a Pratt & Whitney perspective. We anticipate that other organizations have similar views. We focus on customer expectations and the consequent requirements. We enunciate a vision for grid generation, discuss issues that developers must recognize.

Description: This paper presents a perspective on the requirements that Computational Fluid Dynamics (CFD) technology must meet for its effective use in aerospace design. General observations are made on current aerospace design practices and deficiencies are noted that must be rectified for the U.S. aerospace industry to maintain its leadership position in the global marketplace. In order to rectify deficiencies, industry is transitioning to an integrated product and process development (IPPD) environment and design processes are undergoing radical changes. The role of CFD in producing data that design teams need to support flight vehicle development is briefly discussed. An overview of the current state of the art in CFD is given to provide an assessment of strengths and weaknesses of the variety of methods currently available, or under development, to produce aerodynamic data. Effectiveness requirements are examined from a customer/supplier view point with design team as customer and CFD practitioner as supplier. Partnership between the design team and CFD team is identified as an essential requirement for effective use of CFD. Rapid turnaround, reliable accuracy, and affordability are offered as three key requirements that CFD community must address if CFD is to play its rightful role in supporting the IPPD design environment needed to produce high quality yet affordable designs.

Description: The efforts in geometry modeling and grid generation at the NASA Lewis Research Center, as applied to the computational fluid dynamic (CFD) analysis of aeropropulsion systems, are presented. The efforts are mainly characterized by a focus on the analysis of components of an aeropropulsion system, which involve turbulent viscous flow with heat transfer and chemistry. Thus, this discussion will follow that characterization and will sequence through the components of typical propulsion systems consisting of inlets, compressors, combustors, turbines, and nozzles. For each component, some applications of CFD analysis will be presented to show how CFD is used to compute the desired performance information, how geometry modeling and grid generation are performed, and what issues have developed related to geometry modeling and grid generation. The discussion will illustrate the following needs related to geometry modeling and grid generation as observed in aeropropulsion analysis: (1) accurate and efficient resolution of turbulent viscous and chemically-reacting flowfields; (2) easy-to-use interfaces with CAD data for automated grid generation about complex geometries; and (3) automated batch grid generation software for use with design and optimization software.

Description: This viewgraph presentation discusses key features of current combustion system CFD modeling capabilities at GE Aircraft Engines provided by the CONCERT code; CONCERT development history; modeling applied for designing engine combustion systems; modeling applied to improve fundamental understanding; CONCERT3D results for current production combustors; CONCERT3D model of NASA/GE E3 combustor; HYBRID CONCERT CFD/Monte-Carlo modeling approach; and future modeling directions.

Description: A procedure was developed to improve the turn-around time for computational fluid dynamics (CFD) simulations of an inlet-bleed problem involving oblique shock-wave/boundary-layer interactions on a flat plate with bleed into a plenum through one or more circular holes. This procedure is embodied in a preprocessor called AUTOMAT. With AUTOMAT, once data for the geometry and flow conditions have been specified (either interactively or via a namelist), it will automatically generate all input files needed to perform a three-dimensional Navier-Stokes simulation of the prescribed inlet-bleed problem by using the PEGASUS and OVERFLOW codes. The input files automatically generated by AUTOMAT include those for the grid system and those for the initial and boundary conditions. The grid systems automatically generated by AUTOMAT are multi-block structured grids of the overlapping type. Results obtained by using AUTOMAT are presented to illustrate its capability.

Description: The role of Computational Fluid Dynamics (CFD) at Ames Research Center has expanded to address a broad range of aeronautical problems, including wind tunnel support, flight test support, design, and analysis. Balancing the requirements of each new problem against the available resources - software, hardware, time, and expertise - is critical to the effective use of CFD. Several case studies of recent applications highlight the depth of CFD capability at Ames, the tradeoffs involved in various approaches, and lessons learned in the use of CFD as an engineering tool.

Description: When computing the flow around complex three dimensional configurations, the generation of the mesh is the most time consuming part of any calculation. With some meshing technologies this can take of the order of a man month or more. The requirement for a number of design iterations coupled with ever decreasing time allocated for design leads to the need for a significant acceleration of this process. Of the two competing approaches, block-structured and unstructured, only the unstructured approach will allow fully automatic mesh generation directly from a CAD model. Using this approach coupled with the techniques described in this paper, it is possible to reduce the mesh generation time from man months to a few hours on a workstation. The desire to closely couple a CFD code with a design or optimization algorithm requires that the changes to the geometry be performed quickly and in a smooth manner. This need for smoothness necessitates the use of Bezier polynomials in place of the more usual NURBS or cubic splines. A two dimensional Bezier polynomial based design system is described.

Description: The results of a first order perfect gas correction for the effects of the boundary layer formation within expansion tubes with nozzles are presented. The analytical model developed to describe the boundary layer formation within the expansion tube and an expansion nozzle located at the end of the acceleration tube is based on the Karman integral equations. The results of this analytical model are compared with experimental data from an expansion diffuser. The model provides a useful tool for the preliminary design of nozzles for such facilities.

Description: In this paper the capabilities of an automated CFD system which is currently available at NLR are demonstrated. Transonic flow around the AS28G wing/body configuration and hypersonic flow through a generic three-dimensional mixed-compression airbreathing inlet are simulated. An assessment of the level of automation of the current CFD-system is made. The problem-turnaround time lies within the order of a week for both applications.

Description: Detailed simulations of viscous flows in complicated geometries pose a significant challenge to current capabilities of Computational Fluid Dynamics (CFD). To enable routine application of CFD to this class of problems, advanced methodologies are required that employ (a) automated grid generation, (b) adaptivity, (c) accurate discretizations and efficient solvers, and (d) advanced software techniques. Each of these ingredients contributes to increased accuracy, efficiency (in terms of human effort and computer time), and/or reliability of CFD software. In the long run, methodologies employing structured grid systems will remain a viable choice for routine simulation of flows in complex geometries only if genuinely automatic grid generation techniques for structured grids can be developed and if adaptivity is employed more routinely. More research in both these areas is urgently needed.

Description: A discussion of the strengths and weaknesses of overset composite grid and solution technology is given, along with a sampling of current work in the area. Major trends are identified, and the observation is made that generalized and hybridized overset methods provide a natural framework for combining disparate mesh types and physics models. Because of this, the author concludes that overset methods will be the foundation for the general purpose computational fluid dynamics programs of the future.

Description: The status of CFD methods based on the use of block-structured grids for analyzing viscous flows over complex configurations is examined. The objective of the present study is to make a realistic assessment of the usability of such grids for routine computations typically encountered in the aerospace industry. It is recognized at the very outset that the total turnaround time, from the moment the configuration is identified until the computational results have been obtained and postprocessed, is more important than just the computational time. Pertinent examples will be cited to demonstrate the feasibility of solving flow over practical configurations of current interest on block-structured grids.

Description: This viewgraph presentation discusses turbulence modeling requirements, development philosophy, and approach; two major areas of concentration (high speed and low speed turbulence modeling); high speed turbulence modeling; compressibility effects; turbulence models adapted to USA code; M = 9.2 flat plate flow; Mach 7.05 flow over axisymmetric flare; Mach 8.6 flow over cold wall edge; low speed turbulence modeling; turbulence models being assessed; turbulence model deck structure and integration with Navier-Stokes solver; nonlinear algebraic-stress model; rotation modified k-epsilon model; and Reynolds stress model.

Description: This viewgraph presentation covers gas turbine combustor flow physics, turbulence model investigations, turbulent combustion modeling, and present status and future needs.

Description: This viewgraph presentation discusses geometry and flow configuration, effect of y+ on heat transfer computations, standard and extended k-epsilon turbulence model results with wall function, low-Re model results (the Lam-Bremhorst model without wall function), a criterion for flow reversal in a radially rotating square duct, and a summary.

Description: This viewgraph presentation demonstrates that computationally efficient k-l and k-kl turbulence models have been developed and implemented at Lockheed Fort Worth Company. Many years of experience have been gained applying two equation turbulence models to complex three-dimensional flows for design and analysis.

Description: The objective of this viewgraph presentation is to evaluate turbulence models for integrated aircraft components such as the forebody, wing, inlet, diffuser, nozzle, and afterbody. The one-equation models have replaced the algebraic models as the baseline turbulence models. The Spalart-Allmaras one-equation model consistently performs better than the Baldwin-Barth model, particularly in the log-layer and free shear layers. Also, the Sparlart-Allmaras model is not grid dependent like the Baldwin-Barth model. No general turbulence model exists for all engineering applications. The Spalart-Allmaras one-equation model and the Chien k-epsilon models are the preferred turbulence models. Although the two-equation models often better predict the flow field, they may take from two to five times the CPU time. Future directions are in further benchmarking the Menter blended k-w/k-epsilon and algorithmic improvements to reduce CPU time of the two-equation model.

Description: This viewgraph presentation summarizes some CFD experience at GE Aircraft Engines for flows in the primary gaspath of a gas turbine engine and in turbine blade cooling passages. It is concluded that application of the standard k-epsilon turbulence model with wall functions is not adequate for accurate CFD simulation of aerodynamic performance and heat transfer in the primary gas path of a gas turbine engine. New models are required in the near-wall region which include more physics than wall functions. The two-layer modeling approach appears attractive because of its computational complexity. In addition, improved CFD simulation of film cooling and turbine blade internal cooling passages will require anisotropic turbulence models. New turbulence models must be practical in order to have a significant impact on the engine design process. A coordinated turbulence modeling effort between NASA centers would be beneficial to the gas turbine industry.

Description: Program goals at the Center for Modeling of Turbulence and Transition (CMOTT), NASA Lewis Research Center, are (1) to develop reliable turbulence (including bypass transition) and combustion models for complex flows in propulsion systems and (2) to integrate developed models into deliverable CFD tools for propulsion systems in collaboration with industry. This viewgraph presentation covers the following topics: development of turbulence and combustion models; collaboration with industry and technology transfer; isotropic eddy viscosity models; algebraic Reynolds stress models; scalar turbulence models; second order closure models; multiple scale k-epsilon models; and PDF modeling of turbulent reacting flows.

Description: An assessment of two unstructured methods is presented in this paper. A tetrahedral unstructured method USM3D, developed at NASA Langley Research Center is compared to a Cartesian unstructured method, SPLITFLOW, developed at Lockheed Fort Worth Company. USM3D is an upwind finite volume solver that accepts grids generated primarily from the Vgrid grid generator. SPLITFLOW combines an unstructured grid generator with an implicit flow solver in one package. Both methods are exercised on three test cases, a wing, and a wing body, and a fully expanded nozzle. The results for the first two runs are included here and compared to the structured grid method TEAM and to available test data. On each test case, the set up procedure are described, including any difficulties that were encountered. Detailed descriptions of the solvers are not included in this paper.

Description: The past ten years have seen steady progress in surface modeling procedures, and wholesale changes in grid generation technology. Today, it seems fair to state that a satisfactory grid can be developed to model nearly any configuration of interest. The issues at present focus on operational concerns such as cost and quality. Continuing evolution of the engineering process is placing new demands on the technologies of surface modeling and grid generation. In the evolution toward a multidisciplinary analysis-bascd design environment, methods developed for Computational Fluid Dynamics are finding acceptance in many additional applications. These two trends, the normal evolution of the process and a watershed shift toward concurrent and multidisciplinary analysis, will be considered in assessing current capabilities and needed technological improvements.

Description: This viewgraph presentation discusses what models are used in this package and what their advantages and disadvantages are, how the probability density function (PDF) model is implemented and the features of the program, and what can be expected in the future from the NASA Lewis PDF code.

Description: This viewgraph presentation discusses an extension of the probability density function (PDF) method to the modeling of spray flames to evaluate the limitations and capabilities of this method in the modeling of gas-turbine combustor flows. The comparisons show that the general features of the flowfield are correctly predicted by the present solution procedure. The present solution appears to provide a better representation of the temperature field, particularly, in the reverse-velocity zone. The overpredictions in the centerline velocity could be attributed to the following reasons: (1) the use of k-epsilon turbulence model is known to be less precise in highly swirling flows and (2) the swirl number used here is reported to be estimated rather than measured.

Description: Viewgraphs are presented on computation of turbulent combustion, governing equations, closure problem, PDF modeling of turbulent reactive flows, validation cases, current projects, and collaboration with industry and technology transfer.

Description: Viewgraphs are presented on modeling turbulent reacting flows, regimes of turbulent combustion, regimes of premixed and regimes of non-premixed turbulent combustion, chemical closure models, flamelet model, conditional moment closure (CMC), NO(x) emissions from turbulent H2 jet flames, probability density function (PDF), departures from chemical equilibrium, mixing models for PDF methods, comparison of predicted and measured H2O mass fractions in turbulent nonpremixed jet flames, experimental evidence of preferential diffusion in turbulent jet flames, and computation of turbulent reacting flows.

Description: This viewgraph presentation discusses project description, turbulence models, and computational engine and results for second-order closures for compressible turbulence.

Description: This viewgraph presentation gives a profile of Advanced Scientific Computing (ASC) Ltd., applications, clients and clients' needs, ASC's directions, and how the Center for Modeling of Turbulence and Transition (CMOTT) can help.

Description: This viewgraph presentation discusses the following: STAR-CD computational features; STAR-CD turbulence models; common features of industrial complex flows; industry-specific CFD development requirements; applications and experiences of industrial complex flows, including flow in rotating disc cavities, diffusion hole film cooling, internal blade cooling, and external car aerodynamics; and conclusions on turbulence modeling needs.

Description: This viewgraph presentation provides a brief review of two-equation eddy-viscosity models (TEM's) from the perspective of applied CFD. It provides objective assessment of both well-known and newer models, compares model predictions from various TEM's with experiments, identifies sources of modeling error and gives historical perspective of their effects on model performance and assessment, and recommends directions for future research on TEM's.

Description: This viewgraph presentation discusses the following: introduction to CFD Research Corporation; experiences with two-equation models - models used, numerical difficulties, validation and applications, and strengths and weaknesses; and answers to three questions posed by the workshop organizing committee - what are your customers telling you, what are you doing in-house, and how can NASA-CMOTT (Center for Modeling of Turbulence and Transition) help.

Description: This viewgraph presentation concludes that a Monte Carlo probability density function (PDF) solution successfully couples with an existing finite volume code; PDF solution method applied to turbulent reacting flows shows good agreement with data; and PDF methods must be run on parallel machines for practical use.

Description: This viewgraph presentation discusses (1) turbulence modeling: challenges in turbulence modeling, desirable attributes of turbulence models, turbulence models in FLUENT, and examples using FLUENT; and (2) combustion modeling: turbulence-chemistry interaction and FLUENT equilibrium model. As of now, three turbulence models are provided: the conventional k-epsilon model, the renormalization group model, and the Reynolds-stress model. The renormalization group k-epsilon model has broadened the range of applicability of two-equation turbulence models. The Reynolds-stress model has proved useful for strongly anisotropic flows such as those encountered in cyclones, swirlers, and combustors. Issues remain, such as near-wall closure, with all classes of models.

Description: This viewgraph presentation discusses joint velocity-scalar PDF method; turbulent combustion modeling issues for gas turbine combustors; PDF calculations for a recirculating flow; stochastic dissipation model; joint PDF calculations for swirling flows; spray calculations; reduced kinetics/manifold methods; parallel processing; and joint PDF focus areas.

Description: An algorithm has been developed for time-dependent forced convective diffusion-reaction having convection by a recirculating flow field within the drop that is hydrodynamically coupled at the interface with a convective external flow field that at infinity becomes a uniform free-streaming flow. The concentration field inside the droplet is likewise coupled with that outside by boundary conditions at the interface. A chemical reaction can take place either inside or outside the droplet, or reactions can take place in both phases. The algorithm has been implemented, and for comparison results are shown here for the case of no reaction in either phase and for the case of an external first order reaction, both for unsteady behavior. For pure interphase mass transfer, concentration isocontours, local and average Sherwood numbers, and average droplet concentrations have been obtained as a function of the physical properties and external flow field. For mass transfer enhanced by an external reaction, in addition to the above forms of results, we present the enhancement factor, with the results now also depending upon the (dimensionless) rate of reaction.

Description: The application of spectral methods, using a Chebyshev collocation scheme, to solve hydrodynamic stability problems is demonstrated on the Benard problem. Implementation of the Chebyshev collocation formulation is described. The performance of the spectral scheme is compared with that of a 2nd order finite difference scheme. An exact solution to the Marangoni-Benard problem is used to evaluate the performance of both schemes. The error of the spectral scheme is at least seven orders of magnitude smaller than finite difference error for a grid resolution of N = 15 (number of points used). The performance of the spectral formulation far exceeded the performance of the finite difference formulation for this problem. The spectral scheme required only slightly more effort to set up than the 2nd order finite difference scheme. This suggests that the spectral scheme may actually be faster to implement than higher order finite difference schemes.

Description: The effects of different turbulence boundary conditions were examined for two classical flows: a turbulent plane free shear layer and a flat plate turbulent boundary layer with zero pressure gradient. The flow solver used was DTNS, an incompressible Reynolds averaged Navier-Stokes solver with k-epsilon turbulence modeling, developed at the U.S. Navy David Taylor Research Center. Six different combinations of turbulence boundary conditions at the inflow boundary were investigated: In case 1, 'exact' k and epsilon profiles were used; in case 2, the 'exact' k profile was used, and epsilon was extrapolated upstream; in case 3, both k and epsilon were extrapolated; in case 4, the turbulence intensity (I) was 1 percent, and the turbulent viscosity (mu(sub t)) was equal to the laminar viscosity; in case 5, the 'exact' k profile was used and mu(sub t) was equal to the laminar viscosity; in case 6, the I was 1 percent, and epsilon was extrapolated. Comparisons were made with experimental data, direct numerical simulation results, or theoretical predictions as applicable. Results obtained with DTNS showed that turbulence boundary conditions can have significant impacts on the solutions, especially for the free shear layer.

Description: The exchange of models is one of the most serious problems currently encountered in the practice of spacecraft thermal analysis. Essentially, the problem originates in the diversity of computing environments that are used across different sites, and the consequent proliferation of native tool formats. Furthermore, increasing pressure to reduce the development's life cycle time has originated a growing interest in the so-called spacecraft concurrent engineering. In this context, the realization of the interdependencies between different disciplines and the proper communication between them become critical issues. The use of a neutral format represents a step forward in addressing these problems. Such a means of communication is adopted by consensus. A neutral format is not directly tied to any specific tool and it is kept under stringent change control. Currently, most of the groups promoting exchange formats are contributing with their experience to STEP, the Standard for Exchange of Product Model Data, which is being developed under the auspices of the International Standards Organization (ISO 10303). This paper presents the different efforts made in Europe to provide the spacecraft thermal analysis community with a Thermal Neutral Format (TNF) based on STEP. Following an introduction with some background information, the paper presents the characteristics of the STEP standard. Later, the first efforts to produce a STEP Spacecraft Thermal Application Protocol are described. Finally, the paper presents the currently harmonized European activities that follow up and extend earlier work on the area.

Description: Using global interpolation functions (GIF's) boundary element solutions are obtained for two-dimensional laminar flows. Two schemes are proposed for handling the convective terms. The first treats convection as a forcing function, and converts the flow equations to pseudo-Poisson equations. In the second scheme, some convective effect is incorporated into the fundamental solution used in constructing the pertinent integral equations. The lid-driven cavity flow is selected as the benchmark problem.

Description: A two dimensional finite volume method is used to predict the film coefficients in the transitional flow region (laminar or turbulent) for the radiator panel tubes. The code used to perform this analysis is CAST (Computer Aided Simulation of Turbulent Flows). The information gathered from this code is then used to augment a Sinda85 model that predicts overall performance of the radiator. A final comparison is drawn between the results generated with a Sinda85 model using the Sinda85 provided transition region heat transfer correlations and the Sinda85 model using the CAST generated data.

Description: This paper presents a probabilistic one-dimensional finite element model for heat transfer processes in porous heat exchangers. The Galerkin approach is used to develop the finite element matrices. Some of the submatrices are asymmetric due to the presence of the flow term. The Neumann expansion is used to write the temperature distribution as a series of random variables, and the expectation operator is applied to obtain the mean and deviation statistics. To demonstrate the feasibility of the formulation, a one-dimensional model of heat transfer phenomenon in superfluid flow through a porous media is considered. Results of this formulation agree well with the Monte-Carlo simulations and the analytical solutions. Although the numerical experiments are confined to parametric random variables, a formulation is presented to account for the random spatial variations.

Description: This paper presents a numerical scheme, based on the finite element method, to solve strongly coupled fluid flow and heat transfer problems. The surface radiation effect for gray, diffuse and isothermal surfaces is considered. A procedure for obtaining the view factors between the radiating surfaces is discussed. The overall solution strategy is verified by comparing the available results with those obtained using this approach. An analysis of a thermosyphon is undertaken and the effect of considering the surface radiation is clearly explained.

Description: An algorithm has been developed for the forced convective diffusion-reaction problem for convection inside and outside a droplet by a recirculating flow field hydrodynamically coupled at the droplet interface with an external flow field that at infinity becomes a uniform streaming flow. The concentration field inside the droplet is likewise coupled with that outside by boundary conditions at the interface. A chemical reaction can take place either inside or outside the droplet or reactions can take place in both phases. The algorithm has been implemented and results are shown here for the case of no reaction and for the case of an external first order reaction, both for unsteady behavior. For pure interphase mass transfer, concentration isocontours, local and average Sherwood numbers, and average droplet concentrations have been obtained as a function of the physical properties and external flow field. For mass transfer enhanced by an external reaction, in addition to the above forms of results, we present the enhancement factor, with the results now also depending upon the (dimensionless) rate of reaction.

Description: The critical state of vortex cores downstream of vortex breakdown has been studied. Base vortical flows were computed using the Reynolds-averaged, axisymmetric Navier-Stokes equations. Standard K - epsilon, RNG and second-order Reynolds stress models were employed. Results indicate that the return to supercriticality is highly dependent on the turbulence model. The K - epsilon model predicted a rapid return of the vortex to supercritical conditions, the location of which showed little sensitivity to changes in the swirl ratio. The Reynolds stress model predicted that the vortex remains subcritical to the end of the domain for each of the swirl ratios employed, and provided results in qualitative agreement with experimental work. The RNG model produced intermediate results, with a downstream movement in the critical location with increasing swirl. Calculations for which area reductions were introduced at the exit in a subcritical flow were also performed using the Reynolds stress model. The structure of the resulting recirculation zone was altered significantly. However, when area reductions were employed within supercritical flows as predicted using the two-equation models, no significant influence on the recirculation zone was noted.

Description: The view factors which are used in diffuse-gray radiation enclosure calculations are often computed by approximate numerical integrations. These approximately calculated view factors will usually not satisfy the important physical constraints of reciprocity and closure. In this paper several view-factor rectification algorithms are reviewed and a rectification algorithm based on a least-squares numerical filtering scheme is proposed with both weighted and unweighted classes. A Monte-Carlo investigation is undertaken to study the propagation of view-factor and surface-area uncertainties into the heat transfer results of the diffuse-gray enclosure calculations. It is found that the weighted least-squares algorithm is vastly superior to the other rectification schemes for the reduction of the heat-flux sensitivities to view-factor uncertainties. In a sample problem, which has proven to be very sensitive to uncertainties in view factor, the heat transfer calculations with weighted least-squares rectified view factors are very good with an original view-factor matrix computed to only one-digit accuracy. All of the algorithms had roughly equivalent effects on the reduction in sensitivity to area uncertainty in this case study.

Description: A transient, one-dimensional numerical code is developed to model the liquid motion in an axial groove with square cross section. Axial variation in liquid level, shear stress and heat transfer between the groove wall and the liquid, evaporation and transient body forces are accounted for in the model. Dryout and rewet of the groove are allowed; the front location is determined numerically using conservation of mass and linear extrapolation. Several numerical test results are presented and discussed.

Description: A higher-order finite-difference technique is developed to calculate the developing-flow field of steady incompressible laminar flows in the entrance regions of circular pipes. Navier-Stokes equations governing the motion of such a flow field are solved by using this new finite-difference scheme. This new technique can increase the accuracy of the finite-difference approximation, while also providing the option of using unevenly spaced clustered nodes for computation such that relatively fine grids can be adopted for regions with large velocity gradients. The velocity profile at the entrance of the pipe is assumed to be uniform for the computation. The velocity distribution and the surface pressure drop of the developing flow then are calculated and compared to existing experimental measurements reported in the literature. Computational results obtained are found to be in good agreement with existing experimental correlations and therefore, the reliability of the new technique has been successfully tested.

Description: The Systems Improved Numerical Differencing Analyzer and Fluid Integrator (SINDA/FLUINT) code has often been used to determine the transient and steady-state response of various thermal and fluid flow networks. While this code is an often used design and analysis tool, the validation of this program has been limited to a few simple studies. For the current study, the SINDA/FLUINT code was compared to four different analytical solutions. The thermal analyzer portion of the code (conduction and radiative heat transfer, SINDA portion) was first compared to two separate solutions. The first comparison examined a semi-infinite slab with a periodic surface temperature boundary condition. Next, a small, uniform temperature object (lumped capacitance) was allowed to radiate to a fixed temperature sink. The fluid portion of the code (FLUINT) was also compared to two different analytical solutions. The first study examined a tank filling process by an ideal gas in which there is both control volume work and heat transfer. The final comparison considered the flow in a pipe joining two infinite reservoirs of pressure. The results of all these studies showed that for the situations examined here, the SINDA/FLUINT code was able to match the results of the analytical solutions.

Description: An IBM Personal Computer (PC) version of the Groove Analysis program (GAP) was developed to predict the steady state heat transport capability of an axially grooved heat pipe for a specified groove geometry and working fluid. In the model, the capillary limit is determined by the numerical solution of the differential equation for momentum conservation with the appropriate boundary conditions. This governing equation accounts for the hydrodynamic losses due to friction in liquid and vapor flows and due to liquid/vapor shear interaction. Back-pumping in both 0-g and 1-g is accounted for in the boundary condition at the condenser end. Slug formation in 0-g and puddle flow in 1-g are also considered in the model. At the user's discretion, the code will perform the analysis for various fluid inventories (undercharge, nominal charge, overcharge, or a fixed fluid charge) and heat pipe elevations. GAP will also calculate the minimum required heat pipe wall thickness for pressure containment at design temperatures that are greater than or lower than the critical temperature of the working fluid. This paper discusses the theory behind the development of the GAP model. It also presents the many useful and powerful capabilities of the model. Furthermore, a correlation of flight test performance data and the predictions using GAP are presented and discussed.

Description: The paper presents the numerical solution of heat and mass transfer during cross-flow (orthogonal) mixed convection. In this class of flow, a buoyancy-driven transport in the vertical direction and a forced convective flow in the horizontal direction results in a three-dimensional boundary layer structure adjacent to the plate. The rates of heat and mass transfer are determined by a combined influence of the two transport processes. The equations for the conservation of mass, momentum, energy, and species concentration were solved along with appropriate boundary conditions to determine the distributions of velocity components, temperature, and concentration across the thickness of the boundary layer at different locations on the plate. Results were expressed in dimensionless form using Reynolds number, Richardson number for heat transfer, Richardson number for mass transfer, Prandtl number, and Schmidt number as parameters. It was found that the transport is dominated by buoyancy at smaller vertical locations and at larger distances away from the forced convection leading edge. Effects of forced convection appeared to be very strong at smaller horizontal distances from the leading edge. The cross stream forced convection enhanced the rate of heat and mass transfer by a very significant amount.

Description: The paper presents a detailed theoretical analysis of the process of gas absorption to a thin liquid film adjacent to a horizontal rotating disk. The film is formed by the impingement of a controlled liquid jet at the center of the disk and subsequent radial spreading of liquid along the disk. The chemical reaction between the gas and the liquid film can be expressed as a zero-order homogeneous reaction. The process was modeled by establishing equations for the conservation of mass, momentum, and species concentration and solving them analytically. A scaling analysis was used to determine dominant transport processes. Appropriate boundary conditions were used to solve these equations to develop expressions for the local concentration of gas across the thickness of the film and distributions of film height, bulk concentration, and Sherwood number along the radius of the disk. The partial differential equation for species concentration was solved using the separation of variables technique along with the Duhamel's theorem and the final analytical solution was expressed using confluent hypergeometric functions. Tables for eigenvalues and eigenfunctions are presented for a number of reaction rate constants. A parametric study was performed using Reynolds number, Ekman number, and dimensionless reaction rate as parameters. At all radial locations, Sherwood number increased with Reynolds number (flow rate) as well as Ekman number (rate of rotation). The enhancement of mass transfer due to chemical reaction was found to be small when compared to the case of no reaction (pure absorption), but the enhancement factor was very significant when compared to pure absorption in a stagnant liquid film. The zero-order reaction processes considered in the present investigation included the absorption of oxygen in aqueous alkaline solutions of sodiumdithionite and rhodium complex catalyzed carbonylation of methanol. Present analytical results were compared to previous theoretical results for limiting conditions, and were found to have very good agreement.

Description: Forced convective diffusion-reaction is considered for viscous axisymmetric extensional convecting velocity in the neighborhood of a sphere. For Peclet numbers in the range 0.1 less than or equal to Pe less than or equal to 500 and for Damkohler numbers increasing with increasing Pe but in the overall range 0.02 less than or equal to Da less than or equal to 10, average and local Sherwood numbers have been computed. By introducing the eigenfunction expansion c(r, Theta) = Sum of c(n)(r)P(n)(cos Theta) into the forced convective diffusion equation for the concentration of a chemical species undergoing a first order homogeneous reaction and by using properties of the Legendre functions Pn(cos Theta), the variable coefficient PDE can be reduced to a system of N + 1 second order ODEs for the radial functions c(sub n)(r), n = 0, 1, 2,..., N. The adaptive grid algorithm of Pereyra and Lentini can be used to solve the corresponding 2(N + 1) first order differential equations as a two-point boundary value problem on 1 less than or equal to r less than or equal to r(sub infinity). Convergence of the expansion for a specific value of N can thus be established and provides 'spectral' behavior as well as the full concentration field c(r, Theta).

Description: On July 4,1997, the Mars Pathfinder spacecraft lands on Mars and starts conducting technological and scientific experiments. One experiment, the Alpha-Proton-X-ray Spectrometer, uses a sensor head placed against rocks and soil to determine their composition. To guarantee proper placement, a deployment mechanism mounted on the Mars Rover aligns the sensor head to within 20 deg of the rock and soil surfaces. In carrying out its task, the mechanism mimics the action of a human hand and arm. Consisting of a flexible wrist, a parallel link arm, a brush dc motor actuator, and a revolutionary non-pyrotechnic fail-safe release device, the mechanism correctly positions the sensor head on rocks as high as 0.29 m and on targets whose surfaces are tilted as much as 45 deg from the nominal orientation of the sensor head face. The mechanism weighs less than 0.5 kg, can withstand 100 g's, and requires less than 2.8 N x m of actuation torque. The fail-safe coupler utilizes Cerrobend, a metal alloy that melts at 60 C, to fuse the actuator and the rest of the mechanism together. A film heater wrapped around the coupler melts the metal, and Negator springs drive the mechanism into its stowed position. The fail-safe actuates using 6.75 Watts for 5 minutes in the event of an actuator failure.

Description: Annotation is a key activity of data analysis. However, current systems for data analysis focus almost exclusively on visualization. We propose a system which integrates annotations into a visualization system. Annotations are embedded in 3D data space, using the Post-it metaphor. This embedding allows contextual-based information storage and retrieval, and facilitates information sharing in collaborative environments. We provide a traditional database filter and a Magic Lens filter to create specialized views of the data. The system has been customized for fluid flow applications, with features which allow users to store parameters of visualization tools and sketch 3D volumes.

Description: Fluctuating surface pressure measurements have been made to investigate the effectiveness of boundary layer separators (BLS's) in reducing the fluctuating pressure loads produced by separated shock wave turbulent boundary layer interactions. Measurements have been made under unswept and swept compression corner interactions in a Mach 5 flow. BLS's fix the separation location and eliminate the large-amplitude, low-frequency fluctuating pressure loads upstream of the compression corners. The loads on the unswept compression corner face are reduced by as much as 59%. The BLS's also shift the mean pressure distribution on the unswept corner face in the streamwise direction. Results show that the loads on the corner face vary with the BLS height and the distance between the BLS and the compression corner. Suggestions for the optimum placement and the use of the BLS's are also made.

Description: Fluctuating wall pressure measurements have been made in a separated shock wave/turbulent boundary layer interaction produced by an unswept compression corner in a Mach 5 flow. Wheeler doublet vortex generators were placed 15.8 boundary layer thicknesses upstream of the corner to study their effect on the fluctuating pressure loads produced by the translating separation shock. The vortex generators produced significant three-dimensionality in an otherwise two-dimensional interaction. They reduced the upstream influence and the length of the region of shock motion by 60% and 64%, respectively, decreased the maximum wall pressure rms by 23%, and shifted the fluctuations to a higher frequency band. The maximum fraction of energy in the 100-500 Hz frequency band is decreased by 11%. These changes are due to a fuller boundary layer profile, a weaker separation shock, and increased boundary layer turbulence causing increased separation shock jitter.

Description: A semianalytical solution is presented for the unsteady pressure field of a vortical gust interacting with a flat-plate airfoil in subsonic flow. The solution will serve as a benchmark for evaluating the accuracy and efficiency of time dependent numerical schemes. The specific case considered corresponds to the ICASE benchmark problem number 6. The results are compared with those of asymptotic theories for high frequency and show excellent agreement.

Description: A numerical method has been developed in order to address aeroacoustic problems modeled by the linearized Euler equations. A weak formulation of the equations leads to a time-dependent equation for the test functions. The basic solver being one dimensional, two dimensional problems are handled by directional splitting. This method shows low dissipation and dispersion errors.

Description: Accurate capturing of discontinuities within compressible flow computations is achieved by coupling a suitable solver with an automatic adaptive mesh algorithm for unstructured triangular meshes. The mesh adaptation procedures developed rely on non-hierarchical dynamical local refinement/derefinement techniques, which hence enable structural optimization as well as geometrical optimization. The methods described are applied for a number of the ICASE test cases are particularly interesting for unsteady flow simulations.

Description: This paper discusses the issues which arise when combining multigrid strategies with adaptive meshing techniques for solving steady-state problems on unstructured meshes. A basic strategy is described, and demonstrated by solving several inviscid and viscous flow cases. Potential inefficiencies in this basic strategy are exposed, and various alternate approaches are discussed, some of which are demonstrated with an example. Although each particular approach exhibits certain advantages, all methods have particular drawbacks, and the formulation of a completely optimal strategy is considered to be an open problem.

Description: The dynamic solution adaptive grid algorithm, DSAGA3D, is extended to automatically adapt 2-D structured multi-block grids, including adaption of the block boundaries. The extension is general, requiring only input data concerning block structure, connectivity, and boundary conditions. Imbedded grid singular points are permitted, but must be prevented from moving in space. Solutions for workshop cases 1 and 2 are obtained on multi-block grids and illustrate both increased resolution of and alignment with the solution. A mesh quality assessment criteria is proposed to determine how well a given mesh resolves and aligns with the solution obtained upon it. The criteria is used to evaluate the grid quality for solutions of workshop case 6 obtained on both static and dynamically adapted grids. The results indicate that this criteria shows promise as a means of evaluating resolution.

Description: This paper describes adaptive grid methods developed specifically for compressible flow computations. The basic flow solver is a finite-volume implementation of Roe's flux difference splitting scheme or arbitrarily moving unstructured triangular meshes. The grid adaptation is performed according to geometric and flow requirements. Some results are included to illustrate the potential of the methodology.

Description: Solution-adaptive methods based on cutting bodies out of Cartesian grids are gaining popularity now that the ways of circumventing the accuracy problems associated with small cut cells have been developed. Researchers are applying Cartesian-based schemes to a broad class of problems now, and, although there is still development work to be done, it is becoming clearer which problems are best suited to the approach (and which are not). The purpose of this paper is to give a candid assessment, based on applying Cartesian schemes to a variety of problems, of the strengths and weaknesses of the approach as it is currently implemented.

Description: A three dimensional finite volume scheme based on hybrid grids containing both tetrahedral and hexahedral cells is presented. The application to hybrid grids offers the possibility to combine the flexibility of tetrahedral meshes with the accuracy of hexahedral grids. An algorithm to compute a dual mesh for the entire computational domain was developed. The dual mesh technique guarantees conservation in the whole flow field even at interfaces between hexahedral and tetrahedral domains and enables the employment of an accurate upwind flow solver. The hybrid mesh can be adapted to the solution by dividing cells in areas of insufficient resolution. The method is tested on different viscous and inviscid cases for hypersonic, transonic and subsonic flows.

Description: A method for generating high quality unstructured triangular grids for high Reynolds number Navier-Stokes calculations about complex geometries is described. Careful attention is paid in the mesh generation process to resolving efficiently the disparate length scales which arise in these flows. First the surface mesh is constructed in a way which ensures that the geometry is faithfully represented. The volume mesh generation then proceeds in two phases thus allowing the viscous and inviscid regions of the flow to be meshed optimally. A solution-adaptive remeshing procedure which allows the mesh to adapt itself to flow features is also described. The procedure for tracking wakes and refinement criteria appropriate for shock detection are described. Although at present it has only been implemented in two dimensions, the grid generation process has been designed with the extension to three dimensions in mind. An implicit, higher-order, upwind method is also presented for computing compressible turbulent flows on these meshes. Two recently developed one-equation turbulence models have been implemented to simulate the effects of the fluid turbulence. Results for flow about a RAE 2822 airfoil and a Douglas three-element airfoil are presented which clearly show the improved resolution obtainable.

Description: This paper presents an hp-adaptive discontinuous Galerkin method for linear hyperbolic conservation laws. A priori and a posteriori error estimates are derived in mesh-dependent norms which reflect the dependence of the approximate solution on the element size (h) and the degree (p) of the local polynomial approximation. The a posteriori error estimate, based on the element residual method, provides bounds on the actual global error in the approximate solution. The adaptive strategy is designed to deliver an approximate solution with the specified level of error in three steps. The a posteriori estimate is used to assess the accuracy of a given approximate solution and the a priori estimate is used to predict the mesh refinements and polynomial enrichment needed to deliver the desired solution. Numerical examples demonstrate the reliability of the a posteriori error estimates and the effectiveness of the hp-adaptive strategy.

Description: Two symposia related to the exploration of the Venus and Mars atmospheres produced papers which described the results of spacecraft-based and ground-based observations of these two planets. The topics addressed included solar wind interactions, molecular ions, planetary magnetospheres, planetary ionospheres, planetary thermospheres, gravity waves, the homopause, radiative transfer, radiative balance, atmospheric water, planetary lightning, very low frequency bursts, cosmic-ray ionization, planetary loss, and future missions to study Mars aeronomy, energetic particles near Mars, radar exploration of Mars, and the Mars Environmental Survey (MESUR) mission.

Description: Navier-Stokes equation as discretized by new flux conserving method proposed by Chang and Scott results in the system: vector F(vector x) = 0, where F is a vector valued function. The Optimization method we use is based on Quasi-Newton methods: given a nonlinear function vector F(vector x) = 0, we solve, Delta(vector x) = -BF(vector x), where Delta(vector x) is the correction term and B is the inverse Jacobian of F(x). Then, iteratively, vector(x(sub (i+1))) = vector(x (sub i)) + alpha.Delta(vector x(sub i)), where alpha is a line search correction term determined by a line search routine. We use the BFCG's update the Jacobian matrix B(sub k) at each iteration. It is well known that B(sub k) approaches B(*) at the solution X(*). This algorithm has several advantages over the Newton-Raphson method. For example, we do not need to calculate the Jacobian matrix at each iteration which is computationally very expensive.

Description: A newly discovered source of energetic particle streams directed at Venus' nightside may cause several exosphere phenomena. The streams evolve from dispersed cosmoids, a surreptitious population of meteoroids in nearly hyperbolic orbits, measured with three dust experiments on Pioneer 10/11. Loose fragile comet-like ensembles of frozen volatile material, they are orders of magnitude more populous than short-period meteoroids. Dispersion is forced near a planet masking the meteor signature at Earth, but recently terrestrial exosphere interactions have been detected principally from VLF radar returns and as suddenly formed layers of neutral sodium and iron probed by lidar. Venus influx is estimated greater than 10(exp -14)(g)/((sq cm)(s)) with nightside directed kinetic power of greater than 0.3 (erg)/((sq cm)s)). Compared to Earth, Venus, lacking a dipole field with almost no rotation, has a readily recognized set of nightside interaction signatures of these downward energetic particle streams: i.e., density depressions of the neutral thermosphere, hydrogen and helium bulges, a non-disappearing ionosphere, electron holes, nightglows, VLF bursty signals, and local magnetic fields. From Pioneer Venus Orbiter (PVO) measurements, the cosmoid population, its temporal and solar azimuthal variation may be determined.

Description: High-resolution, near-infrared (1.09 to 2.5 micrometers) spectra of the night side of Venus have been obtained in 1990 and 1991 using the Fourier Transform Spectrometer at the 3.6-m Canada-France-Hawaii telescope. Absorptions due to H2O were detected in spectral windows near 2.3, 1.74, and 1.18 micrometers. Our analysis of these absorptions constrains the abundance of water vapor in three different altitude ranges located between the clouds and the surface: 30-40 km, 15-25 km and 0-15 km. A constant water vapor mixing ratio of 30 +/- 15 ppm below the clouds can fit the observations. These values are consistent with recent near-infrared studies of the night side of Venus at lower spectral resolution. The atmosphere of Venus appears to be dryer than originally suggested by the in-situ measurements made by the Pioneer Venus and Venera mass-spectrometers and gas-chromatographs.

Description: Using Pioneer Venus data, hydrogen and deuterium ions are shown to escape from the hydrogen bulge region in the nightside ionosphere. The polarization electric field propels these light ions upward through the ionosphere and into the ion-exosphere, where H(+) and D(+) continue to be accelerated away from Venus and move into the ionotail and beyond. The vertical flow speeds of H(+) and D(+) are found to be about the same; therefore, selective escape between H(+) and D(+) is negligible for this mechanism. Present day planetary loss rates of about 8.6 x 10(exp 25)/s and 3.2 X 10(exp 23)/s were obtained for H(+) and D(+), respectively. Such rates, persisting over a few billion years, should have significantly affected the planetary water budget.

Description: Advanced airbreathing propulsion systems used in Mach 4-6 mission scenarios, usually consist of a single integrated turboramjet or as in this study, a turbojet housed in an upper bay with a separate ramjet housed in a lower bay. As the engines transition from turbojet to ramjet, there is an operational envelope where both engines operate simultaneously. One nozzle concept under consideration has a common nozzle, where the plumes from the turbojet and ramjet interact with one another as they expand to ambient conditions. In this paper, the two plumes interact at the end of a common 2-D cowl, when they both reach an approximate Mach 3.0 condition and then jointly expand to Mach 3.6 at the common nozzle exit plane. At this condition, the turbojet engine operated at a higher NPR than the ramjet, where the turbojet overpowers the ramjet plume, deflecting it approximately 12 degrees downward and in turn the turbojet plume is deflected 6 degrees upward. In the process, shocks were formed at the deflections and a shear layer formed at the confluence of the two jets. This particular case was experimentally tested and the data used to compare with the PARC3D code with k-kl two equation turbulence model. The 2-D and 3-D centerline CFD solutions are in good agreement, but as the CFD solutions approach the outer sidewall, a slight variance occurs. The outer wall boundary layers are thin and do not present much of an interaction, however, where the confluence interaction shocks interact with the thin boundary layer on the outer wall, strong vortices run down each shock causing substantial disturbances in the boundary layer. These disturbances amplify somewhat as they propagate downstream axially from the confluence point. The nozzle coefficient (CFG) is reduced 1/2 percent as a result of this sidewall interaction, from 0.9850 to 0.9807. This three-dimensional reduction is in better agreement with the experimental value of 0.9790.

Description: Wavenumber transition and hysteresis in a highly unstable baroclinic flow are investigated using a high-resolution spectral numerical model. As the flow becomes more supercritical, the dominant wave gradually shifts from the most unstable wave predicted by the linear theory to a longer wave with a larger time-averaged amplitude, while the rectified mean flow attains a stronger shear at the center of the channel. The numerical results display a complex hysteresis behavior, which occurs not only between the states of different dominant wavenumbers, but also between the states of identical dominant wavenumber but of different dynamic characteristics. In a certain parameter range three stable states, each with different dominant wavenumber, are possible, and in another parameter range four stable states are possible, among them three stable states with an identical dominant wave. The numerical results suggest that a multiple weather regime exists even without external forcing in which the flow aperiodically varies between two distinct behaviors. The effects of stable higher harmonics are assessed and it is found that their presence contributes not only to the better approximation of the model solutions but also to the selection of the final equilibrium state, due to the chaotic nature of the initial transient period.

Description: We demonstrate the use of Filtererd Rayleigh Scattering and a 3D reconstruction technique to interrogate the highly three dimensional flow field inside of a supersonic inlet model. A 3 inch by 3 inch by 2.5 inch volume is reconstructed yielding 3D visualizations of the crossing shock waves and of the boundary layer. In this paper we discuss the details of the techniques used, and present the reconstructured 3D images.

Description: Substantial evidence suggests that a UV spectrally Absorbing Material (UV-SAM) exists on Triton's surface. This evidence is found in the positive slope in Triton's spectrum from the UV to the near-IR, and the increasing contrast in Triton's light curve in the blue and UV. Although it is now widely-thought that UV-SAMs exist on Triton, little is known about their distribution and spectral properties. The goal of this NDAP Project is to determine the spatial distribution and geological context of the UV-SaM material. We hope to determine if UV-SAMs on Triton are correlated with geologic wind streaks, craters, calderas, geomorphic/topographic units, regions containing (or lacking) volatile frosts, or some other process (e.g., magnetospheric interactions). Once the location and distribution of UV-SAMs has been determined, further constraints on their composition cable made by analyzing the spectrographic data set. To accomplish these goals, various data sets will be used, including Voyager 2 UV and visible images of Triton's surface, IUE and HST spectra of Triton, and a geologic map of the surface based on voyager 2 and spectrophotometric data. The results of this research will be published in the planetary science literature.

Description: The consistency of second-order closure models with results from hydrodynamic stability theory is analyzed for the simplified case of homogeneous turbulence. In a recent study, Speziale, Gatski, and MacGiolla Mhuiris showed that second-order closures are capable of yielding results that are consistent with hydrodynamic stability theory for the case of homogeneous shear flow in a rotating frame. It is demonstrated in this paper that this success is due to the fact that the stability boundaries for rotating homogeneous shear flow are not dependent on the details of the spatial structure of the disturbances. For those instances where they are -- such as in the case of elliptical flows where the instability mechanism is more subtle -- the results are not so favorable. The origins and extent of this modeling problem are examined in detail along with a possible resolution based on rapid distortion theory (RDT) and its implications for turbulence modeling.

Description: Two linearized solvers (time and frequency domain) based on a high resolution numerical scheme are presented. The basic approach is to linearize the flux vector by expressing it as a sum of a mean and a perturbation. This allows the governing equations to be maintained in conservation law form. A key difference between the time and frequency domain computations is that the frequency domain computations require only one grid block irrespective of the interblade phase angle for which the flow is being computed. As a result of this and due to the fact that the governing equations for this case are steady, frequency domain computations are substantially faster than the corresponding time domain computations. The linearized equations are used to compute flows in turbomachinery blade rows (cascades) arising due to blade vibrations. Numerical solutions are compared to linear theory (where available) and to numerical solutions of the nonlinear Euler equations.