Search

von Neumann modelled the evolution of two-dimensional soap froths as a purely diffusive phenomenon; the area growth of a given cell was found to depend only on the geometry of the bubble lattice. In the model, hexagons are stable, pentagons shrink and heptagons grow. The simplest equivalent to the area growth law is / approximately t(sub beta). The result depends on assuming (1) an incompressible gas; (2) bubble walls which meet at 120 deg and (3) constant wall thickness and curvature. Each assumption is borne out in experiments except the last one: bubble wall thickness between connecting cells varies in unit gravity because of gravity drainage. The bottom part of the soap membrane is thickened, the top part is thinned, such that gas diffusion across the membrane shows a complex dependence on gravity. As a result, experimental tests of von Neumann's law have been influenced by effects of gravity; fluid behavior along cell borders can give non-uniform wall thicknesses and thus alter the effective area and gas diffusion rates between adjacent bubbles. For area plotted as a function of time, Glazier (J.A. Glazier, S.P. Gross, and I. Stavans, Phys. Rev. A. 36, 306 (1987); J. Stavans, J.A, Glazier, Phys. Rev. Lett. 62, 1318 (1989).) suggest that in some cases their failure to observe von Neumann's predicted growth exponent ((sup beta)theor(sup =1; beta)exp(sup =0.70 + 0.10)) may have been the result of such "fluid drainage onto the lower glass plate". Additional experiments which varied plate spacing gave different beta exponents in a fashion consistent with this suggestion. During preliminary long duration experiments (approximately 100 h) aboard Spacelab-J, a low-gravity test of froth coarsening has examined (1) power law scaling of von Neumann's law (beta values) in the appropriate diffusive limits; (2) new bubble lattice dynamics such as greater fluid wetting behavior on froth membranes in low gravity; and (3) explicit relations for the gravity dependence of the second moment (or disorder parameter) governing the geometric spread in cell-sidedness around the mean of perfect hexagonal filling. By reducing the gravity-induced distortion in lattice wall thickness, the diffusion-limited regime of bubble coarsening becomes available for performing critical tests of network dynamics.

The effects of gravity in influencing the theoretical limit for bubble lattice coarsening and aging behavior, otherwise called von Neumann's law, is examined theoretically and experimentally. Preliminary microgravity results will be discussed.

A two component Point Doppler Velocimetry (PDV) system has been developed and tested. Improvements were made to an earlier PDV system, in terms of experimental techniques, as well as the data acquisition and reduction software. Measurements of the streamwise and spanwise mean and fluctuating velocities for flows from a rectangular channel and over an NACA 0012 airfoil were made, and the data were compared against hot wire data. The closest to the airfoil surface that PDV measurements could be made was on the order of 0.005 m(0.2", z/c = 0.0169). When the PDV and hot wire data were compared, the time traces for each appeared similar. The mean velocities agreed to within plus or minus 2 m/sec, while the RMS velocities agreed to plus or minus 0.4 m/sec. While the PDV time autocorrelations agreed with those of the hot wire data, the PDV power spectral densities were noisier above 750 Hz. A major source of error in these experiments was determined to be the drifting of the iodine cell stem temperatures. While the stem temperatures were controlled to within plus or minus 0.1 C, this could lead to a frequency shift of as much as 6 MHz, which translates into an error of 1.6 m/sec for the back scatter channel, and up to 6.9 m/sec for the forward scatter channel. These error estimates are consistent with the observed error magnitudes.

For a typical three dimensional flow in a practical engineering device, the time spent in grid generation can take 70 percent of the total analysis effort, resulting in a serious bottleneck in the design/analysis cycle. The present research attempts to develop a procedure that can considerably reduce the grid generation effort. The DRAGON grid, as a hybrid grid, is created by means of a Direct Replacement of Arbitrary Grid Overlapping by Nonstructured grid. The DRAGON grid scheme is an adaptation to the Chimera thinking. The Chimera grid is a composite structured grid, composing a set of overlapped structured grids, which are independently generated and body-fitted. The grid is of high quality and amenable for efficient solution schemes. However, the interpolation used in the overlapped region between grids introduces error, especially when a sharp-gradient region is encountered. The DRAGON grid scheme is capable of completely eliminating the interpolation and preserving the conservation property. It maximizes the advantages of the Chimera scheme and adapts the strengths of the unstructured and while at the same time keeping its weaknesses minimal. In the present paper, we describe the progress towards extending the DRAGON grid technology into three dimensions. Essential and programming aspects of the extension, and new challenges for the three-dimensional cases, are addressed.

This article describes the application of the Multidisciplinary Analysis (MDA) solver, Spectrum, in analyzing a hydrogen-cooled hypersonic cowl leading-edge structure. Spectrum, a multiphysics simulation code based on the finite element method, addresses compressible and incompressible fluid flow, structural, and thermal modeling, as well as the interactions between these disciplines. Fluid-solid-thermal interactions in a hydrogen impingement-cooled leading edge are predicted using Spectrum. Two- and semi-three-dimensional models are considered for a leading edge impingement coolant, concept under either specified external heat flux or aerothermodynamic heating from a Mach 5 external flow interaction. The solution accuracy is demonstrated from mesh refinement analysis. With active cooling, the leading edge surface temperature is drastically reduced from 1807 K of the adiabatic condition to 418 K. The internal coolant temperature profile exhibits a sharp gradient near channel/solid interface. Results from two different cooling channel configurations are also presented to illustrate the different behavior of alternative active cooling schemes.

In order to expand the predictive capability of single-point turbulence closure models to account for the early-stage transition regime, a methodology for the formulation and calibration of model equations for the ensemble-averaged disturbance kinetic energy and energy dissipation rate is presented. First the decay of laminar disturbances and turbulence in mean shear-free flows is studied. In laminar flows, such disturbances are linear superpositions of modes governed by the Orr-Sommerfeld equation. In turbulent flows, disturbances are described through transport equations for representative mean quantities. The link between a description based on a deterministic evolution equation and a probability based mean transport equation is established. Because an uncertainty in initial conditions exists in the laminar as well as the turbulent regime, a probability distribution must be defined even in the laminar case. Using this probability distribution, it is shown that the exponential decay of the linear modes in the laminar regime can be related to a power law decay of both the (ensemble) mean disturbance kinetic energy and the dissipation rate. The evolution of these mean disturbance quantities is then described by transport equations similar to those for the corresponding turbulent decaying flow. Second, homogeneous shear flow, where disturbances can be described by rapid distortion theory (RDT), is studied. The relationship between RDT and linear stability theory is exploited in order to obtain a closed set of modeled equations. The linear disturbance equations are solved directly so that the numerical simulation yields a database from which the closure coefficients in the ensemble-averaged disturbance equations can be determined.

Since its inception 40 years ago, Marshall Space Flight Center (MSFC) has had the need to maintain and advance state-of-the-art flow analysis and cold-flow testing capability to support its roles and missions. This overview discusses the recent organizational changes that have occurred at MSFC with emphasis on the resulting three groups that form the core of fluid dynamics expertise at MSFC: the Fluid Physics and Dynamics Group, the Applied Fluid Dynamics Analysis Group, and the Experimental Fluid Dynamics Group. Recently completed activities discussed include the analysis and flow testing in support of the Fastrac engine design, the X-33 vehicle design, and the X34 propulsion system design. Ongoing activities include support of the RLV vehicle design, Liquid Fly Back Booster aerodynamic configuration definition, and RLV focused technologies development. Other ongoing activities discussed are efforts sponsored by the Center Director's Discretionary Fund (CDDF) to develop an advanced incompressible flow code and to develop optimization techniques. Recently initiated programs and their anticipated required fluid dynamics support are discussed. Based on recent experiences and on the anticipated program needs, required analytical and experimental technique improvements are presented. Due to anticipated budgetary constraints, there is a strong need to leverage activities and to pursue teaming arrangements in order to advance the state-of-the-art and to adequately support concept development. Throughout this overview there is discussion of the lessons learned and of the capabilities demonstrated and established in support of the hardware development programs.

An evaporative cooling membrane device is disclosed having a flat or pleated plate housing with an enclosed bottom and an exposed top that is covered with at least one sheet of hydrophobic porous material having a thin thickness so as to serve as a membrane. The hydrophobic porous material has pores with predetermined dimensions so as to resist any fluid in its liquid state from passing therethrough but to allow passage of the fluid in its vapor state, thereby, causing the evaporation of the fluid and the cooling of the remaining fluid. The fluid has a predetermined flow rate. The evaporative cooling membrane device has a channel which is sized in cooperation with the predetermined flow rate of the fluid so as to produce laminar flow therein. The evaporative cooling membrane device provides for the convenient control of the evaporation rates of the circulating fluid by adjusting the flow rates of the laminar flowing fluid.

Three problems concerning laminar-turbulent transition are addressed by solving a series of initial value problems. The first problem is the calculation of resonance within the continuous spectrum of the Blasius boundary layer. The second is calculation of the growth of Tollmien-Schlichting waves that are a direct result of disturbances that only lie outside of the boundary layer. And, the third problem is the calculation of non-parallel effects. Together, these problems represent a unified approach to the study of freestream disturbance effects that could lead to transition. Solutions to the temporal, initial-value problem with an inhomogeneous forcing term imposed upon the flow is sought. By solving a series of problems, it is shown that: A transient disturbance lying completely outside of the boundary layer can lead to the growth of an unstable Tollmien-Schlichting wave. A resonance with the continuous spectrum leads to strong amplification that may provide a mechanism for bypass transition once nonlinear effects are considered. A disturbance with a very weak unstable Tollmien-Schlichting wave can lead to a much stronger Tollmien-Schlichting wave downstream, if the original disturbance has a significant portion of its energy in the continuum modes.

A new transport equation for intermittency factor is proposed to model transitional flows. The intermittent behavior of the transitional flows is incorporated into the computations by modifying the eddy viscosity, mu(sub t), obtainable from a turbulence model, with the intermittency factor, gamma: mu(sub t, sup *) = gamma.mu(sub t). In this paper, Menter's SST model (Menter, 1994) is employed to compute mu(sub t) and other turbulent quantities. The proposed intermittency transport equation can be considered as a blending of two models - Steelant and Dick (1996) and Cho and Chung (1992). The former was proposed for near-wall flows and was designed to reproduce the streamwise variation of the intermittency factor in the transition zone following Dhawan and Narasimha correlation (Dhawan and Narasimha, 1958) and the latter was proposed for free shear flows and was used to provide a realistic cross-stream variation of the intermittency profile. The new model was used to predict the T3 series experiments assembled by Savill (1993a, 1993b) including flows with different freestream turbulence intensities and two pressure-gradient cases. For all test cases good agreements between the computed results and the experimental data are observed.

The temporal evolution of combustion flowfields established by the interaction between wedge-shaped bodies and explosive hydrogen-oxygen-nitrogen mixtures accelerated to hypersonic speeds in an expansion tube is investigated. The analysis is carried out using a fully implicit, time-accurate, computational fluid dynamics code that we developed recently for solving the Navier-Stokes equations for a chemically reacting gas mixture. The numerical results are compared with experimental data from the Stanford University expansion tube for two different gas mixtures at Mach numbers of 4.2 and 5.2. The experimental work showed that flow unstart occurred for the Mach 4.2 cases. These results are reproduced by our numerical simulations and, more significantly, the causes for unstart are explained. For the Mach 5.2 mixtures, the experiments and numerical simulations both produced stable combustion. However, the computations indicate that in one case the experimental data were obtained during the transient phase of the flow; that is, before steady state had been attained.

The Geophysical Fluid Flow Cell (GFFC) experiment performed visualizations of thermal convection in a rotating differentially heated spherical shell of fluid. In these experiments dielectric polarization forces are used to generate a radially directed buoyancy force. This enables the laboratory simulation of a number of geophysically and astrophysically important situations in which sphericity and rotation both impose strong constraints on global scale fluid motions. During USML-2 a large set of experiments with spherically symmetric heating were carried out. These enabled the determination of critical points for the transition to various forms of nonaxisymmetric convection and, for highly turbulent flows, the transition latitudes separating the different modes of motion. This paper presents a first analysis of these experiments as well as data on the general performance of the instrument during the USML-2 flight.

The objective of the research is to perform Direct Numerical Simulations (DNS) and Large Eddy Simulations (LES) for film cooling process, and to evaluate and improve advanced forms of the two equation turbulence models for turbine blade surface flow analysis. The DNS/LES were used to resolve the large eddies within the flow field near the coolant jet location. The work involved code development and applications of the codes developed to the film cooling problems. Five different codes were developed and utilized to perform this research. This report presented a summary of the development of the codes and their applications to analyze the turbulence properties at locations near coolant injection holes.

High performance piezoelectric polymers are of interest to NASA as they may be useful for a variety of sensor applications. Over the past few years research on piezoelectric polymers has led to the development of promising high temperature piezoelectric responses in some novel polyimides. In this study, a series of polyimides have been studied with systematic variations in the diamine monomers that comprise the polyimide while holding the dianhydride constant. The effect of structural changes, including variations in the nature and concentration of dipolar groups, on the remanent polarization and piezoelectric coefficient is examined. Fundamental structure-piezoelectric property insight will enable the molecular design of polymers possessing distinct improvements over state-of-the-art piezoelectric polymers including enhanced polarization, polarization stability at elevated temperatures, and improved processability.

The technical details are summarized below: Compressible and incompressible versions of a three-dimensional unstructured mesh Reynolds-averaged Navier-Stokes flow solver have been differentiated and resulting derivatives have been verified by comparisons with finite differences and a complex-variable approach. In this implementation, the turbulence model is fully coupled with the flow equations in order to achieve this consistency. The accuracy demonstrated in the current work represents the first time that such an approach has been successfully implemented. The accuracy of a number of simplifying approximations to the linearizations of the residual have been examined. A first-order approximation to the dependent variables in both the adjoint and design equations has been investigated. The effects of a "frozen" eddy viscosity and the ramifications of neglecting some mesh sensitivity terms were also examined. It has been found that none of the approximations yielded derivatives of acceptable accuracy and were often of incorrect sign. However, numerical experiments indicate that an incomplete convergence of the adjoint system often yield sufficiently accurate derivatives, thereby significantly lowering the time required for computing sensitivity information. The convergence rate of the adjoint solver relative to the flow solver has been examined. Inviscid adjoint solutions typically require one to four times the cost of a flow solution, while for turbulent adjoint computations, this ratio can reach as high as eight to ten. Numerical experiments have shown that the adjoint solver can stall before converging the solution to machine accuracy, particularly for viscous cases. A possible remedy for this phenomenon would be to include the complete higher-order linearization in the preconditioning step, or to employ a simple form of mesh sequencing to obtain better approximations to the solution through the use of coarser meshes. . An efficient surface parameterization based on a free-form deformation technique has been utilized and the resulting codes have been integrated with an optimization package. Lastly, sample optimizations have been shown for inviscid and turbulent flow over an ONERA M6 wing. Drag reductions have been demonstrated by reducing shock strengths across the span of the wing.

Detailed heat transfer measurements and predictions are given for a turbine rotor with 136 deg of turning and an axial chord of 12.7 cm. Data were obtained for inlet Reynolds numbers of 0.5 and 1.0 x 10(exp 6), for isentropic exit Mach numbers of 1.0 and 1.3, and for inlet turbulence intensities of 0.25% and 7.0%. Measurements were made in a linear cascade having a highly three-dimensional flow field resulting from thick inlet boundary layers. The purpose of the work is to provide benchmark quality data for three-dimensional CFD code and model verification. Data were obtained by a steady-state technique using a heated, isothermal blade. Heat fluxes were determined from a calibrated resistance layer in conjunction with a surface temperature measured by calibrated liquid crystals. The results show the effects of strong secondary vortical flows, laminar-to-turbulent transition, shock impingement, and increased inlet turbulence on the surface heat transfer.

The objectives of this research are: (1) to develop and implement a new methodology for large eddy simulation of (LES) of high-speed reacting turbulent flows. (2) To develop algebraic turbulence closures for statistical description of chemically reacting turbulent flows.

Complex streamwise vorticity fields are observed in the evolution of non-circular jets. Generation mechanisms are investigated via Reynolds-averaged (RANS), large-eddy (LES) and direct numerical (DNS) simulations of laminar and turbulent rectangular jets. Complex vortex interactions are found in DNS of laminar jets, but axis-switching is observed only when a single instability mode is present in the incoming mixing layer. With several modes present, the structures are not coherent and no axis-switching occurs, RANS computations also produce no axis-switching. On the other hand, LES of high Reynolds number turbulent jets produce axis-switching even for cases with several instability modes in the mixing layer. Analysis of the source terms of the mean streamwise vorticity equation through post-processing of the instantaneous results shows that, complex interactions of gradients of the normal and shear Reynolds stresses are responsible for the generation of streamwise vorticity which leads to axis-switching. RANS computations confirm these results. k - epsilon turbulence model computations fail to reproduce the phenomenon, whereas algebraic Reynolds stress model (ASM) computations, in which the secondary normal and shear stresses are computed explicitly, succeeded in reproducing the phenomenon accurately.

An unsteady, multiblock, Reynolds Averaged Navier Stokes solver based on Runge-Kutta scheme and Pseudo-time step for turbo-machinery applications was developed. The code was validated and assessed against analytical and experimental data. It was used to study a variety of physical mechanisms of unsteady, three-dimensional, turbulent, transitional, and cooling flows in compressors and turbines. Flow over a cylinder has been used to study effects of numerical aspects on accuracy of prediction of wake decay and transition, and to modify K-epsilon models. The following simulations have been performed: (a) Unsteady flow in a compressor cascade: Three low Reynolds number turbulence models have been assessed and data compared with Euler/boundary layer predictions. Major flow features associated with wake induced transition were predicted and studied; (b) Nozzle wake-rotor interaction in a turbine: Results compared to LDV data in design and off-design conditions, and cause and effect of unsteady flow in turbine rotors were analyzed; (c) Flow in the low-pressure turbine: Assessed capability of the code to predict transitional, attached and separated flows at a wide range of low Reynolds numbers and inlet freestream turbulence intensity. Several turbulence and transition models have been employed and comparisons made to experiments; (d) leading edge film cooling at compound angle: Comparisons were made with experiments, and the flow physics of the associated vortical structures were studied; and (e) Tip leakage flow in a turbine. The physics of the secondary flow in a rotor was studied and sources of loss identified.

The asymptotic solutions, described by Tennekes and Lumley (1972), for surface flows in a channel, pipe or boundary layer at large Reynolds numbers are revisited. These solutions can be extended to more complex flows such as the flows with various pressure gradients, zero wall stress and rough surfaces, etc. In computational fluid dynamics (CFD), these solutions can be used as the boundary conditions to bridge the near-wall region of turbulent flows so that there is no need to have the fine grids near the wall unless the near-wall flow structures are required to resolve. These solutions are referred to as the wall functions. Furthermore, a generalized and unified law of the wall which is valid for whole surface layer (including viscous sublayer, buffer layer and inertial sublayer) is analytically constructed. The generalized law of the wall shows that the effect of both adverse and favorable pressure gradients on the surface flow is very significant. Such as unified wall function will be useful not only in deriving analytic expressions for surface flow properties but also bringing a great convenience for CFD methods to place accurate boundary conditions at any location away from the wall. The extended wall functions introduced in this paper can be used for complex flows with acceleration, deceleration, separation, recirculation and rough surfaces.

Two new calibration algorithms were developed for the calibration of non-nulling multi-hole probes in compressible, subsonic flowfields. The reduction algorithms are robust and able to reduce data from any multi-hole probe inserted into any subsonic flowfield to generate very accurate predictions of the velocity vector, flow direction, total pressure and static pressure. One of the algorithms PROBENET is based on the theory of neural networks, while the other is of a more conventional nature (polynomial approximation technique) and introduces a novel idea of local least-squares fits. Both algorithms have been developed to complete, user-friendly software packages. New technology was developed for the fabrication of miniature multi-hole probes, with probe tip diameters all the way down to 0.035". Several miniature 5- and 7-hole probes, with different probe tip geometries (hemispherical, conical, faceted) and different overall shapes (straight, cobra, elbow probes) were fabricated, calibrated and tested. Emphasis was placed on the development of four stainless-steel conical 7-hole probes, 1/16" in diameter calibrated at NASA Langley for the entire subsonic regime. The developed calibration algorithms were extensively tested with these probes demonstrating excellent prediction capabilities. The probes were used in the "trap wing" wind tunnel tests in the 14'x22' wind tunnel at NASA Langley, providing valuable information on the flowfield over the wing. This report is organized in the following fashion. It consists of a "Technical Achievements" section that summarizes the major achievements, followed by an assembly of journal articles that were produced from this project and ends with two manuals for the two probe calibration algorithms developed.

An experimental examination of crossflow instability and transition on a 45deg swept wing was conducted in the Arizona State University Unsteady Wind Tunnel. The stationary-vortex pattern and transition location are visualized by using both sublimating chemical and liquid-crystal coatings. Extensive hot-wire measurements were obtained at several measurement stations across a single vortex track. The mean and travelling wave disturbances were measured simultaneously. Stationary crossflow disturbance profiles were determined by subtracting either a reference or a span-averaged velocity profile from the mean velocity data. Mean, stationary crossflow, and traveling wave velocity data were presented as local boundary layer profiles and contour plots across a single stationary crossflow vortex track. Disturbance mode profiles and growth rates were determined. The experimental data are compared with predictions from linear stability theory.

Computational Fluid Dynamics (CFD) numerical simulations of low-density shock-wave interactions for an incident shock impinging on a cylinder have been performed. Flow-field density gradient and surface pressure and heating define the type of interference pattern and corresponding perturbations. The maximum pressure and heat transfer level and location for various interaction types (i.e., shock-wave incidence with respect to the cylinder) are presented. A time-accurate solution of the Type IV interference is employed to demonstrate the establishment and the steadiness of the low-density flow interaction.

The Near Infrared Camera and Multi Object Spectrometer (NICMOS) was installed in the Hubble Space Telescope (MST) in February 1997. Shortly thereafter, the instrument experienced a thermal short in its solid nitrogen dewar system which will significantly shorten the instrument's useful life. A reverse Brayton cycle mechanical refrigerator will be installed during the Third Servicing Mission (SM3) to provide cooling for the instrument, and thereby extend its operations. A Capillary Pump Loop (CPL) and radiator system was designed, built and tested to remove up to 500 watts of heat from the mechanical cryocooler and its associated electronics. The HST Orbital Systems Test (HOST) platform was flown on the Space Shuttle Discovery (STS-95) as a flight demonstration of the cryocooler system, CPL control electronics, and the CPL/Radiator. This paper will present the flight test results and thermal performance of the CPL system in detail.

This report summarizes the research undertaken, at Aeronautics Department of the Massachusetts Institute of Technology, during the approximately five year period, February 94 - March 99. This work is part of a larger effort aimed at providing a reliable fast turn around capability for the prediction of hypersonic flows over complete vehicle configurations.

The present work reports a study of bubble generation under reduced gravity conditions for both co-flow and cross-flow configurations. Experiments were performed aboard the DC-9 Reduced Gravity Aircraft at NASA Glenn Research Center, using an air-water system. Three different flow tube diameters were used: 1.27, 1.9, and 2.54 cm. Two different ratios of air injection nozzle to tube diameters were considered: 0.1 and 0.2. Gas and liquid volumetric flow rates were varied from 10 to 200 ml/s. It was experimentally observed that with increasing superficial liquid velocity, the bubbles generated decreased in size. The bubble diameter was shown to increase with increasing air injection nozzle diameters. As the tube diameter was increased, the size of the detached bubbles increased. Likewise, as the superficial liquid velocity was increased, the frequency of bubble formation increased and thus the time to detach forming bubbles decreased. Independent of the flow configuration (for either single nozzle or multiple nozzle gas injection), void fraction and hence flow regime transition can be controlled in a somewhat precise manner by solely varying the gas and liquid volumetric flow rates. On the other hand, it is observed that uniformity of bubble size can be controlled more accurately by using single nozzle gas injection than by using multiple port injection, since this latter system gives rise to unpredictable coalescence of adjacent bubbles. A theoretical model, based on an overall force balance, is employed to study single bubble generation in the dynamic and bubbly flow regime. Under conditions of reduced gravity, the gas momentum flux enhances bubble detachment; however, the surface tension forces at the nozzle tip inhibits bubble detachment. Liquid drag and inertia can act either as attaching or detaching force, depending on the relative velocity of the bubble with respect to the surrounding liquid. Predictions of the theoretical model compare well with performed experiments. However, at higher superficial,liquid velocities, the bubble neck length begins to significantly deviate from the value of the air injection nozzle diameter and thus the theory no longer predicts the experiment behavior. Effects of fluid properties, injection geometry and flow conditions on generated bubble size are investigated using the theoretical model. It is shown that bubble diameter is larger in a reduced gravity environment than in a normal gravity environment at similar flow condition and flow geometry.

This paper summarizes experimental and computational results on the mixing of opposed rows of jets with a confined subsonic crossflow in rectangular ducts. The studies from which these results were excerpted investigated flow and geometric variations typical of the complex three-dimensional flowfield in the combustion chambers in gas turbine engines. The principal observation was that the momentum-flux ratio, J, and the orifice spacing, S/H, were the most significant flow and geometric variables. Jet penetration was critical, and penetration decreased as either momentum-flux ratio or orifice spacing decreased. It also appeared that jet penetration remained similar with variations in orifice size, shape, spacing, and momentum-flux ratio when the orifice spacing was inversely proportional to the square-root of the momentum-flux ratio. It was also seen that planar averages must be considered in context with the distributions. Note also that the mass-flow ratios and the orifices investigated were often very large (jet-to-mainstream mass-flow ratio > 1 and the ratio of orifices-area-to-mainstream- cross-sectional-area up to 0.5, respectively), and the axial planes of interest were often just downstream of the orifice trailing edge. Three-dimensional flow was a key part of efficient mixing and was observed for all configurations.

A multi-block, three-dimensional Navier-Stokes code has been used to compute heat transfer coefficient on the blade, hub and shroud for a rotating high-pressure turbine blade with 172 film-cooling holes in eight rows. Film cooling effectiveness is also computed on the adiabatic blade. Wilcox's k-omega model is used for modeling the turbulence. Of the eight rows of holes, three are staggered on the shower-head with compound-angled holes. With so many holes on the blade it was somewhat of a challenge to get a good quality grid on and around the blade and in the tip clearance region. The final multi-block grid consists of 4784 elementary blocks which were merged into 276 super blocks. The viscous grid has over 2.2 million cells. Each hole exit, in its true oval shape, has 80 cells within it so that coolant velocity, temperature, k and omega distributions can be specified at these hole exits. It is found that for the given parameters, heat transfer coefficient on the cooled, isothermal blade is highest in the leading edge region and in the tip region. Also, the effectiveness over the cooled, adiabatic blade is the lowest in these regions. Results for an uncooled blade are also shown, providing a direct comparison with those for the cooled blade. Also, the heat transfer coefficient is much higher on the shroud as compared to that on the hub for both the cooled and the uncooled cases.

A computational study has been performed to predict the distribution of convective heat transfer coefficient on a simulated blade tip with cooling holes. The purpose of the examination was to assess the ability of a three-dimensional Reynolds-averaged Navier-Stokes solver to predict the rate of tip heat transfer and the distribution of cooling effectiveness. To this end, the simulation of tip clearance flow with blowing of Kim and Metzger was used. The agreement of the computed effectiveness with the data was quite good. The agreement with the heat transfer coefficient was not as good but improved away from the cooling holes. Numerical flow visualization showed that the uniformity of wetting of the surface by the film cooling jet is helped by the reverse flow due to edge separation of the main flow.

An experiment was conducted to study the effect of velocity and sweep angle on the critical distance in ice accretion formation on swept wings at glaze ice conditions. The critical distance is defined as the distance from the attachment line to the beginning of the zone where roughness elements develop into glaze ice feathers. Icing runs were performed on a NACA 00 1 2 swept wing tip at velocities of 75, 100, 150, and 200 miles per hour. At each velocity and tunnel condition, the sweep angle was changed from 0 deg to 45 deg at 5 deg increments. Casting data, ice shape tracings, and close-up photographic data were obtained. The results showed that at given velocity and tunnel conditions, as the sweep angle is increased from 0 deg to 25 deg the critical distance slowly decreases. As the sweep angle is increased past 25 deg, the critical distance starts decreasing more rapidly. For 75 and 100 mph it reaches a value of 0 millimeters at 35 deg. For 150 and 200 mph it reaches a value of 0 millimeters at 40 deg. On the ice accretion, as the sweep angle is increased from 0 deg to 25 deg, the extent of the attachment line zone slowly decreases. In the glaze ice feathers zone, the angle that the preferred direction of growth of the feathers makes with respect to the attachment line direction increases. But overall, the ice accretions remain similar to the 0 deg sweep angle case. As the sweep angle is increased above 25 deg, the extent of the attachment line zone decreases rapidly and complete scallops form at 35 deg sweep angle for 75 and 100 mph, and at 40 deg for 150 and 200 mph.

This paper presents software for solving the non-conforming fluid structure interfaces in aeroelastic simulation. It reviews the algorithm of interpolation and integration, highlights the flexibility and the user-friendly feature that allows the user to select the existing structure and fluid package, like NASTRAN and CLF3D, to perform the simulation. The presented software is validated by computing the High Speed Civil Transport model.

Experimental and theoretical studies of low speed leading edge boundary layer receptivity to free-stream vorticity produced by upstream wires normal to the leading edge are discussed. Data include parametric variations in leading edge configuration and details of the incident disturbance field including single and multiple wakes. The induced disturbance amplitude increases with increases in the leading edge diameter and wake interactions. Measurements agree with the theory of M. E. Goldstein.

Reduction of a large system of equations to a lower-dimensional system of similar dynamics is investigated. For dynamical systems with disparate timescales, a criterion for determining redundant dimensions and a general reduction method based on the minimization of evolution rate are proposed.

The secondary instability of three-dimensional incompressible boundary layers is studied using Floquet theory. Starting from the equilibrium solutions that we obtained from the PSE computations documented in Part 1, we investigate the region where a purely stationary crossflow disturbance saturates for its secondary instability characteristics utilizing developed global and local eigenvalue solvers that are based on the Implicitly Restarted Arnoldi Method, and a Newton-Raphson technique, respectively. The main focuses of this study are on the existence of multiple roots in the eigenvalue spectrum that could explain experimental observations of time-dependent occurrences of an explosive growth of traveling disturbances, on the routes by which high-frequency disturbances enter the boundary layer, as well as on gaining more information about threshold amplitudes for the growth of secondary disturbances.

A comparison of the NPARC, PAB, and WIND (previously known as NASTD) Navier-Stokes solvers is made for two flow cases with turbulent mixing as the dominant flow characteristic, a two-dimensional ejector nozzle and a Mach 1.5 elliptic jet. The objective of the work is to determine if comparable predictions of nozzle flows can be obtained from different Navier-Stokes codes employed in a multiple site research program. A single computational grid was constructed for each of the two flows and used for all of the Navier-Stokes solvers. In addition, similar k-e based turbulence models were employed in each code, and boundary conditions were specified as similarly as possible across the codes. Comparisons of mass flow rates, velocity profiles, and turbulence model quantities are made between the computations and experimental data. The computational cost of obtaining converged solutions with each of the codes is also documented. Results indicate that all of the codes provided similar predictions for the two nozzle flows. Agreement of the Navier-Stokes calculations with experimental data was good for the ejector nozzle. However, for the Mach 1.5 elliptic jet, the calculations were unable to accurately capture the development of the three dimensional elliptic mixing layer.

A three-dimensional viscous flow analysis is performed using a time-marching Reynolds-averaged Navier-Stokes code for a 3:1 rectangular nozzle with two delta tabs located at the nozz1e exit plane to enhance mixing. Two flow configurations, a subsonic jet case and a supersonic jet case using the same rate configuration, which were previously studied experimentally, are computed and compared with the experimental data. The experimental data include streamwise velocity and vorticity distributions for the subsonic case, and Mach number distributions for the supersonic case, at various axial locations downstream of the nozzle exit. The computational results show very good agreement with the experimental data. In addition, the effect of compressibility on vorticity dynamics is examined by comparing the vorticity contours of the subsonic jet case with those of the supersonic jet case which were not measured in the experiment.

A combined experimental and computational study has been performed to investigate the detailed distribution of convective heat transfer coefficients on the first stage blade tip surface for a geometry typical of large power generation turbines (>1OOMW). This paper is concerned with the numerical prediction of the tip surface heat transfer. Good comparison with the experimental measured distribution was achieved through accurate modeling of the most important features of the blade passage and heating arrangement as well as the details of experimental rig likely to affect the tip heat transfer. A sharp edge and a radiused edge tip were considered. The results using the radiused edge tip agreed better with the experimental data. This improved agreement was attributed to the absence of edge separation on the tip of the radiused edge blade.

This paper presents a rational approach to modelling the triple velocity correlations that appear in the transport equations for the Reynolds stresses. All existing models of these correlations have largely been formulated on phenomenological grounds and are defective in one important aspect: they all neglect to allow for the dependence of these correlations on the local gradients of mean velocity. The mathematical necessity for this dependence will be demonstrated in the paper. The present contribution lies in the novel use of Group Representation Theory to determine the most general tensorial form of these correlations in terms of all the second- and third-order tensor quantities that appear in the exact equations that govern their evolution. The requisite representation did not exist in the literature and therefore had to be developed specifically for this purpose by Professor G. F. Smith. The outcome of this work is a mathematical framework for the construction of algebraic, explicit, and rational models for the triple velocity correlations that are theoretically consistent and include all the correct dependencies. Previous models are reviewed, and all are shown to be an incomplete subset of this new representation, even to lowest order.

An experimental study to investigate the effect of steady and oscillatory (with zero net mass flux) blowing on cavity resonance is undertaken. The objective is to study the basic mechanisms of the control of cavity resonance. An actuator is designed and calibrated to generate either steady blowing or oscillatory blowing with A zero net mass flux. The results of the experiment show that both steady and oscillatory blowing are effective, and reduce the amplitude of the dominant resonant mode by 1OdB. The oscillatory blowing is however found to be more superior in that the same effectiveness could be accomplished with a momentum coefficient an order of magnitude smaller than for steady blowing. The experiment also confirms the results of previous computations that suggest the forcing frequency for oscillatory blowing must not be at harmonic frequencies of the cavity resonant modes.

The two equation k-epsilon turbulence model of Chien has been implemented in the WIND Navier-Stokes flow solver. Details of the numerical solution algorithm, initialization procedure, and stability enhancements are described. Results obtained with this version of the model are compared with those from the Chien k-epsilon model in the NPARC Navier-Stokes code and from the WIND SST model for three validation cases: the incompressible flow over a smooth flat plate, the incompressible flow over a backward facing step, and the shock-induced flow separation inside a transonic diffuser. The k-epsilon model results indicate that the WIND model functions very similarly to that in NPARC, though the WIND code appears to he slightly more accurate in the treatment of the near-wall region. Comparisons of the k-epsilon model results with those from the SST model were less definitive, as each model exhibited strengths and weaknesses for each particular case.

The Capillary Pumped Loop 3 (CAPL 3) Experiment is a follow on to the CAPL 1 and CAPL 2 experiments which flew on STS-60 (2/94) and STS-69 (9/95), respectively. CAPL 3 is tentatively scheduled to fly on the Space Shuttle in late 2000 as part of the Hitchhiker Experiments Advancing Technology (HEAT) payload. The experiment is a joint Naval Research Laboratory (NRL)/Goddard Space Flight Center (GSFC) payload which will meet technology objectives for both the Department of Defense and NASA. The primary objective of CAPL 3 is to demonstrate in space a multiple evaporator capillary pumped loop system, capable of reliable start-up, reliable continuous operation, and at least 50% heat load sharing with hardware for a deployable radiator. CAPL 3 is a full scale CPL system with four parallel capillary evaporators. The loop also contains a capillary starter pump, 8 parallel direct condensation condensers with associated flow regulators, a back pressure regulator, a two-phase reservoir, and various headers and transport tubing. A variable conductance heat pipe is located between one of the evaporators and the experiment radiator to provide a cooling source for the demonstration of heat load sharing. The experiment has an operating power range of 100 W to approximately 1400 W. The experiment ammonia charge will cause it to transition to a fixed conductance mode of operation if the radiator usage reaches 85%. Ambient functional tests have been performed on the experiment. Tests performed included start-up, low power, power cycles, high power, heat load sharing, variable/fixed conductance transition, saturation temperature changes, and pressure primes while the system was operating. The majority of the testing was performed at an ammonia saturation temperature of 30C, but a few tests were done at temperatures above and below this. The testing was highly successful. Details of the tests performed and a discussion of the results will be given in the presentation.

The development of high-order accurate numerical discretization techniques for irregular domains and meshes is often cited as one of the remaining challenges facing the field of computational fluid dynamics. In structural mechanics, the advantages of high-order finite element approximation are widely recognized. This is especially true when high-order element approximation is combined with element refinement (h-p refinement). In computational fluid dynamics, high-order discretization methods are infrequently used in the computation of compressible fluid flow. The hyperbolic nature of the governing equations and the presence of solution discontinuities makes high-order accuracy difficult to achieve. Consequently, second-order accurate methods are still predominately used in industrial applications even though evidence suggests that high-order methods may offer a way to significantly improve the resolution and accuracy for these calculations. To address this important topic, a special course was jointly organized by the Applied Vehicle Technology Panel of NATO's Research and Technology Organization (RTO), the von Karman Institute for Fluid Dynamics, and the Numerical Aerospace Simulation Division at the NASA Ames Research Center. The NATO RTO sponsored course entitled "Higher Order Discretization Methods in Computational Fluid Dynamics" was held September 14-18, 1998 at the von Karman Institute for Fluid Dynamics in Belgium and September 21-25, 1998 at the NASA Ames Research Center in the United States. During this special course, lecturers from Europe and the United States gave a series of comprehensive lectures on advanced topics related to the high-order numerical discretization of partial differential equations with primary emphasis given to computational fluid dynamics (CFD). Additional consideration was given to topics in computational physics such as the high-order discretization of the Hamilton-Jacobi, Helmholtz, and elasticity equations. This volume consists of five articles prepared by the special course lecturers. These articles should be of particular relevance to those readers with an interest in numerical discretization techniques which generalize to very high-order accuracy. The articles of Professors Abgrall and Shu consider the mathematical formulation of high-order accurate finite volume schemes utilizing essentially non-oscillatory (ENO) and weighted essentially non-oscillatory (WENO) reconstruction together with upwind flux evaluation. These formulations are particularly effective in computing numerical solutions of conservation laws containing solution discontinuities. Careful attention is given by the authors to implementational issues and techniques for improving the overall efficiency of these methods. The article of Professor Cockburn discusses the discontinuous Galerkin finite element method. This method naturally extends to high-order accuracy and has an interpretation as a finite volume method. Cockburn addresses two important issues associated with the discontinuous Galerkin method: controlling spurious extrema near solution discontinuities via "limiting" and the extension to second order advective-diffusive equations (joint work with Shu). The articles of Dr. Henderson and Professor Schwab consider the mathematical formulation and implementation of the h-p finite element methods using hierarchical basis functions and adaptive mesh refinement. These methods are particularly useful in computing high-order accurate solutions containing perturbative layers and corner singularities. Additional flexibility is obtained using a mortar FEM technique whereby nonconforming elements are interfaced together. Numerous examples are given by Henderson applying the h-p FEM method to the simulation of turbulence and turbulence transition.

The intent of Stanford University's SciVis group is to develop technologies that enabled comparative analysis and visualization techniques for simulated and experimental flow fields. These techniques would then be made available un- der the Joint Research Interchange for potential injection into the DARWIN Workspace Environment (DWE). In the past, we have focused on techniques that exploited feature based comparisons such as shock and vortex extractions. Our current research effort focuses on finding a quantitative comparison of general vector fields based on topological features. Since the method relies on topological information, grid matching an@ vector alignment is not needed in the comparison. This is often a problem with many data comparison techniques. In addition, since only topology based information is stored and compared for each field, there is a significant compression of information that enables large databases to be quickly searched. This report will briefly (1) describe current technologies in the area of comparison techniques, (2) will describe the theory of our new method and finally (3) summarize a few of the results.

The problem of flow past bluff bodies was studied extensively in the past. The problem of drag reduction is very important in many high speed flow applications. Considerable work has been done in this subject area in case of circular cylinders. The present study attempts to investigate the feasibility of drag reduction on a rectangular cylinder by flow injection by flow injection from the rear stagnation region. The physical problem is modeled as two-dimensional body and numerical analysis is carried out with and without trailing jets. A commercial code is used for this purpose. Unsteady computation is performed in case of rectangular cylinders with no trailing jets where as steady state computation is performed when jet is introduced. It is found that drag can be reduced by introducing jets with small intensity in rear stagnation region of the rectangular cylinders.

The intent of Stanford University's SciVis group is to develop technologies that enabled comparative analysis and visualization techniques for simulated and experimental flow fields. These techniques would then be made available under the Joint Research Interchange for potential injection into the DARWIN Workspace Environment (DWE). In the past, we have focused on techniques that exploited feature based comparisons such as shock and vortex extractions. Our current research effort focuses on finding a quantitative comparison of general vector fields based on topological features. Since the method relies on topological information, grid matching and vector alignment is not needed in the comparison. This is often a problem with many data comparison techniques. In addition, since only topology based information is stored and compared for each field, there is a significant compression of information that enables large databases to be quickly searched. This report will briefly (1) describe current technologies in the area of comparison techniques, (2) will describe the theory of our new method and finally (3) summarize a few of the results.

The decay of anomalous effects on shock waves in weakly ionized gases following plasma generator extinction has been measured in the anticipation that the decay time must correlate well with the relaxation time of the mechanism responsible for the anomalous effects. When the relaxation times cannot be measured directly, they are inferred theoretically, usually assuming that the initial state is nearly in thermal equilibrium. In this paper, it is demonstrated that relaxation from any steady state far from equilibrium, including the state of a weakly ionized gas, can proceed much more slowly than arguments based on relaxation from near equilibrium states might suggest. This result justifies a more careful analysis of the relaxation times in weakly ionized gases and suggests that although the experimental measurements of relaxation times did not lead to an unambiguous conclusion, this approach to understanding the anomalous effects may warrant further investigation.

Aeroheating characteristics of the X-33 Rev-F configuration have been experimentally examined in the Langley 20-Inch Mach 6 Air Tunnel (Test 6770). Global surface heat transfer distributions, surface streamline patterns, and shock shapes were measured on a 0.013-scale model at Mach 6 in air. Parametric variations include angles-of-attack of 20-deg, 30-deg, and 40-deg; Reynolds numbers based on model length of 0.9 to 4.9 million; and body-flap deflections of 10-deg and 20-deg. The effects of discrete roughness elements on boundary layer transition, which included trip height, size, and location, both on and off the windward centerline, were investigated. This document is intended to serve as a quick release of preliminary data to the X-33 program; analysis is limited to observations of the experimental trends in order to expedite dissemination.

A thin film of oil on a surface responds primarily to the wall shear stress generated on that surface by a three-dimensional flow. The oil film is also subject to wall pressure gradients, surface tension effects and gravity. The partial differential equation governing the oil film flow is shown to be related to Burgers' equation. Analytical and numerical methods for solving the thin oil film equation are presented. A direct numerical solver is developed where the wall shear stress variation on the surface is known and which solves for the oil film thickness spatial and time variation on the surface. An inverse numerical solver is also developed where the oil film thickness spatial variation over the surface at two discrete times is known and which solves for the wall shear stress variation over the test surface. A One-Time-Level inverse solver is also demonstrated. The inverse numerical solver provides a mathematically rigorous basis for an improved form of a wall shear stress instrument suitable for application to complex three-dimensional flows. To demonstrate the complexity of flows for which these oil film methods are now suitable, extensive examination is accomplished for these analytical and numerical methods as applied to a thin oil film in the vicinity of a three-dimensional saddle of separation.

In this article, we present a reduced order modeling approach suitable for active control of fluid dynamical systems based on proper orthogonal decomposition (POD). The rationale behind the reduced order modeling is that numerical simulation of Navier-Stokes equations is still too costly for the purpose of optimization and control of unsteady flows. We examine the possibility of obtaining reduced order models that reduce computational complexity associated with the Navier-Stokes equations while capturing the essential dynamics by using the POD. The POD allows extraction of certain optimal set of basis functions, perhaps few, from a computational or experimental data-base through an eigenvalue analysis. The solution is then obtained as a linear combination of these optimal set of basis functions by means of Galerkin projection. This makes it attractive for optimal control and estimation of systems governed by partial differential equations. We here use it in active control of fluid flows governed by the Navier-Stokes equations. We show that the resulting reduced order model can be very efficient for the computations of optimization and control problems in unsteady flows. Finally, implementational issues and numerical experiments are presented for simulations and optimal control of fluid flow through channels.