ALBERT

Description: In the design of a combustor, information is necessary for the mixing of the fuel and air in order to determine the optimum combustor length. In scramjet combustors the mixing often takes place in a shear layer that is formed between the fuel and air. This research was an experimental study of shear layers in supersonic flows aimed at determining what mechanisms affect the shear layer so that the mixing could be better predicted. A second goal was to provide sufficient instream information for use in checking existing Computational Fluid Dynamic (CFD) codes. The shear layer between a supersonic two-dimensional air stream (M = 2 or M = 3) was mixed with a near sonic two-dimensional air stream (M = 1.2). Instream measurements of pitot pressure and cone static pressure were used to determine mean velocity profiles at various axial locations. These velocity profiles were used to determine the shear layer spreading rate and are compared with various predictions. Wall measurements of static pressure, temperature and skin friction were also taken and are presented. The instream measurements were also used for comparison with an existing CFD code. The upstream velocity, pressure and temperature profiles were used as a starting profile and the code was used to calculate downstream profiles for comparison with the experimental results. Reasonable agreement between the measured and calculated results was obtained.

Description: This paper presents a set of second-order closure models for low-Reynolds-number turbulence near the wall. Existing closure models for the Reynolds-stress equations were modified to show proper near-wall behavior. A dissipation-rate equation for the turbulent kinetic energy is also reformulated. The proposed models satisfy realizability and will not produce unphysical behavior. Fully developed channel flows are used for model testing. The equations are solved for the mean velocity, the Reynolds stresses, and the dissipation rate of the turbulent kinetic energy. The calculations are compared with both direct numerical simulations and with measurements. It is shown that the present models perform well in predicting the behavior of the turbulence near a wall. Significant improvements over previous models in predicting the components of the Reynolds stress tensor are obtained in the present models.

Description: The low-g fluids management group with the Center for Space Construction is engaged in active research on the following topics: gauging; venting; controlling contamination; sloshing; transfer; acquisition; and two-phase flow. Our basic understanding of each of these topics at present is inadequate to design space structures optimally. A brief report is presented on each topic showing the present status, recent accomplishings by our group and our plans for future research. Reports are presented in graphic and outline form.

Description: Experimental and theoretical studies have been conducted to determine critical parameters at the onset of nonlinear counterflow in He II below the lambda point of He-4. Critical temperature differences have been measured in porous media for zero net mass flow and for Darcy permeabilities in the order of magnitude range from 10 to the -10th to 10 to the -8th sq cm. The normalized critical temperature gradients, which covered the liquid temperature range of 1.5 K to the lambda temperature, are found to vary with T proportional to the ratio of the superfluid density to the normal fluid density. This liquid temperature dependence appears to be consistent with duct data which are limited at low temperature by a Reynolds number criterion.

Description: An analysis is made based upon the concept that the velocity fluctuations, and therefore, the Reynolds stresses, driven by the instability of the original flow grow until a new stable state is approached. The Reynolds stresses incorporated into the Orr-Sommerfeld equation are coupled with the main flow such that all the imaginary parts of the complex eigenvalues vanish, i.e., the original instability is eliminated. Using this stabilization principle, it is possible to find the Reynolds stresses as well as the mean velocity for plane Poiseuille flow with the Reynolds number slightly higher than the critical.

Description: The behavior of the reverse flow ceiling jet against the ventilation flow from 0.58 to 0.87 m/s was investigated in a 1/3 scale model of a wide body aircraft interior. For all tests, strong reverse-flow ceiling jets of hot gases were detected well upstream of the fire. Both thicknesses of the reverse-flow ceiling jet and the smoke layer increased with the fire-crossflow parameter. The thickness of the smoke layer where the smoke flows along the main flow below the reverse-flow ceiling jet was almost twice that of the reverse-flow ceiling jet. Detailed spatial and time-varying temperatures of the gas in the test section were measured, and velocity profiles were also measured using a temperature compensated hot film.

Description: Certain theoretical studies of boundary-layer transition are described, based on high Reynolds numbers and with attention drawn to the various nonlinear interactions and scales present. The article concentrates in particular on theories for which the mean-flow profile is completely altered from its original state. Two- and three-dimensional flow theory and conjectures on turbulent-boundary-layer structures are included. Specific recent findings noted, and in qualitative agreement with experiments, are: nonlinear finite-time break-ups in unsteady interactive boundary layers; strong vortex/wave interactions; and prediction of turbulent boundary-layer displacement- and stress sublayer-thicknesses.

Description: The direct interaction approximation is used to treat turbulence in a compressible medium with constant mean gradients. The set of coupled nonlinear integrodifferential equations is derived that is satisfied by the transverse and longitudinal energy spectral functions, Q(T) and Q(L) and by the transverse and longitudinal response functions, G(T) and G(L). Finally, expressions for the average of the product of pairs of physically relevant fluctuating quantities (velocity, temperature, density) are derived in terms of Q(T) and Q(L).

Description: The conditions under which finite difference schemes for the shallow water equations can conserve both total energy and potential enstrophy are considered. A method of deriving such schemes using operator formalism is developed. Several such schemes are derived for the A-, B- and C-grids. The derived schemes include second-order schemes and pseudo-fourth-order schemes. The simplest B-grid pseudo-fourth-order schemes are presented.

Description: Renormalization-group theory is applied to incompressible three-dimensional Navier-Stokes turbulence so as to eliminate unresolvable small scales. The renormalized Navier-Stokes equation now includes a triple nonlinearity with the eddy viscosity exhibiting a mild cusp behavior, in qualitative agreement with the test-field model results of Kraichnan. For the cusp behavior to arise, not only is the triple nonlinearity necessary but the effects of pressure must be incorporated in the triple term. The renormalized eddy viscosity will not exhibit a cusp behavior if it is assumed that a spectral gap exists between the large and small scales.

Description: The physical mechanism governing the centring of a hollow liquid shell in capillary oscillations, which has been observed in experiments, is investigated theoretically. First, the shell is assumed to be inviscid and to have a thickness that is much less than its spherical radius. A system of one-dimensional nonlinear equations of motion is derived using a thin-sheet model. From a numerical study the nonlinear effects of the wave are found to cause the core to oscillate slowly relative to the shell while the centre of mass of the whole system remains stationary. The effects of small viscosity are then considered in an approximation. Finally the strength of the centring mechanism is compared with that of the decentring effect due to buoyancy. The findings are consistent with the limited experimental information available.

Description: Developing fluid flow in a curved duct of square cross-section is studied numerically by a factored ADI finite-difference method on a staggered grid. A central-difference scheme with primitive variables is used inside the computational domain to reduce numerical diffusion. Two Reynolds numbers, 574 and 790, based upon a bulk velocity and hydraulic diameter are chosen for curvature ratios of 1/6.45 and 1/2.3, respectively. It is found that the secondary flow is far more complicated than expected, with the appearance of at least two pairs of vortices. Main-flow separation is also observed for the higher curvature ratio. Furthermore, it is observed that the flow develops into two quite different states downstream, depending upon the inlet conditions. Solutions of the fully developed Navier-Stokes equations is shown to be not unique beyond a certain critical Reynolds number. Developing flow seems to evolve into the fully developed state along a particular branch into which the fully developed solution bifurcates.

Description: A development history and current status evaluation are presented for the theory of permeability and percolation. The microscale phenomena treated in this field have proven difficult to analyze due both to their tortuous geometry and the influence of capilarity. Capilary effects may be not only important but predominant, and are differentiated into those at the fluid-fluid interface, and those involving the existence of a contact line between the solid substrate and this interface. Percolation theory has been borrowed from physics and adapted to the two-phase engineering context.

Description: The morphological stability of a rotating and solidifying disk is investigated under the assumption that delta, the thickness of the viscous boundary layer, is much larger than delta(c), the thickness of the solute boundary layer. It is found that axisymmetric disturbances with wavelengths comparable to delta respond to nonparallel flow effects and have stability characteristics quite different from disturbances in a parallel flow. These long waves are unstable because of the nonparallel flow and would decay without it. This analysis thus identifies a new mechanism of morphological change induced by flow.

Description: In this paper, the receptivity of a typical free shear layer to pulse-type and periodic excitation is studied. This is done by solving the initial-value problem completely and studing its long-time behavior. This leads to a wave packet for the pulse. By the superposition of many wave packets, a spatial instability mode when the flow is convectively unstable is generated. This establishes a general and simple relationship between the receptivities for pulse-type and sinusoidal excitations. It is found that a shear layer is very receptive to high-frequency disturbances that are generated near the centerline of the layer.

Description: The turbulent boundary layer on a flat plate, with zero pressure gradient, is simulated numerically at four stations between R sub theta = 225 and R sub theta = 1410. The three-dimensional time-dependent Navier-Stokes equations are solved using a spectral method with up to about 10 to the 7th grid points. Periodic spanwise and streamwise conditions are applied, and a multiple-scale procedure is applied to approximate the slow streamwise growth of the boundary layer. The flow is studied, primarily, from a statistical point of view. The solutions are compared with experimental results. The scaling of the mean and turbulent quantities with Reynolds number is compared with accepted laws, and the significant deviations are documented. The turbulence at the highest Reynolds number is studied in detail. The spectra are compared with various theoretical models. Reynolds-stress budget data are provided for turbulence-model testing.

Description: A method is proposed for treating steady-state, thermally driven convection using the full direct interaction approximation to treat the nonlinear energy transfer and a prescribed instability function to model the eenrgy input. The instability function used is the growth rate of the mechanism generating the turbulence. This methodology can be easily generalized and applied to other driving mechanisms. The 1/3 power law form of the N vs. R relation for water is duplicated here and the coefficient is computed, using a two-point closure, to be less than about 0.08.

Description: The relationship between the spectra of the three normal stresses near a turbulent shear layer is derived from the constraint of irrotationality. The relationship implies that the common practice of determining time scales of the turbulence from spectra in the near field is not valid. This implication is reinforced by a discussion of the scaling of irrotational spectra that follows from the formulation of Phillips (1975).

Description: The structure of the vorticity field in the viscous wall layer of a turbulent channel is studied by examining the results of a fully resolved direct numerical simulation. It is shown that this region is dominated by intense three-dimensional shear layers in which the dominant vorticity component is spanwise. The advection and reproduction processes of these structures are examined and shown to be consistent with the classical generation mechanism for two-dimensional Tollmien-Schlichting waves. This process is fundamentally different from the usually accepted mechanism involving hairpin vortices.

Description: An algebraic stress turbulence model and a computational procedure for turbulent boundary layer flows which is based on the semidiscrete Galerkin FEM are discussed. In the algebraic stress turbulence model, the eddy viscosity expression is obtained from the Reynolds stress turbulence model, and the turbulent kinetic energy dissipation rate equation is improved by including a production range time scale. Good agreement with experimental data is found for the examples of a fully developed channel flow, a fully developed pipe flow, a flat plate boundary layer flow, a plane jet exhausting into a moving stream, a circular jet exhausting into a moving stream, and a wall jet flow.

Description: Free-streamline theory is used to analyze the deformation and dislodging by wind pressure of drops of liquid adhered by surface tension to a solid surface. The critical Weber number for droplets to be dislodged is determined as a function of advancing and receding contact angle. Graphical results for drop shape are in good agreement with observation.

Description: This paper presents the application of a class of multi-grid methods to the solution of the Navier-Stokes equations for two-dimensional laminar flow problems. The methods consists of combining the full approximation scheme-full multi-grid technique (FAS-FMG) with point-, line- or plane-relaxation routines for solving the Navier-Stokes equations in primitive variables. The performance of the multi-grid methods is compared to those of several single-grid methods. The results show that much faster convergence can be procured through the use of the multi-grid approach than through the various suggestions for improving single-grid methods. The importance of the choice of relaxation scheme for the multi-grid method is illustrated.

Description: Rapid distortion theory is applied to study distortion of homogeneous turbulence subject to two different axisymmetric strain modes: the axisymmetric contraction (AC, nozzle-type flow), and the axisymmetric expansion (AE, diffuser-type flow). The paper explores the differences in effects of the two axisymmetric strain modes on the anisotropy of correlations and structures of turbulence; examines the effect of dilatation on the distortion of turbulence; and provides a theoretical background for turbulence model development. It is found that velocity and vorticity fluctuations are enhanced more efficiently by contraction than by expansion; contraction produces much higher anisotropy in velocity and vorticity than expansion; root-mean-square pressure is slightly reduced during contraction, whereas it increases rapidly during expansion; and vortical structures of rodlike shape develop in a contraction flow, while disklike structures develop in an expansion flow. A simple model that reflects the dependence of turbulence evolution on structural parameters such as the Reynolds-stress anisotropy and total strain is proposed, and is shown to outperform all other models for all cases examined, regardless of the mean strain rate.

Description: Numerical studies of turbulent flow in an axisymmetric 45-deg-expansion combustor and bifurcated diffuser are presented. The Navier-Stokes equations incorporating a k-epsilon model were solved in a nonorthogonal curvilinear coordinate system. A zonal-grid method, where the flow field was divided into several subsections, was developed. This approach permitted different computational schemes to be used in the various zones. In addition, grid generation was made a more simple task. Boundary overlap and interpolating techniques were used, and an adjustment of the flow variables was required to assure conservation of mass flux. Three finite-differencing methods (hybrid, quadratic upwind, and skew upwind) were used to represent the convection terms. Results were compared with existing experimental data. In general, good agreement between predicted and measured values was obtained.

Description: It has previously been shown that the no-slip boundary conditions leads to a singularity at a moving contact line and that this presumes some form of slip. Present considerations on the energetics of slip due to shear stress lead to a yield stress boundary condition. A model for the distortion of the liquid state near solid boundaries gives a physical basis for this boundary condition. The yield stress condition is illustrated by an analysis of a slender drop rolling down an incline. That analysis provides a formula for the frictional drag resisting the drop movement. With the present boundary condition, the length of the slip region becomes a property of the fluid flow.

Description: The numerical scheme used by the present time-accurate FEM numerical method for incompressible Navier-Stokes equations, using primitive variables as the unknowns, is a Crank-Nicholson implicit treatment of all equation terms with central differencing for space derivatives. The introduction of a continuous auxilliary system in pseudo-time, with artificial compressibility, yields the incompressible solution at the advanced time level; time-accurate solutions are thereby obtained for two-dimensional fluid flows in a square cavity, in the cases of both an impulsively starting lid and an oscillating lid.

Description: The paper describes an investigation of the velocity fluctuations in the free stream above an incompressible turbulent boundary layer developing at constant pressure. Measurements were performed in a wind tunnel with an rms level of axial velocity fluctuations of about 0.2 percent. The possibility of assessing the contributions to the rms level of the velocity fluctuations without using the high-pass filtering technique is demonstrated.

Description: Combinations of filters and subgrid scale stress models for large eddy simulation of the Navier-Stokes equations are examined by a priori tests and numerical simulations. The structure of the subgrid scales is found to depend strongly on the type of filter used, and consistency between model and filter is essential to ensure accurate results. The implementation of consistent combinations of filter and model gives more accurate turbulence statistics than those obtained in previous investigations in which the models were chosen independently from the filter. Results and limitations of the a priori test are discussed. The effect of grid refinement is also examined.

Description: Recent experiments have shown that, in rough-wall turbulent boundary layers, drag varies systematically with the spanwise aspect ratio lambda(z) (span/height) of roughness elements. In this paper, the effect of lambda(z) on turbulence structure has been examined. Based on lambda(z), the roughness in a transversely grooved surface with lambda(z) much greater than 1 is the opposite extreme of model plant canopies with lambda(z) much less than 1, studied in wind tunnels, whereas sandgrain is an intermediate type. Second-, third-, and fourth-order turbulence moments have been measured in turbulent boundary layers over transversely grooved and smooth surfaces and compared with available turbulence structure measurements over other types of surfaces. The near-wall turbulence structure is found to vary with lambda(z). The instantaneous motions involved in the flux of shear stress near the wall in smooth and transversely grooved surfaces are opposite in sign to those in three-dimensional roughness. The former is explained in terms of hairpin vortices alone, while the latter group is modeled to have an additional vortex (the so-called necklace vortex which straddles a three-dimensional roughness element near its base).

Description: The stability of inviscid swirling gas flows to small nonaxisymmetric perturbations is considered. For small Brunt-Vaisala frequencies, the problem reduces to the classical Sturm-Liouville form and the oscillation theorem can be applied. The resulting necessary and sufficient stability condition is compared to various criteria in the literature and a limited numerical study of isothermal rigidly rotating Poiseuille flow. For given azimuthal and axial wavenumbers, it is found numerically that the higher inertial modes become unstable for successively lower Rossby numbers and that this sequence of critical values approaches the theoretical value from above.

Description: The extension of the concept of vortex singularities, developed by Saffman and Meiron (1986) for the case of two-dimensional point vortices in an incompressible vortical flow, to the three-dimensional case of vortex sticks (vortons) is investigated analytically. The derivation of the governing equations is explained, and it is demonstrated that the formulation obtained conserves total vorticity and is a weak solution of the vorticity equation, making it an appropriate means for representing three-dimensional vortical flows with limited numbers of vortex singularities.

Description: Unarmed plastic projectiles can be equipped with small axisymmetric cavities for the generation of intense tones that are useful in training maneuvers. Attention is presently given to the simulation of rainfall in an airstream and the effect of rain droplet impingement on the nose of projectiles, and especially to any penetration or accumulation of water at the base of the cavity that might increase the fundamental cavity frequency and/or reduce the intensity of sound production during rain conditions.

Description: A numerical procedure in which the Navier-Stokes equations are discretized using tightly coupled discretizations of pressure derivatives and continuity equations is used here to extend the range of known terminal velocities of gaseous bubbles in liquids well beyond that in previous investigations. Computations performed for Reynolds numbers up to 2000 and Marangoni numbers up to 1000 show only a modest variation of the scaled bubble velocity between 0.16 and 0.5. The bubble velocity is influenced more by the Marangoni number than by the Reynolds number.

Description: A numerical analysis is performed on thermocapillary buoyancy convection induced by phase change in a liquid droplet. A finite-difference code is developed using an alternating-direction implicit (ADI) scheme. The intercoupling relation between thermocapillary force, buoyancy force, fluid property, heat transfer, and phase change, along with their effects on the induced flow patterns, are disclosed. The flow is classified into three types: thermocapillary, buoyancy, and combined convection. Among the three mechanisms, the combined convection simulates the experimental observations quite well, and the basic mechanism of the observed convection inside evaporating sessile drops is thus identified. It is disclosed that evaporation initiates unstable convection, while condensation always brings about a stable density distribution which eventually damps out all fluid disturbances. Another numerical model is presented to study the effect of boundary recession due to evaporation, and the 'peeling-off' effect (the removal of the surface layer of fluid by evaporation) is shown to be relevant.

Description: Two-dimensional solidification influenced by anisotropic heat conduction has been considered. The interfacial energy balance was derived to account for the heat transfer in one direction (x or y) depending on the temperature gradient in both the x and y directions. A parametric study was made to determine the effects of the Stefan number, aspect ratio, initial superheat, and thermal conductivity ratios on the solidification rate. Because of the imposed boundary conditions, the interface became skewed and sometimes was not a straight line between the interface position at the upper and lower adiabatic walls (spatially nonlinear along the height). This skewness depends on the thermal conductivity ratio k(yy)/k(yx). The nonlinearity of the interface is influenced by the solidification rate, aspect ratio, and k(yy/k(yx).

Description: A multiple-time-scale turbulence model of a single point closure and a simplified split-spectrum method is presented. In the model, the effect of the ratio of the production rate to the dissipation rate on eddy viscosity is modeled by use of the multiple-time-scales and a variable partitioning of the turbulent kinetic energy spectrum. The concept of a variable partitioning of the turbulent kinetic energy spectrum and the rest of the model details are based on the previously reported algebraic stress turbulence model. Example problems considered include: a fully developed channel flow, a plane jet exhausting into a moving stream, a wall jet flow, and a weakly coupled wake-boundary layer interaction flow. The computational results compared favorably with those obtained by using the algebraic stress turbulence model as well as experimental data. The present turbulence model, as well as the algebraic stress turbulence model, yielded significantly improved computational results for the complex turbulent boundary layer flows, such as the wall jet flow and the wake boundary layer interaction flow, compared with available computational results obtained by using the standard kappa-epsilon turbulence model.

Description: The emission from a gray radiating medium is analyzed for transient cooling in surroundings at a low temperature. The medium is rectangular with no variations in the direction normal to the cross section. The integral equation for the transient temperature distribution is solved numerically using a two-dimensional Gaussian integration subroutine. The emissive ability for a rectangle at uniform temperature is compared with that for transient cooling where the temperature distribution of the region has reached a fully developed shape, as shown by a separation of variables solution. The two solutions provide the upper and lower bounds for the emittance of a rectangle during transient cooling. The emittances for various aspect ratios are presented as a function of the mean length of the rectangle and are compared with results for a plane layer and a cylinder.

Description: Calculations for n-decane drops evaporating in a spherical cluster surrounded by unvitiated ambient air at atmospheric pressure were performed using two previously proposed cluster models. Both cluster models predict that turbulent transport effects are more important in the case of small clusters. This is due to the smaller volume to surface ratio and thus to the greater transport of hot unvitiated gas to the drops in order to promote evaporation. The results obtained are compared with those of two turbulent models for each one of the 'trapping factors' and similarity models.

Description: The effects of critical layer nonlinearity are considered on spatially growing oblique instability waves on nominally two-dimensional shear layers between parallel streams. The analysis shows that three-dimensional effects cause nonlinearity to occur at much smaller amplitudes than it does in two-dimensional flows. The nonlinear instability wave amplitude is determined by an integro-differential equation with cubic type nonlinearity. The numerical solutions to this equation are worked out and discussed in some detail. The numerical solutions always end in a singularity at a finite downstream distance.

Description: The effects of critical-layer nonlinearity on spatially growing oblique instability waves on compressible shear layers between two parallel streams are considered. The analysis shows that mean temperature nonuniformities cause nonlinearity to occur at much smaller amplitudes than it does when the flow is isothermal. The nonlinear instability wave growth rate effects are described by an integrodifferential equation which bears some resemblance to the Landau equation, in that it involves a cubic-type nonlinearity. The numerical solutions to this equation are worked out and discussed in some detail. Inviscid solutions always end in a singularity at a finite downstream distance, but viscosity can eliminate this singularity for certain parameter ranges.

Description: A simulation is performed of a passive scalar field convected by a rapidly fluctuating velocity field whose correlation time approaches zero. By using a code proposed in a previous study (Chasnov et al., 1988), the turbulence spectrum of the passive temperature field in the conductive subrange is determined. A theoretical model is proposed which explains the result obtained by representing the transfer of scalar variance by an eddy conductivity, whose correlation time is limited by the correlation time of the velocity field.

Description: Recent work on the equilibrium shapes of a rotating spheroid and drop shape oscillations is reviewed. The related theoretical aspects of the equilibrium shapes of a rotating drop, the stability, shape oscillations, gravitational forces, and drop fission are examined. Experimental findings on figures of a rotating drop in an immiscible system, drop oscillations in an immiscible system, oscillations of a rotating drop, compound drop oscillations, and drop dynamics in space are addressed.

Description: A formulation is presented by which any iteration process for obtaining the entropy increase in the flow of a one-dimensional steady nozzle is eliminated, and the simple solution of a quadratic equation is obtained. The proper parameters are then explicitly seen in the equation, and their effects on the solution are easily determined. Since only one root of the equation is physically admissible, entropy production, and therefore the shock wave, are uniquely determined by this set of parameters.

Description: Long's self-similar vortex is known to have two solutions for each supercritical value of the flow force. Each of those solutions is shown to have a double structure if the flow force is large. The inertial instabilities of one of those large-flow-force limit solutions are investigated, showing that they are related to the instabilities of the Bickley jet in one regime. However, the swirl in the vortex becomes important for long waves, very strongly modifying the sinuous and varicose, Bickley modes. The asymptotic results obtained agree well with the numerical solutions for the sinuous mode, but not for the varicose mode, the difficulty in the latter case being apparently due to mode jumping.

Description: Numerical techniques are developed to solve the Navier-Stokes equations for unsteady incompressible flow. The extension of the finite-difference Galerkin (FDG) method of Stephens et al. (1984) to the continuous-time case in two or three space dimensions is explained, and the numerical implementation of the method is discussed with particular attention to the staggered-MAC-grid primitive-variable discretization, the application of discrete mass balance to avoid problems inherent in FDG schemes, the direct interpretation of the FDG expansion variables as a discrete streamfunction, and a mass-balance approach to two-dimensional problems with throughflow or obstacles. Numerical results are presented graphically for the evolution of asymptotic steady flow in a driven cavity at Reynolds number 400, 1000, or 3200; good agreement with published experimental data is demonstrated, with accurate predictions of secondary-vortex formation from wall bubble recirculations at Reynolds number 1000.

Description: Flow-field measurements of unsteady turbulent flow downstream of a rotating spoked-wheel wake generator were performed in a short-duration light-piston tunnel, and the instantaneous-velocity data were phase averaged based on a signal synchronized with the bar-passing frequency. Mean axial velocities were found to agree well with those obtained from measurements behind a stationary cylinder and to be independent of both Reynolds and bar-passing Strouhal numbers. Reynolds stresses were found to be consistent with related cylinder-wake measurements, but were significantly higher than corresponding measurements obtained in large-scale research turbomachines. Phase-averaged triple velocity correlations were calculated from the digital velocity records, revealing the sign and the magnitude of skewness in the velocity probability density distributions for the two components.

Description: The bifurcation diagram corresponding to the Eckhaus stability curve has been constructed for the one-dimensional Swift-Hohenberg equation in a finite domain. Finite-amplitude solutions with particular spatial wavelength recover linear stability, as predicted by the Eckhaus curve, after a sequence of secondary bifurcations from the branch of solutions with this wavelength. No connectivity between the primary-solution branches is admissible if the stability predicted by this bifurcation diagram is to correspond to the prediction of the Eckhaus analysis. The Eckhaus curve does not exist if nonlinear couplings destroy this pattern. This is demonstrated by analysis of a coupled pair of Swift-Hohenberg equations.

Description: The effects of asymmetry in furnace temperature profile and pulling velocity on the crystal interface shape are demonstrated while neglecting the latent heat of solidification. It is seen that the furnace temperature profile may be varied in order to influence the shape of the melt-crystal interface. An exact thermal analysis is then performed on the Bridgman technique by including the latent heat of solidification as a source term. The exact temperature field required for yielding a flat melt-crystal interface is obtained. The earlier observation regarding the influence of furnace temperature profile on the interface shape is confirmed and a criterion for achieving a flat interface is obtained. Various furnace temperature profiles are selected and their corresponding melt-crystal interface results are presented.

Description: Transient cooling by radiation is analyzed for a cylindrical region filled with axially flowing streams of drops that are becoming solidified. This is of interest for the dissipation of waste heat from orbiting power system in space. The drops absorb, emit, and scatter radiation, and the surroundings are at a lower uniform temperature. The radiative properties are assumed gray, and the scattering is isotropic. The radiating region is a two-phase mixture that remains at the melting temperature of the drops. Its temperature uniformity maintains a high emissive power as energy is lost. This is an advantage over a sensible heat radiator in which the temperature decreases, thereby reducing the emissive power. The results provide the axial length that remains two-phase and the fraction of energy dissipated within this length in which the emissive power has not decreased because of sensible cooling. It is also shown how the radial distribution of the axial velocity of the drops can be modified to increase this energy fraction.

Description: A number of successful applications of a spectral collocation method extended by a multi-domain patching technique are shown. The multi-domain technique can be used to improve resolution for problems with widely disparate scales, and to reduce the ill-conditioning of the spectral operators for problems in which a large number of points are required for distributed resolution. A new nonreflecting outflow boundary treatment for unsteady transition-to-turbulence simulations is also presented, which relies on the multi-domain technique. The role of multi-domain in improving the efficiency of such calculations is discussed.

Description: Spectral element methods are high-order weighted residual techniques based on spectral expansions of variables and geometry for the Navier-Stokes (NS) and transport equations. Here, practical aspects of these methods and their efficient implementation are examined, and several examples of flows in truly complex geometries are presented. The spectral element discretization for NS equations is introduced, and the convergence of the method is addressed. An efficient data management scheme is discussed in the context of parallel processing computations. The method is validated by comparing the spectral element solutions with the exact eigensolutions for the Orr-Sommerfeld equations in two and three dimensions. Computer-aided flow visualizations are presented for an impulsive flow past a sharp edge wedge. Three-dimensional states of channel flow disrupted by an array of cylindrical eddy promoters are studied, and the results of a direct simulation of the turbulent flow in a plane channel are presented.

Description: The properties of the Karhunen-Loeve expansion of a strongly inhomogeneous random process are examined with emphasis on applications to turbulent flow fields. The ability of the KL expansion to represent functions that have both slow and rapid variations in a relatively small number of expansion terms is tested on a one-dimensional model based on the forced Burgers' equation. The rate of the convergence of the expansion is evaluated, and its dependence on the Reynolds number is determined. It is shown that the KL eigenfunctions possess wall boundary layers attached to outer structures that are independent of the Reynolds number (at high Reynolds numbers). It is also shown that the spectrum of eigenvalues is broad at large Reynolds numbers, requiring many terms to represent higher-order derivatives of the function.

Description: The flow in a two-dimensional curved channel driven by an azimuthal pressure gradient can become linearly unstable due to axisymmetric perturbations and/or nonaxisymmetric perturbations depending on the curvature of the channel and the Reynolds number. For a particular small value of curvature, the critical neighborhood of this curvature value and critical Reynolds number, nonlinear interactions occur between these perturbations. The Stuart-Watson approach is used to derive two coupled Landau equations for the amplitudes of these perturbations. The stability of the various possible states of these perturbations is shown through bifurcation diagrams. Emphasis is given to those cases which have relevance to external flows.

Description: To study thermocapillary flows in a two-dimensional thin liquid layer with heat fluxes imposed on the free surface experimentally, a long tray configuration was employed to simulate the infinite layer. The surface temperature distribution due to thermocapillary convection for different flow regimes was measured and compared with theorectical predictions. A short tray configuration was also employed to study the end wall effects (insulating or conducting). The results show that, for a strong convection flow with an insulating wall as the boundary, the surface temperature distribution became quite uniform. Consequently, the thermocapillary driving force was greatly reduced. On the other hand, a strong fluid motion always existed adjacent to the conducting wall because of the large surface temperature gradient near the wall.

Description: The 'transfinite element method' (TFEM) proposed by Tamma and Railkar (1987 and 1988) for the analysis of linear and nonlinear heat-transfer problems is described and demonstrated. The TFEM combines classical Galerkin and transform approaches with state-of-the-art FEMs to obtain a flexible hybrid modeling scheme. The fundamental principles of the TFEM and the derivation of the governing equations are reviewed, and numerical results for sample problems are presented in extensive graphs and briefly characterized. Problems analyzed include a square plate with a hole, a rectangular plate with natural and essential boundary conditions and varying thermal conductivity, the Space Shuttle thermal protection system, a bimaterial plate subjected to step temperature variations, and solidification in a semiinfinite liquid slab.

Description: The budgets for the Reynolds stresses and for the dissipation rate of the turbulence kinetic energy are computed using direct simulation data of a turbulent channel flow. The budget data reveal that all the terms in the budget become important close to the wall. For inhomogeneous pressure boundary conditions, the pressure-strain term is split into a return term, a rapid term, and a Stokes term. The Stokes term is important close to the wall. The rapid and return terms play different roles depending on the component of the term. A split of the velocity pressure-gradient term into a redistributive term and a diffusion term is proposed, which should be simpler to model. The budget data is used to test existing closure models for the pressure-strain term, the dissipation rate, and the transport rate. In general, further work is needed to improve the models.

Description: A cylindrical region filled with absorbing-emitting material is cooled by radiation to surroundings at a much lower temperature. A solution is found showing that, for each set of parameters, the transient radial temperature distribution reaches a fixed shape, although the temperatures are decreasing with time. This 'fully developed' transient region is characterized by having a constant emittance based on instantaneous values of the cylinder heat loss and mean temperature. This emittance depends only on the optical radius of the cylinder and the scattering albedo. The emittance is lower than that for a cylinder at uniform temperature. This arises from the larger local cooling and, hence, reduced temperatures of the outer layers of the cylinder. An examination of this transient emittance provides the ranges of parameters within which the simplification can be made that the cylinder has uniform radial temperature distribution throughout the cooling process.

Description: The nonlinear evolution of periodic disturbances on vortex trails is considered. In addition to following small initial perturbations, large amplitude initial disturbances of the vortex trails are also studied. It is shown that the equations support a rich variety of essentially nonlinear solutions including unbounded and quasisteady ones. These solutions are found to correspond to various modes of vortex clustering in the physical plane. At the close of the paper, comparisons of these results with recent numerical and experimental findings on the wakes behind stationary cylinders, and also transversely oscillating bluff objects, are made.

Description: An infinitely thin shear layer emanating from a semi-infinite flat plate subjected to acoustic excitation is considered. The flow field outside the excited shear layer is computed employing a source distribution approach. Results are given for the region of the velocity field that cannot easily be obtained by analytical approximations.

Description: Many problems in computational fluid dynamics (CFD) involve the calculation of flow fields within or around complex geometric configurations. The flow solution is computed on a computational grid. The construction of the grid is one of the major difficulties in the application of CFD to the analysis of flow about actual aircraft configurations. Due to geometric complexity, the grid has to be constructed in simple subregions and then all of these subgrids have to be pieced together to form a complete grid for the entire flow field. The entire grid, which is called a composite grid because it is formed from many parts, may have common regions. In either case, the computation of the flow field, using any numerical algorithm, will require the transfer of information between individual subgrids. The transfer of information is more difficult with overlapping grids. Algorithms have been developed and tested for automating the transfer of information between two overlapping grids.

Description: An effort directed at developing improved transitional models was initiated. The focus of this work was concentrated on the critical assessment of a popular existing transitional model developed by McDonald and Fish in 1972. The objective of this effort was to identify the shortcomings of the McDonald-Fish model and to use the insights gained to suggest modifications or alterations of the basic model. In order to evaluate the transitional model, a compressible boundary layer code was required. Accordingly, a two-dimensional compressible boundary layer code was developed. The program was based on a three-point fully implicit finite difference algorithm where the equations were solved in an uncoupled manner with second order extrapolation used to evaluate the non-linear coefficients. Iteration was offered as an option if the extrapolation error could not be tolerated. The differencing scheme was arranged to be second order in both spatial directions on an arbitrarily stretched mesh. A variety of boundary condition options were implemented including specification of an external pressure gradient, specification of a wall temperature distribution, and specification of an external temperature distribution. Overall the results of the initial phase of this work indicate that the McDonald-Fish model does a poor job at predicting the details of the turbulent flow structure during the transition region.

Description: A description of the Axial Flow Turbine Research Facility (AFTRF) being built at the Turbomachinery Laboratory of the Pennsylvania State University is presented. The purpose of the research to be performed in this facility is to obtain a better understanding of the rotor/stator interaction, three dimensional viscous flow field in nozzle and rotor blade passages, spanwise mixing and losses in these blade rows, transport of wake through rotor passage, and unsteady aerodynamics and heat transfer of rotor blade row. The experimental results will directly feed and support the analytical and the computational tool development. This large scale low speed facility is heavily instrumented with pressure and temperature probes and has provision for flow visualization and laser Doppler anemometer measurement. The facility design permits extensive use of the high frequency response instrumentation on the stationary vanes and more importantly on the rotating blades. Furthermore it facilitates detailed nozzle wake, rotor wake, and boundary layer surveys. The large size of the rig also has the advantage of operating at Reynolds numbers representative of the engine environment.

Description: A quasi-three-dimensional analysis has been developed for unsteady rotor-stator interaction in turbomachinery. The analysis solves the unsteady Euler or thin-layer Navier-Stokes equations in a body-fitted coordinate system. It accounts for the effects of rotation, radius change, and stress-surface thickness. The Baldwin-Lomax eddy-viscosity model is used for turbulent flows. The equations are integrated in time using an explicit four-stage Runge-Kutta scheme with a constant time step. Implicit residual smoothing is used to increase the stability limit of the time-accurate computations. The scheme is described, and stability and accuracy analyses are given.

Description: A 3-D model was developed for simulating multistage turbomachinery flows using supercomputers. This average passage flow model described the time averaged flow field within a typical passage of a bladed wheel within a multistage configuration. To date, a number of inviscid simulations were executed to assess the resolution capabilities of the model. Recently, the viscous terms associated with the average passage model were incorporated into the inviscid computer code along with an algebraic turbulence model. A simulation of a stage-and-one-half, low speed turbine was executed. The results of this simulation, including a comparison with experimental data, is discussed.

Description: A multiphase turbulence closure model is presented which employs one transport equation, namely the turbulence kinetic energy equation. The proposed form of this equation is different from the earlier formulations in some aspects. The power spectrum of the carrier fluid is divided into two regions, which interact in different ways and at different rates with the suspended particles as a function of the particle-eddy size ratio and density ratio. The length scale is described algebraically. A mass/time averaging procedure for the momentum and kinetic energy equations is adopted. The resulting turbulence correlations are modeled under less retrictive assumptions comparative to previous work. The closures for the momentum and kinetic energy equations are given. Comparisons of the predictions with experimental results on liquid-solid jet and gas-solid pipe flow show satisfactory agreement.

Description: Generally, two types of theory are used to describe the field equations for suspensions. The so-called postulated equations are based on the kinetic theory of mixtures, which logically should give reasonable equations for solutions. The basis for the use of such theory for suspensions is tenuous, though it at least gives a logical path for mathematical arguments. It has the disadvantage that it leads to a system of equations which is underdetermined, in a sense that can be made precise. On the other hand, the so-called averaging theory starts with a determined system, but the very process of averaging renders the resulting system underdetermined. A third type of theory is proposed in which the kinetic theory of gases is used to motivate continuum equations for the suspended particles. This entails an interpretation of the stress in the particles that is different from the usual one. Classical theory is used to describe the motion of the suspending medium. The result is a determined system for a dilute suspension. Extension of the theory to more concentrated systems is discussed.

Description: Detailed measurements were made of energy transfer among the scales of motion in incompressible turbulent fields at low Reynolds numbers generated by direct numerical simulation. It was observed that although the transfer resulted from triad interactions that were non-local in k space, the energy always transferred locally. The results are consistent with the notion of non-uniform advection of small weak eddies by larger and stronger ones, similar to transfer processes in the far dissipation range at high Reynolds numbers.

Description: Employing numerical simulations of 3-D homogeneous shear flow, the associated multifractal spectra of the energy dissipation, scalar dissipation and vorticity fields were calculated. The results for (128) cubed simulations of this flow, and those obtained in recent experiments that analyzed 1- and 2-D intersections of atmospheric and laboratory flows, are in some agreement. A two-scale Cantor set model of the energy cascade process which describes the experimental results from 1-D intersections quite well, describes the 3-D results only marginally.

Description: A detailed study of the intercomponent energy transfer processes by the pressure-strain-rate in homogeneous turbulent shear flow is presented. Probability density functions (pdf's) and contour plots of the rapid and slow pressure-strain-rate show that the energy transfer processes are extremely peaky, with high-magnitude events dominating low-magnitude fluctuations, as reflected by very high flatness factors of the pressure-strain-rate. A concept of the energy transfer class was applied to investigate details of the direction as well as magnitude of the energy transfer processes. In incompressible flow, six disjoint energy transfer classes exist. Examination of contours in instantaneous fields, pdf's and weighted pdf's of the pressure-strain-rate indicates that in the low magnitude regions all six classes play an important role, but in the high magnitude regions four classes of transfer processes, dominate. The contribution to the average slow pressure-strain-rate from the high magnitude fluctuations is only 50 percent or less. The relative significance of high and low magnitude transfer events is discussed.

Description: New definitions of entrainment and mixing based on the passive scalar field in the plane mixing layer are proposed. The definitions distinguish clearly between three fluid states: (1) unmixed fluid, (2) fluid engulfed in the mixing layer, trapped between two scalar contours, and (3) mixed fluid. The difference betwen (2) and (3) is the amount of fluid which has been engulfed during the pairing process, but has not yet mixed. Trends are identified from direct numerical simulations and extensions to high Reynolds number mixing layers are made in terms of the Broadwell-Breidenthal mixing model. In the limit of high Peclet number (Pe = ReSc) it is speculated that engulfed fluid rises in steps associated with pairings, introducing unmixed fluid into the large scale structures, where it is eventually mixed at the Kolmogorov scale. From this viewpoint, pairing is a prerequisite for mixing in the turbulent plane mixing layer.

Description: Several direct numerical simulations were performed and analyzed to study various aspects of the early development of mixing layers. Included are the phase jitter of the large-scale eddies, which was studied using a 2-D spatially-evolving mixing layer simulation; the response of a time developing mixing layer to various spanwise disturbances; and the sound radiation from a 2-D compressible time developing mixing layer.

Description: In the last few years, more and more evidence has emerged suggesting that transition to turbulence may be viewed as a succession of bifurcations to deterministic chaos. Most experimental and numerical observations have been restricted to Rayleigh-Benard convection and Taylor-Couette flow between concentric cylinders. An attempt is made to accurately describe the bifurcation sequence leading to chaos in a 2-D temporal free shear layer on the beta-plane. The beta-plane is a locally Cartesian reduction of the equations describing the dynamicss of a shallow layer of fluid on a rotating spherical planet. It is a valid model for large scale flows of interest in meteorology and oceanography.

Description: The vaporization of a droplet, interacting with its neighbors in a non-dilute spray environment is examined as well as a vaporization scaling law established on the basis of a recently developed theory of renormalized droplet. The interacting droplet consists of a centrally located droplet and its vapor bubble which is surrounded by a cloud of droplets. The distribution of the droplets and the size of the cloud are characterized by a pair-distribution function. The vaporization of a droplet is retarded by the collective thermal quenching, the vapor concentration accumulated in the outer sphere, and by the limited percolative passages for mass, momentum and energy fluxes. The retardation is scaled by the local collective interaction parameters (group combustion number of renormalized droplet, droplet spacing, renormalization number and local ambient conditions). The numerical results of a selected case study reveal that the vaporization correction factor falls from unity monotonically as the group combustion number increases, and saturation is likely to occur when the group combustion number reaches 35 to 40 with interdroplet spacing of 7.5 diameters and an environment temperature of 500 K. The scaling law suggests that dense sprays can be classified into: (1) a diffusively dense cloud characterized by uniform thermal quenching in the cloud; (2) a stratified dense cloud characterized by a radial stratification in temperature by the differential thermal quenching of the cloud; or (3) a sharply dense cloud marked by fine structure in the quasi-droplet cloud and the corresponding variation in the correction factor due to the variation in the topological structure of the cloud characterized by a pair-distribution function of quasi-droplets.

Description: Full turbulence simulations at low Reynolds numbers were made for the single-step, irreversible, bimolecular reaction between non-premixed reactants in isochoric, decaying homogeneous turbulence. Various initial conditions for the scalar field were used in the simulations to control the initial scalar dissipation length scale, and simulations were also made for temperature-dependent reaction rates and for non-stoichiometric and unequal diffusivity conditions. Joint probability density functions (pdf's), conditional pdf's, and various statistical quantities appearing in the moment equations were computed. Preliminary analysis of the results indicates that compressive strain-rate correlates better than other dynamical quantities with local reaction rate, and the locations of peak reaction rates seem to be insensitive to the scalar field initial conditions.

Description: Recent studies of turbulent shear flows have shown that many of their important kinematical and dynamical properties can be more clearly understood by describing the flows in terms of individual events or streamline patterns. These events or flow regions are studied because they are associated with relatively large contributions to certain average properties of the flow, for example kinetic energy, Reynolds stress, or to particular processes in the flow, such as mixing and chemical reactions, which may be concentrated at locations where streamlines converge for fast chemical reactions (referred to as convergence or C regions), or in recirculating eddying regions for slow chemical reactions. The aim of this project was to use the numerical simulations to develop suitable criteria for defining these eddying or vortical zones. The C and streaming (S) zones were defined in order to define the whole flow field. It is concluded that homogeneous and sheared turbulent flow fields are made up of characteristic flow zones: eddy, C, and S zones. A set of objective criteria were found which describe regions in which the streamlines circulate, converge or diverge, and form high streams of high velocity flow.

Description: The subharmonic resonance phenomenon in a spatially evolving mixing layer is studied using direct simulations of the 2-D Navier-Stokes equations. The computational domain extends to + or - infinity in the cross-stream direction with U(1) = 1.25 and U(2) = 0.25 imposed at + infinity and - infinity respectively. The domain is finite in the streamwise direction with inflow and outflow boundary conditions imposed at x/delta sub omega = 0 and 100, respectively. A hyperbolic-tangent mean velocity profile is assumed at the inlet and the Reynolds number based on the inlet vorticity thickness and velocity difference is Re = 600. It is observed that the phase angle between the fundamental and its subharmonic plays a key role in the spatial development of these modes. Contour plots of vorticity show that varying the phase will have a dramatic effect on the dynamics of the vortices. Pairing or shredding is observed depending on the phase. Fourier decomposition of the time traces show that the fundamental grows, saturates and decays with the downstream distance. The subharmonic has a similar behavior. However, the level at which the modes will saturate is affected by the phase. At 0 deg phase, it was found that as the fundamental saturates, the growth rate of the subharmonic is enhanced. At 90 deg phase, it was found that as the fundamental saturates, the growth rate of the subharmonic is inhibited. In the later case, the growth rate of the subharmonic recovers after saturation of the fundamental. These results are in qualitative agreement with experimental data.

Description: The numerical simulation of incompressible spatially-developing shear flows poses a special challenge to computational fluid dynamicists. The Navier-Stokes equations are elliptic and boundary equations need to be specified at the inflow and outflow boundaries in order to compute the fluid properties within the region of interest. It is, however, difficult to choose inflow and outflow conditions corresponding to a given experimental situation. Furthermore the effects that changes in the boundary conditions or in the size of the computational domain may induce on the global dynamics of the flow are presently unknown. These issues are examined in light of recent developments in hydrodynamic stability theory. The particular flow considered is the spatial mixing layer but it was expected that similar phenomena were bound to occur in other cases such as channel flow, the boundary layer, etc. A short summary of local/global and absolute/convective instability concepts is given. The results of numerical simulations are presented which strongly suggest that global resonances may be triggered in domains of finite streamwise extent although the evolution of the perturbation vorticity field is everywhere locally convective. A relationship between finite domains and pressure sources which might help in devising a scheme to eliminate these difficulties is discussed.

Description: A spatially-developing plane mixing layer was analyzed for chaotic behavior. A direct numerical simulation of the Navier-Stokes equations in a 2-D domain infinite in y and having inflow-outflow boundary conditions in x was used for data. Spectra, correlation dimension and the largest Lyapunov exponent were computed as functions of downstream distance x. When forced at a single (fundamental) frequency with maximum amplitude, the flow is periodic at the inflow but becomes aperiodic with increasing x. The aperiodic behavior is caused by the presence of a noisy subharmonic caused by the feedback between the necessarily nonphysical inflow and outflow boundary conditions. In order to overshadow this noise the flow was also studied with the same fundamental forcing and added random forcing of amplitude upsilon prime sub R/delta U = 0.01 at the inlet. Results were qualitatively the same in both cases: for small x, spectral peaks were sharp and dimension was nearly 1, but as x increased a narrowband spectral peak grew, spectra decayed exponentially at high frequencies and dimension increased to greater than 3. Based on these results, the flow appears to exhibit deterministic chaos. However, at no location was the largest Lyapunov exponent found to be significantly greater than zero.

Description: Experiments on wall-bounded shear flows (channel flows and boundary layers) have indicated that the turbulence in the region close to the wall exhibits a characteristic intermittently formed pattern of coherent structures. For a quantitative study of coherent structures it is necessary to make use of conditional sampling. One particularly successful sampling technique is the Variable Integration Time Averaging technique (VITA) first explored by Blackwelder and Kaplan (1976). In this, an event is assumed to occur when the short time variance exceeds a certain threshold multiple of the mean square signal. The analysis presented removes some assumptions in the earlier models in that the effects of pressure and viscosity are taken into account in an approximation based on the assumption that the near-wall structures are highly elongated in the streamwise direction. The appropriateness of this is suggested by the observations but is also self consistent with the results of the model which show that the streamwise dimension of the structure grows with time, so that the approximation should improve with the age of the structure.

Description: Effects of organized turbulent structures on the propagation of an optical beam in a homogeneous shear flow were studied. A passive-scalar field in a computed turbulent shear flow is used to represent index-of-refraction fluctuations, and phase errors induced in a coherent optical beam by turbulent fluctuations are computed. The organized vortical structures produce a scalar distribution with elongated regions of intense fluctuations which have an inclination with respect to the mean flow similar to that of the characteristic hairpin eddies. It is found that r.m.s. phase error is minimized by propagating approximately normal to the inclined vortical structures. Two-point correlations of vorticity and scalar fluctuation suggest that the regions of intense scalar fluctuation are produced primarily by the hairpin eddies.

Description: An analysis is presented of how a plane boundary affects the structure of turbulence in a sheared free stream. A uniform-shear boundary layer (USBL) is formulated with slip velocity condition at the surface, and inhomogeneous rapid distortion theory is applied. The effects of blocking by the surface on the turbulence structure in USBL is compared with those in the shear-free boundary layer (SFBL). Shear produces highly anisotropic eddies elongated in the flow direction. The vertical velocity variance is reduced with shear at all heights, roughly in proportion to the reduction in the homogeneous value, but the shape of the profile remains unchanged only near the surface. The streamwise integral scales increase with shear, indicating elongation of the streamwise extent of eddies.

Description: In wall-bounded shear flows the transition to turbulence through localized disturbances goes through a pattern starting with a development of shear layers. The localized normal velocity fluctuations induce normal vorticity through the lift-up effect. These shear layers become unstable to secondary disturbances, and if the amplitudes of the disturbances are large enough, a turbulent spot develops. Investigations of the spot in boundary layers has shown that the turbulent part of the spot is very similar to a fully developed boundary layer. Wygnanski et al. (1976) showed that the mean profile at the center-symmetry plane has a logarithmic region and Johansson et al. (1987) showed that both the higher-order statistics and flow structures in the spot were the same as in the corresponding fully developed flow. In what respects the turbulence inside the Poiseuille spot is similar to fully developed turbulent channel flow is studied. The numerically simulated spot is used, where the characteristics inside the spot are compared to those of the wave packet in the wingtip area. A recent experimental investigation of the velocity field associated with the Poiseuille spot by Klingmann et al. is used for comparison.

Description: Space-time evolution of near wall flow structures is described by conditional sampling methods, in which conditional averages are formed at various stages of development of shear layer structures. The development of spanwise asymmetry of the structures was found to be important in the creation of the structures and for the process of turbulence production.

Description: It is well known that turbulent mixing layers are dominated by large scale, fairly coherent structures, and that these structures are related to the stability characteristics of the flow. These facts have led researchers to attempt controlling such flows by selectively forcing certain unstable modes, which can in addition have the effect of suppressing other modes. Much of the work on controlling the mixing layer has relied on forcing 2-D instabilities. The results of forcing 3-D instabilities are addressed. The objectives of the work are twofold: to understand how a mixing layer responds to 3-D perturbations, and to test the validity of an amplitude expansion in predicting the mixing layer development. The amplitude expansion could be very useful in understanding and predicting the 3-D response of the flow to a variety of initial conditions.

Description: Data obtained from the numerical simulation of a 2-D mixing layer were used to study the feasibility of using the instantaneous concentration of a passive scalar for detecting the typical coherent structures in the flow. The study showed that this technique works quite satisfactorily and yields results similar to those that can be obtained by using the instantaneous vorticity for structure detection. Using the coherent events educed by the scalar conditioning technique, the contribution of the coherent events to the total turbulent momentum and scalar transport was estimated. It is found that the contribution from the typical coherent events is of the same order as that of the time-mean value. However, the individual contributions become very large during the pairing of these structures. The increase is particularly spectacular in the case of the Reynolds shear stress.

Description: An incompressible, time developing 3-D mixing layer with idealized initial conditions was simulated numerically. Consistent with the suggestions from experimental measurements, the braid region between the dominant spanwise vortices or rolls develops longitudinal vortices or ribs, which are aligned upstream and downstream of a roll and produce spanwise distortion of the rolls. The process by which this distortion occurs is explained by studying a variety of quantities of dynamic importance (e.g., production of enstrophy, vortex stretching). Other quantities of interest (dissipation, helicity density) are also computed and discussed. The currently available simulation only allows the study of the early evolution (before pairing) of the mixing layer. New simulations in progress will relieve this restriction.

Description: To unravel the liquid-jet breakup process in the nondilute region, a newly developed system of real-time X-ray radiography, an advanced digital image processor, and a high-speed video camera were used. Based upon recorded X-ray images, the inner structure of a liquid jet during breakup was observed. The jet divergence angle, jet breakup length, and fraction distributions along the axial and transverse directions of the liquid jets were determined in the near-injector region. Both wall- and free-jet tests were conducted to study the effect of wall friction on the jet breakup process.

Description: A review of the application of single particle hydrodynamics in models for the exchange of interphase momentum in continuum models of multiphase flow is presented. Considered are the equations of motion for a laminar, mechanical two phase flow. Inherent to this theory is a model for the interphase exchange of momentum due to drag between the dispersed particulate and continuous fluid phases. In addition, applications of two phase flow theory to de-mixing flows require the modeling of interphase momentum exchange due to lift forces. The applications of single particle analysis in deriving models for drag and lift are examined.