ALBERT

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: 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: 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: 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: 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: 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: The capability of accurate nonlinear flow analysis of resonance systems is essential in many problems, including combustion instability. Classical numerical schemes are either too diffusive or too dispersive especially for transient problems. In the last few years, significant progress has been made in the numerical methods for flows with shocks. The objective was to assess advanced shock capturing schemes on transient flows. Several numerical schemes were tested including TVD, MUSCL, ENO, FCT, and Riemann Solver Godunov type schemes. A systematic assessment was performed on scalar transport, Burgers' and gas dynamic problems. Several shock capturing schemes are compared on fast transient resonant pipe flow problems. A system of 1-D nonlinear hyperbolic gas dynamics equations is solved to predict propagation of finite amplitude waves, the wave steepening, formation, propagation, and reflection of shocks for several hundred wave cycles. It is shown that high accuracy schemes can be used for direct, exact nonlinear analysis of combustion instability problems, preserving high harmonic energy content for long periods of time.

Description: An explicit multistage Runge-Kutta type of time-stepping scheme is used for solving transonic flow past a transport type wing/fuselage configuration. Solutions for both Euler and Navier-Stokes equations are obtained for quantitative assessment of boundary layer interaction effects. The viscous solutions are obtained on both a medium resolution grid of approximately 270,000 points and a find grid of 460,000 points to assess the effects of grid density on the solution. Computed pressure distributions are compared with the experimental data.

Description: Some considerations toward developing numerical procedures for simulating viscous compressible flows are discussed. Both Navier-Stokes and boundary layer field methods are considered. Because efficient viscous-inviscid interaction methods have been difficult to extend to complex 3-D flow simulations, Navier-Stokes procedures are more frequently being utilized even though they require considerably more work per grid point. It would seem a mistake, however, not to make use of the more efficient approximate methods in those regions in which they are clearly valid. Ideally, a general purpose compressible flow solver that can optionally take advantage of approximate solution methods would suffice, both to improve accuracy and efficiency. Some potentially useful steps toward this goal are described: a generalized 3-D boundary layer formulation and the fortified Navier-Stokes procedure.

Description: The construction and development of the multi-component traversing system and associated control hardware and software are presented. A hydrogen bubble/laser sheet flow visualization technique was developed to visually study the characteristics of the mixing layers. With this technique large-scale rollers arising from the Taylor-Gortler instability and its interaction with the primary Kelvin-Helmholtz structures can be studied.

Description: The sizes and arrangement of the wind tunnel used for the experimentation are described. The specifications for the cold-wire anemometers, hot-wire anemometers, cold-wire rakes, and miniature 3-wire probe are proveded. The results of the experiment are briefly discussed.

Description: A very low Reynolds number turbulent boundary layer subject to an adverse pressure gradient is studied. The aim is to obtain highly accurate mean-flow and turbulence measurements under conditions that can be closely related to the numerical simulations of Philippe Spalart for the purposes of CFD validation. Much of the Boundary Layer Wind Tunnel was completely rebuilt with a new wider contraction and working section which will improve compatibility with the simulations. A unique sophisticated high-speed computer controlled 3-D probe traversing mechanism was integrated into the test section. Construction of the tunnel and traverse is discussed in some detail. The hardware is now complete, and measurements are in progress. The mean-flow data indicate that a suitably two-dimensional base flow was established. Automation of the probe positioning and data acquistion have led to a decreased running time for total pressure measurements. However, the most significant benefits are expected to occur when using hot-wire probes. Calibrations can be performed automatically and there is no need to handle fragile probes when moving between measuring stations. Techniques are being developed which require sampling of the signals from moving hot-wire probes on the basis of their position in the flow. Measurements can be made in high intensity turbulence by flying probes upstream at high speed so that the relative magnitude of the turbulent velocity fluctuations are reduced. In regions, where the turbulence intensity is not too large, the probe can also be repetitively scanned across very dense spatial grids in other directions. With this technique, a complete profile can be measured in about 1/3 the time and with a spatial density about 50 times that obtainable using a stationary probe.

Description: Thermal convection was proposed as a possible mechanism for generation and maintenance of turbulence in the inner accretion disk regime of the primordial solar nebula. It is of fundamental interest to design experiments with the basic physical features of the solar nebula conditions cannot be produced in the laboratory, numerical simulations of hydrodynamic flows, which have been very successful in describing aerodynamic flows, can be suitable modified to provide experimental data for solar nebula modelling. The goals are to modify an extant, proven hydrodynamics code with the most important features of the solar nebula and other thin accretion disks: bouyancy terms to generate convection, internal heating representing the release of gravitational potential energy, a variable gravity linearly proportional the the distance from the vertical midplane due to centrifugal balance, rapid rotation with axis aligned with gravity, and Keplerian rotational shear; to determine the effect that these features have on the turbulent convection by introducing them individually and to determine the cumulative nature of the turbulent convection for accretion disk conditions; and to model the convection and the turbulence. In this manner, prior solar nebula models can be tested and their deficiencies rectified.

Description: It is likely that turbulence played a major role in the evolution of the solar nebula, which is the flattened disk of dust and gas out of which our solar system formed. Relevant turbulent processes include the transport of angular momentum, mass, and heat, which were critically important to the formation of the solar system. This research will break ground in the modeling of compressible turbulence and its effects on the evolution of the solar nebula. The computational techniques which were developed should be of interest to researchers studying other astrophysical disk systems (e.g., active galactic nuclei), as well as turbulence modelers outside the astrophysics community.

Description: A research program for direct numerical simulations of compressible reacting flows is described. Two main research subjects are proposed: the effect of pressure waves on turbulent combustion and the use of direct simulation methods to validate flamelet models for turbulent combustion. The interest of a compressible code to study turbulent combustion is emphasized through examples of reacting shear layer and combustion instabilities studies. The choice of experimental data to compare with direct simulation results is discussed. A tentative program is given and the computation cases to use are described as well as the code validation runs.

Description: Outline of the research program and a recent progress in the studies of sheared turbulence are described. The research program reported is directed at two goals: (1) understanding of fundamental physics of organized structures in turbulent shear flows; and (2) development of phenomenological models of turbulence based on physical arguments. Three projects that were carried out are: (1) structure of sheared turbulence near a plane boundary; (2) distortion of turbulence by axisymmetric strain and dilation; and (3) study of energy transfer in turbulent shear flow.

Description: Vortex interactions and their role in turbulent flow are examined. The objectives are twofold. First, to use the existing axisymmetric code to study the annihilation process of colliding vortex rings and determine the relevance of this problem to similar 3-D phenomena. The second objective is to extend the code to three dimensions. The code under development is unique in that it can compute flows in a truly infinite domain (i.e., without periodic boundary conditions or approximations from truncating the domain). Because of this, the far field sound can be computed, and therefore, contribute to improved models of turbulence generated noise for this class of flows. Issues which can be addressed by the code include: effects of viscosity on mode selection in azimuthal breakdown of vortex rings (i.e., the Widnall instability); reconnection, the associated production of small scales, and the time scale of the process.

Description: The motivation for studying close vortex interactions is briefly discussed in the light of turbulence and coherent structures. Particular attention is given to the interaction known as reconnection. Two reconnection mechanisms are discussed. One is annihilation of vorticity by cross-diffusion, the other is an inviscid head-tail formation. At intermediate Reynolds numbers both mechanisms are operating.

Description: The objective of the present research was to extend the Direct Numerical Simulation (DNS) approach to particle-laden turbulent flows using a simple model of particle/flow interaction. The program addressed the simplest type of flow, homogeneous, isotropic turbulence, and examined interactions between the particles and gas phase turbulence. The specific range of problems examined include those in which the particle is much smaller than the smallest length scales of the turbulence yet heavy enough to slip relative to the flow. The particle mass loading is large enough to have a significant impact on the turbulence, while the volume loading was small enough such that particle-particle interactions could be neglected. Therefore, these simulations are relevant to practical problems involving small, dense particles conveyed by turbulent gas flows at moderate loadings. A sample of the results illustrating modifications of the particle concentration field caused by the turbulence structure is presented and attenuation of turbulence by the particle cloud is also illustrated.

Description: Direct numerical simulations are being performed on two different fluid flows in an attempt to discover the mechanism underlying the transition to turbulence in each. The first system is Taylor-Couette flow; the second, two-dimensional flow over an airfoil. Both flows exhibit a gradual transition to high-dimensional turbulence through low-dimensional chaos. The hope is that the instabilities leading to chaos will be easier to relate to physical processes in this case, and that the understanding of these mechanisms can then be applied to a wider array of turbulent systems.

Description: Cocke (1969) proved that in incompressible, isotropic turbulence the average material line (material surface) elements increase in comparison with their initial values. Good estimates of how much they increase in terms of the eigenvalues of the Green deformation tensor were rigorously obtained.

Description: The objective is to understand and extend a recent theory of turbulence based on dynamic renormalization group (RNG) techniques. The application of RNG methods to hydrodynamic turbulence was explored most extensively by Yakhot and Orszag (1986). An eddy viscosity was calculated which was consistent with the Kolmogorov inertial range by systematic elimination of the small scales in the flow. Further, assumed smallness of the nonlinear terms in the redefined equations for the large scales results in predictions for important flow constants such as the Kolmogorov constant. It is emphasized that no adjustable parameters are needed. The parameterization of the small scales in a self-consistent manner has important implications for sub-grid modeling.

Description: Use of the Smagorinsky eddy-viscosity formulation and related schemes for subgrid-scale parameterization of large eddy simulation models requires specification of a single length scale, earlier related by Lilly to the scale of filtering and/or numerical resolution. An anisotropic integration of the Kolmogoroff enstrophy spectrum allows generalization of that relationship to anisotropic resolution. It is found that the Deardorff assumption is reasonably accurate for small anisotropies and can be simply improved for larger values.

Description: In recent years codes that use the Navier-Stokes equations to compute aerodynamic flows have evolved from computing two-dimensional flows around simple airfoils to computing flows around full scale aircraft configurations. Most flows of engineering interest are turbulent and turbulence models are needed for their prediction. Yet, it is known that present turbulence models are adequate only for simple flows and do poorly in complicated flows such as three-dimensional separation, or large-scale unsteadiness. The same progress that allowed the development of these aerodynamic codes, namely the introduction of supercomputers, has allowed us to compute directly turbulent flows, albeit only for simple flows at moderate Reynolds numbers. These direct turbulence simulations provide us with detailed data that experimentalists were not able to measure. This work is motivated by the fact that data exists for developing better turbulence models and by the need for better models to compute flows of engineering interest. The objective is to develop turbulence models for engineering applications. The model categories that show promise for immediate use are on the two-equation level and the Reynolds-stress level.

Description: In the last two decades there have been extensive developments in computational unsteady transonic aerodynamics. Such developments are essential since the transonic regime plays an important role in the design of modern aircraft. Therefore, there has been a large effort to develop computational tools with which to accurately perform flutter analysis at transonic speeds. In the area of Computational Fluid Dynamics (CFD), unsteady transonic aerodynamics are characterized by the feature of modeling the motion of shock waves over aerodynamic bodies, such as wings. This modeling requires the solution of nonlinear partial differential equations. Most advanced codes such as XTRAN3S use the transonic small perturbation equation. Currently, XTRAN3S is being used for generic research in unsteady aerodynamics and aeroelasticity of almost full aircraft configurations. Use of Euler/Navier Stokes equations for simple typical sections has just begun. A brief history of the development of CFD for aeroelastic applications is summarized. The development of unsteady transonic aerodynamics and aeroelasticity are also summarized.

Description: A finite difference code was implemented for the compressible Navier-Stokes equations on the Connection Machine, a massively parallel computer. The code is based on the ARC2D/ARC3D program and uses the implicit factored algorithm of Beam and Warming. The codes uses odd-even elimination to solve linear systems. Timings and computation rates are given for the code, and a comparison is made with a Cray XMP.

Description: The hardware and software currently used for visualization of fluid dynamics at NASA Ames is described. The software includes programs to create scenes (for example particle traces representing the flow over an aircraft), programs to interactively view the scenes, and programs to control the creation of video tapes and 16mm movies. The hardware includes high performance graphics workstations, a high speed network, digital video equipment, and film recorders.

Description: The present investigation has focused on a computational methodology for the fundamental case of transition in channel flow, in which recently published experimental data are utilized both as a stimulus and as a measure of merit of the method. The research has proceeded along three avenues in parallel. The first task has consisted of the development and verification of a computer code which calculates the mean evolution of flow in a channel similar to the one employed experimentally by Blair and Anderson. An analytical test case was created for the dual purposes of code verification and of highlighting the interactions between the Reynolds stress and the mean velocity profile. This test case generated a Reynolds stress by the residue in the momentum equation which is produced by a typical analytical velocity profile. By a substitution of this Reynolds stress into the appropriate code module, the correctness of the code may be verified, along with the accuracy of the computational method. The second task pursued has involved the development of a triple layer model for the Reynolds stress profile, which was suggested and derived from experimental velocity profiles. It is demonstrated that the innermost length scale is based on the local friction velocity, the intermediate layer corresponds to the usual logarithmic law of the wall region in which the normalized Reynolds stress is approximately unity, and the outermost layer is represented by a closed mathematical form depending explicitly on the velocity profile in the wake region. The third task was comprised of scrutiny of the excellent databases developed by Blair and others, and the planning of its incorporation into the transition analysis. These extensive measurements indicate that turbulent statistics in the transition regime may be considered to alternate between laminar and fully turbulent types, the proportions of which are quantified by a measured intermittency function.

Description: Numerous experimental studies were conducted on the steady, three-dimensional boundary layer over a disk rotating at constant angular speed in an otherwise undisturbed fluid. The subject flow geometry is of interest because it provides a relatively simple way to study the cross-flow instability phenomenon which occurs in three-dimensional boundary layers, as on swept wings. This flow instability results in the formation of a stationary spiral vortex flow field over the disk, as shown by Wilkinson and Malik. Using a hot-wire probe, the spatial wave pattern of stationary vortices, which filled the entire circumference of the disk was mapped. The subject flow instability caused transition-to-turbulent flow as the periphery of the disk was approached. The effect on receptivity and transition of discrete disturbance modes, such as three-dimensional toughness elements and acoustic excitation was investigated. The present study (an extension of the work of Wilkinson and Malik) is focused on the effect of pulsed point suction on flow instability and transition, and consequently, on the classical stationary vortical flow pattern.

Description: In 1978, the Russian mathematician V. Kharitonov published a remarkably simple necessary and sufficient condition in order that a rectangular parallelpiped of polynomials be a stable set. Here, stable is taken to mean that the polynomials have no roots in the closed right-half of the complex plane. The possibility of generalizing this result was studied by numerous authors. A set, Q, of polynomials is given and a necessary and sufficient condition that the set be stable is sought. Perhaps the most general result is due to Barmish who takes for Q a polytope and proceeds to construct a complicated nonlinear function, H, of the points in Q. With the notion of stability which was adopted, Barmish asks that the boundary of the closed right-half plane be swept, that the set G is considered = to (j(omega)(bar) - infinity is less than omega is less than infinity) and for each j(omega)(sigma)G, require H(delta) is greater than 0. Barmish's scheme has the merit that it describes a true generalization of Kharitonov's theorem. On the other hand, even when Q is a polyhedron, the definition of H requires that one do an optimization over the entire set of vertices, and then a subsequent optimization over an auxiliary parameter. In the present work, only the case where Q is a polyhedron is considered and the standard definition of stability described, is used. There are straightforward generalizations of the method to the case of discrete stability or to cases where certain root positions are deemed desirable. The cases where Q is non-polyhedral are less certain as candidates for the method. Essentially, a method of geometric programming was applied to the problem of finding maximum and minimum angular displacements of points in the Nyquist locus (Q(j x omega)(bar) - infinity is less than omega is less than infinity). There is an obvious connection with the boundary sweeping requirement of Barmish.

Description: Recently, NASA, FAA, and other organizations have focused their attention upon the possible effects of rain on airfoil performance. Rhode carried out early experiments and concluded that the rain impacting the aircraft increased the drag. Bergrum made numerical calculation for the rain effects on airfoils. Luers and Haines did an analytic investigation and found that heavy rain induces severe aerodynamic penalties including both a momentum penalty due to the impact of the rain and a drag and lift penalty due to rain roughening of the airfoil and fuselage. More recently, Hansman and Barsotti performed experiments and declared that performance degradation of an airfoil in heavy rain is due to the effective roughening of the surface by the water layer. Hansman and Craig did further experimental research at low Reynolds number. E. Dunham made a critical review for the potential influence of rain on airfoil performance. Dunham et al. carried out experiments for the transport type airfoil and concluded that there is a reduction of maximum lift capability with increase in drag. There is a scarcity of published literature in analytic research of two-phase boundary layer around an airfoil. Analytic research is being improved. The following assumptions are made: the fluid flow is non-steady, viscous, and incompressible; the airfoil is represented by a two-dimensional flat plate; and there is only a laminar boundary layer throughout the flow region. The boundary layer approximation is solved and discussed.

Description: Researchers started their studies on the development and application of computational methods for compressible flows. Particular attention was given to proper numerical treatment of sharp layers occurring in such problems and to general mesh generation capabilities for intricate computational geometries. Mainly finite element methods enhanced with several state-of-the art techniques (such as the streamline-upwind/Petrov-Galerkin, discontinuity capturing, adaptive implicit-explicit, and trouped element-by-element approximate factorization schemes) were employed.

Description: A self-consistent derivation of the conservation laws is given for flows of a fluid-solid mixture. A unified analytical framework for obtaining constitutive relations is provided. This analysis uses a control volume/control surface approach that is widely used in fluid mechanics. All terms in the governing equations and the constitutive relations are written in terms of the mass-weighted averages except solid concentration. It is believed that the mass-weighted average is the natural bridge between micromechanics and constitutive relations. The derived momentum equations contain terms that differ from all existing models except that of Prosperetti and Jones (1984). However, their assumptions are not needed here. Special attention is given to the solid phase pressure. The physical basis of the previously assumed form for this pressure (Givler 1987) becomes clear. A number of related phenomena are also discussed. These include the anti-diffusion and anisotropic normal stresses. The energy equations are also different from existing models.

Description: Some features of two recent approaches of two-phase potential flow are presented. The first approach is based on a set of progressive examples that can be analyzed using common techniques, such as conservation laws, and taken together appear to lead in the direction of a general theory. The second approach is based on variational methods, a classical approach to conservative mechanical systems that has a respectable history of application to single phase flows. This latter approach, exemplified by several recent papers by Geurst, appears generally to be consistent with the former approach, at least in those cases for which it is possible to obtain comparable results. Each approach has a justifiable theoretical base and is self-consistent. Moreover, both approaches appear to give the right prediction for several well-defined situations.

Description: The main characteristics and the potential advantages of generalized drift flux models are presented. In particular it is stressed that the issue on the propagation properties and on the mathematical nature (hyperbolic or not) of the model and the problem of closure are easier to tackle than in two fluid models. The problem of identifying the differential void-drift closure law inherent to generalized drift flux models is then addressed. Such a void-drift closure, based on wave properties, is proposed for bubbly flows. It involves a drift relaxation time which is of the order of 0.25 s. It is observed that, although wave properties provide essential closure validity tests, they do not represent an easily usable source of quantitative information on the closure laws.

Description: Numerical techniques for solving the compressible Euler and Navier-Stokes equations are discussed with an emphasis on characteristic-based schemes. Two popular approaches, flux difference splitting and flux vector splitting, are described in one-dimensional Cartesian coordinates and then extended to three-dimensional generalized coordinates. A technique for increasing the spatial accuracy is presented, followed by a discussion of numerical dissipation mechanisms. An introduction to the use of implicit time integration schemes for accelerating the convergence rate to steady-state solutions including Newton's method, relaxation strategies, and approximate factorization techniques and their implementation on a vector processor concludes the chapter.

Description: Rectangular jets injected from a flat plate into a crossflow at large angles have been studied. Results were obtained as surface pressure distributions, mean velocity vector plots, turbulence intensities, and Reynolds stresses in the jet plume. The length-to-width ratio of the jets was 4, and the jets were aligned streamwise as single and side-by-side dual jets. The jet injection angles were 90 and 60 deg. Surface pressure distribution results were obtained for jet-to-freestream velocity ratios of 2.2, 4, and 8. Mean flow and turbulence flowfield data were obtained for the side-by-side dual jets, mainly for the jet-to-freestream velocity ratio of 4. The jets featured strong negative pressure peaks near the front nozzle corners. The 60-deg jets produced lower magnitude negative pressures, which are distributed over a lesser area when compared to the 90-deg jets.

Description: A comparatively simple procedure is presented for the direct summation of the velocity field introduced by point vortices which significantly reduces the required number of operations by replacing selected partial sums by asymptotic series. Tables are presented which demonstrate the speed of this algorithm in terms of the mere doubling of computational time in dealing with a doubling of the number of vortices; current methods involve a computational time extension by a factor of 4. This procedure need not be restricted to the solution of the Poisson equation, and may be applied to other problems involving groups of points in which the interaction between elements of different groups can be simplified when the distance between groups is sufficiently great.

Description: The phenomenon of hyperbolic heat conduction in contrast to the classical (parabolic) form of Fourier heat conduction involves thermal energy transport that propagates only at finite speeds, as opposed to an infinite speed of thermal energy transport. To accommodate the finite speed of thermal wave propagation, a more precise form of heat flux law is involved, thereby modifying the heat flux originally postulated in the classical theory of heat conduction. As a consequence, for hyperbolic heat conduction problems, the thermal energy propagates with very sharp discontinuities at the wave front. Accurate solutions are found for a class of one-dimensional hyperbolic heat conduction problems involving non-Fourier effects that can be used effectively for representative benchmark tests and for validating alternate schemes. Modeling/analysis formulations via specially tailored hybrid computations are provided for accurately modeling the sharp discontinuities of the propagating thermal wave front. Comparative numerical test models are presented for various hyperbolic heat conduction models involving non-Fourier effects to demonstrate the present formulations.

Description: The spatial distribution of the numerical disturbances that are generated during the numerical solution of a flow is examined. It is shown that the distribution of the disturbances is not uniform. In regions where the structure of a flow is simple, the magnitudes of the generated disturbances is small and their decay is fast. However, in complex flow regions, as in separation and vortical areas, large magnitude disturbances appear and their decay may be very slow. The observed nonuniformity of the numerical disturbances makes possible the reduction of the calculation time by application of what may be called the partial-grid calculation technique, in which a major part of the calculation procedure is applied in selective subregions, where the velocity disturbances are large, and not within the whole grid. This technique is expected to prove beneficial in large-scale calculations such as the flow about complete aircraft configurations at high angle of attack. Also, it has been shown that if the Navier-Stokes equations are written in a generalized coordinate system, then in regions in which the grid is fine, such as near solid boundaries, the norms become infinitesimally small, because in these regions the Jacobian has very large values. Thus, the norms, unless they are unscaled by the Jacobians, reflect only the changes that happen at the outer boundaries of the computation domain, where the value of the Jacobian approaches unity, and not in the whole flow field.

Description: Direct numerical simulations of the unsteady incompressible Navier-Stokes equations have been performed in order to investigate the behavior of passive-scalar fields resulting from mean scalar gradients in each of three orthogonal directions in homogeneous turbulent shear flow. For all orientations of the mean scalar gradient, the sum of the pressure-scalar gradient and velocity gradient-scalar gradient terms in the turbulent scalar flux balance equation are found to be approximately aligned with the scalar flux vector itself. The simulation results are used to obtain dimensionless model coefficients as a function of the turbulence Reynolds and Peclet numbers.

Description: The theoretical model of Lee and Wang (1986) for the instability of an annular jet, in which the jet's liquid layer is treated as a thin liquid sheet, is examined. It is suggested that the model should be altered so that when the envelope is closing its bottleneck during collapse, the new envelope experiences a sharp pressure pulse from its gaseous core, reversing the normal velocity of the sheet enough to maintain continuous constant gas flow. Using this improved version of the model, it is shown that if the liquid velocity is high enough and the gas velocity is greater than the liquid velocity, the bubble-formation frequency varies linearly with the difference between the two velocities, but not with their individual values.

Description: The vortical evolution of mixing layers subject to various types of forcing is numerically simulated using pseudospectral methods. The effect of harmonic forcing and random noise in the initial conditions is examined with some results compared to experimental data. Spanwise forcing is found to enhance streamwise vorticity in a nonlinear process leading to a slow, secondary growth of the shear layer. The effect of forcing on a chemical reaction is favorably compared with experimental data at low Reynolds numbers. Combining harmonic and subharmonic forcing is shown to both augment and later destroy streamwise vorticity.

Description: The nonlinear interaction between sinuous and varicose instability modes in a plane wake is examined in the nonlinear-nonequilibrium critical layer regime. Equations governing the evolution of the instability wave amplitudes and critical layer vorticity distributions are derived. Numerical solutions for these equations are obtained for a number of wake defects and initial amplitude ratios. The results show that the primary effects of the nonlinear interaction are the suppression of the varicose mode and the downstream shift of the peak of the sinuous mode.

Description: The transient cooling of an absorbing, emitting, and scattering cylinder is studied by using the forward integration of energy transfer equations. The analysis presented here provides cooling curves for a cylindrical region, initially at uniform temperature and then suddenly exposed to a much cooler environment. Results are presented for several optical radii and different values of the scattering albedo.

Description: The present treatment of the early stages of boundary layer-transition phenomena, where the unsteady motion is of small amplitude and can be accordingly treated as a small perturbation of an appropriate mean flow, elaborates the Heinrich et al. (1988) discussion of the role played by this 'receptivity' stage: in which the unsteady flow exhibits the same harmonic time-dependence as the externally-imposed forcing. Freestream disturbance wavelengths are noted to often be much longer than the Tollmien-Schlichting wavelength. Attention is given to the variety of wavelength-reduction mechanisms able to couple the long-wavelength, freestream disturbances to the comparatively short Tollmien-Schlichting waves.

Description: Most of the detailed turbulence-structure data available pertain only to the simplest cases, involving zero pressure-gradient boundary layers and free-shear layers, and indicate that each disparate geometry possesses its own set of dominant nonlinear instabilities. Various boundary/input conditions act to modify these instabilities for low input levels; for stronger inputs, the basic instability modes/structures sustaining the turbulence field may be altered. Steady-state inputs are noted to be extremely effective in altering turbulence structures, in the directions of either amplification or diminution.

Description: The liquid drop radiator consists of many directed streams of hot liquid drops that are cooled by passing through space and are then collected for reuse. To facilitate the collection of the cooled liquid, a converging geometry might be useful. As the streams converge toward the collector, the density and optical thickness of the droplet cloud are increased. The increase of optical thickness in the flow direction is shown to reduce the local emittance along the length of the radiator.

Description: It is demonstrated here that inhomogeneous diffusion-limited aggregation (DLA) model can be used to simulate viscous fingering in a medium with inhomogeneous permeability and homogeneous porosity. The medium consists of a pipe-pore square-lattice network in which all pores have equal volume and the pipes have negligible volume. It is shown that fluctuations in a DLA-based growth process may be tuned by noise reduction, and that fluctuations in the velocity of the moving interface are multiplicative in form.

Description: By application of simple computer graphics techniques, the statistical performance of two Monte Carlo methods used in the simulation of rarefied gas flows are assessed. Specifically, two direct simulation Monte Carlo (DSMC) methods developed by Bird and Nanbu are considered. The graphics techniques are found to be of great benefit in the reduction and interpretation of the large volume of data generated, thus enabling important conclusions to be drawn about the simulation results. Hence, it is discovered that the method of Nanbu suffers from increased statistical fluctuations, thereby prohibiting its use in the solution of practical problems.

Description: A highly conducting charged drop that is surrounded by a fluid insulator of another density can be levitated by suitably applying a uniform electric field. Axisymmetric equilibrium shapes and stability of the levitated drop are found by solving simultaneously the augmented Young-Laplace equation for surface shape and the Laplace equation for the elecric field, together with constraints of fixed drop volume, charge, and center of mass. The means are a method of subdomains, finite element basis functions, and Galerkin's method of weighted residuals, all facilitated by a large-scale computer. Shape families of fixed charge are treated systematically by first-order continuation. Previous analyses by Abbas et al. in 1967 and Abbas and Latham in 1969, in which the shapes of levitated drops are approximated as spheroids, are corrected. The new analysis shows that drops charged to less than the Rayleigh limit lose shape stability at turning points, with respect to external field strength, and that the instability seen in experiments of Doyle et al. in 1964 and others is not a bifurcation to a family of two-lobed shapes, but rather is a related imperfect bifurcation.

Description: Axisymmetric equilibrium shapes and stability of isolated charged drops are found by solving simultaneously the Young-Laplace equation for surface shape and the Laplace equation for the electric field. Families of two-, three-, and four-lobed shapes that branch from the trunk family of spheres are treated systematically by means of the Galerkin/finite element method and a tessellation that deforms with the free surface. The results show that at the limit found by Rayleigh in 1882 the spherical family exchanges stability with a family of two-lobed shapes, a transcritically bifurcating family, one arm of which proves to consist of stable shapes. The results are reinforced by those of approximating the stable drop shapes as oblate spheroids. Thus oblate drops carrying charge in excess of the Rayleigh limit ought to be seen in experiments, though none have yet been reported.

Description: Numerical experiments were performed to clarify apparent differences between experimental observations and a theoretical prediction of the secondary instability in plane Poiseuille flow. It is shown that subharmonic breakdown is unlikely in natural transition as a result of the initial growth of what we call the 'minus' modes and consequent forcing of Orr-Sommerfeld modes present in the background noise. Subharmonic breakdown was achieved only when these minus modes were continuously suppressed.

Description: Lumley's proper orthogonal decomposition technique is applied to the turbulent flow in a channel. Coherent structures are extracted by decomposing the velocity field into characteristic eddies with random coefficients. A generalization of the shot-noise expansion is used to determine the characteristic eddies in homogeneous spatial directions. Three different techniques are used to determine the phases of the Fourier coefficients in the expansion: (1) one based on the bispectrum, (2) a spatial compactness requirement, and (3) a functional continuity argument. Similar results are found from each of these techniques.

Description: Reattached turbulent boundary layer relaxation downstream of a wall fence is investigated. An adverse pressure gradient is imposed upon it which is adjusted to bring the boundary layer into equilibrium. The pressure gradient is adjusted so as to bring the Clauser parameter G down to a value of about 11.4 and then maintain it constant. In the region from the reattachment point to 2 or 3 reattachment lengths downstream, the boundary layer recovers from the initial major effects of reattachment. Farther downstream, where G is constant, the pressure-gradient parameter changes very slowly and profiles of non-dimensionalized eddy viscosity appear self-similar. However, pressure gradient and eddy viscosity are both roughly twice as large as expected on the basis of previous equilibrium turbulent boundary layer studies.

Description: The fractal and marginal fractal regimes in surface geometry are studied. The basic notions of fractal geometry are applied to a small surface patch in a developed sea, corresponding to the equilibrium wave number spectrum. Topothesy, outer and inner boundaries of the fractal range, and a cascade pattern in surface geometry are discussed. Theoretical predictions of whitecap and foam coverage are presented. A fractal decomposition for a surface patch is developed based on the Karhunen-Loeve expansion. The resulting series formalizes the cascade process of constructing realization of a Gaussian random patch. The implications of the research for microwave remote sensing signatures are considered.

Description: Various secondary and tertiary instabilities in plane channael flow are explored via time-dependent numerical simulations using the incompressible Navier-Stokes equations. Comparisons are made between transitional flows at Reynolds numbers 1500, 5000, and 8000. The lambda vortex, detached shear layer, and inverted vortex regions are identified, and the origin of the latter is explained. The laminar breakdown of the Re = 1500 flow is computed with high resolution, and the nature of its ensuing hairpin eddies is clarified by numerical particle paths. The potential of center-mode, rather than wall-mode transitions is proposed, and the resulting flow structure is described.

Description: Pressure fluctuations in a turbulent channel flow are investigated by analyzing a database obtained from a direct numerical simulation. Detailed statistics associated with the pressure fluctuations are presented. Characteristics associated with the rapid (linear) and slow (nonlinear) pressure are discussed. It is found that the slow pressure fluctuations are larger than the rapid pressure fluctuations throughout the channel except very near the wall, where they are about the same magnitude. This is contrary to the common belief that the nonlinear source terms are negligible compared to the linear source terms. Probability density distributions, power spectra, and two-point correlations are examined to reveal the characteristics of the pressure fluctuations. The global dependence of the pressure fluctuations and pressure-strain correlations are also examined by evaluating the integral associated with Green's function representations of them. In the wall region where the pressure-strain terms are large, most contributions to the pressure-strain terms are from the wall region (i.e., local), whereas away from the wall where the pressure-strain terms are small, contributions are global. Structures of instantaneous pressure and pressure gradients at the wall and the corresponding vorticity field are examined.

Description: A boundary layer created on an infinite flat plate by a time-dependent freestream velocity vector whose magnitude is independent of time but whose direction changes at a constant angular velocity is theoretically studied using Reynolds-number scaling laws and numerical simulations performed over a range of Reynolds numbers. Results obtained with a higher-order version of existing theories of the Ekman layer are shown to agree well with the numerical results at three Reynolds numbers. The present results can be extrapolated to the case of high Reynolds numbers. The Reynolds-averaged equations reduce to a one-dimensional steady problem, making possible the easy and accurate testing of turbulence models.

Description: Turbulence-producing events in turbulent channel flow were found to be predominantly associated with a symmetric vortical structures rather than pairs of counter-rotating structures. An asymmetry-preserving averaging scheme was devised, allowing a picture of the 'average' structure that more closely resembles the instantaneous one to be obtained. In addition, these structures were found to persist for long distances with little change while convecting downstream.

Description: Plane stagnation-point flow is modulated in the free stream so that the velocity components are proportional to K(H) + K cos omega t. Similarity solutions of the Navier-Stokes equations are examined using high-frequency asymptotics for K and K(H) of unit order. Special attention is focused on the steady streaming generated in this flow with strongly non-parallel streamlines.

Description: This paper presents experimental data taken in the forward region of a separated internal free-shear layer produced in an internal cavity flowfield. It has been found that in the region very near the forward restrictor, experimental velocity profiles agree closely with the exact Stuart instability velocity profile with various values of a steepness parameter. Reynolds shear-stress profiles suggest the presence of counter-rotating longitudinal vortices. Spectral analysis by the maximum entropy method of the time samples within the vortices indicates subharmonic and harmonic components of the fundamental frequency, with a weak indication of the fundamental frequency itself.

Description: The effect of enhancement devices on flow boiling heat transfer in coolant channels, which are heated either from the top side or uniformly was studied. Studies are completed of the variations in the local (axial and circumferential) and mean heat transfer coefficients in horizontal, top-heated coolant channels with smooth walls and internal heat transfer enhancement devices. The working fluid is freon-11. The objectives are to: (1) examine the variations in both the mean and local (axial and circumferential) heat transfer coefficients for a circular coolant channel with either smooth walls or with both a twisted tape and spiral finned walls; (2) examine the effect of channel diameter (and the length-to-diameter aspect ratio) variations for the smooth wall channel; and (3) develop and improved data reduction analysis. The case of the top-heated, horizontal flow channel with smooth wall (1.37 cm inside diameter, and 122 cm heated length) was completed. The data were reduced using a preliminary analysis based on the heated hydraulic diameter. Preliminary examination of the local heat transfer coefficient variations indicated that there are significant axial and circumferential variations. However, it appears that the circumferential variation is more significant than the axial ones. In some cases, the circumferential variations were as much as a factor of ten. The axial variations rarely exceeded a factor of three.

Description: The immediate objective of this research is to measure liquid film thickness from the two equilibrium phases of a monotectic system in order to estimate the film pressure of each phase. Thus liquid film thicknesses on the inside walls of the prism cell above the liquid level have been measured elliposmetrically for the monotectic system of succinonitrile and water. The thickness varies with temperature and composition of each plane. The preliminary results from both layers at 60 deg angle of incidence show nearly uniform thickness from about 21 to 23 C. The thickness increases with temperature but near 30 C the film appears foggy and scatters the laser beam. As the temperature of the cell is raised beyond room temperature it becomes increasingly difficult to equalize the temperature inside and outside the cell. The fogging may also be an indication that solution, not pure water, is adsorbed onto the substrate. Nevertheless, preliminary results suggest that ellipsometric measurement is feasible and necessary to measure more accurately and rapidly the film thickness and to improve thermal control of the prism walls.

Description: Flush-mounted hot-film gages have proved effective in detecting boundary-layer transition and in measuring skin friction but with limited success in detecting laminar separation and reattachment. The development of multielement micro hot-film sensors, and the recent discovery of the phase reversal phenomena associated with low-frequency dynamic shear stress signals across regions of laminar separation and turbulent reattachment, have made it possible to simultaneously and unambiguously detect these surface shear layer characteristics. Experiments were conducted on different airfoils at speeds ranging from low subsonic to transonic speeds to establish the technique for incompressible and compressible flow applications. The multielement dynamic shear stress sensor technique was successfully used to detect laminar separation, turbulent reattachment, as well as, shock induced laminar and turbulent separation.