ALBERT

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Description: The following are discussed in this viewgraph presentation: Glovebox-2 investigations, study of oscillatory thermocapillary flows and validate physical models on the second U.S. Microgravity Laboratory Mission (USML-2), pool boiling experiment, critical fluid viscosity measurement experiment, Hitchhiker-G typical structural configuration, ZENO-critical fluid light scattering experiment, isothermal dendritic growth experiment, microgravity acceleration environment support, and the space acceleration measurement system (SAMS).

Description: The response of a pair of emulsion drops to the imposition of a uniform electric field is examined. The case studied is that of equal-sized drops whose line of centers is parallel to the axis of the applied field. A new boundary integral solution to the governing equations of the leaky dielectric model is developed; the formulation accounts for the electrostatic and hydrodynamic interactions between the drops, as well as their deformations. Numerical calculations show that, after an initial transient during which the drops primarily deform, the pair drift slowly together due to their electrostatic interactions.

Description: A noninvasive fluorescence recovery after photobleaching (FRAP) technique is under development to measure interfacial transport in two phase systems without disturbing the interface. The concentration profiles of a probe solute are measured in both sides of the interface by argon-ion laser, and the system relaxation is then monitored by a microscope-mounted CCD camera.

Description: The effects of hydrodynamic convection on nucleation and broken chiral symmetry have been investigated for a simple inorganic molecule, sodium chlorate (NaClO3). Our experiment suggests that the symmetry breaking is a result of hydrodynamic amplification of rare nucleation events. The effect is more pronounced when the primary nucleation occurs on the solute-vapor interface, where mixing in the surface sublayer becomes important. The transition from the achiral to the chiral states appears to be smooth as the hydrodynamic parameters, such as flow rate, are varied.

Description: It has been proposed to measure the extensional viscosity function of a non-Newtonian polymer solution in a reduced gravity environment as part of the Advanced Fluid Module. In ground-based extensional measurements, the no-sip boundary condition at solid-fluid interfaces always result in appreciable shear gradients in the test fluid; however the removal of gravitational body forces permits controlled extensional deformation of containerless test samples and the first unambiguous measurements of this kind. Imperative to successful implementation of this experiment is the generation and subsequent deformation of a stable cylindrical column of test fluid. A study of the generation and deformation of liquid bridges demonstrates that Newtonian liquid bridges undergo capillary breakup as anticipated when stretched beyond a critical aspect ratio; non-Newtonian liquid bridges, however, are stabilized by the strain-hardening phenomenon exhibited by these materials. Numerical simulations of Newtonian breakup are compared with experimental results, and show that previous ground-based attempts at measuring the extensional viscosity of Newtonian fluids are of limited accuracy.

Description: The physical scales of velocity, length, time, thermal gradient magnitude likely to be involved in gas-solid multiphase flight experiments are assessed for 1-100 micron particles.

Description: An experiment is being conducted to study the effects of buoyancy on planar stratified flows. A wind tunnel has been designed and constructed to generate planar flows with separate heating for the top and bottom planar air jets emerging from slot nozzles separated by an insulating splitter plate. The objective is to generate planar jet flows with well defined and well controlled velocity and temperature profiles. Magnitudes of velocity and temperature will be varied separately in each flow for both laminar and turbulent flow conditions. Both stably and unstably stratified flows will be studied by changing the temperature distributions in each air stream. This paper reports on the design of the apparatus and initial measurements of velocity and turbulence made by laser Doppler velocimetry.

Description: The two-phase bubbly flow and boiling heat transfer in microgravity represents a substantial challenge to scientists and engineers and yet there is an urgent need to seek fundamental understanding in this area for future spacecraft design and space missions. At Washington State University, we have successfully designed, built and tested a 2.1 second drop tower with an innovation airbag deceleration system. Microgravity boiling experiments performed in our 0.6 second Drop Tower produced data flow visualizations that agree with published results and also provide some new understanding concerning flow boiling and microgravity bubble behavior. On the analytical and numerical work, the edge effects of finite divergent electrode plates on the forces experienced by bubbles were investigated. Boiling in a concentric cylinder microgravity and an electric field was numerically predicted. We also completed a feasibility study for microgravity boiling in an acoustic field.

Description: Design of the two-phase systems which are anticipated to be utilized in future spacecraft thermal management systems requires a knowledge of two-phase flow and heat transfer parameters in reduced gravities. A program has been initiated by NASA to design a two-phase test loop and perform a series of experiments to generate the data for the Critical Heat Flux (CHF) and onset of instability under reduced gravities. In addition to low gravity airplane trajectory testing, the experimental program consists of a set of laboratory tests with vertical upflow and downflow configurations. Modularity is considered in the design of this experiment and the test loop in instrumented to provide data for two-phase pressure drop and flow regime behavior. Since the program is in the final stages of the design and construction task, this article is intended to discuss the phenomena, design approach, and the description of the test loop.

Description: A ground based (1-g) experiment is in progress that will measure the turbidity of a density-matched, binary fluid mixture extremely close to the critical point. By covering the range of reduced temperatures t is equivalent to (T-T(sub c))/T(sub c) from 10(exp -8) to 10(exp -2), the turbidity measurements will allow the critical exponent eta to be determined. No experiment has determined a value of the critical exponent eta, yet its value is significant to theorists in critical phenomena. Interpreting the turbidity correctly is important if future NASA flight experiments use turbidity as an indirect measurement of relative temperature in shuttle experiments on critical phenomena in fluids.

Description: Analytical studies along with ground-based experiments are presently being carried out in connection with thermocapillary phenomena associated with drops and bubbles in a containerless environment. The effort here focuses on the thermal and the fluid phenomena associated with the local heating of acoustically levitated drops, both at 1-g and at low-g. In particular, the Marangoni effect on drops under conditions of local spot-heating and other types of heating are being studied. With the experiments conducted to date, fairly stable acoustic levitation of drops has been achieved and successful flow visualization by light scattering from smoke particles has been carried out. The results include situations with and without heating. As a preliminary qualitative interpretation of these experimental results, we consider the external flow pattern as a superposition of three discrete circulation cells operating on different spatial scales. The observations of the flow fields also indicate the existence of a steady state torque induced by the streaming flows. The theoretical studies have been concentrated on the analysis of streaming flows in a gaseous medium with the presence of a spherical particle undergoing periodic heating. A matched asymptotic analysis was carried out for small parameters derived from approximations in the high frequency range. The heating frequency being 'in tune' with the acoustic frequency results in a nonzero time-averaged thermal field. This leads to a steady heat flow across the equatorial plane of the sphere.

Description: An intuitive delineation along with dimensional considerations and experimental evidences are presented to show that in a general case, the evaporation of a liquid droplet undergoes three regimes through the process. Initially, the heat transfer inside the evaporating droplets is conduction controlled; then, in the second stage, convective heat transfer may take over; finally, the convections subside, and the process returns to conduction controlled mode.

Description: We have studied the microscopic aspects of the spreading of liquid drops on a solid surface by molecular dynamics simulations of coexisting three-phase Lennard-Jones systems of liquid, vapor and solid. We consider both spherically symmetric atoms and chain-like molecules, and a range of interaction strengths. As the attraction between liquid and solid increases we observed a smooth transition in spreading regimes, from partial to complete to terraced wetting. In the terraced case, where distinct monomolecular layers spread with different velocities, the layers are ordered but not solid, with qualitative behavior resembling recent experimental findings, but with interesting differences in the spreading rate.

Description: Nonlinear phenomena associated with the dynamics of free drops and bubbles are investigated analytically, numerically and experimentally. Although newly developed levitation and measurement techniques have been implemented, the full experimental validation of theoretical predictions has been hindered by interfering artifacts associated with levitation in the Earth gravitational field. The low gravity environment of orbital space flight has been shown to provide a more quiescent environment which can be utilized to better match the idealized theoretical conditions. The research effort described in this paper is a closely coupled collaboration between predictive and guiding theoretical activities and a unique experimental program involving the ultrasonic and electrostatic levitation of single droplets and bubbles. The goal is to develop and to validate methods based on nonlinear dynamics for the understanding of the large amplitude oscillatory response of single drops and bubbles to both isotropic and asymmetric pressure stimuli. The first specific area on interest has been the resonant coupling between volume and shape oscillatory modes isolated gas or vapor bubbles in a liquid host. The result of multiple time-scale asymptotic treatment, combined with domain perturbation and bifurcation methods, has been the prediction of resonant and near-resonant coupling between volume and shape modes leading to stable as well as chaotic oscillations. Experimental investigations of the large amplitude shape oscillation modes of centimeter-size single bubbles trapped in water at 1 G and under reduced hydrostatic pressure, have suggested the possibility of a low gravity experiment to study the direct coupling between these low frequency shape modes and the volume pulsation, sound-radiating mode. The second subject of interest has involved numerical modeling, using the boundary integral method, of the large amplitude shape oscillations of charged and uncharged drops in the presence of a static or time-varying electric field. Theoretically predicted non linearity in the resonant frequency of the fundamental quadrupole mode has been verified by the accompanying experimental studies. Additional phenomena such as hysteresis in the frequency response of ultrasoncially levitated droplets in the presence of a time varying electric field, and mode coupling in the oscillations of ultrasonically modulated droplets, have also been uncovered. One of the results of this ground-based research has been the identification and characterization of phenomena strictly associated with the influence of the gravitational field. This has also allowed us to identify the specific requirements for potential microgravity investigations yielding new information not obtainable on Earth.

Description: The objectives are to: (1) Understand the influence in low gravity of flow on interface shape. For example, document and control the influence of axial flow on the Plateau-Rayleigh instability of a liquid bridge; and (2) Extend the ground-based density-matching technique of low gravity simulation to situations with flow; that is, develop Plateau chamber experiments for which flow can be controlled. Containerless containment of liquid by surface tension has broad importance in low gravity. For space vehicles, the behavior of liquid/gas interfaces is crucial to successful liquid management systems. In microgravity science, free interfaces are exploited in various applications. Examples include float-zone crystal growth, phase separation near the critical point of liquid mixtures (spinoidal decomposition) and quenching of miscibility gap molten metal alloys. In some cases, it is desired to stabilize the capillary instability while in others it is desired to induce capillary breakup. In all cases, understanding the stability of interface shape in the presence of liquid motion is central.

Description: The hydrodynamics governing the fluid motions on a microscopic scale near moving contact lines are different from those governing motion far from the contact line. We explore these unique hydrodynamics by detailed measurement of the shape of a fluid meniscus very close to a moving contact line. The validity of present models of the hydrodynamics near moving contact lines as well as the dynamic wetting characteristics of a family of polymer liquids are discussed.

Description: The objective of this paper is to conduct ground-based experiments to measure the onset conditions of oscillatory Marangoni flow in laser-heated silicone oil in a cylindrical container. For a single fluid, experimental data are presented using the aspect ratio and the dynamic Bond number. It is found that for a fixed aspect ratio, there seems to be an asymptotic limit of the dynamic Bond number beyond which no onset of flow oscillation could occur. Experimental results also suggested that there could be a lower limit of the aspect ratio below which there is no onset of oscillatory flow.

Description: The present study demonstrated that calculated thermocapillary convection in a non-cylindrical floating zone can now be compared with measured one, by considering the lens effect of the floating zone. Flow visualization and computer simulation of thermocapillary convection in a silicone oil zone and a molten zone in an NaNO3 rod were conducted. The calculated results agree very well with the measured ones, including the free surface shapes, the solid/melt interface shapes and the velocity fields.

Description: Laboratory and numerical experiments are underway to generate, and subsequently suppress, oscillatory thermocapillary convection in thin layer of silicone oil. The laboratory experiments have succeeded in characterizing the flow state in a limited range of Bond number-Marangoni number space of interest, identifying states of: (1) steady, unicellular, thermocapillary convection; (2) steady, multicellular, thermocapillary convection; and (3) oscillatory thermocapillary convection. Comparisons between experimental results and stability computations for a related basic state will be made.

Description: Nonplanar flow oscillations have been shown to be effective in stabilizing buoyancy-induced Rayleigh-Benard convection. The present study was initiated to see if thermocapillary convection of the Marangoni type might also be stabilized by the same means. When surface deflection can be ignored, significant stabilization occurs. However, when the operating parameters are such that surface deflection is nonnegligible, destabilization can occur, in contrast to Rayleigh-Benard convection. Mechanisms for both stabilization and destabilization are discussed.

Description: Some areas of opportunity for future research in microgravity two-phase flow and heat transfer are pointed out. These satisfy the dual requirements of relevance to current and future needs, and scientific/engineering interest.

Description: The probes described herein, in various configurations, permit the measurement of the volume fraction of two or more fluids flowing through a pipe. Each probe measures the instantaneous relative dielectric constant of the fluid in immediate proximity. As long as separation of the relative dielectric constant of the fluid is possible, several or even many fluids can be measured in the same flow stream. By using multiple probes, the velocity of each fluid can generally be determined as well as the distribution of each constituent in the pipe. The values are determined by statistical computation. There are many potential applications for probes of this type in industry and government. Possible NASA applications include measurements of helium/hydrazine flow during rocket tests at White Sands, liquid/gas flow in hydrogen or oxygen lines in Orbiter engines, and liquid/gaseous Freon flow in zero gravity tests with the KS135 aircraft at JSC. Much interest has been shown recently by the oil industry. In this a good method is needed to measure the fractions of oil, water, and natural gas flowing in a pipeline and the velocity of each. This particular problem involves an extension of what has been developed to date and our plans to solve this problem will be discussed herein.

Description: The objective of this research is to obtain an improved understanding of the behavior of droplets in vaporizing sprays, particularly under conditions typical of those in high pressure rocket sprays. Experiments are conducted in a variety of high pressure, high temperature, optically-accessible flow systems, including one which is capable of operation at pressures up to 70 atm, temperatures up to 600 K, gas velocities up to 30 m/sec and turbulence intensities up to 40 percent. Single droplets, 50 to 500 micron in diameter, are produced by an aerodynamic droplet generator and transversely injected into the flow. Measurements are made of the droplet position, size, velocity and temperature and of the droplet's vapor wake from which droplet drag, dispersion, heating, vaporization and breakup are characterized.

Description: The potential significance of drop fragmentation in sprays and other propulsion-related multiphase flows has been noted in the literature. This has motivated recent experimental and theoretical works to: better understand the fundamentals of physics of drop breakup processes, and develop models of drop fragmentation suitable for use in multiphase flow codes. The works summarized below aim to contribute to both sides of this two-pronged attack.

Description: Results of atomization tests of rocket injection elements are correlated in terms of the test fluid properties and the injection element operating conditions. The type of element modeled has a central, swirling liquid oxidizer flow; and a coaxial, axial gaseous fuel flow. Results of the correlations describe the effect of geometry, fluid properties, and operating conditions on liquid atomization. The ranges of test conditions allow reasonable extrapolation to actual high-pressure engine operating conditions. At such conditions, the swirl element should produce exceptionally finely-atomized and well-distributed spray for high combustion performance. United Technologies Research Center (UTRC) completed a test program to measure fluid spray atomization, in terms of droplet size and droplet size distributions, for swirling-flow injection elements. Element sizes and dimensions covered ranges of dimensions representative of typical cryogenic (O2/H2) liquid rocket engine practice. Relevant dimensions include the element diameter, the tangential-entry metering slot length, and the coastal fuel slot radial gap. Test fluids included water, simulated Jet-A, and Freon 113. These fluids cover ranges of properties which include or approach those of liquid oxygen. Gaseous nitrogen simulated the coaxial gaseous hydrogen fuel flow.

Description: Studies of the interaction of unsteady (oscillatory) flows with the growth of a turbulent boundary layer on a flat plate have primarily dealt with an oscillatory component in the primary flow direction. Past studies of the 2-D flow have shown little or no increase in the time averaged heat transfer. The present paper deals with a steady axial and an oscillatory transverse flow. It is shown that for such flows the temporal variation for both the turbulent skin friction and heat transfer are such as to yield increased time averaged values.

Description: Our Sensitivity-Based Minimal Residual (SBMR) method which is based on our earlier Distributed Minimal Residual (DMR) method allows each component of the solution vector in a system of equations to have its own convergence speed. Our global SBMR method was found to consistently outperform the DMR method while requiring considerably less computer memory. Recently, we have developed and tested a new Line SBMR or LSBMR method and a Time-Step-Scaling (TSS) method that are even more robust and computationally efficient than our global SBMR method, especially on highly clustered computational grids in laminar and turbulent flow computations.

Description: The present paper reports the development of the Least-Squares Finite Element Method (LSFEM) for simulating compressible viscous flows at low Mach numbers in which the incompressible flows pose as an extreme. Conventional approach requires special treatments for low-speed flows calculations: finite difference and finite volume methods are based on the use of the staggered grid or the preconditioning technique; and, finite element methods rely on the mixed method and the operator-splitting method. In this paper, however, we show that such difficulty does not exist for the LSFEM and no special treatment is needed. The LSFEM always leads to a symmetric, positive-definite matrix through which the compressible flow equations can be effectively solved. Two numerical examples are included to demonstrate the method: first, driven cavity flows at various Reynolds numbers; and, buoyancy-driven flows with significant density variation. Both examples are calculated by using full compressible flow equations.

Description: A systematic investigation has been conducted to study the effects of ambient flow conditions (i.e. pressure and velocity) on supercritical droplet gasification in a forced-convective environment. The model is based on the time-dependent conservation equations in axisymmetric coordinates, and accommodates thermodynamic nonidealities and transport anomalies. In addition, an efficient scheme for evaluating thermophysical properties over the entire range of fluid thermodynamic states is established. The analysis allows a thorough examination of droplet behavior during its entire lifetime, including transient gasification, dynamic deformation, and shattering. A parametric study of droplet vaporization rate in terms of ambient pressure and Reynolds number is also conducted.

Description: The objective of the study is to assess the capability of two-equation turbulence models for simulating propulsion-related flowfields. The standard kappa-epsilon model with Chien's low Reynolds number formulation for near-wall effects is used as the baseline turbulence model. Several experimental test cases, representative of rocket combustor internal flowfields, are used to catalog the performance of the baseline model. Specific flowfields considered here include recirculating flow behind a backstep, mixing between coaxial jets and planar shear layers. Since turbulence solutions are notoriously dependent on grid and numerical methodology, the effects of grid refinement and artificial dissipation on numerical accuracy are studied. In the latter instance, computational results obtained with several central-differenced and upwind-based formulations are compared. Based on these results, improved turbulence modes such as enhanced kappa-epsilon models as well as other two-equation formulations (e.g., kappa-omega) are being studied. In addition, validation of swirling and reacting flowfields are also currently underway.

Description: The objective of this paper is to review the efforts in spatial direct numerical simulations for transition modeling. Much recent success has been realized in the development of more efficient numerical algorithms as well as a robust downstream boundary condition. Efforts at explaining complex physical phenomena through the use of simulations are reviewed.

Description: A new procedure closely linking dissociation and exchange reactions in air to the vibrational levels of the diatomic molecules has been implemented in both one- and two-dimensional versions of Direct Simulation Monte Carlo (DSMC) programs. The previous modeling of chemical reactions with DSMC was based on the continuum reaction rates for the various possible reactions. The new method is more closely related to the actual physics of dissociation and is more appropriate to the particle nature of DSMC. Two cases are presented: the relaxation to equilibrium of undissociated air initially at 10,000 K, and the axisymmetric calculation of shuttle forebody heating during reentry at 92.35 km and 7500 m/s. Although reaction rates are not used in determining the dissociations or exchange reactions, the new method produces rates which agree astonishingly well with the published rates derived from experiment. The results for gas properties and surface properties also agree well with the results produced by earlier DSMC models, equilibrium air calculations, and experiment.

Description: Computational fluid dynamics is being used increasingly to predict flows for aerospace propulsion applications, yet there is still a need for an easy to use, computationally inexpensive turbulence model capable of accurately predicting a wide range of turbulent flows. The Baldwin-Lomax model is the most widely used algebraic model, even though it has known difficulties calculating flows with strong adverse pressure gradients and large regions of separation. The modified mixing length model (MML) was developed specifically to handle the separation which occurs on airfoils and has given significantly better results than the Baldwin-Lomax model. The success of these calculations warrants further evaluation and development of MML. The objective of this work was to evaluate the performance of MML for zero and adverse pressure gradient flows, and modify it as needed. The Proteus Navier-Stokes code was used for this study and all results were compared with experimental data and with calculations made using the Baldwin-Lomax algebraic model, which is currently available in Proteus. The MML model was first evaluated for zero pressure gradient flow over a flat plate, then modified to produce the proper boundary layer growth. Additional modifications, based on experimental data for three adverse pressure gradient flows, were also implemented. The adapted model, called MMLPG (modified mixing length model for pressure gradient flows), was then evaluated for a typical propulsion flow problem, flow through a transonic diffuser. Three cases were examined: flow with no shock, a weak shock and a strong shock. The results of these calculations indicate that the objectives of this study have been met. Overall, MMLPG is capable of accurately predicting the adverse pressure gradient flows examined in this study, giving generally better agreement with experimental data than the Baldwin-Lomax model.

Description: Measurements of the performance of hot-film wall-shear sensors were performed to aid development of improved sensors. The effect of film size and substrate properties on the sensor performance was quantified through parametric studies carried out both electronically and in a shock tube. The results show that sensor frequency response increases with decreasing sensor size, while at the same time sensitivity decreases. Substrate effects were also studied, through parametric variation of thermal conductivity and heat capacity. Early studies used complex dual-layer substrates, while later studies were designed for both single-layer and dual-layer substrates. Sensor failures and funding limitations have precluded completion of the substrate thermal-property tests.

Description: Time averaged Stanton number and surface-pressure distributions are reported for the first-stage vane row, the first stage blade row, and the second stage vane row of the Rocketdyne Space Shuttle Main Engine two-stage fuel-side turbine. Unsteady pressure envelope measurements for the first blade are also reported. These measurements were made at 10 percent, 50 percent, and 90 percent span on both the pressure and suction surfaces of the first stage components. Additional Stanton number measurements were made on the first stage blade platform blade tip, and shroud, and at 50 percent span on the second vane. A shock tube was used as a short duration source of heated and pressurized air to which the turbine was subjected. Platinum thin-film heat flux gages were used to obtain the heat flux measurements, while miniature silicon-diaphragm flush-mounted pressure transducers were used to obtain the pressure measurements. The first stage vane Stanton number distributions are compared with predictions obtained using a version of STAN5 and a quasi-3D Navier-Stokes solution. This same quasi-3D N-S code was also used to obtain predictions for the first blade and the second vane.

Description: Hypersonic vehicles are exposed to extreme thermal conditions compared to subsonic aircraft; therefore, some level of thermal management is required to protect the materials used. Normally, hypersonic vehicles experience the highest temperatures in the nozzle throat, and aircraft and propulsion system leading edges. Convective heat transfer augmentation techniques can be used in the thermal management system to increase heat transfer of the cooling channels in those areas. The techniques studied in this report are pin-fin, offset-fin, ribbed and straight roughened channel. A smooth straight channel is used as the baseline for comparing the techniques. SINDA '85, a lumped parameter finite difference thermal analyzer, is used to model the channels. Subroutines are added to model the fluid flow assuming steady one dimensional compressible flow with heat addition and friction. Correlations for convective heat transfer and friction are used in conjunction with the fluid flow analysis mentioned. As expected, the pin-fin arrangement has the highest heat transfer coefficient and the largest pressure drop. All the other devices fall in between the pin-fin and smooth straight channel. The selection of the best heat augmentation method depends on the design requirements. A good approach may be a channel using a combination of the techniques. For instance, several rows of pin-fins may be located at the region of highest heat flux, surrounded by some of the other techniques. Thus, the heat transfer coefficient is maximized at the region of highest heat flux while the pressure drop is not excessive.

Description: We report the results of measurements of the disturbance velocity field generated in a boundary layer by a shallow three-dimensional bump oscillating at a very low frequency on the surface of a flat plate. Profiles of the mean velocity, the disturbance velocity at the fundamental frequency and at the first harmonic are presented. These profiles were measured both upstream and downstream of the oscillating bump. Measurements of the disturbance velocity were also made at various spanwise and downstream locations at a fixed distance from the boundary of one displacement thickness. Finally, the spanwise spectrum of the disturbances at three locations downstream of the bump are presented.

Description: A computer model has been developed which evaluates the performance of a heat exchanger. This model is general enough to be used to evaluate many heat exchanger geometries and a number of different operating conditions. The film approach is used to describe condensation in the presence of noncondensables. The model is also easily expanded to include other effects like fog formation or suction.

Description: Experimental studies have been conducted to examine slot film cooling effectiveness and the interaction between the cooling film and an incident planar shock wave in turbulent hypersonic flow. The experimental studies were conducted in the 48-inch shock tunnel at Calspan at a freestream Mach number of close to 6.4 and at a Reynolds number of 35 x 10(exp 6) based on the length of the model at the injection point. The Mach 2.3 planar wall jet was generated from 40 transverse nozzles (with heights of both 0.080 inch and 0.120 inch), producing a film that extended the full width of the model. The nozzles were operated at pressures and velocities close to matching the freestream, as well as at conditions where the nozzle flows were over- and under-expanded. A two-dimensional shock generator was used to generate oblique shocks that deflected the flow through total turnings of 11, 16, and 21 degrees; the flows impinged downstream of the nozzle exits. Detailed measurements of heat transfer and pressure were made both ahead and downstream of the injection station, with the greatest concentration of measurements in the regions of shock-wave/boundary layer interaction. The major objectives of these experimental studies were to explore the effectiveness of film cooling in the presence of regions of shock-wave/boundary layer interaction and, more specifically, to determine how boundary layer separation and the large recompression heating rates were modified by film cooling. Detailed distributions of heat transfer and pressure were obtained in the incident shock/wall-jet interaction region for a series of shock strengths and impingement positions for each of the two nozzle heights. Measurements were also made to examine the effects of nozzle lip thickness on cooling effectiveness. The major conclusion from these studies was that the effect of the cooling film could be readily dispersed by relatively weak incident shocks, so the peak heating in the recompression region was not significantly reduced by even the largest levels of film cooling. For the case studies in the absence of film cooling, the interaction regions were unseparated. However, adding film cooling resulted in regions of boundary layer separation induced in the film cooling layer -- the size of which regions first increased and then decreased with increased film cooling. Surprisingly, the size of the separated regions and the magnitude of the recompression heating were not strongly influenced by the thickness of the cooling film, nor by the point of shock impingement relative to the exit plane of the nozzles. The lip thickness was found to have little effect on cooling effectiveness. Measurements with and in the absence of shock interaction were compared with the results of earlier experimental studies and correlated in terms of the major parameters controlling these flows.

Description: The calculation of engineering flows undergoing laminar-turbulent transition presents special problems. Mean-flow quantities obey neither the fully laminar nor the fully turbulent correlations. In addition, local maxima in skin friction, wall temperature, and heat transfer often occur near the end of the transition region. Traditionally, modeling this region has been important for the design of turbine blades, where the transition region is long in relation to the chord length of the blade. More recently, the need for better transition-region models has been recognized by designers of hypersonic vehicles where the high Mach number, the low Reynolds number, and the low-disturbance flight environment emphasize the importance of the transition region. Needless to say, a model that might work well for the transitional flows typically found in gas turbines will not necessarily work well for the external surface of a hypersonic vehicle. In Section 2 of this report, some of the important flow features that control the transition region will be discussed. In Section 3, different approaches to the modeling problem will be summarized and cataloged. Fully turbulent flow models will be discussed in detail in Section 4; models specifically designed for transitional flow, in Section 5; and the evaluation of models, in Section 6.

Description: The parabolized stability equations (PSE) are a new approach to analyze the streamwise evolution of single or interacting Fourier modes in weakly nonparallel flows such as boundary layers. The concept rests on the decomposition of every mode into a slowly varying amplitude function and a wave function with slowly varying wave number. The neglect of the small second derivatives of the slowly varying functions with respect to the streamwise variable leads to an initial boundary-value problem that can be solved by numerical marching procedures. The PSE approach is valid in convectively unstable flows. The equations for a single mode are closely related to those of the traditional eigenvalue problems for linear stability analysis. However, the PSE approach does not exploit the homogeneity of the problem and, therefore, can be utilized to analyze forced modes and the nonlinear growth and interaction of an initial disturbance field. In contrast to the traditional patching of local solutions, the PSE provide the spatial evolution of modes with proper account for their history. The PSE approach allows studies of secondary instabilities without the constraints of the Floquet analysis and reproduces the established experimental, theoretical, and computational benchmark results on transition up to the breakdown stage. The method matches or exceeds the demonstrated capabilities of current spatial Navier-Stokes solvers at a small fraction of their computational cost. Recent applications include studies on localized or distributed receptivity and prediction of transition in model environments for realistic engineering problems. This report describes the basis, intricacies, and some applications of the PSE methodology.

Description: This paper is concerned with the effect of a weak spanwise-variable mean-flow distortion on the growth of oblique instability waves in a Blasius boundary layer. The streamwise component of the distortion velocity initially grows linearly with increasing streamwise distance, reaches a maximum, and eventually decays through the action of viscosity. This decay occurs slowly and allows the distortion to destabilize the Blasius flow over a relatively large streamwise region. It is shown that even relatively weak distortions can cause certain oblique Rayleigh instability waves to grow much faster than the usual two-dimensional Tollmien-Schlichting waves that would be the dominant instability modes in the absence of the distortion. The oblique instability waves can then become large enough to interact nonlinearly within a common critical layer. It is shown that the resulting nonlinearity is weak and that the common amplitude of the interacting oblique waves is governed by the amplitude evolution equation derived in Goldstein & Choi (1989). The implications of these results for Klebanoff-type transition are discussed.

Description: A new algorithm for identifying and characterizing vortices in complex flows is presented. The scheme uses both the vorticity and pressure fields. A skeleton line along the center of a vortex is produced by a two-step predictor-corrector scheme. The technique uses the vector field to move in the direction of the skeleton line and the scalar field to correct the location in the plane perpendicular to the skeleton line. A general vortex cross section can be concisely defined with five parameters at each point along the skeleton line. The details of the method and examples of its use are discussed.

Description: An experimental study of the flow over a backward-facing step at low Reynolds number was performed for the purpose of validating a direct numerical simulation (DNS) which was performed by the Stanford/NASA Center for Turbulence Research. Previous experimental data on back step flows were conducted at Reynolds numbers and/or expansion ratios which were significantly different from that of the DNS. The geometry of the experiment and the simulation were duplicated precisely, in an effort to perform a rigorous validation of the DNS. The Reynolds number used in the DNS was Re(sub h)=5100 based on step height, h. This was the maximum possible Reynolds number that could be economically simulated. The boundary layer thickness, d, was approximately 1.0 h in the simulation and the expansion ratio was 1.2. The Reynolds number based on the momentum thickness, Re(sub theta), upstream of the step was 610. All of these parameters were matched experimentally. Experimental results are presented in the form of tables, graphs and a floppy disk (for easy access to the data). An LDV instrument was used to measure mean velocity components and three Reynolds stresses components. In addition, surface pressure and skin friction coefficients were measured. LDV measurements were acquired in a measuring domain which included the recirculating flow region.

Description: Bolgiano scaling in Boussinesq turbulence is analyzed using the Yakhot-Orszag renormalization group. For this purpose, an isotropic model is introduced. Scaling exponents are calculated by forcing the temperature equation so that the temperature variance flux is constant in the inertial range. Universal amplitudes associated with the scaling laws are computed by expanding about a logarithmic theory. Connections between this formalism and the direct interaction approximation are discussed. It is suggested that the Yakhot-Orszag theory yields a lowest order approximate solution of a regularized direct interaction approximation which can be corrected by a simple iterative procedure.

Description: Stationary vortex instabilities with wavelengths significantly larger than the thickness of the underlying three-dimensional boundary layer are studied with asymptotic methods. The long-wavelength Rayleigh modes are locally neutral and are aligned with the direction of the local inviscid streamline. For a spanwise wave number Beta much less than 1, the spatial growth rate of these vortices is O(Beta(exp 3/2)). When Beta becomes O(R(exp -1/7)), the viscous correction associated with a thin sublayer near the surface modifies the inviscid growth rate to the leading order. As Beta is further decreased through this regime, viscous effects assume greater significance and dominate the growth-rate behavior. The spatial growth rate becomes comparable to the real part of the wave number when Beta = O(R(exp -1/4)). At this stage, the disturbance structure becomes fully viscous-inviscid interactive and is described by the triple-deck theory. For even smaller values of Beta, the vortex modes become nearly neutral again and align themselves with the direction of the wall-shear stress. Thus, the study explains the progression of the crossflow-vortex structure from the inflectional upper branch mode to nearly neutral long-wavelength modes that are aligned with the wall-shear direction.

Description: This paper is concerned with the continued development of a new numerical method, the space-time solution element (STS) method, for solving conservation laws. The present work focuses on the two-dimensional, steady, incompressible Navier-Stokes equations. Using first an integral approach, and then a differential approach, the discrete flux conservation equations presented in a recent paper are rederived. Here a simpler method for determining the flux expressions at cell interfaces is given; a systematic and rigorous derivation of the conditions used to simulate the differential form of the governing conservation law(s) is provided; necessary and sufficient conditions for a discrete approximation to satisfy a conservation law in E2 are derived; and an estimate of the local truncation error is given. A specific scheme is then constructed for the solution of the thin airfoil boundary layer problem. Numerical results are presented which demonstrate the ability of the scheme to accurately resolve the developing boundary layer and wake regions using grids which are much coarser than those employed by other numerical methods. It is shown that ten cells in the cross-stream direction are sufficient to accurately resolve the developing airfoil boundary layer.

Description: An advanced engineering model has been developed to aid in the analysis and design of hydrogen/oxygen chemical rocket engines. The complete multispecies, chemically reacting and multidiffusing Navier-Stokes equations are modelled, including the Soret thermal diffusion and the Dufour energy transfer terms. In addition to the spectrum of multispecies aspects developed, the model developed in this study is also conservative in axisymmetric flow for both inviscid and viscous flow environments and the boundary conditions employ a viscous, chemically reacting, reference plane characteristics method. Demonstration cases are presented for a 1030:1 area ratio nozzle, a 25 lbf film cooled nozzle, and a transpiration cooled plug and spool rocket engine. The results indicate that the thrust coefficient predictions of the 1030:1 and the 25 lbf film cooled nozzle are within 0.2 to 0.5 percent, respectively, of experimental measurements when all of the chemical reaction and diffusion terms are considered. Further, the model's predictions agree very well with the heat transfer measurements made in all of the nozzle test cases. The Soret thermal diffusion term is demonstrated to have a significant effect on the predicted mass fraction of hydrogen along the wall of the nozzle in both the laminar flow 1030:1 nozzle and the turbulent flow plug and spool nozzle analysis cases performed. Further, the Soret term was shown to represent an important fraction of the diffusion fluxes occurring in a transpiration cooled rocket engine.

Description: Receptivity to a model convected disturbance in the presence of localized and distributed variations in wall geometry and wall-suction velocity is examined. The model free-stream disturbance corresponds to the time-harmonic wake of a vibrating ribbon that is placed at a suitable distance above the surface of a thin airfoil. The advantages of using this disturbance for experiments on receptivity to convected disturbances are outlined. A brief parametric study is presented for a flat-plate boundary layer. The study quantifies the effect of wake position as well as wake width; in addition, it should be helpful in the choice of an optimal setting for a controlled experiment of the above type, which the above parametric study shows as feasible.

Description: In many cases the structure of a fluid flow is well-characterized by its vortices, especially for the purpose of visualization. In this paper we present a new algorithm for identifying vortices in complex flows. The algorithm produces a skeleton line along the center of a vortex by using a two-step predictor-corrector scheme. The vorticity vector field serves as the predictor and the pressure gradient (in the perpendicular plane) serves as the corrector. We describe an economical description of the vortex tube's cross-section: a 5-term truncated Fourier series is generally sufficient, and it compresses the representation of the flow by a factor of 4000 or more. We reconstruct the vortex tubes as generalized cylinders, providing a polygonal mesh suitable for display on a graphics workstation. We show how the reconstructed geometry of vortex tubes can be enhanced to help visualize helical motion in a static image.

Description: The basic objective of this research was to identify, develop and recommend turbulence models which could be incorporated into CFD codes used in the design of the National AeroSpace Plane vehicles. To accomplish this goal, a combined effort consisting of experimental and theoretical phases was undertaken. The experimental phase consisted of a literature survey to collect and assess a database of well documented experimental flows, with emphasis on high speed or hypersonic flows, which could be used to validate turbulence models. Since it was anticipated that this database would be incomplete and would need supplementing, additional experiments in the NASA Ames 3.5-Foot Hypersonic Wind Tunnel (HWT) were also undertaken. The theoretical phase consisted of identifying promising turbulence models through applications to simple flows, and then investigating more promising models in applications to complex flows. The complex flows were selected from the database developed in the first phase of the study. For these flows it was anticipated that model performance would not be entirely satisfactory, so that model improvements or corrections would be required. The primary goals of the investigation were essentially achieved. A large database of flows was collected and assessed, a number of additional hypersonic experiments were conducted in the Ames HWT, and two turbulence models (kappa-epsilon and kappa-omega models with corrections) were determined which gave superior performance for most of the flows studied and are now recommended for NASP applications.

Description: The Goldstein-Ruban theory has been extended within the framework of Zavol'skii et al. to study the acoustic receptivity of compressible boundary layers. We consider the receptivity produced in a region of localized, small-amplitude variation in the surface temperature and compare it with the receptivity that is induced through a similar mechanism by a variation in the suction velocity at the surface. It is found that the orientation of the acoustic wave can have a significant impact on the receptivity process, with the maximum receptivity at a given sound-pressure level being produced by upstream oriented acoustic waves. At sufficiently low Mach numbers, the variation of receptivity with the acoustic-wave orientation can be predicted analytically and is the same for both surface suction and surface heating. However, as a result of the acoustic refraction across the mean boundary layer, the above dependence can become rather complex and, also, dependent on the type of surface nonuniformity. The results also suggest that the receptivity caused by temperature nonuniformities may turn out to be more significant than that produced by the mean-flow perturbations associated with strip suction.

Description: The problems of understanding the origins of turbulent flow and transition to turbulent flow are the most important unsolved problems of fluid mechanics and aerodynamics. It is well known that the stability, transition, and turbulent characteristics of bounded shear layers are fundamentally different from those of free shear layers. Likewise, the stability, transition, and turbulent characteristics of open systems are fundamentally different from those of closed systems. Because of the influence of indigenous disturbances, surface geometry and roughness, sound, heat transfer, and ablation, it is not possible to develop general prediction schemes for transition location and the nature of turbulent structures in boundary-layer flows. At the present time no mathematical model exists that can predict the transition Reynolds number on a flat plate. The recent progress in this area is encouraging, in that a number of distinct transition mechanisms have been found experimentally. The theoretical work finds them to be amplitude and Reynolds-number dependent. The theory remains rather incomplete with regard to predicting transition. Amplitude and spectral characteristics of the disturbances inside the laminar viscous layer strongly influence which type of transition occurs. The major need in this area is to understand how freestream disturbances are entrained into the boundary layer, i.e., to answer the question of receptivity. We refer receptivity to the mechanism(s) that cause freestream disturbances to enter the boundary layer and create the initial amplitudes for unstable waves.

Description: This document is intended to serve as an installation, usage, and basic theory guide for the two dimensional triangulation software 'HARLEY' written for the Silicon Graphics IRIS workstation. This code consists of an incremental triangulation algorithm based on point insertion and local edge swapping. Using this basic strategy, several types of triangulations can be produced depending on user selected options. For example, local edge swapping criteria can be chosen which minimizes the maximum interior angle (a MinMax triangulation) or which maximizes the minimum interior angle (a MaxMin or Delaunay triangulation). It should be noted that the MinMax triangulation is generally only locally optical (not globally optimal) in this measure. The MaxMin triangulation, however, is both locally and globally optical. In addition, Steiner triangulations can be constructed by inserting new sites at triangle circumcenters followed by edge swapping based on the MaxMin criteria. Incremental insertion of sites also provides flexibility in choosing cell refinement criteria. A dynamic heap structure has been implemented in the code so that once a refinement measure is specified (i.e., maximum aspect ratio or some measure of a solution gradient for the solution adaptive grid generation) the cell with the largest value of this measure is continually removed from the top of the heap and refined. The heap refinement strategy allows the user to specify either the number of cells desired or refine the mesh until all cell refinement measures satisfy a user specified tolerance level. Since the dynamic heap structure is constantly updated, the algorithm always refines the particular cell in the mesh with the largest refinement criteria value. The code allows the user to: triangulate a cloud of prespecified points (sites), triangulate a set of prespecified interior points constrained by prespecified boundary curve(s), Steiner triangulate the interior/exterior of prespecified boundary curve(s), refine existing triangulations based on solution error measures, and partition meshes based on the Cuthill-McKee, spectral, and coordinate bisection strategies.

Description: The present invention is directed to a method and apparatus for improving the flow range in centrifugal pumps and compressors. Bleed holes are introduced into a volute tongue of a centrifugal pump or compressor thereby providing a double acting means of boundary layer control at the volute tongue.

Description: An increasing amount of research activity in computational fluid dynamics has been devoted to the development of efficient algorithms for parallel computing systems. The increasing performance to price ratio of engineering workstations has led to research to development procedures for implementing a parallel computing system composed of distributed workstations. This thesis proposal outlines an ongoing research program to develop efficient strategies for performing three-dimensional flow analysis on distributed computing systems. The PVM parallel programming interface was used to modify an existing three-dimensional flow solver, the TEAM code developed by Lockheed for the Air Force, to function as a parallel flow solver on clusters of workstations. Steady flow solutions were generated for three different wing and body geometries to validate the code and evaluate code performance. The proposed research will extend the parallel code development to determine the most efficient strategies for unsteady flow simulations.

Description: Direct numerical simulations of compressible, homogeneous, turbulent shear flows are used to evaluate Reynolds stress models. Three pressure-strain models, which are either linear, quadratic, or cubic in the anisotropy tensor are considered. Dilatational dissipation and pressure-strain correlation models do not correctly capture the compressibility effects seen in the direct simulations. In particular, the increase in the anisotropy of normal stresses and the reduction in the shear stress are not reproduced by any of the models. Also, the use of the incompressible form of the dissipation-rate equation to determine the solenoidal part of the dissipation is found to be questionable.

Description: We present a new approach towards the construction of a genuinely multidimensional high-resolution scheme for computing steady-state solutions of the Euler equations of gas dynamics. The unique advantage of this approach is that the Gauss-Seidel relaxation is stable when applied directly to the high-resolution discrete equations, thus allowing us to construct a very efficient and simple multigrid steady-state solver. This is the only high-resolution scheme known to us that has this property. The two-dimensional scheme is presented in detail. It is formulated on triangular (structured and unstructured) meshes and can be interpreted as a genuinely two-dimensional extension of the Roe scheme. The quality of the solutions obtained using this scheme and the performance of the multigrid algorithm are illustrated by the numerical experiments. Construction of the three dimensional scheme is outlined briefly as well.

Description: The effect of basic flow density gradients on the supersonic wall modes were investigated in Part 1 of this analysis. In that investigation only the 2-D modes were studied. Tam and Hu investigated the 3-D modes in a confined vortex sheet and reported that the first 2-D Class A mode (A01) had the highest growth rate compared to all other 2-D and 3-D modes present in the vortex sheet for that particular set of flow patterns. They also showed that this result also held true for finite thickness shear layers with delta(sub w) less than 0.125. For free shear layers, Sandham and Reynolds showed that the 3-D K-H mode became the dominant mode for M(sub c) greater than 0.6. Jackson and Grosch investigated the effect of crossflow and obliqueness on the slow and fast odes present in a M(sub c) greater than 1 environment and showed that for certain combination of crossflow and wave angles the growth rates could be increased by up to a factor of 2 with respect to the 2-D case. The case studied here is a confined shear layer shown in Part 1. All solution procedures and basic low profiles are the same as in Part 1. The effect of density gradients on the 3-D modes present in the density ratios considered in Part 1 are investigated.

Description: The effect of density gradients on the supersonic wall modes (acoustic modes) of a 2-D confined compressible shear layer were investigated using linear analysis. Due to the inadequacies of the hyperbolic tangent profile, the boundary layer basic profiles were used. First a test case was taken with the same parameters as in Tam and Hu's analysis with convective Mach number M(sub c) = 1.836 and density ratio of 1.398. Three generalized inflection points were found giving rise to three modes. The first two show similar properties to the Class A and B modes, and the third is an 'inner mode' which will be called a Class C mode. As the density ratio is increased, the smallest of the three neutral phase speeds tends towards the speed of the lower velocity stream, and the other two eventually coalesce and then disappear. These two effects lead to a linear resonance between the Class B modes which increases the cutoff frequency and growth rate of the lowest mode. In fact, growth rates of 2-4 times the test case were found as the density ratio was increased to 7. A similar trend is observed for the Class A modes when the density ratio is decreased from the test case, but the growth rate is not changed by much from the test case.

Description: The purpose of the present work was threefold: (1) to determine if a free-stream turbulence length scale existed that would cause the greatest augmentation in stagnation-region heat transfer over laminar levels; (2) to investigate the effect of velocity gradient on stagnation-region heat transfer augmentation by free-stream turbulence; and (3) to develop a prediction tool for stagnation heat transfer in the presence of free-stream turbulence. Heat transfer was measured in the stagnation region of four models with elliptical leading edges that had ratios of major to minor axes of 1:1, 1.5:1, 2.25:1, and 3:1. Five turbulence-generating grids were fabricated; four were square mesh, biplane grids made from square bars. The fifth grid was an array of fine parallel wires that were perpendicular to the model spanwise direction. Heat transfer data were taken at Reynolds numbers ranging from 37 000 to 228 000. Turbulence intensities were in the range of 1.1 to 15.9% while the ratio of integral length scale to leading-edge diameter ranged from 0.05 to 0.30. Stagnation-point velocity gradient was varied by nearly 50%. Stagnation-region heat transfer augmentation was found to increase with decreasing length scale but no optimum length scale was found. Heat transfer augmentation due to turbulence was found to be unaffected by the velocity gradient near the leading edge. A correlation was developed that fit heat transfer data for the square-bar grids to within +/- 4%.

Description: Direct numerical simulations (DNS) of the Navier-Stokes equations are employed to explore the effects of suction holes on transition in a laminar flat plate boundary layer. The Navier-Stokes equations are cast in vorticity-velocity form. Periodicity is imposed in spanwise direction; all other spatial derivatives are discretized with fourth order compact differences. An explicit fourth order Runge-Kutta scheme is employed for the time-integration of the vorticity transport equations. Suction is applied through a row of holes aligned in spanwise direction. For low suction strengths, each hole generates a pair of stable streamwise vortices. When the suction strength exceeds a critical value, the vortices become unstable. For high suction strengths, vortex shedding occurs right at the suction holes. Our numerical findings agree well with experimental observations.

Description: A set of direct numerical simulations of temporally evolving incompressible plane wakes started from a variety of initial conditions have been performed. The intent is to explore the effect of initial conditions on the development of three dimensionality. It was found that the existence of an oblique disturbance at a streamwise wavelength equal to the two-dimensional subharmonic wavelength plays a very important role in the development of the wake. The Reynolds number and the presence or absence of longer wavelength disturbances were also found to be important at late times in the development of the flow. Disturbance phase was found to have an effect on the details of the structure of the wake, but has relatively little effect on the wake growth.

Description: Direct and large eddy simulations of forced and decaying isotropic turbulence have been performed using a pseudospectral and a finite-difference code. Subgrid models that include a one-equation subgrid kinetic energy model with and without a stochastic backscatter forcing term and a new scale similarity model have been analyzed in both Fourier space and physical space. The Fourier space analysis showed that the energy transfer across the cutoff wavenumber k(sub c) is dominated by local interaction. The correlation between the exact and the modeled (by a spectral eddy viscosity) nonlinear terms and the subgrid energy transfer in physical space was found to be quite low. In physical space, a similar correlation analysis was carried out using top hat filtering. Results show that the subgrid stress and the energy flux predicted by the subgrid models correlates very well with the exact data. The scale similarity model showed very high correlation for reasonable grid resolution. However, with decrease in grid resolution, the scale similarity model became more uncorrelated, when compared to the kinetic energy subgrid model. The subgrid models were then used for large-eddy simulations for a range of Reynolds number. It was determined that the dissipation was modeled poorly and that the correlation with the exact results was quite low for all the models. In general, for coarse grid resolution, the scale similarity model consistently showed very low correlation while the kinetic energy model showed a relatively higher correlation. These results suggest that to use the scale similarity model relatively fine grid resolution may be required, whereas, the kinetic energy model could be used even in coarse grid.

Description: Briefs of the 1994 annual progress reports of the Research Fellows and students of the Center for Turbulence Research are presented. Subjects covered include turbulence combustion, large eddy simulation, Reynolds-averaged turbulence modeling, turbulence control, postprocessing, sound generation, and turbulence physics.

Description: The restricted Euler equation captures many important features of the behavior of the velocity gradient tensor observed in direct numerical simulations (DNS) of isotropic turbulence. However, in slightly more complex flows the agreement is not good, especially in regions of low dissipation. In this paper, it is demonstrated that the Reynolds-averaged restricted Euler equation violates the balance of mean momentum for virtually all homogeneous turbulent flows with only two major exceptions: isotropic and homogeneously-sheared turbulence. A new model equation which overcomes this shortcoming and is more widely applicable is suggested. This modele is derived from the Navier-Stokes equation with a restricted Euler type approximation made on the fluctuating velocity gradient field. Analytical solutions of the proposed modified restricted Euler equation appear to be difficult to obtain. Hence, a strategy for numerically calculating the velocity gradient tensor is developed. Preliminary calculations tend to indicate that the modified restricted Euler equation captures many important aspects of the behavior of the fluctuating velocity gradients in anisotropic homogeneous turbulence.

Description: It is known that certain configurations which possess curvature are prone to a class of instabilities which their 'flat' counterparts will not support. The main thrust of the study of these centrifugal instabilities has concentrated on curved solid boundaries and their effect on the fluid motion. In this article attention is shifted towards a fluid-fluid interface observed within a curved incompressible mixing layer. Experimental evidence is available to support the conjecture that this situation may be subjected to centrifugal instabilities. The evolution of modes with wavelengths comparable with the layer's thickness is considered and the high Taylor/Gortler number regime is also discussed which characterizes the ultimate fate of the modes.

Description: The present study examines the wind-tunnel blockage and actuation systems effectiveness in starting and forcibly unstarting a two-dimensional scramjet inlet in the NASA Langley 20-Inch Mach 6 Tunnel. The intent of the overall test program is to study (both experimentally and computationally) the dynamics of the inlet unstart; however, prior to the design and fabrication of an expensive, instrumented wind-tunnel model, it was deemed necessary first to examine potential wind-tunnel blockage issues related to model sizing and to examine the adequacy of the actuation systems in accomplishing the start and unstart. The model is equipped with both a moveable cowl and aft plug. Windows in the inlet sidewalls allow limited optical access to the internal shock structure; schlieren video was used to identify inlet start and unstart. A chronology of each actuation sequence is provided in tabular form along with still frames from the schlieren video. A pitot probe monitored the freestream conditions throughout the start/unstart process to determine if there was a blockage effect due to the model start or unstart. Because the purpose of this report is to make the phase I (blockage and actuation systems) data rapidly available to the community, the data is presented largely without analysis of the internal shock interactions or the unstart process. This series of tests indicated that the model was appropriately sized for this facility and identified operability limits required first to allow the inlet to start and second to force the unstart.

Description: Linear stability theory is used to study the effect of crossflow on Goertler instability in incompressible boundary layers. The results cover a wide range of sweep angle, pressure gradient, and wall curvature parameters. It is shown that the crossflow stabilizes Goertler disturbances by reducing the maximum growth rate and shrinking the unstable band of spanwise wave numbers. On the other hand, the effect of concave wall curvature on crossflow instability is destabilizing. Calculations show that the changeover from Goertler to crossflow instabilities is a function of Goertler number, pressure gradient, and sweep angle. The results demonstrate that Goertler instability may still be relevant in the transition process on swept wings even at large angles of sweep if the pressure gradient is sufficiently small. The influence of pressure gradient and sweep can be combined by defining a crossflow Reynolds number. Thus, the changeover from Goertler to crossflow instability takes place at some critical crossflow Reynolds number whose value increases with Goertler number.

Description: We use numerical and asymptotic techniques to study the stability of a two-phase air/water flow above a flat porous plate. This flow is a model of the boundary layer which forms on a yawed cylinder and can be used as a useful approximation to the air flow over swept wings during heavy rainfall. We show that the interface between the water and air layers can significantly destabilize the flow, leading to traveling wave disturbances which move along the attachment line. This instability occurs for lower Reynolds numbers than in the case of the absence of a water layer. We also investigate the instability of inviscid stationary modes. We calculate the effective wavenumber and orientation of the stationary disturbance when the fluids have identical physical properties. Using perturbation methods we obtain corrections due to a small stratification in viscosity, thus quantifying the interfacial effects. Our analytical results are in agreement with the numerical solution which we obtain for arbitrary fluid properties.

Description: A general class of exact solutions is presented for a time evolving bubble in a two-dimensional slow viscous flow in the presence of surface tension. These solutions can describe a bubble in a linear shear flow as well as an expanding or contracting bubble in an otherwise quiescent flow. In the case of expanding bubbles, the solutions have a simple behavior in the sense that for essentially arbitrary initial shapes the bubble will asymptote an expanding circle. Contracting bubbles, on the other hand, can develop narrow structures ('near-cusps') on the interface and may undergo 'break up' before all the bubble-fluid is completely removed. The mathematical structure underlying the existence of these exact solutions is also investigated.

Description: In 1952, Lighthill developed a theory for determining the sound generated by a turbulent motion of a fluid. With some statistical assumptions, Proudman applied this theory to estimate the acoustic power of isotropic turbulence. Recently, Lighthill established a simple relationship that relates the fourth-order retarded time and space covariance of his stress tensor to the corresponding second-order covariance and the turbulent flatness factor, without making statistical assumptions for a homogeneous turbulence. Lilley revisited Proudman's work and applied the Lighthill relationship to evaluate directly the radiated acoustic power from isotropic turbulence. After choosing the time separation dependence in the two-point velocity time and space covariance based on the insights gained from direct numerical simulations, Lilley concluded that the Proudman constant is determined by the turbulent flatness factor and the second-order spatial velocity covariance. In order to estimate the Proudman constant at high Reynolds numbers, we analyzed a unique data set of measurements in a large wind tunnel and atmospheric surface layer that covers a range of the Taylor microscale based on Reynolds numbers 2.0 x 10(exp 3) less than or equal to R(sub lambda) less than or equal to 12.7 x 10(exp 3). Our measurements demonstrate that the Lighthill relationship is a good approximation, providing additional support to Lilley's approach. The flatness factor is found between 2.7 - 3.3 and the second order spatial velocity covariance is obtained. Based on these experimental data, the Proudman constant is estimated to be 0.68 - 3.68.

Description: A novel experiment to study the interfacial shear instability between two liquids is described. Density-matched immiscible liquids are confined between concentric cylinders such that the interface is parallel to the cylinder walls. Interfacial waves that develop because of viscosity differences between the shearing fluids are studied as a function of rotation rate and depth ratio using optical techniques. Conditions neutral stability and the most unstable wavenumber agree reasonably well with predictions from linear stability analysis of the Navier-Stokes equations. Illumination using laser sheets allows precise measurement of the interface shape. Future experiments will verify the correctness of weakly nonlinear theories that describe energy transfer and saturation of wave growth by nonlinear effects. Measurements of solitary wave shapes, that occur far above neutral stability, will be compared to similar measurements for systems that have gravity as an important force to determine how gravity effects large disturbances. These results will be used to interpret slug and annular flow data that have been obtained in other mu g studies.

Description: Theoretical experimental investigations of the rewetting characteristics of thin liquid films over unheated and heated grooved plates were performed. Studied factors which affected the rewetting characteristics of the plate were mass flow rate of the rewetting liquid on the plate, the orientation of the plate (face up, face down or inclined) and smooth/groove surface conditions. The initial plate temperature was also varied, with experiments being performed between room temperature and 150 C. It was found that the rewetting velocity increased with the initial plate temperature. But when the temperature was increased further above the Leidenfrost temperature (liquid front temperature), the rewetting velocity decreases with the initial plate temperature. Hydrodynamically controlled and conductively controlled rewetting models were presented to explain and to predict the rewetting characteristics in these two distinct regions. Also found was a higher rewetting velocity when the grooved plate was placed in the face up orientation than in the face down orientation.

Description: A series of two-phase gas-liquid flow experiments were developed to study annular flows in microgravity using the NASA Lewis Learjet. A test section was built to measure the liquid film thickness around the perimeter of the tube permitting the three dimensional nature of the gas-liquid interface to be observed. A second test section was used to measure the film thickness, pressure drop and wall shear stress in annular microgravity two-phase flows. Three liquids were studied to determine the effects of liquid viscosity and surface tension. The result of this study provide insight into the wave characteristics, pressure drop and droplet entrainment in microgravity annular flows.

Description: This paper describes the status and plans for the work being performed under NASA NRA contract NASW-4803 so that members of the Microgravity Fluid Dynamics Discipline Working Group are aware of this program. The contract is a cross-disciplinary research program and is administered under the Low Temperature Microgravity Research Program at the Jet Propulsion Laboratory. The purpose of the project is to perform low-gravity verification experiments on the slosh behavior of He II to use in the development of a CFD model that incorporates the two-fluid physics of He II. The two-fluid code predicts a different fluid motion response in low-gravity environment from that predicted by a single-fluid model, while the 1g response is identical for the both types of model.

Description: Ground-based modeling and experiments have been performed on the interaction and coalescence of drops leading to macroscopic phase separation. The focus has been on gravity-induced motion, with research also initiated on thermocapillary motion of drops. The drop size distribution initially shifts toward larger drops with time due to coalescence, and then a back towards smaller drops due to the larger preferentially settling out. As a consequence, the phase separation rate initially increases with time and then decreases.

Description: Wall functions, as used in the typical high Reynolds number k-epsilon turbulence model, can be implemented in various ways. A least disruptive method (to the flow solver) is to directly solve for the flow variables at the grid point next to the wall while prescribing the values of k and epsilon. For the centrally-differenced finite-difference scheme employing artificial viscocity (AV) as a stabilizing mechanism, this methodology proved to be totally useless. This is because the AV gives rise to a large error at the wall due to too steep a velocity gradient resulting from the use of a coarse grid as required by the wall function methodology. This error can be eliminated simply by extrapolating velocities at the wall, instead of using the physical values of the no-slip velocities (i.e. the zero value). The applicability of the technique used in this paper is demonstrated by solving a flow over a flat plate and comparing the results with those of experiments. It was also observed that AV gives rise to a velocity overshoot (about 1 percent) near the edge of the boundary layer. This small velocity error, however, can yield as much as 10 percent error in the momentum thickness. A method which integrates the boundary layer up to only the edge of the boundary (instead of infinity) was proposed and demonstrated to give better results than the standard method.

Description: A method is proposed for determining the response of the steady swirling potential flow past a blade row in an axial turbomachine to the most general type of 'nonacoustic' incident disturbance. The method is based on Goldstein's decomposition of the disturbance velocity and only requires solving a linear inhomogeneous wave equation. It is believed that numerical solutions to this wave equation can be obtained more efficiently and will be more accurate than corresponding solutions to the linearized compressible Euler equations.

Description: A new k-epsilon eddy viscosity model, which consists of a new model dissipation rate equation and a new realizable eddy viscosity formulation, is proposed. The new model dissipation rate equation is based on the dynamic equation of the mean-square vorticity fluctuation at large turbulent Reynolds number. The new eddy viscosity formulation is based on the realizability constraints: the positivity of normal Reynolds stresses and Schwarz' inequality for turbulent shear stresses. We find that the present model with a set of unified model coefficients can perform well for a variety of flows. The flows that are examined include: (1) rotating homogeneous shear flows; (2) boundary-free shear flows including a mixing layer, planar and round jets; (3) a channel flow, and flat plate boundary layers with and without a pressure gradient; and (4) backward facing step separated flows. The model predictions are compared with available experimental data. The results from the standard k-epsilon eddy viscosity model are also included for comparison. It is shown that the present model is a significant improvement over the standard k-epsilon eddy viscosity model.