ALBERT

Description: The SPARK3D and SPARK3D-PNS computer programs were developed to model 3-D supersonic, chemically reacting flow-fields. The SPARK3D code is a full Navier-Stokes solver, and is suitable for use in scramjet combustors and other regions where recirculation may be present. The SPARK3D-PNS is a parabolized Navier-Stokes solver and provides an efficient means of calculating steady-state combustor far-fields and nozzles. Each code has a generalized chemistry package, making modeling of any chemically reacting flow possible. Research activities by the Langley group range from addressing fundamental theoretical issues to simulating problems of practical importance. Algorithmic development includes work on higher order and upwind spatial difference schemes. Direct numerical simulations employ these algorithms to address the fundamental issues of flow stability and transition, and the chemical reaction of supersonic mixing layers and jets. It is believed that this work will lend greater insight into phenomenological model development for simulating supersonic chemically reacting flows in practical combustors. Currently, the SPARK3D and SPARK3D-PNS codes are used to study problems of engineering interest, including various injector designs and 3-D combustor-nozzle configurations. Examples, which demonstrate the capabilities of each code are presented.

Description: A comparative study was made with four different codes for solving the compressible Navier-Stokes equations using three different test problems. The first of these cases was hypersonic flow through the P8 inlet, which represents inlet configurations typical of a hypersonic airbreathing vehicle. The free-stream Mach number in this case was 7.4. This 2-D inlet was designed to provide an internal compression ratio of 8. Initial calculations were made using two state-of-the-art finite-volume upwind codes, CFL3D and USA-PG2, as well as NASCRIN, a code which uses the unsplit finite-difference technique of MacCormack. All of these codes used the same algebraic eddy-viscosity turbulence model. In the experiment, the cowl lip was slightly blunted; however, for the computations, a sharp cowl leading edge was used to simplify the construction of the grid. The second test problem was the supersonic (Mach 3.0) flow in a three-dimensional corner formed by the intersection of two wedges with equal wedge angles of 9.48 degrees. The flow in such a corner is representative of the flow in the corners of a scramjet inlet. Calculations were made for both laminar and turbulent flow and compared with experimental data. The three-dimensional versions of the three codes used for the inlet study (CFL3D, USA-PG3, and SCRAMIN, respectively) were used for this case. For the laminar corner flow, a fourth code, LAURA, which also uses recently-developed upwind technology, was also utilized. The final test case is the two-dimensional hypersonic flow over a compression ramp. The flow is laminar with a free-stream Mach number of 14.1. In the experiment, the ramp angle was varied to change the strength of the ramp shock and the extent of the viscous-inviscid interaction. Calculations were made for the 24-degree ramp configuration which produces a large separated-flow region that extends upstream of the corner.

Description: It is now known that Batchelor's trailing-line vortex is extremely unstable to small amplitude disturbances for swirl numbers in the neighborhood of 0.83. The results of numerical calculations are presented that show the response of the vortex in this range of swirl numbers to finite amplitude, temporal, helical disturbances. Phenomena observed include: (1) ejection of axial vorticity and momentum from the core resulting in the creation of secondary, separate vortices; (2) a great intensification of core axial vorticity and a weakening of core momentum; and (3) the production of azimuthal vorticity in the form of a tightly wrapped spiral wave. The second phenomenon eventually stablizes the vortex, which then smooths and gradually returns to an axisymmetric state. The calculations are mixed spectral-finite-difference, fourth-order accurate, and have been carried out at Reynolds numbers of 1000 to 2000. Some linearized results are also discussed in an attempt to explain the process of vortex intensification.

Description: The NASA-Lewis aircraft icing analysis program is composed of three major sub-programs. These sub-programs are ice accretion simulation, performance degradation evaluation, and ice protection system evaluation. These topics cover all areas of concern related to the simulation of aircraft icing and its consequences. The motivation for these activities is twofold, reduction of time and effort required in experimental programs and the ability to provide reliable information for aircraft certification in icing, over the complete range of environmental conditions. In addition to the analytical activities associated with development of these codes, several experimental programs are underway to provide verification information for existing codes. These experimental programs are also used to investigate the physical processes associated with ice accretion and removal for improvement of present analytical models. The NASA-Lewis icing analysis program is thus striving to provide a full range of analytical tools necessary for evaluation of the consequences of icing and of ice protection systems.

Description: The state-of-the-art in rotor blade drag prediction involves the use of two-dimensional airfoil tables to calculate the drag force on the blade. One of the most serious problems with the current methods is that they cannot be used for airfoils that have yet to be tested. Most of the drag prediction methods also do not take the Reynolds number or the rotational effects of the blade into account, raising doubts about the accuracy of the results. These problems are addressed with the development of an analytical method which includes the shape of airfoil, the effects of Reynolds number, and the rotational motion of the blade.

Description: Interactions between the rotors and wing of a rotary wing aircraft in hover have a significant detrimental effect on its payload performance. The reduction of payload results from the wake of lifting rotors impinging on the wing, which is at 90 deg angle of attack in hover. This vertical drag, often referred as download, can be as large as 15 percent of the total rotor thrust in hover. The rotor wake is a three-dimensional, unsteady flow with concentrated tip vortices. With the rotor tip vortices impinging on the upper surface of the wing, the flow over the wing is not only three-dimensional and unsteady, but also separated from the leading and trailing edges. A simplified two-dimensional model was developed to demonstrate the stability of the methodology. The flow model combines a panel method to represent the rotor and the wing, and a vortex method to track the wing wake. A parametric study of the download on a 20 percent thick elliptical airfoil below a rotor disk of uniform inflow was performed. Comparisons with experimental data are made where the data are available. This approach is now being extended to three-dimensional flows. Preliminary results on a wing at 90 deg angle of attack in free stream is presented.

Description: TranAir is a computer code which solves the full-potential equation for transonic flow about very general and complex configurations. Piecewise flat surface panels are used to describe the surface geometry. This paneled definition is then embedded in an unstructured cartesian flow field grid. Finite elements are used in the discretization of the flow field grid in a manner which is fully conservative and second-order accurate. Since geometries may be defined with relative ease, and since the user is not involved in the generation of the flow field grid, computational results may be generated rather quickly for a wide range of geometries. For transonic cases in the cruise angle-of-attack range, TranAir has generated results which are in generally good agreement with both Euler results and wind tunnel data. A typical transonic case runs in 1 to 2 CPU hours on a Cray X-MP. For subcritical cases, the code runs in 15 to 30 CPU minutes, even for geometries in which several thousand surface panels are used in the definition. This ability to rapidly and accurately provide both subsonic and transonic predictions about very complex aircraft configurations gives TranAir the potential of being a very powerful and widely used design tool.

Description: The objective is developing CFD capabilities to obtain solutions for viscous flows about generic configurations of internally and externally carried stores. The emphasis is placed on the supersonic flow regime with extensions being made to the transonic regime. The project is broken into four steps: (1) Cavity flows for internal carriage configurations; (2) High angle of attack flows, which may be experienced during the separation of the stores: (3) Flows about a body near a flat plate for external carriage configurations; and (4) Flows about a body inside or in the proximity of a cavity. Three-dimensional unsteady cavity flow solutions are obtained by an explicit, MacCormack algorithm, EMCAV3, for open, close, and transitional cavities. High angle of attack flows past cylinders are solved by an implicit, upwind algorithm. All the results compare favorably with the experimental data. For flows about multiple body configurations, the Chimera embedding scheme is modified for finite-volume and multigrid algorithms, MaGGiE. Then a finite volume, implicit, upwind, multigrid Navier-Stokes solver which uses on overlapped/embedded and zonal grids, VUMXZ3, is developed from the CFL3D code. Supersonic flows past a cylinder near a flat plate are computed using this code. The results are compared with the experimental data. Currently the VUMXZ3 code is being modified to accomplish step 4 of this project. Wind tunnel experiments are also being conducted for validation purposes.

Description: Three-dimensional viscous flow computations are presented for the F/A-18 forebody-LEX (Leading Edge EXtensions) geometry. Solutions are obtained from an algorithm for the compressible Navier-Stokes equations which incorporates an upwind-biased, flux-difference-splitting approach along with longitudinally-patched grids. Results are presented for both laminar and fully turbulent flow assumptions and include correlations with wind tunnel as well as flight-test results. A good quantitative agreement for the forebody surface pressure distribution is achieved between the turbulent computations and wind tunnel measurements at Mach number 0.6. The computed turbulent surface flow patterns on the forebody qualitatively agree well with in-flight surface flow patterns obtained on an F/A-18 aircraft at Mach number 0.34.

Description: Fine-grid Navier-Stokes solutions were obtained for flow over the fuselage forebody and wing leading edge extension of the F/A-18 High Alpha Research Vehicle at large incidence. The resulting flows are complex, and exhibit cross flow separation from the sides of the forebody and from the leading edge extension. A well-defined vortex pattern is observed in the leeward-side flow. Results obtained for laminar flow show good agreement with flow visualizations obtained in ground-based experiments. Further, turbulent flows computed at high Reynolds-number flight-test conditions show good agreement with surface and off-surface visualizations obtained in flight.

Description: The objective of the research is to develop and validate accurate, user-oriented viscous CFD codes (with inviscid options) for three-dimensional, unsteady aerodynamic flows about arbitrary rotorcraft configurations. Unsteady, three-dimensional Euler and Navier-Stokes codes are developed, adapted, and extended to rotor-body combinations. Flow solvers are coupled with zonal grid topologies, including rotating and nonrotating blocks. Special grid clustering and wave-fitting techniques were developed to capture low-level radiating acoustic waves. Significant progress was made in computing the propagation of acoustic waves due to the interaction of a concentrated vortex and a helicopter airfoil. The need for higher-order schemes was firmly established in relatively inexpensive two-dimensional calculations. In three dimensions, the number of grid points required to capture the low-level acoustic waves becomes very large, so that large supercomputer memory becomes essential. Good agreement was obtained between the numerical results obtained with a thin-layer Navier-Stokes code and experimental data from a model rotor. In addition, several nonrotating configurations that are sometimes proposed to simulate rotor blade tips in conventional wind tunnels were examined, and the complex flow around the radical tip shape of the world's fastest helicopter is under investigation. These studies demonstrate the flexibility and power of CFD to gain physical insight, study novel ideas, and examine various possibilities that might be difficult or impossible to set up in physical experiments. As a prelude to studies of rotor-body aerodynamic interactions, a preliminary grid topology and moving-interface strategy were developed. A new Euler/Navier-Stokes code using these techniques computes the vortical wake directly, rather than modeling it, as in most previous rotorcraft studies. Several hover cases were run for conventional and advanced-geometry blades. Numerical schemes using multi-zones and/or adaptive grids appear to be necessary to simulate the complex vortical flows in rotor wakes.

Description: The thin layer, Reynolds-averaged, Navier-Stokes equations are used to simulate the transonic viscous flow about the complete F-16A fighter aircraft. These computations demonstrate how computational fluid dynamics can be used to simulate turbulent viscous flow about realistic aircraft geometries. A zonal grid approach is used to provide adequate viscous grid clustering on all aircraft surfaces. Zonal grids extend inside the F-16A inlet and up to the compressor face while power on conditions are modeled by employing a zonal grid extending from the exhaust nozzle to the far field. Computations are compared with existing experimental data and are in fair agreement. Computations for the F-16A in side slip are also presented.

Description: Grid generation and Euler flow about fighter aircraft are described. A fighter aircraft geometry is specified by an area ruled fuselage with an internal duct, cranked delta wing or strake/wing combinations, canard and/or horizontal tail surfaces, and vertical tail surfaces. The initial step before grid generation and flow computation is the determination of a suitable grid topology. The external grid topology that has been applied is called a dual-block topology which is a patched C (exp 1) continuous multiple-block system where inner blocks cover the highly-swept part of a cranked wing or strake, rearward inner-part of the wing, and tail components. Outer-blocks cover the remainder of the fuselage, outer-part of the wing, canards and extend to the far field boundaries. The grid generation is based on transfinite interpolation with Lagrangian blending functions. This procedure has been applied to the Langley experimental fighter configuration and a modified F-18 configuration. Supersonic flow between Mach 1.3 and 2.5 and angles of attack between 0 degrees and 10 degrees have been computed with associated Euler solvers based on the finite-volume approach. When coupling geometric details such as boundary layer diverter regions, duct regions with inlets and outlets, or slots with the general external grid, imposing C (exp 1) continuity can be extremely tedious. The approach taken here is to patch blocks together at common interfaces where there is no grid continuity, but enforce conservation in the finite-volume solution. The key to this technique is how to obtain the information required for a conservative interface. The Ramshaw technique which automates the computation of proportional areas of two overlapping grids on a planar surface and is suitable for coding was used. Researchers generated internal duct grids for the Langley experimental fighter configuration independent of the external grid topology, with a conservative interface at the inlet and outlet.

Description: Information on time dependent incompressible Navier-Stokes equations is given in viewgraph form. Information is given on streamfunction equations for unsteady incompressible flow, the streamfunction algorithm for unsteady incompressible flow, and a multigrid solver for the laminar implicit equations.

Description: The purpose is to describe the development of unstructured grid methods which have several advantages when compared to methods which make use of structured grids. Unstructured grids, for example, easily allow the treatment of complex geometries, allow for general mesh movement for realistic motions and structural deformations of complete aircraft configurations which is important for aeroelastic analysis, and enable adaptive mesh refinement to more accurately resolve the physics of the flow. Steady Euler calculations for a supersonic fighter configuration to demonstrate the complex geometry capability; unsteady Euler calculations for the supersonic fighter undergoing harmonic oscillations in a complete-vehicle bending mode to demonstrate the general mesh movement capability; and vortex-dominated conical-flow calculations for highly-swept delta wings to demonstrate the adaptive mesh refinement capability are discussed. The basic solution algorithm is a multi-stage Runge-Kutta time-stepping scheme with a finite-volume spatial discretization based on an unstructured grid of triangles in 2D or tetrahedra in 3D. The moving mesh capability is a general procedure which models each edge of each triangle (2D) or tetrahedra (3D) with a spring. The resulting static equilibrium equations which result from a summation of forces are then used to move the mesh to allow it to continuously conform to the instantaneous position or shape of the aircraft. The adaptive mesh refinement procedure enriches the unstructured mesh locally to more accurately resolve the vortical flow features. These capabilities are described in detail along with representative results which demonstrate several advantages of unstructured grid methods. The applicability of the unstructured grid methodology to steady and unsteady aerodynamic problems and directions for future work are discussed.

Description: SAGE is a user-friendly, highly efficient, two-dimensional self-adaptive grid code based on Nakahashi and Deiwert's variational principles method. Grid points are redistributed into regions of high flowfield gradients while maintaining smoothness and orthogonality of the grid. Efficiency is obtained by splitting the adaption into 2 directions and applying one-sided torsion control, thus producing a 1-D elliptic system that can be solved as a set of tridiagonal equations.

Description: Strong interactions of structures and fluids are common in many engineering environments. Such interactions can give rise to physically important phenomena such as those occurring for aircraft due to aeroelasticity. Aeroelasticity can significantly influence the safe performance of aircraft. At present exact methods are available for making aeroelastic computations when flows are in either the linear subsonic or supersonic range. However, for complex flows containing shock waves, vortices and flow separations, computational methods are still under development. Several phenomena that can be dangerous and limit the performance of an aircraft occur due to the interaction of these complex flows with flexible aircraft components such as wings. For example, aircraft with highly swept wings experience vortex induced aeroelastic oscillations. Correct understanding of these complex aeroelastic phenomena requires direct coupling of fluids and structural equations. Here, a summary is presented of the development of such coupled methods and applications to aeroelasticity since about 1978 to present. The successful use of the transonic small perturbation theory (TSP) coupled with structures is discussed. This served as a major stepping stone for the current stage of aeroelasticity using computational fluid dynamics. The need for the use of more exact Euler/Navier-Stokes (ENS) equations for aeroelastic problems is explained. The current development of unsteady aerodynamic and aeroelastic procedures based on the ENS equations are discussed. Aeroelastic results computed using both TSP and ENS equations are discussed.

Description: A thin-layer Navier-Stokes has been developed for solving high Reynolds number, turbulent flows past aircraft components under transonic flow conditions. The computer code has been validated through data comparisons for flow past isolated wings, wing-body configurations, prolate spheroids and wings mounted inside wind-tunnels. The basic code employs an explicit Runge-Kutta time-stepping scheme to obtain steady state solution to the unsteady governing equations. Significant gain in the efficiency of the code has been obtained by implementing a multigrid acceleration technique to achieve steady-state solutions. The improved efficiency of the code has made it feasible to conduct grid-refinement and turbulence model studies in a reasonable amount of computer time. The non-equilibrium turbulence model of Johnson and King has been extended to three-dimensional flows and excellent agreement with pressure data has been obtained for transonic separated flow over a transport type of wing.

Description: The INS3D family of computational fluid dynamics computer codes is presented. These codes are used to as tools in developing and assessing algorithms for solving the incompressible Navier-Stokes equations for steady-state and unsteady flow problems. This work involves applying the codes to real-world problems involving complex three-dimensional geometries. The algorithms utilized include the method of pseudocompressibility including both central and upwind differencing, several types of artificial dissipation schemes, approximate factorization, and an implicit line-relaxation scheme. These codes have been validated using a wide range of problems including flow over a backward-facing step, driven cavity flow, flow through various types of ducts, and steady and unsteady flow over a circular cylinder. Many diverse flow applications have been solved using these codes including parts of the Space Shuttle Main Engine, problems in naval hydrodynamics, low-speed aerodynamics, and biomedical fluid flows. The presentation details several of these, including the flow through a Space Shuttle Main Engine inducer, vortex shedding behind a circular cylinder, and flow through an artificial heart.

Description: Significant improvements in predictive accuracies for off-design conditions are achievable through better turbulence modeling; and, without necessarily adding any significant complication to the numerics. One well established fact about turbulence is it is slow to respond to changes in the mean strain field. With the 'equilibrium' algebraic turbulence models no attempt is made to model this characteristic and as a consequence these turbulence models exaggerate the turbulent boundary layer's ability to produce turbulent Reynolds shear stresses in regions of adverse pressure gradient. As a consequence, too little momentum loss within the boundary layer is predicted in the region of the shock wave and along the aft part of the airfoil where the surface pressure undergoes further increases. Recently, a 'nonequilibrium' algebraic turbulence model was formulated which attempts to capture this important characteristic of turbulence. This 'nonequilibrium' algebraic model employs an ordinary differential equation to model the slow response of the turbulence to changes in local flow conditions. In its original form, there was some question as to whether this 'nonequilibrium' model performed as well as the 'equilibrium' models for weak interaction cases. However, this turbulence model has since been further improved wherein it now appears that this turbulence model performs at least as well as the 'equilibrium' models for weak interaction cases and for strong interaction cases represents a very significant improvement. The performance of this turbulence model relative to popular 'equilibrium' models is illustrated for three airfoil test cases of the 1987 AIAA Viscous Transonic Airfoil Workshop, Reno, Nevada. A form of this 'nonequilibrium' turbulence model is currently being applied to wing flows for which similar improvements in predictive accuracy are being realized.

Description: Presented here are several direct simulations of one 2-D second mode perturbation wave, superimposed upon a prescribed mean flow. Periodicity is assumed in the streamwise direction (Fourier) and the variables are expanded in Chebyshev series in the direction normal to the flat plate. The code is fully explicit and is time advanced with a 3rd order Runge-Kutta scheme. The second mode wave (R delta prime = 8000), interacts with itself to generate higher streamwise harmonics. Physical parameters are chosen to maximize the linear growth rate at the prescribed Reynolds number. Initial results indicate that the nonlinear processes begin in the critical layer region and are the result of the cubic interactions in the momentum equations, rather than due to the higher streamwise harmonics. Analysis of the various terms in the momentum equations combined with numerical experiments in which various modes are artificially suppressed, lead to the conclusion that asymptotic methods will produce the saturated state in one or two order of magnitude less computer time than that required by the direct numerical simulations.

Description: Direct simulations consist in solving the full Navier-Stokes equations, without any turbulence model, and describing all the detailed features of the flow. Usually the flows are three-dimensional and time-dependent and contain both coarse and fine structures, which makes the numerical task very challenging in terms of both the algorithm and the computational effort. Most of the work until now has involved spectral methods, which are highly accurate but not very flexible in terms of geometry or complex equations. For that reason, future work will also rely on high-order finite-difference or other methods. Direct simulations complement experimental work, and both contribute to the theory and the empirical knowledge of turbulence. Once such a simulation has been shown to be accurate, the flow field is completely known in three dimensions and time, including the pressure, the vorticity and any other quantity. On the other hand, most simulations to date solved the incompressible equations in rather simple geometries, and direct simulations will always be limited to moderate Reynolds numbers. Extensive simulations have been conducted in homogeneous turbulence, channel flows, boundary layers, and mixing layers. Much effort is devoted to addressing flows with compressibility and chemical reactions, and to new geometries such as a backward-facing step.

Description: Recent applications and development of CFD technology have focused on flow problems that are critically important to the operation and design of space flight vehicles. The main effort is spent on the Space Shuttle in order to provide an understanding of the cryogenic fluid in the duct connecting the External Tank and the Main Engines, the subsonic flow surrounding the Orbiter during crew egress maneuvers, the transonic aerodynamic forces on the Orbiter fuselage and wing, the high angle-of-attack abort flight, and the aerodynamic heating during entry. To provide in-depth analyses for such diverse problems within a timely schedule, matured panel codes and a state-of-the-art incompressible turbulent flow code were adapted. Collaboration with Ames Research Center has resulted in a Shuttle ascent aerodynamic code; and a viscous chemical nonequilibrium code is being developed for predicting Orbiter real-gas aerodynamics and finite-catalytic heating. The remaining activities are devoted to the prediction of the flow environment around the Aeroassist Flight Experiment vehicle at hypersonic speeds and high altitudes. A thermochemical nonequilibrium Navier-Stokes code has been developed on the basis of two- temperature and 11-species models for solving both the shock layer and near wake. After validating the code against wind-tunnel aerodynamic, pressure and heating data, the code is being used to supplement the ground test facilities in predicting a more realistic flight environment. CFD technology is being relied upon by other programs as well in the consideration of candidate configurations.

Description: Several direct simulations of 3-D homogeneous, compressible turbulence are presented with emphasis on the differences with incompressible turbulent simulations. A fully spectral collocation algorithm, periodic in all directions coupled with a 3rd order Runge-Kutta time discretization scheme is sufficient to produce well-resolved flows at Taylor Reynolds numbers below 40 on grids of 128x128x128. A Helmholtz decomposition of velocity is useful to differentiate between the purely compressible effects and those effects solely due to vorticity production. In the context of homogeneous flows, this decomposition in unique. Time-dependent energy and dissipation spectra of the compressible and solenoidal velocity components indicate the presence of localized small scale structures. These structures are strongly a function of the initial conditions. Researchers concentrate on a regime characterized by very small fluctuating Mach numbers Ma (on the order of 0.03) and density and temperature fluctuations much greater than sq Ma. This leads to a state in which more than 70 percent of the kinetic energy is contained in the so-called compressible component of the velocity. Furthermore, these conditions lead to the formation of curved weak shocks (or shocklets) which travel at approximately the sound speed across the physical domain. Various terms in the vorticity and divergence of velocity production equations are plotted versus time to gain some understanding of how small scales are actually formed. Possible links with Burger turbulence are examined. To visualize better the dynamics of the flow, new graphic visualization techniques have been developed. The 3-D structure of the shocks are visualized with the help of volume rendering algorithms developed in-house. A combination of stereographic projection and animation greatly increase the number of visual cues necessary to properly interpret the complex flow.

Description: Computational Fluid Dynamics (CFD) activities at Marshall Space Flight Center (MSFC) have been focused on hardware specific and research applications with strong emphasis upon benchmark validation. The purpose here is to provide insight into the MSFC CFD related goals, objectives, current hardware related CFD activities, propulsion CFD research efforts and validation program, future near-term CFD hardware related programs, and CFD expectations. The current hardware programs where CFD has been successfully applied are the Space Shuttle Main Engines (SSME), Alternate Turbopump Development (ATD), and Aeroassist Flight Experiment (AFE). For the future near-term CFD hardware related activities, plans are being developed that address the implementation of CFD into the early design stages of the Space Transportation Main Engine (STME), Space Transportation Booster Engine (STBE), and the Environmental Control and Life Support System (ECLSS) for the Space Station. Finally, CFD expectations in the design environment will be delineated.

Description: Lewis is a multidisciplinary Center with strong research and development programs in aeronautical and space propulsion, power, space communications, space experiments and materials. Computational fluid dynamics (CFD) is playing an important and growing role in most of these areas. Described here is how CFD is integrated into these programs and highlights elements of the CFD activities. Examples are presented of codes developed to predict flow fields in advanced propulsion systems and several of the code validation experiments are described. The CFD effort at Lewis ranges from basic research on new and improved algorithms through code development to the application of these codes to specific engineering problems. Because of the substantial improvement in CFD's predictive capability, its use at Lewis is on a steep growth path, spreading rapidly into new areas which had not traditionally taken advantage of the techniques of numerical simulation. Multidisciplinary codes and the future direction of CFD at Lewis are discussed.

Description: Information on computational fluid dynamics (CFD) research and applications carried out at the NASA Langley Research Center is given in viewgraph form. The Langley CFD strategy, the five-year plan in CFD and flow physics, 3-block grid topology, the effect of a patching algorithm, F-18 surface flow, entropy and vorticity effects that improve accuracy of unsteady transonic small disturbance theory, and the effects of reduced frequency on first harmonic components of unsteady pressures due to airfoil pitching are among the topics covered.

Description: The Computational Fluid Dynamics (CFD) Program at NASA Ames Research Center is reviewed and discussed. The technical elements of the CFD Program are listed and briefly discussed. These elements include algorithm research, research and pilot code development, scientific visualization, advanced surface representation, volume grid generation, and numerical optimization. Next, the discipline of CFD is briefly discussed and related to other areas of research at NASA Ames including experimental fluid dynamics, computer science research, computational chemistry, and numerical aerodynamic simulation. These areas combine with CFD to form a larger area of research, which might collectively be called computational technology. The ultimate goal of computational technology research at NASA Ames is to increase the physical understanding of the world in which we live, solve problems of national importance, and increase the technical capabilities of the aerospace community. Next, the major programs at NASA Ames that either use CFD technology or perform research in CFD are listed and discussed. Briefly, this list includes turbulent/transition physics and modeling, high-speed real gas flows, interdisciplinary research, turbomachinery demonstration computations, complete aircraft aerodynamics, rotorcraft applications, powered lift flows, high alpha flows, multiple body aerodynamics, and incompressible flow applications. Some of the individual problems actively being worked in each of these areas is listed to help define the breadth or extent of CFD involvement in each of these major programs. State-of-the-art examples of various CFD applications are presented to highlight most of these areas. The main emphasis of this portion of the presentation is on examples which will not otherwise be treated at this conference by the individual presentations. Finally, a list of principal current limitations and expected future directions is given.

Description: An experimental simulation of an in-flight wingtip vortical flow visualization technique uses infrared imaging to observe strong and concentrated vortices. This experiment is phase 1 of a two-phase infrared evaluation program. The system includes a vortex generator (model 320 Vortec Vortex Tube) which generates the required vortex. The mouth of the unit is mounted close to the free end of a half-inch diameter, sixteen and a half foot long stainless steel tubing (sized after tubing currently installed in the wings of an experimental Beechcraft Sundowner 180 aircraft). Dichloro difluoromethane (Freon-12) is entrained into the generated vortex. A breakdown of the vortices is indicated by the rapid diffusion and the resulting pattern is tracked using the infrared imager and video systems. Flow rates (volume and mass) are estimated at the laboratory and proposed flight conditions. The nominal flight altitude is expected to be 2500 feet.

Description: A common airfoil model with the CAST 10-2/DOA-2 profile and 228 mm (9 inches) chord length was tested. The tests performed in NAE covered the Mach numbers from 0.3 to 0.8 and Reynolds numbers from 10 to 30 million. The model was tested with transition free and with transition fixed at 5 percent chord for both the upper and the lower surfaces. The data obtained were analyzed for the effects of Reynolds number, transition fixing and Mach number. The role of the boundary layer on the displacement effect, the interaction with the shock wave and the trailing edge separation are examined. The results are summarized as follows: (1) the airfoil performance depends strongly on Reynolds number and transition fixing; (2) with transition fixed, the aerodynamic quantities such as lift, pitching moment and drag show a monotonic variation with Reynolds number; (3) with transition free, the aerodynamic quantities vary less regularly with Reynolds number and a slight parametric dependency is shown. The weak dependency is due to the compensatory effect of the forward shift of the transition position and the thinning of the turbulent boundary layer as Reynolds number increases; (4) the shock Mach number and the shock position are weakly dependent on Reynolds number; and (5) the long extent of the laminar boundary layer at transonic speeds reduces the drag appreciably at low Reynolds numbers. The drag bucket around the design Mach number can be observed below Reynolds number 15 million.

Description: The aim of this cooperative NASA/DFVLR/ONERA project was to examine the performance of the CAST-10 airfoil in the T2 cryogenic wind tunnel. Tests included general characteristics of the CAST-10 airfoil and fundamental studies on Reynold number effects. Good T2 cryogenic operation was observed. Improvements should be done for moisture elimination and for side wall boundary layer effects.

Description: The development of missiles from early history up to about 1970 is discussed. Early unpowered missiles beyond the rock include the spear, the bow and arrow, the gun and bullet, and the cannon and projectile. Combining gunpowder with projectiles resulted in the first powered missiles. In the early 1900's, the development of guided missiles was begun. Significant advances in missile technology were made by German scientists during World War II. The dispersion of these advances to other countries following the war resulted in accelerating the development of guided missiles. In the late 1940's and early 1950's there was a proliferation in the development of missile systems in many countries. These developments were based primarily on experimental work and on relatively crude analytical techniques. Discussed here are some of the missile systems that were developed up to about 1970; some of the problems encountered; the development of an experimental data base for use with missiles; and early efforts to develop analytical methods applicable to missiles.

Description: Results of the simulation of the viscous flow past the CAST-10 airfoil were shown for different flow conditions. Since the experiments provide only surface pressures and force coefficients, the comparison to the numerical results relies on these. Good agreement of these results is found for the lower Mach number cases except for the shock position. As numerical experiments indicate, this seems to be due to the turbulent shock boundary layer interaction which is not correctly modelled by the algebraic turbulence model employed. For the lower Mach number case the influence of the transition location was investigated, too. Changing the transition location at the lower surface has much more influence on the pressure distribution than changing it on the upper side. For the higher Mach number case the double shock structure found in the experiment for the lower Reynolds numbers was not reproduced by the numerical solutions. The reason for this is unknown though it may be due to the turbulence modelling. For the higher Reynolds number a better resolution of the boundary layer is needed in the computation in order to recover the experimental pressure plateau; but then the shock position is still found downstream of the experimental one.

Description: The objective of the paper is two-fold. First, an upwind/central differencing method for solving the steady Navier-Stokes equations is described. The symmetric line relation method is used to solve the resulting algebraic system to achieve high computational efficiency. The grid spacings used in the calculations are determined from the triple-deck theory, in terms of Mach and Reynolds numbers and other flow parameters. Thus the accuracy of the numerical solutions is improved by comparing them with experimental, analytical, and other computational results. Secondly, the shock wave/boundary layer interactions are studied numerically, with special attention given to the flow separation. The concept of free interaction is confirmed. Although the separated region varies with Mach and Reynolds numbers, it is found that the transverse velocity component behind the incident shock, which has not been identified heretofore, is also an important parameter. A small change of this quantity is sufficient to eliminate the flow separation entirely.

Description: An experimental and computational study of the low-frequency oscillation observed in the flow over an airfoil at the onset of static stall is presented. Wind-tunnel results obtained with two-dimensional airfoil models show that this phenomena takes place only with a transitional state of the separating boundary layer. It is noted that the flowfield does not involve a Karman vortex street. The experimental results agree well with the results of a two-dimensional Navier-Stokes code. The present study demonstrates that the low-frequency oscillations produce intense flow fluctuations which impart much larger unsteady forces to the airfoil than experienced by bluff-body shedding and which may represent the primary aerodynamics of stall flutter of blades and wings.

Description: Thin-element polymeric film riblet surfacing materials previously shown to behave similarly to v-groove riblets in the reduction of aerodynamic drag have been tested to ascertain the effects of groove contamination. Three different levels of contamination were simulated for each of two kinds of contamination. Such contaminants as oil drops and condensates have no significant effect on drag reduction. Atmospheric particulates may, however, degrade riblet performance.

Description: An experimental study has been conducted on the effects of initial swirl and high turbulence in the exhaust of a circular jet injected from a flat plate at a 90-deg angle into a crossflow. The different jet types studied were low-exit turbulence (3 percent), high-exit turbulence (over 10 percent), and 40 and 58 percent swirl. Surface pressure distributions and mean velocity vector plots were obtained for all of these cases. For the surface pressure distribution tests, the jet to crossflow velocity ratios R were 2.2, 4, and 8 for most of the jet types. For the mean velocity vector plots, R = 4 was chosen. Turbulence information in the jet plume was also obtained for the low-exit turbulence case at R = 4. The results showed that the higher-exit turbulence reduced the penetration height, and it also reduced the surface area influenced by negative pressures. The swirl-caused asymmetric pressure distributions and the swirl effects were more pronounced for lower-velocity ratios.

Description: A three-dimensional, Reynolds-averaged Navier-Stokes code has been used to numerically analyze flow through a two-strut, supersonic combustion ramjet (scramjet) inlet configuration. It solves the governing equations in full conservation form using either a fully explicit or explicit-implicit method. An algebraic, two-layer eddy-viscosity model is used for turbulent flow calculations. The analysis allows inclusion of end effects that are caused by the aft placement of the cowl on the underside of the inlet. A special grid has been developed to accommodate the struts embedded in the inlet flowfield. Detailed numerical results are presented here for the two-strut configuration, and a comparison is made with the available experimental results.

Description: Several ways for obtaining aerodynamic parameters of an aircraft from flight data are presented. A brief description of data analysis from steady measurements is given, and a concept of system identification is introduced. This is followed by a mathematical model of an aircraft with aerodynamic forces and moments approximated either by polynomials or splines and a detailed treatment of two often-used techniques for parameter estimation. Because of renewed interest in frequency-domain analysis, one section of the paper is devoted to this problem. The methods are demonstrated in several examples using real flight data.

Description: An effort is made to contribute to the development of CFD by relating the successful use of vortex dynamics in the computation of the pressure drop past a planar array of chevron-shaped obstructions. An ensemble of results was used to compute the loss coefficient k, stimulating an experimental program for the assessment of the measured loss coefficient for the same geometry. The most provocative result of this study has been the representation of kinetic energy production in terms of vorticity source terms.

Description: Computational and experimental performance benefits are presented for a high-aspect-ratio unswept wing configuration with sheared tips. The sheared tip is a highly swept and highly tapered surface located in the same plane as the inboard wing panel to which it is attached. The compuational results were obtained with an inviscid surface panel method that models the nonlinear influence of the trailing wake. Both wind-tunnel and calculated results were obtained for a 12-ft span wing model with various wing-tip configurations. The computational and experimental data are in fair agreement and demonstrate that sheared wing tips can reduce induced drag at cruise and climb lift coefficients. The drag reduction is the result of wake deformation effects and changes in spanwise load distribution. Wind-tunnel measured longitudinal and lateral directional stability characteristics are also presented for the various wing-tip layouts.

Description: Standard techniques used to model chemically-reacting flows require an artificial viscosity for stability in the presence of strong shocks. The resulting shock is smeared over at least three computational cells, so that the thickness of the shock is dictated by the structure of the overall mesh and not the shock physics. A gas passing through a strong shock is thrown into a nonequilibrium state and subsequently relaxes down over some finite distance to an equilibrium end state. The artificial smearing of the shock envelops this relaxation zone which causes the chemical kinetics of the flow to be altered. A method is presented which can investigate these issues by following the chemical kinetics and flow kinetics of a gas passing through a fully resolved shock wave at hypersonic Mach numbers. A nonequilibrium chemistry model for air is incorporated into a spectral multidomain Navier-Stokes solution method. Since no artificial viscosity is needed for stability of the multidomain technique, the precise effect of this artifice on the chemical kinetics and relevant flow features can be determined.

Description: A complex three-dimensional shock-wave/turbulent boundary layer interaction at Mach 4 has been investigated experimentally and computationally, using two turbulence models with substantial refinement. With the use of a fine grid model, secondary flow separation was successfully computed, and grid changes did not improve agreement with experiment for the extent of upstream influence. The use of a non-isotropic turbulence model gave a slight improvement in upstream influence, but the size of the interaction was still significantly less than in the experiment.

Description: In the theory of thin lifting surfaces the minimum drag consistent with a given total lift occurs when the downwash, averaged between forward and reversed motion of the wing, has the same value at all points of the wing planform. In Kogan's theory the conditions for minimum drag are determined on the forward sloping characteristic surface touching the trailing edge of the wing and it is shown that such a surface plays the role of the Trefftz plane familiar in subsonic wing theory. This paper shows how Kogan's theory may be applied to determine the drag of elliptic wings at supersonic speed. It appears that such wings have lower drag than the conventional 'delta' wing.

Description: A class of highly amplified supersonic disturbances are found for high-speed, bounded mixing layers at high values of streamwise wavenumber. Their amplification is an order of magnitude greater than the most amplified modes, which occur at 60-65 deg at low streamwise wavenumber. These disturbances are stabilized by increasing Mach number, viscosity, and sweep; however, the effect of sweep on the most amplified mode is not significant until the wave propagation angle reaches 30 deg. The maximum growth rate of the unstable disturbances decreases as the temperature of the higher Mach number stream is increased. The structure of these disturbances is such that the phase speed with respect to the mean flow is subsonic in a small region in the center of the shear layer, and supersonic on either side of this region.

Description: Numerical solutions of the compressible Navier-Stokes equations assuming a perfect gas and laminar flow are applied to two hypersonic experiments conducted at NASA Langley Research Center. The first is cooling by injection of a gas jet through the nose of an ogive-cone, and the second is the aerothermal environment in the gap formed by the wing and elevon section of a test model of the Space Shuttle. The simulations demonstrate that the model predicts accurate pressures and shock stand-off for the gas jet and accurate pressures for the wing-elevon cove. Experimental heating rates for the gas jet indicate turbulent mixing to be occurring, an effect not modeled in the numerical solution.

Description: Results obtained from iterative solutions of inviscid and boundary-layer equations are presented and compared with experimental values. The calculated results were obtained with an Euler code and a transonic potential code in order to furnish solutions for the inviscid flow; they were interacted with solutions of two-dimensional boundary-layer equations having a strip-theory approximation. Euler code results are found to be in better agreement with the experimental data than with the full potential code, especially in the presence of shock waves, (with the sole exception of the near-tip region).

Description: A series of experiments are described which investigate and quantify the effect of loading on the three-dimensional flow through a subsonic annular cascade of cambered airfoils. At two levels of loading, detailed data quantify the cascade inlet velocity, the intrapassage flow field, the airfoil surface pressure distributions, the exit flow field, and the total pressure loss distributions. Aerodynamic loading is shown to strengthen the radial pressure gradient, the passage vortex structure, the vortex-endwall boundary layer interactions, and the losses.

Description: The problem is to determine the transitional flow aerodynamics and aerothermodynamics, including the base flow characteristics, of the Aeroassist Flight Experiment (AFE). The justification for the computational fluid dynamic (CFD) Application stems from MSFC's system integration responsibility for the AFE. To insure that the AFE objectives are met, MSFC must understand the limitations and uncertainties of the design data. Perhaps the only method capable of handling the complex physics of the rarefied high energy AFE trajectory is Bird's Direct Simulation Monte Carlo (DSMC) technique. The 3-D code used in this analysis is applicable only to the AFE geometry. It uses the Variable Hard Sphere (VHS) collision model and five specie chemistry model available from Langley Research Center. The code is benchmarked against the AFE flight data and used as an Aeroassisted Space Transfer Vehicle (ASTV) design tool. The code is being used to understand the AFE flow field and verify or modify existing design data. Continued application to lower altitudes is testing the capability of the Numerical Aerodynamic Simulation Facility (NASF) to handle 3-D DSMC and its practicality as an ASTV/AFE design tool.

Description: The instability of an attachment-line boundary layer formed on a swept cylinder in a supersonic freestream is considered in the linear regime. The supersonic attachment-line boundary layer is shown to be susceptible to oblique Tollmien-Schlichting wave instability which may be controlled by wall cooling. The critical Reynolds number based upon momentum thickness is found to be about 230. The onset of transition in the attachment-line boundary layer is also studied using the e (sup N) method and results are compared with the experimental data obtained at M (sub infinity) = 3.5 in the absence of any trips or spanwise contamination.

Description: Results from the Viscous Transonic Airfoil Workshop held in January 1987, are compared with each other and with experimental data. Test cases used include attached and separated transonic flows for the NACA 0012 airfoil. A total of 23 sets of numerical results from 15 different author groups are included. The numerical methods used vary widely and include: 16 Navier-Stokes methods, 2 Euler/boundary-layer methods, and 5 potential/boundary-layer methods. The results indicate a high degree of sophistication among the numerical methods with generally good agreement between the various computed and experimental results for attached or moderately separated cases. The agreement for cases with larger separation is only fair and suggests additional work is required in this area.

Description: Theoretical predictions based on the marching technique are compared with experimental observations on an airfoil with a concave region. Theoretical predictions of the wavelength of the most amplified Goertler vortex are in excellent agreement with the experimental observation for the range of chord Reynolds numbers from 1.0 to 3.67 million. In the convex zone, solutions from the marching technique showed that the initial counter-rotating vortex pairs lift off the surface and dissipate while another layer of vortex pairs of opposite rotation develops near the surface. This confirms the experimentally observed double peaked streamwise velocity perturbations. Furthermore, the streamwise velocity perturbations which dominate spanwise variation in the surface shear stress distribution shift by half a wavelength in the convex region. The experimental flow visualization photographs clearly confirm this phenomena.

Description: Instability, turbulence, and relaminarization in the attachment-line region of swept and unswept cylindrical bodies are studied by numerical solution of the full Navier-Stokes equations. The flow is simulated over a strip containing the attachment-line and treated as homogeneous in the spanwise direction; the disturbances decay exponentially upstream. Transpiration through the wall may be prescribed. The new method, which admits completely general disturbance, agrees with published linear-stability results, which were limited to an apparently restrictive form of disturbance. Fully developed turbulent solutions with sweep are generated and compare well with the experiment. The turbulence is subcritical (except for blowing), resulting in large hysteresis loops. By lowering the sweep Reynolds number, or increasing the suction, the turbulent flow is made to relaminarize. The relaminarization Reynolds number is much less sensitive to suction than the linear-stability Reynolds number. Extensive attempts to detect the postulated nonlinear instability of the unswept flow failed, suggesting that this flow is linearly and nonlinearly stable.

Description: Computational experiments have been performed for a few configurations in order to investigate the effects of external flow disturbances on the extent of laminar flow and wake drag. Theoretical results have been compared with experimental data for the AEDC cone, for Mach numbers from subsonic to supersonic, and for both free flight and wind tunnel environments. The comparisons have been found to be very satisfactory, thus establishing the utility of the present method for the design and development of laminar flow configurations and for the assessment of wind tunnel data. In addition, results of calculations concerning the effects of unit Reynolds numbers on transition are presented. In addition to the AEDC cone, computations have been performed for an ogive body of revolution at zero angle of attack and supersonic Mach numbers. Results are presented for transition Reynolds number and wake drag for external disturbances corresponding to free air and the test section of the AEDC-VKF tunnel. These results have been found to compare quite well with wind tunnel data for cases when surface suction is applied as well as when suction is absent.

Description: It is shown that the aerodynamic force and moment can be expressed in terms of vorticity distribution (and entropy variation for compressible flow) on near wake plane, or in terms of boundary vorticity flux on the body surface. Thus the vortical sources of lift and drag are clearly identified, which is the real physical basis of optimal aerodynamic design. Moreover, these sources are highly compact, hence allowing one to concentrate on key local regions of the configuration, which have dominating effect to the lift and drag. A detail knowledge of the vortical low requires measuring or calculating the vorticity and dilatation field, which is however still a challenging task. Nevertheless, this type of formulation has some unique advantages; and how to set up a well-posed problem, in particular how to establish vorticity-dilatation boundary conditions, is addressed.

Description: The dominance of body-like shapes for hypersonic vehicles places emphasis on better understanding of pertinenet body flow physics. Computational methods are implemented to enhance the knowledge of drag components generated by the production of lift forces at hypersonic speeds. Three forebody shapes are examined to identify the effects related to body shape parameters. A Navier-Stokes code and classical Newtonian theory code provide predictions of the polar shapes upon which conclusions are drawn. In particular, it is noted that hypersonic body polar shapes are somewhat irregular; the result of incidence-induced form drag. Test data which supports this finding is identified. The means for reducing hypersonic body lift-induced drag are also identified, along with applied computational schemes that can reduce the cost of a configuration design program.

Description: A fundamentally new approach to the aircraft minimum induced drag problem is presented. The method, a 'viscous lifting line', is based on the minimum entropy production principle and does not require the planar wake assumption. An approximate, closed form solution is obtained for several wing configurations including a comparison of wing extension, winglets, and in-plane wing sweep, with and without a constraint on wing-root bending moment. Like the classical lifting-line theory, this theory predicts that induced drag is proportional to the square of the lift coefficient and inversely proportioinal to the wing aspect ratio. Unlike the classical theory, it predicts that induced drag is Reynolds number dependent and that the optimum spanwise circulation distribution is non-elliptic.

Description: Induced drag is associated with the shedding of vorticity along the span of a finite wing, especially its tip region; for most subsonic aircraft configurations, induced drag constitutes about 50 percent of total aircraft drag throughout the flight envelope. NASA and the U.S. aircraft industry have aggressively studied induced-drag reduction methods. The state-of-the-art CTOL commercial aircraft wing is as a result of these efforts virtually optimal, with a total induced drag lying within a percent of the theoretical minimum. Many of the devices currently under study for induced drag reduction are added to wingtips, yielding benefits through their effects on the wake vortex as well as through forces generated in the flowfield.

Description: A scramjet exhaust simulation technique for hypersonic wind tunnel testing has been developed. Mixtures of Argon and Freon correctly match the inviscid simulation parameters of Mach number, static-pressure ratio, and the ratio of specific heats at the combustor exit location; this simulation is accomplished at significantly reduced temperatures and without combustion. An investigation of nozzle parametrics in a Mach 6 freestream showed that the external nozzle ramp angle, the cowl trailing-edge angle, an external nozzle flow fence and the nozzle static-pressure ratio significantly affected the external nozzle thrust and pitching moment as measured by the integration of surface-pressure data. A comparison of Argon-Freon and air exhaust simulation showed that the external nozzle thrust and pitching moment were in error by roughly a factor of 2 using air due to the incorrect match of the ratio of specific heats. An assessment of two-dimensional Euler and Navier-Stokes codes for predicting external nozzle aerodynamic characteristics was made by comparing computed and experimental results.

Description: The nonlinear interactions of high amplitude instability waves in turbulent jets are described. In plane shear layers Riley and Metcalf (1980) and Monkewitz (1987) have shown that these interactions are dependent, among other parameters, on the phase-difference between the two instability waves. Therefore, here researchers consider the nonlinear development of both the amplitudes and the phase of the instability waves. The development of these waves are also coupled with the development of the mean flow and the background turbulence. In formulating this model it is assumed that each of the flow components can be characterized by conservation equations supplemented by closure models. Results for the interactions between the two instability waves under high-amplitude forcing at fundamental and subharmonic frequencies are presented here. Qualitative agreements are found between the present predictions and available experimental data.

Description: Three-dimensional boundary-layer experiments are currently being conducted on a 45 deg swept wing in the Arizona State University Unsteady Wind Tunnel. Crossflow-dominated transition is produced via a model with contoured end liners to simulate infinite swept-wing flow. Fixed-wavelength stationary and traveling crossflow vortex are observed. The stationary vortex wavelengths vary with Reynolds number as predicted by linear-stability theory, but with observed wavelengths which are about 25 percent smaller than theoretically predicted. The frequencies of the most-amplified moving waves are in agreement with linear-stability theory; traveling waves at higher frequencies than predicted are also observed. These higher-frequency waves may be harmonics of the primary crossflow waves generated by a parametric resonance phenomena. Boundary-layer profiles measured at several spanwise locations show streamwise disturbance profiles characteristic of the crossflow instability. Near the transition location, severe distortions of the laminar boundary-layer profiles are observed. Mean-profile measurements are compared to theoretical predictions.

Description: Advanced Short Takeoff/Vertical Landing (STOVL) aircraft which can operate from remote locations, damaged runways, and small air capable ships are being pursued for deployment around the turn of the century. To achieve this goal, a cooperative program has been defined for testing in the NASA Lewis 9- by 15-foot Low Speed Wind Tunnel (LSWT) to establish a database for hot gas ingestion, one of the technologies critical to STOVL. This paper presents results showing the engine inlet distortions (both temperature and pressure) in a 9.2 percent scale Vectored Thrust STOVL model in ground effects. Results are shown for the forward nozzle splay angles of 0, -6, and 18 deg. The model support system had 4 deg of freedom, heated high pressure air for nozzle flow, and a suction system exhaust for inlet flow. The headwind (freestream) velocity was varied from 8 to 23 kn.

Description: This work outlines a procedure for simulating the flow field within multistage turbomachinery which includes the effects of unsteadiness, compressibility, and viscosity. The associated modeling equations are the average passage equation system which governs the time-averaged flow field within a typical passage of a blade row embedded within a multistage configuration. The results from a simulation of a low aspect ratio stage and a one-half turbine will be presented and compared with experimental measurements. It will be shown that the secondary flow field generated by the rotor causes the aerodynamic performance of the downstream vane to be significantly different from that of an isolated blade row.

Description: The effect of periodic rotor wakes on the secondary flow structure in a turbine stator cascade was investigated. A mechanism simulated the wakes shed from rotor blades by passing cylindrical rods across the inlet to a linear cascade installed in a recirculating water flow loop. Velocity measurements showed a passage vortex, similar to that seen in steady flow, during the time associated with undisturbed fluid. However, as the rotor wake passed through the blade row, a large crossflow toward the suction surface was observed in the midspan region. This caused the development of two large areas of circulation between the midspan and endwall regions, significantly distorting and weakening the passage vortices.

Description: A finite-volume lower-upper (LU) implicit scheme is used to simulate an inviscid flow in a turbine cascade. This approximate factorization scheme requires only the inversion of sparse lower and upper triangular matrices, which can be done efficiently without extensive storage. As an implicit scheme it allows a large time step to reach the steady state. An interactive grid generation program (TURBO), which is being developed, is used to generate grids. This program uses the control point form of algebraic grid generation which uses a sparse collection of control points from which the shape and position of coordinate curves can be adjusted. A distinct advantage of TURBO compared with other grid generation programs is that it allows the easy change of local mesh structure without affecting the grid outside the domain of dependence. Sample grids are generated by TURBO for a compressor rotor blade and a turbine cascade. The turbine cascade flow is simulated by using the LU implicit scheme on the grid generated by TURBO.

Description: Performance data of 0.17-scale model counterrotation pusher propeller configurations were taken in the NASA Lewis 8- by 6-Foot Supersonic Wind Tunnel at Mach numbers of 0.66, 0.71, 0.75, and 0.79. These tests investigated the aerodynamic performance of the unducted fan (UDF) demonstrator propeller engine developed in a joint program by General Electric and NASA. Data were recorded to show the effect on counterrotation propeller cruise efficiency of two takeoff noise-reduction concepts. These two concepts are reduced aft blade diameter and increased forward blade number. The four configurations tested were a baseline (F1/A1 8/8) configuration, a reduced aft diameter (F1/A3 8/8) configuration, an increase forward blade number (F1/A1 9/8) configuration, and a combination of the latter two (F1/A3 9/8) configurations. Data were collected with a complex counterrotation propeller test rig via rotating thrust and torque balances and pressure instrumentation. Data comparisons documented the power differences between the baseline and the reduced aft diameter concepts. Performance comparisons to the baseline configuration showed that reducing the aft blade diameter reduced the net efficiency, and adding a blade to the front rotor increased the net efficiency. The combination of the two concepts showed only slightly lower net efficiency than the baseline configuration. It was also found that the counterrotation demonstrator propeller model (F7/A7 8/8) configuration outperformed the baseline (F1/A1 8/8) configuration.

Description: The state of statistical turbulence modeling is reviewed with emphasis on compressible aerodynamic flows. A three-stage modeling process is suggested to act as a frame-work for improving turbulence models for complex flows. This process is then used to guide assessments of the comparative performance of various levels of turbulence models for several examples of flow fields of aerodynamic interest. The examples cover flows in two and three dimensions and with a broad range of Mach numbers. In addition to considerations of individual flow fields, some generalizations are drawn from the aggregate behavior of the comparisons.

Description: A cooperative program was undertaken by research organizations in England, France, Australia and the U.S. to study the capabilities of computational fluid dynamics codes (CFD) to predict the aerodynamic loading on helicopter rotor blades. The program goal is to compare predictions with experimental data for flight tests of a research Puma helicopter with rectangular and swept tip blades. Two topics are studied. First, computed results from three CFD codes are compared for flight test cases where all three codes use the same partial inflow-angle boundary conditions. Second, one of the CFD codes (FPR) is iteratively coupled with the CAMRAD/JA heilcopter performance code. These results are compared with experimental data and with an uncoupled CAMRAD/JA solution. The influence of flow field unsteadiness is found to play an important role in the blade aerodynamics. Alternate boundary conditions are suggested in order to properly model this unsteadiness in the CFD codes.

Description: Stability and transition experiments are conducted on a 45 degree sweep airfoil in the ASU Unsteady Wind Tunnel. Combined flow-visualization techniques and hot-wire measurements are used in conjuction with stability-code calculations to map out transition behavior for this flow geometry. Both steady and unsteady crossflow vortices are observed and at this time there is no contradiction with the theoretical predictions.

Description: Wind tunnel experiments were performed on a 70 deg sweep delta wing to determine the effect of a sinusoidal pitching motion on the pressure field on the suction side of the wing. Twelve pressure taps were placed from 35 to 90 percent of the chord, at 60 percent of the local semi-span. Pressure coefficients were measured as a function of Reynolds number and pitch rate. The pressure coefficient was seen to vary at approximately the same frequency as the pitching frequency. The relative pressure variation at each chord location was comparable for each case. The average pressure distribution through each periodic motion was near the static distribution for the average angle of attack. Upon comparing the upstroke and downstroke pressures for a specific angle of attack, the downstroke pressures were slightly larger. Vortex breakdown was seen to have the most significant effect at the 40 to 45 percent chord location, where a decrease in pressure was apparent.

Description: Goertler vortices are thought to be the cause of transition in many fluid flows of practical importance. A review of the different stages of vortex growth is given. In the linear regime, nonparallel effects completely govern this growth, and parallel flow theories do not capture the essential features of the development of the vortices. A detailed comparison between the parallel and nonparallel theories is given and it is shown that at small vortex wavelengths, the parallel flow theories have some validity; otherwise nonparallel effects are dominant. New results for the receptivity problem for Goertler vortices are given; in particular vortices induced by free stream perturbations impinging on the leading edge of the walls are considered. It is found that the most dangerous mode of this type can be isolated and it's neutral curve is determined. This curve agrees very closely with the available experimental data. A discussion of the different regimes of growth of nonlinear vortices is also given. Again it is shown that, unless the vortex wavelength is small, nonparallel effects are dominant. Some new results for nonlinear vortices of 0(1) wavelengths are given and compared to experimental observations.

Description: The large eddy interaction model (LEIM) is a statistical model of turbulence based on the interaction of selected eddies with the mean flow and all of the eddies in a turbulent shear flow. It can be utilized as the starting point for obtaining physical structures in the flow. The possible application of the LEIM to a mixing layer formed between two parallel, incompressible flows with a small temperature difference is developed by invoking a detailed similarity between the spectra of velocity and temperature.

Description: An investigation was conducted in the Langley Spin Tunnel to measure the pressures on the surface of the horizontal and vertical tail and the aft fuselage of an aircraft model. The pressures were measured on a model of a proposed Australian Primary Trainer airplane configuration while the model was rotating at spinning attitudes. The test results indicate that the presence of the horizontal tail adversely modifies the surface pressure on the vertical tail. The presence of the wing also adversely modifies the pressures on the horizontal and vertical tails.

Description: Recent developments in the state of the art in missile aerodynamics are reviewed. Among the subjects covered are: (1) Tri-service/NASA data base, (2) wing-body interference, (3) nonlinear controls, (4) hypersonic transition, (5) vortex interference, (6) airbreathers, supersonic inlets, (7) store separation problems, (8) correlation of missile data, (9) CFD codes for complete configurations, (10) engineering prediction methods, and (11) future configurations. Suggestions are made for future research and development to advance the state of the art of missile aerodynamics.

Description: A three-dimensional viscous-inviscid interaction analysis has been developed to predict the performance of rotors in hover and forward flight at subsonic and transonic tip speeds. The analysis solves the full-potential and boundary-layer equations by finite-difference numerical procedures. Calculations were made for several different model rotor configurations in hover and forward flight at subsonic and transonic tip speeds. The results were compared with predictions from a two-dimensional integral method and with experimental data. The comparisons show good agreement between test data and predictions.

Description: Flight boundary-layer transition experiments were conducted on a 30-degree swept wing with a perforated leading-edge suction panel. The transition location on the panel was changed by systematically varying the location and amount of suction. Transition from laminar to turbulent flow was due to leading-edge turbulence contamination or crossflow disturbance growth and/or Tollmien-Schlichting disturbance growth, depending on flight condition and suction variation. Amplification factor correlations with transition location were made for various suction configurations using a state-of-the-art linear stability theory which accounts for body and streamline curvature and compressibility.

Description: Navier-Stokes computations of three-dimensional vortical flows over a round-edged double-delta wing and a tangent-ogive forebody are performed using an implicit upwind finite-difference scheme. Computed particle traces are compared with experimental oil-flow streaks.

Description: Turbulence modeling for high speed compressible flows is described and discussed. Starting with the compressible Navier-Stokes equations, methods of statistical averaging are described by means of which the Reynolds-averaged Navier-Stokes equations are developed. Unknown averages in these equations are approximated using various closure concepts. Zero-, one-, and two-equation eddy viscosity models, algebraic stress models and Reynolds stress transport models are discussed. Computations of supersonic and hypersonic flows obtained using several of the models are discussed and compared with experimental results. Specific examples include attached boundary layer flows, shock wave boundary layer interactions and compressible shear layers. From these examples, conclusions regarding the status of modeling and recommendations for future studies are discussed.