ALBERT

Description: The low-g fluids management group with the Center for Space Construction is engaged in active research on the following topics: gauging; venting; controlling contamination; sloshing; transfer; acquisition; and two-phase flow. Our basic understanding of each of these topics at present is inadequate to design space structures optimally. A brief report is presented on each topic showing the present status, recent accomplishings by our group and our plans for future research. Reports are presented in graphic and outline form.

Description: Through the Surveyor 3 and 7, and Apollo 11-17 missions a knowledge of the mechanical properties of Lunar regolith were gained. These properties, including material cohesion, friction, in-situ density, grain-size distribution and shape, and porosity, were determined by indirect means of trenching, penetration, and vane shear testing. Several of these properties were shown to be significantly different from those of terrestrial soils, such as an interlocking cohesion and tensile strength formed in the absence of moisture and particle cementation. To characterize the strength and deformation properties of Lunar regolith experiments have been conducted on a lunar soil simulant at various initial densities, fabric arrangements, and composition. These experiments included conventional triaxial compression and extension, direct tension, and combined tension-shear. Experiments have been conducted at low levels of effective confining stress. External conditions such as membrane induced confining stresses, end platten friction and material self weight have been shown to have a dramatic effect on the strength properties at low levels of confining stress. The solution has been to treat these external conditions and the specimen as a full-fledged boundary value problem rather than the idealized elemental cube of mechanics. Centrifuge modeling allows for the study of Lunar soil-structure interaction problems. In recent years centrifuge modeling has become an important tool for modeling processes that are dominated by gravity and for verifying analysis procedures and studying deformation and failure modes. Centrifuge modeling is well established for terrestrial enginering and applies equally as well to Lunar engineering. A brief review of the experiments is presented in graphic and outline form.

Description: Experimental and theoretical studies have been conducted to determine critical parameters at the onset of nonlinear counterflow in He II below the lambda point of He-4. Critical temperature differences have been measured in porous media for zero net mass flow and for Darcy permeabilities in the order of magnitude range from 10 to the -10th to 10 to the -8th sq cm. The normalized critical temperature gradients, which covered the liquid temperature range of 1.5 K to the lambda temperature, are found to vary with T proportional to the ratio of the superfluid density to the normal fluid density. This liquid temperature dependence appears to be consistent with duct data which are limited at low temperature by a Reynolds number criterion.

Description: Nonthermal radio emissions from earth, Jupiter, Saturn, and Uranus are reviewed. The dominant source of emission at each planet appears to be AKR-like auroral emission in the X-mode. O-mode emissions are substantially responsible. There is a remarkably constant scaling factor relating the total solar wind input power into each planetary system and the AKR-like auroral emissions.

Description: This article presents the results of a fit of a model of the Martian satellite orbits to earthbased and spacecraft-based observations. An assessment of the orbit accuracies is given and the orbits are compared with those obtained by previous investigators.

Description: The petrology and chronology of early lunar crust is examined using the least equivocal of the available petrographic and age data on lunar rock samples, and the possible processes which produced the lunar crust are discussed. The results suggest that the lunar anorthositic crust was formed by about 120 Ma after the primary accretion of the moon at 4.56 Ga. At least some members of the diverse Mg-suites of rocks, such as norites, troctolites, and dunites, crystallized within a very few 100s of Ma after 4.56 Ga. A trace-element-rich material (KREEP) was formed by about 4.3 Ga ago, and this residue was subsequently reworked in melting and impact processes such that most samples which contain it have ages around 3.9-4.0 Ga. The findings also suggest that the onset of ferrous mare basalt volcanism began about 4.33 Ga, much earlier than was once assumed, and was still in process before the end of the most intense period of bombardment (3.9-4.0 Ga ago).

Description: An analysis is made based upon the concept that the velocity fluctuations, and therefore, the Reynolds stresses, driven by the instability of the original flow grow until a new stable state is approached. The Reynolds stresses incorporated into the Orr-Sommerfeld equation are coupled with the main flow such that all the imaginary parts of the complex eigenvalues vanish, i.e., the original instability is eliminated. Using this stabilization principle, it is possible to find the Reynolds stresses as well as the mean velocity for plane Poiseuille flow with the Reynolds number slightly higher than the critical.

Description: The behavior of the reverse flow ceiling jet against the ventilation flow from 0.58 to 0.87 m/s was investigated in a 1/3 scale model of a wide body aircraft interior. For all tests, strong reverse-flow ceiling jets of hot gases were detected well upstream of the fire. Both thicknesses of the reverse-flow ceiling jet and the smoke layer increased with the fire-crossflow parameter. The thickness of the smoke layer where the smoke flows along the main flow below the reverse-flow ceiling jet was almost twice that of the reverse-flow ceiling jet. Detailed spatial and time-varying temperatures of the gas in the test section were measured, and velocity profiles were also measured using a temperature compensated hot film.

Description: Certain theoretical studies of boundary-layer transition are described, based on high Reynolds numbers and with attention drawn to the various nonlinear interactions and scales present. The article concentrates in particular on theories for which the mean-flow profile is completely altered from its original state. Two- and three-dimensional flow theory and conjectures on turbulent-boundary-layer structures are included. Specific recent findings noted, and in qualitative agreement with experiments, are: nonlinear finite-time break-ups in unsteady interactive boundary layers; strong vortex/wave interactions; and prediction of turbulent boundary-layer displacement- and stress sublayer-thicknesses.

Description: The findings made by the Voyager 2 Neptune encounter are reviewed. Data on the bowshock, magnetic field, magnetosphere, rings, plasma sheet, aurora, moons, and dust of Neptune are discussed. Findings made concerning Triton are summarized.

Description: The composition of lunar regolith and its attendant properties are discussed. Tables are provided listing lunar minerals, the abundance of plagioclase feldspar, pyroxene, olivine, and ilmenite in lunar materials, typical compositions of common lunar minerals, and cumulative grain-size distribution for a large number of lunar soils. Also provided are charts on the chemistry of breccias, the chemistry of lunar glass, and the comparative chemistry of surface soils for the Apollo sites. Lunar agglutinates, constructional particles made of lithic, mineral, and glass fragments welded together by a glassy matrix containing extremely fine-grained metallic iron and formed by micrometeoric impacts at the lunar surface, are discussed. Crystalline, igneous rock fragments, breccias, and lunar glass are examined. Volatiles implanted in lunar materials and regolith maturity are also addressed.

Description: Theoretical models of planetary-ring dynamics are discussed in a detailed analytical review and illustrated with graphs and diagrams. The streamline concept is introduced, and the phenomena associated with the transport of angular momentum are described. Particular attention is then given to (1) broad rings like those of Saturn (shepherding, density-wave excitation, gaps, bending-wave excitation, multiringlet structures, inner-edge shepherding, and the possibility of polar rings around Neptune), (2) narrow rings like those of Uranus (shepherding, ring shapes, and a self-gravity model of rigid precession), and (3) ring arcs like those seen in stellar-occultation observations of Neptune.

Description: Images of Neptune obtained by the narrow-angle camera of the Voyager 2 spacecraft reveal large-scale cloud features that persist for several months or longer. The features' periods of rotation about the planetary axis range from 15.8 to 18.4 hours. The atmosphere equatorward of -53 deg rotates with periods longer than the 16.05-hour period deduced from Voyager's planetary radio astronomy experiment (presumably the planet's internal rotation period). The wind speeds computed with respect to this radio period range from 20 meters per second eastward to 325 meters per second westward. Thus, the cloud-top wind speeds are roughly the same for all the planets ranging from Venus to Neptune, even though the solar energy inputs to the atmospheres vary by a factor of 1000.

Description: The theory of turbulent plumes maintained above steady line sources of buoyancy is worked out in detail within the limitations of Taylor's entrainment assumption. It is applied to the structure of a pure plume injected into a stably stratified atmosphere. Volcanic basalt eruptions that develop from long, narrow vents create line source plumes, which rise well above the magmatic fire fountains playing near the ground level. The eruption of Laki in 1783 may provide an example of this style of eruption. Flood basalts are more ancient examples. Evidence of enormous fissure eruptions that occurred in the past on Mars and Venus also exists. Owing to the different properties of the atmospheres on these two planets from those on the earth, heights of line source plumes are expected to vary in the ratios 1:6:0.6 (earth:Mars:Venus). It is very unlikely that the observed increase of sulfur dioxide above the Venusian cloud deck in 1978 could have been due to a line source volcanic eruption, even if it had been a flood basalt eruption.

Description: An analytical and numerical study of the stability of polar rings around Neptune is presented. The stability proofs are based on various methods used to study gas disks in galaxies. It is shown that stable polar rings can exist despite energy dissipation by collisions between particles. Also, four equilibrium orientations which pass nearly over the pole of Neptune are found, two of which are stable in the presence of dissipation.

Description: Issues related to the establishment of lunar observatories are briefly addressed. The advantages of placing an observatory in a crater at one of the poles, where permanent darkness exists, are pointed out, and the methods required to emplace and operate such an observatory are considered. Planning for the installation of the first set of observatory instruments is discussed.

Description: This paper presents the application of a class of multi-grid methods to the solution of the Navier-Stokes equations for two-dimensional laminar flow problems. The methods consists of combining the full approximation scheme-full multi-grid technique (FAS-FMG) with point-, line- or plane-relaxation routines for solving the Navier-Stokes equations in primitive variables. The performance of the multi-grid methods is compared to those of several single-grid methods. The results show that much faster convergence can be procured through the use of the multi-grid approach than through the various suggestions for improving single-grid methods. The importance of the choice of relaxation scheme for the multi-grid method is illustrated.

Description: Rapid distortion theory is applied to study distortion of homogeneous turbulence subject to two different axisymmetric strain modes: the axisymmetric contraction (AC, nozzle-type flow), and the axisymmetric expansion (AE, diffuser-type flow). The paper explores the differences in effects of the two axisymmetric strain modes on the anisotropy of correlations and structures of turbulence; examines the effect of dilatation on the distortion of turbulence; and provides a theoretical background for turbulence model development. It is found that velocity and vorticity fluctuations are enhanced more efficiently by contraction than by expansion; contraction produces much higher anisotropy in velocity and vorticity than expansion; root-mean-square pressure is slightly reduced during contraction, whereas it increases rapidly during expansion; and vortical structures of rodlike shape develop in a contraction flow, while disklike structures develop in an expansion flow. A simple model that reflects the dependence of turbulence evolution on structural parameters such as the Reynolds-stress anisotropy and total strain is proposed, and is shown to outperform all other models for all cases examined, regardless of the mean strain rate.

Description: The phenomena of the Saturn system discovered by the Voyager missions are addressed. The characteristics of the rings, including grooves, spikes, and warps, are described, and the discovery of bombarded rings and young rings is discussed. The unique and unexplained properties of Enceladus are summarized, and the possibility that some moons of Saturn have been reassembled is addressed. The interaction of moons with the rings is examined, and the possibility that Titan's atmosphere may contain complex molecules that are the precursors of life is considered.

Description: Spectral data from earth observations have indicated the presence of N2 and CH4 on Triton. This paper outlines the use of the 1-D radiative-convective model developed for Titan to calculate the current pressure of N2 and CH4 on Triton. The production of haze material is obtained by scaling down from the Titan value. Results and predictions for the Voyager Triton encounter are as follows: A N2-CH4 atmosphere on Triton is thermodynamically self consistent and would have a surface pressure of approximately 50 millibar; due to the chemically produced haze, Triton has a hot atmosphere with a temperature of approximately 130 K; Triton's troposphere is a region of saturation of the major constituent of the atmosphere, N2.

Description: The process by which ambient gases can be implanted into silicates by shocks was investigated by analyzing the noble-gas content of several experimentally and naturally shocked silicate samples. The retentivity of shock-implanted gas during stepwise heating in the laboratory was defined in terms of two parameters, namely, the activation energy for diffusion and the extraction temperature at which 50 percent of the gas is released, both of which correlate with the shock pressure. The experiments indicate that, with increasing shock pressure, gas implantation occurs through an increasing production of microcracks/defects in the silicate lattice. The degree of annealing of these defects control the degree of diffusive loss of implanted gas.

Description: The long-term modulation of Saturn's nonthermal radio emission in the kilometric wavelength range has been studied based upon data obtained by Voyagers 1 and 2. A comparison of the ballistic and hydrodynamic propagation of solar wind features from the spacecraft to Saturn allows the uncertainty inherent in the projection to be determined. The results confirm the previous suggestion that momentum, ram pressure, and kinetic energy flux are the primary solar wind parameters that drive the nonthermal radio emission. It is suggested that, under certain conditions and for limited periods of time, the magnetic properties and time derivatives of the solar wind have increased importance.

Description: It has been observed from the plasma experiments on the Pioneer Venus Orbiter that the altitude of the upper boundary of the ionosphere decreases in response to increasing solar wind dynamic pressure. However, at pressures above about 2.5 x 10 to the -8th dynes/sq cm, the further decrease in the ionopause height is rather small. Following the model of Cloutier et al. (1969), it is suggested that during high solar wind conditions, when the ionopause is formed at lower altitudes, the solar wind induces vertical and horizontal flows which sweep away the ionospheric plasma that is produced locally by photoionization. As a result, a disturbed photodynamical ionosphere is formed which has the scale height of the ionizable neutral constituent. It is shown that such a photodynamical ionosphere is observed at the subsolar ionopause under these conditions. As a consequence of this interaction, the ionopause altitude is observed to follow the small-scale height of the ionizable species, atomic oxygen, showing only small changes with solar wind pressure.

Description: Asteroids appear in light of telescopic and meteority studies to be the most accessible repositories of early solar system history available. In the cooler regions of the outer asteroid belt, apparently unaffected by severe heating, the C, P, and D populations appear to harbor significant inventories of volatiles; the larger primordial belt population may have had an even greater percentage of volatile-rich, low-albedo asteroids, constituting a potent asteroid for veneering early terrestrial planet atmospheres. The volatile-rich asteroids contain carbon, structurally bound and adsorbed water, as well as remnants of interstellar material predating the solar system.

Description: The development and implementation of the Voyager missions are reviewed. The interplanetary missions preceding Voyager are discussed, focusing on the technological development leading up to the Voyager spacecraft. The main results from Voyager observations of Jupiter, Saturn, and Uranus are outlined. Also, consideration is given to the prospects for observations of Neptune.

Description: Isotopic analysis of nesquehonite recovered from the surface of the LEW 85320 H5 ordinary chondrite shows that the delta C-13 and delta O-18 values of the two generations of bicarbonate (Antarctic and Texas) are different: delta C-13 = + 7.9 per mil and + 4.2 per mil; delta O-18 = + 17.9 per mil and + 12.1 per mil, respectively. Carbon isotopic compositions are consistent with equilibrium formation from atmospheric carbon dioxide at - 2 + or - 4 C (Antarctic) and + 16 + or - 4 C (Texas). Oxygen isotopic data imply that the water required for nesquehonite precipitation was derived from atmospheric water vapor or glacial meltwater which had locally exchanged with silicates, either in the meteorite or in underlying bedrock. Although carbonates with similar delta C-13 values have been identified in the SNC meteorites EETA 79001 and Nakhla, petrographic and temperature constraints argue against their simply being terrestrial weathering products.

Description: Numerical studies of turbulent flow in an axisymmetric 45-deg-expansion combustor and bifurcated diffuser are presented. The Navier-Stokes equations incorporating a k-epsilon model were solved in a nonorthogonal curvilinear coordinate system. A zonal-grid method, where the flow field was divided into several subsections, was developed. This approach permitted different computational schemes to be used in the various zones. In addition, grid generation was made a more simple task. Boundary overlap and interpolating techniques were used, and an adjustment of the flow variables was required to assure conservation of mass flux. Three finite-differencing methods (hybrid, quadratic upwind, and skew upwind) were used to represent the convection terms. Results were compared with existing experimental data. In general, good agreement between predicted and measured values was obtained.

Description: A primary concern in the formation of a multinational Mars-exploration program is that of technology transfer, which would be most acute in the U.S./U.S.S.R. case but would exist in lesser degree among other nations. Another concern is that of management-complexity, which could inflate total costs and substantially counteract the anticipated benefits of spreading program costs among a number of nations. It is presently suggested that these problems can be substantially reduced by having each nation become responsible for one (or more) complete flight system. Each nation's role in mission operations must be clear, and the cleanliness of interfaces among flight systems must carry over into mission operations.

Description: Observations of the stratospheric IR emission structure on Saturn are presented. The high-spatial-resolution global images show a variety of new features, including a narrow equatorial belt of enhanced emission at 7.8 micron, a prominent symmetrical north polar hotspot at all three wavelengths, and a midlatitude structure which is asymmetrically brightened at the east limb. The results confirm the polar brightening and reversal in position predicted by recent models for seasonal thermal variations of Saturn's stratosphere.

Description: An overview of the Voyager 2 encounter with Neptune is presented, including a brief discussion of the trajectory, the planned observations, and highlights of the results described in the 11 companion papers. Neptune's blue atmosphere has storm systems reminiscent of those in Jupiter's atmosphere. An optically thin methane ice cloud exists near the 1.5-bar pressure level, and an optically thick cloud exists below 3 bars. Neptune's magnetic field is highly tilted and offset from the planet's center; it rotates with a period of 16.11 hours. Two narrow and two broad rings circle the planet; the outermost of these rings has three optically thicker arc segments. Six new moons were discovered in circular prograde orbits, all well inside Triton's retrograde orbit. Triton has a highly reflective and geologically young surface, a thin nitrogen atmosphere, and at least two active geyser-like plumes.

Description: A numerical procedure in which the Navier-Stokes equations are discretized using tightly coupled discretizations of pressure derivatives and continuity equations is used here to extend the range of known terminal velocities of gaseous bubbles in liquids well beyond that in previous investigations. Computations performed for Reynolds numbers up to 2000 and Marangoni numbers up to 1000 show only a modest variation of the scaled bubble velocity between 0.16 and 0.5. The bubble velocity is influenced more by the Marangoni number than by the Reynolds number.

Description: A numerical analysis is performed on thermocapillary buoyancy convection induced by phase change in a liquid droplet. A finite-difference code is developed using an alternating-direction implicit (ADI) scheme. The intercoupling relation between thermocapillary force, buoyancy force, fluid property, heat transfer, and phase change, along with their effects on the induced flow patterns, are disclosed. The flow is classified into three types: thermocapillary, buoyancy, and combined convection. Among the three mechanisms, the combined convection simulates the experimental observations quite well, and the basic mechanism of the observed convection inside evaporating sessile drops is thus identified. It is disclosed that evaporation initiates unstable convection, while condensation always brings about a stable density distribution which eventually damps out all fluid disturbances. Another numerical model is presented to study the effect of boundary recession due to evaporation, and the 'peeling-off' effect (the removal of the surface layer of fluid by evaporation) is shown to be relevant.

Description: Two-dimensional solidification influenced by anisotropic heat conduction has been considered. The interfacial energy balance was derived to account for the heat transfer in one direction (x or y) depending on the temperature gradient in both the x and y directions. A parametric study was made to determine the effects of the Stefan number, aspect ratio, initial superheat, and thermal conductivity ratios on the solidification rate. Because of the imposed boundary conditions, the interface became skewed and sometimes was not a straight line between the interface position at the upper and lower adiabatic walls (spatially nonlinear along the height). This skewness depends on the thermal conductivity ratio k(yy)/k(yx). The nonlinearity of the interface is influenced by the solidification rate, aspect ratio, and k(yy/k(yx).

Description: A multiple-time-scale turbulence model of a single point closure and a simplified split-spectrum method is presented. In the model, the effect of the ratio of the production rate to the dissipation rate on eddy viscosity is modeled by use of the multiple-time-scales and a variable partitioning of the turbulent kinetic energy spectrum. The concept of a variable partitioning of the turbulent kinetic energy spectrum and the rest of the model details are based on the previously reported algebraic stress turbulence model. Example problems considered include: a fully developed channel flow, a plane jet exhausting into a moving stream, a wall jet flow, and a weakly coupled wake-boundary layer interaction flow. The computational results compared favorably with those obtained by using the algebraic stress turbulence model as well as experimental data. The present turbulence model, as well as the algebraic stress turbulence model, yielded significantly improved computational results for the complex turbulent boundary layer flows, such as the wall jet flow and the wake boundary layer interaction flow, compared with available computational results obtained by using the standard kappa-epsilon turbulence model.

Description: The emission from a gray radiating medium is analyzed for transient cooling in surroundings at a low temperature. The medium is rectangular with no variations in the direction normal to the cross section. The integral equation for the transient temperature distribution is solved numerically using a two-dimensional Gaussian integration subroutine. The emissive ability for a rectangle at uniform temperature is compared with that for transient cooling where the temperature distribution of the region has reached a fully developed shape, as shown by a separation of variables solution. The two solutions provide the upper and lower bounds for the emittance of a rectangle during transient cooling. The emittances for various aspect ratios are presented as a function of the mean length of the rectangle and are compared with results for a plane layer and a cylinder.

Description: Calculations for n-decane drops evaporating in a spherical cluster surrounded by unvitiated ambient air at atmospheric pressure were performed using two previously proposed cluster models. Both cluster models predict that turbulent transport effects are more important in the case of small clusters. This is due to the smaller volume to surface ratio and thus to the greater transport of hot unvitiated gas to the drops in order to promote evaporation. The results obtained are compared with those of two turbulent models for each one of the 'trapping factors' and similarity models.

Description: The effects of critical layer nonlinearity are considered on spatially growing oblique instability waves on nominally two-dimensional shear layers between parallel streams. The analysis shows that three-dimensional effects cause nonlinearity to occur at much smaller amplitudes than it does in two-dimensional flows. The nonlinear instability wave amplitude is determined by an integro-differential equation with cubic type nonlinearity. The numerical solutions to this equation are worked out and discussed in some detail. The numerical solutions always end in a singularity at a finite downstream distance.

Description: The effects of critical-layer nonlinearity on spatially growing oblique instability waves on compressible shear layers between two parallel streams are considered. The analysis shows that mean temperature nonuniformities cause nonlinearity to occur at much smaller amplitudes than it does when the flow is isothermal. The nonlinear instability wave growth rate effects are described by an integrodifferential equation which bears some resemblance to the Landau equation, in that it involves a cubic-type nonlinearity. The numerical solutions to this equation are worked out and discussed in some detail. Inviscid solutions always end in a singularity at a finite downstream distance, but viscosity can eliminate this singularity for certain parameter ranges.

Description: A simulation is performed of a passive scalar field convected by a rapidly fluctuating velocity field whose correlation time approaches zero. By using a code proposed in a previous study (Chasnov et al., 1988), the turbulence spectrum of the passive temperature field in the conductive subrange is determined. A theoretical model is proposed which explains the result obtained by representing the transfer of scalar variance by an eddy conductivity, whose correlation time is limited by the correlation time of the velocity field.

Description: Recent models of the internal structure of Pluto and Charon, made possible by analysis of the Pluto/Charon mutual events are reviewed. At a mean density of just over 2 g/cu cm and a predicted rock mass fraction of roughly 0.7, the Pluto/Charon system is significantly rockier than the icy satellites of the giant planets, a contrast which may reflect its formation in a CO-rich outer solar nebula rather than a circumplanetary nebula. Pluto and Charon may in fact be so rocky that they lost volatiles early in their history (possibly during a Charon-forming impact event), although this is still an open issue.

Description: Preliminary results from electric field measurements in the environment of Mars using the plasma-wave system on board Phobos 2 are reported. Electron-plasma oscillations observed upstream of the bow shock correspond to a solar-wind density of 2/cu cm. The shock-foot boundary was crossed up to three times on each orbit. The shock ramp was detected at altitudes between 0.45 and 0.75 Mars radii R(M) above the planetary surface. The density increased by about a factor of two at the ramp. The shock position, although variable, seems to be consistent with previous measurements. The downstream magnetosheath contained broadband electric-field noise below the plasma frequency. The boundary of th obstacle, or plasmapause, was crossed at altitudes of the order of 0.28 R(M); the cold plasma density was highly variable within the planetopause and reached the unexpected value of 700/cu cm on the third orbit, at 0.25 R(M) altitude. Bursts of waves with frequencies below the electron cyclotron frequency occur within the planetopause.

Description: The magnetic fields of Mars have been measured from Phobos 2 with high temporal resolution in the tail and down to an 850-km altitude. During four successive highly elliptical orbits, the position of the bow shock as well as that of a transition layer, the 'planetopause', were identified. Subsequent circular orbits at 6000-km altitude provided the first high-resolution data in the planetary tail and indicate that the interplanetary magnetic field mainly controls the magnetic tail. Magnetic turbulence was also detected when the spacecraft crossed the orbit of Phobos, indicating the possible existence of a torus near the orbit of this moon.

Description: The formation of the ice-rock cores of the giant planets by density wave-assisted accretion is outlined. The process could be rapid (100,000-1,000,000 yr) and completed within the probable lifetime of the solar nebula. The mechanism works for both Jupiter and Saturn and does not require a large excess of mass over that believed present in their cores.

Description: This paper identifies and documents several well-preserved examples of Martian strike-slip faults and examines their relationships to wrinkle-ridges. The strike-slip faulting predates or overlaps periods of wrinkle-ridge growth southeast of Valles Marineris, and some wrinkle ridges may have nucleated and grown as a result of strike-slip displacements along the echelon fault arrays. Lateral displacements of several km inferred along these arrays may be related to tectonism in Tharsis.

Description: Results are presented on the long-term behavior of the main peak electron density in the Venus ionosphere during the solar cycle 21, based on 104 radio occultation measurements of the vertical electron density profile in the dayside ionosphere of Venus carried out aboard the Pioneer Venus Orbiter spacecraft (along with published data on 11 Venera 9-10 measurements). The equation representing the electron temperature at h = 140 km is presented. The results imply that the electron temperature at h = 140 km decreased by about 25 percent from solar maximum to solar minimum, compared to a decrease of 50-75 percent above 200 km found by Kliore and Mullen (1989).

Description: Theoretical models of early Martian atmospheric evolution describe the maintenance of a dense CO2 atmosphere and a warm, wet climate until the end of the heavy-bombardment phase of impacting. However, the presence of very young, earthlike fluvial valleys on the northern flank of Alba Patera conflicts with this scenario. Whereas the widespread ancient Martian valleys generally have morphologies indicative of sapping erosion by the slow outflow of subsurface water, the local Alba valleys were probably formed by surface-runoff processes. Because subsurface water flow might be maintained by hydrothermal energy inputs and because surface-runoff valleys developed late in Martian history, it is not necessary to invoke drastically different planet-wide climatic conditions to explain valley development on Mars. The Alba fluvial valleys can be explained by hydrothermal activity or outflow-channel discharges that locally modified the atmosphere, including precipitation and local overland flow on low-permeability volcanic ash.

Description: On the Mars rover sample-return mission, the rover vehicle will collect and select samples from different locations on the Martian surface to be brought back to earth for laboratory studies. It is anticipated that an in situ energy-dispersive X-ray fluorescence (XRF) spectrometer will be on board the rover. On such a mission, sample selection is of higher priority than in situ quantitative chemical anlaysis. With this in mind, a pattern recognition technique is proposed as a simple, direct, and speedy alternative to detailed chemical analysis of the XRF spectra. The validity and efficacy of the pattern recognition technique are demonstrated by the analyses of laboratory XRF spectra obtained from a series of geological samples, in the form both of standardized pressed pellets and as unprepared rocks. It is found that pattern recognition techniques applied to the raw XRF spectra can provide for the same discrimination among samples as a knowledge of their actual chemical composition.

Description: Long's self-similar vortex is known to have two solutions for each supercritical value of the flow force. Each of those solutions is shown to have a double structure if the flow force is large. The inertial instabilities of one of those large-flow-force limit solutions are investigated, showing that they are related to the instabilities of the Bickley jet in one regime. However, the swirl in the vortex becomes important for long waves, very strongly modifying the sinuous and varicose, Bickley modes. The asymptotic results obtained agree well with the numerical solutions for the sinuous mode, but not for the varicose mode, the difficulty in the latter case being apparently due to mode jumping.

Description: Numerical techniques are developed to solve the Navier-Stokes equations for unsteady incompressible flow. The extension of the finite-difference Galerkin (FDG) method of Stephens et al. (1984) to the continuous-time case in two or three space dimensions is explained, and the numerical implementation of the method is discussed with particular attention to the staggered-MAC-grid primitive-variable discretization, the application of discrete mass balance to avoid problems inherent in FDG schemes, the direct interpretation of the FDG expansion variables as a discrete streamfunction, and a mass-balance approach to two-dimensional problems with throughflow or obstacles. Numerical results are presented graphically for the evolution of asymptotic steady flow in a driven cavity at Reynolds number 400, 1000, or 3200; good agreement with published experimental data is demonstrated, with accurate predictions of secondary-vortex formation from wall bubble recirculations at Reynolds number 1000.

Description: The present photoelectric lightcurve observations of 44 Nysa at one aspect, and 64 Angelina at both polar aspects, indicate the presence of an opposition spike or brightening of approximately 0.25 mag within several deg of zero-phase angle. Since the three curves are identical within the bounds of observational scatter, the opposition spike cannot be judged anomalous within the asteroid taxonomic class E. The general similarity of the observed phase curves to those of the Uranian satellites and the rings of Saturn supports the status of the spike as an ordinary property of moderate-to-high albedo atmosphereless surfaces.

Description: In order to ascertain the features of organic compound-production in planetary atmospheres under the effects of plasmas and shocks, various mixtures of N2, CH4, and H2 modeling the atmosphere of Titan were subjected to discrete sparks, laser-induced plasmas, and UV radiation. The experimental results obtained suggest that UV photolysis from the plasma is an important organic compound synthesis process, as confirmed by the photolysis of gas samples that were exposed to the light but not to the shock waves emitted by the sparks. The thermodynamic equilibrium theory is therefore incomplete in the absence of photolysis.

Description: An account is given of the results of photometric light curve observations in the V band for asteroids, yielding estimates of the mean and maximum reduced magnitudes for each object. On the basis of fits of the H-G magnitude relation for 33 objects, the mean values of the slope parameter are examined for different taxonomic classes; these values are then applied to the analysis of less complete data sets. While in the case of the moderate albedo S and M class asteroids the H-G relation appears to fit the available data, the relation for dark asteroids appears to predict more of an opposition than is characteristically present.

Description: Flow-field measurements of unsteady turbulent flow downstream of a rotating spoked-wheel wake generator were performed in a short-duration light-piston tunnel, and the instantaneous-velocity data were phase averaged based on a signal synchronized with the bar-passing frequency. Mean axial velocities were found to agree well with those obtained from measurements behind a stationary cylinder and to be independent of both Reynolds and bar-passing Strouhal numbers. Reynolds stresses were found to be consistent with related cylinder-wake measurements, but were significantly higher than corresponding measurements obtained in large-scale research turbomachines. Phase-averaged triple velocity correlations were calculated from the digital velocity records, revealing the sign and the magnitude of skewness in the velocity probability density distributions for the two components.

Description: The bifurcation diagram corresponding to the Eckhaus stability curve has been constructed for the one-dimensional Swift-Hohenberg equation in a finite domain. Finite-amplitude solutions with particular spatial wavelength recover linear stability, as predicted by the Eckhaus curve, after a sequence of secondary bifurcations from the branch of solutions with this wavelength. No connectivity between the primary-solution branches is admissible if the stability predicted by this bifurcation diagram is to correspond to the prediction of the Eckhaus analysis. The Eckhaus curve does not exist if nonlinear couplings destroy this pattern. This is demonstrated by analysis of a coupled pair of Swift-Hohenberg equations.

Description: Voyager imaging, radio occultation, and stellar occultation data for the regular structure of Saturn's inner Cassini Division are presently analyzed. The regular optical depth variation observed by the radio occultation experiment scan and the feature noted in Voyager images is the same structure, namely the gravitational wakes of two 10-km radius satellites orbiting within the division. The structure is azimuthally symmetric, and is judged to rule out the possibility that large moonlets may be responsible for the observed structure.

Description: TEM and IR spectroscopy investigations of the interplanetary dust particles (IDPs) collected in the stratosphere have shown the majority of IDPs in the layer-lattice silicate and pyroxene classes to not have been heated to temperatures above 600 C during atmospheric entry. This implies that they arrive at the upper atmosphere with low geocentric encounter velocities, and limits the possible encounter trajectories for these particles to relatively circular prograde orbits. On this basis, it is judged unlikely that these IDPs are from earth-crossing comets or asteroids; collected IDPs dominated by olivine include a larger portion of above-600 C-heated particles, suggesting their capture from more eccentric orbits.

Description: The present examination of the west Deuteronilus Mensae region of Mars notes the changes in fretted terrain across the gradational boundary from uplands to lowlands to include a reduction of canyon wall slopes and depths, so that the fretted terrain north of the gradational boundary appears matted, but not obscured. The two process-classes that may account for the lateral overlap are the erosion of stratified upland terrain, and the deposition of plains materials onto the sloping upland margin and fretted terrain. A variety of plains-emplacement mechanisms is considered.

Description: The compositional contrast between the giant-planet satellites and the significantly rockier Pluto/Charon system is indicative of different formation mechanisms; cosmic abundance calculations, in conjunction with an assumption of the Pluto/Charon system's direct formation from solar nebula condensates, strongly suggest that most of the carbon in the outer solar nebula was in CO form, in keeping with both the inheritance from the dense molecular clouds in the interstellar medium, and/or the Lewis and Prinn (1980) kinetic-inhibition model of solar nebula chemistry. Laboratory studies of carbonaceous chondrites and Comet Halley flyby studies suggest that condensed organic material, rather than elemental carbon, is the most likely candidate for the small percentage of the carbon-bearing solid in the outer solar nebula.

Description: An attempt is made to provide a constraint on the combined mass of Janus and Epimetheus from an analysis of Voyager I and Voyager 2 data and ground-based observations obtained during the 1966 and 1980 ring plane crossings. The results of the analysis presented here suggest that the total mass is 2.59 + or - 0.26 x 10 to the 21st g, the mass ratio is 3.61 + or - 0.01, and Janus' density is 0.67 + or - 0.10 g/cu cm. The low density of Janus is attributed to its porosity rather than composition.

Description: A coordinated multinational Mars Surveyor Program involving a series of robotic missions is proposed in order to select worthwhile human landing sites for discoveries, for safe operations, and for testing and proving technologies for making human flight to Mars possible. Some characteristics motivating exploration of the planet are briefly discussed, including the possibility of life, geological features, and meteorological conditions. The necessity for preliminary exploration of Mars by robots prior to human exploration is discussed, and the rationale behind a multinational approach for a Mars Surveyor Program is presented.

Description: The effects of asymmetry in furnace temperature profile and pulling velocity on the crystal interface shape are demonstrated while neglecting the latent heat of solidification. It is seen that the furnace temperature profile may be varied in order to influence the shape of the melt-crystal interface. An exact thermal analysis is then performed on the Bridgman technique by including the latent heat of solidification as a source term. The exact temperature field required for yielding a flat melt-crystal interface is obtained. The earlier observation regarding the influence of furnace temperature profile on the interface shape is confirmed and a criterion for achieving a flat interface is obtained. Various furnace temperature profiles are selected and their corresponding melt-crystal interface results are presented.

Description: Arecibo high-resolution (1.5 to 2 km) radar data of Venus for the area extending from Beta Regio to western Eisila Regio provide strong evidence that the mountains in Beta and Eisila Regiones and plains in and adjacent to Guinevere Planitia are of volcanic origin. Recognized styles of volcanism include large volcanic edifices on the Beta and Eisila rises related to regional structural trends, plains with multiple source vents and a mottled appearance due to the ponding of volcanic flows, and plains with bright features surrounded by extensive quasi-circular radar-dark halos. The high density of volcanic vents in the plains suggests that heat loss by abundant and widely distributed plains volcanism may be more significant than previously recognized. The low density of impact craters greater than 15 km in diameter in this region compared to the average density for the higher northern latitudes suggests that the plains have a younger age.

Description: In February and March 1989 the Phobos 2 spacecraft took 37 television images of Phobos from a distance of 190-1100 km. The data are being used to update the three-dimensional model of Phobos, to provide improved determinations of its density and orbital dynamics, and to study its surface color, composition, and texture. Preliminary findings are presented here which include different integrated photometric behavior in visible and near-infrared bands, observation of a region immediately west of Stickney which is relatively free of large grooves, the prevalence of bright rims on grooves and younger craters, and low bulk density.

Description: Transient cooling by radiation is analyzed for a cylindrical region filled with axially flowing streams of drops that are becoming solidified. This is of interest for the dissipation of waste heat from orbiting power system in space. The drops absorb, emit, and scatter radiation, and the surroundings are at a lower uniform temperature. The radiative properties are assumed gray, and the scattering is isotropic. The radiating region is a two-phase mixture that remains at the melting temperature of the drops. Its temperature uniformity maintains a high emissive power as energy is lost. This is an advantage over a sensible heat radiator in which the temperature decreases, thereby reducing the emissive power. The results provide the axial length that remains two-phase and the fraction of energy dissipated within this length in which the emissive power has not decreased because of sensible cooling. It is also shown how the radial distribution of the axial velocity of the drops can be modified to increase this energy fraction.

Description: A number of successful applications of a spectral collocation method extended by a multi-domain patching technique are shown. The multi-domain technique can be used to improve resolution for problems with widely disparate scales, and to reduce the ill-conditioning of the spectral operators for problems in which a large number of points are required for distributed resolution. A new nonreflecting outflow boundary treatment for unsteady transition-to-turbulence simulations is also presented, which relies on the multi-domain technique. The role of multi-domain in improving the efficiency of such calculations is discussed.

Description: Spectral element methods are high-order weighted residual techniques based on spectral expansions of variables and geometry for the Navier-Stokes (NS) and transport equations. Here, practical aspects of these methods and their efficient implementation are examined, and several examples of flows in truly complex geometries are presented. The spectral element discretization for NS equations is introduced, and the convergence of the method is addressed. An efficient data management scheme is discussed in the context of parallel processing computations. The method is validated by comparing the spectral element solutions with the exact eigensolutions for the Orr-Sommerfeld equations in two and three dimensions. Computer-aided flow visualizations are presented for an impulsive flow past a sharp edge wedge. Three-dimensional states of channel flow disrupted by an array of cylindrical eddy promoters are studied, and the results of a direct simulation of the turbulent flow in a plane channel are presented.

Description: The nearly ubiquitous presence of a negative Eu anomaly in the mare basalts has been suggested to indicate prior separation and flotation of plagioclase from the basalt source region during its crystallization from a lunar magma ocean (LMO). Are there any mare basalts derived from a mantle source which did not experience prior plagioclase separation? Crystal chemical rationale for REE substitution in pyroxene suggests that the combination of REE size and charge, M2 site characteristics of pyroxene, fO2, magma chemistry, and temperature may account for the negative Eu anomaly in the source region of some types of primitive, low TiO2 mare basalts. This origin for the negative Eu anomaly does not preclude the possibility of the LMO as many mare basalts still require prior plagioclase crystallization and separation and/or hybridization involving a KREEP component.

Description: Pioneer Venus ion composition measurements are used to study the Venus ionosphere during solar minimum. It is suggested that the topside electron density profile at Venus during solar minimum has two distinct regimes. One beween 140 and 180 km is dominated by O2(+) ions which are in photochemical equilibrium. The other regime is above 180 km and is dominated by O(+) ions which are disturbed by the solar wind induced plasma transport. For Pioneer Venus, Mariner 10, and Venera 9 and 10 data, it is found that Venus exhibits a photodynamical type of ionopause during solar minimum.

Description: Voyager 2's IR observations of Neptune encompass thermal emissions and broadband radiometer measurements of reflected solar radiation. Temperature maps were obtained for the planet between 80 deg S and 30 deg N for two atmospheric layers, one in the lower stratosphere and the other in the troposphere. The relatively warm pole and equator, with cooler midlatitudes, are qualitatively similar to Uranus, despite the two planets' very different obliquities and internal heat fluxes. Powerful wavelike longitudinal thermal structure is noted, of which some appears to be associated with the Great Dark Spot; a localized cold region uncorrelated with any visible feature is found in the lower stratosphere.

Description: Detection of very intense short radio bursts from Neptune was possible as early as 30 days before closest approach and at least 22 days after closest approach. The bursts lay at frequencies in the range 100 to 1300 kilohertz, were narrowband and strongly polarized, and presumably originated in southern polar regions of the planet. Episodes of smooth emissions in the frequency range from 20 to 865 kilohertz were detected during an interval of at least 10 days around closest approach. The bursts and the smooth emissions can be described in terms of rotation in a period of 16.11 + or - 0.05 hours. The bursts came at regular intervals throughout the encounter, including episodes both before and after closest approach. The smooth emissions showed a half-cycle phase shift between the five episodes before and after closest approach. This experiment detected the foreshock of Neptune's magnetosphere and the impacts of dust at the times of ring-plane crossings and also near the time of closest approach. Finally, there is no evidence for Neptunian electrostatic discharges.

Description: Based on the success of several two-dimensional (latitude, longitude) linear barotropic instability models at matching some of the observed characteristics of the cloud level, polar region of the Venus atmosphere, a more realistic, linear, three-dimensional (height, latitude and longitude) model has been developed to further test the hypothesis that the observed features can be described by linear instability theory. The approach taken is to vary the model input parameters to see whether it is possible to produce modes that resemble the observations of wave activity and to compare those input parameters with other observations of the mean state. Sensitivity studies show that in addition to a well-documented dependence on the mean zonal wind, the growth and propagation of unstable modes depends on the latitude variation of the mean temperature (and hence static stability). These studies have led to the specification of a model basic state wind and temperature field that produces modes which are matched to observations of spatial structure, preferred wavenumber and phase speed of the polar disturbances. Wavenumber 2 is found to have the shortest growth time and unlike the two-dimensional results, wavenumbers 1-3 share a nearly common period of about 3 days. The derived basic state has a temperature field that is quite similar to Pioneer Venus observations; however, in some regions the model basic state wind field departs from cyclostrophic values based on temperature observations.

Description: A description of the Axial Flow Turbine Research Facility (AFTRF) being built at the Turbomachinery Laboratory of the Pennsylvania State University is presented. The purpose of the research to be performed in this facility is to obtain a better understanding of the rotor/stator interaction, three dimensional viscous flow field in nozzle and rotor blade passages, spanwise mixing and losses in these blade rows, transport of wake through rotor passage, and unsteady aerodynamics and heat transfer of rotor blade row. The experimental results will directly feed and support the analytical and the computational tool development. This large scale low speed facility is heavily instrumented with pressure and temperature probes and has provision for flow visualization and laser Doppler anemometer measurement. The facility design permits extensive use of the high frequency response instrumentation on the stationary vanes and more importantly on the rotating blades. Furthermore it facilitates detailed nozzle wake, rotor wake, and boundary layer surveys. The large size of the rig also has the advantage of operating at Reynolds numbers representative of the engine environment.

Description: A quasi-three-dimensional analysis has been developed for unsteady rotor-stator interaction in turbomachinery. The analysis solves the unsteady Euler or thin-layer Navier-Stokes equations in a body-fitted coordinate system. It accounts for the effects of rotation, radius change, and stress-surface thickness. The Baldwin-Lomax eddy-viscosity model is used for turbulent flows. The equations are integrated in time using an explicit four-stage Runge-Kutta scheme with a constant time step. Implicit residual smoothing is used to increase the stability limit of the time-accurate computations. The scheme is described, and stability and accuracy analyses are given.

Description: A 3-D model was developed for simulating multistage turbomachinery flows using supercomputers. This average passage flow model described the time averaged flow field within a typical passage of a bladed wheel within a multistage configuration. To date, a number of inviscid simulations were executed to assess the resolution capabilities of the model. Recently, the viscous terms associated with the average passage model were incorporated into the inviscid computer code along with an algebraic turbulence model. A simulation of a stage-and-one-half, low speed turbine was executed. The results of this simulation, including a comparison with experimental data, is discussed.

Description: A multiphase turbulence closure model is presented which employs one transport equation, namely the turbulence kinetic energy equation. The proposed form of this equation is different from the earlier formulations in some aspects. The power spectrum of the carrier fluid is divided into two regions, which interact in different ways and at different rates with the suspended particles as a function of the particle-eddy size ratio and density ratio. The length scale is described algebraically. A mass/time averaging procedure for the momentum and kinetic energy equations is adopted. The resulting turbulence correlations are modeled under less retrictive assumptions comparative to previous work. The closures for the momentum and kinetic energy equations are given. Comparisons of the predictions with experimental results on liquid-solid jet and gas-solid pipe flow show satisfactory agreement.

Description: Generally, two types of theory are used to describe the field equations for suspensions. The so-called postulated equations are based on the kinetic theory of mixtures, which logically should give reasonable equations for solutions. The basis for the use of such theory for suspensions is tenuous, though it at least gives a logical path for mathematical arguments. It has the disadvantage that it leads to a system of equations which is underdetermined, in a sense that can be made precise. On the other hand, the so-called averaging theory starts with a determined system, but the very process of averaging renders the resulting system underdetermined. A third type of theory is proposed in which the kinetic theory of gases is used to motivate continuum equations for the suspended particles. This entails an interpretation of the stress in the particles that is different from the usual one. Classical theory is used to describe the motion of the suspending medium. The result is a determined system for a dilute suspension. Extension of the theory to more concentrated systems is discussed.

Description: The vaporization of a droplet, interacting with its neighbors in a non-dilute spray environment is examined as well as a vaporization scaling law established on the basis of a recently developed theory of renormalized droplet. The interacting droplet consists of a centrally located droplet and its vapor bubble which is surrounded by a cloud of droplets. The distribution of the droplets and the size of the cloud are characterized by a pair-distribution function. The vaporization of a droplet is retarded by the collective thermal quenching, the vapor concentration accumulated in the outer sphere, and by the limited percolative passages for mass, momentum and energy fluxes. The retardation is scaled by the local collective interaction parameters (group combustion number of renormalized droplet, droplet spacing, renormalization number and local ambient conditions). The numerical results of a selected case study reveal that the vaporization correction factor falls from unity monotonically as the group combustion number increases, and saturation is likely to occur when the group combustion number reaches 35 to 40 with interdroplet spacing of 7.5 diameters and an environment temperature of 500 K. The scaling law suggests that dense sprays can be classified into: (1) a diffusively dense cloud characterized by uniform thermal quenching in the cloud; (2) a stratified dense cloud characterized by a radial stratification in temperature by the differential thermal quenching of the cloud; or (3) a sharply dense cloud marked by fine structure in the quasi-droplet cloud and the corresponding variation in the correction factor due to the variation in the topological structure of the cloud characterized by a pair-distribution function of quasi-droplets.

Description: Particles in the mass range from 10 to the minus 7th power to 10 to the minus 3rd power grams contribute 80 percent of the total mass influx of meteoritic material in the 10 to the minus 13th power to 10 to the 6th power gram mass range at Earth (Hughes, 1978). On Earth atmospheric entry, all but the smallest particles in the 10 to the minus 7th power to 10 to the minus 3rd power gram mass range, about 60 to 1200 micrometers in diameter, are heated sufficiently to melt and vaporize. Mars, because of its lower escape velocity and larger atmospheric scale height, is a much more favorable site for unmelted survival of micrometeorites on atmospheric deceleration. Researchers calculate that a significant fraction of particles throughout the 60 to 1200 micrometer diameter range will survive atmospheric entry unmelted. Thus returned Mars soils may offer a resource for sampling micrometeorites in a size range uncollectable in unaltered form at Earth. The addition of meteoritic material to the Mars soils should perturb their chemical composition, as has been detected using the soils on the Moon (Anders, et al., 1973). Using measured mass influx at Earth and estimates of the Mars/Earth flux ratio, researchers estimate a mass influx at Mars of between 2,700 and 202,000 metric tons per year.

Description: The capability of accurate nonlinear flow analysis of resonance systems is essential in many problems, including combustion instability. Classical numerical schemes are either too diffusive or too dispersive especially for transient problems. In the last few years, significant progress has been made in the numerical methods for flows with shocks. The objective was to assess advanced shock capturing schemes on transient flows. Several numerical schemes were tested including TVD, MUSCL, ENO, FCT, and Riemann Solver Godunov type schemes. A systematic assessment was performed on scalar transport, Burgers' and gas dynamic problems. Several shock capturing schemes are compared on fast transient resonant pipe flow problems. A system of 1-D nonlinear hyperbolic gas dynamics equations is solved to predict propagation of finite amplitude waves, the wave steepening, formation, propagation, and reflection of shocks for several hundred wave cycles. It is shown that high accuracy schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time.

Description: An explicit multistage Runge-Kutta type of time-stepping scheme is used for solving transonic flow past a transport type wing/fuselage configuration. Solutions for both Euler and Navier-Stokes equations are obtained for quantitative assessment of boundary layer interaction effects. The viscous solutions are obtained on both a medium resolution grid of approximately 270,000 points and a find grid of 460,000 points to assess the effects of grid density on the solution. Computed pressure distributions are compared with the experimental data.

Description: Some considerations toward developing numerical procedures for simulating viscous compressible flows are discussed. Both Navier-Stokes and boundary layer field methods are considered. Because efficient viscous-inviscid interaction methods have been difficult to extend to complex 3-D flow simulations, Navier-Stokes procedures are more frequently being utilized even though they require considerably more work per grid point. It would seem a mistake, however, not to make use of the more efficient approximate methods in those regions in which they are clearly valid. Ideally, a general purpose compressible flow solver that can optionally take advantage of approximate solution methods would suffice, both to improve accuracy and efficiency. Some potentially useful steps toward this goal are described: a generalized 3-D boundary layer formulation and the fortified Navier-Stokes procedure.

Description: The construction and development of the multi-component traversing system and associated control hardware and software are presented. A hydrogen bubble/laser sheet flow visualization technique was developed to visually study the characteristics of the mixing layers. With this technique large-scale rollers arising from the Taylor-Gortler instability and its interaction with the primary Kelvin-Helmholtz structures can be studied.

Description: The sizes and arrangement of the wind tunnel used for the experimentation are described. The specifications for the cold-wire anemometers, hot-wire anemometers, cold-wire rakes, and miniature 3-wire probe are proveded. The results of the experiment are briefly discussed.

Description: A very low Reynolds number turbulent boundary layer subject to an adverse pressure gradient is studied. The aim is to obtain highly accurate mean-flow and turbulence measurements under conditions that can be closely related to the numerical simulations of Philippe Spalart for the purposes of CFD validation. Much of the Boundary Layer Wind Tunnel was completely rebuilt with a new wider contraction and working section which will improve compatibility with the simulations. A unique sophisticated high-speed computer controlled 3-D probe traversing mechanism was integrated into the test section. Construction of the tunnel and traverse is discussed in some detail. The hardware is now complete, and measurements are in progress. The mean-flow data indicate that a suitably two-dimensional base flow was established. Automation of the probe positioning and data acquistion have led to a decreased running time for total pressure measurements. However, the most significant benefits are expected to occur when using hot-wire probes. Calibrations can be performed automatically and there is no need to handle fragile probes when moving between measuring stations. Techniques are being developed which require sampling of the signals from moving hot-wire probes on the basis of their position in the flow. Measurements can be made in high intensity turbulence by flying probes upstream at high speed so that the relative magnitude of the turbulent velocity fluctuations are reduced. In regions, where the turbulence intensity is not too large, the probe can also be repetitively scanned across very dense spatial grids in other directions. With this technique, a complete profile can be measured in about 1/3 the time and with a spatial density about 50 times that obtainable using a stationary probe.

Description: Thermal convection was proposed as a possible mechanism for generation and maintenance of turbulence in the inner accretion disk regime of the primordial solar nebula. It is of fundamental interest to design experiments with the basic physical features of the solar nebula conditions cannot be produced in the laboratory, numerical simulations of hydrodynamic flows, which have been very successful in describing aerodynamic flows, can be suitable modified to provide experimental data for solar nebula modelling. The goals are to modify an extant, proven hydrodynamics code with the most important features of the solar nebula and other thin accretion disks: bouyancy terms to generate convection, internal heating representing the release of gravitational potential energy, a variable gravity linearly proportional the the distance from the vertical midplane due to centrifugal balance, rapid rotation with axis aligned with gravity, and Keplerian rotational shear; to determine the effect that these features have on the turbulent convection by introducing them individually and to determine the cumulative nature of the turbulent convection for accretion disk conditions; and to model the convection and the turbulence. In this manner, prior solar nebula models can be tested and their deficiencies rectified.

Description: It is likely that turbulence played a major role in the evolution of the solar nebula, which is the flattened disk of dust and gas out of which our solar system formed. Relevant turbulent processes include the transport of angular momentum, mass, and heat, which were critically important to the formation of the solar system. This research will break ground in the modeling of compressible turbulence and its effects on the evolution of the solar nebula. The computational techniques which were developed should be of interest to researchers studying other astrophysical disk systems (e.g., active galactic nuclei), as well as turbulence modelers outside the astrophysics community.

Description: A research program for direct numerical simulations of compressible reacting flows is described. Two main research subjects are proposed: the effect of pressure waves on turbulent combustion and the use of direct simulation methods to validate flamelet models for turbulent combustion. The interest of a compressible code to study turbulent combustion is emphasized through examples of reacting shear layer and combustion instabilities studies. The choice of experimental data to compare with direct simulation results is discussed. A tentative program is given and the computation cases to use are described as well as the code validation runs.

Description: Outline of the research program and a recent progress in the studies of sheared turbulence are described. The research program reported is directed at two goals: (1) understanding of fundamental physics of organized structures in turbulent shear flows; and (2) development of phenomenological models of turbulence based on physical arguments. Three projects that were carried out are: (1) structure of sheared turbulence near a plane boundary; (2) distortion of turbulence by axisymmetric strain and dilation; and (3) study of energy transfer in turbulent shear flow.

Description: Vortex interactions and their role in turbulent flow are examined. The objectives are twofold. First, to use the existing axisymmetric code to study the annihilation process of colliding vortex rings and determine the relevance of this problem to similar 3-D phenomena. The second objective is to extend the code to three dimensions. The code under development is unique in that it can compute flows in a truly infinite domain (i.e., without periodic boundary conditions or approximations from truncating the domain). Because of this, the far field sound can be computed, and therefore, contribute to improved models of turbulence generated noise for this class of flows. Issues which can be addressed by the code include: effects of viscosity on mode selection in azimuthal breakdown of vortex rings (i.e., the Widnall instability); reconnection, the associated production of small scales, and the time scale of the process.

Description: The motivation for studying close vortex interactions is briefly discussed in the light of turbulence and coherent structures. Particular attention is given to the interaction known as reconnection. Two reconnection mechanisms are discussed. One is annihilation of vorticity by cross-diffusion, the other is an inviscid head-tail formation. At intermediate Reynolds numbers both mechanisms are operating.

Description: The objective of the present research was to extend the Direct Numerical Simulation (DNS) approach to particle-laden turbulent flows using a simple model of particle/flow interaction. The program addressed the simplest type of flow, homogeneous, isotropic turbulence, and examined interactions between the particles and gas phase turbulence. The specific range of problems examined include those in which the particle is much smaller than the smallest length scales of the turbulence yet heavy enough to slip relative to the flow. The particle mass loading is large enough to have a significant impact on the turbulence, while the volume loading was small enough such that particle-particle interactions could be neglected. Therefore, these simulations are relevant to practical problems involving small, dense particles conveyed by turbulent gas flows at moderate loadings. A sample of the results illustrating modifications of the particle concentration field caused by the turbulence structure is presented and attenuation of turbulence by the particle cloud is also illustrated.

Description: Direct numerical simulations are being performed on two different fluid flows in an attempt to discover the mechanism underlying the transition to turbulence in each. The first system is Taylor-Couette flow; the second, two-dimensional flow over an airfoil. Both flows exhibit a gradual transition to high-dimensional turbulence through low-dimensional chaos. The hope is that the instabilities leading to chaos will be easier to relate to physical processes in this case, and that the understanding of these mechanisms can then be applied to a wider array of turbulent systems.

Description: Cocke (1969) proved that in incompressible, isotropic turbulence the average material line (material surface) elements increase in comparison with their initial values. Good estimates of how much they increase in terms of the eigenvalues of the Green deformation tensor were rigorously obtained.