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: The use of satellite-derived imagery for measuring subresolution horizontal terrain displacments associated with present-day earthquakes is discussed with reference to data from the French SPOT satellite whose sensor array provides 10-m panchromatic imagery. The measured terrain displacements can be up to several meters, but usually no more than 6-8 m even for major earthquakes. The general approach is to spatially match the after image to the before image at each point on a half-kilometer grid by iteratively interpolating one and testing its correlation with the other. The discussion covers the basic algorithm, results of initial tests, error types and limitations, and future work.

Description: A new approach to stratigraphic analysis is described which uses photogeologic and spectral interpretation of multispectral remote sensing data combined with topographic information to determine the attitude, thickness, and lithology of strata exposed at the surface. The new stratigraphic procedure is illustrated by examples in the literature. The published results demonstrate the potential of spectral stratigraphy for mapping strata, determining dip and strike, measuring and correlating stratigraphic sequences, defining lithofacies, mapping biofacies, and interpreting geological structures.

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: Column ozone measurements obtained from a Spin-scan Ozone Imager (SOI) aboard the Dynamics Explorer I (DE-I) satellite are discussed. The orbital geometry has many features similar to the Molniya satellite orbits. The ozone imagery obtained from DE-I are unique in being synoptic views of a region of continental scale for a period of a few hours. A brief description of the instrument, data, and orbit is followed by a discussion of some applications of DE-I data and the potential usefulness of an ozone imager aboard a satellite in a Molniya-type orbit.

Description: The Colima Volcanic Complex at the western end of the Mexican Volcanic Belt is the most active andesitic volcano in Mexico. Short-wavelength infrared data from the Landsat Thematic Mapper satellite were used to determine the temperature and fractional area of radiant picture elements for two January data acquisitions in 1985 and 1986. The 1986 data showed four 28.5 m by 28.5 m pixels (picture elements) whose hot subpixel components had temperatures ranging from 511-774 C and areas of 1.8-13 sq m. The 1985 data had no radiating areas above background temperatures. Ground observations and measurements in November 1985 and February 1986 reported the presence of hot fumaroles at the summit with temperatures of 135-895 C. This study demonstrates the utility of satellite data for monitoring volcanic activity.

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: Present understanding of the earth's atmosphere is briefly reviewed. The structure and composition of the atmosphere are described. The origin of the atmosphere and the factors involved in global atmospheric change are addressed.

Description: Ground-based and airborne measurements of biomass-burning smoke particle optical properties, obtained with a view to aerosol-absorption properties, are presented as a function of time and atmospheric height. The wavelength dependence of the optical thickness can be explained by a log-normal size distribution, with particles' effective radius varying between 0.1 and 0.2 microns. The strong correlation noted between aerosol particle profile and CO profile indicates that smoke particulates constitute a good tracer for emission trace gases from tropical biomass burning.

Description: Results are presented from airborne measurements of aerosol mass loading, size distribution, and elemental composition obtained in a smoke plume from a fire at a Florida wildlife refuge. It is found that there was a high concentration of carbon-containing aerosols and salt crystals in the 0.1-0.2 micron size range; the composition and morphology of the aerosols differed with size. Aerosol mass concentrations are used in conjunction with those obtained for CO2 in order to calculate ratios of aerosol and CO2. The ratios are noted to be higher for the smoldering phase of the fire than for its flaming phase.

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: Results derived from a spectroscopic analysis of the nu5 band R19 transition of C2H2 observed in solar spectra recorded at the Jungfraujoch station, Switzerland, between June 1986 and April 1991 are reported. A least-squares sine fit to the data reveals a seasonal variability with an amplitude of about +/-40 percent of the mean; the maximum occurs during mid-winter and the minimum in summer. In general, the observed seasonal observations are comparable with those derived from airborne in-situ measurements and those reported from ground-based stations.

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: The objectives of this study are the following: (1) to assess the quality of the sea level observed by altimetry on a global scale, especially in tropical regions where the atmospheric effects are the most critical, and to prepare the altimeter data for their assimilation into tropical Oceanic General Circulation Models; and (2) to validate and calibrate the altimeter-derived sea state parameters, to assimilate these parameters into numerical models, and to estimate the altimeter's sea state and electromagnetic biases.

Description: In the proposed research, TOPEX/POSEIDON altimeter data will be used with Geosat and European Remote Sensing satellite (ERS-1) data to compute global, region, and local oceanic geoid surfaces. These observations will then be analyzed to conduct geophysical studies relative to the structure of the oceanic lithosphere and mantle.

Description: The primary objective of our research program is the application of TOPEX/POSEIDON oceanic data to the study of various geophysical observations, thereby improving our knowledge of the Earth's interior. Observations of the Earth's solid tides, rotation, time-dependent gravity, and crustal deformation provide potentially useful constraints on properties of the solid Earth and its liquid core, but all these parameters are also sensitive to oceanic effects. Depending on the type of observation, those effects are usually either removed by modeling (often using particularly simple assumptions about the dynamics of the ocean) or ignored entirely. Useful oceanic data are rarely available. In many cases, the uncertainty in the oceanic correction is the main source of error in the geophysical interpretation of the results. We describe, below several of these geophysical observations, how their interpretation in terms of the Earth's interior is affected by the oceans, and how data from the TOPEX/POSEIDON mission can be used to improve this situation.

Description: The goal of the proposal is to determine the present motion of the main tectonic plates from the Doppler data of the Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS) orbitography system, which includes in its final configuration about 50 tracking stations with a world-wide distribution.

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 recently reprocessed Nimbus-7 TOMS total ozone data from 1979 to 1989 are analyzed to assess the global impact of the El Chichon volcanic eruption on stratospheric total ozone. An apparent decrease in total ozone of 5 to 6 percent is present during the winter of 1982-1983 following the eruption of El Chichon. A regression analysis of total ozone with the equatorial zonal wind at 30 mb indicates that response to the quasi-biennial oscillation can explain much of the observed ozone anomaly, and that the total ozone decrease which can be attributed to El Chichon is at most 2 to 4 percent. This study also suggests that the interannual variability caused by the quasi-biennial oscillation and planetary wave activity may introduce apparent seasonal trends in total ozone and temperature in the lower stratosphere, particularly if the length of the data record is not very long. Such trends may affect the assessment of total ozone changes caused by chemical perturbations.

Description: The present study compares the diurnal variation of NO2 measured near 30 km by the BLISS in situ laser spectrometer with calculations from a photochemical model that includes a detailed description of multiple scattering. Even better agreement is found between the data and the model, both at sunset and during the day. The conclusions of an earlier study that the high-resolution in situ measurements of NO2 facilitated validation of the understanding of the diurnal chemistry of NO2 are confirmed.

Description: Variations in crustal magnetization along a seismic section across the Archean Yilgarn block of Western Australia inferred from Magsat data are interpreted as a subtle thermal effect arising from variations in depth to the Curie isotherm. The isotherm lies deep within the mantle of the eastern part of the province, but transects the crust-mantle transition and rises well into the crust on the western side. The model is consistent with heat flow variations along the section line. The mean crustal magnetization implied by the model is approximately 2 A/m. The temperature variation implied by the model is consistent with the hypothesis that the crust-mantle transition seen seismically corresponds to the mafic granulite-eclogite phase transition within a zone of igneous crustal underplating.

Description: Results are presented of solar irradiance measurements in the spectral range 160-400 nm at approximately 0.15-0.20-nm intervals and at 1-nm resolution performed continually since November 1978. Solar irradiance data from the Nimbus-7 SBUV satellite instrument, the SBUV/2 instruments on the NOAA-9 and NOAA-11 satellites, and the October 1989 flight of the Shuttle SBUV instrument are presented and compared. Uncertainties in the instruments' absolute and long-term radiometric calibrations, which vary among the four instruments, are discussed. Comparisons of the initial solar spectra from the four instruments show agreement to within approximately 10 percent, with spectral biases on the order of +/-4 percent. Irradiances measured by the two NOAA instruments and SSBUV agree to within about 5 percent overall from 270 to 360 nm, with spectral biases on the order of about +/-2 percent. The Nimbus-7 SBUV irradiances are an additional 5-10 percent lower in this region than those measured by the other three instruments.

Description: Stratospheric concentrations of HO2, H2O2, and OH have been retrieved simultaneously from the far-infrared emission spectra obtained with a balloon-borne Fourier transform spectrometer in June 1983 at 32 deg N latitude. Retrieved concentrations of HO2 and H2O2 are reported, along with vertical distributions of OH which were reported in an earlier paper for the afternoon, sunset, and nighttime periods for altitudes from 26 to 38 km. HO2 distributions are obtained with uncertainties that are about the same as OH for the same vertical range and for the afternoon and sunset periods. H2O2 concentration is obtained at an altitude of 30 km for the period that covers afternoon and sunset hours. The retrieved concentrations of these HO(x) species agree well with other individually measured results and the steady state photochemical predictions. The ratio HO2/OH at around 32 km seems to increase from the afternoon period to the sunset period.

Description: An aerosol size distribution model for the stratosphere is inferred based on 5 years of Stratospheric Aerosol and Gas Experiment (SAGE) II measurements of multispectral aerosol and water vapor extinction. The SAGE II aerosol and water vapor extinction data strongly suggest that there is a critical particle radius below which there is a relatively weak dependence of particle number density with size and above which there are few, if any, particles. A segmented power law model, as a simple representation of this dependence, is used in theoretical calculations and intercomparisons with a variety of aerosol measurements including dustsondes, longwave lidar, and wire impactors and shows a consistently good agreement.

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: Several characteristic transition heights can be found in ionospheric ion-composition profiles. These are the boundaries between the regions dominated by light ions, atomic oxygen ions, molecular ions, and cluster ions, respectively. Ion-composition modeling can benefit from the use of these transition heights. Special emphasis is given to the 'upper' transition height H sub T (light ions to atomic oxygen) and to the 'lower' transition height h sub t (atomic oxygen to molecular ions). Transition-height models deduced from rocket and satellite measurements are compared with the heights predicted by the International Reference Ionosphere. Considerable discrepancies are found between the different models as well as between the models and independent measurements.

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: 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: The application of van Leer's scheme, a monotonic, upstream-biased differencing scheme, to three-dimensional constituent transport calculations is shown. The major disadvantage of the scheme is shown to be a self-limiting diffusion. A major advantage of the scheme is shown to be its ability to maintain constituent correlations. The scheme is adapted for a spherical coordinate system with a hybrid sigma-pressure coordinate in the vertical. Special consideration is given to cross-polar flow. The vertical wind calculation is shown to be extremely sensitive to the method of calculating the divergence. This sensitivity implies that a vertical wind formulation consistent with the transport scheme is essential for accurate transport calculations. The computational savings of the time-splitting method used to solve this equation are shown. Finally, the capabilities of this scheme are illustrated by an ozone transport and chemistry model simulation.

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 unified approach for combining active and passive microwave measurements for remote sensing applications is described. A synergic inversion technique has been developed and applied to the retrieval of geophysical parameters of the ocean surface and of the atmosphere. It is based on the combination of radiometric and radar measurements at the electromagnetic and cell level and not only on the correction of radar measurements by radiometric measurements, or conversely. Such a combination is performed through a common quantity: the bistatic scattering coefficient of the observed surface. This is used in a direct model to simulate combined measurements from active and passive sensors. It requires a rather complete and accurate calculation of the scattering of microwaves by the rough sea surface. Comparisons with experimental models derived for scatterometers and radar altimeters or with data derived from radiometric measurements show that the results of this unified approach are consistent with those from separate microwave sensors, at least for a foam-free sea surface.

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: Ozone mixing ratios were measured by ozonesondes on board balloons launched from Esrange, near Kiruna, Sweden (68 deg N, 20 deg E) from January 11 to February 9, 1990. The data obtained prior to a sudden warming on February 7, 1990 show that at potential temperatures between 460 and 640 K, the ozone mixing ratio just inside the polar vortex was systematically smaller than that outside, the largest difference being 29 percent at around 525 K. The ozone mixing ratio at 525 K inside the vortex decreased at a rate of about 1.5 percent per day between January 26 and February 4. The temperatures simultaneously observed were quite often low enough to allow for formation of nitric acid trihydrate particles around this altitude. Depletion of ozone due to highly perturbed chemical conditions in late January and early February is strongly suggested.

Description: A new cryogenic collection system has been flown on board a balloon gondola to obtain separate samples of ozone and carbon dioxide without entrapping major atmospheric gases. Precision laboratory isotopic analysis of CO2 samples collected between 26 and 35.5 km show a mass-independent enrichment in both O-17 and O-18 of about 11 per mil above tropospheric values. Ozone enrichment in its heavy isotopes was 9 to 16 percent in O3-50 and 8 to 11 percent in O3-49, respectively (Schueler et al., 1990). A mechanism to explain the isotope enrichment in CO2 has been recently proposed by Yung et al. (1991). The model is based on the isotope exchange between CO2 and O3 via O(1D), resulting in a transfer of the ozone isotope enrichment to carbon dioxide. Predicted enrichment and measured values agree well.

Description: The 1990 Antarctic ozone hole matched the record 1987 ozone hole in depth, duration, and area. During the formation phase of the hole (August), total ozone values were the lowest yet recorded. The decline rate approximately matched the record 1987 decline and reached a minimum of 125 Dobson Units on October 4, 1990. October total ozone averages were marginally higher that 1987. As during 1987, the 1990 total ozone values within the hole slowly and steadily increased during the mid-October through November period. The ozone hole breakup was the latest yet recorded (early December), with low ozone values persisting over the pole through December, setting a record low for December average polar ozone. Temperatures were near average during the early spring, but were below normal for the late spring. Temperatures in the early spring of 1990 were substantially warmer than those observed in the early spring of 1987.

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: 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: Geodetic measurements with Rogue GPS receivers from sites in the California Permanent GPS geodetic Array (PGGA) have been analyzed using the GIPSY orbit-determination and baseline-estimation software. Based on an unbiased selection of 23 daily measurements spanning 8 months, the LF contributions to the long-term repeatabilities of baseline measurements are approximately 5, 3, and 8 mm for the east, north, and vertical components. Short-term contributions to the long-term repeatabilities were evaluated by examining data from the week of October 21, 1990, which showed the lowest short-term scatter. For this week, daily repeatabilities of 2-3 mm in the horizontal and 4 mm in the vertical have been achieved for the 172-km JPL-Pinyon baseline, consistent with carrier phase date noise of about 6 mm. High quality (less than about 5 mm) repeatabilities have been achieved for all components of the other baselines as well.

Description: Thermal ion energy distribution functions and local electric and magnetic fields were directly measured for the first time in the ionospheric E region. Measured ion distribution functions were fitted to shifted Maxwellian distributions, and their resulting ion drift velocities were compared with E x B/B-squared velocities from the double-probe electric field observations. The results show that the ion drift direction rotates with respect to the local electric field direction and that the ratio of the magnitudes of the ion velocity to the E x B/B-squared velocity decreases with decreasing altitudes. Using these observations, the quiet time ion-neutral collison frequencies and neutral wind velocities were estimated and found to be consistent with theoretical estimates. However, significant discrepancies between observations and theory are found in the disturbed E region near auroral particle precipitation regions. These data indicate that the auroral atmosphere is significantly perturbed due to Joule as well as particle heating effects.

Description: Magnetic properties of the high-latitude tail boundary are examined with IMP 8 magnetometer data. The high-latitude tail boundary separates the tail lobe from the magnetosheath. Magnetic fields are stable in the tail lobe, but very irregular in the magnetosheath. Boundary crossings are marked by the rotation of magnetic components parallel to the boundary plane. The magnetic component normal to the boundary, if any, is very small in comparison to this rotational change. Despite large magnetic fluctuations in the magnetosheath, the magnetosheath-side field orientation is consistent with the draping of the IMF against the magnetotail. The boundary current has a component parallel to the lobe field (tail-aligned current), as well as a circumferential component. The IMF orientation controls the sheath-side B(Y), while the lobe-side field has a more rigid configuration flaring antisunward.

Description: Electron densities from the Langmuir probes on the Atmospheric Explorer C and Dynamics Explorer 2 are used for analyzing the behavior of the high-altitude night-side F region polar hole as a function of solar and magnetic activity and of universal time (UT). The polar region of invariant latitude from 70 deg to 80 deg and MLT from 22 to 03 hours is examined. The strongest dependencies are observed in F10.7 and UT; a strong hemispherical difference due to the offset of the magnetic poles from the earth's rotation axis is observed in the UT dependence of the ionization hole. A seasonal variation in the dependence of ion density on solar flux is indicated, and an overall asymmetry in the density level between hemispheres is revealed, with the winter-hole density about a factor of 10 greater in the north than in the south.

Description: This paper reexamines arguments formulated and restated by Heikkila to the effect that a magnetosheath plasma cloud having 'excess' momentum which impinges on the dayside magnetopause boundary is able to flow continuously through the boundary onto both open and closed flux tubes in the interior. It is shown that the argument used to arrive at this conclusion is not correct. The error in the argument concerns the nature of the flow which is associated with the induction electric field produced by the perturbed current layer, which was assumed by Heikkila to be such as to keep the plasma jet just moving with the boundary. Heikkila's argument, correctly applied, does not lead to 'impulsive transport' of plasma through the magnetopause.

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: Laboratory experiments document that liquid iron reacts chemically with silicates at high pressures (above 2.4 x 10 to the 10th Pa) and temperatures. In particular, (Mg,Fe)SiO3 perovskite, the most abundant mineral of earth's lower mantle, is expected to react with liquid iron to produce metallic alloys (FeO and FeSi) and nonmetallic silicates (SiO2 stishovite and MgSiO3 perovskite) at the pressures of the core-mantle boundary, 14 x 10 to the 10th Pa. The experimental observations, in conjunction with seismological data, suggest that the lowermost 200 to 300 km of earth's mantle, the D-double-prime layer, may be an extremely heterogeneous region as a result of chemical reactions between the silicate mantle and the liquid iron alloy of earth's core. The combined thermal-chemical-electrical boundary layer resulting from such reactions offers a plausible explanation for the complex behavior of seismic waves near the core-mantle boundary and could influence earth's magnetic field observed at the surface.

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: For more than 50 years, observations of earthquakes to depths of 100 to 650 kilometers inside earth have been enigmatic: at these depths, rocks are expected to deform by ductile flow rather than brittle fracturing or frictional sliding on fault surfaces. Laboratory experiments and detailed calculations of the pressures and temperatures in seismically active subduction zones indicate that this deep-focus seismicity could originate from dehydration and high-pressure structural instabilities occurring in the hydrated part of the lithosphere that sinks into the upper mantle. Thus, seismologists may be mapping the recirculation of water from the oceans back into the deep interior of the planet.

Description: Measurements of water diffusivity in a basaltic liquid are reported. The concentration-dependent total water diffusivities in the basaltic melt at 1300-1500 C are 30-50 times as large as those in rhyolitic melts and are greater than the total CO2 diffusivity in basaltic melts, contrary to previous expectations. These results suggest that diffusive fractionation would increase the ratio of water to CO2 in growing bubbles relative to equilibrium partitioning and decrease the ratio in interface melts near an advancing anhydrous phenocryst.

Description: An airborne continuous-wave (CW) focused CO2 Doppler lidar and a ground-based pulsed CO2 Doppler lidar were to obtain seven pairs of comparative measurements of tropospheric aerosol backscatter profiles at 10.6-micron wavelength, near Denver, Colorado, during a 20-day period in July 1982. In regions of uniform backscatter, the two lidars show good agreement, with differences usually less than about 50 percent near 8-km altitude and less than a factor of 2 or 3 elsewhere but with the pulsed lidar often lower than the CW lidar. Near sharp backscatter gradients, the two lidars show poorer agreement, with the pulsed lidar usually higher than the CW lidar. Most discrepancies arise from a combination of atmospheric factors and instrument factors, particularly small-scale areal and temporal backscatter heterogeneity above the planetary boundary layer, unusual large-scale vertical backscatter structure in the upper troposphere and lower stratosphere, and differences in the spatial resolution, detection threshold, and noise estimation for the two lidars.

Description: Since 1981 the Royal Signals and Radar Establishment and the Royal Aircraft Establishment, United Kindom, have made vertical and horizontal sounding measurements of aerosol backscatter coefficients at 10.6 microns, using an airborne continuous-wave-focused CO2 Doppler lidar, the Laser True Airspeed System (LATAS). In this paper, the heterodyne signal from the LATAS detector is spectrally analyzed. Then, in conjunction with aircraft flight parameters, the data are processed in a six-stage computer algorithm: set search window, search for peak signal, test peak signal, measure total signal, calculate signal-to-noise ratio, and calculate backscatter coefficient.