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: In June 1978, the Seasat satellite was launched carrying, among other instruments, the Seasat-A scatterometer system (SASS), which produced ambiguous wind speed and direction data at the ocean surface. A fifteen day subset of dealiased wind vector data with the inherent ambiguities removed was produced for the period of September 6-20, 1978. On September 8, SASS began to observe a development of frontal cyclogenesis in the South Pacific off the east coast of New Zealand, in an area of few surface observations. A large mature cyclone contained weak warm and cold fronts and an occlusion with a strong horizontal wind shear. Satellite imagery shows that a strong upper-level jet streak was moving rapidly over the area of the surface frontal occlusion and as the jet passed over this area a new vortex formed. This cyclogenesis event was studied using 50-km resolution scatterometer surface wind data. High-resolution fields of wind vectors, divergence and vorticity are computed and plotted from the scatterometer data to study the structure and development of the newly formed cyclonic vortex, not otherwise possible using conventional observations.

Description: The NASA scatterometer (NSCAT) is a spaceborne scatterometer scheduled to be deployed in the mid-1990s. An analysis of the wind retrieval error distribution for wind estimates based on backscatter measurements made by the NSCAT instrument is presented. The results are based on an end-to-end simulation of the scatterometer instrument and data processing. In general, the distribution of the wind speed error, when normalized, is independent of the true wind speed and direction. The wind speed error can be characterized by a normal distribution. The wind direction error is independent of the true wind speed, but depends on the true wind direction. Details for wind vectors with true wind speeds from 3 m/s to 33 m/s and true wind directions from 0 to 360 deg are presented.

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: Modulation of a rain wave pattern by longer waves has been studied. An analytical model taking into account capillarity effects and obliquity of short waves has been developed. Modulation rates in wave number and amplitude have been computed. Experiments were carried out in a wave tank. First results agree with theoretical models, but higher values of modulation rates are measured. These results could be taken into account for understanding the radar response from the sea surface during rain.

Description: Rainfall modification of directional scatterometer response from the sea surface was simulated in wind-wave tank experiments. Data show that for the range of conditions in laboratory experiments, rain enhances radar cross section for all azimuthal angles relative to wind direction. This result broadens previous measurements, which showed that scatterometer response increases with increasing rainfall for radars pointing upwind. But more to the point, the data also show that the directional dynamic-range of scatterometry diminishes rapidly as rainfall rate increases. Thus, while it may be possible to determine wind speed and direction during rain, it will require adequate system sensitivity.

Description: An account of the construction of surface pressure fields from Seasat-A satellite scatterometer (SASS) winds as carried out by different methods, and the comparison of these pressure fields with those derived from in situ ship observations is presented. On the assumption that the pressure adjusts itself instantaneously to the motion field, it may be computed by various methods. One of these makes use of planetary boundary theory, and of the possible techniques in this category a two-layer iterative scheme admitting of the parametrization of diabatic and baroclinic effects and of secondary flow was chosen. A second method involves the assumption of zero two-dimensional divergence, leading to a Laplace's equation (the balance equation) in pressure, with the wind field serving as a forcing function. This method does not accommodate adiabatic or baroclinic effects, and requires a knowledge of the pressure at all boundary points. Two comparison fields are used for validation: the conventional operational analyses of the US National Meteorological Center (NMC), and the special analyses of the Gulf of Alaska Experiment (GOASEX), which were done by hand. The results of the computations were as follows: (1) The pressure fields, as computed from the SASS winds alone, closely approximated the NMC fields in regions where reasonable in situ coverage was available (typically, one or two mb differences over most of the chart, three to four mb in extreme cases); (2) In some cases the SASS-derived pressure fields displayed high-resolution phenomena not detected by the NMC fields, but evident in the GOASEX data; and, (3) As expected, the pressure fields derived from the balance equation were much smoother and less well resolved than the SASS-derived or NMC fields. The divergence as measured from the SASS winds is smaller than, but of the same order of magnitude as, the vorticity.

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: The accuracy of temperature, pressure, potential temperature, and horizontal wind measurements is discussed in connection with the use of Meteorological Measurement System data in the AAOE. The vertical distribution of temperature measurements and latitudinal variations of the zonal wind for 12 flights over Antarctica during the 1987 AAOE campaign are summarized. Model atmospheres from 0 to 32 km at 70 deg and 55 deg S for the August-September period are constructed. Above the 420 K isentropic surface, the polar vortex remains strong throughout August and September of 1987.

Description: The formulation and evolution of polar stratospheric ice clouds are simulated using a one-dimensional model of cloud microphysics. It is found that the optical thickness and particle size of ice clouds depend on the cooling rate of the air in which the cloud formed. It is necessary that there be an energy barrier to ice nucleation upon the preexisting aerosols in order to account for the cooling rate dependence of the cloud properties.

Description: The history of minimum temperatures at 50 and 70 mb is examined from NMC, UK Met O and ECMWF analyses. MSU channel 24 data are similarly inspected. South Pole sonde data are used to calculate saturation humidity mixing ratio as a function of altitude and time throughout 1987. Saturation with respect to ice could be maintained for water mixing ratios of 3.5 ppmv for a period of about 80 days from mid-June to mid-September. Dehydration to mixing ratios of 1 ppmv or less was possible sporadically. Data from the ER-2 flights between 53 S and 72 S are used in conjunction with particle size measurements and air parcel trajectories to demonstrate the dehydration occurring over Antarctica. Water mixing ratios at the latitude of Punta Arens (53 S), in conjunction with tracer measurements and trajectory analysis, show that at potential temperatures from about 325 to 400 K, the dryness (less than 3 ppmv) had its origin over Antarctica rather than in the tropics. Water mixing ratios within the Antarctic vortex varied from 1.5 to 3.8 ppmv, with a strong isentropic gradient being evident in the region of high potential vorticity gradients.

Description: Measurements of total reactive nitrogen, NOy, total water vapor, and aerosols were made as part of the Airborne Antarctic Ozone Experiment. The measurements were made using instruments located onboard the NASA ER-2 aircrafts which conducted twelve flights over the Antarctic continent reaching altitudes of 18 km at 72 S latitude. Each instrument utilized an ambient air sample and provided a measurement up to 1 Hz or every 200 m of flight path. The data presented focus on the flights of Aug. 17th and 18th during which Polar Stratospheric Clouds (PSCs) were encountered containing concentrations of 0.5 to 1.0 micron diameter aerosols greater than 1 cm/cu. The temperature pressure during these events ranged as low as 184 K near 75 mb pressure, with water values near 3.5 ppm by volume (ppmv). With the exception of two short periods, the PSC activity was observed at temperatures above the frost point of water over ice. The data gathered during these flights are analyzed and presented.

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: A multichannel statistical approach is used to retrieve rainfall rates from the brightness temperature T(B) observed by passive microwave radiometers flown on a high-altitude NASA aircraft. T(B) statistics are based upon data generated by a cloud radiative model. This model simulates variabilities in the underlying geophysical parameters of interest, and computes their associated T(B) in each of the available channels. By further imposing the requirement that the observed T(B) agree with the T(B) values corresponding to the retrieved parameters through the cloud radiative transfer model, the results can be made to agree quite well with coincident radar-derived rainfall rates. Some information regarding the cloud vertical structure is also obtained by such an added requirement. The applicability of this technique to satellite retrievals is also investigated. Data which might be observed by satellite-borne radiometers, including the effects of nonuniformly filled footprints, are simulated by the cloud radiative model for this purpose.

Description: This paper presents the results of a field program using a ground-based Raman lidar system to observe changes in moisture profiles as a cold and a warm front passed over the NASA/Goddard Space Flight Center in Greenbelt, Maryland. The lidar operating only during darkness is capable of providing continuous high vertical resolution profiles of water vapor mixing ratio and aerosol scattering ratio from near the surface to about 7 km altitude. The lidar data acquired on three consecutive nights from shortly after sunset to shortly before sunrise, along with upper air data from specially launched rawinsondes, have provided a unique visualization of the detailed structure of the two fronts.

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: An improved version of the GISS Model II cumulus parameterization designed for long-term climate integrations is used to study the effects of entrainment and multiple cloud types on the January climate simulation. Instead of prescribing convective mass as a fixed fraction of the cloud base grid-box mass, it is calculated based on the closure assumption that the cumulus convection restores the atmosphere to a neutral moist convective state at cloud base. This change alone significantly improves the distribution of precipitation, convective mass exchanges, and frequencies in the January climate. The vertical structure of the tropical atmosphere exhibits quasi-equilibrium behavior when this closure is used, even though there is no explicit constraint applied above cloud base.

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: The joint frequency distribution technique was used to analyze buoyancy fluxes in the marine atmospheric boundary layer (MABL) for the cloud street regime noted during the Genesis of Atlantic Lows Experiment. It is found that for the lower half of the MABL, the buoyancy flux is mainly generated by the rising thermals and the sinking compensating ambient air, and is mainly consumed by the entrainment and detrainment of thermals, penetrative convection, and the entrainment from the MABL top. If the buoyancy flux is primarily driven by the temperature flux, these buoyancy-flux generating processes should be the same for the lower boundary layers over land and ocean. The results of the scale analysis of the buoyancy flux agree well with those obtained for mesoscale cellular convection during the Air-Mass Transformation Experiment.

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: Results are presented on numerical experiments that were carried out to investigate the mechanisms of the observed variabilities in wind and convection associated with supercloud clusters (SCCs), westerly wind bursts, and 30-60 day oscillations in the western Pacific region. It is shown that the generation of a 30-60 day eastward propagating precipitation pattern in the Lau and Peng (1987) model, which can be identified as SCC, is accompanied by convective clusters coming in opposite direction to that of the SCC itself. The results suggest that the westward propagating cloud clusters are produced at the initial stage of the 30-60 day disturbance due to mutual adjustment of the large-scale flow and heating.

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: Multidecadal time series of surface winds from central tropical Pacific islands are used to compute trends in the trade winds between the end of WWII and 1985. Over this period, averaged over the whole region, there is no statistically significant trend in speed or zonal or meridional wind (or pseudostress). However, there is some tendency, within a few degrees of the equator, toward weakening of the easterlies and increased meridional flow toward the equator. Anomalous conditions subsequent to the 1972-73 ENSO event make a considerable contribution to the long-term trends. The period 1974-80 has been noted previously to have been anomalous, and trends over that period are sharply greater than those over the longer records.

Description: This paper presents an analysis of the uncertainties expected in vertical velocities using a vertically pointing airborne Doppler radar which has a nadir or zenith-pointing beam. To examine the expected uncertainty, the Doppler velocity equation for a moving platform is derived and it is applied to cases of nadir-fixed and stabilized beams. The main emphasis of the paper is on the effect of platform stability on the deduced vertical air motions and it is shown that the antenna must be stabilized to obtain desired accuracy in the vertical velocity measurements.

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: Monthly fields of shortwave radiation (SR) and latent heat flux (LE) over the central and eastern tropical Pacific between 1980 and 1983 have been computed using satellite data. They are the dominant variable components of surface thermal forcing on the ocean in this time scale. During the 1982-1983 ENSO episode, surface-wind convergence and cloudiness associated with the displacement of equatorial organized convection caused a reduction in both the SR into the ocean and the LE out of the ocean. The lag-correlation coefficients between the forcing (SR-LE) and the sea surface temperature are found to be significantly high outside the equatorial region, showing that surface thermal forcing is the dominant factor in sea surface temperature change. In the narrow equatorial wave guide, ocean dynamics play a more important role, and surface heat flux is a consequence rather than the cause of sea surface temperature change.

Description: An aircraft experiment has been conducted with a dual-frequency (10 GHz and 35 GHz) radar/radiometer system and an 18-GHz radiometer to test various rain-rate retrieval algorithms from space. In the experiment, which took place in the fall of 1988 at the NASA Wallops Flight Facility, VA, both stratiform and convective storms were observed. A ground-based radar and rain gauges were also used to obtain truth data. An external radar calibration is made with rain gauge data, thereby enabling quantitative reflectivity measurements. Comparisons between path attenuations derived from the surface return and from the radar reflectivity profile are made to test the feasibility of a technique to estimate the raindrop size distribution from simultaneous radar and path-attenuation measurements.

Description: Azimuthal response of a scatterometer to radiation scattered by the sea surface was studied in a wind-wave tank. The variation of the normalized radar cross section with the azimuth angle is fitted by a three-term series. Results show that the upwind-downwind asymmetry decreases as the wind speed increases. The crosswind modulation depends on the wind velocity. The results show that the evolution of the long-wind-crosswind ratio evolves with wind speed in a manner similar to the evolution of the isotropy of short capillary-gravity waves. The maximum of the isotropy of the short wind waves is obtained for wind velocities close to 4 m/s. For the same value of the velocity, the variations of radar response between long-wind and crosswind directions is minimum. For lower or higher values of wind velocities the directional accuracy of the radar increases, since the wind-wave field tends to align in the wind direction.

Description: The present interpretation of the radar cross section sigma exp 0 measured by satellite altimeters implies that the rms wave slope gamma is controlled solely by the local wind. However, parameters of wave spectra, including the exponent in the power law for the equilibrium range, depend on sea maturity. The latter is characterized by the nondimensional fetch, x = gX/U-squared. Consequently, gamma and sigma exp 0 are controlled by both U and the wind fetch X. Geosat data for one year are used jointly with in-situ wind and wave observations to assess the fetch-related error trend in altimeter wind speeds. This trend results in overestimated winds in the regions and seasons characterized by a high x, and vice versa. A procedure for wind speed retrieval based on processing sigma exp 0 jointly with the significant wave height information contained in the altimeter wave forms is proposed.

Description: Seasat scatterometer data over the Arabian Sea are used to build wind-stress fields during July and August 1978. They are first compared with 3-day wind analyses from ship data along the Somali coast. Seasat scatterometer specifications of 2-m/s and 20-deg accuracy are fulfilled in almost all cases. The exceptions are for winds stronger than 14 m/s, which are underestimated by the scatterometer by 15 percent. Wind stress is derived from these wind data using a bulk formula with a drag coefficient depending on the wind intensity. A successive-correction objective analysis is used to build the wind-stress field over the Arabian Sea with 2 x 2-deg and 6-day resolution. The final wind-stress fields are not significantly dependent on the objective analysis because of the dense coverage of the scatterometer. The combination of scatterometer and coastal ship data gives the best coverage to resolve monsoon wind structures even close to the coast. The final wind stress fields show wind features consistent with other monthly mean wind stress field. However, a high variability is observed on the 6-day time scale.

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: Cumulus processes involved in the interaction and merging of clouds under the influence of different imposed conditions (including large-scale lifting forcing, environmental wind shear, and cloud microphysical processes) were studied using simulations with a three-dimensional model. The design of the study was to generate several convective clouds randomly inside the model domain, and then to observe and analyze the interactions and merging between the simulated clouds. Ten merged clouds were identified. Seven of these, each involving two previously separated clouds, generally lie along a line parallel to the initial environmental wind shear vector, while one (also a two-cloud system) lies along a line perpendicular to the wind shear vector prior to merging. The remaining two merging systems involve three parent clouds each; they are a combination of parallel and perpendicular cells. The merging mechanisms associated with three-cloud merging cases are studied by examining the temperature, pressure, and wind fields prior to, during, and following the merging.

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: The possibility of initiating the growth of ice sheets by solar insolation variations is examined. The study is conducted using a climate model with three different orbital configurations corresponding to 116,000 and 106,000 yr before the present and a modified insolation field with greater reductions in summer insolation at high northern latitudes. Despite the reduced summer and fall insolation, the model fails to maintain snow cover through the summer at locations of suspected ice sheet initiation. The results suggest that there is a discrepancy between the model's response to Milankovitch perturbations and the geophysical evidence of ice sheet initiation. If the model results are correct, the growth of ice shown by geophysical evidence would have occurred in an extremely ablative environment, demanding a complicated strategy.

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: Thirty-three Landsat TM scenes of California stratocumulus cloud fields were acquired as part of the FIRE Marine Stratocumulus Intensive Field Observations in July 1987. They exhibit a wide variety of stratocumulus structures. Analysis has so far focused upon the July 7 scene, in which aircraft from NASA, NCAR, and the British Meteorological Office repeatedly gathered data across a stratocumulus-fair weather cumulus transition. The aircraft soundings validate the cloud base temperature threshold determined by spatial coherence analysis of the TM thermal band. Brightness variations in the stratocumulus region exhibit a -5/3 power-law decrease of the wavenumber spectra for scales larger than the cloud thickness, about 200 m, changing to a -3 power at smaller scales. Observations by an upward-looking three-channel microwave radiometer on San Nicolas Island also show the -5/3 power-law in total integrated liquid water, suggesting that the largest-scale TM brightness variations are primarily due to variations in the liquid water. The Kolmogorov 5/3 power suggests that for some purposes liquid water in turbulent stratocumulus clouds may be treated as a passive scalar, simply reflecting variations in vertical velocity. This may be tested using the velocities measured by the aircraft.

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: New nonlinear partial differential equations containing terrain curvature and its rate of change are derived that describe the flow of an atmospheric density current. Unlike the classical hydraulic-type equations for density currents, the new equations are valid for two-dimensional, gradually varied flow over highly curved terrain, hence suitable for computing unsteady (or steady) flows over arbitrary mountain/valley profiles. The model assumes the atmosphere above the density current exerts a known arbitrary variable pressure upon the unknown interface. Later this is specialized to the varying hydrostatic pressure of the atmosphere above. The new equations yield the variable velocity distribution, the interface position, and the pressure distribution that contains a centrifugal component, often significantly larger than its hydrostatic component. These partial differential equations are hyperbolic, and the characteristic equations and characteristic directions are derived. Using these to form a characteristic mesh, a hypothetical unsteady curved-flow problem is calculated, not based upon observed data, merely as an example to illustrate the simplicity of their application to unsteady flows over mountains.

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 sources of sub-Saharan precipitation are studied using diagnostic procedures integrated into the code of the GISS climate model. Water vapor evaporating from defined source regions is 'tagged', allowing the determination of the relative contributions of each evaporative source to the simulated July rainfall in the Sahel. Two June-July simulations are studied to compare the moisture sources, moisture convergence patterns and the spatial variations of precipitation for rainy and drought conditions. Results for this case study indicate that patterns of moisture convergence and divergence over northern Africa had a stronger influence on model rainfall over the sub-Sahara than did evaporation rates over the adjacent oceans or moisture advection from ocean to continent. While local continental evaporation contributed significant amounts of water to Sahelian precipitation in the'rainy' simulation, moisture from the Indian Ocean did not precipitate over the Sahel in either case.

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.