ALBERT

Description: Basic algorithms for unstructured mesh generation and fluid flow calculation are discussed. In particular the following are addressed: preliminaries of graphs and meshes; duality and data structures; basic graph operations important in CFD (Computational Fluid Dynamics); triangulation methods, including Varonoi diagrams and Delaunay triangulation; maximum principle analysis; finite volume schemes for scalar conservation law equations; finite volume schemes for the Euler and Navier-Stokes equations; and convergence acceleration for steady state calculations.

Description: Transient solutions were obtained for a square region of heat conducting semitransparent material cooling by thermal radiation. The region is in a vacuum environment, so energy is dissipated only by radiation from within the medium leaving through its boundaries. The effect of heat conduction during the transient is to partially equalize the internal temperature distribution. As the optical thickness of the region is increased, the temperature gradients increase near the boundaries and corners, unless heat conduction is large. The solution procedure must provide accurate temperature distributions in these regions to prevent error in the calculated radiation losses. Two-dimensional numerical Gaussian integration is used to obtain the local radiative source term. A finite difference procedure with variable space and time increments is used to solve the transient energy equation. Variable spacing was used to concentrate grid points in regions with large temperature gradients.

Description: The 3.0-micrometers water of hydration absorption feature observed in the IR photometry of many low-albedo and some medium-albedo asteroids strongly correlates with the 0.7-micrometers Fe(+2) to Fe(+3) oxidized iron absorption feature observed in narrowband spectrophotometry of these asteroids. Using this relationship, an empirical algorithm for predicting the presence of water of hydration in the surface material of a Solar System body using photometry obtained through the Eight-Color Asteroid Survey nu (0.550 micrometers), w (0.701 micrometers), and x (0.853 micrometers) filters was developed and applied to the ECAS photometry of asteroids and outer planet satellites. The percentage of objects in low-albedo, outer main-belt asteroid classes that test positively for water of hydration increases from P to B to C to G class and correlates linearly with the increasing mean albedos of those objects testing positively. The medium-albedo M-class asteroids do not test positively in large number using this algorithm. Aqueously altered asteroids dominate the Solar System population between heliocentric distances of 2.6 to 3.5 AU, bracketing the Solar System region where the aqueous alteration mechanism operated most strongly. One jovian satellite, J VI Himalia, and one saturnian satellite. Phoebe, tested positively for water of hydration, supporting the hypothesis that these may be captured C-class asteroids from a postaccretional dispersion. The proposed testing technique could be applied to an Earth-based survey of asteroids or a space-probe study of an asteroid's surface characteristic in order to identify a potential water source.

Description: It is shown that to satisfy the general accepted compressible law of the wall derived from the Van Driest transformation, turbulence modeling coefficients must actually be functions of density gradients. The transformed velocity profiles obtained by using standard turbulence model constants have too small a value of the effective von Karman constant kappa in the log-law region (inner layer). Thus, if the model is otherwise accurate, the wake component is overpredicted and the predicted skin friction is lower than the expected value.

Description: The infrared transmission spectra and photochemical behavior of various organic compounds isolated in solid N2 ices, appropriate for applications to Triton ad Pluto, are presented. It is shown that excess absorption in the surface spectra of Triton and Pluto, i.e., absorption not explained by present models incorporating molecules already identified on these bodies (N2, CH4, CO, and CO2), that starts near 4450/cm (2.25 microns) and extends to lower frequencies, may be due to alkanes (C(n)H(2n+2)) and related molecules frozen in the nitrogen. Branched and linear alkanes may be responsible. Experiments in which the photochemstry of N2: CH4 and N2: CH4: CO ices was explored demonsrtrate that the surface ices of Triton and Pluto may contain a wide variety of additional species containing H, C, O, and N. Of these, the reactive molecule diazomethane, CH2N2, is particularly important since it may be largely responsible for the synthesis of larger alkanes from CH4 and other small alkanes. Diazomethane would also be expected to drive chemical reactions involving organics in the surface ices of Triton and Pluto toward saturation, i.e., to reduce multiple CC bonds. The positions and intrinsic strengths (A values) of many of the infrared absorption bands of N2 matrix-isolated molecules of relevance to Triton and Pluto have also been determined. These can be used to aid in their search and to place constraints on their abundances.

Description: The multigrid method has been applied to an existing three-dimensional compressible Euler solver to accelerate the convergence of the implicit symmetric relaxation scheme. This lower-upper symmetric Gauss-Seidel implicit scheme is shown to be an effective multigrid driver in three dimensions. A grid refinement study is performed including the effects of large cell aspect ratio meshes. Performance figures of the present multigrid code on Cray computers including the new C90 are presented. A reduction of three orders of magnitude in the residual for a three-dimensional transonic inviscid flow using 920 k grid points is obtained in less than 4 min on a Cray C90.

Description: The present paper explores the use of large-eddy simulations as a tool for predicting noise from first principles. A high-order numerical scheme is used to perform large-eddy simulations of a supersonic jet flow with emphasis on capturing the time-dependent flow structure representating the sound source. The wavelike nature of this structure under random inflow disturbances is demonstrated. This wavelike structure is then enhanced by taking the inflow disturbances to be purely harmonic. Application of Lighthill's theory to calculate the far-field noise, with the sound source obtained from the calculated time-dependent near field, is demonstrated. Alternative approaches to coupling the near-field sound source to the far-field sound are discussed.

Description: This report presents the most recent spherical harmonic topography model of Venus developed at Jet Propulsion Laboratory. It was produced by a spherical harmonic analysis of the most complete set of Magellan altimetry data, augmented by Pioneer Venus and Venera data. The harmonic coefficients of the topography were computed to degree and order 360. Compared to previous topography models, this one has the highest correlation with the gravity field of Venus.

Description: The 500-Myr average crater retention age for Venus has raised questions about the present-day level of tectonic activity. In this study we examine the relationship between the gravity and topography of four large volcanic swells, Beta, Atla, Bell, and Western Eistla Regiones, for clues about their stage evolution. The Magellan line-of-sight gravity data are inverted using a point mass model of the anomalous mass to solve for the local vertical gravity field. Spectral admittance calculated from both the local gravity inversions and a spherical harmonic model is compared to three models of compensation: local compensation, a 'flexural' model with local and regional compensation of surface and subsurface loads, and a 'hotspot' model of compensation that includes top loading by volcanoes and subsurface loading due to a deep, low density mass anomaly. The coherence is also calculated in each region, but yields an elastic thickness estimate only at Bell Regio. In all models, the long wavelengths are compensated locally. Our results may indicate a relatively old, possibly inactive plume.

Description: The steady state solution of the system of equations consisting of the full Navier-Stokes equations and two turbulence equations has been obtained using a multigrid strategy of unstructured meshes. The flow equations and turbulence equations are solved in a loosely coupled manner. The flow equations are advanced in time using a multistage Runge-Kutta time-stepping scheme with a stability-bound local time step, while turbulence equations are advanced in a point-implicit scheme with a time step which guarantees stability and positivity. Low-Reynolds-number modifications to the original two-equation model are incorporated in a manner which results in well-behaved equations for arbitrarily small wall distances. A variety of aerodynamic flows are solved, initializing all quantities with uniform freestream values. Rapid and uniform convergence rates for the flow and turbulence equations are observed.

Description: Recent mapping studies west of Elysium Mons, Mars, have pinpointed subice features that suggest the existence of a frozen paleolake in Utopia Planitia as recently as 1.8 billion years ago. The subice features are interpreted to be hyaloclastic ridges and hills, table moutains, associated joekulhalaup deposits, and fluvial channels. Photoclinometric studies of these features and of a basal scarp around the northwest flank of Elysium Mons interpreted to have been an ice-sheet boundary indicate that the maximum thickness of ice within the basin may have been about 180 m. This thickness of ice during a relatively late stage of Martian geologic history would have important implications concerning the atmospheric, the climatic, and possibly the exobiologic history of the planet.

Description: During the final, low solar activity phase of the Pioneer Venus (PV) mission, the Orbiter Ion Mass Spectrometer (OIMS) measurements found all ion species, in the midnight-dusk sector, reduced in concentration relative to that observed at solar maximum. Molecular ion species comprised a greater part of the total ion concentration as O(+) and H(+) had the greatest depletions. The nightside ionospheric states were strikingly similar to the isolated solar maximum 'disappearing' ionospheres. Both are very dynamic states characterized by a rapidly drifting plasma and 30-100 eV superthermal O(+) ions.

Description: In October, 1992 the Pioneer Venus Orbiter entered the atmosphere of Venus, ending nearly 14 years of observations at Venus. Prior to the entry into the atmosphere and subsequent loss of the spacecraft careful management of spacecraft resources had allowed the acquisition of much low altitude data over the nightside of the planet. The long duration of the Pioneer Venus mission has enabled us to study the ionosphere and atmosphere of Venus under different levels of solar activity.

Description: Atmospheric drag measurements from the orbital decay of the Pioneer Venus Orbiter and Magellan spacecraft have recently been obtained of the Venus dayside and nightside atmosphere between 130 and 210 km during a period of low solar activity. These new measurements, combined with the earlier Pioneer Venus drag measurements (1978-80) obtained near the maximum of the 11-year solar cycle, have allowed the detection of the detailed response of temperature, atomic oxygen and carbon dioxide to solar variations. We have found a weak but detectable temperature response on the dayside which is in accord with the response predicted by Keating and Bougher when they assumed very strong CO2 radiative cooling resulting from atomic oxygen exciting CO2 into 15 micron emission. This same radiative process may cause strong cooling in the Earth's upper atmosphere with the doubling of CO2 in the future. With decreasing solar activity, the O/CO2 ratio in the lower thermosphere is found to decrease, apparently due to decreased photodissociation of CO2 and lower temperatures. The percent decrease in atomic oxygen with decreasing solar activity on the dayside is found to be approximately the same as the percent decreases of atomic oxygen transported to the nightside. A very weak response of nightside temperatures to solar activity variations has also been detected.

Description: The ion and electron momentum equations, along with Ampere's law, are solved for the ion and electron drift velocities and the electric field in the subsolar Venus ionosphere, assuming a partially ionized gas and a single ion species having the ion mean mass. All collision terms among the ions, electrons and neutral particles are retained in the equations. A general expression for the evolution of the magnetic field is derived and compared with earlier expressions. Subsolar region data in the altitude range 150-300 km from the Pioneer Venus Orbiter are used to calculate altitude profiles of the components of the current due to the electric field, gradients of pressure, and gravity. Altitude profiles of the ion and electron velocities as well as the electric field, electrodynamic heating, and the energy density are determined. Only orbits having a complete set of measured plasma temperatures and densities, neutral densities, and magnetic field were considered for analysis; the results are shown only for orbit 202. The vertical velocity at altitudes above 220 km is upgoing for orbit 202. This result is consistent with observations of molecular ions at high altitudes and of plasma flow to the nightside, both of which require upward velocity of ions from the dayside ionosphere. Above about 230 km the momentum equations are extremely sensitive to the altitude profiles of density, temperature, and magnetic field.

Description: We report here analyses of olivines and pyroxenes, and petrofabrics of 27 chondritic interplanetary dust particles (IDPs), comparing those from anhydrous and hydrous types. Approximately 40% of the hydrous particles contain diopside, a probable indicator of parent body thermal metamorphism, while this mineral is rarely present in the anhydrous particles. Based on this evidence, we find that hydrous and anhydrous IDPs are, in general, not directly related, and we conclude that olivine and pyroxene major-element compositions can be used to help discriminate between IDPs that are (1) predominantly nebular condensates, and lately resided in anhydrous or icy (no liquids) primitive parent bodies, and (2) those originating from more geochemically active parent bodies (probably hydrous and anhydrous asteroids).

Description: In the present work, the high Reynolds number flow past an inclined plate with a splitter plate placed in its wake is considered numerically. A numerical conformal mapping technique is employed to transform the two-plate system into the same number of cylinders: the flow field is assumed to be two-dimensional. The vortex shedding from the inclined plate is modelled using the discrete vortex method. It is shown that the splitter plate has a profound effect on the development of the flow over a range of values of a suitably defined offset parameter and for a range of positions of the leading edge of the splitter plate. The acoustic field is also calculated and the spectrum reflects the flow results.

Description: This paper reports on a magnetic field phenomenon, hereafter referred to as null fields, which were discovered during the inbound pass of the recent flyby of Jupiter by the Ulysses spacecraft. These null fields which were observed in the outer dayside magnetosphere are characterised by brief but sharp decreases of the field magnitude to values less than 1 nT. The nulls are distinguished from the current sheet signatures characteristic of the middle magnetosphere by the fact that the field does not reverse across the event. A field configuration is suggested that accounts for the observed features of the events.

Description: Geochemical profiles of surface units, impact, and volcanic features are studied in detail to determine the underlying structure in an area of extensive mare/highland interface, Sinus Amoris. This study region includes and surrounds the northeastern embayment of Mare Tranquillitatis. The concentrations of two major rock-forming elements (Mg and Al), which were derived from the Apollo 15 orbital geochemical measurements, were used in this study. Mapped units and deposits associated with craters in the northwestern part of the region tend to have correlated low Mg and Al concentrations, indicating the presence of Potassium (K)-Rare Earth Elements (REE)-Phosphorus (P) (KREEP)-enriched basalt. Found along the northeastern rim of Tranquillitatis were areas with correlated high Mg and Al concentration, indicating the presence of troctolite. Distinctive west/east and north/south trends were observed in the concentrations of Mg and Al, and, by implication, in the distribution of major rock components on the surface. Evidence for a systematic geochemical transition in highland or basin-forming units may be observed here in the form of distinctive differences in chemistry in otherwise similar units in the western and eastern portions of the study region.

Description: Numerical results obtained with direct simulation Monte Carlo and Navier-Stokes methods are presented for a Mach-20 nitrogen flow about a 70-deg blunted cone. The flow conditions simuulated are those that can be obtained in existing low-density hypersonic wind tunnels. Three sets of flow conditions are considered with freestream Knudsen numbers ranging from 0.03 to 0.001. The focus is on the wake structure: how the wake structure changes as a function of rarefaction, what the afterbody levels of heating are, and to what limits the continuum models are realistic as rarefaction in the wake is progressively increased. Calculations are made with and without an afterbody sting. Results for the after body sting are emphasizes in anticipation of an experimental study for the current flow conditions and model configuration. The Navier-Stokes calculations were made with and without slip boundary conditions. Comparisons of the results obtained with the two simulation methodologies are made for both flowfield structure and surface quantities.

Description: The aerobraking orbital activities of Magelland during the gravity mapping of Venus are discussed. The goal of aerobraking was to circularize Magellan's orbit. By aerobraking the spacecraft into a nearly circula orbit, the Magellan team was able to provide scientists with a different data set to deepen their understanding of what is going on beneath Venus' surface. Before undertaking its gravity-mapping mission, Magellan completed three cycles of radar mapping. This repeated coverage allowed the spacecraft to see some of Venus' geologic features from different viewing angles. Various aspects of the mission are discussed, and maps of Venus are presented.

Description: We have combined the most recent Pioneer Venus Orbiter (PVO) and Magellan (MGN) data with the earlier 1978-1982 PVO data set to obtain a new 60th degree and order spherical harmonic gravity model and a 120th degree and order spherical harmonic topography model. Free-air gravity maps are shown over regions where the most marked improvement has been obtained (Ishtar-Terra, Alpha, Bell and Artemis). Gravity versus topography relationships are presented as correlations per degree and axes orientation.

Description: We report infrared heterodyne spectroscopy (lambda/delta lambda is approximately 10(exp 6)) of C2H6 emission at 11.9 microns from the northern Jovian auroral region, in observations conducted over December 2-7, 1989. Accurately measured line shapes provide information on C2H6 abundance as well as temperature and permit retrieval of the source pressure region. Enhanced emission was observed in the longitude range approximately 150-180 deg at approximately 60 deg north latitude, approximately corresponding to the CH4 7.8 micron hot spot and the region of brightest UV aurora. Significant brightness variations were observed in the hot spot emissions on a time scale of approximately 20 hours. Analysis of the brightest hot spot spectra indicates C2H6 mole fractions of approximately (6.3-6.8) x 10(exp -6) at temperatures of approximately 182-184 K at 1 mbar, compared to mole fractions of (3.8 +/- 1.4) x 10(exp -6) averaged over spectra outside the hot spot at a temperature of approximately 172 K at the same pressure. Fixing the mole fraction to the lower limit retrieved in the quiescent (non-hot spot) region allows the temperature at 1 mbar to be as high as approximately 200 K within the hot spot. These results provide upper limits to the temperature increase near the source of the C2H6 thermal infrared emission. Combined with results from similar measurements of ethylene emission probing the approximately 10-microbar region (Kostiuk et al., this issue), altitude information on the thermal structure of the Jovian auroral stratosphere can be obtained for the first time.

Description: Doppler tracking data of three orbiting spacecraft have been reanalyzed to develop a new gravitational field model for the planet Mars, Goddard Mars Model 1 (GMM-1). This model employs nearly all available data, consisting of approximately 1100 days of S band tracking data collected by NASA's Deep Space Network from the Mariner 9 and Viking 1 and Viking 2 spacecraft, in seven different orbits, between 1971 and 1979. GMM-1 is complete to spherical harmonic degree and order 50, which corresponds to a half-wavelength spatial resolution of 200-300 km where the data permit. GMM-1 represents satellite orbits with considerably better accuracy than previous Mars gravity models and shows greater resolution of identifiable geological structures. The notable improvement in GMM-1 over previous models is a consequence of several factors: improved computational capabilities, the use of otpimum weighting and least squares collocation solution techniques which stabilized the behavior of the solution at high degree and order, and the use of longer satellite arcs than employed in previous solutions that were made possible by improved force and measurement models. The inclusion of X band tracking data from the 379-km altitude, nnear-polar orbiting Mars Observer spacecraft should provide a significant improvement over GMM-1, particularly at high latitudes where current data poorly resolve the gravitational signature of the planet.

Description: A series of surface-modified clays containing nanophase (np) iron/oxyhydroxides of extremely small particle sizes, with total iron contents as high as found in Mars soil, were prepared by iron deposition on the clay surface from ferrous chloride solution. Comprehensive studies of the iron mineralogy in these 'Mars-soil analogs' were conducted using chemical extractions, solubility analyses, pH and redox, x ray and electron diffractometry, electron microscopic imaging specific surface area and particle size determinations, differential thermal analyses, magnetic properties characterization, spectral reflectance, and Viking biology simulation experiments. The clay matrix and the procedure used for synthesis produced nanophase iron oxides containing a certain proportion of divalent iron, which slowly converts to more stable, fully oxidized iron minerals. The noncrystalline nature of the iron compounds precipitated on the surface of the clay was verified by their complete extractability in oxalate. Lepidocrocite (gamma-FeOOH) was detected by selected area electron diffraction. It is formed from a double iron Fe(II)/Fe(III) hydroxyl mineral such as 'green rust', or ferrosic hydroxide. Magnetic measurements suggested that lepidocrocite converted to the more stable meaghemite (gamma-Fe203) by mild heat treatment and then to nanophase hematite (aplha-Fe203) by extensive heat treatment. Their chemical reactivity offers a plausible mechanism for the somewhat puzzling observations of the Viking biology experiments. Their unique chemical reactivities are attributed to the combined catalytic effects of the iron oxide/oxyhydroxide and silicate phase surfaces. The mode of formation of these (nanophase) iron oxides on Mars is still unknown.

Description: Spatial correlation among densely packed particles can substantially change their single-scattering properties, thus making questionable the applicability of the independent scattering approximation in calculations of light scattering by planetary regoliths. The same problem arises in geophysics in light scattering computations for snow, frosts, and bare soil. In this paper, we use a dense-medium light-scattering theory based on the introduction of the static structure factor to calculate asymmetry parameters of the phase function for densely packed particles with real refractive indices 1.31 and 1.66, approximating water ice and soil particles, respectively, and imaginary refractive indices 0, 0.01, and 0.3. For sparsely distributed, independently scattering grains, the calculated asymmetry parameters are always positive and always larger than those for densely packed particles. For densely packed grains, the asymmetry parameters may be negative but only for radius-to-wavelength ratios from about 0.1 to about 0.4. With decreasing particle size, the calculated asymmetry parameters tend to zero independently of the compaction state. In the geometrical optics regime, the asymmetry parameters for densely packed scatterers are positive and very close to those for independently scattering grains. These results may have important implications for remote sensing of the Earth and solid planetary surfaces. In particular, it is demonstrated that negative asymmetry parameters derived with some approximate multiple-scattering theories may be physically irrelevant and can be the result of using an inaccurate bidirectional reflection function combined with the ill-conditionally of the inverse scattering problem.

Description: The 3-5 micrometer thermal emission of the nightside of Venus, recorded by the Near-Infrared Mapping Spectrometer (NIMS) instrument at the time of the Galileo flyby of Venus, is analysed to infer the properties of the upper cloud boundary. From the global maps of Venus at fixed wavelengths, the limb darkening of the flux is measured at several latitudes, within each infrared channel. By using the nominal Pioneer Venus thermal profile, these data give access to two parameters: the cloud deck temperature and the cloud scale height. It is verified independently, from the NIMS spectra, that this thermal profile is consistent with all the NIMS observations, and that the thermal structure does not vary significantly in the latitude range (25 deg S, 30 deg N). Within this range, the cloud scale height is found to be constant with latitude, and is H = 5.2 km, with an accuracy of about 15%, taking into account the various sources of theoretical and observational uncertainties. At higher latitudes, the temperature profile becomes more isothermal and the presented method to retrieve H is no longer valid.

Description: A large number of i.r. spectra of Venus was obtained using the Near-Infrared Mapping Spectrometer (NIMS) on the Galileo spacecraft, during the February 1990 encounter. Preliminary results show an apparent increase in the tropospheric CO volume mixing ratio (vmr) in the northern polar region. Other possible explanations of the observations are examined and rejected and an increase of the CO abundance north of 47 deg N of (35 +/- 15)% is inferred. Some possible causes of this enhancement are suggested.

Description: The spectroscopic data of the Near-Infrared Mapping Spectrometer (NIMS), recorded during the Galileo flyby of Venus, are analysed to retrieve the water vapor abundance variations in the lower atmosphere of Venus at night. The 1.18 micrometer spectral window, which probes altitude levels below 20 km, is used for this purpose. Constraints on the CO2 continuum and far-wing opacity from existing ground-based high-resolution observations are included in the modelling of the NIMS spectra. The NIMS measurements can be fitted with a water vapor mixing ratio of 30 +/- 15 ppm, in agreement with analyses of ground-based nightside observations. The water vapor abundance shows no horizontal variations exceeding 20% over a wide latitude range (40 deg S, 50 deg N) on the nightside of Venus. Within the same selection of NIMS spectra, a large enhancement in the O2 fluorescence emission at 1.27 micrometer is observed at a latitude of 40 deg S, over a spatial area about 100 km wide.

Description: Using Venus nightside data obtained by the Galileo Near-Infrared Mapping Spectrometer (NIMS), we have studied the correlation of 1.74 and 2.30 micrometer radiation which is transmitted through the clouds. Since the scattering and absorption properties of the cloud particles are different at these two wavelengths, one can distinguish between abundance variations and variations in the properties of the cloud particles themselves. The correlation of intensities shows a clustering of data into five distinct branches. Using radiative transfer calculations, we interpret these branches as regions of distinct but different mixes of Mode 2' and 3 particles. The data and calculations indicate large differences in these modal ratios, the active cloud regions varying in content from nearly pure Mode 2' particles to almost wholly Mode 3. The spatial distribution of these branches shows large scale sizes and both hemispheric symmetries and asymmetries. High-latitude concentrations of large particles are seen in both hemispheres and there is banded structure of small particles seen in both the North and South which may be related. The mean particle size in the Northern Hemisphere is greater than found in the South. If these different branch regions are due to mixing of vertically stratified source regions (e.g. photochemical and condensation source mechanisms), then the mixing must be coherent over very large spatial scales.

Description: The Radio and Plasma Wave Science (RPWS) experiment being built for the Cassini spacecraft will study a wide range of plasma and radio wave phenomena in the magnetosphere of Saturn and will also make valuable measurements during the cruise phase and at other encounters. A feature of data from wave receivers is the capability of producing vastly more data than the spacecraft telemetry link is capable of transmitting back to the Earth. Thus, techniques of on-board data compression and data reduction are important. The RPWS instrument has one processor dedicated to data compression tasks.

Description: The first comprehensive discussion of the south seasonal polar cap spectra obtained by the Mariner 7 infrared spectrometer in the short-wavelength region (2-4 microns) is presented. The infrared spectra is correlated with images acquired by the wide-angle camera. Significant spectral variation is noted in the cap interior and regions of varying water frost abundance, CO2 ice/frost cover, and CO2-ice path length can be distinguished. Many of these spectral variations correlate with heterogeneity noted in the camera images, but certain significant infrared spectral variations are not discernible in the visible. Simple reflectance models are used to classify the observed spectral variations into four regions. Region I is at the cap edge, where there is enhanced absorption beyond 3 microns inferred to be caused by an increased abundance of water frost. The increase in water abundance over that in the interior is on the level of a few parts per thousand or less. Region II is the typical cap interior characterized by spectral features of CO2 ice at grain sizes of several millimeters to centimeters. These spectra also indicate the presence of water frost at the parts per thousand level. A third, unusual region (III), is defined by three spectra in which weak CO2 absorption features are as much as twice as strong as in the average cap spectra and are assumed to be caused by an increased path length in the CO2. Such large paths are inconsistent with the high reflectance in the visible and at 2.2 microns and suggest layered structures or deposition conditions that are not accounted for in current reflectance models. The final region (IV) is an area of thinning frost coverage or transparent ice well in the interior of the seasonal cap. These spectra are a combination of CO2 and ground signatures.

Description: The transformation validity question utilizing resulting data from direct numerical simulations (DNS) of supersonic, isothermal cold wall channel flow was investigated. The DNS results stood for a wide scope of parameter and were suitable for the purpose of examining the generality of Van Driest transformation. The Van Driest law of the wall can be obtained from the inner-layer similarity arguments. It was demonstrated that the Van Driest transformation cannot be incorporated to collapse the sublayer and log-layer velocity profiles simultaneously. Velocity and temperature predictions according to the preceding composite mixing-length model were presented. Despite satisfactory congruity with the DNS data, the model must be perceived as an engineering guide and not as a rigorous analysis.

Description: A Monte Carlo solution technique has been formulated to predict the radiative heat transfer in three-dimensional, inhomogeneous participating media which exhibit spectrally dependent emission and absorption and anisotropic scattering. Details of the technique and selected numerical sensitivities are discussed. The technique was applied to a problem involving a medium composed of a gas mixture of carbon dioxide and nitrogen and suspended carbon particles. A homogeneous medium was modeled to examine the effect of total pressure and carbon-particle concentration on radiative heat transfer. Variation in total pressure, over the range studied, had minimal effect on the amount of heat radiated to the enclosure walls and on the radiative-flux distribution within the medium. Increases in the carbon particle concentration produced significantly higher heat fluxes at the boundaries and altered the radiative flux distribution. The technique was then applied to an inhomogeneous medium to examine effects of specific temperature and carbon particle concentration distributions on radiative heat transfer. For the inhomogeneous conditions examined, the largest radiative flux divergence occurs near the center of the medium and the regions near some enclosure walls act as energy sinks.

Description: We report new measurements of the sodium emission intensity seen in a line of sight just above the surface of the Moon. These data show a strong dependence on lunar phase. The emission intensity decreases from a maximum around first quarter (phase angle 90 deg) to very small values near full Moon (phase angle 0 deg). This suggests that the rate of sodium vapor production from the lunar surface is largest at the subsolar point and becomes small near the terminator. However, the sodium emission near full Moon falls below that which would be expected for solar photon-driven processes. Since the solar wind flux decreases substantially when the Moon enters the Earth's magnetotail near full Moon, while the global solar photon flux is undiminished, we suggest that solar wind sputtering is the dominant process for sodium production.

Description: Numerical results obtained with direct simulation Monte Carlo and Navier-Stokes methods are presented for a Mach-20 nitrogen flow about a 70-deg blunted cone. The flow conditions simulated are those that can be obtained in existing low-density hypersonic wind tunnels. Three sets of flow conditions are considered with freestream Knudsen numbers ranging from 0.03 to 0.001. The focus is on the wake structure: how the wake structure changes as a function of rare faction, what the afterbody levels of heating are, and to what limits the continuum models are realistic as rarefunction in the wake is progressively increased. Calculations are made with and without an afterbody sting. Results for the afterbody sting are emphasized in anticipation of an experimental study for the current flow conditions and model configuration. The Navier-Stokes calculations were made with and without slip boundary conditions. Comparisons of the results obtained with the two simulation methodologies are made for both flowfield structure and surface quantities.

Description: Infrared diffuse reflectance spectra (2.53-25 microns) of some carbonaceous (C) chondrites were measured. The integrated intensity of the absorption bands near 3 microns caused by hydrous minerals were compared with the modal content of hydrous minerals for the meteorites. The CM and CI chondrites show larger values of the intergated intensity than those of the unique C chondrites Y82162, Y86720 and B7904, suggesting that the amount of hydrous minerals in the CM and CI chondrites is larger, which supports the contention that hydrous minerals were dehydrated by thermal metamorphism in the unique chondrites. Orgueil (CI) has the largest value of the integrated intensity among the C chondrites we measured and shows a sharp absorption band at 3685/cm (2.71 microns) that is not seen in the spectra of the CM chondrites. There is an excellent correlation between the observed hydrogen content in C chondrites and the integrated intensity. The CM chondrites show a wide variation in the strength of absorption bands at 1470/cm (6.8 microns), despite the similarity in absorption features near 3 micron for all CM chondites. The 1470/cm band could be due to the presence of some hydrocarbons but may also be a result of terrestrial alteration processes.

Description: Oxygen production from a lunar rock has been experimentally demonstrated for the first time. A 10 g sample of high-Ti basalt 70035 was reduced with hydrogen in seven experiments at temperatures of 900-1050 C and pressures of 14.7-150 psia. In all experiments, water evolution began almost immediately and was essentially complete in tens of minutes. Oxygen yields ranged from 2.93 to 4.61% of the starting sample weight, and showed weak dependence on temperature and pressure. Analysis of the solid samples demonstrated total reduction of Fe(2+) in ilmenite and small degrees of reduction in olivine and pyroxene. Ti O2 was also partially reduced to one or more suboxides. Data from these experiments provide a basis for predicting the yield of oxygen from lunar basalt as well as new constraints on natural reduction in the lunar regolith.

Description: The Simplified Shuttle Payload Thermal Analyzer program (SSPTA) was developed to aid in the evaluation of thermal design concepts of instruments to be flown in the Space Shuttle cargo bay. SSPTA consists of a collection of programs that are currently used in the thermal analysis of spacecraft and have been modified for quick, preliminary analysis of payloads. SSPTA includes a reduced math model of the Shuttle cargo bay to simplify use of the program for payload analysis. One of the prime objectives in developing SSPTA was to create a program which was easy to use. With SSPTA, the user required input is simple and the user is free from many of the concerns of computer usage such as disk space handling, tape usage, and complicated program control. Although SSPTA was designed primarily to analyze Shuttle payloads, it can easily be used to perform thermal analysis in other situations. SSPTA is comprised of a system of data files called 'bins', a master program, and a set of thermal subprograms. The bin system is a collection of disk files which contain data required by or computed by the thermal subprograms. SSPTA currently has the capability of handling 50 bins. The master program serves primarily as a manager for the bin system and its interaction with the thermal subprograms. Input to the master program consists of simple user commands which direct the data manipulation procedures, prepare the data for these procedures, and call the appropriate thermal subprograms. The subprograms of SSPTA are all based on programs which have been used extensively in the analysis of orbiting spacecraft and space hardware. Subprogram CONSHAD uses the user supplied geometric radiation model to compute black body view factors, shadow factors, and a description of the surface model. The subprogram WORKSHEET uses the surface model description, optical property data, and node assignment data to prepare input for SCRIPTF. Subprogram SCRIPTF computes the inverses of the infrared (IR) and ultraviolet (UV) radiation transfer equations; it also computes the radiation coupling between nodes in the thermal model. Subprogram ORBITAL uses the shadow tables to compute incident flux intensities on each surface in the geometric model. Subprogram ABSORB uses these flux intensities combined with the IR and UV inverses to compute the IR and UV fluxes absorbed by each surface. The radiation couplings from SCRIPTF and the absorbed fluxes from ABSORB are used by subprogram TTA to compute the temperature and power balance for each node in the thermal model. Output consists of tabulated data from each of the subprograms executed during a particular analysis. Due to the modular form of SSPTA, analyses may be run in whole or in part, and new subprograms may be added by the user. SSPTA is written in FORTRAN for use on a DEC VAX-11/780. SSPTA was originally developed in 1977 for use on IBM 370 series computers. This version is an update which was ported to the VAX in 1980.

Description: Remote Interactive Particle-tracing (RIP) is a distributed-graphics program which computes particle traces for computational fluid dynamics (CFD) solution data sets. A particle trace is a line which shows the path a massless particle in a fluid will take; it is a visual image of where the fluid is going. The program is able to compute and display particle traces at a speed of about one trace per second because it runs on two machines concurrently. The data used by the program is contained in two files. The solution file contains data on density, momentum and energy quantities of a flow field at discrete points in three-dimensional space, while the grid file contains the physical coordinates of each of the discrete points. RIP requires two computers. A local graphics workstation interfaces with the user for program control and graphics manipulation, and a remote machine interfaces with the solution data set and performs time-intensive computations. The program utilizes two machines in a distributed mode for two reasons. First, the data to be used by the program is usually generated on the supercomputer. RIP avoids having to convert and transfer the data, eliminating any memory limitations of the local machine. Second, as computing the particle traces can be computationally expensive, RIP utilizes the power of the supercomputer for this task. Although the remote site code was developed on a CRAY, it is possible to port this to any supercomputer class machine with a UNIX-like operating system. Integration of a velocity field from a starting physical location produces the particle trace. The remote machine computes the particle traces using the particle-tracing subroutines from PLOT3D/AMES, a CFD post-processing graphics program available from COSMIC (ARC-12779). These routines use a second-order predictor-corrector method to integrate the velocity field. Then the remote program sends graphics tokens to the local machine via a remote-graphics library. The local machine interprets the graphics tokens and draws the particle traces. The program is menu driven. RIP is implemented on the silicon graphics IRIS 3000 (local workstation) with an IRIX operating system and on the CRAY2 (remote station) with a UNICOS 1.0 or 2.0 operating system. The IRIS 4D can be used in place of the IRIS 3000. The program is written in C (67%) and FORTRAN 77 (43%) and has an IRIS memory requirement of 4 MB. The remote and local stations must use the same user ID. PLOT3D/AMES unformatted data sets are required for the remote machine. The program was developed in 1988.

Description: This software package includes two programs, the KPD12 and the KPD12P. Both programs utilizes the vortex-blob method to simulate flow around solid bodies, in an unbounded domain using the KPD12, with periodicity in one direction using the KPD12P. The main advantage of the vortex-blob method is the ability to handle situations involving arbitrary shapes including multiple bodies. The user just supplies points on the solid boundaries; there is no grid. The KPD12 program has worked successfully on bluff bodies, stalled wings, and multiple-element airfoils. The KPD12P program has been used successfully on high-solidity separated cascades and on cases of rotating stall in cascades of thin airfoils. However, they do not capture subtle viscous effects such as incipient separation and friction drag. The KPD12 and the KPD12P programs apply the vortex-blob method to time-dependent, high-Reynolds-number flows around solid bodies. Both programs solve the two-dimensional incompressible Navier-Stokes equations, neglecting the viscous effects away from the walls. By creating new vortices along the wall at every time step, they treat the no-penetration and no-slip boundary conditions while using an influence matrix. The code automatically controls the number of vortices. Furthermore, the code has the option of treating the boundary layers by simple integral methods to determine the separation points. The KPD12 outputs forces, moments, and pressure distributions on the bodies. The KPD12P also outputs the turning angle and loss of total pressure. The source code is in Cray FORTRAN and contains a few calls to Cray vector functions which are vectorized with the Cray compiler. However, substitutes for these vector functions are provided. The code is set up to plot the bodies, vortex positions, and streamlines using the DISSPLA graphics software. The software requires a mainframe computer with at least 589k of memory available running under COS 1.16. KPD12 was developed in 1988.

Description: The Steady State Thermal Analysis Program (STEADY) provides the thermal designer with a quick and convenient method for calculating heat loads and temperatures. STEADY can be used on small nodal networks for conceptual or preliminary thermal design and analysis. STEADY will accept up to 20 nodes of fixed or variable temperature, with constant or temperature-dependent thermal conductivities, and any set of consistent units. In a steady state thermal network, the heat balance on each variable temperature node must sum to zero. The general heat transfer equations are solved with a Newton-Raphson technique and refined by a fourth order quartic solution. Input data includes the number of nodes, number of boundary nodes, the fixed temperatures at all boundary nodes, initial temperature guesses for variable nodes, impressed heat loads, conduction and radiation coefficients, and control parameters such as convergence criteria, maximum iterations, and damping factors. The output is stored in a print file and tabulates final temperatures and heat flows for all nodes. STEADY is menu driven and allows the user to save files for future modification. STEADY is written in FORTRAN 77 (Ryan McFarland's RMFORTRAN) for interactive execution and has been implemented on the IBM PC computer series under DOS with a central memory requirement of approximately 92K of 8 bit bytes using a math coprocessor, and 103K bytes without the coprocessor. This program was developed in 1987.

Description: The ORBSIM program was developed for the accurate extraction of geophysical model parameters from Doppler radio tracking data acquired from orbiting planetary spacecraft. The model of the proposed planetary structure is used in a numerical integration of the spacecraft along simulated trajectories around the primary body. Using line of sight (LOS) Doppler residuals, ORBSIM applies fast and efficient modelling and optimization procedures which avoid the traditional complex dynamic reduction of data. ORBSIM produces quantitative geophysical results such as size, depth, and mass. ORBSIM has been used extensively to investigate topographic features on the Moon, Mars, and Venus. The program has proven particulary suitable for modelling gravitational anomalies and mascons. The basic observable for spacecraft-based gravity data is the Doppler frequency shift of a transponded radio signal. The time derivative of this signal carries information regarding the gravity field acting on the spacecraft in the LOS direction (the LOS direction being the path between the spacecraft and the receiving station, either Earth or another satellite). There are many dynamic factors taken into account: earth rotation, solar radiation, acceleration from planetary bodies, tracking station time and location adjustments, etc. The actual trajectories of the spacecraft are simulated using least squares fitted to conic motion. The theoretical Doppler readings from the simulated orbits are compared to actual Doppler observations and another least squares adjustment is made. ORBSIM has three modes of operation: trajectory simulation, optimization, and gravity modelling. In all cases, an initial gravity model of curved and/or flat disks, harmonics, and/or a force table are required input. ORBSIM is written in FORTRAN 77 for batch execution and has been implemented on a DEC VAX 11/780 computer operating under VMS. This program was released in 1985.

Description: The Thermal Radiation Analyzer System, TRASYS, is a computer software system with generalized capability to solve the radiation related aspects of thermal analysis problems. TRASYS computes the total thermal radiation environment for a spacecraft in orbit. The software calculates internode radiation interchange data as well as incident and absorbed heat rate data originating from environmental radiant heat sources. TRASYS provides data of both types in a format directly usable by such thermal analyzer programs as SINDA/FLUINT (available from COSMIC, program number MSC-21528). One primary feature of TRASYS is that it allows users to write their own driver programs to organize and direct the preprocessor and processor library routines in solving specific thermal radiation problems. The preprocessor first reads and converts the user's geometry input data into the form used by the processor library routines. Then, the preprocessor accepts the user's driving logic, written in the TRASYS modified FORTRAN language. In many cases, the user has a choice of routines to solve a given problem. Users may also provide their own routines where desirable. In particular, the user may write output routines to provide for an interface between TRASYS and any thermal analyzer program using the R-C network concept. Input to the TRASYS program consists of Options and Edit data, Model data, and Logic Flow and Operations data. Options and Edit data provide for basic program control and user edit capability. The Model data describe the problem in terms of geometry and other properties. This information includes surface geometry data, documentation data, nodal data, block coordinate system data, form factor data, and flux data. Logic Flow and Operations data house the user's driver logic, including the sequence of subroutine calls and the subroutine library. Output from TRASYS consists of two basic types of data: internode radiation interchange data, and incident and absorbed heat rate data. The flexible structure of TRASYS allows considerable freedom in the definition and choice of solution method for a thermal radiation problem. The program's flexible structure has also allowed TRASYS to retain the same basic input structure as the authors update it in order to keep up with changing requirements. Among its other important features are the following: 1) up to 3200 node problem size capability with shadowing by intervening opaque or semi-transparent surfaces; 2) choice of diffuse, specular, or diffuse/specular radiant interchange solutions; 3) a restart capability that minimizes recomputing; 4) macroinstructions that automatically provide the executive logic for orbit generation that optimizes the use of previously completed computations; 5) a time variable geometry package that provides automatic pointing of the various parts of an articulated spacecraft and an automatic look-back feature that eliminates redundant form factor calculations; 6) capability to specify submodel names to identify sets of surfaces or components as an entity; and 7) subroutines to perform functions which save and recall the internodal and/or space form factors in subsequent steps for nodes with fixed geometry during a variable geometry run. There are two machine versions of TRASYS v27: a DEC VAX version and a Cray UNICOS version. Both versions require installation of the NASADIG library (MSC-21801 for DEC VAX or COS-10049 for CRAY), which is available from COSMIC either separately or bundled with TRASYS. The NASADIG (NASA Device Independent Graphics Library) plot package provides a pictorial representation of input geometry, orbital/orientation parameters, and heating rate output as a function of time. NASADIG supports Tektronix terminals. The CRAY version of TRASYS v27 is written in FORTRAN 77 for batch or interactive execution and has been implemented on CRAY X-MP and CRAY Y-MP series computers running UNICOS. The standard distribution medium for MSC-21959 (CRAY version without NASADIG) is a 1600 BPI 9-track magnetic tape in UNIX tar format. The standard distribution medium for COS-10040 (CRAY version with NASADIG) is a set of two 6250 BPI 9-track magnetic tapes in UNIX tar format. Alternate distribution media and formats are available upon request. The DEC VAX version of TRASYS v27 is written in FORTRAN 77 for batch execution (only the plotting driver program is interactive) and has been implemented on a DEC VAX 8650 computer under VMS. Since the source codes for MSC-21030 and COS-10026 are in VAX/VMS text library files and DEC Command Language files, COSMIC will only provide these programs in the following formats: MSC-21030, TRASYS (DEC VAX version without NASADIG) is available on a 1600 BPI 9-track magnetic tape in VAX BACKUP format (standard distribution medium) or in VAX BACKUP format on a TK50 tape cartridge; COS-10026, TRASYS (DEC VAX version with NASADIG), is available in VAX BACKUP format on a set of three 6250 BPI 9-track magnetic tapes (standard distribution medium) or a set of three TK50 tape cartridges in VAX BACKUP format. TRASYS was last updated in 1993.

Description: TDIGG is a fast and versatile program for generating two-dimensional computational grids for use with finite-difference flow-solvers. Both algebraic and elliptic grid generation systems are included. The method for grid generation by algebraic transformation is based on an interpolation algorithm and the elliptic grid generation is established by solving the partial differential equation (PDE). Non-uniform grid distributions are carried out using a hyperbolic tangent stretching function. For algebraic grid systems, interpolations in one direction (univariate) and two directions (bivariate) are considered. These interpolations are associated with linear or cubic Lagrangian/Hermite/Bezier polynomial functions. The algebraic grids can subsequently be smoothed using an elliptic solver. For elliptic grid systems, the PDE can be in the form of Laplace (zero forcing function) or Poisson. The forcing functions in the Poisson equation come from the boundary or the entire domain of the initial algebraic grids. A graphics interface procedure using the Silicon Graphics (GL) Library is included to allow users to visualize the grid variations at each iteration. This will allow users to interactively modify the grid to match their applications. TDIGG is written in FORTRAN 77 for Silicon Graphics IRIS series computers running IRIX. This package requires either MIT's X Window System, Version 11 Revision 4 or SGI (Motif) Window System. A sample executable is provided on the distribution medium. It requires 148K of RAM for execution. The standard distribution medium is a .25 inch streaming magnetic IRIX tape cartridge in UNIX tar format. This program was developed in 1992.

Description: The Systems Improved Numerical Fluids Analysis Code, SINFAC, consists of additional routines added to the April 1983 revision of SINDA, a general thermal analyzer program. The purpose of the additional routines is to allow for the modeling of active heat transfer loops. The modeler can simulate the steady-state and pseudo-transient operations of 16 different heat transfer loop components including radiators, evaporators, condensers, mechanical pumps, reservoirs and many types of valves and fittings. In addition, the program contains a property analysis routine that can be used to compute the thermodynamic properties of 20 different refrigerants. SINFAC can simulate the response to transient boundary conditions. SINFAC was first developed as a method for computing the steady-state performance of two phase systems. It was then modified using CNFRWD, SINDA's explicit time-integration scheme, to accommodate transient thermal models. However, SINFAC cannot simulate pressure drops due to time-dependent fluid acceleration, transient boil-out, or transient fill-up, except in the accumulator. SINFAC also requires the user to be familiar with SINDA. The solution procedure used by SINFAC is similar to that which an engineer would use to solve a system manually. The solution to a system requires the determination of all of the outlet conditions of each component such as the flow rate, pressure, and enthalpy. To obtain these values, the user first estimates the inlet conditions to the first component of the system, then computes the outlet conditions from the data supplied by the manufacturer of the first component. The user then estimates the temperature at the outlet of the third component and computes the corresponding flow resistance of the second component. With the flow resistance of the second component, the user computes the conditions down stream, namely the inlet conditions of the third. The computations follow for the rest of the system, back to the first component. On the first pass, the user finds that the calculated outlet conditions of the last component do not match the estimated inlet conditions of the first. The user then modifies the estimated inlet conditions of the first component in an attempt to match the calculated values. The user estimated values are called State Variables. The differences between the user estimated values and calculated values are called the Error Variables. The procedure systematically changes the State Variables until all of the Error Variables are less than the user-specified iteration limits. The solution procedure is referred to as SCX. It consists of two phases, the Systems phase and the Controller phase. The X is to imply experimental. SCX computes each next set of State Variables in two phases. In the first phase, SCX fixes the controller positions and modifies the other State Variables by the Newton-Raphson method. This first phase is the Systems phase. Once the Newton-Raphson method has solved the problem for the fixed controller positions, SCX next calculates new controller positions based on Newton's method while treating each sensor-controller pair independently but allowing all to change in one iteration. This phase is the Controller phase. SINFAC is available by license for a period of ten (10) years to approved licensees. The licenced program product includes the source code for the additional routines to SINDA, the SINDA object code, command procedures, sample data and supporting documentation. Additional documentation may be purchased at the price below. SINFAC was created for use on a DEC VAX under VMS. Source code is written in FORTRAN 77, requires 180k of memory, and should be fully transportable. The program was developed in 1988.

Description: The Thermal Radiation Analyzer System, TRASYS, is a computer software system with generalized capability to solve the radiation related aspects of thermal analysis problems. TRASYS computes the total thermal radiation environment for a spacecraft in orbit. The software calculates internode radiation interchange data as well as incident and absorbed heat rate data originating from environmental radiant heat sources. TRASYS provides data of both types in a format directly usable by such thermal analyzer programs as SINDA/FLUINT (available from COSMIC, program number MSC-21528). One primary feature of TRASYS is that it allows users to write their own driver programs to organize and direct the preprocessor and processor library routines in solving specific thermal radiation problems. The preprocessor first reads and converts the user's geometry input data into the form used by the processor library routines. Then, the preprocessor accepts the user's driving logic, written in the TRASYS modified FORTRAN language. In many cases, the user has a choice of routines to solve a given problem. Users may also provide their own routines where desirable. In particular, the user may write output routines to provide for an interface between TRASYS and any thermal analyzer program using the R-C network concept. Input to the TRASYS program consists of Options and Edit data, Model data, and Logic Flow and Operations data. Options and Edit data provide for basic program control and user edit capability. The Model data describe the problem in terms of geometry and other properties. This information includes surface geometry data, documentation data, nodal data, block coordinate system data, form factor data, and flux data. Logic Flow and Operations data house the user's driver logic, including the sequence of subroutine calls and the subroutine library. Output from TRASYS consists of two basic types of data: internode radiation interchange data, and incident and absorbed heat rate data. The flexible structure of TRASYS allows considerable freedom in the definition and choice of solution method for a thermal radiation problem. The program's flexible structure has also allowed TRASYS to retain the same basic input structure as the authors update it in order to keep up with changing requirements. Among its other important features are the following: 1) up to 3200 node problem size capability with shadowing by intervening opaque or semi-transparent surfaces; 2) choice of diffuse, specular, or diffuse/specular radiant interchange solutions; 3) a restart capability that minimizes recomputing; 4) macroinstructions that automatically provide the executive logic for orbit generation that optimizes the use of previously completed computations; 5) a time variable geometry package that provides automatic pointing of the various parts of an articulated spacecraft and an automatic look-back feature that eliminates redundant form factor calculations; 6) capability to specify submodel names to identify sets of surfaces or components as an entity; and 7) subroutines to perform functions which save and recall the internodal and/or space form factors in subsequent steps for nodes with fixed geometry during a variable geometry run. There are two machine versions of TRASYS v27: a DEC VAX version and a Cray UNICOS version. Both versions require installation of the NASADIG library (MSC-21801 for DEC VAX or COS-10049 for CRAY), which is available from COSMIC either separately or bundled with TRASYS. The NASADIG (NASA Device Independent Graphics Library) plot package provides a pictorial representation of input geometry, orbital/orientation parameters, and heating rate output as a function of time. NASADIG supports Tektronix terminals. The CRAY version of TRASYS v27 is written in FORTRAN 77 for batch or interactive execution and has been implemented on CRAY X-MP and CRAY Y-MP series computers running UNICOS. The standard distribution medium for MSC-21959 (CRAY version without NASADIG) is a 1600 BPI 9-track magnetic tape in UNIX tar format. The standard distribution medium for COS-10040 (CRAY version with NASADIG) is a set of two 6250 BPI 9-track magnetic tapes in UNIX tar format. Alternate distribution media and formats are available upon request. The DEC VAX version of TRASYS v27 is written in FORTRAN 77 for batch execution (only the plotting driver program is interactive) and has been implemented on a DEC VAX 8650 computer under VMS. Since the source codes for MSC-21030 and COS-10026 are in VAX/VMS text library files and DEC Command Language files, COSMIC will only provide these programs in the following formats: MSC-21030, TRASYS (DEC VAX version without NASADIG) is available on a 1600 BPI 9-track magnetic tape in VAX BACKUP format (standard distribution medium) or in VAX BACKUP format on a TK50 tape cartridge; COS-10026, TRASYS (DEC VAX version with NASADIG), is available in VAX BACKUP format on a set of three 6250 BPI 9-track magnetic tapes (standard distribution medium) or a set of three TK50 tape cartridges in VAX BACKUP format. TRASYS was last updated in 1993.

Description: INS3D computes steady-state solutions to the incompressible Navier-Stokes equations. The INS3D approach utilizes pseudo-compressibility combined with an approximate factorization scheme. This computational fluid dynamics (CFD) code has been verified on problems such as flow through a channel, flow over a backwardfacing step and flow over a circular cylinder. Three dimensional cases include flow over an ogive cylinder, flow through a rectangular duct, wind tunnel inlet flow, cylinder-wall juncture flow and flow through multiple posts mounted between two plates. INS3D uses a pseudo-compressibility approach in which a time derivative of pressure is added to the continuity equation, which together with the momentum equations form a set of four equations with pressure and velocity as the dependent variables. The equations' coordinates are transformed for general three dimensional applications. The equations are advanced in time by the implicit, non-iterative, approximately-factored, finite-difference scheme of Beam and Warming. The numerical stability of the scheme depends on the use of higher-order smoothing terms to damp out higher-frequency oscillations caused by second-order central differencing. The artificial compressibility introduces pressure (sound) waves of finite speed (whereas the speed of sound would be infinite in an incompressible fluid). As the solution converges, these pressure waves die out, causing the derivation of pressure with respect to time to approach zero. Thus, continuity is satisfied for the incompressible fluid in the steady state. Computational efficiency is achieved using a diagonal algorithm. A block tri-diagonal option is also available. When a steady-state solution is reached, the modified continuity equation will satisfy the divergence-free velocity field condition. INS3D is capable of handling several different types of boundaries encountered in numerical simulations, including solid-surface, inflow and outflow, and far-field boundaries. Three machine versions of INS3D are available. INS3D for the CRAY is written in CRAY FORTRAN for execution on a CRAY X-MP under COS, INS3D for the IBM is written in FORTRAN 77 for execution on an IBM 3090 under the VM or MVS operating system, and INS3D for DEC RISC-based systems is written in RISC FORTRAN for execution on a DEC workstation running RISC ULTRIX 3.1 or later. The CRAY version has a central memory requirement of 730279 words. The central memory requirement for the IBM is 150Mb. The memory requirement for the DEC RISC ULTRIX version is 3Mb of main memory. INS3D was developed in 1987. The port to the IBM was done in 1990. The port to the DECstation 3100 was done in 1991. CRAY is a registered trademark of Cray Research Inc. IBM is a registered trademark of International Business Machines. DEC, DECstation, and ULTRIX are trademarks of the Digital Equipment Corporation.

Description: SAGE, Self Adaptive Grid codE, is a flexible tool for adapting and restructuring both 2D and 3D grids. Solution-adaptive grid methods are useful tools for efficient and accurate flow predictions. In supersonic and hypersonic flows, strong gradient regions such as shocks, contact discontinuities, shear layers, etc., require careful distribution of grid points to minimize grid error and produce accurate flow-field predictions. SAGE helps the user obtain more accurate solutions by intelligently redistributing (i.e. adapting) the original grid points based on an initial or interim flow-field solution. The user then computes a new solution using the adapted grid as input to the flow solver. The adaptive-grid methodology poses the problem in an algebraic, unidirectional manner for multi-dimensional adaptations. The procedure is analogous to applying tension and torsion spring forces proportional to the local flow gradient at every grid point and finding the equilibrium position of the resulting system of grid points. The multi-dimensional problem of grid adaption is split into a series of one-dimensional problems along the computational coordinate lines. The reduced one dimensional problem then requires a tridiagonal solver to find the location of grid points along a coordinate line. Multi-directional adaption is achieved by the sequential application of the method in each coordinate direction. The tension forces direct the redistribution of points to the strong gradient region. To maintain smoothness and a measure of orthogonality of grid lines, torsional forces are introduced that relate information between the family of lines adjacent to one another. The smoothness and orthogonality constraints are direction-dependent, since they relate only the coordinate lines that are being adapted to the neighboring lines that have already been adapted. Therefore the solutions are non-unique and depend on the order and direction of adaption. Non-uniqueness of the adapted grid is acceptable since it makes possible an overall and local error reduction through grid redistribution. SAGE includes the ability to modify the adaption techniques in boundary regions, which substantially improves the flexibility of the adaptive scheme. The vectorial approach used in the analysis also provides flexibility. The user has complete choice of adaption direction and order of sequential adaptions without concern for the computational data structure. Multiple passes are available with no restraint on stepping directions; for each adaptive pass the user can choose a completely new set of adaptive parameters. This facility, combined with the capability of edge boundary control, enables the code to individually adapt multi-dimensional multiple grids. Zonal grids can be adapted while maintaining continuity along the common boundaries. For patched grids, the multiple-pass capability enables complete adaption. SAGE is written in FORTRAN 77 and is intended to be machine independent; however, it requires a FORTRAN compiler which supports NAMELIST input. It has been successfully implemented on Sun series computers, SGI IRIS's, DEC MicroVAX computers, HP series computers, the Cray YMP, and IBM PC compatibles. Source code is provided, but no sample input and output files are provided. The code reads three datafiles: one that contains the initial grid coordinates (x,y,z), one that contains corresponding flow-field variables, and one that contains the user control parameters. It is assumed that the first two datasets are formatted as defined in the plotting software package PLOT3D. Several machine versions of PLOT3D are available from COSMIC. The amount of main memory is dependent on the size of the matrix. The standard distribution medium for SAGE is a 5.25 inch 360K MS-DOS format diskette. It is also available on a .25 inch streaming magnetic tape cartridge in UNIX tar format or on a 9-track 1600 BPI ASCII CARD IMAGE format magnetic tape. SAGE was developed in 1989, first released as a 2D version in 1991 and updated to 3D in 1993.

Description: INS3D computes steady-state solutions to the incompressible Navier-Stokes equations. The INS3D approach utilizes pseudo-compressibility combined with an approximate factorization scheme. This computational fluid dynamics (CFD) code has been verified on problems such as flow through a channel, flow over a backwardfacing step and flow over a circular cylinder. Three dimensional cases include flow over an ogive cylinder, flow through a rectangular duct, wind tunnel inlet flow, cylinder-wall juncture flow and flow through multiple posts mounted between two plates. INS3D uses a pseudo-compressibility approach in which a time derivative of pressure is added to the continuity equation, which together with the momentum equations form a set of four equations with pressure and velocity as the dependent variables. The equations' coordinates are transformed for general three dimensional applications. The equations are advanced in time by the implicit, non-iterative, approximately-factored, finite-difference scheme of Beam and Warming. The numerical stability of the scheme depends on the use of higher-order smoothing terms to damp out higher-frequency oscillations caused by second-order central differencing. The artificial compressibility introduces pressure (sound) waves of finite speed (whereas the speed of sound would be infinite in an incompressible fluid). As the solution converges, these pressure waves die out, causing the derivation of pressure with respect to time to approach zero. Thus, continuity is satisfied for the incompressible fluid in the steady state. Computational efficiency is achieved using a diagonal algorithm. A block tri-diagonal option is also available. When a steady-state solution is reached, the modified continuity equation will satisfy the divergence-free velocity field condition. INS3D is capable of handling several different types of boundaries encountered in numerical simulations, including solid-surface, inflow and outflow, and far-field boundaries. Three machine versions of INS3D are available. INS3D for the CRAY is written in CRAY FORTRAN for execution on a CRAY X-MP under COS, INS3D for the IBM is written in FORTRAN 77 for execution on an IBM 3090 under the VM or MVS operating system, and INS3D for DEC RISC-based systems is written in RISC FORTRAN for execution on a DEC workstation running RISC ULTRIX 3.1 or later. The CRAY version has a central memory requirement of 730279 words. The central memory requirement for the IBM is 150Mb. The memory requirement for the DEC RISC ULTRIX version is 3Mb of main memory. INS3D was developed in 1987. The port to the IBM was done in 1990. The port to the DECstation 3100 was done in 1991. CRAY is a registered trademark of Cray Research Inc. IBM is a registered trademark of International Business Machines. DEC, DECstation, and ULTRIX are trademarks of the Digital Equipment Corporation.

Description: The ability to treat arbitrary boundary shapes is one of the most desirable characteristics of a method for generating grids, including those about airfoils. In a grid used for computing aerodynamic flow over an airfoil, or any other body shape, the surface of the body is usually treated as an inner boundary and often cannot be easily represented as an analytic function. The GRAPE computer program was developed to incorporate a method for generating two-dimensional finite-difference grids about airfoils and other shapes by the use of the Poisson differential equation. GRAPE can be used with any boundary shape, even one specified by tabulated points and including a limited number of sharp corners. The GRAPE program has been developed to be numerically stable and computationally fast. GRAPE can provide the aerodynamic analyst with an efficient and consistent means of grid generation. The GRAPE procedure generates a grid between an inner and an outer boundary by utilizing an iterative procedure to solve the Poisson differential equation subject to geometrical restraints. In this method, the inhomogeneous terms of the equation are automatically chosen such that two important effects are imposed on the grid. The first effect is control of the spacing between mesh points along mesh lines intersecting the boundaries. The second effect is control of the angles with which mesh lines intersect the boundaries. Along with the iterative solution to Poisson's equation, a technique of coarse-fine sequencing is employed to accelerate numerical convergence. GRAPE program control cards and input data are entered via the NAMELIST feature. Each variable has a default value such that user supplied data is kept to a minimum. Basic input data consists of the boundary specification, mesh point spacings on the boundaries, and mesh line angles at the boundaries. Output consists of a dataset containing the grid data and, if requested, a plot of the generated mesh. The GRAPE program is written in FORTRAN IV for batch execution and has been implemented on a CDC 6000 series computer with a central memory requirement of approximately 135K (octal) of 60 bit words. For plotted output the commercially available DISSPLA graphics software package is required. The GRAPE program was developed in 1980.

Description: The Thermal Radiation Analyzer System, TRASYS, is a computer software system with generalized capability to solve the radiation related aspects of thermal analysis problems. TRASYS computes the total thermal radiation environment for a spacecraft in orbit. The software calculates internode radiation interchange data as well as incident and absorbed heat rate data originating from environmental radiant heat sources. TRASYS provides data of both types in a format directly usable by such thermal analyzer programs as SINDA/FLUINT (available from COSMIC, program number MSC-21528). One primary feature of TRASYS is that it allows users to write their own driver programs to organize and direct the preprocessor and processor library routines in solving specific thermal radiation problems. The preprocessor first reads and converts the user's geometry input data into the form used by the processor library routines. Then, the preprocessor accepts the user's driving logic, written in the TRASYS modified FORTRAN language. In many cases, the user has a choice of routines to solve a given problem. Users may also provide their own routines where desirable. In particular, the user may write output routines to provide for an interface between TRASYS and any thermal analyzer program using the R-C network concept. Input to the TRASYS program consists of Options and Edit data, Model data, and Logic Flow and Operations data. Options and Edit data provide for basic program control and user edit capability. The Model data describe the problem in terms of geometry and other properties. This information includes surface geometry data, documentation data, nodal data, block coordinate system data, form factor data, and flux data. Logic Flow and Operations data house the user's driver logic, including the sequence of subroutine calls and the subroutine library. Output from TRASYS consists of two basic types of data: internode radiation interchange data, and incident and absorbed heat rate data. The flexible structure of TRASYS allows considerable freedom in the definition and choice of solution method for a thermal radiation problem. The program's flexible structure has also allowed TRASYS to retain the same basic input structure as the authors update it in order to keep up with changing requirements. Among its other important features are the following: 1) up to 3200 node problem size capability with shadowing by intervening opaque or semi-transparent surfaces; 2) choice of diffuse, specular, or diffuse/specular radiant interchange solutions; 3) a restart capability that minimizes recomputing; 4) macroinstructions that automatically provide the executive logic for orbit generation that optimizes the use of previously completed computations; 5) a time variable geometry package that provides automatic pointing of the various parts of an articulated spacecraft and an automatic look-back feature that eliminates redundant form factor calculations; 6) capability to specify submodel names to identify sets of surfaces or components as an entity; and 7) subroutines to perform functions which save and recall the internodal and/or space form factors in subsequent steps for nodes with fixed geometry during a variable geometry run. There are two machine versions of TRASYS v27: a DEC VAX version and a Cray UNICOS version. Both versions require installation of the NASADIG library (MSC-21801 for DEC VAX or COS-10049 for CRAY), which is available from COSMIC either separately or bundled with TRASYS. The NASADIG (NASA Device Independent Graphics Library) plot package provides a pictorial representation of input geometry, orbital/orientation parameters, and heating rate output as a function of time. NASADIG supports Tektronix terminals. The CRAY version of TRASYS v27 is written in FORTRAN 77 for batch or interactive execution and has been implemented on CRAY X-MP and CRAY Y-MP series computers running UNICOS. The standard distribution medium for MSC-21959 (CRAY version without NASADIG) is a 1600 BPI 9-track magnetic tape in UNIX tar format. The standard distribution medium for COS-10040 (CRAY version with NASADIG) is a set of two 6250 BPI 9-track magnetic tapes in UNIX tar format. Alternate distribution media and formats are available upon request. The DEC VAX version of TRASYS v27 is written in FORTRAN 77 for batch execution (only the plotting driver program is interactive) and has been implemented on a DEC VAX 8650 computer under VMS. Since the source codes for MSC-21030 and COS-10026 are in VAX/VMS text library files and DEC Command Language files, COSMIC will only provide these programs in the following formats: MSC-21030, TRASYS (DEC VAX version without NASADIG) is available on a 1600 BPI 9-track magnetic tape in VAX BACKUP format (standard distribution medium) or in VAX BACKUP format on a TK50 tape cartridge; COS-10026, TRASYS (DEC VAX version with NASADIG), is available in VAX BACKUP format on a set of three 6250 BPI 9-track magnetic tapes (standard distribution medium) or a set of three TK50 tape cartridges in VAX BACKUP format. TRASYS was last updated in 1993.

Description: INS3D computes steady-state solutions to the incompressible Navier-Stokes equations. The INS3D approach utilizes pseudo-compressibility combined with an approximate factorization scheme. This computational fluid dynamics (CFD) code has been verified on problems such as flow through a channel, flow over a backwardfacing step and flow over a circular cylinder. Three dimensional cases include flow over an ogive cylinder, flow through a rectangular duct, wind tunnel inlet flow, cylinder-wall juncture flow and flow through multiple posts mounted between two plates. INS3D uses a pseudo-compressibility approach in which a time derivative of pressure is added to the continuity equation, which together with the momentum equations form a set of four equations with pressure and velocity as the dependent variables. The equations' coordinates are transformed for general three dimensional applications. The equations are advanced in time by the implicit, non-iterative, approximately-factored, finite-difference scheme of Beam and Warming. The numerical stability of the scheme depends on the use of higher-order smoothing terms to damp out higher-frequency oscillations caused by second-order central differencing. The artificial compressibility introduces pressure (sound) waves of finite speed (whereas the speed of sound would be infinite in an incompressible fluid). As the solution converges, these pressure waves die out, causing the derivation of pressure with respect to time to approach zero. Thus, continuity is satisfied for the incompressible fluid in the steady state. Computational efficiency is achieved using a diagonal algorithm. A block tri-diagonal option is also available. When a steady-state solution is reached, the modified continuity equation will satisfy the divergence-free velocity field condition. INS3D is capable of handling several different types of boundaries encountered in numerical simulations, including solid-surface, inflow and outflow, and far-field boundaries. Three machine versions of INS3D are available. INS3D for the CRAY is written in CRAY FORTRAN for execution on a CRAY X-MP under COS, INS3D for the IBM is written in FORTRAN 77 for execution on an IBM 3090 under the VM or MVS operating system, and INS3D for DEC RISC-based systems is written in RISC FORTRAN for execution on a DEC workstation running RISC ULTRIX 3.1 or later. The CRAY version has a central memory requirement of 730279 words. The central memory requirement for the IBM is 150Mb. The memory requirement for the DEC RISC ULTRIX version is 3Mb of main memory. INS3D was developed in 1987. The port to the IBM was done in 1990. The port to the DECstation 3100 was done in 1991. CRAY is a registered trademark of Cray Research Inc. IBM is a registered trademark of International Business Machines. DEC, DECstation, and ULTRIX are trademarks of the Digital Equipment Corporation.

Description: The ability to treat arbitrary boundary shapes is one of the most desirable characteristics of a method for generating grids. 3DGRAPE is designed to make computational grids in or about almost any shape. These grids are generated by the solution of Poisson's differential equations in three dimensions. The program automatically finds its own values for inhomogeneous terms which give near-orthogonality and controlled grid cell height at boundaries. Grids generated by 3DGRAPE have been applied to both viscous and inviscid aerodynamic problems, and to problems in other fluid-dynamic areas. 3DGRAPE uses zones to solve the problem of warping one cube into the physical domain in real-world computational fluid dynamics problems. In a zonal approach, a physical domain is divided into regions, each of which maps into its own computational cube. It is believed that even the most complicated physical region can be divided into zones, and since it is possible to warp a cube into each zone, a grid generator which is oriented to zones and allows communication across zonal boundaries (where appropriate) solves the problem of topological complexity. 3DGRAPE expects to read in already-distributed x,y,z coordinates on the bodies of interest, coordinates which will remain fixed during the entire grid-generation process. The 3DGRAPE code makes no attempt to fit given body shapes and redistribute points thereon. Body-fitting is a formidable problem in itself. The user must either be working with some simple analytical body shape, upon which a simple analytical distribution can be easily effected, or must have available some sophisticated stand-alone body-fitting software. 3DGRAPE does not require the user to supply the block-to-block boundaries nor the shapes of the distribution of points. 3DGRAPE will typically supply those block-to-block boundaries simply as surfaces in the elliptic grid. Thus at block-to-block boundaries the following conditions are obtained: (1) grids lines will match up as they approach the block-to-block boundary from either side, (2) grid lines will cross the boundary with no slope discontinuity, (3) the spacing of points along the line piercing the boundary will be continuous, (4) the shape of the boundary will be consistent with the surrounding grid, and (5) the distribution of points on the boundary will be reasonable in view of the surrounding grid. 3DGRAPE offers a powerful building-block approach to complex 3-D grid generation, but is a low-level tool. Users may build each face of each block as they wish, from a wide variety of resources. 3DGRAPE uses point-successive-over-relaxation (point-SOR) to solve the Poisson equations. This method is slow, although it does vectorize nicely. Any number of sophisticated graphics programs may be used on the stored output file of 3DGRAPE though it lacks interactive graphics. Versatility was a prominent consideration in developing the code. The block structure allows a great latitude in the problems it can treat. As the acronym implies, this program should be able to handle just about any physical region into which a computational cube or cubes can be warped. 3DGRAPE was written in FORTRAN 77 and should be machine independent. It was originally developed on a Cray under COS and tested on a MicroVAX 3200 under VMS 5.1.

Description: In the present paper, direct numerical methods by which to simulate the spatially developing free shear flows in the transitional region are described and the numerical results of a spatially developing plane wake are presented. The incompressible time-dependent Navier-Stokes equations were solved using Pade finite difference approximations in the streamwise direction, a mapped pseudospectral Fourier method in the cross-stream direction, and a third-order compact Runge-Kutta scheme for time advancement. The unstable modes of the Orr-Sommerfeld equations were used to perturb the inlet of the wake. Statistical analyses were performed and some numerical results were compared with experimental measurements. When only the fundamental mode is forced, the energy spectra show amplification of the fundamental and its higher harmonics. In this case, unperturbed alternate vortices develop in the saturation region of the wake. The phase jitter around the fundamental frequency plays a critical role in generating vortices of random shape and spacing. Large- and small-scale distortions of the fundamental structure are observed. Pairing of vortices of the same sign is observed, as well as vortex coupling of vortices of the opposite sign.

Description: This book first reviews the overall aspects and background information related to thermal radiation heat transfer and incorporates new general information, advances in analytical and computational techniques, and new reference material. Coverage focuses on radiation from opaque surfaces, radiation interchange between various types of surfaces enclosing a vacuum or transparent medium, and radiation including the effects of partially transmitting media, such as combustion gases, soot, or windows. Boundary conditions and multiple layers are discussed with information on radiation in materials with nonunity refractive indices.

Description: The present approach to the prediction of instability generation that is due to the interaction of freestream disturbances with regions of subscale variations in surface boundary conditions can account for the finite Reynolds number effects, while furnishing a framework for the study of receptivity in compressible flow and in 3D boundary layers. The approach is illustrated for the case of Tollmien-Schlichting wave generation in a Blasius boundary layer, due to the interaction of a freestream acoustic wave with a localized wall inhomogeneity. Results are presented for the generation of viscous and inviscid instabilities in adverse pressure-gradient boundary layers, supersonic boundary layer instabilities, and cross-flow vortex instabilities.

Description: A novel modeling technique is developed and applied to the long-wavelength gravity, topography, and internal density structure of Venus. The focus of this study is to employ data on the gravity and topography of Venus, which are obtained from the Pioneer Venus Orbiter mission, to enhance understanding some key aspects of the internal structure of the planet. A modeling strategy utilizing 'stokeslets' or basic units of slow viscous flow (governed by the Stokes equations) to duplicate the observed gravity and topography is implemented.

Description: The growth and development of a horseshoe vortex system in an incompressible, three-dimensional turbulent junction flow were investigated experimentally. A streamlined cylinder mounted with its axis normal to a flat surface was used to generate the junction vortex flow. The flow environment was characterized by a body Reynolds number of 183,000, based on the leading edge diameter of the streamlined cylinder. The study included surface flow visualizations, surface pressure measurements, and mean flow measurements of total pressure, static pressure, and velocity distributions in three planes around the base of the streamlined cylinder, and in two planes in the wake flow. Some characterizations of vortex properties based on the measured mean cross-flow velocity components are presented. The results show the presence of a single large, dominant vortex, with strong evidence of a very small corner vortex in the junction between the cylinder and the flat surface. The center of the dominant vortex drifts away from both the body and the flat surface as the flow develops along and downstream of the body. The growth and development of the core of the large, dominant vortex are documented.

Description: The nonlinear resonant-triad interaction, proposed by Raetz (1959), Craik (1971), and others for a Blasius boundary layer, is analyzed here for an adverse-pressure-gradient boundary layer. We assume that the adverse pressure gradient is in some sense weak and, therefore, that the instability growth rate is small. This ensures that there is a well-defined critical layer located somewhere within the flow and that the nonlinear interaction is effectively confined to that layer. The initial interaction is of the parametric resonance type, even when the modal amplitudes are all of the same order. This means that the oblique instability waves exhibit faster than exponential growth and that the growth rate of the two-dimensional mode remains linear. However, the interaction and the resulting growth rates become fully coupled, once oblique-mode amplitudes become sufficiently large, but the coupling terms are now quartic, rather than quadratic as in the Craik (1971) analysis. More importantly, however, new nonlinear interactions, which were not present in the Craik-type analyses, now come into play. These interactions eventually have a dominant effect on the instability wave development.

Description: This algorithm has been developed for calculating both the quantity of compressor bleed flow required to cool a turbine and the resulting decrease in efficiency due to cooling air injected into the gas stream. Because of the trend toward higher turbine inlet temperatures, it is important to accurately predict the required cooling flow. This program is intended for use with axial flow, air-breathing jet propulsion engines with a variety of airfoil cooling configurations. The algorithm results have compared extremely well with figures given by major engine manufacturers for given bulk metal temperatures and cooling configurations. The program calculates the required cooling flow and corresponding decrease in stage efficiency for each row of airfoils throughout the turbine. These values are combined with the thermodynamic efficiency of the uncooled turbine to predict the total bleed airflow required and the altered turbine efficiency. There are ten airfoil cooling configurations and the algorithm allows a different option for each row of cooled airfoils. Materials technology is incorporated and requires the date of the first year of service for the turbine stator vane and rotor blade. The user must specify pressure, temperatures, and gas flows into the turbine. This program is written in FORTRAN IV for batch execution and has been implemented on an IBM 3080 series computer with a central memory requirement of approximately 61K of 8 bit bytes. This program was developed in 1980.

Description: A fast algorithm has been developed for accurately generating boundary-conforming, three-dimensional consecutively refined computational grids applicable to arbitrary wing-body and axial turbomachinery geometries. This algorithm has been incorporated into the GRID3O computer program. The method employed in GRID3O is based on using an analytic function to generate two-dimensional grids on a number of coaxial axisymmetric surfaces positioned between the centerbody and the outer radial boundary. These grids are of the O-type and are characterized by quasi-orthogonality, geometric periodicity, and an adequate resolution throughout the flow field. Because the built-in nonorthogonal coordinate stretching and shearing cause the grid lines leaving the blade or wing trailing-edge to end at downstream infinity, use of the generated grid simplifies the numerical treatment of three-dimensional trailing vortex sheets. The GRID3O program is written in FORTRAN IV for batch execution and has been implemented on an IBM 370 series computer with a central memory requirement of approximately 450K of 8 bit bytes. The GRID3O program was developed in 1981.

Description: As turbine-engine core operating conditions become more severe, designers must develop more effective means of cooling blades and vanes. In order to design reliable, cooled turbine blades, advanced transient thermal calculation techniques are required. The TACT1 computer program was developed to perform transient and steady-state heat-transfer and coolant-flow analyses for cooled blades, given the outside hot-gas boundary condition, the coolant inlet conditions, the geometry of the blade shell, and the cooling configuration. TACT1 can analyze turbine blades, or vanes, equipped with a central coolant-plenum insert from which coolant-air impinges on the inner surface of the blade shell. Coolant-side heat-transfer coefficients are calculated with the heat transfer mode at each station being user specified as either impingement with crossflow, forced convection channel flow, or forced convection over pin fins. A limited capability to handle film cooling is also available in the program. The TACT1 program solves for the blade temperature distribution using a transient energy equation for each node. The nodal energy balances are linearized, one-dimensional, heat-conduction equations which are applied at the wall-outer-surface node, at the junction of the cladding and the metal node, and at the wall-inner-surface node. At the mid-metal node a linear, three-dimensional, heat-conduction equation is used. Similarly, the coolant pressure distribution is determined by solving the set of transfer momentum equations for the one-dimensional flow between adjacent fluid nodes. In the coolant channel, energy and momentum equations for one-dimensional compressible flow, including friction and heat transfer, are used for the elemental channel length between two coolant nodes. The TACT1 program first obtains a steady-state solution using iterative calculations to obtain convergence of stable temperatures, pressures, coolant-flow split, and overall coolant mass balance. Transient calculations are based on the steady-state solutions obtained. Input to the TACT1 program includes a geometrical description of the blade and insert, the nodal spacing to be used, and the boundary conditions describing the outside hot-gas and the coolant-inlet conditions. The program output includes the value of nodal temperatures and pressures at each iteration. The final solution output includes the temperature at each coolant node, and the coolant flow rates and Reynolds numbers. This program is written in FORTRAN IV for batch execution and has been implemented on an IBM 360 computer with a central memory requirement of approximately 480K of 8 bit bytes. The TACT1 program was developed in 1978.

Description: This is a finite-difference program for calculating the viscous compressible boundary layer flow over either planar or axisymmetric surfaces. The flow may be initially laminar and progress through a transitional zone to a fully turbulent flow, or it may remain laminar, depending on the imposed boundary conditions, laws of viscosity, and numerical solution of the momentum and energy equations. The flow may also be forced into a turbulent flow at a chosen spot by the data input. The input may contain factors of arbitrary Reynolds number, free-stream Mach number, free stream turbulence, wall heating or cooling, longitudinal wall curvature, wall suction or blowing, and wall roughness. The solution may start from an initial Falkner-Skan similarity profile, an approximate equilibrium turbulent profile, or an initial arbitrary input profile. This program has been implemented on the IBM 7094/7044 Direct Couple System. This program is written in FORTRAN IV and was developed in 1974.

Description: A computer program has been developed for the design of sharp-edged throat supersonic nozzles where losses are accounted for by correcting the ideal nozzle geometry for boundary layer displacement thickness. The ideal nozzle is designed by the method of characteristics to produce uniform parallel flow at the nozzle exit in the smallest possible distance. Boundary-layer parameters (displacement and momentum thicknesses) are calculated for the ideal nozzle, and the final nozzle geometry is obtained by adding the displacement thickness to the ideal nozzle coordinates. The boundary layer parameters are also used to calculate the aftermixing conditions downstream of the nozzle assuming the flow mixes to a uniform state. The computer program input consists essentially of the nozzle-exit Mach number, specific-heat ratio, nozzle angle, throat half-height, nozzle subsonic section coordinates and corresponding pressure ratios, total temperature and pressure, gas constant, and initial momentum or displacement thickness. The program gas properties are set up for air; for other gases, changes are required to the program. The computer program output consists of the corrected nozzle coordinates, the principal boundary-layer parameters, and the aftermixing conditions. This program has been implemented on the IBM 7094/7044 Direct Couple System.

Description: A computer program has been developed which analyzes by means of mathematical models the temperature profiles in the contents of a filled propellant tank. In designing space vehicles using cryogenic liquid propellants, it is necessary to know how heat transferred from the tank walls and heat absorbed internally affect the temperature distribution with the tank contents. The mathematical flow model is based on results from small-scale experiments. The results showed that when a subcooled fluid is subject to both nonuniform internal heating and wall heating, two distinct temperature regions are developed. In the lower region, the fluid is thoroughly mixed and maintains a uniform temperature profile. In the upper region, a stratified layer develops, and a temperature gradient is formed from the accumulation of warm fluid from the boundary layer along the tank walls; it also indicated that the temperature profiles in the stratified layer exhibited similarity. This concept was developed primarily for internal heating caused by nuclear radiation. However, the theory and computer program are applicable for any form of internal or bulk heating. This program is written in FORTRAN IV for batch execution and has been implemented on the IBM 7094. This program was developed in 1970.

Description: This program solves the two-dimensional, compressible laminar or turbulent boundary-layer equations in an arbitrary pressure gradient. Cohen and Reshotko's method is used for the laminar boundary layer, Sasman and Cresci's method for the turbulent boundary layer, and the Schlichting-Ulrich-Granville method to predict transition. Transition may also be forced at any point by the user. Separation, if it occurs, is predicted for both laminar and turbulent flow. The user may begin values for displacement thickness and momentum thickness in either laminar or turbulent flow. This program was implemented on the IBM 7094.

Description: Reacting free shear layers are of fundamental importance in many industrial systems including gas turbine combustors and rockets. Efficient propulsion systems are essential for air breathing supersonic ramjets in the high Mach number range. A limiting factor in these engines is the time for fuel and oxidizer to mix in the combustion chamber; for fast mixing, the flow must be vigorously turbulent which requires the laminar flow to be unstable. Understanding the stability characteristics of compressible reacting free shear layers is, therefore, very important and may allow one to control the flow. Low speed shear layers are highly unstable but, as chemical reaction and compressibility effects tend to stabilize them, it is important to investigate the stability of high speed reacting mixing layers. The latter consists of two fluid streams containing fuel and oxidizer respectively, and the conclusions are expected to apply, with quantitative modifications, to other shear flows, e.g., jets. Since low speed reacting cases have been studied earlier, we concentrate on the effects of Mach number and heat release. We are primarily interested in solving the stability problem over a large range of Mach number and heat release. In order to understand the effect of the heat release on the stability of this flow, one must first study the characteristics of the non-reacting flow. Inviscid theory is a reliable guide for understanding stability of compressible shear flows at moderate and large Reynolds numbers and is the basis for this work.

Description: Several direct numerical simulations of high-speed turbulent Couette flow were performed with a new spectral code. Mach numbers up to three and a Reynolds number of 3000 were used. A new time-integration scheme was developed to handle Mach numbers above 1.5, which require greater accuracy and stability than lower Mach numbers. At low Mach number, the large streamwise eddies found by M. J. Lee in high incompressible Couette flow simulations were reproduced. At higher Mach numbers these structures still exist, but they become considerably less organized (although the disorganization may be a function of the spanwise box size). While the same types of vortical structures seen in the incompressible flow are observed at higher Mach numbers, a new structure involving the divergence of the velocity is also observed. This structure is generally associated with low shear areas next to the walls, but it has not been determined whether it is a cause or an effect of the low shear. A 'nonphysical' simulation was performed to determine by what mechanism the Mach number affects the flow. It appears that pressure gradient (acoustic) effects are more important than variable viscosity effects in determining the wall shear, but the size of vortical structures is determined more by the local kinematic viscosity. Low-order mean statistics are provided to help quantify these effects.

Description: Many of the turbulent layers encountered in practical flows develop in adverse pressure gradients; hence, the dynamics of the thickening and possible separation of the boundary layer has important implications for design practices. What are the key physical processes that govern how a turbulent boundary layer responds to an adverse pressure gradient, and how should these processes be modeled? Despite the ubiquity of such flows in engineering and nature, these equations remain largely unanswered. The turbulence closure models presently used to describe these flows commonly use 'wall functions' that have ad hoc corrections for the effects of pressure gradients. There is, therefore, a practical and theoretical need to examine the effects of adverse pressure gradients on wall bounded turbulent flows in order to develop models based on sound physical principle. The evolution of a turbulent boundary layer on a flat wall with an externally imposed pressure gradient is studied.

Description: Advancing the knowledge and understanding of turbulence theory is addressed. Specific problems to be addressed will include studies of subgrid models to understand the effects of unresolved small scale dynamics on the large scale motion which, if successful, might substantially reduce the number of degrees of freedom that need to be computed in turbulence simulation.

Description: The increase in the range of length scales with increasing Reynolds number limits the direct simulation of turbulent flows to relatively simple geometries and low Reynolds numbers. However, since most flows of engineering interest occur at much higher Reynolds number than is currently within the capabilities of full simulation, prediction of these flow fields can only be obtained by solving some suitably-averaged set of governing equations. In the traditional Reynolds-averaged approach, the Navier-Stokes equations are averaged over time. This in turn yields correlations between various turbulence fluctuations. It is these terms, e.g. the Reynolds stresses, for which a turbulence model must be derived. Turbulence modeling of incompressible flows has received a great amount of attention in the literature. An area of research that has received comparatively less attention is the modeling of compressible turbulent flows. An approach to simulating compressible turbulence at high Reynolds numbers is through the use of Large-Eddy Simulation (LES). In LES the dependent variables are decomposed into a large-scale (resolved) component and a sub-grid scale component. It is the small-scale components of the velocity field which are presumably more homogeneous than the large scales and, therefore, more easily modeled. Thus, it seems plausible that simpler models, which should be more universal in character than those employed in second-order closure schemes, may be developed for LES of compressible turbulence. The objective of the present research, therefore, is to explore models for the Large-Eddy Simulation of compressible turbulent flows. Given the recent successes of Zeman in second order closure modeling of compressible turbulence, model development was guided by principals employed in second-order closures.

Description: With the recent revitalization of high speed flow research, compressibility presents a new set of challenging problems to turbulence researchers. Questions arise as to what extent compressibility affects turbulence dynamics, structures, the Reynolds stress-mean velocity (constitutive) relation, and the accompanying processes of heat transfer and mixing. In astrophysical applications, compressible turbulence is believed to play an important role in intergalactic gas cloud dynamics and in accretion disk convection. Understanding and modeling of the compressibility effects in free shear flows, boundary layers, and boundary layer/shock interactions is discussed.

Description: Energy of a high harmonic standing wave tends to be distributed equally over the whole wave even in a stratified medium where the wave's peak amplitude can be much larger near the upper boundary than the lower one. This fact is generalized to the many diverse physical problems which solve second-order differential equations of Sturm-Liouville type. For any such solution y(z) whose sign fluctuates along the z-axis, quantities are found which have the same value between any two neighboring zeros of y. One of the equidistributed quantities for an oscillating fluid sphere is similar to kinetic energy but is identical only in limiting cases. The acoustic midpoint of a cavity can be a unique place where some nonlinear perturbations have extra strength. This may apply to the puzzling solar phenomenon called supergranulation.

Description: A computational code for the vorticity-potential method is developed for a three-dimensional bounded vorticity field. The evaluation of the boundary data for the vector potential in the code is improved so that the numerical solution simulates that in an unbounded domain to a high order. The time evolution of two vortex rings and that of an elliptic ring are investigated with this code. The cut-and-connect phenomena of vortex rings are successfully captured. The results are compared with those of asymptotic theory and the experiment. They also highlight the need for additional theoretical and numerical investigations.

Description: The interaction of the thermosphere and ionosphere is largely governed by collisions between ions and neutral particles. On Venus and the Earth, O(+) is a dominant ion, and atomic O dominates throughout much of the thermosphere; therefore an accurate O(+)-O cross section is an important prerequisite for understanding the dynamics of planetary upper atmospheres. The cross section and momentum transfer collision frequency are calculated with a quantum mechanical code which includes resonance charge exchange, polarization, and charge-quadrupole effects. Our results agree well with earlier calculations of Stubbe (1968) and Stallcop et al. (1991).

Description: Synthetic aperture radar (SAR) focusing can be achieved either based on accurate ephemeris data or on an autofocusing process. For the Magellan project, such a decision must be made in the early phase of Magellan SAR system design. The analysis of the emphemeris requirement is complicated. The analysis given by the author leads to the conclusion that empheris data obtained from the Magellan navigation system provide sufficient accuracy to meet the Magellan image resolution requirement.

Description: An overview of the history and current status of research on planetary noble gases is presented. The discovery that neon and argon are vastly more abundant on Venus than on earth points to the solar wind rather than condensation as the fundamental process for placing noble gases in the atmospheres of the terrestrial planets; however, solar wind implantation may not be able to fully reproduce the observed gradient, nor does it obviously account for similar planetary Ne/Ar ratios and dissimilar planetary Ar/Kr ratios. More recent studies have emphasized escape rather than accretion. Hydrodynamic escape, which is fractionating, readily accounts for the difference between atmospheric neon and isotopically light mantle neon. Atmospheric cratering, which is nearly nonfractionating, can account for the extreme scarcity of nonradiogenic noble gases (and other volatiles) on Mars.

Description: An account is given of interface-driven motions of drops and bubbles. It is shown that even in the simplest cases, theory predicts exotic flow topologies. Attention is given to several unsolved problems that must be addressed both theoretically and experimentally.

Description: The study points out three difficulties with the hypothesis of Hartung (1993) that the Corvid meteor system, observed only once in 1937, may be the return of the clump of ejecta from the formation of a lunar crater, specifically, an event recorded in the chronicles of Gervase of Canterbury on June 25, 1178, which Hartung (1976) previously suggested may be an eyewitness account of the formation of the lunar crater Giordano Bruno. On the basis of this, he predicts that another Corvid shower may be observed in 2003 or 2006. This hypothesis is rejected on the contention that it is implausible that a clump of ejecta could be launched into heliocentric orbit with a low enough dispersion in velocity among separate pieces that it would produce a meteor shower in just one year and not others. Subsequent perturbations by the earth on parts of the clump passing near the earth but not impacting would destroy the coherence of the clump on a time scale much shorter than the 759-yr interval proposed. There are so many orbits that could yield a shower after 759 yr that it is unlikely that a prediction of a return in a specific year would be correct.

Description: A new gravity field determination that has been produced combines both the Pioneer Venus Orbiter (PVO) and the Magellan Doppler radio data. Comparisonsbetween this estimate, a spherical harmonic model of degree and order 21, and previous models show that significant improvements have been made. Results are displayed as gravity contours overlaying a topographic map. We also calculate a new spherical harmonic model of topography based on Magellan altimetry, with PVO altimetry included where gaps exist in the Magellan data. This model is also of degree and order 21, so in conjunction with the gravity model, Bouguer and isostatic anomaly maps can be produced. These results are very consistent with previous results, but reveal more spatial resolution in the higher latitudes.

Description: Results of a comprehensive morphological analysis of the dayside superthermal ion presence in the dayside ionopause region of Venus made using the Orbiter Ion Mass Spectrometer measurements are presented. There is a strong inbound-outbound difference in both the number and intensity of the superthermal ions with most events detected on the inbound crossing of the ionopause. The ambient energy of the detected ions is comparable to the ram energy of cold ions relative to the moving spacecraft. A comparison of the locations of the superthermal O(+) ion occurrences with respect to other parameters shows that the superthermals are a characteristic of the ionopause transition separating the essentially stationary ionosphere plasma from the flowing ionosheath plasma as deduced by Taylor et al. (1980) from individual examples.

Description: Compared to most other Yamato polymict eucrites, Yamato Y792769 eucrite includes fewer and smaller eucritic clasts with homogenized pyroxenes, and its fine-grained matrix is shock-compacted and sintered. In this work, the relationships between the Antarctic eucrite Y792769, monomict eucrites, polymict eucrites, and isotopic ages are investigated, using results of Ar-39/Ar-40 method to date the time of the major thermal event on the Y792769 body and the Rb-Sr and Sm-Nd methods to determine whether relict older ages might have been preserved in some of the breccia materials. The Ar-39/Ar-40 time of the last thermal event which produced the Y792769 texture is 3.99 +/- 0.04 Ga. The complete resetting of the Ar-39/Ar-40 age is consistent with the texture of Y792769 observed in TEM, suggesting that shock compaction converted part of the matrix plagioclase to maskelynite. The Sm-Nd data give an age of 4.23 +/- 0.12 Ga, reflecting partial resetting of the Sm-Nd system during breccia formation. The 3.9 Ga Ar-39/Ar-40 age probably reflects a period of intense meteoroid bompardment which affected the entire inner solar system.

Description: In three-dimensional numerical simulations of a rapidly rotating Boussinesq fluid shell, thermally driven convection in the form of columns parallel to the rotation axis generates an alternately directed mean zonal flow with a cylindrical structure. The mean structure at the outer spherical surface consists of a broad eastward flow at the equator and alternating bands of westward and eastward flows at higher latitudes in both hemispheres. The banded structure persists even though the underlying convective motions are time-dependent. These results, although still far from the actual motions seen on Jupiter and Saturn, provide support for theoretical suggestions that thermal convection can account for the remarkable banded flow structures on these planets.

Description: A review of the planetary system is presented from a uniform perspective, taking advantage of the data obtained during the planetary space missions. A comparative approach is used, which makes it possible to consider each planet or a satellite in context, focusing both on their common origin and on the processes that have influenced their evolution. The discussion focuses on the Jovian planets, the moon and the Mercury, Venus and Earth, and small bodies in the outer solar system. Attention is also given to comets, asteroids, and cosmic impacts.

Description: An electron probe microanalyzer is used to measure the Mn, Fe, and oxygen zoning profiles of olivines in the ALHA 77257 ureilite. This is done to study the effects of reduction on the Mn-Fe value, as ureilite olivines exhibit thin reduced rims. Since the Mn content gradually increases toward the rim of ureilite olivines, while the Fa (= 100 x Fe/(Mg + Fe), mol percent) component decreases, the Mn-Fe content of olivine is likely related to redox conditions. The results of melting experiments suggest that the Mn-Fe positive correlation is related to temperature and that the negative correlation of Mn-Fe in olivine and low-Ca pyroxene is related to reduction.

Description: Voyager IRIS observations of the NEB reveal longitudinal variability of 5-micron brightness temperatures of order 100 C. An anisotropic multiple scattering radiative transfer model is used to calculate synthetic spectra for comparison with the IRIS observations. Mie theory is used to model the spectral dependence of cloud extinction from 180 to 2300/cm. Cloud base locations within the model vary with assumed gas abundances according to thermochemical equilibrium. It is found that spatial variations in the abundance profiles of the condensible species, parahydrogen profiles and cloud optical depths can be used as tracers of the local and large-scale dynamics. The variation of cloud opacity is strongly correlated with the variation of relative humidity, which suggests that dynamic depletion of water vapor above the cloud forming level is the most plausible model to explain the spatial variation in the water profile within the NEB.

Description: In the present grid-independent approximate Riemann solver for 2D and 3D flows that are governed by the Euler or Navier-Stokes equations, fluxes on grid faces are obtained by wave decomposition; the assumption of information-propagation in the velocity-difference directions leads to a more accurate resolution of shear and shock waves, when these are are oblique to the grid. The model, which yields significantly greater accuracy in both supersonic and subsonic first-order spatially accurate computations, describes the difference in states at each grid interface by the action of five waves.

Description: Io's brightness was monitored following emergence from eclipse by Jupiter on 14 occasions during the Jupiter apparitions from 1981 to 1989 and no instance is found of what has been called posteclipse brightening. If all the observations are averaged, a 2 percent effect cannot be ruled out; however, this effect is also the size of the rms errors of the summed data set. If condensation followed by sublimation of SO2 frost is hypothesized to be the mechanism which causes posteclipse brightening, then solar insolation alone may not be sufficient to remove an optically thick layer of SO2 frost in the time required to conform with the positive reports of posteclipse brightening in the literature. An additional source of energy is required, such as warming caused by the impact of magnetospheric particles as they dissipate the energy of their motion into Io's surface.

Description: For a fixed heat flow, the surface flow velocity of a convecting layer is not strongly sensitive to the variation of viscosity as a function of depth. Thus, the inferred absence of a low viscosity asthenosphere on Venus can not account for the limited surface motions there. The surface velocity is dependent on the convective geometry. Cartesian geometry convection can produce large surface velocities if the high viscosity surface layer is broken in places by weak zones. On the other hand, a high viscosity surface layer may inhibit the development of large surface velocities in axisymmetric convection.

Description: The paper characterizes the mineralogy and spectral properties of three slightly palagonitized basaltic tephra samples collected near the summit of Mauna Kea in order to contribute to the basis for inferring mineralogy and processes for Martian surface materials. The mineralogy of size fractions of these samples is examined by diffuse reflectance and FIR spectroscopy, optical microscopy, X-ray diffraction, Moessbauer spectroscopy, magnetic analysis, EMPA, TEM, and SEM. For the 20-1000 micron size fraction, sample HWMK11 (red) is essentially completely oxidized and has a hematite (Ti-hematite) pigment dispersed throughout the silicate matrix. Sample HWMK12 (black) has the lowest proportion of ferric-bearing phases and is thus least weathered. For HWMK11, the amount of hematite is essentially constant, and mica is present only in the coarse clay-sized fraction; smectites are low in structural Fe.

Description: The action of weak, streamwise vortices on a plane, incompressible, steady mixing layer is examined in the large Reynolds number limit. The outer, inviscid region is bounded by a vortex sheet to which the viscous region is confined. It is shown that the local linear analysis becomes invalid at streamwise distances O(epsilon sup -1), where (epsilon much less than 1) is the crossflow amplitude, and a new nonlinear analysis is constructed for this region. Numerical solutions of the nonlinear problem show that the vortex sheet undergoes an O(1) change in position and that the solution is ultimately terminated by a breakdown in the numerical procedure. The corresponding viscous layer shows downstream thickening, but appears to remain well behaved up to the terminal location.

Description: This paper investigates the linear stability of confined mixing layers with special emphasis on the effects of heat release and compressibility. The results show that reflection of supersonic disturbances by the walls makes the confined supersonic mixing layer more unstable than the unconfined free shear layer. Decreasing the distance between the walls makes the flow more unstable. However, subsonic disturbances are relatively unaffected by the walls. Heat release and Mach number hardly change the growth rates of supersonic disturbances. The most unstable supersonic disturbances are two-dimensional in rectangular channel flows, but three-dimensional in partially confined flows. Finally, the reactants are not strongly mixed by supersonic instabilities, which mainly disturb one side of the layer.

Description: Recent laboratory polarimetric measurements of light scattering by layers of alumina grains of different size are used to examine the relevance of the Wolff & Dollfus (1990) theory of negative polarization, which has been developed to describe light scattering by regolithic grains of arbitrary size. It is demonstrated that, while reproducing satisfactorily the main polarization characteristics of light scattering by grains much greater than the wavelength, this theory does not describe the negative polarization of light measured by Geake & Geake (1990) for subwavelength-sized grains. The smallest grains exhibit a polarization feature that may be called the 'polarization opposition effect'. This feature requires an alternative explanation and may be due to the so-called weak localization of photons in discrete random media.

Description: It is shown that N2 may be present in the troposphere of Neptune in an amount difficult to evaluate but which could easily be as high as 0.003, while there is no evidence that it is present in the atmosphere of Uranus. The estimate of the helium abundance depends on the assumed value for N2. If there is no N2 in the observed region of the atmosphere of Uranus and an N2 mole fraction of 0.003 on Neptune, the central value of the estimates of the helium abundance are equal to 0.26 by mass in both planets, which is close to the protosolar value of 0.28. This would imply that the He/H2 ratios measured in the outer atmospheres of Uranus and Neptune are representative of the ratio in the primitive solar nebula and thus were not modified during planetary formation.

Description: The GISS GCM is used here to examine the hypothesis that equatorial superrotation on slowly rotating planets is sensitive to the nature of the vertical radiative heating profile and can exist in the absence of diurnally varying forcing. The general circulation, diabatic heating, and thermal structure produced in the experiments are described and the heat and angular momentum budgets and energy cycles are analyzed to understand the factors conducive to equatorial superrotation. The implications of the results for future planetary missions are addressed.

Description: The spins of the terrestrial planets likely arose as the planets formed by the accretion of planetesimals. Depending on the masses of the impactors, the planet's final spin can either be imparted by many small bodies (ordered accretion), in which case the spin is determined by the mean angular momentum of the impactors, or by a few large bodies (stochastic accretion), in which case the spin is a random variable whose distribution is determined by the rms angular momentum of the impactors. In the case of ordered accretion, the planet's obliquity is expected to be near 0 or 180 deg, whereas, if accretion is stochastic, there should be a wide range of obliquities. Analytic arguments and extensive orbital integrations are used to calculate the expected distributions of spin rate and obliquity as a function of the planetesinal mass and velocity distributions. The results imply that the spins of the terrestrial planets are determined by stochastic accretion.

Description: Direct numerical simulations were made to examine the local structure of the reaction zone for a moderately fast reaction between unmixed species in decaying, homogeneous turbulence and in a homogeneous turbulent shear flow. Pseudospectral techniques were used in domains of 64 exp 3 and higher wavenumbers. A finite-rate, single step reaction between non-premixed reactants was considered, and in one case temperature-dependent Arrhenius kinetics was assumed. Locally intense reaction rates that tend to persist throughout the simulations occur in locations where the reactant concentration gradients are large and are amplified by the local rate of strain. The reaction zones are more organized in the case of a uniform mean shear than in isotropic turbulence, and regions of intense reaction rate appear to be associated with vortex structures such as horseshoe vortices and fingers seen in mixing layers. Concentration gradients tend to align with the direction of the most compressive principal strain rate, more so in the isotropic case.