ALBERT

Description: This paper presents new trends in Particle Image Velocimetry and practical aspects relevant to the application of the technique to large scale wind tunnel testing. The various problems and their solutions to the operation of PIV in large scale wind tunnels are discussed. Application of the technique in mapping complex flows are also presented.

Description: Pressure and temperature sensitive paints have been utilized for the measurement of blade surface pressure and temperature distributions in a high speed axial compressor and an Allied Signal F109 gas turbine engine. Alternate blades were painted with temperature sensitive paints and then pressure sensitive paint. This combination allows temperature distributions to be accounted for when determining the blade suction surface pressure distribution. Measurements were taken and pressure maps on the suction surface of a blade were obtained over a range of rotational speeds. Pressure maps of the suction surface show-strong shock waves at the higher speeds.

Description: The oil film interferometry skin-friction technique is described and applied to flows in some of the NASA Ames large wind tunnel facilities. Various schemes for applying the technique are discussed. Results are shown for tests in several wind tunnels which illustrate the oil film's ability to measure a variety of flow features such as shock waves separation, and 3D flow.

Description: Minuscule amounts (e.g., 1 part in 10,000) of a surface-active material in a liquid can drastically affect the surface behavior of the liquid, influencing how the material flows and mixes with other liquid and solid materials. In many respects, the science of surfactants has been empirical, with trial and error dominating over the ability to predict how surfactant type and concentration influence surface behavior. A program for the modeling of surfactant behavior has been established at Yale. This program combines experimental work performed both on the ground and in space, and theoretical and numerical modeling. By levitating a drop of liquid in air, away from solid container surfaces, and by manipulating the drop with acoustic radiation forces, we have been able to establish idealized conditions for surface behavior studies. The primary experiments involve the study of the free oscillations of initially deformed drops. In STS-73, the USML-2 mission of the Space Shuttle, we performed the following measurements: 1) the oscillation of a spherical drop in its quadrupole mode; 2) the oscillation of a drop about a deformed (oblate) shape; 3) the slow static squeezing of the drop from spherical to nearly flat; and 4) the superoscillations of drops when the radiation forces maintaining the drop in a flattened state are suddenly reduced. Analytic and numerical studies have enabled us to understand the physics of these oscillations and to extract material properties such as the dynamic surface tension and the surface viscosities (shear and dilatational). The relation to ground-based studies is essential, because the knowledge and understanding gleaned from our space studies enable us to interpret ground-based data.

Description: While turbulence models are being developed which show promise for simulating the transition region on a turbine blade or vane, it is believed that the best approach with the greatest potential for practical use is the use of models which incorporate the physics of turbulent spots present in the transition region. This type of modeling results in the prediction of transition region intermittency which when incorporated in turbulence models give a good to excellent prediction of the transition region heat transfer. Some models are presented which show how turbulent spot characteristics and behavior can be employed to predict the effect of pressure gradient and Mach number on the transition region. The models predict the spot formation rate which is needed, in addition to the transition onset location, in the Narasimha concentrated breakdown intermittency equation. A simplified approach is taken for modeling turbulent spot growth and interaction in the transition region which utilizes the turbulent spot variables governing transition length and spot generation rate. The models are expressed in terms of spot spreading angle, dimensionless spot velocity, dimensionless spot area, disturbance frequency and Mach number. The models are used in conjunction with a computer code to predict the effects of pressure gradient and Mach number on the transition region and compared with VKI experimental turbine data.

Description: In recent years the techniques of computational fluid dynamics (CFD) have been used to compute flows associated with geometrically complex configurations. However, success in terms of accuracy and reliability has been limited to cases where the effects of turbulence and transition could be modeled in a straightforward manner. Even in simple flows, the accurate computation of skin friction and heat transfer using existing turbulence models has proved to be a difficult task, one that has required extensive fine-tuning of the turbulence models used. In more complex flows (for example, in turbomachinery flows in which vortices and wakes impinge on airfoil surfaces causing periodic transitions from laminar to turbulent flow) the development of a model that accounts for all scales of turbulence and predicts the onset of transition is an extremely difficult task. Fortunately, current trends in computing suggest that it may be possible to perform direct simulations of turbulence and transition at moderate Reynolds numbers in some complex cases in the near future. This presentation will focus on direct simulations of transition and turbulence using high-order accurate finite-difference methods. The advantage of the finite-difference approach over spectral methods is that complex geometries can be treated in a straightforward manner. Additionally, finite-difference techniques are the prevailing methods in existing application codes. An application of high-order-accurate finite-difference methods to direct simulations of transition and turbulence in a spatially evolving boundary layer subjected to high levels of freestream turbulence will be presented.

Description: During hypersonic flight, high temperatures and high heat fluxes are generated. The Flight Loads Laboratory (FLL) at Dryden Flight Research Center (DFRC) is equipped to calibrate high heat fluxes up to 1100 kW/sq m. There are numerous uncertainties associated with these heat flux calibrations, as the process is transient, there are expected to be interactions between transient conduction, natural and forced convection, radiation, and possibly an insignificant degree of oxidation of the graphite cavity. Better understanding, of these mechanisms during the calibration process, will provide more reliable heat transfer data during either ground testing or flight testing of hypersonic vehicles.

Description: 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 non-axisymmetric 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.

Description: An overview of the NASA Lewis Low-Pressure Turbine (LPT) Flow Physics Program will be presented. The program was established in response to the aero-engine industry's need for improved LPT efficiency and designs. Modern jet engines have four to seven LPT stages, significantly contributing to engine weight. In addition, there is a significant efficiency degradation between takeoff and cruise conditions, of up to 2 points. Reducing the weight and part count of the LPT and minimizing the efficiency degradation will translate into fuel savings. Accurate prediction methods of LPT flows and losses are needed to accomplish those improvements. The flow in LPT passages is at low Reynolds number, and is dominated by interplay of three basic mechanisms: transition, separation and wake interaction. The affecting parameters traditionally considered are Reynolds number, freestream turbulence intensity, wake frequency parameter, and the pressure distribution (loading). Three-dimensional effects and additional parameters, particularly turbulence characteristics like length scales, spectra and other statistics, as well as wake turbulence intensity and properties also play a role. The flow of most interest is on the suction surface, where large losses are generated as the flow tends to separate at the low Reynolds numbers. Ignoring wakes, a common flow scenario, there is laminar separation, followed by transition on the separation bubble and turbulent reattachment. If transition starts earlier the separation will be eliminated and the boundary layer will be attached leading to the well known bypass transition issues. In contrast, transition over a separation bubble is closer to free shear layer transition and was not investigated as well, particularly in the turbine environment. Unsteadiness created by wakes complicates the picture. Wakes induce earlier transition, and the calmed regions trailing the induced turbulent spots can delay or eliminate separation via shear stress modification. Three-dimensional flow physics and geometry will have strong effects. Altogether a very complex and challenging problem emerges. The objective of the program is to provide improved models and physical understanding of the complex flow, which are essential for accurate prediction of flow and losses in the LPT. Experimental, computational and analytical work as complementing and augmenting approaches are used. The program involves industry, universities and research institutes, and other government laboratories. It is characterized by strong interaction among participants, quick dissemination of results, and responsiveness to industry's needs. The presentation will describe the work elements. Highlighting some activities in progress are experiments on simulated blade suction surface in low-speed wind tunnels, on curved wall, and on a flat-plate, both with pressure gradient. In the area of computation, assessment of existing models is performed using RANS (Reynolds Averaged Navier Stokes) simulations. Laminar flow DNS was completed. Analytical studies of instability and receptivity in attached and separated flows were started. In the near future the program is moving to include wake effects and development of improved modeling. Experimental work in preparation stages are: (1) Addition of wakes to the curved tunnel experiment; (2) Low-speed rotating rig experiment on GE90 engine LPT; and (3) Transonic cascade. In the area of computation, it is expected to move from model assessment towards development of improved models. In addition, a new project of Large Eddy Simulation (LES) of LPT is to begin and will provide numerical data bases. It is planned to implement the emerging improved models in a multistage turbomachinery code and to validate against the GE90 engine LPT.

Description: Rarefaction effects are important for hypersonic applications for a wide spectrum of conditions ranging from low-density (high altitude) situations to relatively high-density flows where the characteristic dimension is small. The present chapter concentrates on two hypersonic flow problems at flow conditions that produce a significant range of rarefaction effects: corner flow with jet interaction and blunt body flow with special emphasis on the near wake, These problems were chosen because they involve complex flow interactions that have significant implications for both spacecraft and re-entry vehicles. In an effort to clarify issues associated with these two general flow problems and to enhance their respective databases, both experimental and computational contributions were executed by an international group of researchers. In some cases, multiple data sources for both experimental and computational contributions are achieved.

Description: In this chapter recent activity in real-gas database definition and code validation will be summarized. In the Phase I report of the Working Group (WG) 181, aerothermodynamic problems were classified, for purpose of discussion, into seven types: aerodynamic parameters, viscous/shock interaction, boundary-layer transition, forebody-heating/heat-transfer, radiation and ablation, lee and base-region flow, and low-density flow. Several of these problem types were the subject of various chapters of the Phase 1 report describing real-gas effects and ground test facility issues. In this chapter some background and objectives outlined in the real-Gas effects Chapter V of the Phase 1 report will be reviewed. The results of the blunt cone test campaign developed under the auspices of the WG18 activity to study real-gas phenomena will be summarized, including the experimental and computational programs, issues and questions, and recommendations. Further, recent progress in other real-gas areas beyond the blunt cone test campaign will be discussed. Finally, a summary in which the present status of our understanding of real-gas issues will be presented.

Description: The operation of Capillary Pumped Loops (CPL's) in low gravity has generally been unable to match ground-based performance. The reason for this poorer performance has been elusive. In order to investigate the behavior of a CPL in low-gravity, an idealized, glass CPL experiment was constructed. This experiment, known as the Capillary-driven Heat Transfer (CHT) experiment, was flown on board the Space Shuttle Columbia in July 1997 during the Microgravity Science Laboratory mission. During the conduct of the CHT experiment an unexpected failure mode was observed. This failure mode was a result of liquid collecting and then eventually bridging the vapor return line. With the vapor return line blocked, the condensate was unable to return to the evaporator and dry-out subsequently followed. The mechanism for this collection and bridging has been associated with long wavelength instabilities of the liquid film forming in the vapor return line. Analysis has shown that vapor line blockage in present generation CPL devices is inevitable. Additionally, previous low-gravity CPL tests have reported the presence of relatively low frequency pressure oscillations during erratic system performance. Analysis reveals that these pressure oscillations are in part a result of long wavelength instabilities present in the evaporator pores, which likewise lead to liquid bridging and vapor entrapment in the porous media. Subsequent evaporation to the trapped vapor increases the vapor pressure. Eventually the vapor pressure causes ejection of the bridged liquid. Recoil stresses depress the meniscus, the vapor pressure rapidly increases, and the heated surface cools. The process then repeats with regularity.

Description: Heinzl et al. (1985) reports that experiments in ink-jets to produce drawings or signals occurred as early as 1930. Various companies such as IBM and Pitney-Bowes have conducted extensive studies on these devices for many years. Many such reports are available in such journals as the IBM Journal of Research and Development. While numerous articles have been published on the jetting characteristics of ink and water, the literature is rather limited on fluids such as waxes (Gao & Sonin 1994) or non-water based fluids (Passow, et al. 1993). This present study extends the knowledge base to determine the performance of molten waxes in "ink-jet" type printers for rapid prototyping. The purpose of this research was to qualitatively and quantitatively study the droplet formation of a drop-on-demand ink-jet type nozzle system for rapid prototyping.

Description: Capillary Pumped Loops (CPL's) and Loop Heat Pipes (LHP's) are finding increased acceptance on upcoming NASA spacecraft missions, as well as military and commercial applications, The transition from a research and development tool to an "off the shelf' system is underway. The state of the art of two phase systems (TPS) is reviewed and applications on various NASA missions are discussed, with particular emphasis on new technology developments. Upcoming research areas and flight experiments are also addressed, along with recommendations for future activities.

Description: This paper details the flight configuration and pre-flight performance test results of the fifth generation cryogenic capillary pumped loop (CCPL-5). This device will fly on STS-95 in October 1998 as part of the CRYOTSU Flight Experiment. This flight represents the first in-space demonstration of a CCPL, a miniaturized two-phase fluid circulator for thermally linking cryogenic cooling sources to remote cryogenic components. CCPL-5 utilizes N2 as the working fluid and has a practical operating range of 75-110 K. Test results indicate that CCPL-5, which weighs about 200 grams, can transport over 10 W of cooling a distance of 0.25 m (or more) with less than a 5 K temperature drop.

Description: Capillary Pumped Loops (CPL's) and Loop Heat Pipes (LHP's) are finding increased acceptance on upcoming NASA spacecraft missions, as well as military and commercial applications. The transition from a research and development tool to an "off the shelf' system is underway. The state of the art of two phase systems (TPS) is reviewed and applications on various NASA missions are discussed, with particular emphasis on new technology developments. Upcoming research areas and flight experiments are also addressed, along with recommendations for future activities.

Description: This paper details the flight configuration and pre-flight performance test results of the fifth generation cryogenic capillary pumped loop (CCPL-5). This device will fly on STS-95 in October 1998 as part of the CRYOTSU Flight Experiment. This flight represents the first in-space demonstration of a CCPL; a miniaturized two-phase fluid circulator for thermally linking cryogenic components. CCPL-5 utilizes N2 as the working fluid and has a practical operating range of 75-110 K. Test results indicate that CCPL-5, which weighs about 200 grams, can transport over 10 W of cooling a distance of 0.25 m (or more) with less than a 5 K temperature drop.

Description: Planar Doppler Velocimetry (PDV) concepts using a pulsed laser are described and the obtainable minimum resolved velocities in large-scale wind tunnels are evaluated. Velocity-field measurements are shown to be possible at ranges of tens of meters and with single pulse resolutions as low as 2 m/s. Velocity measurements in the flow of a low-speed, turbulent jet are reported that demonstrate the ability of PDV to acquire both average velocity fields and their fluctuation amplitudes, using procedures that are compatible with large-scale facility operations. The advantages of PDV over current Laser Doppler Anemometry and Particle Image Velocimetry techniques appear to be significant for applications to large facilities.

Description: Experimental results for the spreading and centerline velocity decay rates for round, compressible jets, from a convergent and a convergent-divergent nozzle, are presented. The spreading rate is determined from the variation of streamwise mass flux obtained from Pitot probe surveys. Results for the far asymptotic region show that both spreading and centerline velocity decay rates, when nondimensionalized by parameters at the nozzle exit, decrease with increasing "jet Mach number" M(sub J). Dimensional analysis with the assumption of momentum conservation, together with compressible flow calculations for the conditions at the nozzle exit, predict this Mach number dependence well. The analysis also demonstrates that an increase in the "potential core length" of the jet occurring with increasing M(sub J), a commonly observed trend, is largely accounted for simply by the variations in the density and static pressure at the nozzle exit. The effect of decreasing mixing efficiency with increasing compressibility is inferred to contribute only partially to the latter trend.

Description: A three-dimensional Navier-Stokes code has been used to compare the heat transfer coefficient on a film-cooled, rotating turbine blade. The blade chosen is the ACE rotor with five rows containing 93 film cooling holes covering the entire span. This is the only film-cooled rotating blade over which experimental data is available for comparison. Over 2.278 million grid points are used to compute the flow over the blade including the tip clearance region, using Coakley's q-omega turbulence model. Results are also compared with those obtained by Garg and Abhari (1997) using the zero-equation Baldwin-Lomax (B-L) model. A reasonably good comparison with the experimental data is obtained on the suction surface for both the turbulence models. At the leading edge, the B-L model yields a better comparison than the q-omega model. On the pressure surface, however, the comparison between the experimental data and the prediction from either turbulence model is poor. A potential reason for the discrepancy on the pressure surface could be the presence of unsteady effects due to stator-rotor interaction in the experiments which are not modeled in the present computations. Prediction using the two-equation model is in general poorer than that using the zero-equation model, while the former requires at least 40% more computational resources.

Description: Stationary onset of convection due to surface tension variation in an unbounded multicomponent fluid layer is considered. Surface deformation is included and general flux boundary conditions are imposed on the stratifying agencies (temperature/composition) disturbance equations. Exact solutions are obtained to the general N-component problem for both finite and infinitesimal wavenumbers. Long wavelength instability may coexist with a finite wavelength instability for certain sets of parameter values, often referred to as frontier points. For an impermeable/insulated upper boundary and a permeable/conductive lower boundary, frontier boundaries are computed in the space of Bond number, Bo, versus Crispation number, Cr, over the range 5 x 10(exp -7) less than or equal to Bo less than or equal to 1. The loci of frontier points in (Bo, Cr) space for different values of N, diffusivity ratios, and, Marangoni numbers, collapsed to a single curve in (Bo, D(dimensional variable)Cr) space, where D(dimensional variable) is a Marangoni number weighted diffusivity ratio.

Description: Radiative equilibrium surface temperatures and surface heating rates from a combined inviscid-boundary layer method are presented for the X-34 Reusable Launch Vehicle for several points along the hypersonic descent portion of its trajectory. Inviscid, perfect-gas solutions are generated with the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) and the Data-Parallel Lower-Upper Relaxation (DPLUR) code. Surface temperatures and heating rates are then computed using the Langley Approximate Three-Dimensional Convective Heating (LATCH) engineering code employing both laminar and turbulent flow models. The combined inviscid-boundary layer method provides accurate predictions of surface temperatures over most of the vehicle and requires much less computational effort than a Navier-Stokes code. This enables the generation of a more thorough aerothermal database which is necessary to design the thermal protection system and specify the vehicle's flight limits.

Description: A higher order accurate numerical procedure has been developed for solving incompressible Navier-Stokes equations for 2D or 3D fluid flow problems. It is based on low-storage Runge-Kutta schemes for temporal discretization and fourth and sixth order compact finite-difference schemes for spatial discretization. The particular difficulty of satisfying the divergence-free velocity field required in incompressible fluid flow is resolved by solving a Poisson equation for pressure. It is demonstrated that for consistent global accuracy, it is necessary to employ the same order of accuracy in the discretization of the Poisson equation. Special care is also required to achieve the formal temporal accuracy of the Runge-Kutta schemes. The accuracy of the present procedure is demonstrated by application to several pertinent benchmark problems.

Description: An experimental study was made to obtain quantitative information on heat transfer, flow, and pressure distribution in a branched duct test section that had several significant features of an internal cooling passage of a turbine blade. The objective of this study was to generate a set of experimental data that could be used for validation of computer codes that would be used to model internal cooling. Surface heat transfer coefficients and entrance flow conditions were measured at nominal entrance Reynolds numbers of 45,000, 335,000, and 726,000. Heat transfer data were obtained by using a steady-state technique in which an Inconel heater sheet is attached to the surface and coated with liquid crystals. Visual and quantitative flow-field data from particle image velocimetry measurements for a plane at midchannel height for a Reynolds number of 45,000 were also obtained. The flow was seeded with polystyrene particles and illuminated by a laser light sheet. Pressure distribution measurements were made both on the surface with discrete holes and in the flow field with a total pressure probe. The flow-field measurements yielded flow-field velocities at selected locations. A relatively new method, pressure sensitive paint, was also used to measure surface pressure distribution. The pressure paint data obtained at Reynolds numbers of 335,000 and 726,000 compared well with the more standard method of measuring pressures by using discrete holes.

Description: This project focused on the generation and suppression of oscillatory thermocapillary convection in a thin liquid layer. The bulk of the research was experimental in nature, some theoretical work was also done. ne first phase of this research generated, for the first time, the hydrothermal-wave instability predicted by Smith and Davis in 1983. In addition, the behavior of the fluid layer under a number of conditions was investigated and catalogued. A transition map for the instability of buoyancy-thermocapillary convection was prepared which presented results in terms of apparatus-dependent and apparatus-independent parameters, for ease of comparison with theoretical results. The second phase of this research demonstrated the suppression of these hydrothermal waves through an active, feed-forward control strategy employing a CO2 laser to selectively heat lines of negative disturbance temperature on the free surface of the liquid layer. An initial attempt at this control was only partially successful, employing a thermocouple inserted slightly below the free surface of the liquid to generate the control scheme. Subsequent efforts, however, were completely successful in suppressing oscillations in a portion of the layer by utilizing data from an infrared image of the free surface to compute hydrothermal-wave phase speeds and, using these, to tailor the control scheme to each passing wave.

Description: Analytical studies of compressible turbulence have found that compressible velocity fluctuations create both effective fluid transport properties and an effective equation of state. This paper investigates the effects of rotation on compressible turbulence. It is shown that rotation modifies the transport properties of compressible turbulence by replacing the turbulence time scale by a rotational time scale, much as rotation modifies the transport properties of incompressible turbulence. But thermal equilibrium properties are modified in a more complex manner. Two regimes are possible: one dominated by incompressible fluctuations, in which the sound speed is modified as it is in non-rotating turbulence, and a rotation dominated regime in which the sound speed enhancement is rotation dependent. The dimensionless parameter which discriminates between regimes is identified. In general, rotation is found to suppress the effects of compressibility. A novel feature of the present analysis is the use of a non-Kolmogorov steady state as the reference state of turbulence. introduction of such steady states expands the power and utility of analytical turbulence closures to a wider range of problems.

Description: The structure of a numerically simulated turbulent boundary layer over a flat plate at Re(theta) = 670 was studied using the invariants of the velocity gradient tensor (Q and R) and a related scalar quantity, the cubic discriminant (D = 27R(exp 2)/4 + Q(exp 3)). These invariants have previously been used to study the properties of the small-scale motions responsible for the dissipation of turbulent kinetic energy. In addition, these scalar quantities allow the local flow patterns to be unambiguously classified according to the terminology proposed by Chong et al. The use of the discriminant as a marker of coherent motions reveals complex, large-scale flow structures that are shown to be associated with the generation of Reynolds shear stress -u'v'(bar). These motions are characterized by high spatial gradients of the discriminant and are believed to be an important part of the mechanism that sustains turbulence in the near-wall region.

Description: Direct Numerical Simulation (DNS) is used to study the mechanism of generation and evolution of turbulence structures in a temporally evolving supersonic swirling round jet and also to examine the resulting acoustic radiations. Fourier spectral expansions are used in the streamwise and azimuthal directions and a 1-D b-spline Galerkin representation is used in the radial direction. Spectral-like accuracy is achieved using this numerical scheme. Direct numerical simulations, using the b-spline spectral method, are carried out starting from mean flow initial conditions which are perturbed by the most unstable linear stability eigenfunctions. It is observed that the initial helical instability waves evolve into helical vortices which eventually breakdown into smaller scales of turbulence. 'Rib' structures similar to those seen in incompressible mixing layer flow of Rogers and Moserl are observed. The jet core breakdown stage exhibits increased acoustic radiations.

Description: The start up of a liquid metal heat pipe from the frozen state was evaluated experimentally with a Nb-1%Zr heat pipe with potassium as the working fluid. The heat pipe was fabricated and tested at Los Alamos National Laboratory. RF induction heating was used to heat 13 cm of the 1-m-long heat pipe. The heat pipe and test conditions are well characterized so that the test data may be used for comparison with numerical analyses. An attempt was made during steady state tests to calibrate the heat input so that the heat input would be known during the transient cases. The heat pipe was heated to 675 C with a throughput of 600 W and an input heat flux of 6 W/cm(exp 2). Steady state tests, start up from the frozen state, and transient variations from steady state were performed.

Description: Refractory-composite/heat-pipe-cooled wing an tail leading edges are being considered for use on hypersonic vehicles to limit maximum temperatures to values below material reuse limits and to eliminate the need to actively cool the leading edges. The development of a refractory-composite/heat-pipe-cooled leading edge has evolved from the design stage to the fabrication and testing of heat pipes embedded in carbon/carbon (C/C). A three-foot-long, molybdenum-rhenium heat pipe with a lithium working fluid was fabricated and tested at an operating temperature of 2460 F to verify the individual heat-pipe design. Following the fabrication of this heat pipe, three additional heat pipes were fabricated and embedded in C/C. The C/C heat-pipe test article was successfully tested using quartz lamps in a vacuum chamber in both a horizontal and vertical orientation. Start up and steady state data are presented for the C/C heat-pipe test article. Radiography and eddy current evaluations were performed on the test article.

Description: In support of the wake vortex effect of the Terminal Area Productivity program, we have put forward four tasks to be accomplished in our proposal. The first task is validation of two-dimensional wake vortex-turbulence interaction. The second task is investigation of three-dimensional interaction between wake vortices and atmospheric boundary layer (ABL) turbulence. The third task is ABL studies. The, fourth task is addition of a Klemp-Durran condition at the top boundary for TASS model. The accomplishment of these tasks will increase our understanding of the dynamics of wake vortex and improve forecasting systems responsible for air safety and efficiency. The first two tasks include following three parts: (a) Determine significant length scale for vortex decay and transport, especially the length scales associated with the onset of Crow instability (Crow, 1970); (b) Study the effects of atmospheric turbulence on the decay of the wake vortices; and (c) Determine the relationships between decay rate, transport properties and atmospheric parameters based on large eddy simulation (LES) results and the observational data. These parameters may include turbulence kinetic energy, dissipation rate, wind shear and atmospheric stratification. The ABL studies cover LES modeling of turbulence structure within planetary boundary layer under transition and stable stratification conditions. Evidences have shown that the turbulence in the stable boundary layer can be highly intermittent and the length scales of eddies are very small compared to those in convective case. We proposed to develop a nesting grid mesh scheme and a modified Klemp-Durran conditions (Klemp and Wilhelmson, 1978) at the top boundary for TASS model to simulate planetary boundary layer under stable stratification conditions. During the past year, our group has made great efforts to carry out the above mentioned four tasks simultaneously. The work accomplished in the last year will be described in the next section.

Description: For weakly compressible turbulent fluctuations there exists a small parameter, the square of the fluctuating Mach number, that allows an investigation using a perturbative treatment. The consequences of such a perturbative analysis in three different subject areas are described: 1) initial conditions in direct numerical simulations, 2) an explanation for the oscillations seen in the compressible pressure in the direct numerical simulations of homogeneous shear, and 3) for turbulence closures accounting for the compressibility of velocity fluctuations. Initial conditions consistent with small turbulent Mach number asymptotics are constructed. The importance of consistent initial conditions in the direct numerical simulation of compressible turbulence is dramatically illustrated: spurious oscillations associated with inconsistent initial conditions are avoided, and the fluctuating dilatational field is some two orders of magnitude smaller for a compressible isotropic turbulence. For the isotropic decay it is shown that the choice of initial conditions can change the scaling law for the compressible dissipation. A two-time expansion of the Navier-Stokes equations is used to distinguish compressible acoustic and compressible advective modes. A simple conceptual model for weakly compressible turbulence - a forced linear oscillator is described. It is shown that the evolution equations for the compressible portions of turbulence can be understood as a forced wave equation with refraction. Acoustic modes of the flow can be amplified by refraction and are able to manifest themselves in large fluctuations of the compressible pressure.

Description: A simple kinematically consistent expression for the longitudinal two-point correlation function related to both the integral length scale and the Taylor microscale is obtained. On the inner scale, in a region of width inversely proportional to the turbulent Reynolds number, the function has the appropriate curvature at the origin. The expression for two-point correlation is related to the nonlinear cascade rate, or dissipation epsilon, a quantity that is carried as part of a typical single-point turbulence closure simulation. Constructing an expression for the two-point correlation whose curvature at the origin is the Taylor microscale incorporates one of the fundamental quantities characterizing turbulence, epsilon, into a model for the two-point correlation function. The integral of the function also gives, as is required, an outer integral length scale of the turbulence independent of viscosity. The proposed expression is obtained by kinematic arguments; the intention is to produce a practically applicable expression in terms of simple elementary functions that allow an analytical evaluation, by asymptotic methods, of diverse functionals relevant to single-point turbulence closures. Using the expression devised an example of the asymptotic method by which functionals of the two-point correlation can be evaluated is given.

Description: A nearly homogeneous nearly isotropic compressible turbulent flow interacting with a normal shock wave has been studied experimentally in a large shock tube facility. Spatial resolution of the order of 8 Kolmogorov viscous length scales was achieved in the measurements of turbulence. A variety of turbulence generating grids provide a wide range of turbulence scales. Integral length scales were found to substantially decrease through the interaction with the shock wave in all investigated cases with flow Mach numbers ranging from 0.3 to 0.7 and shock Mach numbers from 1.2 to 1.6. The outcome of the interaction depends strongly on the state of compressibility of the incoming turbulence. The length scales in the lateral direction are amplified at small Mach numbers and attenuated at large Mach numbers. Even at large Mach numbers amplification of lateral length scales has been observed in the case of fine grids. In addition to the interaction with the shock the present work has documented substantial compressibility effects in the incoming homogeneous and isotropic turbulent flow. The decay of Mach number fluctuations was found to follow a power law similar to that describing the decay of incompressible isotropic turbulence. It was found that the decay coefficient and the decay exponent decrease with increasing Mach number while the virtual origin increases with increasing Mach number. A mechanism possibly responsible for these effects appears to be the inherently low growth rate of compressible shear layers emanating from the cylindrical rods of the grid.

Description: Radiative equilibrium surface temperatures, heating rates, streamlines, surface pressures, and flow-field features as predicted by the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) are presented for the X-34 Technology Demonstrator. Results for two trajectory points corresponding to entry peak heating and two control surface deflections are discussed. This data is also discussed in the context of Thermal Protection System (TPS) design issues. The work presented in this report is part of a larger effort to define the X-34 aerothermal environment, including the application of engineering codes and wind-tunnel studies.

Description: We studied the interaction of a vortical density inhomogeneity (VDI) with shock and expansion waves. We call the VDI the region of concentrated vorticity (vortex) with a density different from that of ambiance. Non-parallel directions of the density gradient normal to the VDI surface and the pressure gradient across a shock wave results in an additional vorticity. The roll-up of the initial round VDI towards a non-symmetrical shape is studied numerically. Numerical modeling of this interaction is performed by a 2-D Euler code. The use of an adaptive unstructured numerical grid makes it possible to obtain high accuracy and capture regions of induced vorticity with a moderate overall number of mesh points. For the validation of the code, the computational results are compared with available experimental results and good agreement is obtained. The interaction of the VDI with a propagating shock wave is studied for a range of initial and induced circulations and obtained flow patterns are presented. The splitting of the VDI develops into the formation of a non-symmetrical vortex pair and not in a set of vortices. A method for the analytical computation of an overall induced circulation Gamma(sub 1) as a result of the interaction of a moving VDI with a number of waves is proposed. Simplified, approximated, expressions for Gamma(sub 1) are derived and their accuracy is discussed. The splitting of the VDI passing through the Prandtl-Meyer expansion wave is studied numerically. The obtained VDI patterns are compared to those for the interaction of the VDI with a propagating shock wave for the same values of initial and induced circulations. These patterns have similar shapes for corresponding time moments.

Description: The topics contained in this progress report are direct numerical simulation of turbulent non-premixed combustion with realistic chemistry; LES of non-premixed turbulent reacting flows with conditional source term estimation; measurements of the three-dimensional scalar dissipation rate in gas-phase planar turbulent jets; direct simulation of a jet diffusion flame; on the use of interpolating wavelets in the direct numerical simulation of combustion; on the use of a dynamically adaptive wavelet collocation algorithm in DNS (direct numerical simulation) of non-premixed turbulent combustion; 2D simulations of Hall thrusters; computation of trailing-edge noise at low mach number using LES and acoustic analogy; weakly nonlinear modeling of the early stages of bypass transition; interactions between freestream turbulence and boundary layers; interfaces at the outer boundaries of turbulent motions; largest scales of turbulent wall flows; the instability of streaks in near-wall turbulence; an implementation of the v(sup 2) - f model with application to transonic flows; heat transfer predictions in cavities; a structure-based model with stropholysis effects; modeling a confined swirling coaxial jet; subgrid-scale models based on incremental unknowns for large eddy simulations; subgrid scale modeling taking the numerical error into consideration; towards a near-wall model for LES of a separated diffuser flow; on the feasibility of merging LES with RANS (Reynolds Averaging Numerical simulation) for the near-wall region of attached turbulent flows; large-eddy simulation of a separated boundary layer; numerical study of a channel flow with variable properties; on the construction of high order finite difference schemes on non-uniform meshes with good conservation properties; development of immersed boundary methods for complex geometries; and particle methods for micro and macroscale flow simulations.

Description: The Strutjet approach to Rocket Based Combined Cycle (RBCC) propulsion depends upon fuel-rich flows from the rocket nozzles and turbine exhaust products mixing with the ingested air for successful operation in the ramjet and scramjet modes. It is desirable to delay this mixing process in the air-augmented mode of operation present during low speed flight. A model of the Strutjet device has been built and is undergoing test to investigate the mixing of the streams as a function of distance from the Strutjet exit plane during simulated low speed flight conditions. Cold flow testing of a 1/6 scale Strutjet model is underway and nearing completion. Planar Laser Induced Fluorescence (PLIF) diagnostic methods are being employed to observe the mixing of the turbine exhaust gas with the gases from both the primary rockets and the ingested air simulating low speed, air augmented operation of the RBCC. The ratio of the pressure in the turbine exhaust duct to that in the rocket nozzle wall at the point of their intersection is the independent variable in these experiments. Tests were accomplished at values of 1.0, 1.5 and 2.0 for this parameter. Qualitative results illustrate the development of the mixing zone from the exit plane of the model to a distance of about 19 equivalent rocket nozzle exit diameters downstream. These data show the mixing to be confined in the vertical plane for all cases, The lateral expansion is more pronounced at a pressure ratio of 1.0 and suggests that mixing with the ingested flow would be likely beginning at a distance of 7 nozzle exit diameters downstream of the nozzle exit plane.

Description: A theoretical explanation of some experimentally observed phenomena associated with the so-called Klebanoff modes is obtained by analyzing the flow over a finite thickness flat plate resulting from a small-amplitude distortion imposed on the upstream mean flow. The analysis shows (among other things) how the stretching of the vortex lines around the plate leads to streamwise vorticity at the plate surface, which then produces a streamwise velocity perturbation within the boundary layer that can be related to the experimentally observed Klebanoff mode. The complete evolution of this flow must be found by solving the boundary-region equations of Kemp (1951) and Davis and Rubin (1980), but a limiting analytical solution can also be obtained. Since the initial growth of the boundary-layer disturbance is nearly algebraic, our results demonstrate how the algebraically growing disturbances promoted by Landahl and others can be generated by a realistic external-disturbance environment. The relationship between these results and various bypass transition mechanisms is discussed.

Description: Six months of funding was received for the proposed three year research program (funding for the period from March 1, 1997 to August 31, 1997). Although the official starting date for the project was March 1, 1997, no funding for the project was received until July 1997. In the funded research period, considerable progress was made on Phase I of the proposed research program. The initial research efforts concentrated on applying the 10-, 20-, and 35-moment Gaussian-based closures to a series of standard two-dimensional non-reacting single species test flow problems, such as the flat plate, couette, channel, and rearward facing step flows, and to some other two-dimensional flows having geometries similar to those encountered in chemical-vapor deposition (CVD) reactors. Eigensystem analyses for these systems for the case of two spatial dimensions was carried out and efficient formulations of approximate Riemann solvers have been formulated using these eigenstructures. Formulations to include rotational non-equilibrium effects into the moment closure models for the treatment of polyatomic gases were explored, as the original formulations of the closure models were developed strictly for gases composed of monatomic molecules. The development of a software library and computer code for solving relaxing hyperbolic systems in two spatial dimensions of the type arising from the closure models was also initiated. The software makes use of high-resolution upwind finite-volumes schemes, multi-stage point implicit time stepping, and automatic adaptive mesh refinement (AMR) to solve the governing conservation equations for the moment closures. The initial phase of the code development was completed and a numerical investigation of the solutions of the 10-moment closure model for the simple two-dimensional test cases mentioned above was initiated. Predictions of the 10-moment model were compared to available theoretical solutions and the results of direct-simulation Monte Carlo (DSMC) calculations. The first results of this study were presented at a meeting last year.

Description: The Marshall Space Flight Center (MSFC) is the Center of Excellence for Space Propulsion Technology within NASA. One of MSFC's goals is to lower the cost of access to space. MSFC's Fluid Dynamics Division, a part of the Structures and Dynamics Lab, assists in accomplishing this goal. This division is divided into two branches: experimental and analytical, whose work is complimentary and integrated. In the analysis branch, the emphasis is to apply, rather than develop computational fluid dynamics (CFD), to reduce the cost of design and development of: (1) pump and turbine feedlines and feed manifolds; (2) all primary and secondary flow paths in pumps; (3) all primary and secondary flow paths in turbines; (4) combustion chambers; (5) rocket engine nozzles and plumes; and (6) launch vehicles ascent. The broad range of applications imposes unique constraints on the software. The use of CFD in evaluating rocket engine hardware, specifically, is of value due to several reasons: (1) Rocket engine hardware is extremely expensive to build and test; (2) high engineering and development cost is due to limited production; (3) the cost of one engine test is greater than the average yearly hardware and software cost; (4) CFD represents a relatively small investment relative to the impact it has on the engine development; and (5) computing costs are continuously decreasing. In terms of accuracy the predictable error is more useful than the inconsistent accuracy. The achievable accuracy varies from application to application. The lack of accuracy is often caused by improper analyses: such as improperly applied or located boundary conditions, or attempting steady analysis of unsteady flows. This report is in viewgraph format.

Description: The heat transport capability of a capillary pumped loop (CPL) is limited by the pressure drop that its evaporator wick can sustain. The pressure drop in a CPL is not constant even under seemingly steady operation, but rather exhibits an oscillatory behavior. A hydrodynamic theory based on a mass-spring-dashpot model was previously developed to predict the pressure oscillation in a CPL with a single evaporator and a single condenser. The theory states that the pressure oscillation is a function of physical dimensions of the CPL components and operating conditions. Experimental data agreed very well with theoretical predictions. The hydrodynamic stability theory has recently been extended to predict the pressure oscillations in CPLs with multiple evaporators and multiple condensers. Concurrently, an experimental study was conducted to verify the theory and to investigate the effects of various parameters on the pressure oscillation. Four evaporators with different wick properties were tested using a test loop containing two condenser plates. The test loop allowed the four evaporators to be tested in a single-pump, two-pump or four-pump configuration, and the two condenser plates to be plumbed either in parallel or in series. Test conditions included varying the power input, the reservoir set point temperature, the condenser sink temperature, and the flow resistance between the reservoir and the loop. Experimental results agreed well with theoretical predictions.

Description: The Liquid Motion Experiment (LME), designed to study the effects of liquid motion in rotating tanks, was flown on STS 84. LME was essentially a spin table that created a realistic nutation motion of scale-model tanks containing liquid. TWo spherical and two cylindrical transparent tanks were tested simultaneously, and three sets of such tanks were employed to vary liquid viscosity, fill level, and propellant management device (PMD) design. All the tanks were approximately 4.5 inches diameter. The primary test measurements were the radial and tangential torques exerted on the tanks by the liquid. Resonant frequencies and damping of the liquid oscillations were determined by sine sweep tests. For a given tank shape, the resonant frequency depended on fill level. For the cylindrical tanks, the resonances had somewhat different frequencies for the tangential axis (0.55 to 0.75 times spin rate) and the radial axis (0.73 to 0.78 times spin rate), and the tangential axis resonance agreed more closely with available analytical models. For the spherical tanks, the resonant frequencies were between 0.74 to 0.77 times the spin rate and were the same for the tangential and radial axes. The damping coefficients varied from about I% to 3% of critical, depending on tank shape, fill level, and liquid viscosity. 'Me viscous energy dissipation rates of the liquid oscillations were determined from sine dwell tests. The LME energy dissipation rates varied from 0.3 to 0.5 times the estimates obtained from scaling previous ground tests and spacecraft flight data. The PNDs sometimes enhanced the resonances and energy dissipation rates and sometimes decreased them, which points out the need to understand better the effects of PMD on liquid motion as a function of PMD and tank design.

Description: A tab placed suitably on a nozzle that produces a jet in a cross-flow can reduce the penetration of the jet. This effect, achieved when the tab is placed on the windward side of the nozzle relative to the cross flow, may be of interest in film cooling applications. Wind tunnel experiments are carried out, in the momentum ratio (J) range of 10-90, to investigate the tab geometry that would maximize this effect. The preliminary results show that a 'delta tab' having a base width approximately fifty percent of the nozzle diameter may be considered optimum. With a given tab size, the effect is more pronounced at higher J. Reduction in jet penetration by as much as 40% is observed. Comparable reduction in jet penetration is also obtained when a triangular shaped tab is placed flush with the tunnel wall or with its apex tilted down into the jet nozzle (the 'delta tab' being the configuration in which the apex is tilted up). However, the delta tab involves the least flow blockage and pressure loss. Relative to the baseline case, the lateral spreading of the jet is found to be more with the delta tab but less with other orientations of the tab.

Description: A theoretical explanation of some experimentally observed phenomena associated with the so-called Klebanoff modes is obtained by analyzing the flow over a finite thickness flat plate resulting from a small-amplitude distortion imposed on the upstream mean flow. The analysis shows (among other things) how the stretching of the vortex lines around the plate leads to streamwise vorticity at the plate surface, which then produces a streamwise velocity perturbation within the boundary layer that can be related to the experimentally observed Klebanoff mode. The complete evolution of this flow must be found by solving the boundary-region equations of Kemp (1951) and Davis and Rubin (1980), but a limiting analytical solution can also be obtained. Since the initial growth of the boundary-layer disturbance is nearly algebraic, our results demonstrate how the algebraically growing disturbances promoted by Landahl and others can be generated by a realistic external-disturbance environment. The relationship between these results and various bypass transition mechanisms is discussed.

Description: Wave packets produced by localized disturbances play an important role in transition in three-dimensional boundary layers, such as that on a swept wing. Starting with the receptivity process, we show the effects of wave-space energy distribution on the development of packets and other three-dimensional disturbance patterns. Nonlinearity in the receptivity process is specifically addressed, including demonstration of an effect which can enhance receptivity of traveling crossflow disturbances. An efficient spatial numerical simulation method is allowing most of the simulations presented to be carried out on a workstation.

Description: The objective of this study was to analyze jet plume mass flow entrainment rates associated with the introduction of counter-rotating streamwise vorticity by prism shaped devices (tabs) located at the lip of the nozzle. We have examined the resulting mixing process through coordinated experimental tests and numerical simulations of the supersonic flow from a model axisymmetric nozzle. In the numerical simulations, the total induced vorticity was held constant while varying the distribution of counter-rotating vorticity around the nozzle lip training edge. In the experiment, the number of tabs applied was varied while holding the total projected area constant. Evaluations were also conducted on initial vortex strength. The results of this work show that the initial growth rate of the jet shear layer is increasingly enhanced as more tabs are added, but that the lowest tab count results in the largest entrained mass flow. The numerical simulations confirm these results.

Description: The main goal of this project was efficient distributed parallel and workstation cluster implementations of Newton-Krylov-Schwarz (NKS) solvers for implicit Computational Fluid Dynamics (CFD.) "Newton" refers to a quadratically convergent nonlinear iteration using gradient information based on the true residual, "Krylov" to an inner linear iteration that accesses the Jacobian matrix only through highly parallelizable sparse matrix-vector products, and "Schwarz" to a domain decomposition form of preconditioning the inner Krylov iterations with primarily neighbor-only exchange of data between the processors. Prior experience has established that Newton-Krylov methods are competitive solvers in the CFD context and that Krylov-Schwarz methods port well to distributed memory computers. The combination of the techniques into Newton-Krylov-Schwarz was implemented on 2D and 3D unstructured Euler codes on the parallel testbeds that used to be at LaRC and on several other parallel computers operated by other agencies or made available by the vendors. Early implementations were made directly in Massively Parallel Integration (MPI) with parallel solvers we adapted from legacy NASA codes and enhanced for full NKS functionality. Later implementations were made in the framework of the PETSC library from Argonne National Laboratory, which now includes pseudo-transient continuation Newton-Krylov-Schwarz solver capability (as a result of demands we made upon PETSC during our early porting experiences). A secondary project pursued with funding from this contract was parallel implicit solvers in acoustics, specifically in the Helmholtz formulation. A 2D acoustic inverse problem has been solved in parallel within the PETSC framework.

Description: We have developed a code where the nonlinear terms are treated implicitly. The equations are discretized using the two-point fourth order compact scheme in the y-direction and the backward Euler method in the x-direction. We investigated the transition process in a Blasius boundary layer due to fundamental type breakdown. With 8 modes in the w and 3 planes, we could compute the evolution of disturbances up to Re(x)=910, which is well into the strongly nonlinear region. The transition onset point is located around Re(x)=850. The comparison with the measurements and with the DNS computations are very good up to Re(x)=880.

Description: The Strutjet approach to Rocket Based Combined Cycle (RBCC) propulsion depends upon fuel-rich flows from the rocket nozzles and turbine exhaust products mixing with the ingested air for successful operation in the ramjet and scramjet modes. It is desirable to delay this mixing process in the air-augmented mode of operation present during low speed flight. A model of the Strutjet device has been built and is undergoing test to investigate the mixing of the streams as a function of distance from the Strutjet exit plane during simulate low speed flight conditions. Cold flow testing of a 1/6 scale Strutjet model is underway and nearing completion. Planar Laser Induced Fluorescence (PLIF) diagnostic methods are being employed to observe the mixing of the turbine exhaust gas with the gases from both the primary rockets and the ingested air simulating low speed, air augmented operation of the RBCC. These diagnostics are supplemented by pressure profile measurements across the various jet streams present in the flow duct downstream of the model exit plane. The results of these experiments will be presented.

Description: This paper studies concepts for the Hydrogen On-Orbit Storage and Supply Experiment (HOSS). HOSS is a space flight experiment whose objectives are: Show stable gas supply for storage and direct gain solar-thermal thruster designs; and evaluate and compare low-g performance of active and passive pressure control via a thermodynamic vent system (TVS) suitable for solar-thermal upper stages. This paper shows that the necessary experimental equipment for HOSS can be accommodated in a small hydrogen dewar of 36 to 80 liter. Thermal designs for these dewars which meet the on-orbit storage requirements can be achieved. Furthermore ground hold insulation and shielding concepts are achieved which enable storing initially subcooled liquid hydrogen in these small dewars without venting in excess of 144 hours.

Description: The objective of this project is to conduct measurements of elementary reaction rate constants and photochemical parameters for processes of importance in the atmosphere. These measurements are being carried out under temperature and pressure conditions covering those applicable to the stratosphere and upper troposphere, using the chemical ionization mass spectrometry turbulent flow technique developed in our laboratory. The next section summarizes our research activities during the first year of the project, and the section that follows consists of the statement of work for the third year. Additional details concerning the projects listed in the statement of work were described in our original proposal.

Description: The non-equilibrium critical-layer analysis of a system of frequency-detuned resonant-triads is presented using the generalized scaling of Lee. It is shown that resonant-triads can interact nonlinearly within the common critical layer when their (fundamental) Strouhal numbers are different by a factor whose magnitude is of the order of the growth rate multiplied by the wavenumber of the instability wave. Since the growth rates of the instability modes become larger and the critical layers become thicker as the instability waves propagate downstream, the frequency-detuned resonant-triads that grow independently of each other in the upstream region can interact nonlinearly in the later downstream stage. In the final stage of the non-equilibrium critical-layer evolution, a wide range of instability waves with the scaled frequencies differing by almost an Order of (l) can nonlinearly interact. Low-frequency modes are also generated by the nonlinear interaction between oblique waves in the critical layer. The system of partial differential critical-layer equations along with the jump equations are presented here. The amplitude equations with their numerical solutions are given in Part 2. The nonlinearly generated low-frequency components are also investigated in Part 2.

Description: An approach which closely maintains the non-dissipative nature of classical fourth or higher- order spatial differencing away from shock waves and steep gradient regions while being capable of accurately capturing discontinuities, steep gradient and fine scale turbulent structures in a stable and efficient manner is described. The approach is a generalization of the method of Gustafsson and Olsson and the artificial compression method (ACM) of Harten. Spatially non-dissipative fourth or higher-order compact and non-compact spatial differencings are used as the base schemes. Instead of applying a scalar filter as in Gustafsson and Olsson, an ACM like term is used to signal the appropriate amount of second or third-order TVD or ENO types of characteristic based numerical dissipation. This term acts as a characteristic filter to minimize numerical dissipation for the overall scheme. For time-accurate computations, time discretizations with low dissipation are used. Numerical experiments on 2-D vortical flows, vortex-shock interactions and compressible spatially and temporally evolving mixing layers showed that the proposed schemes have the desired property with only a 10% increase in operations count over standard second-order TVD schemes. Aside from the ability to accurately capture shock-turbulence interaction flows, this approach is also capable of accurately preserving vortex convection. Higher accuracy is achieved with fewer grid points when compared to that of standard second-order TVD or ENO schemes. To demonstrate the applicability of these schemes in sustaining turbulence where shock waves are absent, a simulation of 3-D compressible turbulent channel flow in a small domain is conducted.

Description: The objective of the research was to develop a capability in the NPARC computational fluid dynamics (CFD) code to efficiently solve for unsteady airflows with moving geometry and grids. The application of interest was the unsteady flow in a high-speed aircraft inlet operating at the supercritical condition in which a terminal shock resides within the diffuser.

Description: A technique is implemented for computing hypersonic aeroheating, shear stress, and other flow properties on the windward side of a three-dimensional (3D) blunt body. The technique uses a 2D/axisymmetric flow solver modified by scale factors for a, corresponding equivalent axisymmetric body. Examples are given in which a 2D solver is used to calculate the flow at selected meridional planes on elliptic paraboloids in reentry flight. The report describes the equations and the codes used to convert the body surface parameters into input used to scale the 2D viscous shock layer equations in the axisymmetric viscous shock layer code. Very good agreement is obtained with solutions to finite rate chemistry 3D thin viscous shock layer equations for a finite rate catalytic body.

Description: Direct Numerical Simulation (DNS) is used to study the mechanism of generation and evolution of turbulence structures in a temporally evolving supersonic swirling round jet and also to examine the resulting acoustic radiations. Fourier spectral expansions are used in the streamwise and azimuthal directions and a 1-D b-spline Galerkin representation is used in the radial direction. Spectral-like accuracy is achieved using this numerical scheme. Direct numerical simulations, using the b-spline spectral method, are carried out starting from mean flow initial conditions which are perturbed by the most unstable linear stability eigenfunctions. It is observed that the initial helical instability waves evolve into helical vortices which eventually breakdown into smaller scales of turbulence. 'Rib' structures similar to those seen in incompressible mixing layer flow of Rogers and Moserl are observed. The jet core breakdown stage exhibits increased acoustic radiations.

Description: This is a secondary analysis of a portion of the Halstead data. The hot-film traces from an embedded stage of a low pressure turbine have been extensively analyzed by Halstead et al. In this project, wavelet analysis is used to develop the quantitative characterization of individual coherent structures in terms of size, amplitude, phase, convection speed, etc., as well as phase-averaged time scales. The purposes of the study are (1) to extract information about turbulent time scales for comparison with unsteady model results (e.g. k/epsilon). Phase-averaged maps of dominant time scales will be presented; and (2) to evaluate any differences between wake-induced and natural spots that might affect model performance. Preliminary results, subject to verification with data at higher frequency resolution, indicate that spot properties are independent of their phase relative to the wake footprints: therefore requirements for the physical content of models are kept relatively simple. Incidentally, we also observed that spot substructures can be traced over several stations; further study will examine their possible impact.

Description: The objective of the research reported in this presentation is to develop computational techniques for the prediction of unsteady transitional flows associated with the rotor stator interaction in turbomachinery. Three low-Reynolds number turbulence models are incorporated in two unsteady Navier-Stokes codes (one is pressure based and the other is time marching with Runge-Kutta time stepping) and evaluated for accuracy in predicting the onset and the end of unsteady transitional patches due to wake passing. The best model is then used for modification and improvement for the leading edge effect. An existing steady Navier-Stokes code was modified to include pseudo-time stepping, which provided acceleration from 5 to 25 times that of the original code. A systematic validation procedure was implemented to assess the effects of the grid, artificial dissipation, physical, and the pseudo-time step for an accurate prediction of transitional flows resulting from the rotor-stator interaction. The ability of the Navier-Stokes code to predict the unsteady transitional flow on a turbomachinery blade is demonstrated. The unsteady pressure and velocity fields are in good agreement with the experimental data and the prediction from the Euler/boundary layer approach. The numerical solver was able to capture all zones (wake induced transitional strip, wake induced turbulent strip, calmed region, etc.) associated with wake induced transition in a compressor cascade. Another significant step is the assessment of k-epsilon turbulence models, including the leading edge modifications. Best results were obtained from the FLB model. The LB model predicted earlier inception of the transition and shorter transition length. Modification of the k-epsilon model was found to be essential for an accurate prediction of the unsteady transitional flow in a compressor cascade. The CH model failed to predict the unsteady transitional flow. Predicted boundary layer was turbulent from the leading edge, even with the modification of the k-epsilon model near the stagnation point. A comparison of the instantaneous shape factor, the skin friction coefficient, and the momentum thickness indicates that the Navier-Stokes predictions are reasonably good. Interaction between the upstream wake and stator wake results in shedding of unsteady vortices from the trailing edge and increased dissipation in the stator wake and, as a consequence, increased rate of decay of the stator wake. The procedure developed at Penn State should aid the designers of turbomachinery in allowing for the unsteady transitional flows due to rotor-stator interaction.

Description: We have further developed our capabilities to analyze transition in turbine boundary layers from first principles by integrating the nonlinear parabolized stability equations (PSE) with improved initial and boundary conditions. With modified iteration schemes, we are able to proceed deeper into the transition region where skin friction coefficient and heat transfer coefficient significantly increase. Initial and boundary conditions at elevated turbulence levels can be derived by receptivity analysis. Test runs for ERCOFTAC test case T3A at 2.4\% turbulence level provide results in good agreement with the experimental data. The sharper minimum of the skin coefficient also shown by DNS results is likely due to the missing intermittency. The method has been applied to various experimentally studied turbine blades (UTRC, VKI, Zierke, Langston, Hippensteele, and others). The PSE results, though physically reasonable, do not agree as well as expected with the experimental findings. We have, therefore, performed an extensive search for the reasons of the seemingly systematic deviations. A first source of uncertainty has been found in the often insufficient documentation of the experiments (e.g. on blockage by end-wall boundary layers). However, variation of the relevant parameters does not lead to more satisfactory agreement. A second reason has been found in the "standard procedure" which considers a 2D flow at midspan and uses a panel code and subsequent boundary-layer code to obtain the laminar basic flow for the transition analysis. Comparison with the pressure distribution obtained with a 3D design code (RVC3D) shows significant three-dimensionality of the flow (e.g. in the UTRC experiments). The spanwise variation has been neglected in our original PSE code. To overcome this problem, we have developed the PSE/3D for fully 3D boundary layers to account for streamwise and spanwise variations. Since the design code does not provide the boundary-layer flow with sufficient resolution, we have generated the Euler solution and employed a 3D boundary-layer code to obtain the viscous basic flow. Although only the linear stability level of PSE/3D has been implemented so far, the discrepancies with the experiments change but do not disappear. We still find deviations between the computed and experimental variations of C(sub f), and St along the blade for laminar flow. The main reason can be seen by comparing the solution of the boundary-layer code with the viscous results of the design code. The conventional boundary-layer solution exhibits an asymptotic behavior appropriate in external aerodynamics but does not match the steep gradients of the inviscid flow through the passage and consequently provides biased results for C(sub f), and St. An attempt is currently being made to correct this deficiency. Before attempting to perform the transition analysis for the viscous flow provided by the design code, we have analyzed the implementation and "best possible" results. Code and results exhibit flaws that may negatively affect the design and are intolerable for transition analysis. Therefore, we have decided to develop a new code to obtain a reliable basis for stability and transition studies. We expect to report improved results by the time of the meeting.

Description: Current turbulence models tend to give too early and too short a length of flow transition to turbulence, and hence fail to predict flow separation induced by the adverse pressure gradients and streamline flow curvatures. Our discussion will focus on the development and validation of transition models. The baseline data for model comparisons are the T3 series, which include a range of free-stream turbulence intensity and cover zero-pressure gradient to aft-loaded turbine pressure gradient flows. The method will be based on the conditioned N-S equations and a transport equation for the intermittency factor. First, several of the most popular 2-equation models in predicting flow transition are examined: k-e [Launder-Sharina], k-w [Wilcox], Lien-Leschiziner and SST [Menter] models. All models fail to predict the onset and the length of transition, even for the simplest flat plate with zero-pressure gradient(T3A). Although the predicted onset position of transition can be varied by providing different inlet turbulent energy dissipation rates, the appropriate inlet conditions for turbulence quantities should be adjusted to match the decay of the free-stream turbulence. Arguably, one may adjust the low-Reynolds-number part of the model to predict transition. This approach has so far not been very successful. However, we have found that the low-Reynolds-number model of Launder and Sharma [1974], which is an improved version of Jones and Launder [1972] gave the best overall performance. The Launder and Sharma model was designed to capture flow re-laminarization (a reverse of flow transition), but tends to give rise to a too early and too fast transition in comparison with the physical transition. The three test cases were for flows with zero pressure gradient but with different free-stream turbulent intensities. The same can be said about the model when considering flows subject to pressure gradient(T3C1). To capture the effects of transition using existing turbulence models, one approach is to make use of the concept of the intermittency to predict the flow transition. It was originally based on the intermittency distribution of Narasimha [1957], and then gradually evolved into a transport equation for the intermittency factor. Gostelow and associates [1994,1995] have made some improvements to Narasimha's method in an attempt to account for both favorable and adverse pressure gradients. Their approach is based on a linear, explicit combination of laminar and turbulent solutions. This approach fails to predict the overshoot of the skin friction on a flat plate near the end of transition zone, even though the length of transition is well predicted. The major flaw of Gostelow's approach is that it assumes the non-turbulent part being the laminar solution and the turbulent part being the turbulent solution and they do not interact across the transitional region. The technique in condition averaging the flow equations in intermittent flows was first introduced by Libby [1975] and Dopazo [1977] and further refined by Dick and associates [1988, 1996]. This approach employs two sets of transport equations for the non-turbulent part and the other for the turbulent part. The advantage of this approach is that it allows the interaction of non-turbulent and turbulent velocities through the introduction of additional source terms in the continuity and momentum equations for the non-turbulent and turbulent velocities. However, the strong coupling of the two sets of equations has caused some numerical difficulties, which requires special attention. The prediction of the skin friction can be improved by this approach via the implicit coupling of non-turbulent and turbulent velocity flelds. Another improvement of the interrmittency model can be further made by allowing the intermittency to vary in the cross-stream direction. This is one step prior to testing any proposal for the transport equation for the intermittency factor. Instead of solving the transport equation for the intermittency factor, the distribution for the intermittency factor is prescribed by Klebanoff's empirical formula [1955]. The skin friction is very well predicted by this new modification, including the overshoot of the profile near the end of the transition zone. The outcome of this study is very encouraging since it indicates that the proper description of the intermittency distribution is the key to the success of the model prediction. This study will be used to guide us on the modelling of the intermittency transport equation.

Description: Transition plays an important role in the prediction of losses and performance in low-pressure turbines. The transition location on a turbine blade may vary significantly because of the wakes from upstream blade rows, and intermittent flow separation can also affect the transition process in an unsteady flow environment. In the present investigation, an unsteady Navier-Stokes analysis is used to predict transition in a low-pressure turbine cascade and a low-pressure turbine stage. The numerical flow analysis is third-order spatially accurate and second-order temporally accurate, and the equations of motion are integrated using an implicit time-marching procedure. The Baldwin-Lomax and k-epsilon turbulence models, in conjunction with several algebraic transition models, have been used to predict the location of transition. Predicted results include unsteady blade loadings, time-histories of the pressure, transition locations and boundary layer quantities, as well as performance quantities and comparisons with the available experimental/design data.

Description: Results of an experimental study of flow separation and transition in either attached boundary layers or separated shear layers over the suction surface of a simulation of a low-pressure turbine airfoil flow are presented. Detailed velocity profiles were measured with the hot-wire technique. Static pressure distributions are also presented. Flow transition is documented using measured intermittency distributions in the boundary layer and the separated shear layer. Cases for Reynolds numbers of 50,000, 100,000, 200,000 and 300,000 are reported. These Reynolds numbers are based on suction surface length and exit velocity. Three Free Stream Turbulence Intensity values, 0.5%, 2.5% and 10%, are represented. Flow separation is observed for all the low-FSTI cases. Of these, the lowest Reynolds number case was not able to complete transition of the shear layer and the separation bubble persisted over the entire blade surface. For the other low-FSTI cases, transition is observed in the shear layer over the separation bubble. This transition proceeded quickly, spreading rapidly toward the wall. Elevated FSTI drives an earlier transition than in the low-FSTI cases and the separation bubbles are smaller. For the highest Reynolds number cases with 2.5% and 10% FSTI, transition is of the attached boundary layer and no separation exists. Flow separation with shear flow transition is observed for the lower-Re cases. Models for intermittency and transition length and location from the modern literature are assessed.

Description: When a shear-sensitive liquid crystal coating is illuminated from the normal direction by white light and observed from an oblique above-plane view angle, its color-change response to shear depends on both shear stress vector magnitude and the direction of the applied shear vector relative to the observer's in-plane line of sight. At any point, the maximum color change is always seen or measured when the local shear vector is aligned with, and directed away from, the observer; the magnitude of the color change at this vector/observer aligned orientation scales directly with shear stress magnitude. Conversely, any point exposed to a shear vector with a component directed toward the observer exhibits a noncolor-change response, always characterized by a rusty red or brown color, independent of both shear magnitude and direction. Based on this knowledge, full-surface shear stress vector visualization and measurement methodologies were formulated and successfully demonstrated. The present paper reviews the observations and measurements that led to the development of these methodologies and applications of both are discussed.

Description: Lie group methods are used to find both exact and numerical similarity solutions for compressible perturbations to all incompressible, two-dimensional, axisymmetric vortex reference flow. The reference flow vorticity satisfies an eigenvalue problem for which the solutions are a set of two-dimensional, self-similar, incompressible vortices. These solutions are augmented by deriving a conserved quantity for each eigenvalue, and identifying a Lie group which leaves the reference flow equations invariant. The partial differential equations governing the compressible perturbations to these reference flows are also invariant under the action of the same group. The similarity variables found with this group are used to determine the decay rates of the velocities and thermodynamic variables in the self-similar flows, and to reduce the governing partial differential equations to a set of ordinary differential equations. The ODE's are solved analytically and numerically for a Taylor vortex reference flow, and numerically for an Oseen vortex reference flow. The solutions are used to examine the dependencies of the temperature, density, entropy, dissipation and radial velocity on the Prandtl number. Also, experimental data on compressible free vortex flow are compared to the analytical results, the evolution of vortices from initial states which are not self-similar is discussed, and the energy transfer in a slightly-compressible vortex is considered.

Description: Results are reported from the Second Surface Tension Driven Convection Experiment (STDCE-2) performed aboard the Second United States Microgravity Laboratory (USML-2), which flew as Space Shuttle mission STS-73 and launched on October 20, 1995. Oscillatory thermocapillary flows were investigated in open cylindrical containers filled with 2 centistokes kinematic viscosity (Pr=27 at 25 C) silicone oil. Two different heating modes were investigated in detail. The main objectives of the experiments were to determine the onset of oscillatory thermocapillary flow under highly reduced buoyancy and gravity conditions and to study the important features of the oscillatory flow. The onset conditions were determined in three different size containers (1.2-, 2-, and 3-cm diameter) and for various free surface shapes (flat and curved). Numerical and scaling analyses were also performed to understand the basic steady flows. The analysis shows that the main flow is viscous-dominated near the onset of oscillations. The onset conditions determined in the present experiments together with our earlier ground-based data show that one could not describe the oscillation phenomenon if the fluid free surface is assumed to be rigid. Therefore, a parameter representing free surface deformation was derived, and the data are shown to be correlated well by that parameter. The oscillation patterns and frequencies are also presented.

Description: The effect of a Mach number correction on a model for predicting the length of transition was investigated. The transition length decreases as the turbulent spot production rate increases. Much of the data for predicting the spot production rate comes from low speed flow experiments. Recent data and analysis showed that the spot production rate is affected by Mach number. The degree of agreement between analysis and data for turbine blade heat transfer without film cooling is strongly dependent of accurately predicting the length of transition. Consequently, turbine blade heat transfer data sets were used to validate a transition length turbulence model. A method for modifying models for the length of transition to account for Mach number effects is presented. The modification was made to two transition length models. The modified models were incorporated into the two-dimensional Navier-Stokes code, RVCQ3D. Comparisons were made between predicted and measured midspan surface heat transfer for stator and rotor turbine blades. The results showed that accounting for Mach number effects significantly improved the agreement with the experimental data.

Description: A flat-plate heat pipe (FPHP) experiment has been set up for micro-gravity tests on a NASA supplied aircraft. This report presents an analysis on various components of the experimental setup to certify that it will satisfy the flight safety and operation requirements.

Description: Oscillatory thermocapillary flow experiments were performed in the Glovebox aboard the USML-2 Spacelab which was launched on October 20, 1995. Open cylindrical containers of 1.2 and 2.0 cm in diameter were used. The ratio of container depth to radius (aspect ratio) was set at 0.5 and 2. Silicone oil of 2 centistokes kinematic viscosity was the test fluid. The fluid was heated by a cylindrical heater placed along the center axis of the container. The fluid motion was studied by flow visualization. The effect of aspect ratio on the onset of oscillations and on the oscillatory flow was investigated. It was found that the onset of oscillations was delayed when the container was made more shallow.

Description: Numerical modeling of the nonreactive mixing processes associated with a lobed fuel injector in a coflowing air stream is presented. The lobed fuel injector is a device which generates strong streamwise vorticity, producing locally high strain rates which can enhance the molecular mixing of reactants while delaying ignition in a controlled manner. Vortex element modeling is used to simulate flow field evolution and fuel element mixing characteristics for this lobed fuel injector. Quantitative predictions for vorticity generation and qualitative results for streamwise rollup compare well qualitatively with recent experimental investigations of this flow field [Smith et al, Phys. Fluids 9, 667 (1997)]. Parametric studies of the effects of lobe amplitude-to-wavelength ratio, lobe angle, and lobe shape for given flow conditions suggest that geometrical features may be optimized to enhance mixing and control reaction processes.

Description: Under this grant, the flutter versions of the three-dimensional Navier-Stokes aeroelastic code (TURBO-AE) have been developed. The TURBO-AE aeroelastic code allows flutter calculations for a fan, compressor, or turbine blade row. This code models a vibrating three-dimensional bladed disk geometry and the associated unsteady flow (including shocks, and viscous effects) to calculate the aeroelastic instability using a work-per-cycle or eigenvalue approach. Steady and unsteady calculations have been performed with the 3D Navier-Stokes (TURBO) code to gain familiarity with the code. Calculations were done with the NASA Rotor-67 configuration which has been used as a standard test case by turbomachinery researchers. Additional calculations have been performed for the NASA/GE EEE fan and a proprietary configuration.

Description: A quasitransient in-ground effect method is developed to study the effect of vertical landing on a launch vehicle base-heating environment. This computational methodology is based on a three-dimensional, pressure-based, viscous flow, chemically reacting, computational fluid dynamics formulation. Important in-ground base-flow physics such as the fountain-jet formation, plume growth, air entrainment, and plume afterburning are captured with the present methodology. Convective and radiative base-heat fluxes are computed for comparison with those of a flight test. The influence of the laminar Prandtl number on the convective heat flux is included in this study. A radiative direction-dependency test is conducted using both the discrete ordinate and finite volume methods. Treatment of the plume afterburning is found to be very important for accurate prediction of the base-heat fluxes. Convective and radiative base-heat fluxes predicted by the model using a finite rate chemistry option compared reasonably well with flight-test data.

Description: A Direct Numerical Simulation of a coolant jet injected normally into a periodic crossflow is presented. The physical situation simulated represents a periodic module in a coolant hole array with a heated crossflow. A collocated finite difference scheme is used which is fifth-order accurate spatially and second-order accurate temporally. The scheme is based on a fractional step approach and requires the solution of a pressure-Poisson equation. The simulations are obtained for a blowing ratio of 0.25 and a channel Reynolds number of 5600. The simulations reveal the dynamics of several large scale structures including the Counter-rotating Vortex Pair (CVP), the horse-shoe vortex, the shear layer vortex, the wall vortex and the wake vortex. The origins and the interactions of these vortical structures are identified and explored. Also presented are the turbulence statistics and how they relate to the flow structures.

Description: The effects of rotation on turbulent, compressible convection within stellar envelopes are studied through three-dimensional numerical simulations conducted within a local f-plane model. This work seeks to understand the types of differential rotation that can be established in convective envelopes of stars like the Sun, for which recent helioseismic observations suggest an angular velocity profile with depth and latitude at variance with many theoretical predictions. This paper analyzes the mechanisms that are responsible for the mean (horizontally averaged) zonal and meridional flows that are produced by convection influenced by Coriolis forces. The compressible convection is considered for a range of Rayleigh, Taylor, and Prandtl (and thus Rossby) numbers encompassing both laminar and turbulent flow conditions under weak and strong rotational constraints. When the nonlinearities are moderate, the effects of rotation on the resulting laminar cellular convection leads to distinctive tilts of the cell boundaries away from the vertical. These yield correlations between vertical and horizontal motions that generate Reynolds stresses that can drive mean flows, interpretable as differential rotation and meridional circulations. Under more vigorous forcing, the resulting turbulent convection involves complicated and contorted fluid particle trajectories, with few clear correlations between vertical and horizontal motions, punctuated by an evolving and intricate downflow network that can extend over much of the depth of the layer. Within such networks are some coherent structures of vortical downflow that tend to align with the rotation axis. These yield a novel turbulent alignment mechanism, distinct from the laminar tilting of cellular boundaries, that can provide the principal correlated motions and thus Reynolds stresses and subsequently mean flows. The emergence of such coherent structures that can persist amidst more random motions is a characteristic of turbulence with symmetries broken by rotation and stratification. Such structure is here found to play a crucial role in defining the mean zonal and meridional flows that coexist with the convection. Though they are subject to strong inertial oscillations, the strength and type of the mean flows are determined by a combination of the laminar tilting and the turbulent alignment mechanisms. Varying the parameters produces a wide range of mean motions. Among these, some turbulent solutions exhibit a mean zonal velocity profile that is nearly constant with depth, much as deduced by helioseismology at midlatitudes within the Sun. The solutions exhibit a definite handedness, with the direction of the persistent mean flows often prescribing a spiral with depth near the boundaries, also in accord with helioseismic deductions. The mean helicity has a profile that is positive in the upper portion of the domain and negative in the lower portion, a property bearing on magnetic dynamo processes that may be realized within such rotating layers of turbulent convection.

Description: There is a continuing need to understand the fluid physics occurring under low gravity conditions in processes such as crystal growth, materials processing, and the movement of bubbles or droplets. The fluid flow in such situations is often caused by a gradient in interfacial tension. If a temperature gradient is created due to a heat source, the resulting flow is called thermocapillary flow, a special case of Marangoni Convection. In this study, an experimental investigation was conducted using silicone oil in cylindrical containers with a laser heat source at the free surface. It was desired to determine the conditions under which steady, axisymmetrical thermocapillary flow becomes unstable and oscillatory three-dimensional flow states develop. The critical Marangoni number for each observed oscillatory state was measured as a function of the container aspect ratio and the dynamic Bond number, a measure of buoyant force versus ii thermocapillary force. Various oscillatory modes were observed during three- dimensional convection, and chaotic flow was reached in one test condition. The critical Marangoni numbers are compared with those measured in previous studies, and the power spectra and phase trajectories of the instantaneous surface temperature distributions are used to characterize the routes of transitions to the chaotic flow state. Results show that only superharmonic modes appear in the routes to chaos while infinite number of subharmonic modes occur in flow transitions for pure Rayleigh convection.

Description: This paper describes a general purpose computer program for analyzing steady state and transient flow in a complex network. The program is capable of modeling phase changes, compressibility, mixture thermodynamics and external body forces such as gravity and centrifugal. The program's preprocessor allows the user to interactively develop a fluid network simulation consisting of nodes and branches. Mass, energy and specie conservation equations are solved at the nodes; the momentum conservation equations are solved in the branches. The program contains subroutines for computing "real fluid" thermodynamic and thermophysical properties for 33 fluids. The fluids are: helium, methane, neon, nitrogen, carbon monoxide, oxygen, argon, carbon dioxide, fluorine, hydrogen, parahydrogen, water, kerosene (RP-1), isobutane, butane, deuterium, ethane, ethylene, hydrogen sulfide, krypton, propane, xenon, R-11, R-12, R-22, R-32, R-123, R-124, R-125, R-134A, R-152A, nitrogen trifluoride and ammonia. The program also provides the options of using any incompressible fluid with constant density and viscosity or ideal gas. Seventeen different resistance/source options are provided for modeling momentum sources or sinks in the branches. These options include: pipe flow, flow through a restriction, non-circular duct, pipe flow with entrance and/or exit losses, thin sharp orifice, thick orifice, square edge reduction, square edge expansion, rotating annular duct, rotating radial duct, labyrinth seal, parallel plates, common fittings and valves, pump characteristics, pump power, valve with a given loss coefficient, and a Joule-Thompson device. The system of equations describing the fluid network is solved by a hybrid numerical method that is a combination of the Newton-Raphson and successive substitution methods. This paper also illustrates the application and verification of the code by comparison with Hardy Cross method for steady state flow and analytical solution for unsteady flow.

Description: The results of the phase 2 effort conducted under AGARD Working Group 18 on Hypersonic Experimental and Computational Capability, Improvement and Validation are presented in this report. The first volume, published in May 1996, mainly focused on the design methodology, plans and some initial results of experiments that had been conducted to serve as validation benchmarks. The current volume presents the detailed experimental and computational data base developed during this effort.

Description: Dryout occurs readily in certain slug and annular two-phase flows for systems that exhibit partial wetting. The mechanism for the ultimate rupture of the film is attributed to van der Waals forces, but the pace towards rupture is quickened by the surface tension instability (Rayleigh-type) of the annular film left by the advancing slug and by the many perturbations of the free surface present in the Re(sub g) approximately 0(10(exp 3)), Re(sub l) approximately 0(10(exp 4)), and Ca approximately 0(10(exp -1) flows. Results from low-gravity experiments using three different test fluids are presented and discussed. For the range of tests conducted, the effect of increasing viscosity is shown to eliminate the film rupture while the decrease of surface tension via a surfactant additive is shown to dramatically enhance it. Laboratory measurements using capillary tubes are presented which reveal the sensitivity of the dryout phenomena to particulate and surfactant contamination. Rom such observations, dryout due to the hydrodynamic-van der Waals instability can be expected in a certain range of flow parameters in the absence of heat transfer. The addition of heat transfer may only exacerbate the problem by producing thermal transport lines replete with "hot spots." A caution to this effect is issued to future space systems designers concerning the use of partially wetting working fluids.

Description: Propulsion performance and operability are key factors in the development of a successful aircraft. For high-speed supersonic aircraft, mixed-compression inlets offer high performance but are susceptible to an instability referred to as unstart. An unstart occurs when a disturbance originating in the atmosphere or the engine causes the shock system to be expelled from the inlet. This event can have adverse effects on control of the aircraft, which is unacceptable for a passenger plane such as the high speed civil transport (HSCT). The ability to predict the transient response of such inlets to flow perturbations is, therefore, important to the proper design of the inlet and the control measures used to prevent unstart. Computational fluid dynamics (CFD) is having an increasing role in the analysis of individual propulsion components. isolated inlet studies are relatively easy to perform, but a major uncertainty is the boundary condition used at the inlet exit to represent the engine - the so-called compressor face boundary condition. A one-dimensional (I-D) Euler inlet simulation (ref. 1) showed that the predicted inlet unstart tolerance to free-stream pressure perturbations can vary by as much as a factor of about six, depending on the boundary condition used. Obviously, a thorough understanding of dynamic interactions between inlets and compressors/fans is required to provide the proper boundary condition.

Description: Abstract Accurate and efficient direct numerical simulation of turbulence in the presence of shock waves represents a significant challenge for numerical methods. The objective of this paper is to evaluate the performance of high order compact and non-compact central spatial differencing employing total variation diminishing (TVD) shock-capturing dissipations as characteristic based filters for two model problems combining shock wave and shear layer phenomena. A vortex pairing model evaluates the ability of the schemes to cope with shear layer instability and eddy shock waves, while a shock wave impingement on a spatially-evolving mixing layer model studies the accuracy of computation of vortices passing through a sequence of shock and expansion waves. A drastic increase in accuracy is observed if a suitable artificial compression formulation is applied to the TVD dissipations. With this modification to the filter step the fourth-order non-compact scheme shows improved results in comparison to second-order methods, while retaining the good shock resolution of the basic TVD scheme. For this characteristic based filter approach, however, the benefits of compact schemes or schemes with higher than fourth order are not sufficient to justify the higher complexity near the boundary and/or the additional computational cost.

Description: A videotape presentation of flow in a packed bed of spheres is provided. The flow experiment consisted of three principal elements: (1) an oil tunnel 76.2 mm by 76.2 mm in cross section, (2) a packed bed of spheres in regular and irregular arrays, and (3) a flow characterization methodology, either (a) full flow field tracking (FFFT) or (b) computational fluid dynamic (CFD) simulation. The refraction indices of the oil and the test array of spheres were closely matched, and the flow was seeded with aluminum oxide particles. Planar laser light provided a two-dimensional projection of the flow field, and a traverse simulated a three-dimensional image of the entire flow field. Light focusing and reflection rendered the spheres black, permitting visualization of the planar circular interfaces in both the axial and transverse directions. Flows were observed near the wall-sphere interface and within the set of spheres. The CFD model required that a representative section of a packed bed be formed and gridded, enclosing and cutting six spheres so that symmetry conditions could be imposed at all cross-boundaries. Simulations had to be made with the flow direction at right angles to that used in the experiments, however, to take advantage of flow symmetry. Careful attention to detail was required for proper gridding. The flow field was three-dimensional and complex to describe, yet the most prominent finding was flow threads, as computed in the representative 'cube' of spheres with face symmetry and conclusively demonstrated experimentally herein. Random packing and bed voids tended to disrupt the laminar flow, creating vortices.

Description: Midspan aerodynamic measurements for a three vane-four passage linear turbine vane cascade are given. The vane axial chord was 4.45 cm. Surface pressures and loss coefficients were measured at exit Mach numbers of 0.3, 0.7, and 0.9. Reynolds number was varied by a factor of six at the two highest Mach numbers, and by a factor of ten at the lowest Mach number. Measurements were made with and without a turbulence grid. Inlet turbulence intensities were less than I% and greater than IO%. Length scales were also measured. Pressurized air fed the test section, and exited to a low pressure exhaust system. Maximum inlet pressure was two atmospheres. The minimum inlet pressure for an exit Mach number of 0.9 was one-third of an atmosphere, and at a Mach number of 0.3, the minimum pressure was half this value. The purpose of the test was to provide data for verification of turbine vane aerodynamic analyses, especially at low Reynolds numbers. Predictions obtained using a Navier-Stokes analysis with an algebraic turbulence model are also given.

Description: A simple mathematical model is described to account for the losses incurred when the flow in the duct (port) of a wave rotor is not aligned with the passages. The model, specifically for wave rotors with axially aligned passages, describes a loss mechanism which is sensitive to incident flow angle and Mach number. Implementation of the model in a one-dimensional CFD based wave rotor simulation is presented. Comparisons with limited experimental results are consistent with the model. Sensitivity studies are presented which highlight the significance of the incidence loss relative to other loss mechanisms in the wave rotor.

Description: Gravity modulation of an unbounded fluid layer with surface tension variations along its free surface is investigated. In parameter space of (wavenumber, Marangoni number) modulation has a destabilizing effect on the unmodulated neutral stability curve for large Prandtl number, Pr, and small modulation frequency, Omega, while a stabilizing effect is observed for small Pr and large Omega. As Omega yields infinity, the modulated neutral stability curves approach the unmodulated neutral stability curve. At certain values of Pr and L2 multiple minima are observed and the neutral stability curves become highly distorted. Closed regions of subharmonic instability are also observed. Alternating regions of synchronous and subharmonic instability separated by very thin stable regions are observed in (1/Omega,g(sub 1)) space for the singly diffusive cases. Quasiperiodic behavior in addition to the synchronous and subharmonic responses, are observed for the case of a double diffusive fluid layer. Minimum acceleration amplitudes were observed to closely correspond with a subharmonic response, Lambda(sub im) = Omega/2 .

Description: A metallic 'glove' structure has been built and attached to the wing of the Pegasus(trademark) space booster. An experiment on the upper surface of the glove has been designed to help validate boundary-layer stability codes in a free-flight environment. Three-dimensional thermal analyses have been performed to ensure that the glove structure design would be within allowable temperature limits in the experiment test section of the upper skin of the glove. Temperature results obtained from the design-case analysis show a peak temperature at the leading edge of 490 F. For the upper surface of the glove, approximately 3 in. back from the leading edge, temperature calculations indicate transition occurs at approximately 45 sec into the flight profile. A worst-case heating analysis has also been performed to ensure that the glove structure would not have any detrimental effects on the primary objective of the Pegasus a launch. A peak temperature of 805 F has been calculated on the leading edge of the glove structure. The temperatures predicted from the design case are well within the temperature limits of the glove structure, and the worst-case heating analysis temperature results are acceptable for the mission objectives.

Description: Accurate and efficient direct numerical simulation of turbulence in the presence of shock waves represents a significant challenge for numerical methods. The objective of this paper is to evaluate the performance of high order compact and non-compact central spatial differencing employing total variation diminishing (TVD) shock-capturing dissipations as characteristic based filters for two model problems combining shock wave and shear layer phenomena. A vortex pairing model evaluates the ability of the schemes to cope with shear layer instability and eddy shock waves, while a shock wave impingement on a spatially-evolving mixing layer model studies the accuracy of computation of vortices passing through a sequence of shock and expansion waves. A drastic increase in accuracy is observed if a suitable artificial compression formulation is applied to the TVD dissipations. With this modification to the filter step the fourth-order non-compact scheme shows improved results in comparison to second-order methods, while retaining the good shock resolution of the basic TVD scheme. For this characteristic based filter approach, however, the benefits of compact schemes or schemes with higher than fourth order are not sufficient to justify the higher complexity near the boundary and/or the additional computational cost.

Description: The ability to predict the dynamics of integrated inlet/compressor systems is an important part of designing high-speed propulsion systems. The boundaries of the performance envelope are often defined by undesirable transient phenomena in the inlet (unstart, buzz, etc.) in response to disturbances originated either in the engine or in the atmosphere. Stability margins used to compensate for the inability to accurately predict such processes lead to weight and performance penalties, which translate into a reduction in vehicle range. The prediction of transients in an inlet/compressor system requires either the coupling of two complex, unsteady codes (one for the inlet and one for the engine) or else a reliable characterization of the inlet/compressor interface, by specifying a boundary condition. In the context of engineering development programs, only the second option is viable economically. Computations of unsteady inlet flows invariably rely on simple compressor-face boundary conditions (CFBC's). Currently, customary conditions include choked flow, constant static pressure, constant axial velocity, constant Mach number or constant mass flow per unit area. These conditions are straightforward extensions of practices that are valid for and work well with steady inlet flows. Unfortunately, it is not at all likely that any flow property would stay constant during a complex system transient. At the start of this effort, no experimental observation existed that could be used to formulate of verify any of the CFBC'S. This lack of hard information represented a risk for a development program that has been recognized to be unacceptably large. The goal of the present effort was to generate such data. Disturbances reaching the compressor face in flight may have complex spatial structures and temporal histories. Small amplitude disturbances may be decomposed into acoustic, vorticity and entropy contributions that are uncoupled if the undisturbed flow is uniform. This study is focused on the response of an inlet/compressor system to acoustic disturbances. From the viewpoint of inlet computations, acoustic disturbances are clearly the most important, since they are the only ones capable of moving upstream. Convective and entropy disturbances may also produce upstream-moving acoustic waves, but such processes are outside the scope of the present study.

Description: This paper presents application of a Generalized Fluid System Simulation Program (GFSSP) to model fluid flow in a very complex secondary flow circuit of a rocket engine turbopump, SIMPLEX. SIMPLEX was a low cost turbopump designed, manufactured and tested to demonstrate the ability to reduce the overall cost and design cycle time of turbomachinery. GFSSP is a general purpose computer program for analyzing flowrate, pressure, temperature and concentration distribution in a complex flow network. The program employs a finite volume formulation of mass, momentum and energy conservation equations in conjunction with thermodynamic equation of state of a real fluid. The secondary flow passages modeled include flow between the impeller shroud and housing, impeller back face and housing, flow through labyrinth seal and bearing, return line flow to inducer and flow adjacent to front and back face of the turbine. Axial load on the bearing are calculated from the predicted pressure on all axial faces of the rotor. The predicted pressure and temperature distributions compared well with the test data.

Description: A multi-block, three-dimensional Navier-Stokes code has been used to study the within-hole and near-hole physics in relation to heat transfer on a film-cooled blade. The flow domain consists of the coolant flow through the plenum and hole-pipes for the three staggered rows of shower-head holes on the VK1 rotor, and the main flow over the blade. A multi-block grid is generated that is nearly orthogonal to the various surfaces. It may be noted that for the VK1 rotor the shower-head holes are inclined at 30 deg. to the spanwise direction, and are normal to the streamwise direction on the blade. Wilcox's k-omega turbulence model is used. The present study provides a much better comparison for the heat transfer coefficient at the blade mid-span with the experimental data than an earlier analysis wherein coolant velocity and temperature distributions were specified at the hole exits rather than extending the computational domain into the hole-pipe and plenum. Details of the distributions of coolant velocity, temperature, k and omega at the hole exits are also presented.

Description: A recently reported solution for stationary stability of a thermosolutal system with Soret diffusion is re-derived and examined using a symbolic computational package. Symbolic computational languages are well suited for such an analysis and facilitate a pragmatic approach that is adaptable to similar problems. Linearization of the equations, normal mode analysis, and extraction of the final solution are performed in a Mathematica notebook format. An exact solution is obtained for stationary stability in the limit of zero gravity. A closed form expression is also obtained for the location of asymptotes in relevant parameter, (Sm(sub c), Mac(sub c)), space. The stationary stability behavior is conveniently examined within the symbolic language environment. An abbreviated version of the Mathematica notebook is given in the Appendix.

Description: A program is underway to improve the efficiency of a three-dimensional Navier-Stokes code and generalize it for nozzle and turbopump geometries. Code modifications have included the implementation of parallel processing software, incorporation of new physical models and generalization of the multiblock capability. The final report contains details of code modifications, numerical results for several nozzle and turbopump geometries, and the implementation of the parallelization software.

Description: Recently, perfectly matched layer (PML) as an absorbing boundary condition has widespread applications. The idea was first introduced by Berenger for electromagnetic waves computations. In this paper, it is shown that the PML equations for the linearized Euler equations support unstable solutions when the mean flow has a component normal to the layer. To suppress such unstable solutions so as to render the PML concept useful for this class of problems, it is proposed that artificial selective damping terms be added to the discretized PML equations. It is demonstrated that with a proper choice of artificial mesh Reynolds number, the PML equations can be made stable. Numerical examples are provided to illustrate that the stabilized PML performs well as an absorbing boundary condition. In a ducted environment, the wave mode are dispersive. It will be shown that the group velocity and phase velocity of these modes can have opposite signs. This results in a confined environment, PML may not be suitable as an absorbing boundary condition.

Description: Two-dimensional thermosolutal convection with no-slip boundary conditions is studied using numerical simulations in a periodic domain. The domain is large enough to follow the evolution of phase instabilities of fully nonlinear traveling waves. In the parameter regime studied these instabilities evolve, without loss of phase or hysteresis, into a series of confined states or pulses characterized by locally enhanced heat and solute transport. The wavelength and phase velocity of the traveling rolls within a pulse differ substantially from those in the background. The pulses drift in the same direction as the convection rolls on which they ride but more slowly, and are characterized by an exponential leading front and an oscillatory trailing end. Multiple, apparently stable, states are found for identical parameter values. The qualitative properties of the pulses are in good agreement with the predictions of a third-order phase equation which accounts for the relation between wave number and phase velocity, the oscillatory tails and the multiplicity of states. These properties of the pulses are shown to be a consequence of Shil'nikov dynamics in the spatial domain.

Description: Two-phase flows are present in a wide variety of applications for spacecraft thermal control systems design. Bubble formation and detachment is an integral part of the two-phase flow science. The authors objective is to experimentally investigate the effects of liquid cross velocity, gas velocity, and nozzle diameter on the bubble diameter at detachment under reduced and normal gravity and under relatively low gas flow rates. Results from ground (1 g) and reduced gravity experiments will be presented in this paper. For the 1 g experiment, a flow loop was designed and built to accommodate a range of liquid and gas flow rates. The reduced gravity experiment was conducted on the NASA DC-9 reduced gravity platform using the two-phase flow loop qualified for operation on the low-gravity platform. Flow visualization is accomplished using a high speed 500 frames/s camera. The results suggest that the existence of buoyancy force contributes to the faster detachment of bubbles. Buoyancy helps the detachment process which results in smaller bubbles being formed. In reduced gravity, although drag force is present, the virtual non-existence of buoyant force results in larger bubbles and longer times for detachment. Theoretical predictions are also presented in this paper and seem to agree with the experimental results.

Description: This paper describes the single passage model based on CFD to analyze the flow in blade passages of a centrifugal pump. The model consists of the flow passage between two impeller blades and the spaces in the inlet eye as well as in the volute. The incompressible Navier-Stokes equations in the conservation form are solved by a finite difference method. The code is designed to investigate the velocity and pressure distributions and intended to investigate how the pump design affects fluid flow through the rotor as well as the pump performance. An early part of the paper investigates the behavior of the model as well as validity of the assumptions made in the model. Then, applications to a rotodynamic heart pump are presented.

Description: Aircraft trailing vortices can be influenced significantly by atmospheric conditions such as crosswind, turbulence, and stratification. According to the NASA 1994 and 1995 field measurement program in Memphis, Tennessee, the descending aircraft wake vortices could stall or be deflected at the top of low-level temperature inversions that usually produce pronounced shear zones. Numerical simulations of vortex/shear interactions with ground effects have been performed by several groups. Burnham used a series of evenly spaced line vortices at a particular altitude to model the ground shear layer of the cross- wind. He found that the wind shear was swept up around the downwind vortex and caused the downwind vortex to move upward, and claimed that the effect was actually produced by the vertical gradient in the wind shear rather than by the wind shear directly, because uniformly distributed wind-shear vortices would have no effect on the trailing vortex vertical motion. Recently, Proctor et al. numerically tested the effects of narrow shear zones on the behavior of the vortex pair, motivated by the observation of the Memphis field data. The shear-layer sensitivity tests indicated that the downwind vortex was more sensitive and deflected to a higher altitude than its upwind counterpart. The downstream vortex contained vorticity of opposite sign to that of the shear. There was no detectable preference for the downwind vortex (or upwind vortex) to weaken (or strengthen) at a greater rate.

Description: The Strutjet approach to Rocket Based Combined Cycle (RBCC) propulsion depends upon fuel-rich flows from the rocket nozzles and turbine exhaust products mixing with the ingested air for successful operation in the ramjet and scramjet modes. It is desirable to delay this mixing process in the air-augmented mode of operation present during take-off and low speed flight. A scale model of the Strutjet device was built and tested to investigate the mixing of the streams as a function of distance from the Strut exit plane in simulated sea level take-off conditions. The Planar Laser Induced Fluorescence (PLIF) diagnostic method has been employed to observe the mixing of the turbine exhaust gas with the gases from both the primary rockets and the ingested air. The ratio of the pressure in the turbine exhaust to that in the rocket nozzle wall at the point where the two jets meet, is the independent variable in these experiments. Tests were accomplished at values of 1.0 (the original design point), 1.5 and 2.0 for this parameter at 8 locations downstream of the rocket nozzle exit. The results illustrate the development of the mixing zone from the exit plane of the strut to a distance of about 18 equivalent rocket nozzle exit diameters downstream (18"). These images show the turbine exhaust to be confined until a short distance downstream. The expansion into the ingested air is more pronounced at a pressure ratio of 1.0 and 1.5 and shows that mixing with this air would likely begin at a distance of 2" downstream of the nozzle exit plane. Of the pressure ratios tested in this research, 2.0 is the best value for delaying the mixing at the operating conditions considered.

Description: The Seventh Summer Program of the Center for Turbulence Research took place in the four-week period, July 5 to July 31, 1998. This was the largest CTR Summer Program to date, involving thirty-six participants from the U. S. and nine other countries. Thirty-one Stanford and NASA-Ames staff members facilitated and contributed to most of the Summer projects. A new feature, and perhaps a preview of the future programs, was that many of the projects were executed on non-NASA computers. These included supercomputers located in Europe as well as those operated by the Departments of Defense and Energy in the United States. In addition, several simulation programs developed by the visiting participants at their home institutions were used. Another new feature was the prevalence of lap-top personal computers which were used by several participants to carry out some of the work that in the past were performed on desk-top workstations. We expect these trends to continue as computing power is enhanced and as more researchers (many of whom CTR alumni) use numerical simulations to study turbulent flows. CTR's main role continues to be in providing a forum for the study of turbulence for engineering analysis and in facilitating intellectual exchange among the leading researchers in the field. Once again the combustion group was the largest. Turbulent combustion has enjoyed remarkable progress in using simulations to address increasingly complex and practically more relevant questions. The combustion group's studies included such challenging topics as fuel evaporation, soot chemistry, and thermonuclear reactions. The latter study was one of three projects related to the Department of Energy's ASCI Program (www.llnl.gov/asci); the other two (rocket propulsion and fire safety) were carried out in the turbulence modeling group. The flow control and acoustics group demonstrated a successful application of the so-called evolution algorithms which actually led to a previously unknown forcing strategy for jets yielding increased spreading rate. A very efficient algorithm for flow in complex geometries with moving boundaries based on the immersed boundary forcing technique was tested with very encouraging results. Also a new strategy for the destruction of aircraft trailing vortices was introduced and tested. The Reynolds Averaged Modeling (RANS) group demonstrated that the elliptic relaxation concept for RANS calculations is also applicable to transonic flows with shocks; however, prediction of laminar/turbulent transition remains an important pacing item. A large fraction of the LES effort was devoted to the development and testing of a new algorithmic procedure (as opposed to phenomenological model) for subgrid scale modeling based on regularized de-filtering of the flow variables. This appears to be a very promising approach, and a significant effort is currently underway to assess its robustness in high Reynolds number flows and in conjunction with numerical methods for complex flows. As part of the Summer Program two review tutorials were given on Turbulent structures in hydrocarbon pool fires (Sheldon Tieszen), and Turbulent combustion modeling: from RANS to LES via DNS (Luc Vervisch); and two seminars entitled Assessment of turbulence models for engineering applications (Paul Durbin) and Subgrid-scale modeling for non-premixed, turbulent reacting flows (James Riley) were presented. A number of colleagues from universities, government agencies, and industry attended the final presentations of the participants on July 31 and participated in the discussions. There are twenty-six papers in this volume grouped in five areas. Each group is preceded with an overview by its coordinator.

Description: A pre-release version of the Navier-Stokes solver (TURBO) was obtained from MSU. Along with Dr. Milind Bakhle of the University of Toledo, subroutines for aeroelastic analysis were developed and added to the TURBO code to develop versions 1 and 2 of the TURBO-AE code. For specified mode shape, frequency and inter-blade phase angle the code calculates the work done by the fluid on the rotor for a prescribed sinusoidal motion. Positive work on the rotor indicates instability of the rotor. The version 1 of the code calculates the work for in-phase blade motions only. In version 2 of the code, the capability for analyzing all possible inter-blade phase angles, was added. The version 2 of TURBO-AE code was validated and delivered to NASA and the industry partners of the AST project. The capabilities and the features of the code are summarized in Refs. [1] & [2]. To release the version 2 of TURBO-AE, a workshop was organized at NASA Lewis, by Dr. Srivastava and Dr. M. A. Bakhle, both of the University of Toledo, in October of 1996 for the industry partners of NASA Lewis. The workshop provided the potential users of TURBO-AE, all the relevant information required in preparing the input data, executing the code, interpreting the results and bench marking the code on their computer systems. After the code was delivered to the industry partners, user support was also provided. A new version of the Navier-Stokes solver (TURBO) was later released by MSU. This version had significant changes and upgrades over the previous version. This new version was merged with the TURBO-AE code. Also, new boundary conditions for 3-D unsteady non-reflecting boundaries, were developed by researchers from UTRC, Ref. [3]. Time was spent on understanding, familiarizing, executing and implementing the new boundary conditions into the TURBO-AE code. Work was started on the phase lagged (time-shifted) boundary condition version (version 4) of the code. This will allow the users to calculate non-zero interblade phase angles using, only one blade passage for analysis.

Description: Mathematical theory predicts that small changes in container shape or in contact angle can give rise to large shifts of liquid in a microgravity environment. For studying such shifts, container shapes are described that were flown on board the Space Shuttle USML-2 mission as part of the Glovebox Interface Configuration Experiment. These containers are in the form of a circular cylinder with two diametrically opposed "canonical proboscis" protrusions. The containers were designed with the goal of having two desirable properties -that sufficient liquid would participate in the shift to permit easy observation, and that the change would be abrupt enough to allow accurate determination of critical contact angle. The observed behavior for these vessels is depicted, along with behavior for a movable wedge vessel, which also formed part of the experiment. The experimental results support the validity of the concept of macroscopic contact angle, basic to the theory, and thereby its use in predicting fluid behavior under reduced gravity. The results indicate, as well, the role of hysteresis in impeding orientation to equilibrium.