ALBERT

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

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

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

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

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

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

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

Description: The objective of this work is to develop a space-time accurate numerical method for the solution of incompressible Navier-Stokes equations in generalized coordinates. The resulting code is to be used for direct and large-eddy simulation of turbulence in complex geometries. In a previous paper, the system of Navier-Stokes equations in general curvilinear coordinates was solved by a second-order accurate finite-difference scheme. Satisfactory results were obtained for several flows in two and three dimensions. The system of Navier-Stokes for the fluxes are given in Orlandi (1989). The main deficiency of the numerical scheme was the large CPU time required for the solution of the Poisson equation for the 'pressure' field. The point SOR relaxation, in conjunction with a multigrid scheme, was used for the Poisson equation. In some cases, particularly with very fine grids, it was impossible to obtain a divergent-free flow. A preliminary attempt is made to compute the spatially evolving flow of Swearingen & Blackwelder. To reduce the streamwise distance, the inflow was at a distance x = 60 cm from the leading edge.

Description: Vortex filaments in superfluids such as helium 2 may provide new insights into very high Reynolds number flows. The behavior of a superfluid vortex ring interacting with a normal fluid shear flow, specifically channel flow, is simulated. The vortex ring evolves into a stable horseshoe configuration which propagates without further change of form. In this simulation, a boundary layer behavior in a superfluid through the coupling of the superfluid and the normal fluid is demonstrated.

Description: Explicit solutions of the stationary Hopf equation are discussed and their computational possibilities are explored. The motivation is to circumvent the infinite hierarchy of coupled equations for the velocity moments and obtain an exact closure of the steady-state 3D Navier-Stokes equations, without modeling assumptions or truncation. The Hopf formulation of the Navier-Stokes equation is reviewed. A stationary homogeneous solution for 2D flow is displayed and discussed. It is shown how depletion of nonlinearity may arise for 3D forced homogeneous flow. The general 3D forced case is considered and a method for closing the 3D unforced equations with arbitrary boundary conditions is derived.

Description: The solar nebula, from which the planets in our solar system formed, featured a disk of gas and dust grains in rapid, differential rotation, and at some stage was likely to have been unstable to thermal convection. This situation is suspected by many to lead to significant turbulent Reynolds stress production and angular momentum transport in such systems, and estimates of transport rates have been attempted from unsubstantiated phenomenological models. In order to determine the circumstances and physical conditions under which our own planetary system formed and to explain recent observations of young stellar systems, it is necessary to develop realistic models of heat and angular momentum transport for such flows. Developing an understanding of complicated flows featuring thermal convection, rotation, and shear is also of wide interest in stellar astrophysics and in planetary and terrestrial atmospheric studies. The ultimate objective is to develop workable models based on the numerical simulations for constructing global solar nebula models; viz., relatively simple prescriptions for heat and angular momentum fluxes from given system parameters (e.g., ratios of rotation, shear, and convective lapse rates) are characterized, quantified, and developed. Toward this end, our program has been to attempt to understand the behavior of the direct numerical simulations of Boussinesq convection, which, despite the complexity of the results, is still an overly simplified approximation to the real system and should be more amenable to analysis. These results are also intended to be tested against turbulence models, especially those designed for atmospheric boundary layers, and may provide a basis for subgrid-scale models. In order to make the numerical simulations more realistic with regard to the solar nebula problem, a fully compressible code that will allow incorporation of large density stratifications and realistic thermodynamic and radiative properties is developed. In order to explore the properties of these flows at the very high values of Re found in natural systems and the very low values of Pr found in most astrophysical contexts, we will need to employ large-eddy simulations for which we want to determine the most appropriate subgrid-scale model to incorporate.

Description: Turbulent plane Couette flow was numerically simulated at a Reynolds number (U(sub w)h/nu) of 6000, where U(sub w) is the relative wall speed and h is half the channel-height. Unlike in Poiseuille flow, where the mean shear rate changes its sign at the centerline, the sign of mean shear rate in plane Couette flow remains the same across the whole channel. This difference is expected to yield several differences between the two flows, especially in the core region. The most significant and dramatic difference observed was the existence of large-scale structures in the core region of the plane Couette flow. The large eddies are extremely long in the flow direction and fill the entire channel (i.e., their vertical extent is 2h). The large-scale structures have the largest contribution from the wavenumber (k(sub x)h,k(sub z)h) = (0, plus or minus 1.5), corresponding to a wavelength lambda(sub z)/h is approximately equal to 4. The secondary motion associated with the k(sub x)h = 0 mode consists of the large-scale vortices. The large eddies contribute about 30 percent of turbulent kinetic energy.

Description: The National Aerospace Plane (NASP) will require thermal insulation systems which are consistent with cryogenic fluids, high thermal loads, and design restrictions such as weight and volume. Test sections of the proposed system have been constructed and evaluated. In this paper we discuss the components of the insulation system, the application of the insulation system to the NASP liquid hydrogen fuel tank system, and thermal conductivity measurements performed on test sections of the system. Both steady-state and transient thermal measurements are presented.

Description: A model is developed for describing the interaction of vortex-drop clusters in a flowing gaseous jet, convecting downstream from an injection location. Results are presented for a stationary case representing the situation when identical clusters are continuously injected and the injection rate is constant. The results indicate that, in a rich mixture high-drop-number density regime, the mass evaporated from the drops controls the velocity of the cluster-in-vortex as it propagates downstream.

Description: An algorithm is presented for unsteady two-dimensional incompressible Navier-Stokes calculations. This algorithm is based on the fourth order partial differential equation for incompressible fluid flow which uses the streamfunction as the only dependent variable. The algorithm is second order accurate in both time and space. It uses a multigrid solver at each time step. It is extremely efficient with respect to the use of both CPU time and physical memory. It is extremely robust with respect to Reynolds number.

Description: A Navier-Stokes equations solver based on a pressure correction method with a pressure-staggered mesh and calculations of separated three-dimensional flows are presented. It is shown that the velocity pressure decoupling, which occurs when various pressure correction algorithms are used for pressure-staggered meshes, is caused by the ill-conditioned discrete pressure correction equation. The use of a partial differential equation for the incremental pressure eliminates the velocity pressure decoupling mechanism by itself and yields accurate numerical results. Example flows considered are a three-dimensional lid driven cavity flow and a laminar flow through a 90 degree bend square duct. For the lid driven cavity flow, the present numerical results compare more favorably with the measured data than those obtained using a formally third order accurate quadratic upwind interpolation scheme. For the curved duct flow, the present numerical method yields a grid independent solution with a very small number of grid points. The calculated velocity profiles are in good agreement with the measured data.

Description: The least squares (L sub 2) finite element method is introduced for 2-D steady state pure convection problems with smooth solutions. It is proven that the L sub 2 method has the same stability estimate as the original equation, i.e., the L sub 2 method has better control of the streamline derivative. Numerical convergence rates are given to show that the L sub 2 method is almost optimal. This L sub 2 method was then used as a framework to develop an iteratively reweighted L sub 2 finite element method to obtain a least absolute residual (L sub 1) solution for problems with discontinuous solutions. This L sub 1 finite element method produces a nonoscillatory, nondiffusive and highly accurate numerical solution that has a sharp discontinuity in one element on both coarse and fine meshes. A robust reweighting strategy was also devised to obtain the L sub 1 solution in a few iterations. A number of examples solved by using triangle and bilinear elements are presented.

Description: An implicit numerical algorithm for the time accurate solution of the compressible Navier-Stokes equations is described. Results for steady flow past a finite flat plate are presented, together with preliminary results for the temporal simulation of second mode instability in a flat plate boundary layer at Mach 4.5.

Description: The main objective of this research is to address two important but unresolved problems: (1) the measurement of vertical and transverse length scales via space correlations for all Reynolds stress components and velocity-temperature correlations, both in the free stream and within the boundary layer using the existing triple and quad-wire probes; and (2) to relate the character of the free stream turbulence to the character of the turbulence within the boundary layer in order to determine the effect on surface heat transfer.

Description: An experimental research project aimed at obtaining quantitative data on the behavior of the secondary vortex structure in a turbulent mixing layer at moderate Reynolds numbers (Re(sub delta) = 2.9 x 10(exp 4)) is discussed. This project was terminated before all the contemplated measurements could be made, and data were obtained only on the spatially stationary part of the secondary structure. Nonetheless, these data reveal some interesting facets of mixing layer behavior which are discussed.

Description: The wavelet-transformed Navier-Stokes equations are used to define quantities such as the transfer of kinetic energy and the flux of kinetic energy through scale r at position x. Direct numerical simulations of turbulent shear flow reveal that although their mean spatial values agree with their traditional counterparts in Fourier space, their spatial variability at every scale is very large, exhibiting non-Gaussian statistics. The local flux of energy involving scales smaller than some r also exhibits large spatial intermittency, and it is negative quite often, indicative of local inverse cascades.

Description: Large-scale structures in turbulent and transitional wall-bounded flows make a significant contribution to the Reynolds stress and turbulent energy. The behavior of these structures is examined. Small perturbations are introduced into a laminar and a turbulent boundary layer to trigger the formation of large-scale features. Both flows use the same inlet unit Reynolds number, and they experience the same pressure gradient history, i.e. a favorable pressure gradient (FPG) followed by an adverse pressure gradient (APG). The perturbation consists of a small short duration flow repetitively introduced through a hole in the wall located at the C(sub p) minimum. Hot-wire data are averaged on the basis of the phase of the disturbance, and automation of the experiment was used to obtain measurements on large spatially dense grids. In the turbulent boundary, the perturbation evolves into a vortex loop which retains its identity for a considerable streamwise distance. In the laminar layer, the perturbation decays to a very small magnitude before growing rapidly and triggering the transition process in the APG. The 'time-like' animations of the phase-averaged data are used to gain insight into the naturally occurring physical mechanisms in each flow.

Description: The dependence of product generation on Peclet and Reynolds numbers in a numerically simulated, reacting, two dimensional, temporally growing mixing layer is related theoretically to the fractal dimension of the passive scalar interfaces. This relation is verified using product generation measurements and dimensions derived from a standard box counting technique. A transition from a low initial dimension to a higher one of approximately 5/3 is identified and shown to be associated to the kinematic distortion on the flow field during the first pairing interaction. It is suggested that the structures responsible for this transition are non-deterministic, non-random, inhomogeneous fractals. Only the large scales are involved. No further transitions, either in the spectra of the vorticity field or in the mixing behavior, are found for Reynolds numbers up to 90,000.

Description: Studies of supersonic mixing which were accomplished over the last year with Center for Turbulence Research (CTR) support are described. During this period, a Nd:YAG laser, optical components, and data acquisition computer were obtained. This allowed detailed visualizations of the flow structure to be performed at a rapid rate, representing a significant improvement over our previous attempts. Aspects of the flow structure are described below. In addition, preliminary findings on a possible mixing enhancement strategy are also shown using the flow visualization technique.

Description: A formulation has been developed to describe the evaporation of dense or dilute clusters of binary-fuel drops. The binary fuel is assumed to be made of a solute and a solvent whose volatility is much lower than that of the solute. Convective flow effects, inducing a circulatory motion inside the drops, are taken into account, as well as turbulence external to the cluster volume. Results obtained with this model show that, similar to the conclusions for single isolated drops, the evaporation of the volatile is controlled by liquid mass diffusion when the cluster is dilute. In contrast, when the cluster is dense, the evaporation of the volatile is controlled by surface layer stripping, that is, by the regression rate of the drop, which is in fact controlled by the evaporation rate of the solvent. These conclusions are in agreement with existing experimental observations. Parametric studies show that these conclusions remain valid with changes in ambient temperature, initial slip velocity between drops and gas, initial drop size, initial cluster size, initial liquid mass fraction of the solute, and various combinations of solvent and solute. The implications of these results for computationally intensive combustor calculations are discussed.

Description: It is found through two-dimensional temporal simulations of high-speed free shear layers that mean flow distortion is significantly increased when supersonic disturbances are introduced as initial conditions. The shear layer exhibits no subharmonic growth or roll-up, but rather a spectral broadening as energy is distributed into higher harmonics. Increasing the velocity of one side of the mixing layer (u2) to roughly 1/5 the speed of the high speed side (u1), allows a slight subharmonic growth at a very slow rate for two-dimensional modes. A first look at three-dimensional free shear flows is also presented for M = 2. No effect is seen for incompressible flow; however, stabilization is seen with respect to maximum temporal growth rates as the transverse velocity is increased. A much stronger, but similar effect is seen if u2 is increased. The wave direction of maximum growth for u2 is greater than 0.4 is found to be the direction of the faster stream (u1) over a broad range of transverse velocities.

Description: The effects of rotation on grid turbulence were studied in a low speed tunnel. Results are presented for both axisymmetric and non-axisymmetric flows. As observed by previous workers, the cascade process was seen to be effectively blocked by rotation leading to reduced dissipation. In the non-axisymmetric case, rotation caused the flow to tend rapidly towards axisymmetry; however, for the length studied, the flow tended to move away from isotropy. Based on the results, some recommendations are made on the design of a larger facility for studying similar flows.

Description: It is shown that the ideas of selective amplification and direct resonance, based on linear theory, can not provide an explanation for the well-defined streak spacing of about 100 wall units (referred to as 100(+) hereafter) in wall-bounded turbulent shear flows. In addition, for the direct resonance theory, the streaks are created by the non-linear self-interaction of the vertical velocity rather than of the directly forced vertical vorticity. In view of the failure of these approaches, it is then proposed that the selection mechanism must be inherently non-linear and correspond to a self-sustaining mechanism. The 100(+) value should thus be considered as a critical Reynolds number for that mechanism. Indeed, in the case of Poiseuille flow, this 100(+) criterion for transition to turbulence corresponds to the usually quoted value of 1000 based on the half-width and the centerline velocity. In Couette flow, it corresponds to a critical Reynolds number of about 400 based on the half width and half velocity difference.

Description: The development of an automated technique for the eduction of 3-D spatial patterns in vector or scalar diagnostics was completed. The method is based on an iterative convolution between a trial pattern and the data field. It was applied to the analysis of low Reynolds number turbulent channel flow and homogeneous shear flow. The results yielded new information on the dominant flow structures in these flows, particularly with respect to the spatial relationships between various forms of organized motion. A particular application of the pattern eduction method, which is tentatively referred to as an 'adaptive wavelet transformation', is proposed with the objective of investigating the way turbulence structure changes with scale. Preliminary results using data from homogeneous turbulent shear flow simulations are presented. At the low Reynolds numbers of the simulations, there is no evidence of scale similarity. The small scales appear to be associated with the edges of the larger scale vortical structures.

Description: The existence and importance of large-scale spanwise vortical structures for 2-D straight mixing layers has been well documented in the last decade. Computer models and simulations have sought to reproduce these vortical structures associated with the Kelvin-Helmholtz (K-H) instability mode which is due to the shear per se. Secondary streamwise vortical structures for the same flows were also seen experimentally and have recently been given importance in computational efforts. Curved mixing layers can be characterized as stable (the high-speed stream is placed on the outside of the longitudinal bend), leading to a suppression of the Taylor-Gortler (T-G) instability, and unstable (high-speed stream on the inside of the bend), leading to an enhancement of the T-G instability. The T-G instability is associated with the centripetal acceleration that the curvature imparts. Thus, curvature superimposed on 2-D shear layer flows provides a way for studying the importance of streamwise vorticity, its competition with spanwise vorticity, and changes to entrainment and mixing. Furthermore, the outcome of the competition of a relatively enhanced or suppressed T-G instability with the K-H instability offers the possibility of achieving passive mixing enhancement. As a first step in understanding the competition between the K-H and the T-G instabilities and the resulting changes to the structure of the flow, highly resolved visualizations of the flow structure for the stable and the unstable configurations are provided. The straight layer is also visualized for comparison with earlier works.

Description: Computer flow simulation aided by dynamical systems analysis is used to investigate the kinematics of time-periodic vortex shedding past a two-dimensional circular cylinder in the context of the following general questions: (1) Is a dynamical systems viewpoint useful in the understanding of this and similar problems involving time-periodic shedding behind bluff bodies; and (2) Is it indeed possible, by adopting such a point of view, to complement previous analyses or to understand kinematical aspects of the vortex shedding process that somehow remained hidden in previous approaches. We argue that the answers to these questions are positive. Results are described.

Description: An ongoing research effort designed to apply the best possible second-moment-closure model to simulate complex hypersonic flows is presented. The baseline model under consideration is the Launder-Reece-Rodi Reynolds stress transport turbulence model. Two add-ons accounting for wall effects, namely, the Launder-Shima low-Reynolds-number model and the compressible wall-function technique, are tested. Results are reported for flow over a flat plate, both adiabatic-wall and cooled-wall cases. It has been found that further improvements of the existing models are necessary to achieve accurate prediction in high Mach number flow range.

Description: Scalar fields undergoing random advection have attracted much attention from researchers in both the theoretical and practical sectors. Research interest spans from the study of the small scale structures of turbulent scalar fields to the modeling and simulations of turbulent reacting flows. The probability density function (PDF) method is an effective tool in the study of turbulent scalar fields, especially for those which involve chemical reactions. It has been argued that a one-point, joint PDF approach is the one to choose from among many simulation and closure methods for turbulent combustion and chemically reacting flows based on its practical feasibility in the foreseeable future for multiple reactants. Instead of the multi-point PDF, the joint PDF of a scalar and its gradient which represents the roles of both scalar and scalar diffusion is introduced. A proper closure model for the molecular diffusion term in the PDF equation is investigated. Another direction in this research is to study the mapping closure method that has been recently proposed to deal with the PDF's in turbulent fields. This method seems to have captured the physics correctly when applied to diffusion problems. However, if the turbulent stretching is included, the amplitude mapping has to be supplemented by either adjusting the parameters representing turbulent stretching at each time step or by introducing the coordinate mapping. This technique is still under development and seems to be quite promising. The final objective of this project is to understand some fundamental properties of the turbulent scalar fields and to develop practical numerical schemes that are capable of handling turbulent reacting flows.

Description: The immediate goal is to understand and validate the Yakhot-Orszag model of the velocity-derivative skewness and model equation for the rate of energy dissipation epsilon. A summary of a more detailed manuscript in preparation is presented. The purpose is to clarify some limitations of the theory by careful examination of key assumptions and approximations, and thereby to encourage its improvement.

Description: Investigations carried out for revisiting homogeneous turbulent flows in the presence of mean shear, rotation, or external compression are summarized. The simplest and most concise RDT (Rapid Distortion Theory) formulation, which includes a comprehensive linear stability analysis, is used for this purpose. Such a linear approach could be extended by a generalized EDQNM (Eddy Damped Quasi-Normal Markovian) to two point closure in order to model non-linear interactions, especially when pure Coriolis effects are present. The results are discussed in connection with databases obtained by DNS (Direct Numerical Simulations), including previous Center for Turbulence Research (CTR) results and new calculations in progress. The main goal is to contribute to and significantly improve on the rational one-point closure models in progress at the CTR and at Ecole Centrale de Lyon (ECL).

Description: The near-wall region plays an essential role in turbulent boundary layers: it is a region of high shear; the peak rate of production and peak intensity of turbulence occurs there; and the peak rate of dissipation occurs right at the wall. Nevertheless, this region has received less attention from modelers than have more nearly homogeneous flows. One reason for this is that when the boundary layer is near equilibrium, experimental data can be used to prescribe the flow in the wall layer. Another reason is that most turbulence models are developed under assumptions of near homogeneity. This is a poor approximation in the wall region. A single-point moment closure model for the strongly non-homogeneous A turbulent flow near a rigid boundary is developed.

Description: A prototype ceramic fabric/titanium water heat pipe has been constructed and tested; it transported 25 to 80 W of power at 423 K. Component development and testing is continuing with the aim of providing an improved prototype, with a 38 micron stainless steel liner covered by a biaxially-braided Nextel (trademark) sleeve that is approximately 300 microns thick. This fabric has been tested to 800 K, and its emittance is about 0.5 at that temperature. Advanced versions of the water heat pipe will probably require a coating over the ceramic fabric in order to increase this emittance to the 0.8 to 0.9 range.

Description: Preliminary results are presented of observations of Mars using three telescopes and a CCD camera, grating-array spectrometers, and a near-IR array camera. Martian albedo features are shown to be similar in the visible and 1.3-3.5-micron range, and reflectances of nearly zero are noted in the image at 3.02 microns. The results demonstrate that these techniques provide effective maps of spatial variations in spectral absorption features caused by Martian atmospheric volatiles and by surface minerals.

Description: Recent advances in the methodology for direct numerical simulation of turbulent flows and some of the current applications are reviewed. It is argued that high-order finite difference schemes yield solutions with comparable accuracy to the spectral methods with the same number of degrees of freedom. The effects of random inflow conditions on the downstream evolution of turbulence are discussed.

Description: A variable explicit time integration algorithm is developed for unsteady diffusion problems. The algorithm uses nodal partitioning and allows the nodal groups to be updated with different time steps. The stability of the algorithm is analyzed using energy methods and critical time steps are found in terms of element eigenvalues with no restrictions on element types. Several numerical examples are given to illustrate the accuracy of the method.

Description: A methodology for correcting raw pressure-drop data on the influence of acceleration on the instrumentation in two-phase flow thermal management systems is described. Such systems are now being considered as an alternative to conventional single-phase systems in future space missions because of the potential to reduce overall system mass, size, and pumping power requirements. Corrected pressure-drop measurements are presented and compared with predictions of two-phase flow pressure-drop models. A set of flow and acceleration conditions is defined for which the frictional pressure drop does not increase upon entry into 0-'g' pressure drop.

Description: This article describes the results of numerical simulations of oscillating wall-bounded developing flows. The full phase-averaged Navier-Stokes equations are solved. The application of quasi-steady turbulence modeling to unsteady flows is demonstrated using an unsteady version of the k-epsilon model. The effects of unsteadiness on the mean flow and turbulence are studied. Critical evaluation of the applicability of the quasi-steady approach to turbulence modeling is presented. Suggestions are given for the future efforts in turbulence modeling of unsteady flows.

Description: Propulsion systems planned for use late in this century and beyond will require appropriate physical models for describing supersonic combustion and numerical techniques for solving the model governing equations. A computer program to study these flows is reported which considers the multicomponent diffusion and convection of important chemical species, the finite-rate reaction of these species, and the resulting interaction of the field mechanics and the chemistry. The application of the program to a spatially developing and reacting mixing layer, which serves an an excellent physical model for the mixing and reaction processes that take place in a scramjet combustor, is reported. Several techniques to enhance the fuel-air mixing and growth of that layer and improve its overall combustion efficiency are considered.

Description: A three-component laser velocimeter is used to investigate the flow over a backward-facing step. The backward-facing step had an expansion ratio of 2, a boundary layer height to step height ratio of 0.34 and a Reynolds number based on step height of 19,000. Results from three-component velocimeter surveys of the flow over the backward-facing step are presented with comparisons of the current experiment with previous experiments and computational results. The present results compared well with previous experiments with the exception of the reattachment length. The short reattachment length was due to the short length of the channel downstream. The measurement of the lateral velocity component showed that there is a mean flow in and out of the centerline plane as high as 7 percent of the freestream velocity. However, the shear stresses show no correlation between the lateral fluctuations and the longitudinal and vertical fluctuations, indicating that the flow is 2D in terms of the turbulence quantities.

Description: A method presented by Schlosser et al. (1989) for analyzing the pressure dependence of experimental melting-temperature data is applied to rare-gas solids. The plots of the logarithm of the reduced melting temperature vs that of the reduced pressure are straight lines in the absence of phase transitions. The plots of the reduced melting temperatures for Ar, Kr, and Xe are shown to be approximately straight lines.

Description: The physical mechanism driving the weakly chaotic Taylor-Couette flow is investigated using the short-time Liapunov exponent analysis. In this procedure, the transition from quasi-periodicity to chaos is studied using direct numerical 3D simulations of axially periodic Taylor-Couette flow, and a partial Liapunov exponent spectrum for the flow is computed by simultaneously advancing the full solution and a set of perturbations. It is shown that the short-time Liapunov exponent analysis yields more information on the exponents and dimension than that obtained from the common Liapunov exponent calculations. Results show that the chaotic state studied here is caused by a Kelvin-Helmholtz-type instability of the outflow boundary jet of Taylor vortices.

Description: The complex nature of the Uranus radio emissions, both magnetospheric and atmospheric, is reviewed, with emphasis on the identification of distinct components and the determination of their source locations. Seven radii components were discovered in addition to the RF signature of lightning in the planet's atmosphere. Six of the seven magnetospheric components are freely propagating emissions; one component, the nonthermal continuum, is trapped in the density cavity between the magnetopause and the dense inner magnetosphere. The radio components are divided into two types according to their emission signature: bursty emission and smooth emission. The inferred source location for the dominant nightside emission is above the nightside magnetic pole, largely overlapping the UV auroral region and the magnetic polar cap. The N-burst component appears to be associated with solar-wind enhancements at Uranus, consistent with the idea that the solar wind was triggering magnetospheric substormlike activity during the encounter.

Description: The energetic particle measurements by the low-energy charged-particle and cosmic-ray instruments on the Voyager 2 spacecraft in the magnetosphere of Uranus are reviewed. Upstream events were observed outside the Uranian bow shock, probably produced by ion escape from the magnetosphere. Evidence of earthlike substorm activity was discovered within the Uranian magnetosphere. A proton injection event was observed within the orbit of Umbriel and proton events were observed in the magnetotail plasma-sheet boundary layer that are diagnostic of earthlike substorms. The magnetospheric composition is totally dominated by protons, with only a trace abundance of H(2+) and no evidence for He or heavy ions; the Uranian atmophere is argued to be the principal plasma source. Phase-space densities of medium energy protons show inward radial diffusion and are quantitatively similar to those observed at the earth, Jupiter, and Saturn. These findings and plasma wave data suggest the existence of structures analogous to the earth's plasmasphere and plasmapause.

Description: An overview of the observational results on the plasma environment at Uranus is given, and the implications of these observations for magnetospheric physics at Uranus are discussed. During the Voyager 2 encounter with Uranus, an extended magnetosphere filled with a tenuous plasma was detected. This low-energy plasma was found to consist of protons and electrons, with no significant heavy ion contribution, and with a density in the regions sampled by the spacecraft of at most three electrons per cubic centimeter. The plasma electrons and ions exhibit both a thermal component (with temperatures of tens of eV) and a hot component (with temperatures of a few keV). The thermal ion component is observed both inside and outside an L-shell value near 5, whereas the hot ion and electron component is excluded from the region inside of that L-shell. The source of the thermal component of the plasma is either the planetary ionosphere or the neutral hydrogen corona surrounding Uranus, whereas the hot component is convected in from the magnetotail, with probably an ionospheric source.

Description: Recent telescopic observations have led to the identification of cyanogroup-containing molecules in the dark surface solids of several D-class asteroids, cometary dusts, and the rings of Uranus, as well as the low-albedo atmosphere of Iapetus. The occurrence of the 2.2-micron overtone of C triple-bond N's stretching fundamental mode in all four classes of small solar system bodies is presently suggested to serve as a diagnostic of both exposure duration and degree of modification of surface materials.

Description: Results are presented from experimental and numerical-modeling studies carried out on a rotating thermally driven fluid system in a cylindrical annulus with the horizontal gradient imposed upon the lower horizontal surface. The rotation rate and the temperature difference were varied to construct a regime diagram in thermal Rossby-number/Taylor-number space. It is shown that the curve separating the axisymmetric flow from the wave flow is 'knee-shaped', similar to the conventional side-heated and side-cooled baroclinic annulus. The experimentally observed transition curve agrees well with numerical calculations and the agreement between the predicted and measured wavenumbers is good both near the transition and within the wave regime. Away from the transition curve and well within the wave regime, there is an evolution from larger to smaller wavenumbers as the flow equilibrates, which is well captured by the numerical model for these cases.

Description: The fine-scale structure of turbulence in a fully developed turbulent duct flow is examined by considering the 3D velocity derivative field obtained from direct numerical simulations at two relatively small Reynolds numbers. The magnitudes of all mean-square derivatives (normalized by wall variables) increase with the Reynolds number, the increase being largest at the wall. These magnitudes are not consistent with the assumption of local isotropy except perhaps near the duct center-line. When the assumption of local isotropy is relaxed to one of local axisymmetry, or invariance with respect to rotation about a coordinate axis (here chosen in the streamwise direction), satisfactory agreement is indicated by the data outside the wall region. Support for axisymmetry is demonstrated by anisotropy invariant maps of the dissipation and vorticity tensors.

Description: A quantitative analysis of error sources in 1D planetary photoclinometry is presented. The technique is affected by error sources arising from the spacecraft image, the planetary body, and the scan line orientation. Slope errors are calculated for each of these sources, using examples of Voyager imaging of Ganymede and Viking orbiting imaging of Mars. Slope errors are investigated for a variety of viewing and lighting geometries, slope angles, and slope orientations. The results are broken down into nonsystematic and systematic errors. Derivations that allow the calculation of photoclinometric slope errors for any photometric function are presented, and the implications of these results for 2D photoclinometric techniques are discussed.

Description: This paper presents a general solution algorithm for the set of difference equations that arise when two-point central differences are used to approximate the flux difference terms in systems of hyperbolic differential equations. The general algorithm eliminates the weak points associated with the nonstandard algorithm reported by Wornom and Hafez (1986). The disadvantages of their algorithm relate to its implementation. It consists of separate algorithms for subsonic, supersonic, sonic and shock cells, applied individually, which presents a major bookkeeping problem when multiple sonic and shock cells are present. The general algorithm eliminates this problem and introduces an improved shock treatment which produces shocks with at most one interior shock point.

Description: A morphometric analysis of north-south trending graben at 35 deg N in northern Tharsis is used to provide the first quantitative measurements of crustal extension in and strain for areas of the Martian crust that have undergone both extreme and mild deformation. These estimates also demonstrate the viability of the Viking data for making such quantitative assessment. Graben in the region accommodate about 8 km of post Early Amazonian east-west extension. This extension corresponds to a net regional strain of 0.45-0.61 percent; local strains vary from less than 1 percent to 5 percent. Extension is nonuniform and localized near 110 deg (south of Alba Patera and north of Ceraunius Fossae) because of reactivation of preexisting fractures in underlying basement and possibly due to superposition of stress fields from Tharsis and Alba Patera. A major mechanical discontinuity in the shallow crust, probably representing the base of or a strength contrast within the megaregolith, is observed. It occurs at a depth of 2 km on the profile flanks and is depressed to about 7.5 km near 110 deg due to loading of the crust by the volcano.

Description: The dynamic subgrid-scale (SGS) model of Germano et al. (1991) is generalized for the large eddy simulation (LES) of compressible flows and transport of a scalar. The model was applied to the LES of decaying isotropic turbulence, and the results are in excellent agreement with experimental data and direct numerical simulations. The expression for the SGS turbulent Prandtl number was evaluated using direct numerical simulation (DNS) data in isotropic turbulence, homogeneous shear flow, and turbulent channel flow. The qualitative behavior of the model for turbulent Prandtl number and its dependence on molecular Prandtl number, direction of scalar gradient, and distance from the wall are in accordance with the total turbulent Prandtl number from the DNS data.

Description: The momentumless coalescence of drops of the same liquid, separated by an immiscible host, is studied experimentally. Observations show that for low-viscosity drops of unequal sizes, there is considerable mixing following coalescence, with the smaller drop penetrating the larger drop as a vortex. The extreme case of coalescence of a small drop with the bulk of the same liquid at a flat interface with an immiscible liquid is studied in detail. The penetration depths of small drops (1-5 mm) following coalescence are measured and correlated with theoretical predictions. It is found that in the range of the investigation, the penetration depth is proportional to the 5/4 power of drop diameter and inversely proportional to the square root of the drop viscosity.

Description: Dynamical aspects of a drop drastically flattened by acoustic radiation stress are considered. Its static equilibrium has been studied, starting with a dislike shape and modeling the sound field and the associated radiation stress according to this geometry. It is suggested that, at low viscosity, the ripples are capillary waves generated by the parametric instability excited by the membrane vibration, which is driven by the sound pressure. Atomization occurs whenever the membrane becomes so thin that the vibration is sufficiently intense. Buckling occurs when an existent equilibrium is unstable to a radial oscillation of the membrane because of the Bernoulli effect. The radiation stress at the rim of the flattened drop is also destabilizing and leads to horizontal expansion and subsequent breakup.

Description: A time-accurate numerical solution was carried out for transient radiative cooling of a gray emitting and absorbing medium in a square two-dimensional region. The integro-differential energy equation for transient temperature distributions was solved in two stages. At each time increment, the local radiative source term was obtained by numerical integration of the temperature field using two-dimensional Gaussian integration over rectangular subregions. Then the differential portion of the equation was integrated forward in time by use of the local first and second time derivatives. The results were compared with available limiting case, and excellent agreement was obtained. Transient results are given for a wide range of optical thicknesses of the region. Optimum transient cooling is obtained when the optical side length is about 4.

Description: A variety of numerical simulations of transition and turbulence in incompressible flow are presented to compare the commonly used rotation form with the skew-symmetric (and other) forms of the nonlinear terms. The results indicate that the rotation form is much less accurate than the other forms for spectral algorithms which include aliasing errors. For de-aliased methods the difference is minimal.

Description: Three-component laser velocimetry techniques are applied to characterize the fluid dynamics of an MOCVD reactor. These methods provide three-dimensional quantitative measurements of the gas velocities inside the reactor at typical growth temperatures. The effects of buoyancy-induced convection are examined by comparing data from the hot reactor with cold reactor results. It was found that thermal convection dominated the fluid dynamics of the reactor at growth temperatures. Cold reactor tests showed unstable flow patterns that were subject to fluidic switching effects.

Description: This paper presents numerical solutions of jet-induced mixing in a partially full cryogenic tank. An axisymmetric laminar jet is discharged from the central part of the tank bottom toward the liquid-vapor interface. Liquid is withdrawn at the same volume flow rate from the outer part of the tank. The jet is at a temperature lower than the interface, which is maintained at a certain saturation temperature. The interface is assumed to be flat and shear free and the condensation-induced velocity is assumed to be negligibly small compared with radial interface velocity. Finite-difference method is used to slove the nondimensional form of steady-state continuity, momentum, and energy equations. Calculations are conducted for jet Reynolds numbers ranging from 150 to 600 and Prandtl numbers ranging from 0.85 to 2.65. The effects of previously stated parameters on the condensation Nusselt and Stanton numbers that characterize the steady-state interface condensation process are investigated. Detailed analysis is performed to gain a better understanding of the fundamentals of fluid mixing and interface condensation.

Description: The Martian southern hemisphere atmospheric water vapor column abundance measurements reported agree with Viking Orbiter atmospheric water detectors during early southern spring and southern autumnal equinox; profiles obtained in southern mid- and late summer, however, indicate the presence of twice as much water both in the southern hemisphere and planetwide. This discrepancy is accounted for by the high optical depths created by two global dust storms during the Viking year, while the present observations were obtained in the case of the relatively dust-free atmosphere of the 1988-1989 opposition.

Description: Based on multilevel partitioning, this paper develops a structural parallelizable solution methodology that enables a significant reduction in computational effort and memory requirements for very large scale linear and nonlinear steady and transient thermal (heat conduction) models. Due to the generality of the formulation of the scheme, both finite element and finite difference simulations can be treated. Diverse model topologies can thus be handled, including both simply and multiply connected (branched/perforated) geometries. To verify the methodology, analytical and numerical benchmark trends are verified in both sequential and parallel computer environments.

Description: An algorithm for the solution of the incompressible Navier-Stokes equations in three-dimensional generalized curvilinear coordinates is presented. The algorithm can be used to compute both steady-state and time-dependent flow problems. The algorithm is based on the method of artificial compressibility and uses a third-order flux-difference splitting technique for the convective terms and the second-order central difference for the viscous terms. The accuracy is obtained in the numerical solutions by subiterating the equations in pseudotime for each physical time step. The equations are solved with a line-relaxation scheme that allows the use of very large pseudotime steps leading to fast convergence for steady-state problems as well as for the subiterations of time-dependent problems. The steady-state solution of flow through a square duct with a 90-deg bend is computed, and the results are compared with experimental data. Good agreement is observed. Computations of unsteady flow over a circular cylinder are presented and compared to other experimental and computational results. Finally, the flow through an artificial heart configuration with moving boundaries is calculated and presented.

Description: Near-infrared spectra of a bright and a dark thermal emission feature on the night side of Venus have been obtained from 2.2 to 2.5 microns at a spectral resolution of 1200 to 1500. Both bright and dark features show numerous weak absorption bands produced by CO2, CO, water vapor, and other gases. The bright feature emits more radiation than the dark feature throughout this spectral region, but the largest contrasts occur between 2.21 s 2.32 microns, where H2SO4 clouds and a weak CO2 band provide the only known sources of extinction. The contrast decreases by 55 to 65 percent at wavelengths longer than 2.34 microns, where CO, clouds, and water vapor also absorb and scatter upwelling radiation. This contrast reduction may provide direct spectroscopic evidence for horizontal variations in the water vapor concentrations in the Venus atmosphere at levels below the cloud tops.

Description: Magellan images confirm that volcanism is widespread and has been fundamentally important in the formation and evolution of the crust of Venus. High-resolution imaging data reveal evidence for intrusion (dike formation and cryptodomes) and extrusion (a wide range of lava flows). Also observed are thousands of small shield volcanoes, larger edifices up to several hundred kilometers in diameter, massive outpourings of lavas, and local pyroclastic deposits. Although most features are consistent with basaltic compositions, a number of large pancake-like domes are morphologically similar to rhyolite-dacite domes on earth. Flows and sinuous channels with lengths of many hundreds of kilometers suggest that extremely high effusion rates or very fluid magmas (perhaps komantiites) may be present. Volcanism is evident in various tectonic settings (coronae, linear extensional and compressional zones, mountain belts, upland rises, highland plateaus, and tesserae). Volcanic resurfacing rates appear to be low (less than 2 cu km/yr) but the significance of dike formation and intrusions, and the mode of crustal formation and loss remain to be established.

Description: Mars exploration crews will be exposed to such high radiation dosages in route from earth that as to sharply reduce the allowable dose they should receive while on the Martian surface. An account is presently given of the possibility of using Martian regolith as crew shielding to maintain very low short-term dose limits. NASA-Langley's nucleon and heavy-ion transport computer codes are used to predict the propagation and interaction of the free-space nucleons and heavy ions through the Martian atmosphere and then through various thicknesses of regolith. The largest reduction in dose occurs in the first 15-20 cm of regolith material.

Description: Accelerator mass spectrometry is utilized to determine the half-life of Ca-41 from the decrease of its concentration with terrestrial age in five Antarctic meteorites and a recent fall. The meteorites were selected on the basis of their Cl-36 concentrations, which showed a span of terrestrial ages of about 600 ka, and on the basis of other cosmogenic nuclide concentrations which indicated that the meteorites had small preatmospheric sizes, and sufficiently long irradiation times in space that the concentrations of Ca-41 and Cl-36 were in secular equilibrium prior to the meteorites' fall to earth. The half-life of Ca-41 is determined at 103 + or - 7 ka. Topics discussed include the effects of undersaturation (short exposure time in space), shielding (the samples are from the interior of a large meteorite), and weathering on the cosmogenic nuclide concentrations in meteorites.

Description: Analytic expressions are derived to model the surface topography and the normal to the surface of Phobos. The analytic expressions are comprised of a spherical harmonic expansion for the global figure of Phobos, augmented by addition terms for the large crater Stickney and other craters. Over 300 craters were measured in more than 100 Viking Orbiter images to produce the model. In general, the largest craters were measured since they have a significant effect on topography. The topographic model derived has a global spatial and topographic accuracy ranging from about 100 m in areas having the highest resolution and convergent, stereo coverage, up to 500 m in the poorest areas.

Description: From February to March 1989 the Phobos 2 spacecraft took 37 TV images of Phobos at a distance of 190-1100 km. These images complement Mariner-9 and Viking data by providing higher-resolution coverage of a laarge region west of the crater Stickney (40-160 deg W) and by providing disk-resolved measurements of surface brightness at a greater range of wavelengths and additional phase angles. These images have supported updated mapping and characterization of large craters and grooves, and have provided additional observations of craters' and grooves' bright rims. Variations in surface visible/near-infrared color ratio of almost a factor of 2 have been recognized; these variations appear to be associated with the ejecta of specific large impact craters. Updated determinations of satellite mass and volume allow calculation of a more accurate value of bulk density, 1.90 + or - 0.1 g/cu cm. This is sigificantly lower than the density of meteoritic analogs to Phobos' surface, suggesting a porous interior perhaps containing interstitial ice.

Description: The unsteady incompressible Navier-Stokes equations have been accurately solved for the laminar flow past a circular cylinder in the Reynolds number range 50-200. A direct elliptic solver called the SEVP is used to rapidly advance the streamfunction in time, facilitating the overall convergence to the fully periodic or quasi-steady state. A new integral-series method is developed for the far-field streamfunction condition on a finite two-dimensional computational domain. The use of fourth-order Hermitain relations for the convection terms in the conservation-form vorticity transport equation has also contributed to the good comparison of the present results with the earlier experimental data. The vortex-shedding patterns visualized by the experimentalist are numerically reproduced here in the given Reynolds number range. Discussions that may be helpful in interpreting the behavior of the shedding frequency are presented in the main text.

Description: The morphology and temporal variability of the Venus ionosphere are characterized, reviewing the results of recent theoretical investigations, observations, and in situ measurements, especially by the Pioneer Venus Orbiter (PVO). Consideration is given to the Pioneer mission and orbit evolution, early radio occultation profiles of ionospheric N(e), the mean structure and thermal balance of the ionosphere, the ion composition and its dawn-dusk asymmetry, the small-scale spatial structure on the nightside, latitudinal and seasonal variations, solar-cycle effects, suprathermal electrons and superthermal ions, and the global configuration and stability of the ionopause. Data from a single PVO passage through the ionotail are discussed in detail, examining the implications for ion escape and the solar-cycle and short-term variability. The differences among the terrestrial, Martian, and Venusian ionospheres are outlined; the PVO data base is described; and the sources of measurement error are indicated.

Description: Internal heat flow from radioactive decay in Triton's interior along with absorbed thermal energy from Neptune total 5 to 20 percent of the isolation absorbed by Triton, thus comprising a significant fraction of Triton's surface energy balance. These additional energy inputs can raise Triton's surface temperature between about 0.5 and 1.5 K above that possible with absorbed sunlight alone, resulting in an increase of about a factor of about 1.5 to 2.5 in Triton's basal atmospheric pressure. If Triton's internal heat flow is concentrated in some areas, as is likely, local effects such as enhanced sublimation with subsequent modification of albedo could be quite large. Furthermore, indications of recent global albedo change on Triton suggest that Triton's surface temperature and pressure may not now be in steady state, further suggesting that atmospheric pressure on Triton was as much as ten times higher in the recent past.

Description: A turbulence model to compute bulk properties is presented in which a prescribed source function is used to provide the rate of energy input into the turbulent flow, and the EDQNM model is used to treat the nonlinear transfer in the Navier-Stokes equations. The predictions of the model are tested against (1) the measured Nusselt number versus Rayleigh number relation in turbulent laboratory convection, and (2) the measured bulk kinetic energies and dissipation rates in turbulent channel flow for Reynolds numbers Re = 12,300 and 30,800. In addition, a sensitivity study is performed with respect to the choice of Kolmogorov constant. Generally, the predictions of the model are in reasonable accord with the available experimental data.

Description: Model spectra have been constructed for the CH4 lines of the polar regions of Jupiter in such as way as to encompass the hydrogen pressure-induced opacity and haze layers revealed by high-resolution spectroscopy. Comparisons between observation and model spectra indicate that, for optically thin haze models, the polar haze is primarily located between the 70 and 5 mbar pressure levels; for opaque haze models, the top of the haze layers lies in the vicinity of the 15-mbar level. It is concluded that the polar haze is widely distributed over both polar regions and exhibits no strong correlation with auroral activities.

Description: Since the vapor phase composition of Titan's methane-nitrogen lower atmosphere is uniquely determined as a function of the Gibbs phase rule, these data are presently computed via integration of the Gibbs-Duhem equation. The thermodynamic consistency of published measurements and calculations of the vapor phase composition is then examined, and the saturated mole fraction of gaseous methane is computed as a function of altitude up to the 700-mbar level. The mole fraction is found to lie approximately halfway between that computed from Raoult's law, for a gas in equilibrium with an ideal solution of liquid nitrogen and methane, and that for a gas in equilibrium with pure liquid methane.

Description: Exospheric temperatures, T(ex), inferred from in situ measurements of neutral densities by the Pioneer Venus Orbiter have been analyzed to examine any possible dependence on solar wind dynamic pressure, P(sw). No linear relationship was found between P(sw) and T(ex). However, on one occasion when P(sw) was unusually high for a prolonged period, an unmistakable temperature enhancement was observed. Minor enhancement of T(ex) were also observed in few other cases, with sustained but moderately high P(sw).

Description: A 0.3 - 0.6 micron UV-visible spectrophotometer and a 5 - 50 micron radiometer in the KRFM experiment on Phobos 2 measured two groundtracks in the equatorial region of Phobos. Preliminary results indicate that three surface units can be recognized on the basis of differing UV-visible spectral reflectance properties. One of the units is most comparable spectrally to optically darkened mafic material, and a second is comparable either to anhydrous carbonaceous chondrite or to blackened mafic material. Spectral properties of the third unit do not resemble those of known meteorite types. Brightness temperatures measured by the radiometer are consistent with a typical surface thermal inertia of 1 - 3 x 10 to the -3 cal/(sq cm deg s exp 1/2), as suggested by previous investigations, implying a lunar-like regolith texture. At least one area of possibly higher thermal inertia has been tentatively identified, where a large degraded crater is crossed by several grooves. These results indicate significant lateral heterogeneity in the optical and textural properties of Phobos' surface.

Description: Numerical solutions of the Navier-Stokes equations governing shear flow in a trapezoidal enclosure have been obtained at fine grid resolutions using an efficient multigrid calculation procedure. The effects of flow Reynolds number and wall boundary conditions are studied in detail. Interesting flow patterns are observed to develop, especially at high Reynolds numbers.

Description: Reference is made to the work of Shah (1979) which demonstrated the possibility of partially integrating the radiative equations analytically to obtain an 'exact' solution. Shah's solution was given as a double integration of the modified Bessel function of order zero. Here, it is shown that the 'exact' solution for a rectangular region radiating to cold black walls can be conveniently derived, and expressed in simple form, by using an integral function, Sn, analogous to the exponential integral function appearing in plane-layer solutions.

Description: The modeling of the dynamics of particle collisions within planetary rings is discussed. Particles in the rings collide with one another because they have small random motions in addition to their orbital velocity. The orbital speed is roughly 10 km/s, while the random motions have an average speed of about a tenth of a millimeter per second. As a result, the particle collisions are very gentle. Numerical analysis and simulation of the ring dynamics, performed with the aid of a supercomputer, is outlined.

Description: The kinetic theory of planetary rings developed by Araki and Tremaine (1986) and Araki (1988) is extended and refined, with a focus on the implications of finite particle size: (1) nonlocal collisions and (2) finite filling factors. Consideration is given to the derivation of the equations for the local steady state, the low-optical-depth limit, and the steady state at finite filling factors (including the effects of collision inelasticity, spin degrees of freedom, and self-gravity). Numerical results are presented in extensive graphs and characterized in detail. The importance of distinguishing effects (1) and (2) at low optical depths is stressed, and the existence of vertical density profiles with layered structures at high filling factors is demonstrated.

Description: The surface of the Jovian satellite Europa is characterized on the basis of an analysis of ground photoelectric photometry at 470 and 550 nm and Voyager images. The data are presented in extensive tables and graphs and discussed in detail. At 550 nm, Europa has single-scattering albedo 0.964, opposition-effect amplitude 0.5, opposition-effect width 0.0016, double-lobed Henyey-Greenstein factors b = -0.429 and c = 0.113, and mean roughness angle 10 deg (much lower than on other solar-system objects). From the small roughness and the 96-percent porosity implied by the narrow opposition peak, it is concluded that the surface was formed mainly by endogenic processes. It is also noted that only one of three observational criteria for preferential ion bombardment of the trailing hemisphere are met in Europa.

Description: The surfaces of the Uranian satellites Ariel, Miranda, Oberon, Titania, and Umbriel are characterized on the basis of Voyager observations. Tables of spectrophotometric data and maps of normal reflectances, green/violet ratios, and possible geological formations are presented and discussed in detail. Variations in albedo are found to be associated with impact features, and it is inferred from color differences that the upper surface of Ariel contains a higher proportion of redder material (tentatively identified as accreted low-albedo meteoritic dust) than those of the other moons.

Description: Direct numerical simulations of unsteady channel flow were performed at low to moderate Reynolds numbers on computational boxes chosen small enough so that the flow consists of a doubly periodic (in x and z) array of identical structures. The goal is to isolate the basic flow unit, to study its morphology and dynamics, and to evaluate its contribution to turbulence in fully developed channels. For boxes wider than approximately 100 wall units in the spanwise direction, the flow is turbulent, and the low-order turbulence statistics are in good agreement with experiments in the near-wall region. For a narrow range of widths below that threshold, the flow near only one wall remains turbulent, but its statistics are still in fairly good agreement with experimental data when scaled with the local wall stress. For narrower boxes only laminar solutions are found. In all cases, the elementary box contains a single low-velocity streak, consisting of a longitudinal strip on which a thin layer of spanwise vorticity is lifted away from the wall.

Description: Near-wall flow structures in turbulent shear flows are analyzed, with particular emphasis on the study of their space-time evolution and connection to turbulence production. The results are obtained from investigation of a database generated from direct numerical simulation of turbulent channel flow at a Reynolds number of 180 based on half-channel width and friction velocity. New light is shed on problems associated with conditional sampling techniques, together with methods to improve these techniques, for use both in physical and numerical experiments. The results clearly indicate that earlier conceptual models of the processes associated with near-wall turbulence production, based on flow visualization and probe measurements need to be modified. For instance, the development of asymmetry in the spanwise direction seems to be an important element in the evolution of near-wall structures in general, and for shear layers in particular. The inhibition of spanwise motion of the near-wall streaky pattern may be the primary reason for the ability of small longitudinal riblets to reduce turbulent skin friction below the value for a flat surface.

Description: The effects of a nonlinear-nonequilibrium-viscous critical layer on the spatial evolution of subsonic and supersonic instability modes on a compressible free shear layer is considered. It is shown that the instability wave amplitude is governed by an integrodifferential equation with cubic-type nonlinearity. Numerical and asymptotic solutions to this equation show that the amplitude either ends in a singularity at a finite downstream distance or reaches an equilibrium value, depending on the Prandtl number, viscosity law, viscous parameter and a real parameter which is determined by the linear inviscid stability theory. A necessary condition for the existence of the equilibrium solution is derived, and whether or not this condition is met is determined numerically for a wide range of physical parameters including both subsonic and supersonic disturbances. it is found that no equilibrium solution exists for the subsonic modes unless the temperature ratio of the low-to-high-speed streams exceeds a critical value, while equilibrium solutions for the most rapidly growing supersonic mode exist over most of the parameter range examined.

Description: Novel approaches to the human exploration of Mars are considered with emphasis on a space suit design, extraterrestrial surface mobility, and water supply. A possible way of transporting personnel on the surface of Mars uses a suborbital rocket that will hop from one site to the next, refuelling each time it lands and giving the Martian explorers effective global mobility. Telepresence could be used to avoid limiting the people on Mars to a small exploration area as a result of a lack of transportation infrastructure. Drawings and photographs are included.

Description: A new type of meteoritic material, intermediate in size between meteorites and interplanetary dust particles (IDPs), is described. Melting and filtering of about 100 tons of blue ice near Cap Prudhomme, Antarctica, yielded 7500 or more irregular, friable particles and about 1500 melted spherules, about 100 microns in size, both showing a 'chondritic' composition suggestive of an extraterrestrial origin. Analyzed irregular particles appear to be unmelted and have similarities with the fine-grained matrix of primitive carbonaceous chondrites, but are extremely diverse in composition. Isotopic analysis of trapped neon confirms an extraterrestrial origin for 16 of 47 irregular particles and 2 of 19 spherules studied and strongly suggests that they were exposed in space as micrometeoroids. These large Antarctic micrometeorites constitute a new family, or at least a new population, of solar system objects, in a mass range corresponding to the bulk of extraterrestrial material accreted by the earth today.

Description: Research throughout the Caribbean suggests that the geophysical anomalies in the Yucatan first noted by Penfield and Camargo (1981) and called the Chicxulub crater could be the site of the impact purported to have caused the K/T extinctions. A semicircular ring of sink holes, known locally as cenotes, which correlates with the geophysical anomalies has been identified, and it is argued that the origin of the cenote ring is related to postimpact subsidence of the Chicxulub crater rim. If there is indeed a crater, the region within the cenote ring corresponds to its floor and the crater rim diameter is probably larger than 200 km. If confirmed as a site of impact, the Chicxulub crater would be the largest terrestrial impact crater known, which is consistent with the uniqueness of the K/T global catastrophe.

Description: An analysis of nonlinear wave-wave interactions in turbulent jets based on the integrated energy of each scale of motion in a cross section of the jet shows that two frequency components in the axisymmetric mode can interact with other mackground frequencies in that mode, thereby amplifying many other frequencies. The present computations produce several features consistent with experimental observations on two-frequency excitation, such as the dependence of the interaction on the initial phase differences between the waves, the enhancement of the momentum thickness under multifrequency forcing, and the increase in background turbulence under forcing. Mixing enhancement is found to be due to turbulence enhancement, rather than the simple amplification of forced-wave components.