ALBERT

Description: Periodic Comet Shoemaker/Levy 9 will impact Jupiter in late July 1994. The comet, which broke into more than 20 telescopically detectable fragments when it passed with the Roche lobe of Jupiter on July 8, 1992, is captured in a highly eccentric orbit about Jupiter. The 21 recognized nuclei will be spread out in a train of the order 7 x 10(exp 6) km long at the time of impact, and the impacts will be spread in time over about 5 1/2 days centered on about July 21.2 UT. In addition to the train of recognized bright nuclei, the comet consists of 'wings' of unresolved bodies that are the source of a very broad composite dust tail. The linear extent of the wings is about an order of magnitude greater than that of the train of recognized discrete nuclei. Collision of the wings will be spread in time over several months. Thus the impact of P/S-L 9 with Jupiter will be an event of appreciable duration.

Description: Approximately 950 impact craters have been identified on the surface of Venus, mainly in Magellan radar images. From a combination of Earth-based Arecibo, Venera 15/1, and Magellan radar images, we have interpreted 72 as unequivocal peak-ring craters and four as multiringed basins. The morphological and structural preservation of these craters is high owing to the low level of geologic activity on the venusian surface (which is in some ways similar to the terrestrial benthic environment). Thus these craters should prove crucial to understanding the mechanics of ringed crater formation. They are also the most direct analogs for craters formed on the Earth in Phanerozoic time, such as Chicxulub. We summarize our findings to date concerning these structures.

Description: We propose sending a balloon-borne UV photometer sensor package to measure atmospheric ozone on Mars, and this package could be a Discovery Program sensor candidate. Past measurements of ozone on Mars are highly uncertain, perhaps a factor of 3 or so uncertain, due primarily to interference and masking by cloud and dust. In-situ balloon measurements would avoid these problems, and would provide 'ground truth' which remote sensing techniques cannot. We have explored this approach to measure ozone abundance in the terrestrial stratosphere with a balloon-borne UV absorption photometer. Atmospheric pressures and temperatures and ozone concentrations near the surface of Mars are similar to those in the terrestrial stratosphere.

Description: The Hydroblast Research Cell at Marshall Space Flight Center is used to investigate the use of high pressure waterjets to strip paint, grease, adhesive and thermal spray coatings from various substrates. Current methods of cleaning often use ozone depleting chemicals (ODC) such as chlorinated solvents. High pressure waterjet cleaning has proven to be a viable alternative to the use of solvents. A popular method of waterjet cleaning involves the use of a rotating, multijet, high pressure water nozzle which is robotically controlled. This method enables rapid cleaning of a large area, but problems such as incomplete coverage and damage to the substrate from the waterjet have been observed. This report summarizes research consisting of identifying and investigating the basic properties of rotating, multijet, high pressure water nozzles, and how particular designs and modes of operation affect such things as stripping rate, standoff distance and completeness of coverage. The study involved computer simulations, an extensive literature review, and experimental studies of different nozzle designs.

Description: An integrated moon program has often been proposed as a logical next step for today's space efforts. In the context of preparing for the possibility of launching a moon program, the European Space Agency is currently conducting an internal study effort which is focusing on the assessment of key technologies. Current thinking has this moon program organized into four phases. Phase 1 will deal with lunar resource exploration. The goal would be to produce a complete chemical inventory of the moon, including oxygen, water, other volatiles, carbon, silicon, and other resources. Phase 2 will establish a permanent robotic presence on the moon via a number of landers and surface rovers. Phase 3 will extend the second phase and concentrate on the use and exploitation of local lunar resources. Phase 4 will be the establishment of a first human outpost. Some preliminary work such as the building of the outpost and the installation of scientific equipment will be done by unmanned systems before a human crew is sent to the moon.

Description: A new model for the two-point vector stream function correlation has been developed using tensor invariant arguments and evaluated by the comparison of model predictions with DNS data for incompressible homogeneous turbulent shear flow. This two-point vector stream function model correlation can then be used to calculate the two-point velocity correlation function and other quantities useful in turbulence modeling. The model assumes that the two-point vector stream function correlation can be written in terms of the separation vector and a new tensor function that depends only on the magnitude of the separation vector. The model has a single free model coefficient, which has been chosen by comparison with the DNS data. The relative error of the model predictions of the two-point vector stream function correlation is only a few percent for a broad range of the model coefficient. Predictions of the derivatives of this correlation, which are of interest in turbulence modeling, may not be this accurate.

Description: A non-local model is presented for approximating the pressure diffusion in calculations of turbulent free shear and boundary layer flows. It is based on the solution of an elliptic relaxation equation which enables local diffusion sources to be distributed over lengths of the order of the integral scale. The pressure diffusion model was implemented in a boundary layer code within the framework of turbulence models based on both the kappa-epsilon-(bar)upsilon(exp 2) system of equations and the full Reynolds stress equations. Model computations were performed for mixing layers and boundary layer flows. In each case, the pressure diffusion model enabled the well-known free-stream edge singularity problem to be eliminated. There was little effect on near-wall properties. Computed results agreed very well with experimental and DNS data for the mean flow velocity, the turbulent kinetic energy, and the skin-friction coefficient.

Description: The structure of the separated flow region over a 20 degree compression corner at a free-stream Mach number of 2.84 is investigated computationally using a Reynolds averaged Navier Stokes (R.A.N.S.) solver and kappa-epsilon model. At this Mach number and ramp angle, a steady-state recirculation region of order delta(sub o) is observed, with onset of a 'plateau' in the wall pressure distribution near the corner. At lower ramp angles, separation is negligible, while at an angle of 24 degrees, separation regions of length 2 delta(sub o) are expected. Of interest here is the response of the mathematical model to inclusion of the pressure dilatation term for turbulent kinetic energy. Compared with the experimental data of Smits and Muck (1987), steady-state computations show improvement when the pressure dilatation term is included. Unsteady computations, using both unforced and then forced inlet conditions, did not predict the oscillatory motion of the separation bubble as observed in laboratory experiments. An analysis of the separation bubble oscillation and the turbulent boundary layer (T.B.L.) frequencies for this flow suggests that the bubble oscillations are of nearly the same order as the turbulent frequencies, and therefore difficult for the model to separate and resolve.

Description: A direct numerical approach has been developed to simulate supersonic turbulent boundary layers. The mean flow quantities are obtained by solving the parabolized Reynolds-averaged Navier-Stokes equations (globally). Fluctuating quantities are computed locally with a temporal direct numerical simulation approach, in which nonparallel effects of boundary layers are partially modeled. Preliminary numerical results obtained at the free-stream Mach numbers 3, 4.5, and 6 with hot-wall conditions are presented. Approximately 5 million grid points are used in all three cases. The numerical results indicate that compressibility effects on turbulent kinetic energy, in terms of dilatational dissipation and pressure-dilatation correlation, are small. Due to the hot-wall conditions the results show significant low Reynolds number effects and large streamwise streaks. Further simulations with a bigger computational box or a cold-wall condition are desirable.

Description: A major problem in modeling of turbulent supersonic flows is the correct assessment of viscous-inviscid interaction problems. Of particular interest is the interaction of boundary layers with shocks. Present turbulence models give in most cases unsatisfactory results in the region of rapid distortion and in the separation region (if one is present) in particular with regard to mean flow profiles and turbulence quantities. The objective of the present work is the direct numerical simulation of shock boundary layer interaction. This report summarizes the first phase during which a numerical method suitable for this problem has been developed and a computer code has been written and tested.

Description: Since high resolution numerical simulations of turbulent flows, or experiments at high Reynolds numbers, represent a substantial investment in resources, the estimation of the minimum resolution required for the study of a given property has been the subject of continued interest. Different properties require, in general, different resolutions, and the present paper is dedicated to the requirements for the measurement of the probability distribution functions of the velocity gradients and, in particular, of their low order moments. The deviation of these quantities from the values corresponding to a Gaussian distribution was one of the first indications of the presence of Reynolds number-dependent intermittency and has been the object of recent interest as numerical simulations have become able to explore the distribution of gradients in the low Reynolds number range, while new experiments have extended the range to increasingly high Reynolds numbers. We will use progressive filtering of the results of numerical simulations of isotropic turbulence as a model for the effect of a sensor of finite size. The numerical issues will be addressed first to ensure that the simulations are fully resolved from the point of view of the velocity gradients. This will also give us an estimate for the numerical resolution required for the different quantities.

Description: The MSU compact four-wire transverse vorticity probe permits omega(sub z)(t) measurements in a nominally 1 sq mm domain. Note that a conventional coordinate system is used with x and y in the streamwise and normal directions respectively. The purpose of this investigation was to acquire time series data in the same access port at the ceiling of the 80 ft x 120 ft wind tunnel (NASA Ames Research Center) as earlier used by the Wallace group from the University of Maryland and to compare the present results with those of the three-component vorticity probe used in that earlier study.

Description: Boundary layer flow transition has long been suggested as a potential noise source in both marine (sonar-dome self noise) and aeronautical (aircraft cabin noise) applications, owing to the highly transient nature of process. The design of effective noise control strategies relies upon a clear understanding of the source mechanisms associated with the unsteady flow dynamics during transition. Due to formidable mathematical difficulties, theoretical predictions either are limited to early linear and weakly nonlinear stages of transition, or employ acoustic analogy theories based on approximate source field data, often in the form of empirical correlation. In the present work, an approach which combines direct numerical simulation of the source field with the Lighthill acoustic analogy is utilized. This approach takes advantage of the recent advancement in computational capabilities to obtain detailed information about the flow-induced acoustic sources. The transitional boundary layer flow is computed by solving the incompressible Navier-Stokes equations without model assumptions, thus allowing a direct evaluation of the pseudosound as well as source functions, including the Lighthill stress tensor and the wall shear stress. The latter are used for calculating the radiated pressure field based on the Curle-Powell solution of the Lighthill equation. This procedure allows a quantitative assessment of noise source mechanisms and the associated radiation characteristics during transition from primary instability up to the laminar breakdown stage. In particular, one is interested in comparing the roles played by the fluctuating volume Reynolds stress and the wall-shear-stresses, and in identifying specific flow processes and structures that are effective noise generators.

Description: Turbulence modeling continues to be one of the most difficult problems in fluid mechanics. Existing prediction methods are well developed for certain classes of simple equilibrium flows, but are still not entirely satisfactory for a large category of complex non-equilibrium flows found in engineering practice. Direct and large-eddy simulation (LES) approaches have long been believed to have great potential for the accurate prediction of difficult turbulent flows, but the associated computational cost has been prohibitive for practical problems. This remains true for direct simulation but is no longer clear for large-eddy simulation. Advances in computer hardware, numerical methods, and subgrid-scale modeling have made it possible to conduct LES for flows or practical interest at Reynolds numbers in the range of laboratory experiments. The objective of this work is to apply ES and the dynamic subgrid-scale model to the flow of a boundary layer over a concave surface.

Description: Previously the theoretical development leading to the dynamic localization model (DLM) for large-eddy simulation (LES) was presented. The method has been successfully applied to isotropic turbulence, channel flow, and the flow over a backward-facing step. Here we apply the model to the computation of the temporally developing place wake. The two main objectives of this project are: (1) Use the model to perform an LES of a time developing plane wake and compare the results with direction numerical simulation (DNS) data to see if important statistical measures can be readily predicted, and to provide a relative evaluation of the several versions of the model in terms of predictive capability and cost; and (2) If the tests in (1) show that the model generates reliable predictions, then use the LES to study various aspects of the physics of turbulent wakes and mixing layers.

Description: The modeling of subgrid scales in large-eddy simulation (LES) has been rationalized by the introduction of the dynamic localization procedure. This method allows one to compute rather than prescribe the unknown coefficients in the subgrid-scale model. Formally, the LES equations are supposed to be obtained by applying to the Navier-Stokes equations a 'grid filter' operation. Though the subgrid stress itself is unknown, an identity between subgrid stresses generated by different filters has been derived. Although preliminary tests of the Dynamic Localization Model (DLM) with k-equation have been satisfactory, the use of a negative eddy viscosity to describe backscatter is probably a crude representation of the physics of reverse transfer of energy. Indeed, the model is fully deterministic. Knowing the filtered velocity field and the subgrid-scale energy, the subgrid stress is automatically determined. We know that the LES equations cannot be fully deterministic since the small scales are not resolved. This stems from an important distinction between equilibrium hydrodynamics and turbulence. In equilibrium hydrodynamics, the molecular motions are also not resolved. However, there is a clear separation of scale between these unresolved motions and the relevant hydrodynamic scales. The result of molecular motions can then be separated into an average effect (the molecular viscosity) and some fluctuations. Due to the large number of molecules present in a box with size of the order of the hydrodynamic scale, the ratio between fluctuations and the average effect should be very small (as a result of the 'law of large numbers'). For that reason, the hydrodynamic balance equations are usually purely deterministic. In turbulence, however, there is no clear separation of scale between small and large eddies. In that case, the fluctuations around a deterministic eddy viscosity term could be significant. An eddy noise would then appear through a stochastic term in the subgrid-scale model and could be the source of backscatter.

Description: The elliptic relaxation model of Durbin (1993) for modeling near-wall turbulence using second moment closures (SMC) is compared to DNS data for a channel flow at Re(sub t) = 395. The agreement for second order statistics and even the terms in their balance equation is quite satisfactory, confirming that very little viscous effects (via Kolmogoroff scales) need to be added to the high Reynolds versions of SMC for near-wall-turbulence. The essential near-wall feature is thus the kinematic blocking effect that a solid wall exerts on the turbulence through the fluctuating pressure, which is best modeled by an elliptic operator. Above the transition layer, the effect of the original elliptic operator decays rapidly, and it is suggested that the log-layer is better reproduced by adding a non-homogeneous reduction of the return to isotropy, the gradient of the turbulent length scale being used as a measure of the inhomogeneity of the log-layer. The elliptic operator was quite easily applied to the non-linear Craft & Launder pressure-strain model yielding an improved distinction between the spanwise and wall normal stresses, although at higher Reynolds number (Re) and away from the wall, the streamwise component is severely underpredicted, as well as the transition in the mean velocity from the log to the wake profiles. In this area a significant change of behavior was observed in the DNS pressure-strain term, entirely ignored in the models.

Description: This work carries earlier finite-difference calculations of the Widnall instability of vortex rings into the late non-linear stage. Plots of energy in azimuthal Fourier modes indicate that low-order modes dominate at large times; their structure and dynamics remain unexplored, however. An attempt was made to calculate the acoustic signal using the theory of Mohring (1978), valid for unbounded flow. This theory shows that only low-order azimuthal modes contribute to the sound. As a check on the effects of axial periodicity and a slip wall at large radius imposed by the numerical scheme, the acoustic integrals were also computed in a truncated region. Half of the terms contributing to the sound have large differences between the two regions, and the results are therefore unreliable. The error is less severe for a contribution involving only the m = 2 mode, and its low frequency is consistent with a free elliptic bending wave on a thin ring.

Description: The spectra of the pressure, and other higher-order quantities including the dissipation, the enstrophy, and the square of the longitudinal velocity derivative are computed using data obtained from direct numerical simulation of homogeneous isotropic turbulence at Taylor-Reynolds numbers R(sub lambda) in the range 38 - 170. For the pressure spectra we find reasonable collapse in the dissipation range (of the velocity spectrum) when scaled in Kolmogorov variables and some evidence, which is not conclusive, for the existence of a k(exp -7/3) inertial range where k = absolute value of K, is the modulus of the wavenumber. The power spectra of the dissipation, the enstrophy, and the square of the longitudinal velocity derivative separate in the dissipation range, but appear to converge together in the short inertial range of the simulations. A least-squares curve-fit in the dissipation range for one value of R(sub lambda) = 96 gives a form for the spectrum of the dissipation as k(exp 0)exp(-Ck eta), for k(eta) greater than 0.2, where eta is the Kolmogorov length and C is approximately equal to 2.5.

Description: Panel methods are commonly used in computational fluid dynamics for the solution of potential flow problems. The methods are a numerical technique based on the surface distribution of singularity elements. The solution is the process of finding the strength of the singularity elements distributed over the body's surface. This process involves the solution of the matrix problem Pq = p' for a set of unknowns q. The Fast Multipole Method is used to directly compute q without using matrix solvers. The algorithm works in O(N) time for N points, a great improvement over standard matrix solvers. In panel methods, the surface of a body is divided into a series of quadrilateral panels. The methods involve the computation of the influence of all other panels on each individual panel. The influence is based on the surface distribution, though this can be approximated by the area for distant panels. An alternative approximation, though with arbitrary accuracy, is to develop a multipole expansion about the center of the panel to describe the effect of a given panel on distant points in space. The expansion is based on the moments of the panel, thus allow the use of various surface distributions without changing the basic algorithm, just the computation of the various moments. The expansions are then manipulated in a tree walk to develop Taylor series expansions about a point in space which describe the effect of all distant panels on any point within a volume of convergence. The effect of near panels then needs to be computed directly, but the effect of all distant panels can be computed by simply evaluating the resulting expansion. The Fast Multipole Method has been applied to panel methods for the solution of source and doublet distributions. A major feature of the algorithm is that the algorithm does not change to derive the potential and velocity for sources and doublets. The same expansions can be used for both sources and doublets. Since the velocity is related to the potential, and the doublet potential is related to the z-component of the source velocity, all values can be derived from the same expansion by taking a series of partial derivatives. This requires more expansion terms to be kept since terms are lost in the process of taking partial derivatives. Thus to maintain accuracy for the doublet computation, more terms are required than if just evaluating for sources. The resulting Fast Multipole code should then parallelize better than classical panel methods due to the locality of data dependencies found in the Fast Multipole Method. Theoretically the parallelized code should execute in O(log N) time with O(N) processors, though this is not practical. Ongoing work includes implementing the parallel accelerated panel method, including methods to improve the load balancing of the problem by taking advantage of the known geometry of panels, and to encorporate sensitivity analysis into the algorithm.

Description: Asymptotic similarity states at large Reynolds numbers and small Rossby numbers in rotating homogeneous turbulence are investigated using the database obtained from large-eddy simulations of the incompressible Navier-Stokes equations. Previous work has shown that the turbulence kinetic energy and integral length scales are accurately described by simple scaling laws based on the low wavenumbers part of the three-dimensional energy spectrum. The primary interest of the present study is to search for spectrum similarity in the asymptotic state. Four independent energy spectra are defined. It is shown that rescaling of these energy spectra in the asymptotic regime will collapse three out of the four spectra. The spectrum which does not collapse is a function only of the vertical wavenumber and corresponds to two-component motions in the plane normal to the rotation axis. Detailed investigation of the cause of this anomalous behavior reveals the existence of a strong reverse cascade of energy from small-to-large scales of the two-dimensional, two-component motions. This feature of the rotating flow is presumably linked to the lack of a complete similarity state, though further study of this issue is required.

Description: This paper describes a path planning method for planetary rovers to search for paths on planetary surfaces. The planetary rover is required to travel safely over a long distance for many days over unfamiliar terrain. Hence it is very important how planetary rovers process sensory information in order to understand the planetary environment and to make decisions based on that information. As a new data structure for informational mapping, an extended elevation map (EEM) has been introduced, which includes the effect of the size of the rover. The proposed path planning can be conducted in such a way as if the rover were a point while the size of the rover is automatically taken into account. The validity of the proposed methods is verified by computer simulations.

Description: The long-term goal of the current research is the study of Large-Eddy Simulation (LES) as a tool for aeroacoustics. New algorithms and developments in computer hardware are making possible a new generation of tools for aeroacoustic predictions, which rely on the physics of the flow rather than empirical knowledge. LES, in conjunction with an acoustic analogy, holds the promise of predicting the statistics of noise radiated to the far-field of a turbulent flow. LES's predictive ability will be tested through extensive comparison of acoustic predictions based on a Direct Numerical Simulation (DNS) and LES of the same flow, as well as a priori testing of DNS results. The method presented here is aimed at allowing simulation of a turbulent flow field that is both simple and amenable to acoustic predictions. A free shear flow is homogeneous in both the streamwise and spanwise directions and which is statistically stationary will be simulated using equations based on the Navier-Stokes equations with a small number of added terms. Studying a free shear flow eliminates the need to consider flow-surface interactions as an acoustic source. The homogeneous directions and the flow's statistically stationary nature greatly simplify the application of an acoustic analogy.

Description: Detailed three-dimensional three-component phase averaged measurements of the spanwise and streamwise vorticity formation and evolution in acoustically forced plane free-shear flows have been obtained. For the first time, phase-averaged measurements of all three velocity components have been obtained in both a mixing layer and a wake on three-dimensional grids, yielding the spanwise and streamwise vorticity distributions without invoking Taylor's hypothesis. Initially, two-frequency forcing was used to phase-lock the roll-up and first pairing of the spanwise vortical structures in a plane mixing layer. The objective of this study was to measure the near-field vortical structure morphology in a mixing layer with 'natural' laminar initial boundary layers. For the second experiment the second and third subharmonics of the fundamental roll-up frequency were added to the previous two-frequency forcing in order to phase-lock the roll-up and first three pairings of the spanwise rollers in the mixing layer. The objective of this study was to determine the details of spanwise scale changes observed in previous time-averaged measurements and flow visualization of unforced mixing layers. For the final experiment, single-frequency forcing was used to phase-lock the Karman vortex street in a plane wake developing from nominally two-dimensional laminar initial boundary layers. The objective of this study was to compare measurements of the three-dimensional structure in a wake developing from 'natural' initial boundary layers to existing models of wake vortical structure.

Description: It was previously shown that the strong vorticity in isotropic turbulence is organized into tubular vortices ('worms') whose properties were characterized through the use of full numerical simulations at several Reynolds numbers. At the time most of the observations were kinematic, and several scaling laws were discovered for which there was no theoretical explanation. In the meantime, further analysis of the same fields yielded new information on the generation of the vortices, and it was realized that even if they had to be formed by stretching, they were at any given moment actually compressed at many points of their axes. This apparent contradiction was partially explained by postulating axial inertial waves induced by the nonuniformity of the vortex cores, which helped to 'spread' the axial strain and allowed the vortices to remain compact even if not uniformly stretched. The existence of such solutions was recently proved numerically. The present report discusses a set of new numerical simulations of isotropic turbulence, and a reanalysis of the old ones, in an effort to prove or disprove the presence of these waves in actual turbulent flows and to understand the dynamics, as opposed to the kinematics, of the vortices.

Description: A method for active fluid flow control based on control theory is discussed. Dynamic programming and fixed point successive approximations are used to accommodate the nonlinear control problem. The long-term goal of this project is to establish an effective method applicable to complex flows such as turbulence and jets. However, in this report, the method is applied to stochastic Burgers equation as an intermediate step towards this goal. Numerical results are compared with those obtained by gradient search methods.

Description: The objective of the present research is the development and application of efficient adaptive numerical algorithms for the study, via direct numerical simulations, of active vortex generators. We are using innovative computational schemes to investigate flows past complex configurations undergoing arbitrary motions. Some of the questions we try to answer are: Can and how may we control the dynamics of the wake? What is the importance of body shape and motion in the active control of the flow? What is the effect of three-dimensionality in laboratory experiments? We are interested not only in coupling our results to ongoing, related experimental work, but furthermore to develop an extensive database relating the above mechanisms to the vortical wake structures with the long-range objective of developing feedback control mechanisms. This technology is very important to aircraft, ship, automotive, and other industries that require predictive capability for fluid mechanical problems. The results would have an impact in high angle of attack aerodynamics and help design ways to improve the efficiency of ships and submarines (maneuverability, vortex induced vibration, and noise).

Description: A challenge for traditional turbulence modeling, based on the Reynolds averaged Navier-Stokes equations, remains the accurate prediction of 'mild', adverse pressure-gradient driven separation from a smooth surface. With this study we want to explore the capability of large-eddy simulation to predict the separation which occurs on the deflected wall of an asymmetric, plane diffuser with opening angle of 10 deg. The flow through the plane diffuser exhibits some additional interesting physical phenomena which make it a challenging test case. In addition to 'mild' separation about halfway down the deflected ramp, the flow is characterized by a small backflow zone with stalled fluid in the rear part of the expanding section. The turbulent flow entering the diffuser is subject to combined adverse and radial pressure gradients stemming from the convex curvature. Finally the flow recovers into a developed, turbulent channel flow in the outlet section. Obi et al. provide measurements of mean flow, Reynolds stresses, and pressure recovery, which were obtained by means of LDV in a wind tunnel. The objective of this study is to investigate whether LES with the standard dynamic model is able to accurately predict the flow in the one-sided diffuser and to explore the resolution requirements and associated costs.

Description: The dynamic subgrid-scale (SGS) model has given good results in the large-eddy simulation (LES) of homogeneous isotropic or shear flow, and in the LES of channel flow, using averaging in two or three homogeneous directions (the DA model). In order to simulate flows in general, complex geometries (with few or no homogeneous directions), the dynamic SGS model needs to be applied at a local level in a numerically stable way. Channel flow, which is inhomogeneous and wall-bounded flow in only one direction, provides a good initial test for local SGS models. Tests of the dynamic localization model were performed previously in channel flow using a pseudospectral code and good results were obtained. Numerical instability due to persistently negative eddy viscosity was avoided by either constraining the eddy viscosity to be positive or by limiting the time that eddy viscosities could remain negative by co-evolving the SGS kinetic energy (the DLk model). The DLk model, however, was too expensive to run in the pseudospectral code due to a large near-wall term in the auxiliary SGS kinetic energy (k) equation. One objective was then to implement the DLk model in a second-order central finite difference channel code, in which the auxiliary k equation could be integrated implicitly in time at great reduction in cost, and to assess its performance in comparison with the plane-averaged dynamic model or with no model at all, and with direct numerical simulation (DNS) and/or experimental data. Other local dynamic SGS models have been proposed recently, e.g., constrained dynamic models with random backscatter, and with eddy viscosity terms that are averaged in time over material path lines rather than in space. Another objective was to incorporate and test these models in channel flow.

Description: Two projects are described in this report. The first involves assessing turbulence models in separated flow. The second addresses the anomalous behavior of certain turbulence models in stagnation point flow. The primary motivation for developing turbulent transport models is to provide tools for computing non-equilibrium, or complex, turbulent flows. Simple flows can be analyzed using data correlations or algebraic eddy viscosities, but in more complicated flows such as a massively separated boundary layer, a more elaborate level of modeling is required. It is widely believed that at least a two-equation transport model is required in such cases. The transport equations determine the evolution of suitable velocity and time-scales of the turbulence. The present study included assessment of second-moment closures in several separated flows, including sharp edge separation; smooth wall, pressure driven separation; and unsteady vortex shedding. Flows with mean swirl are of interest for their role in enhancing mixing both by turbulent and mean motion. The swirl can have a stabilizing effect on the turbulence. An axi-symmetric extension to the INS-2D computer program was written adding the capability of computing swirling flow. High swirl can produce vortex breakdown on the centerline of the jet and it occurs in various combustors.

Description: A model for the effects of rotation on turbulence is proposed and tested. These effects which influence mainly the rate of turbulence decay are modeled in a modified turbulent energy dissipation rate equation that has explicit dependence on the mean rotation rate. An appropriate definition of the rotation rate derived from critical point theory and based on the invariants of the deformation tensor is proposed. The modeled dissipation rate equation is numerically well behaved and can be used in conjunction with any level of turbulence closure. The model is applied to the two-equation kappa-epsilon turbulence model and is used to compute separated flows in a backward-facing step and an axisymmetric swirling coaxial jets into a sudden expansion. In general, the rotation modified dissipation rate model shows some improvements over the standard kappa-epsilon model.

Description: The long-time, asymptotic state of rotating homogeneous shallow-water equations is investigated. Our analysis is based on long-time averaged rotating shallow-water equations describing interactions of large-scale, horizontal, two-dimensional motions with surface inertial-gravity waves field for a shallow, uniformly rotating fluid layer. These equations are obtained in two steps: first by introducing a Poincare/Kelvin linear propagator directly into classical shallow-water equations, then by averaging. The averaged equations describe interaction of wave fields with large-scale motions on time scales long compared to the time scale 1/f(sub o) introduced by rotation (f(sub o)/2-angular velocity of background rotation). The present analysis is similar to the one presented by Waleffe (1991) for 3D Euler equations in a rotating frame. However, since three-wave interactions in rotating shallow-water equations are forbidden, the final equations describing the asymptotic state are simplified considerably. Special emphasis is given to a new conservation law found in the asymptotic state and decoupling of the dynamics of the divergence free part of the velocity field. The possible rising of a decoupled dynamics in the asymptotic state is also investigated for homogeneous turbulence subjected to a background rotation. In our analysis we use long-time expansion, where the velocity field is decomposed into the 'slow manifold' part (the manifold which is unaffected by the linear 'rapid' effects of rotation or the inertial waves) and a formal 3D disturbance. We derive the physical space version of the long-time averaged equations and consider an invariant, basis-free derivation. This formulation can be used to generalize Waleffe's (1991) helical decomposition to viscous inhomogeneous flows (e.g. problems in cylindrical geometry with no-slip boundary conditions on the cylinder surface and homogeneous in the vertical direction).

Description: Understanding of turbulent free shear flows (wakes, jets, and mixing layers) is important, not only for scientific interest, but also because of their appearance in numerous practical applications. Turbulent wakes, in particular, have recently received increased attention by researchers at NASA Langley. The turbulent wake generated by a two-dimensional airfoil has been selected as the test-case for detailed high-resolution particle image velocimetry (PIV) experiments. This same wake has also been chosen to enhance NASA's turbulence modeling efforts. Over the past year, the author has completed several wake computations, while visiting NASA through the 1993 and 1994 ASEE summer programs, and also while on sabbatical leave during the 1993-94 academic year. These calculations have included two-equation (K-omega and K-epsilon) models, algebraic stress models (ASM), full Reynolds stress closure models, and direct numerical simulations (DNS). Recently, there has been mutually beneficial collaboration of the experimental and computational efforts. In fact, these projects have been chosen for joint presentation at the NASA Turbulence Peer Review, scheduled for September 1994. DNS calculations are presently underway for a turbulent wake at Re(sub theta) = 1000 and at a Mach number of 0.20. (Theta is the momentum thickness, which remains constant in the wake of a two dimensional body.) These calculations utilize a compressible DNS code written by M. M. Rai of NASA Ames, and modified for the wake by J. Cimbala. The code employs fifth-order accurate upwind-biased finite differencing for the convective terms, fourth-order accurate central differencing for the viscous terms, and an iterative-implicit time-integration scheme. The computational domain for these calculations starts at x/theta = 10, and extends to x/theta = 610. Fully developed turbulent wake profiles, obtained from experimental data from several wake generators, are supplied at the computational inlet, along with appropriate noise. After some adjustment period, the flow downstream of the inlet develops into a fully three-dimensional turbulent wake. Of particular interest in the present study is the far wake spreading rate and the self-similar mean and turbulence profiles. At the time of this writing, grid resolution studies are underway, and a code is being written to calculate turbulence statistics from these wake calculations; the statistics will be compared to those from the ongoing PIV wake measurements, those of previous experiments, and those predicted by the various turbulence models. These calculations will lead to significant long-term benefits for the turbulence modeling effort. In particular, quantities such as the pressure-strain correlation and the dissipation rate tensor can be easily calculated from the DNS results, whereas these quantities are nearly impossible to measure experimentally. Improvements to existing turbulence models (and development of new models) require knowledge about flow quantities such as these. Present turbulence models do a very good job at prediction of the shape of the mean velocity and Reynolds stress profiles in a turbulent wake, but significantly underpredict the magnitude of the stresses and the spreading rate of the wake. Thus, the turbulent wake is an ideal flow for turbulence modeling research. By careful comparison and analysis of each term in the modeled Reynolds stress equations, the DNS data can show where deficiencies in the models exist; improvements to the models can then be attempted.

Description: A new formulation of the dynamic subgrid-scale model is tested in which the error associated with the Germano identity is minimized over flow pathlines rather than over directions of statistical homogeneity. This procedure allows the application of the dynamic model with averaging to flows in complex geometries that do not possess homogeneous directions. The characteristic Lagrangian time scale over which the averaging is performed is chosen such that the model is purely dissipative, guaranteeing numerical stability when coupled with the Smagorinsky model. The formulation is tested successfully in forced and decaying isotropic turbulence and in fully developed and transitional channel flow. In homogeneous flows, the results are similar to those of the volume-averaged dynamic model, while in channel flow, the predictions are superior to those of the plane-averaged dynamic model. The relationship between the averaged terms in the model and vortical structures (worms) that appear in the LES is investigated. Computational overhead is kept small (about 10 percent above the CPU requirements of the volume or plane-averaged dynamic model) by using an approximate scheme to advance the Lagrangian tracking through first-order Euler time integration and linear interpolation in space.

Description: We present analysis and numerical experiments on the instability of streamwise vortices in 'minimal channel' flows and argue that this instability is a key feature in the observed intermittent cycle of formation, break-up, and re-formation of these structures. The base flow is a three-component, two-dimensional pair of counter-rotating rolls with axes aligned along the direction of the mean shear. While it is not a steady solution to the Navier-Stokes equations, we show numerically that this flow is unstable on a fast time scale to a secondary, three-dimensional Floquet mode. The growth of the secondary instability does not saturate in a new equilibrium, but continues until highly unstable local shear layers form and the entire flow breaks down into turbulence. Our analysis is motivated in part by the strong similarities between the intermittent turbulent cycle in minimal channel flows and one studied, both experimentally and in computations, in Couette-Taylor flow.

Description: The decay of a homogeneous turbulence generated by an axisymmetric distribution of random impulsive forces acting at the initial instant is considered. The impulsive forces may be either parallel or perpendicular to the symmetry axis. For impulsive forces which result in a kappa(exp 4) low wavenumber energy spectrum of the turbulence, it is determined that the flow approaches isotropy on all scales of motion at long-times provided the Reynolds number is large. However, for the type of impulsive forces originally proposed by Saffman in which a sq kappa low wavenumber energy spectrum is produced, the turbulence approaches isotropy only at the smallest scales and remains significantly anisotropic at the largest and energy-containing scales. Nevertheless, a similarity state of the flow field establishes itself asymptotically, in which the kinetic energy per unit mass of the turbulence decays as t(exp -6/5).

Description: The statistics of pressure and pressure gradient in stationary isotropic turbulence are measured within direct numerical simulations at low to moderate Reynolds numbers. It is found that the one-point pdf of the pressure is highly skewed and that the pdf of the pressure gradient is of stretched exponential form. The power spectrum of the pressure P(k) is found to be larger than the corresponding spectrum P(sub G)(k) computed from a Gaussian velocity field having the same energy spectrum as that of the DNS field. The ratio P(k)/P(sub G)(k), a measure of the pressure-field intermittence, grows with wavenumber and Reynolds number as -R(sub lambda)(exp 1/2)log(k/k(sub d)) for k less than k(sub d)/2 where k(sub d) is the Kolmogorov wavenumber. The Lagrangian correlations of pressure gradient and velocity are compared and the Lagrangian time scale of the pressure gradient is observed to be much shorter than that of the velocity.

Description: The summer assignment to study sheared turbulent flow was divided into three phases which were: (1) literature survey, (2) computational familiarization, and (3) pilot computational studies. The governing equations of fluid dynamics or Navier-Stokes equations describe the velocity, pressure, and density as functions of position and time. In principle, when combined with conservation equations for mass, energy, and thermodynamic state of the fluid a determinate system could be obtained. In practice the Navier-Stokes equations have not been solved due to the nonlinear nature and complexity of these equations. Consequently, the importance of experiments in gaining insight for understanding the physics of the problem has been an ongoing process. Reasonable computer simulations of the problem have occured as the computational speed and storage of computers has evolved. The importance of the microstructure of the turbulence dictates the need for high resolution grids in extracting solutions which contain the physical mechanisms which are essential to a successful simulation. The recognized breakthrough occurred as a result of the pioneering work of Orzag and Patterson in which the Navier-Stokes equations were solved numerically utilizing a time saving toggling technique between physical and wave space, known as a spectral method. An equally analytically unsolvable problem, containing the same quasi-chaotic nature as turbulence, is known as the three body problem which was studied computationally as a first step this summer. This study was followed by computations of a two dimensional (2D) free shear layer.

Description: The Center for Turbulent Research (CTR) post-processing effort aims to make turbulence simulations and data more readily and usefully available to the research and industrial communities. The Tensoral language, introduced in this document and currently existing in prototype form, is the foundation of this effort. Tensoral provides a convenient and powerful protocol to connect users who wish to analyze fluids databases with the authors who generate them. In this document we introduce Tensoral and its prototype implementation in the form of a user's guide. This guide focuses on use of Tensoral for post-processing turbulence databases. The corresponding document - the Tensoral 'author's guide' - which focuses on how authors can make databases available to users via the Tensoral system - is currently unwritten. Section 1 of this user's guide defines Tensoral's basic notions: we explain the class of problems at hand and how Tensoral abstracts them. Section 2 defines Tensoral syntax for mathematical expressions. Section 3 shows how these expressions make up Tensoral statements. Section 4 shows how Tensoral statements and expressions are embedded into other computer languages (such as C or Vectoral) to make Tensoral programs. We conclude with a complete example program.

Description: In the present work, we propose to conduct direct numerical simulations (DNS) of incompressible turbulent axisymmetric jets and wakes. The objectives of the study are to understand the fundamental behavior of axisymmetric jets and wakes, which are perhaps the most technologically relevant free shear flows (e.g. combuster injectors, propulsion jet). Among the data to be generated are various statistical quantities of importance in turbulence modeling, like the mean velocity, turbulent stresses, and all the terms in the Reynolds-stress balance equations. In addition, we will be interested in the evolution of large-scale structures that are common in free shear flow. The axisymmetric jet or wake is also a good problem in which to try the newly developed b-spline numerical method. Using b-splines as interpolating functions in the non-periodic direction offers many advantages. B-splines have local support, which leads to sparse matrices that can be efficiently stored and solved. Also, they offer spectral-like accuracy that are C(exp O-1) continuous, where O is the order of the spline used; this means that derivatives of the velocity such as the vorticity are smoothly and accurately represented. For purposes of validation against existing results, the present code will also be able to simulate internal flows (ones that require a no-slip boundary condition). Implementation of no-slip boundary condition is trivial in the context of the b-splines.

Description: During the last three years we have conducted high- and low-Reynolds-number experiments, including hot-wire measurements of the velocity fluctuations, in the test-section-ceiling boundary layer of the 80- by 120-foot Full-Scale Aerodynamics Facility at NASA Ames Research Center, to test the local-isotropy predictions of Kolmogorov's universal equilibrium theory. This hypothesis, which states that at sufficiently high Reynolds numbers the small-scale structures of turbulent motions are independent of large-scale structures and mean deformations, has been used in theoretical studies of turbulence and computational methods such as large-eddy simulation; however, its range of validity in shear flows has been a subject of controversy. The present experiments were planned to enhance our understanding of the local-isotropy hypothesis. Our experiments were divided into two sets. First, measurements were taken at different Reynolds numbers in a plane boundary layer, which is a 'simple' shear flow. Second, experiments were designed to address this question: will our criteria for the existence of local isotropy hold for 'complex' nonequilibrium flows in which extra rates of mean strain are added to the basic mean shear?

Description: previously a description was given of an active control scheme using wall transpiration that leads to a 15% reduction in surface skin friction beneath a turbulent boundary layer, according to direct numerical simulation. In this research brief further details of that scheme and its variants are given together with some suggestions as to how sensor/actuator arrays could be configured to reduce surface drag. The research which is summarized here was performed during the first half of 1994. This research is motivated by the need to understand better how the dynamics of near-wall turbulent flow can be modified so that skin friction is reduced. The reduction of turbulent skin friction is highly desirable in many engineering applications. Experiments and direct numerical simulations have led to an increased understanding of the cycle of turbulence production and transport in the boundary layer and raised awareness of the possibility of disrupting the process with a subsequent reduction in turbulent skin friction. The implementation of active feedback control in a computational setting is a viable approach for the investigation of the modifications to the flow physics that can be achieved. Bewley et al. and Hill describe how ideas from optimal control theory are employed to give 'sub-optimal' drag reduction schemes. The objectives of the work reported here is to investigate in greater detail the assumptions implicit within such schemes and their limitations. It is also our objective to describe how an array of sensors and actuators could be arranged and interconnected to form a 'smart' surface which has low skin friction.

Description: Conventional vortex generators as found on many civil aircrafts are mainly for off-design conditions - e.g. suppression of separation or loss of aileron power when the Mach number accidentally rises above the design (cruise) value. In normal conditions they perform no useful function and exert a significant drag penalty. Recently there have been advances in new designs for passive vortex generators and boundary layer control. While traditionally the generators heights were of the order of the boundary layer thickness (delta), recent advances have been made where generators of the order of delta/4 have been shown to be effective. The advancement of MIcro-Electro-Mechanical (MEM) devices has prompted several efforts in exploring the possibility of using such devices in turbulence control. These new devices offer the possibility of boundary layer manipulation through the production of vortices, momentum jets, or other features in the flow. However, the energy output of each device is low in general, but they can be used in large numbers. Therefore, the possibility of moving from passive vortex generators to active (on-demand) devices becomes of interest. Replacement of fixed rectangular or delta-wing generators by devices that could be activated when needed would produce substantial economies. Our proposed application is not strictly 'active' control: the vortex generators would simply be switched on, all together, when needed (e.g. when the aircraft Mach number exceeded a certain limit). To this extent our scheme is simpler; however, to promote mixing and suppress separation we desire to deposit longitudinal vortices into the outer layer of the boundary layer as in conventional vortex generators. This requires a larger device although an alternative might be an array of smaller devices, for example, a longitudinal row with phase differences in the modulation signals so that the periodic vortices join up. The vortex pair with common flow up has the advantage that it will naturally drift away from the surface, but the disadvantage is that the net vorticity is zero so that the pair is eventually obliterated by turbulent mixing, rather than simply being diffused as in the case of a single vortex. It should be possible to devise alternative shapes of cavity wall so that the jet emerges obliquely and produces net longitudinal vorticity.

Description: This work investigates the turbulent constitutive relation when turbulence is subjected to solid body rotation. Laws regarding spectra and asymptotic decay of rotating homogeneous turbulence were confirmed through large-eddy simulation (LES) computations. Rotating turbulent flows exist in many industrial, geophysical, and astrophysical applications. From Lagrangian analysis a relation between turbulent stress and strain in rotating homogeneous turbulence was inferred. This relation was used to derive the spectral energy flux and, ultimately, the energy spectrum form. If the rotation wavenumber k(sub Omega) lies in the inertial subrange, then for wavenumbers less than k(sub Omega) the turbulence motions are affected by rotation and the energy spectrum slope is modified. Energy decay laws inferred in other reports and the present results suggest a modification of the epsilon model equation and eddy viscosity in k-epsilon models.

Description: Large-eddy simulation results are used to investigate the development of anisotropies and the possible transition towards a quasi two-dimensional state in rotating turbulence at high Reynolds number. The present study demonstrates the existence of two transitions that are identified by two Rossby numbers. The first transition marks the onset of anisotropic effects and corresponds to a macro Rossby number Ro(sup L) (based on a longitudinal integral length scale) near unity. A second transition can be defined in terms of a lower bound of micro-Rossby number Ro(sup w) also near unity (defined in this work as the ratio of the rms fluctuating vorticity to background vorticity) and corresponds to a continued development of anisotropy but with an increasing emergence of those indicators based on the pure two-dimensional component of the flow, e.g., integral length scales measured along the rotation axis. Investigation of the vorticity structure shows that the second transition is also characterized by an increasing tendency for alignment between the fluctuating vorticity vector and the basic angular velocity vector with a preference for corotative vorticity.

Description: Direct numerical simulations are performed for homogeneous turbulence with a mean flow having elliptic streamlines. This flow combines the effects of rotation and strain on the turbulence. Qualitative comparisons are made with linear theory for cases with high Rossby number. The nonlinear transfer process is monitored using a generalized skewness. In general, rotation turns off the nonlinear cascade; however, for moderate ellipticities and rotation rates the nonlinear cascade is turned off and then reestablished. Turbulence statistics of interest in turbulence modeling are calculated, including full Reynolds stress budgets.

Description: NASA/Marshall Space Flight Center (MSFC) has an ongoing effort to transfer to industry the technologies developed at MSFC for rocket propulsion systems. The Technology Utilization (TU) Office at MSFC promotes these efforts and accepts requests for assistance from industry. One such solicitation involves a request from North American Marine Jet, Inc. (NAMJ) for assistance in the design of a water-jet-drive system to fill a gap in NAMJ's product line. NAMJ provided MSFC with a baseline axial flow impeller design as well as the relevant working parameters (rpm, flow rate, etc.). This baseline design was analyzed using CFD, and significant deficiencies identified. Four additional analyses were performed involving MSFC changes to the geometric and operational parameters of the baseline case. Subsequently, the impeller was redesigned by NAMJ and analyzed by MSFC. This new configuration performs significantly better than the baseline design. Similar cooperative activities are planned for the design of the jet-drive inlet.

Description: The development of fractures at regular length scales is a widespread feature of Venusian tectonics. Models of lithospheric deformation under extension based on non-Newtonian viscous flow and brittle-plastic flow develop localized failure at preferred wavelengths that depend on lithospheric thickness and stratification. The characteristic wavelengths seen in rift zones and tessera can therefore provide constraints on crustal and thermal structure. Analytic solutions were obtained for growth rates in infinitesimal perturbations imposed on a one-dimensional, layered rheology. Brittle layers were approximated by perfectly-plastic, uniform strength, overlying ductile layers exhibiting thermally-activated power-law creep. This study investigates the formation of faults under finite amounts of extension, employing a finite-element approach. Our model incorporates non-linear viscous rheology and a Coulomb failure envelope. An initial perturbation in crustal thickness gives rise to necking instabilities. A small amount of velocity weakening serves to localize deformation into planar regions of high strain rate. Such planes are analogous to normal faults seen in terrestrial rift zones. These 'faults' evolve to low angle under finite extension. Fault spacing, orientation and location, and the depth to the brittle-ductile transition, depend in a complex way on lateral variations in crustal thickness. In general, we find that multiple wavelengths of deformation can arise from the interaction of crustal and mantle lithosphere.

Description: A coupled problem of diffusion and condensation is solved for the H2SO4-H2O system in Venus' cloud layer. The position of the lower cloud boundary, profiles of the H2O and H2SO vapor mixing ratios, and of the H2O/H2SO4 ratio of sulfuric acid aerosol and its flux are calculated as functions of the column photochemical production rate of sulfuric acid. Variations of the lower cloud boundary are considered. Our basic model, which is constrained to yield f(sub H2O)(30km) = 30 ppm, predicts the position of the lower cloud boundary at 48.4 km coinciding with the mean Pioneer Venus value, the peak H2SO4 mixing ratio of 5.4 ppm, and the H2SO4 production rate phi(sub H2SO4) = 2.2 x 10(exp 12) cm(exp -2)s(exp -1).

Description: Based on the conservation of chemical elements in chemical reactions, a rule is proved that the number of boundary conditions given by densities and/or non-zero velocities should not be less than the number of chemical elements in the system, and the components given by densities and velocities should include all elements in the system. Applications of this rule to Mars are considered. It is proved that a problem of CO2-H2O chemistry in the lower and middle atmosphere of Mars, say, in the range of 0-80 km does not have an unique solution, if only CO2 and H2O densities are given at the lower boundary, while all other boundary conditions are fluxes. Two models of this type are discussed. These models fit the same boundary conditions, are balanced with a relative uncertainty of 10(exp -4) for H2, and predict the O2, CO, and H2 mixing ratios which differ by order of magnitude. One more species density, e.g. that of O2, should be specified at the boundary to obtain the unique solution. The situation is better if the upper boundary is extended to the exobase where thermal escape velocities of H and H2 can be specified. However, in this case, either oxygen nonthermal escape rate or the O2 density at the surface should be given as the boundary condition. Two models of Mars' photochemistry, with and without nitrogen chemistry, are considered. The oxygen nonthermal escape rate of 1.2 x 10(exp 8) cm(exp -2) s(exp -1) is given at 240 km and is balanced with the total hydrogen escape rate within uncertainty of 1 percent for both models. Both models fit the measured O2 and CO mixing ratios, the O3 line absorption at 9.6 microns, and the O2 1.27 microns dayglow within the uncertainties of the measured values; although, the model without nitrogen chemistry fits better.

Description: Mid-infrared spectroscopic observations of the surface of Mercury are reported for the wavelength range 7.3 to 13.5 microns. The observed spectral radiance emanated from equatorial and low latitude regions between 110-130 deg Mercurian longitude. The area is primarily an intercrater plain. The spectra show distinct and recognizable features, the principal Christiansen emission peak being the most prominent. The Christiansen feature strongly suggests the presence of plagioclase (Ca,Na)(Al,Si)AlSi2O8, (in particular labradorite: Ab(50) - Ab(30)). In addition we have studied the effects of thermal gradients to gain insight into the effects of thermal conditions on the spectral radiance of rock samples. This simulates the thermophysical effects as the rotating surface of Mercury is alternately heated and cooled. The spectral features of the samples are retained; however, the relative and absolute amplitudes vary as illustrated by laboratory reflectance and emittance spectra from quartzite.

Description: Given the absence of ground truth information on seismic structure, heat flow, and rock strength, or short wavelength gravity or magnetic data for Venus, information on the thermal, mechanical and compositional nature of the shallow interior must be obtained by indirect methods. Using pre-Magellan data, theoretical models constrained by the depths of impact craters and the length scales of tectonic features yielded estimates on the thickness of Venus' brittle-elastic lithosphere and the allowable range of crustal thickness and surface thermal gradient. The purpose of this study is to revisit the question of the shallow structure of Venus based on Magellan observations of the surface and recent experiments that address Venus' crustal rheology.

Description: An outstanding question relevant to understanding the tectonics of Venus is the mechanism of formation of fold and thrust belts, such as the mountain belts that surround Lakshmi Planum in western Ishtar Terra. These structures are typically long (hundreds of km) and narrow (many tens of km), and are often located at the margins of relatively high (km-scale) topographic rises. Previous studies have attempted to explain fold and thrust belts in various areas of Venus in the context of viscous and brittle wedge theory. However, while wedge theory can explain the change in elevation from the rise to the adjacent lowland, it fails to account for a fundamental aspect of the deformation, i.e., the topographic high at the edge of the rise. In this study we quantitatively explore the hypothesis that fold and thrust belt morphology on Venus can alternatively be explained by horizontal shortening of a lithosphere that is laterally heterogeneous, due either to a change in thickness of the lithosphere or the crust. Lateral heterogeneities in lithosphere structure may arise in response to thermal thinning or extensive faulting, while variations in crustal thickness may arise due to either spatially variable melting of mantle material or by horizontal shortening of the crust. In a variable thickness lithosphere or crust that is horizontally shortened, deformation will tend to localize in the vicinity of thickness heterogeneity, resulting in a higher component of dynamic topography there as compared to elsewhere in the shortening lithosphere. This mechanism may thus provide a simple explanation for the topographic high at the edge of the rise.

Description: The Clementine spacecraft is to be launched into Earth orbit in late January for subsequent insertion into lunar orbit in late February, 1994. There, its primary mission is to produce -- over a period of about two months -- a new photographic map of the entire surface of the Moon; this will be done, in a variety of wavelengths and spatial resolutions, in a manner greatly superior to that previously accomplished for the whole Moon. It will then go on to fly by and photograph the asteroid Geographos. A secondary goal that has been accepted for this mission is to take a series of photographs designed to capture images of, and determine the brightness and extent of, the Lunar Horizon Glow (LHG). One form of LHG is caused by the solar stimulation of emission from Na and K atoms in the lunar exosphere. The scale height of this exosphere is of the order of 100 km. There are also brighter LHG components, with much smaller scale heights, that appear to be caused by scattered off of an exospheric lunar dust cloud.

Description: The Martian regolith is the most substantial volatile reservoir on the planet; it holds CO2 as adsorbate, and can exchange that CO2 with the atmosphere-cap system over timescales of 10(exp 5) to 10(exp 6) years. The climatic response to insolation changes caused by obliquity and eccentricity variations depends in part on the total reservoir of adsorbed CO2. Previous estimates of the adsorbate inventory have been made by measuring the adsorptive behavior of one or more Mars-analyog materials, and deriving an empirical equation that described that adsorption as a function of the partial pressure of CO2 and the temperature of the regolith. The current CO2 inventory is that which satisfies adsorptive equilibrium, observed atmospheric pressure, and no permanent CO2 caps. There is laboratory evidence that H2O poisons the CO2 adsorptive capacity of most materials. No consideration of CO2 - H2O co-adsorption was given in previous estimates of the Martian CO2 inventory, although H2O is present in the vapor phase, and so as adsorbate, throughout the regolith.

Description: The rate at which the wind can redistribute sedimentary material is an important part of any planet's sedimentologic cycle, particularly for Mars, where the competing effects of other gradational processes are less than on Earth. The aeolian drift potential (DP) is a measure of the amount of material capable of being moved through a unit length by the wind for a given period of time. DP is a useful measure of the potential redistribution rate of windblown material on regional scales. The Martian aeolian DP was calculated from laboratory studies of sand movement conducted at Martian atmospheric densities and from surface stress, temperature, and pressure values for that region as determined from the Mars General (Atmospheric) Circulation Model (GCM) developed at the NASA/Ames Research Center. In our simulations for Mars, DP changes in both magnitude (as expected) and direction if the saltation threshold is altered.

Description: Gravitational wave astronomy continues to be one of the exploration concepts under consideration in NASA's strategy for conducting physics and astrophysics from the lunar surface. As with other proposals for new concepts in science and astronomy from the Moon, this one has a number of very interesting features which need to be developed further in order to assess them adequately. The possibility of robotic deployment of a gravitational wave antenna on the Moon in a triangular configuration and the question of closure on the third interferometer leg are discussed here.

Description: With Hapke scattering theory and absorption coefficients derived from our laboratory measurements of solid N2 we have modeled the spectrum of Triton. By comparing a Hapke scattering model to the measured spectrum from Triton, we determined the temperature of the N2 on the satellite's surface to be 38 (+2, -1) K which is in accord with the measurements of Voyager 2. Applying this technique to Pluto we find that the temperature of N2 on that body is 40 +/- 2 K. Other aspects of this investigation are discussed.

Description: The two-dimensional chemical mapping analysis (CMA) techniques of EPMA and XRF were applied to a new polished thin section (PTS) of EET84302,28, Acapulco and a 5x3 cm slice of Caddo County to find heterogeneous regional distribution of low temperature fractions in the lodranite-acapulcoite groups and silicate inclusions in the IAB irons. A region richer in metal-plagioclase was found in EET84302,28 and Caddo County. The mineralogy of EET84302,28 is not much different from coarse-grained, metal-rich acapulcoite-like mineral assemblage in EET84302,19, which has chromite-orthopyroxene segregation. Nearly uniform Mg/Fe ratios of silicates modified by reduction at regional oxygen fugacity and large difference in modal abundances of minerals in this meteorite group can be explained by regional concentration of materials when the source materials were partly melted.

Description: We propose a new center (35.5 W, 32.5 N) and ring assignment for the original Chryse impact basin based on photogeologic mapping and re-examination of the published geology. Noachian features in the Chryse Planitia area are the best indicators of the original ancient multiringed impact structure. While other workers have centered the Chryse impact on the topographic low associated with Hesperian volcanic and fluvial deposits, we suggest that the center of the original Noachian-age excavation cavity was located 800 km farther NE, and the basin topography was significantly modified over time.

Description: Our present knowledge of the topography of Mars is completely inadequate for addressing a wide range of geophysical, geological, and atmospheric problems. The data acquired to date by several techniques has not provided us with reliable and consistent values for even the equatorial and polar radii, and large uncertainties exist in the altitudes of many of the major volcanic constructs. While much of this can be blamed on a lack of reliable data, we feel that much more could be done to improve the analysis of the present data in a consistent global system that would necessarily involve the re-analysis of early spacecraft and Earth-based data in conjunction with more recent models of Mars' gravity field.

Description: We have estimated the mass of Phobos, Deimos, and Mars using the Viking Orbiter and Mariner 9 tracking data. We divided the data into 282 arcs and sorted the data by periapse height, by inclination, and by satellite. The data were processed with the GEODYN/SOLVE orbit determination programs, which have previously been used to analyze planetary tracking data. The a priori Mars gravity field applied in this study was the 50th degree and order GMM-1 (Goddard Mars Model-1) model. The subsets of data were carefully edited to remove any arcs with close encounters of less than 500 km with either Phobos or Deimos. Whereas previous investigators have used close flybys (less than 500 km) to estimate the satellite masses, we have attempted to estimate the masses of Phobos and Deimos from multiday arcs which only included more distant encounters. The subsets of data were further edited to eliminate spurious data near solar conjunction (Nov.-Dec. 1976 and January 1979). In addition, the Viking-1 data from Oct. through Dec. 1978 were also excluded because of the low periapse altitude (as low as 232 km) and thus high sensitivity to atmospheric drag.

Description: Apollo 16 double drive tube 68002/68001 has been opened, and information on its contents and structure is now available. The sample was taken in Apr. 1972 at Station 8, a gently undulating surface of South Ray Crater ejecta material. The tube was driven into the surface about 2 m from the edge of a 10-15 m diameter crater in an area where small (less than 0.5 m diameter) craters are common. Fragments of the less than 1 cm diameter size are most abundant in this area. The top section, 68002, and the lower section, 68001, have a combined length of about 61 cm.

Description: The regolith breccias from the Luna 24 core were analyzed. The Luna 24 regolith is a mixture of fine and coarse grain materials. The comparable analysis of the grain size distributions, the modal and chemical compositions of the breccias, and the regolith from the same levels show that the friable slightly litificated breccia with a friable fine grain matrix may be a source of fine grain material of the Luna 24 present day regolith.

Description: The distribution of ancient massifs and old cratered terrain in the southern Daedalia region indicate the presence of at least two and probably three impact basins of large size. One of these is located near where Craddock et al. placed their center for the Daedalia Basin, but it has very different ring diameters. These basins have rings exceeding 1000 km diameter and overlap significantly with centers separated by 500 to 600 km at nearly identical latitudes of -26 to -29 deg. The smaller westernmost basin appears slightly better preserved, but there is little evidence for obvious superposition that might imply a temporal sequence. Recognizing the improbability of random impacts producing aligned, nearly contemporaneous features, we suggest these basins may have resulted from clustered impactors.

Description: Although present Martian surface conditions appear unfavorable for life as we know it, there is compelling geological evidence that the climate of early Mars was much more Earth-like, with a denser atmosphere and abundant surface water. The fact that life developed on the Earth within the first billion years of its history makes it quite plausible that life may have also developed on Mars. If life did develop on Mars, it is likely to have left behind a fossil record. This has led to the development of a new subdiscipline of paleontology, herein termed 'exopaleontology', which deals with the exploration for fossils on other planets. The most important factor enhancing microbial fossilization is the rapid entombment of microorganisms by fine-grained, stable mineral phases, such as silica, phosphate, or carbonate. The oldest body fossils on Earth are preserved in this way, occurring as permineralized cells in fine-grained siliceous sediments (cherts) associated with ancient volcanic terranes in Australia and South Africa. Modern terrestrial environments where minerals may precipitate in the presence of microorganisms include subaerial thermal springs and shallow hydrothermal systems, sub-lacustrine springs and evaporitic alkaline lakes, zones of mineralization within soils where 'hardpans' (e.g. calcretes, silcretes) form, and high latitude frozen soils or ground ice.

Description: Geologic mapping of the fretted terrain of the Nilosyrtis Mensae region of Mars has revealed geomorphic evidence that the breakup of the plateau units to the south of Nilosyrtis occurred well before the plains units to the north were emplaced in the late Hesperian time. The plains units were deposited against the fretted terrain which has undergone some modification by mass wasting but not significant backwasting. The morphology observed at the contact between plains and the fretted terrain is consistent with that expected where the edge of a pile of sedimentary debris has undergone mass wasting and other erosion.

Description: Many ancient, highly degraded large impact basins exist on Mars. In many cases these basins overlap or are overlapped by more easily observed, presumably younger, impact basins. While impact basin overlap is becoming more recognized, the effects of such overlap have only occassionally been described. Such effects will depend on a variety of factors including the absolute and relative size of the basins, the degree of overlap, the state of the lithosphere and its thermal gradient at the time of impact, and the time between impacts. There now exists enough evidence for overlapping basins of different sizes that some of these can be discussed. This paper highlights some examples of the obvious effects of basin overlap.

Description: In efforts to understand regolith evolution on airless bodies, increasing attention is now being payed to those processes and events that alter or 'weather' the surfaces of regolith grains. This attention has developed partly out of the ongoing need to optimize models of planetary reflectance spectra and the growing recognition that diverse types of grain coatings and surface alterations occur which can strongly influence mineral reflectance properties. In addition to their implications for optical properties, surface features on regolith grains have provided useful clues to the basic thermal, chemical, and radiation history of regoliths.

Description: Enlargements of Lunar-Orbiter photography were used in conjunction with a digitizing tablet to collect the locations and dimensions of blocks surrounding the Surveyor 1, 3, 6, and 7 landing sites. Data were reduced to the location and the major axis of the visible portion of each block. Shadows sometimes made it difficult to assess whether the visible major axis corresponded with the actual principal dimension. These data were then correlated with the locations of major craters in the study areas, thus subdividing the data set into blocks obviously associated with craters and those in intercrater areas. A block was arbitrarily defined to be associated with a crater when its location was within 1.1 crater radii of the crater's center. Since this study was commissioned for the ultimate purpose of determining hazards to landing spacecraft, such a definition was deemed appropriate in defining block-related hazards associated with craters. Size distributions of smaller fragments as determined from Surveyor photography were obtained as measurements from graphical data. Basic comparisons were performed through use of cumulative frequency distributions identical to those applied to studies of crater-count data.

Description: According to exobiology site selection criteria for Mars, the search for potential extinct/extant water dependent life should focus on sites were water flowed and ponded. The Ma'adim Vallis/Gusev crater system is of high priority for exobiology research, because it appears to have involved long term flooding, different periods and rates of sedimentation, and probable episodic ponding. The topics covered include the following: evidence of nonuniform fluvial processes and early overflooding of the plateau and ponding.

Description: The rocks and soils of the Moon provide raw materials essential to the successful establishment of a lunar base. Efficient exploitation of these resources requires accurate characterization of mineral abundances, sizes/shapes, and association of 'ore' and 'gangue' phases, as well as the technology to generate high-yield/high-grade feedstocks. Only recently have x-ray mapping and digital imaging techniques been applied to lunar resource evaluation. The topics covered include inherent differences between lunar basalts and soils and quantitative comparison of rock-derived and soil-derived ilmenite concentrates. It is concluded that x-ray digital-imaging characterization of lunar raw materials provides a quantitative comparison that is unattainable by traditional petrographic techniques. These data are necessary for accurately determining mineral distributions of soil and crushed rock material. Application of these techniques will provide an important link to choosing the best raw material for mineral beneficiation.

Description: The climate of Venus is to a large degree controlled by the radiative properties of its massive atmosphere. In addition, outgassing due to volcanic activity, exospheric escape processes, and surface/atmosphere interactions may all be important in moderating the abundances of atmospheric CO2 and other volatiles. We have developed an evolutionary climate model for Venus using a systems approach that emphasizes feedbacks between elements in the climate system. Modules for atmospheric radiative transfer, surface/atmosphere interactions, tropospheric chemistry, and exospheric escape processes have so far been developed. Climate feedback loops result from interconnections between modules, in the form of the environmental parameters pressure, temperature, and atmospheric mixing ratios. The radiative transfer module has been implemented by using Rosseland mean opacities in a one dimensional grey radiative-convective model. The model has been solved for the static (time independent) case to determine climate equilibrium points. The dynamics of the model have also been explored by employing reaction/diffusion kinetics for possible surface atmosphere heterogeneous reactions over geologic timescales. It was found that under current conditions, the model predicts that the climate of Venus is at or near an unstable equilibrium point. The effects of constant rate volcanism and corresponding exsolution of volatiles on the stability of the climate model were also explored.

Description: Ureilites Novo Urei, Havero, and Kenna show strong evidence of one or more Ar-40 degassing events in the time period of 3.3-4.1 Ga ago. These ages may be compared to current interpretations of ureilite chronology. These include the suggestion of metasomatic activity on the parent body 3.7 Ga ago that reset some Sm-Nd ages and the suggestion that ureilites have experienced terrestrial contamination of several trace elements (including Pb and LREE), which makes suspect ages younger than approximately 4.5 Ga. Because the K-Ar chronometer can be sensitive to metamorphic events, we made Ar-39-Ar-40 determinations on bulk samples (0.12-0.14 g each) of four ureilites. The Ar-39-Ar-40 age spectra and K/Ca ratios as a function of cumulative Ar release from stepwise temperature extractions for the four ureilites analyzed are shown. Because Ar-39-Ar-40 ages shown by low and high temperature extractions may be suspect, we examined the intermediate temperature extractions. Although interpretation of these spectra is obviously uncertain, we believe that the most recent times of Ar degassing can be roughly inferred. These times are approximately 3.3 Ga for Havero, 3.3-3.7 Ga for Novo Urei, and approximately 4.1 Ga for Kenna, for which Ar degassing may not have been complete. The indication of Ar-39-Ar-40 degassing ages of 3.3-4.1 Ga for three ureilites that also contain an enhanced LREE component and (excepting Havero) produce a 3.74 Ga Sm-Nd age, suggests that both chronometers may have responded to the same parent body event. On the other hand, it is also possible that the Ar data reflect one or more separate events that did not strongly affect the Sm-Nd system, a situation that commonly occurs in eucrites. Thus the existence of reset Ar ages does not require similarly reset Sm-Nd ages.

Description: All types of coarse-grained components in CM chondrites are surrounded by fine-grained dust coatings, but the origin of these rims is not yet clear. Although a strictly nebular origin seems likely for rims in the relatively unaltered type 3 chondrites, the rims in CM chondrites are dominated by secondary alteration phases. It has been argued that either the coarse-grained cores accreted altered rim materials while still in the nebula or that alteration of primary rim phases occurred on the CM parent body. To constrain the origin of alteration phases in rim material, we have analyzed the textures and mineral associations from 10 CM chondritic falls by optical and scanning electron microscopy. Our results indicate that the secondary phases in CM chondritic rims were produced by parent body fluid-rock interactions which redefined some primary rim textures and may have produced, in some cases, both coarse-grained components and the rims that surround them. Textural features demonstrate the interactive exchange of alteration fluids between rims, matrix, and chondrules on the CM parent body. For example, most matrix-rim contacts are gradational, suggesting the synchronous alteration of both components. Several observations suggest the possibility of in situ rim production. For example, tochilinite and phyllosilicates commonly form rims around matrix carbonates, which are generally believed to have precipitated from alteration fluids on the CM parent body. This suggests that the rims surrounding matrix carbonates may also have been produced by alteration processes. Partially replaced chondrule olivines bear a striking resemblance to many rimmed olivines in the matrix which suggests, by analogy, that site-specific precipitation of S-bearing phases may also be responsible for the occurrence of many tochilinite-rich rims around isolated matrix olivines. Non-silicate rims precipitate around olivines of any composition, but the process is most effective for fayalitic olivines. Most of the remaining olivines in CM chondrites are relatively Mg-rich, which suggests that the precipitation of S-bearing rims on olivines may not have been an important process in the aqueous alteration of CM chondrites. We conclude that: (1) precursor rim materials in CM chondrites were subjected to pervasive aqueous alteration on the CM parent body; and (2) textures and mineral associations observed in CM chondrites also suggest the possibility of in situ rim production.

Description: We have used digitized scanning electron micrographs and computer image analysis programs to measure track densities in lunar soil grains. Tracks were formed by highly ionizing solar energetic particles and cosmic rays. Back-scattered electron images produced suitable high contrast images for analysis. We used computer counting and measurement of area to obtain track densities. We found an excellent correlation with manual measurements for track densities below 1x10(exp 8) cm(exp -2). For track densities between 1x10(exp 8) to 1x10(exp 9) cm(exp -2) we found that a regression formula using the percentage area covered by tracks gave good agreement with manual measurements. Measurement of tract densities in lunar samples has been a very rewarding technique for measuring exposure ages and soil maturation processes. We have shown that we can reliably measure track densities in lunar grains using image analysis techniques. Automating track counting may allow application of this technique to important problems in regolith dynamics including the ratio of radiation exposure to reworking in various surface and core samples and in regolith breccias.

Description: The mineralogy of a planetary surface can be used to identify the provenance of soil or sediment and reveal the volcanic, metamorphic and/or sedimentological history of a particular region. We have discussed elsewhere the applications and the instrument design of possible X-ray diffraction and X-ray fluorescence (XRD/XRF) devices for the mineralogical characterization of planetary surfaces. In this abstract we evaluate some aspects of sample-detector geometry and sample collection strategies.

Description: 1994 marks the 25th anniversary of the Mariner 6 and 7 flyby missions to Mars. Despite its age, the Mariner 6,7 Infrared Spectrometer (IRS) data are a unique set of measurements that can provide important information about the Martian surface, atmospheric, and atmospheric aerosol composition. For certain mid-IR wavelengths, the IRS spectra are the only such spacecraft data obtained for Mars. At other wavelengths, IRS measured surface regions different from those measured by Mariner 9 or Phobos 2 and under different dust opacity conditions. We are interested in examining the IRS reflectance data in the 1.8 to 3.0 micron region because there are numerous diagnostic absorption features at these wavelengths that could be indicative of hydrated silicate minerals or of carbonate- or sulfate-bearing minerals. Groundbased telescopic data and recent Phobos ISM measurements have provided controversial and somewhat contradictory evidence for the existence of mineralogic absorption features at these wavelengths. Our goal is to determine whether any such features can be seen in the IRS data and to use their presence or absence to re-assess the quality and interpretations of previous telescopic and spacecraft measurements.

Description: We obtained new high resolution reflectance spectra of Mars during the 1993 opposition from Mauna Kea Observatory using the UKIRT CGS4 spectrometer. Fifty spectra of 1600-2000 km surface regions and a number of standard star spectra were obtained in the 2.04 to 2.44 micron wavelength region on 4 February 1993 UT. Near-simultaneous observations of bright standard stars were used to perform terrestrial atmospheric corrections and an absolute flux calibration. Using the known magnitude of the stars and assuming blackbody continuum behavior, the flux from Mars could be derived. A radiative transfer model and the HITRAN spectral line data base were used to compute atmospheric transmission spectra for Mars and the Earth in order to simulate the contributions of these atmospheres to our observed data. Also, we examined the ATMOS solar spectrum in the near-IR to try to identify absorption features in the spectrum of the Sun that could be misinterpreted as Mars features. Eleven absorption features were detected in our Mars spectra. Our data provide no conclusive identification of the mineralogy responsible for the absorption features we detected. However, examination of terrestrial spectral libraries and previous high spectral resolution mineral studies indicates that the most likely origin of these features is either CO3(sup 2-), HCO3(-), or HSO4(-) anions in framework silicates or possibly (Fe, Mg)-OH bonds in sheet silicates.

Description: The region of the moon near the border between Mare Serenitatis and Mare Tranquillitatis is one of the most geologically and compositionally complex areas of the nearside. The geologic history of this region has been shaped by impacts of widely-varying spatial scale and temporal occurrence, by volcanism of variable style and composition with time, and by limited tectonism. We have been studying this region as part of a larger multi remote sensing technique effort to understand the composition, morphology, geology, and stratigraphy of the moon at spatial scales of 2 km or less. The effort has been aided by the proximity of this area to the Apollo 11, 15, and 17 landing sites and by the occurrence of one of the primary lunar spectroscopic 'standard areas' within our scene (MS2). Here, some of the findings from the multispectral imaging and spectroscopy part of this effort are reported.

Description: The 0.4 to 1.0 micron spectrum of Mars is dominated by a steep red, relatively featureless spectral slope. Earlier lower spectral observations interpreted the red color and the lack of absorption features in the spectra as evidence of poorly crystalline ferric oxide minerals. More recent higher spectral resolution observations and reinterpretations of older data sets have revealed measureable spectral structure, however. For example, absorption features near 0.65 and 0.86 micron were detected and spatially mapped in data obtained during the 1988 opposition. These absorptions were interpreted as evidence for crystalline hematite on Mars, occuring as an accessory phase in abundances of 3 to 6 percent in the soil. We are attempting to verify the existence of these subtle crystalline Fe(3+) absorption features and to map their spatial distribution in regions of the planet not imaged in 1988. During the 1990 opposition, we obtained imaging spectroscopic data of Mars from the University of Hawaii 2.24 m telescope at Mauna Kea Observatory. The data were obtained with the Wide Field Grism Spectrograph (WFGS), which uses an 800 x 800 CCD and a transmission grating ruled on a prism. We used a grating blazed at 4800 A in first order to obtain data from 0.50 to 0.94 micron at a spectral resolution of R = 200 to 350. The moon/Mars slit design used had projected dimensions of 0.29 x 153 inches, allowing for high spectral resolution and adequate cross-slit spatial sampling of the Martian disk.

Description: Previous workers on single agglutinates have variously interpreted the composition of agglutinitic glass to represent impact melts of (1) bulk soil, (2) mixed components in finer sizes, and (3) microtargets. Separately, Papike has argued in favor of fusion of the finest fraction of bulk soils. Thirty-four single agglutinates were hand-picked from the mature Apollo 16 soil 61181 (I(sub s)/FeO = 82) and the FMR and chemical composition (INAA for Fe, Sc, Sm, Co, Ni, and Cr) of each agglutinate particle were measured. Thirteen of these single agglutinates were selected for electron beam microanalysis and imaging. Less than 1 micron spots were analyzed (for Na, Mg, Al, Si, P, S, K, Ca, Ti, Cr, Mn, Fe, Ni, and Ba) on pure glassy areas (approximately ten in each particle) selected on the basis of optical and BSE images (avoiding all clasts and inclusions) with an electron microprobe to obtain average glass compositions of each single agglutinate.

Description: The impact crater Stickney dominates one hemisphere of the Martian moon Phobos; its diameter (11 km) is about half the size of the body (19 x 22 x 27 km). Besides demarking a threshold between cratering and catastrophic disruption, this impact reveals a great deal about the target's interior. Because Phobos has an unusually low density yet exhibits no direct evidence for volatiles such as water ice, it has been supposed that it sequesters volatiles in the deep interior, or that it is made of some exotic substance, or that it is a loosely-aggregated rubble pile. The network of fracture grooves created by the Stickney impact constrain which, if any, of these models accord with observation.

Description: Thermodynamic properties of the alpha/beta phase transformation in terrestrial troilite (FeS), as measured by differential scanning calorimetry (DSC), vary systematically with prior thermal history of the troilite, as imposed under laboratory conditions. Both the transition temperature and enthalpy change for the alpha/beta transformation decrease with increasing maximum temperature of prior heat treatment. DSC measurements on troilite from various meteorites indicate clear differences in the alpha/beta thermodynamic properties that are consistent with differences in the natural thermal histories of the meteorites.

Description: Mare soil and orange and black pyroclastic glass were reduced in hydrogen gas at temperatures of 900 to 1100 C. The experiments support studies of regolith maturation, lunar volcanism, and the production of oxygen on the moon. The most reactive component in the high-Ti soil was FeO in ilmenite, which was completely reduced to iron metal at all temperatures. Vitreous orange glass crystallized and was partially reduced to iron metal, pyroxene, and minor olivine. Initially devitrified black glass was similarly reduced and converted to iron metal, pyroxene, and minor olivine. The degree of reaction in both glasses increased with temperature.

Description: Pyroclastic glass particles of diverse compositions are sparsely distributed in the lunar soil. Numerous suspected pyroclastic deposits have been pinpointed on the nearside. However, a dearth of pure samples large enough for spectral study, coupled with the kilometer spatial resolution of earth-based telescope spectrometry, makes it difficult to remotely analyze these deposits or tie most glasses to specific sources. Two developments should soon improve this situation. First, the technique of microscope photometry combined with microprobe analysis allows the correlation of the reflectance spectra and chemical compositions of individual glass particles. Second, the Clementine spacecraft will provide multispectral images of pyroclastic deposits at much higher spatial resolution than is currently achievable. These developments, combined with traditional laboratory and telescope studies, should allow compositions of many pyroclastic deposits and sources of many soil particles to be determined.

Description: The temperature and internal energy fields calculated by Takata et al. in the plume are used to calculate the greybody thermal radiation emitted versus wavelength to predict what might be observed by several spectral sensors operating from different platforms when fragments of Comet Shoemaker-Levy 9 (SL-9) impact Jupiter in July 1994. A SPH code was used by Takata et al. to calculate the full three dimensional flow and thermodynamic fields in the comet fragment and the atmosphere of Jupiter. We determined the fragment penetration depth, energy partitioning between the atmosphere and the impactor, and energy density deposited per unit length over the trajectory. Once the impactor had disintegrated and stopped, and the strong atmospheric shock decayed, the flow is driven by buoyancy effects. We then used our SPH code to calculate the flow and thermodynamic fields: pressure, article velocity, temperature, and internal energy distributions in the plume. The calculations for 2 and 10 km cometary fragments yield maximum deposition depths of approximately 175 and 525 km, respectively (1 bar = 0 km depth). We also calculated that 0.7 and 0.6 of the initial kinetic energy of the 10 and 2 km bolides, respectively, are deposited as internal energy in Jupiter's atmosphere.

Description: Fluid whirl and fluid whip are rotor self-excited, lateral vibrations which occur due to rotor interactions with the surrounding fluid. There exist various modes of fluid whirl and fluid whip. These modes are close to rotor modes corresponding to free vibrations (based on the linear model). Small differences are due to nonlinearities in the system. This paper presents experimental and analytical results on the lowest modes of fluid whirls and fluid whip. Examples of rotors supported in fluid lubricated bearings show the variations of rotor deflection amplitudes and phases in the whirl and whip modes with changes of rotative speeds and/or changes in lumped mass locations along the shaft.

Description: The beneficial effects claimed for whirl control devices demonstrate that the dynamic behavior of rotors is influenced by the fluid whirl in shaft and balance drum seals. The present paper reports results from two series of experiments, the first on the factors affecting the whirl at the seal inlet, and the second on the variation of whirl velocity along the seal. In both cases the LDA measurement technique required the clearance between the fixed and rotating parts of the models to be substantially greater than occurs in real machines, but the results are indicative nevertheless. Experimental and theoretical results are given for the radial distribution of whirl velocity in the gap between impeller shroud and pump casing. Results of tests with modified stator surfaces are also shown. This work leads naturally into the second series of experiments where some preliminary measurements of velocity distribution in the clearance between a fixed stator and a rotating shaft are reported for a range of inlet whirl conditions.

Description: Ozone is a key to understanding atmospheric chemistry on Mars. Over 20 photochemical models of the martian atmosphere have been published, and O3 is often used as a benchmark for these models. O3 abundance has been inferred from instrumentation on several spacecraft, with the most complete coverage provided by Mariner 9. The Mariner 9 UV spectrometer scanned from 2100 to 3500 Angstroms with a spectral resolution of 15 Angstroms and an effective field-of-view of approximately 300 sq km. The only atmospheric absorption in the 2000 to 3000 A wavelength region was previously assumed to come from the Hartley band system of ozone, which has an opacity of order unity at winter polar latitudes. Therefore, the amount of ozone was inferred by fitting this absorption feature with laboratory data of ozone absorption. Mars O3 shows strong seasonal and latitudinal variation, with column abundances ranging from 0.2 micron-atm at equatorial latitudes to 60 micron-atm over northern winter polar latitudes (1 micron-atm is a column abundance of 2.689 x 10 exp 15 molecules cm(exp -2)). However, the O3 abundance is never great enough to significantly affect atmospheric temperatures or surface temperatures and frost amounts. Some of the previously-inferred O3 abundances are shown. A radiative transfer computer model is used to re-examine the Mariner 9 UV spectra. Assuming a constant mixing ratio for O3 and no chemical or radiative interaction between O3 and clouds/dust, it is shown that when typical amounts of dust and cloud are present that significant underestimation of O3 abundance occurs. A factor of 3 times as much O3 is needed to generate the same spectrum the spacecraft would measure for a cloudy, dusty atmosphere as for a clear atmosphere. If the scattering properties of martian clouds and dust were well known, then their appearance would not be a problem, as a model would be capable of retrieving the O3 abundance. However, these properties are not well known, which raises doubts about the effectiveness of the UV reflectance spectroscopy technique for measuring O3 abundance on Mars.

Description: The long-time, asymptotic state of rotating homogeneous turbulence at high Reynolds numbers has been examined using large-eddy simulation of the incompressible Navier-Stokes equations. The simulations were carried out using 128 x 128 x 512 collocation points in a computational domain that is four times longer along the rotation axis than in the other directions. Subgrid-scale motions in the simulations were parameterized using a spectral eddy viscosity modified for system rotation. Simulation results show that in the asymptotic state the turbulence kinetic energy undergoes a power-law decay with an exponent which is independent of rotation rate, depending only on the low-wavenumber form of the initial energy spectrum. Integral lengthscale growth in the simulations is also characterized by power-law growth; the correlation length of transverse velocities exhibiting much more rapid growth than observed in non-rotating turbulence.

Description: The objective of this paper is to review our efforts in spatial direct numerical simulations for modeling leading-edge receptivity to freestream sound and vorticity. These results begin to provide the link between the freestream and the initial boundary-layer response and can provide the upstream conditions for further simulations marching through the transition process toward turbulence.

Description: The spatial evolution of disturbances in plane Poiseuille flow and zero pressure gradient boundary layer flow is considered. For disturbances governed by the linearized equations, potential for significant transient growth of the amplitude is demonstrated. The maximum amplification occurs for disturbances with zero or near zero frequencies. Spatial numerical simulations of the transition scenario involving a pair of oblique waves has been conducted for both flows. A fully spectral solver using a simple but efficient fringe region technique allowed the flows to be computed with high resolution into the fully turbulent domain. A modal decomposition of the simulation results indicates that non-linear excitation of the transient growth is responsible for the rapid emergence of low-frequency structures. Physically, this corresponds to streaky flow structures, as seen from the results of a numerical amplitude expansion. Thus, this spatial transition scenario has been found to be similar to the corresponding temporal one. In the boundary layer simulations the streaks are seen to break down from what appears to be a secondary instability.

Description: After publication of the Kolmogorov refined similarity hypotheses, the small-scale intermittency of the energy dissipation field became a central problem in fully developed turbulence (FDT). This phenomenon has been studied in many different ways, e.g. by searching for corrections to scaling exponents in the inertial range velocity structure functions. A direct measure of this intermittency is, however, available by studying the local rate of energy dissipation, and it may be quantitatively characterized by the intermittency exponent mu. As far as we know, nobody has posed an obvious question: Is the intermittency exponent mu a unique constant, i.e., are the values mu(sub kappa), mu(sub epsilon), mu(sub r), mu(sub b), and mu(sub e) the same at high Reynolds numbers, or do they create a set of different (and perhaps independent) exponents? This paper addresses the above question using the high Reynolds number experiments.

Description: Direct simulations of two time-developing turbulent wakes have been performed. Initial conditions for the simulations were obtained from two realizations of a direct simulation of a turbulent boundary layer at momentum thickness Reynolds number 670. In addition, extra two-dimensional disturbances were added in one of the cases to mimic two-dimensional forcing. The unforced wake is allowed to evolve long enough to attain self-similarity. The mass-flux Reynolds number (equivalent to the momentum thickness Reynolds number in spatially developing wakes) is 2000, which is high enough for a short kappa(exp -5/3) range to be evident in the streamwise one-dimensional velocity spectrum. Several turbulence statistics have been computed by averaging in space and over the self-similar period in time. The growth rate in the unforced flow is low compared to experiments, but when this growth-rate difference is accounted for, the statistics of the unforced case are in reasonable agreement with experiments. However, the forced case is significantly different. The growth rate, turbulence Reynolds number, and turbulence intensities are as much as ten times larger in the forced case. In addition, the forced flow exhibits large-scale structures similar to those observed in transitional wakes, while the unforced flow does not.

Description: The present contribution reviews some of the recent progress obtained at our group in the direct numerical simulation (DNS) of compressible boundary layer transition. Elements of the different simulation approaches and numerical techniques employed are surveyed. Temporal and spatial simulations, as well as comparisons with results obtained from Parabolized Stability Equations, are discussed. DNS results are given for flat plate boundary layers in the Mach number range 1.6 to 4.5. A temporal DNS at Mach 4.5 has been continued through breakdown all the way to the turbulent stage. In addition results obtained with a recently developed extended temporal DNS approach are presented, which takes into account some nonparallel effects of a growing boundary layer. Results from this approach are quite close to those of spatial DNS, while preserving the efficiency of the temporal DNS.

Description: The dynamic modeling procedure for large eddy simulation of turbulent flows is reviewed and recent developments in the theoretical aspects and applications are described. Methods for inclusion of backscatter of energy from small to large scale motions are presented. New formulations of the dynamic procedure are proposed which are optimized based on the subgrid scale flux vector or the energy dissipation rate instead of the subgrid scale stress tensor. Recent results from application of the model to forced isotropic turbulence with an inertial subrange, flow over a backward facing step at Reynolds number of 28000, and flow over a concave curved surface are presented.

Description: The energy balance of a lamp varies with the thermal and optical characteristics of the reflector. The photosynthetic radiation efficiency of lamps, defined as input power divided by photosynthetically active radiation (PAR, 400-700 nm) emitted from the lamp ranges between 0.17 and 0.26. The rest of the energy input is wasted as longwave (3000 nm and over) and non-PAR shortwave radiation (from 700 nm to 3000 nm), convective, and conductive heat from the lamp, reflector, and ballast, and simply for increasing the cooling load. Furthermore, some portion of the PAR is uselessly absorbed by the inner walls, shelves, vessels, etc. and some portion of the PAR received by the plantlets is converted into sensible and latent heat. More than 98% of the energy input is probably converted into heat, with only less than 2% of the energy input being converted into chemical energy as carbohydrates by photosynthesis. Therefore, it is essential to reduce the generation of heat in the culture room in order to reduce the cooling load. Through use of a water-cooled reflector, the generation of convective and conductive heat and longwave radiation from the reflector can be reduced, without reduction of PAR.

Description: The use of plate glass-water thermal barriers in controlled environment facilities effectively reduces the thermal load within the plant growth chamber. This allows high PPFs to be provided for plant growth and development studies, adequate simulation of daily light integrals, and simulation of peak PPF's. Further, substantial amounts of incandescent lamp supplementation can be used to achieve simulation of daylight R:FR ratios which are needed to ensure adequate stem development in some species. While the focus in this paper is on the use of entire thermal barriers which separate the lighting enclosure from the plant growth chamber, the same principles apply to the use of water jackets for cooling individual lamps (such as can occur with xenon-arc lamps). In this instance, the barrier separating the lamps from the plant chamber can be much simpler (e.g., plexiglass) as the main function of the barrier is to separate the air ventilation of the lamp enclosure from the air system within the plant growth chamber. The main advantage of water as a thermal barrier is the negligible absorption of radiation in the photosynthetically-active and near infrared wavebands. Consequently, plate glass-water barriers typically allow transmission of approximately 90% of radiation in these regions. While ventilated double and triple glazing systems appear to be attractive alternative to water barriers from an operating standpoint, their significant absorption in the biologically-important wavebands (7 - 12%) with each glass layer and longer-wave cut-offs (typically 2500 - 4000 nm) makes them a much less attractive alternative. The data presented here demonstrate clearly that measurement of PPF alone is not an adequate representation of the radiation environment being used in a controlled environment study. The amounts and proportions of long-wave and short-wave radiation in a plant growth chamber are dependent on lamp type, lamp combination, presence of a thermal barrier, the type of thermal barrier between the lamps and the plant growing area and the overall construction and design of the chamber. It is important, therefore, in reporting results of controlled environment studies, to adequately describe both the details of the lighting system used and the characteristics of the radiation produced by that system, so results of different studies can be adequately evaluated and compared.

Description: Automatic differentiation (AD) is a powerful computational method that provides for computing exact sensitivity derivatives (SD) from existing computer programs for multidisciplinary design optimization (MDO) or in sensitivity analysis. A pre-compiler AD tool for FORTRAN programs called ADIFOR has been developed. The ADIFOR tool has been easily and quickly applied by NASA Langley researchers to assess the feasibility and computational impact of AD in MDO with several different FORTRAN programs. These include a state-of-the-art three dimensional multigrid Navier-Stokes flow solver for wings or aircraft configurations in transonic turbulent flow. With ADIFOR the user specifies sets of independent and dependent variables with an existing computer code. ADIFOR then traces the dependency path throughout the code, applies the chain rule to formulate derivative expressions, and generates new code to compute the required SD matrix. The resulting codes have been verified to compute exact non-geometric and geometric SD for a variety of cases. in less time than is required to compute the SD matrix using centered divided differences.