ALBERT

Description: A methodology for the simulation of turbine cooling flows is being developed. The methodology seeks to combine numerical techniques that optimize both accuracy and computational efficiency. Key components of the methodology include the use of multiblock grid systems for modeling complex geometries, and multigrid convergence acceleration for enhancing computational efficiency in highly resolved fluid flow simulations. The use of the methodology has been demonstrated in several turbo machinery flow and heat transfer studies. Ongoing and future work involves implementing additional turbulence models, improving computational efficiency, adding AMR.

Description: The following are conclusions and recomendations from the study. Primary wake effect is linear reduction in (eta) with St. Secondary wake effect is skewing of suction/pressure side cooling. Steady computations match experimental Nu, but overpredict (eta). Unsteady computations elucidate wake/film interaction. Model may be used to estimate wake passing effect. Need boundary layer and full stage experiments. Need resolved film hole and full stage unsteady computations. Need validated turbulence models for film cooling.

Description: The following topics are discussed: (1) The Transient Liquid-Crystal Heat-Transfer Technique; (2) 2-D Film-Cooling Heat-Transfer on an AlliedSignal Vane; and (3) Effects of Tab Vortex Generators on Surface Heat Transfer. Downstream of a Jet in Crossflow.

Description: A tab placed on the leeward side of the nozzle was expected to increase jet penetration into the cross-flow. An experiment at UTRC showed insignificant effect. The primary objective of the present study was to confirm and explain the ineffectiveness. The overall approach of the study was to conduct experiments in a low-speed wind tunnel and to conduct hot-wire measurements for mean velocity and streamwise vorticity fields.

Description: Testing is currently in progress to acquire simple geometry surface heat transfer data for internal channels with trips and bleed holes which can be used in the development and validation of models. The transient liquid crystal technique is used on a simple multipass model with rectangular channels and normal ribs. Normal bleed holes are located on the floor of the model in the first channel. Each hole is attached to a flow meter, allowing various bleed flow rates to simulate external pressures on the blade.

Description: Reasonable heat transfer prediction can be achieved in complex geometries. Multi-block grid allows efficient placement of grid points, and efficient use of computer resources. Wilcox k-(omega) turbulence model predicts heat transfer well, and has good numerical behavior.

Description: The SPAcecraft SIMulator (SPASIM) simulates the functions and resources of a spacecraft to quickly perform conceptual design (Phase A) trade-off and sensitivity analyses and uncover any operational bottlenecks during any part of the mission. Failure modes and operational contingencies can be evaluated allowing operational planning (what-if scenarios) and optimization for a range of mission scenarios. The payloads and subsystems are simulated, using a hierarchy of graphical models, in terms of how their functions affect resources such as propellant, power, and data. Any of the inputs and outputs of the payloads and subsystems can be plotted during the simulation or stored in a file so they can be used by other programs. Most trade-off analyses, including those that compare current versus advanced technology, can be performed by changing values in the parameter menus. However, when a component is replaced by one with a different functional architecture, its graphical model can also be modified or replaced by drawing from a component library. SPASIM has been validated using several spacecraft designs that were at least at the Critical Design Review level. The user and programmer guide, including figures, is available on line as a hypertext document. This is an easy-to-use and expandable tool which is based on MATLAB(R) and SIMULINK(R). It runs on Silicon Graphics Inc. workstations and personal computers with Windows 95(TM) or NT(TM).

Description: Detailed flow field measurements have been carried out for a turbulent circular jet perturbed by tabs and artificial excitation. Two "delta tabs" were placed at the nozzle exit at diametrically opposite locations. The excitation condition involved subharmonic resonance that manifested in a periodic vortex pairing in the near flow field. While the excitation and the tabs independently increased jet spreading, a combination of the two diminished the effect. The jet spreading was most pronounced with the tabs but was reduced when excitation was applied to the tabbed jet. The tabs generated streamwise vortex pairs that caused a lateral spreading of the jet in a direction perpendicular to the plane containing the tabs. The excitation, on the other band, organized the azimuthal vorticity into coherent ring structures whose evolution and pairing also increased entrainment by the jet. In the tabbed case, the excitation produced coherent azimuthal structures that were distorted and asymmetric in shape. The self-induction of these structures produced an effect that opposed the tendency for the lateral spreading of the streamwise vortex pairs. The passage of the distorted vortices, and their pairing, also had a cancellation effect on the time-averaged streamwise vorticity field. These led to the reduction in jet spreading.

Description: The interaction between the main passage flow through a high pressure turbine and the secondary flows due to leakage through the wheelspace rim seals is reviewed. Various experimental and analytical studies of these interactions are outlined and a numerical investigation of the interaction between the main passage flow and a wheelspace cavity seal flow is described. The numerical investigation uses a structured grid method to study the overall interaction between the turbine stage components and the wheelspace seal flow, and an unstructured grid method to resolve the detailed flow features within the geometrically complex cavity seal. The numerical results agree with various observations from experimental studies under similar flow conditions. As the flow rate through the rim cavity seal is increased, the ingestion of fluid from the main passage flow into the rim seal area decreases rapidly. A small amount of main flow gas is ingested into the rim seal area, even at high seal flow rates.

Description: Cartesian methods for computational fluid dynamics are considered to offer an accurate and robust approach for the simulation of aerodynamic flows around geometrically complex bodies. A part of the ongoing research in this domain is reviewed with the aim of providing insight into the fundamental challenges faced by the practitioners of this approach, and a guide to further research. The integration schemes used in Cartesian solvers are similar to those used in other approaches. The emphasis is on the geometric algorithms, surface modeling and boundary conditions required to design a successful Cartesian mesh scheme.

Description: In the past forty years much progress has been made in computational methods applied to the solution of problems in spacecraft hypervelocity flow and heat transfer. Although the basic thermochemical and physical modeling techniques have changed little in this time, several orders of magnitude increase in the speed of numerically solving the Navier-Stokes and associated energy equations have been achieved. The extent to which this computational power can be applied to the design of spacecraft heat shields is dependent on the proper coupling of the external flow equations to the boundary conditions and governing equations representing the thermal protection system in-depth conduction, pyrolysis and surface ablation phenomena. A discussion of the techniques used to do this in past problems as well as the current state-of-art is provided. Specific examples, including past missions such as Galileo, together with the more recent case studies of ESA/Rosetta Sample Comet Return, Mars Pathfinder and X-33 will be discussed. Modeling assumptions, design approach and computational methods and results are presented.

Description: A loosely coupled two-phase vacuum water plume model has been developed. This model consists of a source flow model to describe the expansion of water vapor, and the Lagrangian equations of motion for particle trajectories. Gas/Particle interaction is modeled through the drag force induced by the relative velocities. Particles are assumed traveling along streamlines. The equations of motion are integrated to obtain particle velocity along the streamline. This model has been used to predict the mass flux in a 5 meter radius hemispherical domain resulting from the burst of a water jet of 1.5 mm in diameter, mass flow rate of 24.2 g/s, and stagnation pressure of 21.0 psia, which is the nominal Orbiter water dump condition. The result is compared with an empirical water plume model deduced from a video image of the STS-29 water dump. To further improve the model, work has begun to numerically simulate the bubble formation and bursting present in a liquid stream injected into a vacuum. The technique of smoothed particle hydrodynamics was used to formulate this simulation. A status and results of the on-going effort are presented and compared to results from the literature.

Description: A distributed algorithm for a high-order-accurate finite-difference approach to the direct numerical simulation (DNS) of transition and turbulence in compressible flows is described. This work has two major objectives. The first objective is to demonstrate that parallel and distributed-memory machines can be successfully and efficiently used to solve computationally intensive and input/output intensive algorithms of the DNS class. The second objective is to show that the computational complexity involved in solving the tridiagonal systems inherent in the DNS algorithm can be reduced by algorithm innovations that obviate the need to use a parallelized tridiagonal solver.

Description: The behavior of velocity fluctuations subjected to rapid rotation is examined. The rapid rotation considered is any arbitrary combination of two basic forms of rotation, reference frame rotation and mean flow rotation. It is recognized that the two types of rotating flows differ in the manner in which the fluctuating fields are advected. The first category is comprised of flows in rotating systems of which synoptic scale geophysical flows are a good example. In this class of flows the fluctuating velocity field advects and rotates with the mean flow. In the rapid rotation limit, the Taylor-Proudman theorem describes the behavior of this class of fluctuations. Velocity fluctuations that are advected without rotation by the mean flow constitute the second category which includes vortical flows of aerodynamic interest. The Taylor-Proudman theorem is not pertinent to I his class flows and a new result appropriate to this second category of fluctuations is derived. The present development demonstrates that the fluctuating velocity fields are rendered two-dimensional and horizontally non-divergent in the limit of any large combination of reference frame rotation and mean-flow rotation. The concommitant 'geostrophic' balance of the momentum equation is, however, dependent upon the form of rapid rotation. It is also demonstrated that the evolution equations of a two-dimensional fluctuating velocity fields are frame-indifferent with any imposed mean-flow rotation. The analyses and results of this paper highlight many fundamental aspects of rotating flows and have important consequences for their turbulence closures in inertial and non-inertial frames.

Description: A free-floating droplet in microgravity is ideal for scientific observation since it is free of confounding factors such as wetting and nonsymmetrical heat transfer introduced by contact with surfaces. However, the technology to reliably deploy in microgravity has not yet been developed. In some recent fluid deployment experiments, droplets are either shaken off the dispenser or the dispenser is quickly retracted from the droplet. These solutions impart random residual motion to deployed droplet, which can be undesirable for certain investigations. In the present study, two new types of droplet injectors were built and tested. Testing of the droplet injectors consisted of neutral buoyancy tank tests, 5-sec drop tower tests at the NASA Lewis Zero Gravity Facility, and DC-9 tests. One type, the concentric injector, worked well in the neutral buoyancy tank but did not do well in low-gravity. However, it appeared that it makes a fine apparatus for constructing bubbles in low-gravity conditions. The other type, the T-injector, showed the most promise for future development. In both neutral buoyancy and DC-9 tests, water droplets were formed and deployed with some control and repeatability, although in low-gravity the residual velocities were higher than desirable. Based on our observations, further refinements are suggested for future development work.

Description: The term renormalization and renormalization group are explained by reference to various physical systems. The extension of renormalization group to turbulence is then discussed; first as a comprehensive review and second concentrating on the technical details of a few selected approaches. We conclude with a discussion of the relevance and application of renormalization group to turbulence modelling.

Description: This report describes the use of novel techniques for heat containment that could be used to build a high temperature transparent furnace. The primary objective of the work was to experimentally demonstrate transparent furnace operation at 1200 C. Secondary objectives were to understand furnace operation and furnace component specification to enable the design and construction of a low power prototype furnace for delivery to NASA in a follow-up project. The basic approach of the research was to couple high temperature component design with simple concept demonstration experiments that modify a commercially available transparent furnace rated at lower temperature. A detailed energy balance of the operating transparent furnace was performed, calculating heat losses through the furnace components as a result of conduction, radiation, and convection. The transparent furnace shells and furnace components were redesigned to permit furnace operation at at least 1200 C. Techniques were developed that are expected to lead to significantly improved heat containment compared with current transparent furnaces. The design of a thermal profile in a multizone high temperature transparent furnace design was also addressed. Experiments were performed to verify the energy balance analysis, to demonstrate some of the major furnace improvement techniques developed, and to demonstrate the overall feasibility of a high temperature transparent furnace. The important objective of the research was achieved: to demonstrate the feasibility of operating a transparent furnace at 1200 C.

Description: Turbulent wall-bounded complex flows are commonly encountered in engineering practice and are of considerable interest in a variety of industrial applications. The presence of a wall significantly affects turbulence characteristics. In addition to the wall effects, turbulent wall-bounded flows become more complicated by the presence of additional body forces (e.g. centrifugal force and Coriolis force) and complex geometry. Most near-wall Reynolds stress models are developed from a high-Reynolds-number model which assumes turbulence is homogenous (or quasi-homogenous). Near-wall modifications are proposed to include wall effects in near-wall regions. In this process, wall normals are introduced. Good predictions could be obtained by Reynolds stress models with wall normals. However, ambiguity arises when the models are applied in flows with multiple walls. Many models have been proposed to model turbulent flows. Among them, Reynolds stress models, in which turbulent stresses are obtained by solving the Reynolds stress transport equations, have been proved to be the most successful ones. To apply the Reynolds stress models to wall-bounded flows, near-wall corrections accounting for the wall effects are needed, and the resulting models are called near-wall Reynolds stress models. In most of the existing near-wall models, the near-wall corrections invoke wall normals. These wall-dependent near-wall models are difficult to implement for turbulent flows with complex geometry and may give inaccurate predictions due to the ambiguity of wall normals at corners connecting multiple walls. The objective of this study is to develop a more general and flexible near-wall Reynolds stress model without using any wall-dependent variable for wall-bounded turbulent flows. With the aid of near-wall asymptotic analysis and results of direct numerical simulation, a new near-wall Reynolds stress model (NNWRS) is formulated based on Speziale et al.'s high-Reynolds-stress model with wall-independent near-wall corrections. Moreover, only one damping function is used for flows with a wide range of Reynolds numbers to ensure that the near-wall modifications diminish away from the walls.

Description: Despite significant achievements in computational fluid dynamics, there still remain many fluid flow phenomena not well understood. For example, the prediction of temperature distributions is inaccurate when temperature gradients are high, particularly in shock wave turbulent boundary layer interactions close to the wall. Complexities of fluid flow phenomena include transition to turbulence, relaminarization separated flows, transition between viscous and inviscid incompressible and compressible flows, among others, in all speed regimes. The purpose of this paper is to introduce a new approach, called the Flowfield-Dependent Mixed Explicit-Implicit (FDMEI) method, in an attempt to resolve these difficult issues in Computational Fluid Dynamics (CFD). In this process, a total of six implicitness parameters characteristic of the current flowfield are introduced. They are calculated from the current flowfield or changes of Mach numbers, Reynolds numbers, Peclet numbers, and Damkoehler numbers (if reacting) at each nodal point and time step. This implies that every nodal point or element is provided with different or unique numerical scheme according to their current flowfield situations, whether compressible, incompressible, viscous, inviscid, laminar, turbulent, reacting, or nonreacting. In this procedure, discontinuities or fluctuations of an variables between adjacent nodal points are determined accurately. If these implicitness parameters are fixed to certain numbers instead of being calculated from the flowfield information, then practically all currently available schemes of finite differences or finite elements arise as special cases. Some benchmark problems to be presented in this paper will show the validity, accuracy, and efficiency of the proposed methodology.

Description: Three equilibrium-air numerical solutions are presented for the Reentry-F flight-test vehicle at Mach 20, 80,000 Ft. conditions, including turbulent flow predictions. The three solutions are from a thin-layer Navier-Stokes code, coupled thin-layer and parabolized Navier-Stokes codes, and an approximate viscous shock-layer code. Boundary-layer and shock-layer profiles are presented and compared between the three solutions, revealing close agreement between the three solution methods. Notable exceptions to the close agreement, with 7-10 percent discrepancies, occur in the density profiles at the boundary-layer edge, in the boundary-layer velocity profiles, and in the shock-layer profiles in regions influenced by the nose bluntness.

Description: In this paper, we examine the effectiveness of absorbing layers as non-reflecting computational boundaries for the Euler equations. The absorbing-layer equations are simply obtained by splitting the governing equations in the coordinate directions and introducing absorption coefficients in each split equation. This methodology is similar to that used by Berenger for the numerical solutions of Maxwell's equations. Specifically, we apply this methodology to three physical problems shock-vortex interactions, a plane free shear flow and an axisymmetric jet- with emphasis on acoustic wave propagation. Our numerical results indicate that the use of absorbing layers effectively minimizes numerical reflection in all three problems considered.

Description: The objective of this investigation is experimental and computational study of three dimensional viscous flow field in the nozzle passage of an axial flow turbine stage. The nozzle passage flow field has been measured using a two sensor hot-wire probe at various axial and radial stations. In addition, two component LDV measurements at one axial station (x/c(sum m) = 0.56) were performed to measure the velocity field. Static pressure measurements and flow visualization, using a fluorescent oil technique, were also performed to obtain the location of transition and the endwall limiting streamlines. A three dimensional boundary layer code, with a simple intermittency transition model, was used to predict the viscous layers along the blade and endwall surfaces. The boundary layers on the blade surface were found to be very thin and mostly laminar, except on the suction surface downstream of 70% axial chord. Strong radial pressure gradient, especially close to the suction surface, induces strong cross flow components in the trailing edge regions of the blade. On the end-walls the boundary layers were much thicker, especially near the suction corner of the casing surface, caused by secondary flow. The secondary flow region near the suction-casing surface corner indicates the presence of the passage vortex detached from the blade surface. The corner vortex is found to be very weak. The presence of a closely spaced rotor downstream (20% of the nozzle vane chord) introduces unsteadiness in the blade passage. The measured instantaneous velocity signal was filtered using FFT square window to remove the periodic unsteadiness introduced by the downstream rotor and fans. The filtering decreased the free stream turbulence level from 2.1% to 0.9% but had no influence on the computed turbulence length scale. The computation of the three dimensional boundary layers is found to be accurate on the nozzle passage blade surfaces, away from the end-walls and the secondary flow region. On the nozzle passage endwall surfaces the presence of strong pressure gradients and secondary flow limit the validity of the boundary layer code.

Description: The interaction between a shock wave (attached to a wedge) and small amplitude, three-dimensional disturbances of a uniform, supersonic, freestream flow are investigated. The paper extends the two-dimensional study of Duck et al, through the use of vector potentials, which render the problem tractable by the same techniques as in the two-dimensional case, in particular by expansion of the solution by means of a Fourier-Bessel series, in appropriately chosen coordinates. Results are presented for specific classes of freestream disturbances, and the study shows conclusively that the shock is stable to all classes of disturbances (i.e. time periodic perturbations to the shock do not grow downstream), provided the flow downstream of the shock is supersonic (loosely corresponding to the weak shock solution). This is shown from our numerical results and also by asymptotic analysis of the Fourier-Bessel series, valid far downstream of the shock.

Description: Measurements were made in a large scale single stage turbine facility. Within the nozzle passage measurements were made using a five hole probe, a two-component Laser Doppler Velocimeter (LDV), and a single sensor hot wire probe. These measurements showed weak secondary flows at midchord, and two secondary flow loss cores at the nozzle exit. The casing vortex loss core was the larger of the two. At the exit radial inward flow was found over the entire passage, and was more pronounced in the wake. Nozzle wake decay was found to be more rapid than for an isolated vane row due to the rotor's presence. The midspan rotor flow field was measured using a two-component LDV. Measurements were made from upstream of the rotor to a chord behind the rotor. The distortion of the nozzle wake as it passed through the rotor blade row was determined. The unsteadiness in the rotor flow field was determined. The decay of the rotor wake was also characterized.

Description: Three-dimensional, incompressible turbulent jets with rectangular and elliptical cross-sections are simulated with a finite-difference numerical method. The full Navier- Stokes equations are solved at low Reynolds numbers, whereas at high Reynolds numbers filtered forms of the equations are solved along with a sub-grid scale model to approximate the effects of the unresolved scales. A 2-N storage, third-order Runge-Kutta scheme is used for temporary discretization and a fourth-order compact scheme is used for spatial discretization. Although such methods are widely used in the simulation of compressible flows, the lack of an evolution equation for pressure or density presents particular difficulty in incompressible flows. The pressure-velocity coupling must be established indirectly. It is achieved, in this study, through a Poisson equation which is solved by a compact scheme of the same order of accuracy. The numerical formulation is validated and the dispersion and dissipation errors are documented by the solution of a wide range of benchmark problems. Three-dimensional computations are performed for different inlet conditions which model the naturally developing and forced jets. The experimentally observed phenomenon of axis-switching is captured in the numerical simulation, and it is confirmed through flow visualization that this is based on self-induction of the vorticity field. Statistical quantities such as mean velocity, mean pressure, two-point velocity spatial correlations and Reynolds stresses are presented. Detailed budgets of the mean momentum and Reynolds stresses are presented. Detailed budgets of the mean momentum and Reynolds stress equations are presented to aid in the turbulence modeling of complex jets. Simulations of circular jets are used to quantify the effect of the non-uniform curvature of the non-circular jets.

Description: The flow characteristics in the base region of a two-dimensional supersonic compression ramp are investigated. A stream-wise oriented air jet, M = 1.75, is injected through a thin horizontal slot into a supersonic air main flow, M = 2.3, at the end of a two-dimensional compression ramp. The velocity profile and basic characteristics of the flow in the base region immediately following the ramp are determined. Visualization of the flowfield for qualitative observations is accomplished via Dark Central Ground Interferometry (DCGI). Two-dimensional velocity profiles are obtained using Laser Doppler Velocimetry (LDV). The study is the initial phase of a four-year investigation of base flow mixing. The current study is to provide more details of the flowfield.

Description: The current compute environment that most researchers are using for the calculation of 3D unsteady Computational Fluid Dynamic (CFD) results is a super-computer class machine. The Massively Parallel Processors (MPP's) such as the 160 node IBM SP2 at NAS and clusters of workstations acting as a single MPP (like NAS's SGI Power-Challenge array and the J90 cluster) provide the required computation bandwidth for CFD calculations of transient problems. If we follow the traditional computational analysis steps for CFD (and we wish to construct an interactive visualizer) we need to be aware of the following: (1) Disk space requirements. A single snap-shot must contain at least the values (primitive variables) stored at the appropriate locations within the mesh. For most simple 3D Euler solvers that means 5 floating point words. Navier-Stokes solutions with turbulence models may contain 7 state-variables. (2) Disk speed vs. Computational speeds. The time required to read the complete solution of a saved time frame from disk is now longer than the compute time for a set number of iterations from an explicit solver. Depending, on the hardware and solver an iteration of an implicit code may also take less time than reading the solution from disk. If one examines the performance improvements in the last decade or two, it is easy to see that depending on disk performance (vs. CPU improvement) may not be the best method for enhancing interactivity. (3) Cluster and Parallel Machine I/O problems. Disk access time is much worse within current parallel machines and cluster of workstations that are acting in concert to solve a single problem. In this case we are not trying to read the volume of data, but are running the solver and the solver outputs the solution. These traditional network interfaces must be used for the file system. (4) Numerics of particle traces. Most visualization tools can work upon a single snap shot of the data but some visualization tools for transient problems require dealing with time.

Description: A code has been developed to automatically simplify full chemical mechanisms. The method employed is based on the Intrinsic Low Dimensional Manifold (ILDM) method of Maas and Pope. The ILDM method is a dynamical systems approach to the simplification of large chemical kinetic mechanisms. By identifying low-dimensional attracting manifolds, the method allows complex full mechanisms to be parameterized by just a few variables; in effect, generating reduced chemical mechanisms by an automatic procedure. These resulting mechanisms however, still retain all the species used in the full mechanism. Full and skeletal mechanisms for various fuels are simplified to a two dimensional manifold, and the resulting mechanisms are found to compare well with the full mechanisms, and show significant improvement over global one step mechanisms, such as those by Westbrook and Dryer. In addition, by using an ILDM reaction mechanism in a CID code, a considerable improvement in turn-around time can be achieved.

Description: The viscous driven-cavity problem is solved using a stream-function and vorticity formulation for the incompressible Navier-Stokes equations. This report provides the user's manual and FORTRAN code for the set of governing equations presented in NASA TM-110262.

Description: Researchers are developing the technology of 'Ballistic Particle Manufacturing' (BPM) in which individual drops are precisely layered onto a substrate, and the drops are deposited so as to prevent splatting. These individual drops will ultimately be combined to form a net-shape, three-dimensional object. Our understanding of controlled drop deposition as applied to BPM is far from complete. Process parameters include the size and temperature of the liquid metal drop, its impact velocity and trajectory, and the condition and temperature of the substrate. Quantitative knowledge of the fluid mechanics and heat transfer of drop deposition and solidification are necessary to fully optimize the manufacturing process and to control the material microstructure of the final part. The object of this study is to examine the dynamics of liquid metal drops as they impinge upon a solid surface and solidify under conditions consistent with BPM (i.e. conditions which produce non-splatting drops). A program of both numerical simulations and experiments will be conducted. Questions this study will address include the following: How do the deformation and solidification of the drop depend on the properties of the fluid drop and the solid substrate? How does the presence of previously deposited drops affect the impingement and solidification process? How does the impingement of the new drop affect already deposited material? How does the cooling rate and solidification of the drops influence the material microstructure?

Description: The current study expands the application of computational fluid dynamics to three-dimensional multi-element high-lift systems by investigating the flow dynamics created by a slat edge. Flow is computed over a three-element high-lift configuration using an incompressible Navier-Stokes solver with structured, overset grids processed assuming full turbulence with the one-equation Baldwin-Barth turbulence model. The geometry consists of an unswept wing, which spans the wind tunnel test section, a single element half-span Fowler flap, and a three-quarter span slat. Results are presented for the wing configured for landing with a chord based Reynolds number of 3.7 million. Results for the three-quarter span slat case are compared to the full-span slat and two-dimensional investigations.

Description: This document is a manual describing how to use the Volume Grid Manipulation (VGM) software. The code is specifically designed to alter or manipulate existing surface and volume structured grids to improve grid quality through the reduction of grid line skewness, removal of negative volumes, and adaption of surface and volume grids to flow field gradients. The software uses a command language to perform all manipulations thereby offering the capability of executing multiple manipulations on a single grid during an execution of the code. The command language can be input to the VGM code by a UNIX style redirected file, or interactively while the code is executing. The manual consists of 14 sections. The first is an introduction to grid manipulation; where it is most applicable and where the strengths of such software can be utilized. The next two sections describe the memory management and the manipulation command language. The following 8 sections describe simple and complex manipulations that can be used in conjunction with one another to smooth, adapt, and reuse existing grids for various computations. These are accompanied by a tutorial section that describes how to use the commands and manipulations to solve actual grid generation problems. The last two sections are a command reference guide and trouble shooting sections to aid in the use of the code as well as describe problems associated with generated scripts for manipulation control.

Description: Jet noise is a major concern in the design of commercial aircraft. Studies by various researchers suggest that aerodynamic noise is a major contributor to jet noise. Some of these studies indicate that most of the aerodynamic jet noise due to turbulent mixing occurs when there is a rapid variation in turbulent structure, i.e. rapidly growing or decaying vortices. The objective of this research was to simulate a compressible round jet to study the non-linear evolution of vortices and the resulting acoustic radiations. In particular, to understand the effect of turbulence structure on the noise. An ideal technique to study this problem is Direct Numerical Simulations(DNS), because it provides precise control on the initial and boundary conditions that lead to the turbulent structures studied. It also provides complete 3-dimensional time dependent data. Since the dynamics of a temporally evolving jet are not greatly different from those, of a spatially evolving jet, a temporal jet problem was solved, using periodicity ill the direction of the jet axis. This enables the application of Fourier spectral methods in the streamwise direction. Physically this means that turbulent structures in the jet are repeated in successive downstream cells instead of being gradually modified downstream into a jet plume. The DNS jet simulation helps us understand the various turbulent scales and mechanisms of turbulence generation in the evolution of a compressible round jet. These accurate flow solutions will be used in future research to estimate near-field acoustic radiation by computing the total outward flux across a surface and determine how it is related to the evolution of the turbulent solutions. Furthermore, these simulations allow us to investigate the sensitivity of acoustic radiations to inlet/boundary conditions, with possible application to active noise suppression. In addition, the data generated can be used to compute various turbulence quantities such as mean velocities, turbulent stresses, etc. which will aid in turbulence modeling. This report will be presented in two chapters. The first chapter describes some work on the linear stability of a supersonic round jet and the implications of this for the jet noise problem. The second chapter is an extensive discussion of numerical work using the spectral method which we use to solve the compressible Navier-Stokes equations to study turbulent jet flows. The method uses Fourier expansions in the azimuthal and streamwise direction and a 1-D B-spline basis representation in the radial direction. The B-spline basis is locally supported and this ensures block diagonal matrix equations which can be solved in O(N) steps. This is a modification of a boundary layer code developed by Robert Moser. A very accurate highly resolved Direct Numerical Simulation (DNS) of a turbulent jet flow is produced.

Description: Jet noise is a major concern in the design of commercial aircraft. Studies by various researchers suggest that aerodynamic noise is a major contributor to jet noise. Some of these studies indicate that most of the aerodynamic jet noise due to turbulent mixing occurs when there is a rapid variation in turbulent structure, i.e. rapidly growing or decaying vortices. The objective of this research was to simulate a compressible round jet to study the non-linear evolution of vortices and the resulting acoustic radiations. In particular, to understand the effect of turbulence structure on the noise. An ideal technique to study this problem is Direct Numerical Simulations (DNS), because it provides precise control on the initial and boundary conditions that lead to the turbulent structures studied. It also provides complete 3-dimensional time dependent data. Since the dynamics of a temporally evolving jet are not greatly different from those of a spatially evolving jet, a temporal jet problem was solved, using periodicity in the direction of the jet axis. This enables the application of Fourier spectral methods in the streamwise direction. Physically this means that turbulent structures in the jet are repeated in successive downstream cells instead of being gradually modified downstream into a jet plume. The DNS jet simulation helps us understand the various turbulent scales and mechanisms of turbulence generation in the evolution of a compressible round jet. These accurate flow solutions will be used in future research to estimate near-field acoustic radiation by computing the total outward flux across a surface and determine how it is related to the evolution of the turbulent solutions. Furthermore, these simulations allow us to investigate the sensitivity of acoustic radiations to inlet/boundary conditions, with possible appli(,a- tion to active noise suppression. In addition, the data generated can be used to compute, various turbulence quantities such as mean velocities, turbulent stresses, etc. which will aid in turbulence modeling. This report will be presented in two chapters. The first chapter describes some work on the linear stability of a supersonic round jet and the implications of this for the jet noise problem. The second chapter is an extensive discussion of numerical work using the spectral method which we use to solve the compressible Navier-Stokes equations to study turbulent jet flows. The method uses Fourier expansions in the azimuthal and streamwise direction and a 1-D B-spline basis representation in the radial direction. The B-spline basis is locally supported and this ensures block diagonal matrix equations which can be solved in O(N) steps. This is a modification of a boundary layer code developed by Robert Moser. A very accurate highly resolved DNS of a turbulent jet flow is produced.

Description: A numerical method for the convective heat transfer problem is developed for low speed flow at mild temperatures. A simplified energy equation is added to the incompressible Navier-Stokes formulation by using Boussinesq approximation to account for the buoyancy force. A pseudocompressibility method is used to solve the resulting set of equations for steady-state solutions in conjunction with an approximate factorization scheme. A Neumann-type pressure boundary condition is devised to account for the interaction between pressure and temperature terms, especially near a heated or cooled solid boundary. It is shown that the present method is capable of predicting the temperature field in an incompressible flow.

Description: Boundary element algorithms for the solution of steady-state and transient heat conduction are presented. The algorithms are designed for efficient coupling with computational fluid dynamic discretizations and feature piecewise linear elements with offset nodal points. The steady-state algorithm employs the fundamental solution approach; the integration kernels are computed analytically based on linear shape functions, linear elements, and variably offset nodal points. The analytic expressions for both singular and nonsingular integrands are presented. The transient algorithm employs the transient fundamental solution; the temporal integration is performed analytically and the nonsingular spatial integration is performed numerically using Gaussian quadrature. A series solution to the integration is derived for the instance of a singular integrand. The boundary-only character of the algorithm is maintained by integrating the influence coefficients from initial time. Numerical results are compared to analytical solutions to verify the current boundary element algorithms. The steady-state and transient algorithms are numerically shown to be second-order accurate in space and time, respectively.

Description: A numerical method based on b-spline polynomials was developed to study incompressible flows in cylindrical geometries. A b-spline method has the advantages of possessing spectral accuracy and the flexibility of standard finite element methods. Using this method it was possible to ensure regularity of the solution near the origin, i.e. smoothness and boundedness. Because b-splines have compact support, it is also possible to remove b-splines near the center to alleviate the constraint placed on the time step by an overly fine grid. Using the natural periodicity in the azimuthal direction and approximating the streamwise direction as periodic, so-called time evolving flow, greatly reduced the cost and complexity of the computations. A direct numerical simulation of pipe flow was carried out using the method described above at a Reynolds number of 5600 based on diameter and bulk velocity. General knowledge of pipe flow and the availability of experimental measurements make pipe flow the ideal test case with which to validate the numerical method. Results indicated that high flatness levels of the radial component of velocity in the near wall region are physical; regions of high radial velocity were detected and appear to be related to high speed streaks in the boundary layer. Budgets of Reynolds stress transport equations showed close similarity with those of channel flow. However contrary to channel flow, the log layer of pipe flow is not homogeneous for the present Reynolds number. A topological method based on a classification of the invariants of the velocity gradient tensor was used. Plotting iso-surfaces of the discriminant of the invariants proved to be a good method for identifying vortical eddies in the flow field.

Description: A number of wavelet-based techniques for the analysis of experimental data are developed and illustrated. A multiscale analysis based on the Mexican hat wavelet is demonstrated as a tool for acquiring physical and quantitative information not obtainable by standard signal analysis methods. Experimental data for the analysis came from simultaneous hot-wire velocity traces in a bypass transition of the boundary layer on a heated flat plate. A pair of traces (two components of velocity) at one location was excerpted. A number of ensemble and conditional statistics related to dominant time scales for energy and momentum transport were calculated. The analysis revealed a lack of energy-dominant time scales inside turbulent spots but identified transport-dominant scales inside spots that account for the largest part of the Reynolds stress. Momentum transport was much more intermittent than were energetic fluctuations. This work is the first step in a continuing study of the spatial evolution of these scale-related statistics, the goal being to apply the multiscale analysis results to improve the modeling of transitional and turbulent industrial flows.

Description: Loop heat pipes (LHPs) can transport very large thermal power loads over long distances, through flexible, small diameter tubes against gravitational heads. In order to overcome the evaporator limit of LHPs, which is of about 0.07 MW/sq m, work was carried out to improve the efficiency by threefold to tenfold. The vapor passage geometry for the high heat flux conditions is shown. A bidisperse wick material within the circumferential vapor passages was used. Along with heat flux enhancement, several underlying issues were demonstrated, including the fabrication of bidisperse powder with controlled properties and the fabrication of a device geometry capable of replacing vapor passages with bidisperse powder.

Description: A multiblock approach is presented for solving two-dimensional incompressible turbulent flows on unstructured grids. The artificial compressibility form of the governing equations is solved by a vertex-centered, finite-volume implicit scheme which uses a backward Euler time discretization. Point Gauss-Seidel relaxations are used to solve the linear system of equations at each time step. This work introduces a multiblock strategy to the solution procedure, which greatly improves the efficiency of the algorithm by significantly reducing the memory requirements while not increasing the CPU time. Results presented in this work shows that the current multiblock algorithm requires 70% less memory than the single block algorithm.

Description: The final report consists of three papers which outline and demonstrate the new method for determining transition onset. The procedure developed under this grant requires specification of the instability mechanism, i.e., Tollmien-Schlichting or crossflow, that leads to transition. The attached papers are entitled: 'An Alternative to the e(sup n) Method for Determining Onset of Transition', 'Transition Model for Swept Wing Flows', and 'A Transition Closure Model for Predicting Transition Onset'.

Description: The authors review the development of diffuse-interface models of hydrodynamics and their application to a wide variety of interfacial phenomena. The authors discuss the issues involved in formulating diffuse-interface models for single-component and binary fluids. Recent applications and computations using these models are discussed in each case. Further, the authors address issues including sharp-interface analyses that relate these models to the classical free-boundary problem, related computational approaches to describe interfacial phenomena, and related approaches describing fully-miscible fluids.

Description: Previous work at NASA LeRC has shown that flow distortions in aircraft engine inlet ducts can be significantly reduced by mounting vortex generators, or small wing sections, on the inside surface of the engine inlet. The placement of the vortex generators is an important factor in obtaining the optimal effect over a wide operating envelope. In this regard, the only alternative to a long and expensive test program which would search out this optimal configuration is a good prediction procedure which could narrow the field of search. Such a procedure has been developed in collaboration with NASA LeRC, and results obtained by NASA personnel indicate that it shows considerable promise for predicting the viscous turbulent flow in engine inlet ducts in the presence of vortex generators. The prediction tool is a computer code which numerically solves the reduced Navier-Stokes equations and so is commonly referred to as RNS3D. Obvious deficiencies in RNS3D have been addressed in previous work. Primarily, it is known that the predictions of the mean velocity field of a turbulent boundary layer flow approaching separation are not in good agreement with data. It was suggested that the use of an algebraic mixing-length turbulence model in RNS3D is at least partly to blame for this. Additionally, the current turbulence model includes an assumption of isotropy which will ultimately fail to capture turbulence-driven secondary flow known to exist in noncircular ducts.

Description: The multiblock reacting Navier-Stokes flow solver RPLUS2D was modified for parallel implementation. Results for non-reacting flow calculations of this code indicate parallelization efficiencies greater than 84% are possible for a typical test problem. Results tend to improve as the size of the problem increases. The convergence rate of the scheme is degraded slightly when additional artificial block boundaries are included for the purpose of parallelization. However, this degradation virtually disappears if the solution is converged near to machine zero. Recommendations are made for further code improvements to increase efficiency, correct bugs in the original version, and study decomposition effectiveness.

Description: The results of an asymptotic theory for statistical closures for compressible turbulence are explored and validated with the direct numerical simulation of the isotropic decay and the homogeneous shear. An excellent collapse of the data is seen. The slow portion is found to scale, as predicted by the theory, with the quantity M(sub t)(sup 2) and epsilon(sub s). The rapid portion has an unambiguous scaling with alpha(sup 2)M(sub t)(sup s)epsilon(sub s)[P(sub k)/epsilon - l](Sk/epsilon)(sup 2). Implicit in the scaling is a dependence, as has been noted by others, on the gradient Mach number. A new feature of the effects of compressibility, that of the Kolmogorov scaling coefficient, alpha, is discussed. It is suggested that alpha may contain flow specific physics associated with large scales that might provide further insight into the structural effects of compressibility.

Description: This is a user's manual of the CMOTT turbulence module, version 2.0, developed for the NPARC code. The module is written in a self-contained manner so that the user can use any turbulence model in the module without concern as to how it is implemented and solved. Three two-equation turbulence models have been built into the module: Chien, Shih-Lumley and CMOTT models, and all of them have both the low Reynolds number and wall function options. Unlike Chien's model, both the Shih-Lumley and CMOTT models do not involve the dimensionless wall distance y(sup +) in the low Reynolds number approach, an advantage for separated flow calculations. The Van Driest transformation is used so that the wall functions can be applied to both incompressible and compressible flows. The manual gives the details of the turbulence models used and their numerical implementation. It also gives two application examples, one for subsonic and the other for transonic flow, for demonstration. The module can be easily linked to the NPARC code for practical applications.

Description: We present a detailed analysis of a recently proposed perfectly matched layer (PML) method for the absorption of acoustic waves. The split set of equations is shown to be only weakly well-posed, and ill-posed under small low order perturbations. This analysis provides the explanation for the stability problems associated with the split field formulation and illustrates why applying a filter has a stabilizing effect. Utilizing recent results obtained within the context of electromagnetics, we develop strongly well-posed absorbing layers for the linearized Euler equations. The schemes are shown to be perfectly absorbing independent of frequency and angle of incidence of the wave in the case of a non-convecting mean flow. In the general case of a convecting mean flow, a number of techniques is combined to obtain a absorbing layers exhibiting PML-like behavior. The efficacy of the proposed absorbing layers is illustrated though computation of benchmark problems in aero-acoustics.

Description: In this paper, we establish the causality between the model coefficients in the standard pressure-strain correlation model and the predicted equilibrium states for homogeneous turbulence. We accomplish this by performing a comprehensive fixed point analysis of the modeled Reynolds stress and dissipation rate equations. The results from this analysis will be very useful for developing improved pressure-strain correlation models to yield observed equilibrium behavior.

Description: This report details the development of a new two-equation turbulence closure model based on the exact turbulent kinetic energy k and the variance of vorticity, zeta. The model, which is applicable to three dimensional flowfields, employs one set of model constants and does not use damping or wall functions, or geometric factors.

Description: We proposed to create a single computational code incorporating methods that can model both rarefied and continuum flow to enable the efficient simulation of flow about space craft and high altitude hypersonic aerospace vehicles. The code was to use a single grid structure that permits a smooth transition between the continuum and rarefied portions of the flow. Developing an appropriate computational boundary between the two regions represented a major challenge. The primary approach chosen involves coupling a four-speed Lattice Boltzmann model for the continuum flow with the DSMC method in the rarefied regime. We also explored the possibility of using a standard finite difference Navier Stokes solver for the continuum flow. With the resulting code we will ultimately investigate three-dimensional plume impingement effects, a subject of critical importance to NASA and related to the work of Drs. Forrest Lumpkin, Steve Fitzgerald and Jay Le Beau at Johnson Space Center. Below is a brief background on the project and a summary of the results as of the end of the grant.

Description: With the aid of near-wall asymptotic analysis and results of direct numerical simulation, a new near-wall Reynolds stress model (NNWRS) is formulated based on the SSG high-Reynolds-stress model with wall-independent near-wall corrections. Only one damping function is used for flows with a wide range of Reynolds numbers to ensure that the near-wall modifications diminish away from the walls. The model is able to reproduce complicated flow phenomena induced by complex geometry, such as flow recirculation, reattachment and boundary-layer redevelopment in backward-facing step flow and secondary flow in three-dimensional square duct flow. In simple flows, including fully developed channel/pipe flow, Couette flow and boundary-layer flow, the wall effects are dominant, and the NNWRS model predicts less degree of turbulent anisotropy in the near-wall region compared with a wall-dependent near-wall Reynolds Stress model (NWRS) developed by So and colleagues. The comparison of the predictions given by the two models rectifies the misconception that the overshooting of skin friction coefficient in backward-facing step flow prevalent in those near-wall, models with wall normal is caused by he use of wall normal.

Description: Asymptotic regimes of geophysical dynamics are described for different Burger number limits. Rotating Boussinesq equations are analyzed in the asymptotic limit, of strong stratification in the Burger number of order one situation as well as in the asymptotic regime of strong stratification and weak rotation. It is shown that in both regimes horizontally averaged buoyancy variable is an adiabatic invariant for the full Boussinesq system. Spectral phase shift corrections to the buoyancy time scale associated with vertical shearing of this invariant are deduced. Statistical dephasing effects induced by turbulent processes on inertial-gravity waves are evidenced. The 'split' of the energy transfer of the vortical and the wave components is established in the Craya-Herring cyclic basis. As the Burger number increases from zero to infinity, we demonstrate gradual unfreezing of energy cascades for ageostrophic dynamics. The energy spectrum and the anisotropic spectral eddy viscosity are deduced with an explicit dependence on the anisotropic rotation/stratification time scale which depends on the vertical aspect ratio parameter. Intermediate asymptotic regime corresponding to strong stratification and weak rotation is analyzed where the effects of weak rotation are accounted for by an asymptotic expansion with full control (saturation) of vertical shearing. The regularizing effect of weak rotation differs from regularizations based on vertical viscosity. Two scalar prognostic equations for ageostrophic components (divergent velocity potential and geostrophic departure ) are obtained.

Description: A versatile and efficient multi-block method is presented for the simulation of both steady and unsteady flow, as well as aerodynamic design optimization of complete aircraft configurations. The compressible Euler and Reynolds Averaged Navier-Stokes (RANS) equations are discretized using a high resolution scheme on body-fitted structured meshes. An efficient multigrid implicit scheme is implemented for time-accurate flow calculations. Optimum aerodynamic shape design is achieved at very low cost using an adjoint formulation. The method is implemented on parallel computing systems using the MPI message passing interface standard to ensure portability. The results demonstrate that, by combining highly efficient algorithms with parallel computing, it is possible to perform detailed steady and unsteady analysis as well as automatic design for complex configurations using the present generation of parallel computers.

Description: Three-dimensional turbulent jets with rectangular cross-section are simulated with a finite-difference numerical method. The full Navier-Stokes equations are solved at low Reynolds numbers, whereas at the high Reynolds numbers filtered forms of the equations are solved along with a sub-grid scale model to approximate effects of the unresolved scales. A 2-N storage, third-order Runge-Kutta scheme is used for temporal discretization and a fourth-order compact scheme is used for spatial discretization. Computations are performed for different inlet conditions which represent different types of jet forcing. The phenomenon of axis-switching is observed, and it is confirmed that this is based on self-induction of the vorticity field. Budgets of the mean streamwise velocity show that convection is balanced by gradients of the Reynolds stresses and the pressure.

Description: A series of experiments were performed to investigate the effects of Mach number variation on the characteristics of the unsteady shock wave/turbulent boundary layer interaction generated by a blunt fin. A single blunt fin hemicylindrical leading edge diameter size was used in all of the experiments which covered the Mach number range from 2.0 to 5.0. The measurements in this investigation included surface flow visualization, static and dynamic pressure measurements, both on centerline and off-centerline of the blunt fin axis. Surface flow visualization and static pressure measurements showed that the spatial extent of the shock wave/turbulent boundary layer interaction increased with increasing Mach number. The maximum static pressure, normalized by the incoming static pressure, measured at the peak location in the separated flow region ahead of the blunt fin was found to increase with increasing Mach number. The mean and standard deviations of the fluctuating pressure signals from the dynamic pressure transducers were found to collapse to self-similar distributions as a function of the distance perpendicular to the separation line. The standard deviation of the pressure signals showed initial peaked distribution, with the maximum standard deviation point corresponding to the location of the separation line at Mach number 3.0 to 5.0. At Mach 2.0 the maximum standard deviation point was found to occur significantly upstream of the separation line. The intermittency distributions of the separation shock wave motion were found to be self-similar profiles for all Mach numbers. The intermittent region length was found to increase with Mach number and decrease with interaction sweepback angle. For Mach numbers 3.0 to 5.0 the separation line was found to correspond to high intermittencies or equivalently to the downstream locus of the separation shock wave motion. The Mach 2.0 tests, however, showed that the intermittent region occurs significantly upstream of the separation line. Power spectral densities measured in the intermittent regions were found to have self-similar frequency distributions when compared as functions of a Strouhal number for all Mach numbers and interaction sweepback angles. The maximum zero-crossing frequencies were found to correspond with the peak frequencies in the power spectra measured in the intermittent region.

Description: Coupling of sloshing dynamics within a partially filled rotating dewar of superfluid helium 2 with spacecraft dynamics are investigated in response to the realistic environmental disturbance forces and torques acting on the spacecraft during normal operation. This study investigates: (1) the rotating bubble of superfluid helium 2 reacting to combined environmental disturbances, including gravity gradient, aerodynamic, and magnetic forces and torques; (2) characteristics of slosh reaction forces and torques coupling with spacecraft dynamics; (3) the contribution of slosh dynamics to over-all spacecraft dynamics; and (4) activating of attitude and translation control system. The numerical computation of sloshing dynamics is based on the rotational frame, while the spacecraft dynamics is associated with non-rotational frame. Results show that the contributions of spacecraft dynamics are driven by the environmental disturbances coupling with slosh dynamics. Without considering the effects of environmental disturbances-driven slosh dynamics acting on spacecraft coupling with the spacecraft dynamics may lead to the wrong results for the development of spacecraft system guidance and attitude control techniques.

Description: A series of electrically heated tube experiments was performed to investigate the effect of high aspect ratio on curvature heat transfer enhancement in uniformly heated rectangular cooling passages. Three hardware geometries were tested: a baseline straight aspect ratio 10 tube, an aspect ratio 1 (square) tube with a 45 deg. curve, and an aspect ratio 10 tube with a 45 deg. curve. Gaseous nitrogen with the following properties was used as the coolant: ambient inlet temperature, pressures to 8.3 MPa, wall-to-bulk temperature ratios less than two, and Reynolds numbers based on hydraulic diameter ranging from 250,000 to 1,600,000. The measured curvature enhancement factors were compared to values predicted by three previously published models which had been developed for low aspect ratio tubes. The models were shown to be valid for the high aspect ratio tube as well the low aspect ratio tube, indicating that aspect ratio had little impact on the curvature heat transfer enhancement in these tests.

Description: The primary objective of this work is to provide accurate numerical solutions for selected flow fields and to compare and evaluate the performance of selected turbulence models with experimental results. Four popular turbulence models have been tested and validated against experimental data often turbulent flows. The models are: (1) the two-equation k-epsilon model of Wilcox, (2) the two-equation k-epsilon model of Launder and Sharma, (3) the two-equation k-omega/k-epsilon SST model of Menter, and (4) the one-equation model of Spalart and Allmaras. The flows investigated are five free shear flows consisting of a mixing layer, a round jet, a plane jet, a plane wake, and a compressible mixing layer; and five boundary layer flows consisting of an incompressible flat plate, a Mach 5 adiabatic flat plate, a separated boundary layer, an axisymmetric shock-wave/boundary layer interaction, and an RAE 2822 transonic airfoil. The experimental data for these flows are well established and have been extensively used in model developments. The results are shown in the following four sections: Part A describes the equations of motion and boundary conditions; Part B describes the model equations, constants, parameters, boundary conditions, and numerical implementation; and Parts C and D describe the experimental data and the performance of the models in the free-shear flows and the boundary layer flows, respectively.

Description: Method and apparatus for making measurements on fluids related to their complex permeability are disclosed. A microwave probe is provided for exposure to the fluids. The probe can be non-intrusive or can also be positioned at the location where measurements are to be made. The impedance of the probe is determined. in part. by the complex dielectric constant of the fluids at the probe. A radio frequency signal is transmitted to the probe and the reflected signal is phase and amplitude detected at a rapid rate for the purpose of identifying the fluids. Multiple probes may be selectively positioned to monitor the behavior of the fluids including their flow rate. Fluids may be identified as between two or more different fluids as well as multiple phases of the same fluid based on differences between their complex permittivities.

Description: The transpiration cooling method was considered for a scram-jet engine to accommodate thermally the situation where a very high heat flux (200 Btu/sq. ft sec) from hydrogen fuel combustion process is imposed to the engine walls. In a scram-jet engine, a small portion of hydrogen fuel passes through the porous walls of the engine combustor to cool the engine walls and at the same time the rest passes along combustion chamber walls and is preheated. Such a regenerative system promises simultaneously cooling of engine combustor and preheating the cryogenic fuel. In the experiment, an optical heating method was used to provide a heat flux of 200 Btu/sq. ft sec to the cylindrical surface of a porous stainless steel specimen which carried helium gas. The cooling efficiencies by transpiration were studied for specimens with various porosity. The experiments of various test specimens under high heat flux have revealed a phenomenon that chokes the medium flow when passing through a porous structure. This research includes the analysis of the system and a scaling conversion study that interprets the results from helium into the case when hydrogen medium is used.

Description: The analytic expression of the time evolution of the Reynolds stress anisotropy tensor in all planar homogeneous flows is obtained by exact integration of the modeled differential Reynolds stress equations. The procedure is based on results of tensor representation theory, is applicable for general pressure-strain correlation tensors, and can account for any additional turbulence anisotropy effects included in the closure. An explicit solution of the resulting system of scalar ordinary differential equations is obtained for the case of a linear pressure-strain correlation tensor. The properties of this solution are discussed, and the dynamic behavior of the Reynolds stresses is studied, including limit cycles and sensitivity to initial anisotropies.

Description: Work is devoted to study of free-stream vorticity normal to leading edge interaction with boundary layer over plate and resulting flow distortion influence on laminar-turbulent transition. In experiments made the wake behind the vertically stretched wire was used as a source of vortical disturbances and its effect on the boundary layer over the horizontally mounted plate with various leading edge shapes was investigated. The purpose of experiments was to check the predictions of theoretical works of M.E. Goldstein, et. al. This theory shows that small free-stream inhomogeneity interacting with leading edge produces considerable distortion of boundary layer flow. In general, results obtained confirms predictions of Goldstein's theory, i.e., the amplification of steady vortical disturbances in boundary layer caused by vortex lines stretching was observed. Experimental results fully coincide with predictions of theory for large Reynolds number, relatively sharp leading edge and small disturbances. For large enough disturbances the flow distortion caused by symmetric wake unexpectedly becomes antisymmetric in spanwise direction. If the leading edge is too blunt the maximal distortion takes place immediately at the nose and no further amplification was observed. All these conditions and results are beyond the scope of Goldstein's theory.

Description: The dissipation rate transport equation remains the most uncertain part of turbulence modeling. The difficulties arc increased when external agencies like rotation prevent straightforward dimensional analysis from determining the correct form of the modelled equation. In this work, the dissipation rate transport equation and subgrid scale models for rotating turbulence are derived from an analytical statistical theory of rotating turbulence. In the strong rotation limit, the theory predicts a turbulent steady state in which the inertial range energy spectrum scales as k(sup -2) and the turbulent time scale is the inverse rotation rate. This scaling has been derived previously by heuristic arguments.

Description: Aircraft engine combustors generally involve turbulent swirling flows in order to enhance fuel-air mixing and flame stabilization. It has long been recognized that eddy viscosity turbulence models are unable to appropriately model swirling flows. Therefore, it has been suggested that, for the modeling of these flows, a second order closure scheme should be considered because of its ability in the modeling of rotational and curvature effects. However, this scheme will require solution of many complicated second moment transport equations (six Reynolds stresses plus other scalar fluxes and variances), which is a difficult task for any CFD implementations. Also, this scheme will require a large amount of computer resources for a general combustor swirling flow. This report is devoted to the development of a cubic Reynolds stress-strain model for turbulent swirling flows, and was inspired by the work of Launder's group at UMIST. Using this type of model, one only needs to solve two turbulence equations, one for the turbulent kinetic energy k and the other for the dissipation rate epsilon. The cubic model developed in this report is based on a general Reynolds stress-strain relationship. Two flows have been chosen for model evaluation. One is a fully developed rotating pipe flow, and the other is a more complex flow with swirl and recirculation.

Description: A new post-processing technique has been developed, based on the Intrinsic Low Dimensional Manifold (ILDM) method of Maas and Pope. The ILDM method is a dynamical systems approach to the simplification of large chemical kinetic mechanisms. By identifying low-dimensional attracting manifolds, the method allows complex full mechanisms to be parameterized by just a few variables: In effect, generating reduced chemical mechanisms by an automatic procedure. These resulting mechanisms however, still retain all the species used in the full mechanism. The NO(x) post-processor takes an ILDM reduced mechanism and attempts to map this mechanism to the results of a CFD calculation. This mapping allows the NO(x) concentrations at each grid node to be obtained from the ILDM reduced mechanism, as well as other trace species of interest. Because a mapping procedure is used, this method is very fast, being able to process one million node calculations in just a few minutes.

Description: Radiative heat transfer through fibrous insulation used in thermal protection systems (TPS) is significant at high temperatures (1200 C). Decreasing the radiative heat transfer through the fibrous insulation can thus have a major impact on the insulating ability of the TPS. Reflective coatings applied directly to the individual fibers in fibrous insulation should decrease the radiative heat transfer leading to an insulation with decreased effective thermal conductivity. Coatings with high infrared reflectance have been developed using sol-gel techniques. Using this technique, uniform coatings can be applied to fibrous insulation without an appreciable increase in insulation weight or density. Scanning electron microscopy, Fourier Transform infrared spectroscopy, and ellipsometry have been performed to evaluate coating performance.

Description: The aeroheating characteristics of the X-38 Revision 3.1 lifting-body configuration have been experimentally examined in the Langley 20-inch Mach 6 Tunnel. Global surface heat transfer distributions, surface streamline patterns, and shock shapes were measured on a 0.0362-scale model of a proposed Space Station Crew Return Vehicle at Mach 6 in air. Parametric variations include angles-of-attack of 20 deg, 30 deg, and 40 deg; Reynolds numbers based on model length of 0.9 to 3.7 million; and body-flap deflections of O deg, 20 deg, 25 deg, and 30 deg. The effects of discrete roughness elements, which included trip height, location, size, and orientation, as well as multiple-trip parametrics, were investigated. This document is intended to serve as a quick release of preliminary data to the X-38 program; analysis is limited to observations of the experimental trends in order to expedite dissemination.

Description: This report discusses some analytical procedures to enhance the real time solutions of PMARC matrices applicable to the Wall Interference Correction Scheme (WICS) currently being implemented at the 12 foot Pressure Tunell. WICS calculations involve solving large linear systems in a reasonably speedy manner necessitating exploring further improvement in solution time. This paper therefore presents some of the associated theory of the solution of linear systems. Then it discusses a geometrical interpretation of the residual correction schemes. Finally, some results of the current investigation are presented.

Description: The global nonlinear behavior of finite discretizations for constant time steps and fixed or adaptive grid spacings is studied using tools from dynamical systems theory. Detailed analysis of commonly used temporal and spatial discretizations for simple model problems is presented. The role of dynamics in the understanding of long time behavior of numerical integration and the nonlinear stability, convergence, and reliability of using time-marching approaches for obtaining steady-state numerical solutions in computational fluid dynamics (CFD) is explored. The study is complemented with examples of spurious behavior observed in steady and unsteady CFD computations. The CFD examples were chosen to illustrate non-apparent spurious behavior that was difficult to detect without extensive grid and temporal refinement studies and some knowledge from dynamical systems theory. Studies revealed the various possible dangers of misinterpreting numerical simulation of realistic complex flows that are constrained by available computing power. In large scale computations where the physics of the problem under study is not well understood and numerical simulations are the only viable means of solution, extreme care must be taken in both computation and interpretation of the numerical data. The goal of this paper is to explore the important role that dynamical systems theory can play in the understanding of the global nonlinear behavior of numerical algorithms and to aid the identification of the sources of numerical uncertainties in CFD.

Description: Three types of turbulence models which account for rotational effects in noninertial frames of reference are evaluated for the case of incompressible, fully developed rotating turbulent channel flow. The different types of models are a Coriolis-modified eddy-viscosity model, a realizable algebraic stress model, and an algebraic stress model which accounts for dissipation rate anisotropies. A direct numerical simulation of a rotating channel flow is used for the turbulent model validation. This simulation differs from previous studies in that significantly higher rotation numbers are investigated. Flows at these higher rotation numbers are characterized by a relaminarization on the cyclonic or suction side of the channel, and a linear velocity profile on the anticyclonic or pressure side of the channel. The predictive performance of the three types of models are examined in detail, and formulation deficiencies are identified which cause poor predictive performance for some of the models. Criteria are identified which allow for accurate prediction of such flows by algebraic stress models and their corresponding Reynolds stress formulations.

Description: Theories of turbulent time correlations are applied to compute frequency spectra of sound radiated by isotropic turbulence and by turbulent shear flows. The hypothesis that Eulerian time correlations are dominated by the sweeping action of the most energetic scales implies that the frequency spectrum of the sound radiated by isotropic turbulence scales as omega(exp 4) for low frequencies and as omega(exp -3/4) for high frequencies. The sweeping hypothesis is applied to an approximate theory of jet noise. The high frequency noise again scales as omega(exp -3/4), but the low frequency spectrum scales as omega(exp 2). In comparison, a classical theory of jet noise based on dimensional analysis gives omega(exp -2) and omega(exp 2) scaling for these frequency ranges. It is shown that the omega(exp -2) scaling is obtained by simplifying the description of turbulent time correlations. An approximate theory of the effect of shear on turbulent time correlations is developed and applied to the frequency spectrum of sound radiated by shear turbulence. The predicted steepening of the shear dominated spectrum appears to be consistent with jet noise measurements.

Description: Direct numerical simulations have been used to examine the effect of the initial disturbance field on the development of three-dimensionality and the transition to turbulence in the incompressible plane wake. The simulations were performed using a new numerical method for solving the time-dependent, three-dimensional, incompressible Navier-Stokes equations in flows with one infinite and two periodic directions. The method uses standard Fast Fourier Transforms and is applicable to cases where the vorticity field is compact in the infinite direction. Initial disturbances fields examined were combinations of two-dimensional waves and symmetric pairs of 60 deg oblique waves at the fundamental, subharmonic, and sub-subharmonic wavelengths. The results of these simulations indicate that the presence of 60 deg disturbances at the subharmonic streamwise wavelength results in the development of strong coherent three-dimensional structures. The resulting strong three-dimensional rate-of-strain triggers the growth of intense fine scale motions. Wakes initiated with 60 deg disturbances at the fundamental streamwise wavelength develop weak coherent streamwise structures, and do not develop significant fine scale motions, even at high Reynolds numbers. The wakes which develop strong three-dimensional structures exhibit growth rates on par with experimentally observed turbulent plane wakes. Wakes which develop only weak three-dimensional structures exhibit significantly lower late time growth rates. Preliminary studies of wakes initiated with an oblique fundamental and a two-dimensional subharmonic, which develop asymmetric coherent oblique structures at the subharmonic wavelength, indicate that significant fine scale motions only develop if the resulting oblique structures are above an angle of approximately 45 deg.

Description: A review of recent literature has been made in the areas of unsteady measurement in turbomachinery and LP turbine design. This has included a detailed review of methods for extracting quantitative shear stress data from surface mounted hot-film arrays. Research goals were identified on the basis of the literature survey and following interaction with the GE and NASA contracts. Progress has also been made in the detailed planning of the research program. An array of hot-film sensors has been mounted at mid-span on a second stage nozzle blade ready for installation in the current turbine buildup. A new technique of applying the hot-film array using double sided adhesive is being tested. The temperature in the vicinity of the hot-films will be monitored using a thermocouple attached to the blade pressure surface. This will increase the accuracy of the data reduction process compared with the methods applied by Halstead et at. (1995) and others. Data obtained from this first set of sensors, along with data from the hot-wire investigation will help to determine the best locations for surface hot-film instrumentation for subsequent builds of the turbine. Other LSRT instrumentation issues which have been addressed include calibration of the torque measurement system and procurement of custom design hot-wire sensors required to compensate for the slope flow path. Preliminary attempts at modeling the boundary layer transition process have been completed. A steady, coupled viscous-inviscid flow solver (the MISES code of Drela and Youngren) was used to calculate the steady-flow boundary layer transition location.

Description: An experimental investigation of the non-reactive mixing processes associated with a lobed fuel injector in a coflowing air stream is presented. The lobed fuel injector is a device which generates streamwise vorticity, producing high strain rates which can enhance the mixing of reactants while delaying ignition in a controlled manner. The lobed injectors examined in the present study consist of two corrugated plates between which a fuel surrogate, CO2, is injected into coflowing air. Acetone is seeded in the CO2 supply as a fuel marker. Comparison of two alternative lobed injector geometries is made with a straight fuel injector to determine net differences in mixing and strain fields due to streamwise vorticity generation. Planar laser-induced fluorescence (PLIF) of the seeded acetone yields two-dimensional images of the scalar concentration field at various downstream locations, from which local mixing and scalar dissipation rates are computed. It is found that the lobed injector geometry can enhance molecular mixing and create a highly strained flowfield, and that the strain rates generated by scalar energy dissipation can potentially delay ignition in a reacting flowfield.

Description: A three-dimensional Navier-Stokes code has been used to compute the heat transfer coefficient on two film-cooled turbine blades, namely, the VKI rotor with six rows of cooling holes, including three rows on the shower head and the C3X vane with nine rows of holes, including five rows on the shower head. Predictions of heat transfer coefficient at the blade surface using three two-equation turbulence model specifically, Coakley's q-omega model, Chien's k-epsilon model and Wilcox's k-omega model with Menter's modifications, have been compared with the experimental data of Camci and Arts for the VKI rotor, and of Hylton et al. for the C3X vane along with predictions using the Baldwin-Lomar (B-L) model taken from Garg and Gaugler. It is found that for the cases considered here the two equation models predict the blade heat transfer somewhat better than the B-L model except immediately downstream of the film-cooled holes on the suction surface of the VKI rotor, and over most of the suction surface of the C3X vane. However, all two-equation models require 40% more computer core than the B-L model for solution, and while the q-omega and k-epsilon models need 40% more computer time than the B-L model the k-omega model requires at least 65% more time because of the slower rate of convergence. It is found that the heat transfer coefficient exhibit a strong spanwise as well as streamwise variation for both blades and all turbulence models.

Description: We have found that it is possible to have steady thermocapillary flows in any of these three systems provided the liquid layer is not too thin or the surface tension is not too small. If these conditions are not met, then no steady solutions are possible and the film thins and approaches a condition of rupture near the hot end. Future work on these problems will include the effect of liquid inertia in the calculations for the cylindrical cavity. The stabilization effect of inertia may offset the destabilization effect inherent in the cylindrical geometry and result in a stable film flow on the cylinder. This type of behavior could be very significant in terms of understanding the stability of thermocapillary flow in liquid bridges and in the float-zone crystal growth process.

Description: Due to noise considerations, ultra high bypass ducted fans have become a more viable design. These ducted fans typically consist of a rotor stage containing a wide chord fan and a stator stage. One of the concerns for this design is the classical flutter that keeps occurring in various unducted fan blade designs. These flutter are catastrophic and are to be avoided in the flight envelope of the engine. Some numerical investigations by Williams, Cho and Dalton, have suggested that a duct around a propeller makes it more unstable. This needs to be further investigated. In order to design an engine to safely perform a set of desired tasks, accurate information of the stresses on the blade during the entire cycle of blade motion is required. This requirement in turn demands that accurate knowledge of steady and unsteady blade loading be available. Aerodynamic solvers based on unsteady three-dimensional analysis will provide accurate and fast solutions and are best suited for aeroelastic analysis. The Euler solvers capture significant physics of the flowfield and are reasonably fast. An aerodynamic solver Ref. based on Euler equations had been developed under a separate grant from NASA Lewis in the past. Under the current grant, this solver has been modified to calculate the aeroelastic characteristics of unducted and ducted rotors. Even though, the aeroelastic solver based on three-dimensional Euler equations is computationally efficient, it is still very expensive to investigate the effects of multiple stages on the aeroelastic characteristics. In order to investigate the effects of multiple stages, a two-dimensional multi stage aeroelastic solver was also developed under this task, in collaboration with Dr. T. S. R. Reddy of the University of Toledo. Both of these solvers were applied to several test cases and validated against experimental data, where available.

Description: This document contains papers presented at the Eighth Annual Thermal and Fluids Analysis Workshop (TFAWS) on Spacecraft Analysis and Design hosted by the NASA/Johnson Space Center (JSC) on September 8-11, 1997, and held at the University of Houston - Clear Lake (UHCL) in the Bayou Building. The Workshop was sponsored by NASA/JSC. Seminars were hosted and technical papers were provided in fluid and thermal dynamics. Seminars were given in GASP, SINDA, SINAPS Plus, TSS, and PHOENICS. Seventeen papers were presented.

Description: A series of uni-element rocket injector studies were completed to provide benchmark quality data needed to validate computational fluid dynamic models. A shear coaxial injector geometry was selected as the primary injector for study using gaseous hydrogen/oxygen and gaseous hydrogen/liquid oxygen propellants. Emphasis was placed on the use of nonintrusive diagnostic techniques to characterize the flowfields inside an optically-accessible rocket chamber. Measurements of the velocity and species fields were obtained using laser velocimetry and Raman spectroscopy, respectively. Qualitative flame shape information was also obtained using laser-induced fluorescence excited from OH radicals and laser light scattering studies of aluminum oxide particle seeded combusting flows. The gaseous hydrogen/liquid oxygen propellant studies for the shear coaxial injector focused on breakup mechanisms associated with the liquid oxygen jet under subcritical pressure conditions. Laser sheet illumination techniques were used to visualize the core region of the jet and a Phase Doppler Particle Analyzer was utilized for drop velocity, size and size distribution characterization. The results of these studies indicated that the shear coaxial geometry configuration was a relatively poor injector in terms of mixing. The oxygen core was observed to extend well downstream of the injector and a significant fraction of the mixing occurred in the near nozzle region where measurements were not possible to obtain. Detailed velocity and species measurements were obtained to allow CFD model validation and this set of benchmark data represents the most comprehensive data set available to date. As an extension of the investigation, a series of gas/gas injector studies were conducted in support of the X-33 Reusable Launch Vehicle program. A Gas/Gas Injector Technology team was formed consisting of the Marshall Space Flight Center, the NASA Lewis Research Center, Rocketdyne and Penn State. Injector geometries studied under this task included shear and swirl coaxial configurations as well as an impinging jet injector.

Description: Applications are described of high-performance computing methods to the numerical simulation of complete jet engines. The methodology focuses on the partitioned analysis of the interaction of the gas flow with a flexible structure and with the fluid mesh motion driven by structural displacements. The latter is treated by a ALE technique that models the fluid mesh motion as that of a fictitious mechanical network laid along the edges of near-field elements. New partitioned analysis procedures to treat this coupled three-component problem were developed. These procedures involved delayed corrections and subcycling, and have been successfully tested on several massively parallel computers, including the iPSC-860, Paragon XP/S and the IBM SP2. The NASA-sponsored ENG10 program was used for the global steady state analysis of the whole engine. This program uses a regular FV-multiblock-grid discretization in conjunction with circumferential averaging to include effects of blade forces, loss, combustor heat addition, blockage, bleeds and convective mixing. A load-balancing preprocessor for parallel versions of ENG10 was developed as well as the capability for the first full 3D aeroelastic simulation of a multirow engine stage. This capability was tested on the IBM SP2 parallel supercomputer at NASA Ames.

Description: This report describes a concurrent Euler flow solver for flows around complex 3-D bodies. The solver is based on a cell-centered finite volume methodology on 3-D unstructured tetrahedral grids. In this algorithm, spatial discretization for the inviscid convective term is accomplished using an upwind scheme. A localized reconstruction is done for flow variables which is second order accurate. Evolution in time is accomplished using an explicit three-stage Runge-Kutta method which has second order temporal accuracy. This is adapted for concurrent execution using another proven methodology based on concurrent graph abstraction. This solver operates on heterogeneous network architectures. These architectures may include a broad variety of UNIX workstations and PCs running Windows NT, symmetric multiprocessors and distributed-memory multi-computers. The unstructured grid is generated using commercial grid generation tools. The grid is automatically partitioned using a concurrent algorithm based on heat diffusion. This results in memory requirements that are inversely proportional to the number of processors. The solver uses automatic granularity control and resource management techniques both to balance load and communication requirements, and deal with differing memory constraints. These ideas are again based on heat diffusion. Results are subsequently combined for visualization and analysis using commercial CFD tools. Flow simulation results are demonstrated for a constant section wing at subsonic, transonic, and a supersonic case. These results are compared with experimental data and numerical results of other researchers. Performance results are under way for a variety of network topologies.

Description: A model of sound generated in a high subsonic (Mach 0.9) circular jet is solved numerically in cylindrical coordinates for nonaxisymmetric disturbances. The jet is excited by transient mass injection by a finite duration pulse via a modulated ring source. The nonaxisymmetric solution is computed for long times after the initial disturbance has exited the computational domain. The long time behavior of the jet is dominated by vorticity and pressure disturbances generated at the nozzle lip and growing as they convect down-stream in the jet. These disturbances generate sound as they propagate. The primary non-axisymmetric effect that we simulate is that of a flapping mode where regions of high and low pressure alternate on opposite sides of the jet. The predominant feature of this mode is the appearance of relatively large deviations of the pressure from the ambient pressure on opposite sides of the jet and the convection of these regions downstream. We illustrate flow field, near field and far field data. Important nonaxisymmetric characteristics of the near and flow field disturbances include roughly periodic pressure elevations and depressions at opposite values of the azimuthal angle psi. These correspond to pressure disturbances propagating in the axial direction. The azimuthal velocity exhibits a sinusoidal dependence on psi with similar roughly periodic disturbances. For every azimuthal angle psi, the jet radiation peaks about 30 deg. from the jet axis, however there is now a pronounced dependence of the far field radiation pattern on psi.

Description: This report describes the final results of the project, which are reported in the reference in Appendix A. This reference contains the technical details, and will constitute the main body of the report. The paper in Appendix A was first presented at the 8th International Symposium on Applications of Laser Techniques to Fluid Mechanics, July, 1996. This paper was one of 30 papers selected from more than 200 papers presented at the symposium for publication in a book by Springer-Verlag. It will become available during 1997.

Description: The pressure disturbances generated by an instability wave in the shear layer of a supersonic jet are studied for an axisymmetric jet inside a lined cylindrical duct. For the supersonic jet, locally linear stability analysis with duct wall boundary conditions is used to calculate the eigenvalues and the eigenfunctions. These values are used to determine the growth rates and phase velocities of the instability waves and the radial pressure disturbance patterns. The study is confined to the dominant Kelvin-Helmholtz instability mode and to the region just downstream of the nozzle exit where the shear layer is growing but is still small in size compared to the radius of the duct. Numerical results are used to study the effects of changes in the outer flow, growth in the shear layer thickness, wall distance, wall impedance, and frequency. Results indicate that the effects of the duct wall on shear layer growth rates diminish as the outer flow increases. Also, wall reflections cause variations in growth rates depending on wall height and Strouhal number. These variations are due to the phase relationship between the outgoing and the reflected incoming pressure disturbances at the shear layer. The growth rate variations can be reduced and the maximum growth rate minimized by keeping the imaginary part of the impedance negative.

Description: A linear instability model for multiple spatially periodic supersonic rectangular jets is solved using Floquet-Bloch theory. The disturbance environment is investigated using a two dimensional perturbation of a mean flow. For all cases large temporal growth rates are found. This work is motivated by an increase in mixing found in experimental measurements of spatially periodic supersonic rectangular jets with phase-locked screech. The results obtained in this paper suggests that phase-locked screech or edge tones may produce correlated spatially periodic jet flow downstream of the nozzles which creates a large span wise multi-nozzle region where a disturbance can propagate. The large temporal growth rates for eddies obtained by model calculation herein are related to the increased mixing since eddies are the primary mechanism that transfer energy from the mean flow to the large turbulent structures. Calculations of growth rates are presented for a range of Mach numbers and nozzle spacings corresponding to experimental test conditions where screech synchronized phase locking was observed. The model may be of significant scientific and engineering value in the quest to understand and construct supersonic mixer-ejector nozzles which provide increased mixing and reduced noise.

Description: We studied the functionality of the IFA300 Hot Wire Anemometry system and its feasibility for the calibration of Direct Estimator Velocimetry (DEV) using Laser Induced Florescence (LIF). The experimental setup consisted of a single-wire hot wire probe measuring one component of velocity, the anemometry computational hardware and software, an HP oscilloscope, and a table fan to produce a simple flow with fluctuations. Measurements were taken at several points in the stream wise and transverse directions in the flow and various parameters recorded such as mean velocities temperature, turbulence intensities, skewness coefficients and flatness coefficients. The IFA300 software also allowed us to perform statistical manipulations such as spectrum analysis on velocities samples and correlation. Utilization of data files, also produced by the anemometry software, and post analysis were performed to produce graphical representations of turbulent intensity versus probe position and a flow field velocity profile. We concluded that the IFA300 Hot Wire Anemometry system is a reliable and functional method for calibration of DEV using LIF. Our future intentions are to set up a test chamber such that both velocity measurement techniques can be applied simultaneously, thus the calibration.

Description: The potential for using a novel diamond-shaped nozzle which may allow for superior mixing characteristics of supersonic jets without significant thrust losses is explored. The results of flow visualization and pressure measurements indicate the presence of distinct structures in the shear layers, not normally observed in shear layers of axisymmetric and rectangular jets. As characteristics of these features suggests that they are a manifestation of significant streamwise vorticity in the shear layers. Despite the distinct nature of the flowfield structure of the present shear layer, the global growth rates of this shear layer were found to be very similar to its two-dimensional and axisymmetric counterparts. These and other observations suggest that the presence of streamwise vorticity may not play a significant role in the global development of a compressible shear layer.

Description: Reports by the research staff and graduate students of the Mechanical Engineering Department at the University of Queensland are collected and presented. These reports cover various studies related to the advancement of scramjet technology and the operation of advanced hypervelocity shock-expansion tubes. The report topics include the experimental studies of mixing and combustion in a scramjet flow path, the measurement of integrated thrust and skin friction, and the development of a free-piston-driven expansion tunnel capable of delivering a test gas at superorbital velocities.

Description: It has been shown that turbulence and temperature influence the droplet sizes in expansion fan induced condensation by studying the Rayleigh scattering from one port in our shock tube's test section. We have modified our set-up so as to allow, using two ports, the real time measurement of the influence of turbulence and temperature on the rate at which these droplets grow. To do this, we looked at the Rayleigh scattering from two different ports for ten Reynolds numbers at five different temperatures. We modeled the time of flight of droplets, using the equations of one-dimensional gas dynamics and the measured shock wave speed in shock tube's driven section.

Description: An experimental method is described to measure and compare the convective heat transfer coefficient of natural and simulated ice accretion roughness and to provide a rational means for determining accretion-related enhanced heat transfer coefficients. The natural ice accretion roughness was a sample casting made from accretions at the NASA Lewis Icing Research Tunnel (IRT). One of these castings was modeled using a Spectral Estimation Technique (SET) to produce three roughness elements patterns that simulate the actual accretion. All four samples were tested in a flat-plate boundary layer at angle of attack in a "dry" wind tunnel test. The convective heat transfer coefficient was measured using infrared thermography. It is shown that, dispite some problems in the current data set, the method does show considerable promise in determining roughness-induced heat transfer coefficients, and that, in addition to the roughness height and spacing in the flow direction, the concentration and spacing of elements in the spanwise direction are important parameters.

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 ambient air for successful operation in the ramjet and scramjet modes. 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. The modeling basis was centered on using convective Mach Number as the similarity parameter to establish correlation between subscale, cold flow tests and full scale, hot firing modes. This parameter has been used successfully to correlate supersonic shear layer growth rates. The experiment design includes hot (600 R) air as the rocket exhaust simulant and hot (760 R) carbon dioxide as the turbine exhaust gas simulant. The combination of gases and their elevated temperatures was required to achieve a convective Mach Number which matched the fall scale item design conditions. The carbon dioxide is seeded with Acetone to permit tracing of the mixing processes through Laser Induced Fluorescence (LIF) techniques. The experiment and its design will be discussed in detail. Both the rocket and turbine exhaust duct nozzles are of unique (square and rectangular) shape and the turbine exhaust e)dt intersects the rocket nozzle wall upstream of the exit. Cold flow testing with the individual nozzles has been conducted to ascertain their behavior in comparison to conventional flow theory. These data are presented.

Description: A fully threaded tree (FTT) for adaptive refinement of regular meshes is described. By using a tree threaded at all levels, tree traversals for finding nearest neighbors are avoided. All operations on a tree including tree modifications are O(N), where N is a number of cells, and are performed in parallel. An efficient implementation of the tree is described that requires 2N words of memory. A filtering algorithm for removing high frequency noise during mesh refinement is described. A FTT can be used in various numerical applications. In this paper, it is applied to the integration of the Euler equations of fluid dynamics. An adaptive mesh time stepping algorithm is described in which different time steps are used at different l evels of the tree. Time stepping and mesh refinement are interleaved to avoid extensive buffer layers of fine mesh which were otherwise required ahead of moving shocks. Test examples are presented, and the FTT performance is evaluated. The three dimensional simulation of the interaction of a shock wave and a spherical bubble is carried out that shows the development of azimuthal perturbations on the bubble surface.

Description: In the second year, we continued to built upon and improve our scanline-based direct volume renderer that we developed in the first year of this grant. This extremely general rendering approach can handle regular or irregular grids, including overlapping multiple grids, and polygon mesh surfaces. It runs in parallel on multi-processors. It can also be used in conjunction with a k-d tree hierarchy, where approximate models and error terms are stored in the nodes of the tree, and approximate fast renderings can be created. We have extended our software to handle time-varying data where the data changes but the grid does not. We are now working on extending it to handle more general time-varying data. We have also developed a new extension of our direct volume renderer that uses automatic decimation of the 3D grid, as opposed to an explicit hierarchy. We explored this alternative approach as being more appropriate for very large data sets, where the extra expense of a tree may be unacceptable. We also describe a new approach to direct volume rendering using hardware 3D textures and incorporates lighting effects. Volume rendering using hardware 3D textures is extremely fast, and machines capable of using this technique are becoming more moderately priced. While this technique, at present, is limited to use with regular grids, we are pursuing possible algorithms extending the approach to more general grid types. We have also begun to explore a new method for determining the accuracy of approximate models based on the light field method described at ACM SIGGRAPH '96. In our initial implementation, we automatically image the volume from 32 equi-distant positions on the surface of an enclosing tessellated sphere. We then calculate differences between these images under different conditions of volume approximation or decimation. We are studying whether this will give a quantitative measure of the effects of approximation. We have created new tools for exploring the differences between images produced by various rendering methods. Images created by our software can be stored in the SGI RGB format. Our idtools software reads in pair of images and compares them using various metrics. The differences of the images using the RGB, HSV, and HSL color models can be calculated and shown. We can also calculate the auto-correlation function and the Fourier transform of the image and image differences. We will explore how these image differences compare in order to find useful metrics for quantifying the success of various visualization approaches. In general, progress was consistent with our research plan for the second year of the grant.

Description: Using a two-point closure theory, the Eddy-Damped-Quasi-Normal-Markovian (EDQNM) approximation, we have investigated the energy transfer process and triadic interactions of compressible turbulence. In order to analyze the compressible mode directly, the Helmholtz decomposition is used. The following issues were addressed: (1) What is the mechanism of energy exchange between the solenoidal and compressible modes, and (2) Is there an energy cascade in the compressible energy transfer process? It is concluded that the compressible energy is transferred locally from the solenoidal part to the compressible part. It is also found that there is an energy cascade of the compressible mode for high turbulent Mach number (M(sub t) greater than or equal to 0.5). Since we assume that the compressibility is weak, the magnitude of the compressible (radiative or cascade) transfer is much smaller than that of solenoidal cascade. These results are further confirmed by studying the triadic energy transfer function, the most fundamental building block of the energy transfer.

Description: Reports by the research staff and graduate students of the Mechanical Engineering Department at the University of Queensland are collected and presented. These reports cover various studies related to the advancement of scramjet technology and the operation of advanced hypervelocity shock-expansion tubes. The report topics include the experimental studies of mixing and combustion in a scramjet flow path, the measurement of integrated thrust and skin friction, and the development of a free-piston-driven expansion tunnel capable delivering a test gas at super orbital velocities.

Description: Results of an experimental investigation are presented in which the use of porous and microporous honeycomb composite materials is evaluated as an alternate to perforated solid plates for boundary-layer bleed in supersonic aircraft inlets. The terms "porous" and "microporous," respectively, refer to bleed orifice diameters roughly equal to and much less than the displacement thickness of the approach boundary-layer. A Baseline porous solid plate, two porous honeycomb, and three microporous honeycomb configurations are evaluated. The performance of the plates is characterized by the flow coefficient and relative change in boundary-layer profile parameters across the bleed region. The tests were conducted at Mach numbers of 1.27 and 1.98. The results show the porous honeycomb is not as efficient at removing mass compared to the baseline. The microporous plates were about equal to the baseline with one plate demonstrating a significantly higher efficiency. The microporous plates produced significantly fuller boundary-layer profiles downstream of the bleed region for a given mass flow removal rate than either the baseline or the porous honeycomb plates.

Description: Parallel computation of unsteady flows requires significant computational resources. The utilization of a network of workstations seems an efficient solution to the problem where large problems can be treated at a reasonable cost. This approach requires the solution of several problems: 1) the partitioning and distribution of the problem over a network of workstation, 2) efficient communication tools, 3) managing the system efficiently for a given problem. Of course, there is the question of the efficiency of any given numerical algorithm to such a computing system. NPARC code was chosen as a sample for the application. For the explicit version of the NPARC code both two- and three-dimensional problems were studied. Again both steady and unsteady problems were investigated. The issues studied as a part of the research program were: 1) how to distribute the data between the workstations, 2) how to compute and how to communicate at each node efficiently, 3) how to balance the load distribution. In the following, a summary of these activities is presented. Details of the work have been presented and published as referenced.

Description: Distributed-memory parallel computers dominate today's parallel computing arena. These machines, such as Intel Paragon, IBM SP2, and Cray Origin2OO, have successfully delivered high performance computing power for solving some of the so-called "grand-challenge" problems. Despite initial success, parallel machines have not been widely accepted in production engineering environments due to the complexity of parallel programming. On a parallel computing system, a task has to be partitioned and distributed appropriately among processors to reduce communication cost and to attain load balance. More importantly, even with careful partitioning and mapping, the performance of an algorithm may still be unsatisfactory, since conventional sequential algorithms may be serial in nature and may not be implemented efficiently on parallel machines. In many cases, new algorithms have to be introduced to increase parallel performance. In order to achieve optimal performance, in addition to partitioning and mapping, a careful performance study should be conducted for a given application to find a good algorithm-machine combination. This process, however, is usually painful and elusive. The goal of this project is to design and develop efficient parallel algorithms for highly accurate Computational Fluid Dynamics (CFD) simulations and other engineering applications. The work plan is 1) developing highly accurate parallel numerical algorithms, 2) conduct preliminary testing to verify the effectiveness and potential of these algorithms, 3) incorporate newly developed algorithms into actual simulation packages. The work plan has well achieved. Two highly accurate, efficient Poisson solvers have been developed and tested based on two different approaches: (1) Adopting a mathematical geometry which has a better capacity to describe the fluid, (2) Using compact scheme to gain high order accuracy in numerical discretization. The previously developed Parallel Diagonal Dominant (PDD) algorithm and Reduced Parallel Diagonal Dominant (RPDD) algorithm have been carefully studied on different parallel platforms for different applications, and a NASA simulation code developed by Man M. Rai and his colleagues has been parallelized and implemented based on data dependency analysis. These achievements are addressed in detail in the paper.

Description: This Wu Chang-Hua lecture is concerned with the development of analysis and computational capability for turbomachinery flows which is based on detailed flow field physics. A brief review of the work of Professor Wu is presented as well as a summary of the current NASA aeropropulsion programs. Two major areas of research are described in order to determine our predictive capabilities using modern day computational tools evolved from the work of Professor Wu. In one of these areas, namely transonic rotor flow, it is demonstrated that a high level of accuracy is obtainable provided sufficient geometric detail is simulated. In the second case, namely turbine heat transfer, our capability is lacking for rotating blade rows and experimental correlations will provide needed information in the near term. It is believed that continuing progress will allow us to realize the full computational potential and its impact on design time and cost.

Description: The turbulence module recently developed for the NPARC code has been extended to include wall functions. The Van Driest transformation is used so that the wall functions can be applied to both incompressible and compressible flows. The module is equipped with three two-equation K-epsilon turbulence models: Chien, Shih-Lumley and CMOTR models. Details of the wall functions as well as their numerical implementation are reported. It is shown that the inappropriate artificial viscosity in the near-wall region has a big influence on the solution of the wall function approach. A simple way to eliminate this influence is proposed, which gives satisfactory results during the code validation. The module can be easily linked to the NPARC code for practical applications.