ALBERT

Description: The Earth's atmosphere from 90 km to 200 km provides the last aerothermodynamics frontier. Present NASA programs which require but also can provide an understanding of the aerodynamics and aerothermodynamics of the free molecule and transition flows that exist at these altitudes are the Aeroassisted OTV, Entry Research Vehicle and the Tethered Satellite. Each of these programs provides a unique opportunity to do flight research in the rarefied upper atmosphere. However, the Tethered Satellite Program provides, because of its capability to obtain global, in-situ, steady state data, the greatest potential to: (1)define the performance of aerodynamic shapes as a function of environmental characteristics (free molecule, transition, slip flow regimes); (2)define the characteristics of the upper atmosphere and the global variability of properties such as composition temperature, pressure and density. Such data are required to accomplish the systematic development and verification of analytical prediction techniques required to support advance configuration designs.

Description: The prediction of missile aerodynamic characteristics is presently undertaken through the application of supersonic paneling methods and nonlinear corrections to the prediction of missile aerodynamic characteristics. Attention is given to supersonic panel methods and line-singularity methods for the modeling of axisymmetric bodies, in combination with corrections for nonlinear flow phenomena, which are applied to complete missile, inlets, and wing-body combinations. The LRCDM2 computer program is used as an example of the methods presented.

Description: The discrete vortex cloud approach models a missile airframe's vortex wake by means that are capable of treating a variety of configurations over a range of flow conditions. Attention is given to the sheets of vorticity formed on the lee side of a missile at moderate angles-of-attack. While three-dimensional attached flow models are used to represent the missile body, two-dimensional, incompressible, separated flow models are used to represent the separated vortex wake. The predicted pressure distribution of the body under the influence of the freestream and the separation vortex wake are used to calculate aerodynamic loads on the body. The separation vortex wake is represented by clouds of discrete vortices in cross flow planes normal to the body axis.

Description: Applications of laser velocimetry to the measurement of turbulent flow properties of strong transonic viscous-inviscid interactions are reviewed. The data resulting from these studies are then discussed in relation to their importance in the development of improved viscous-flow calculation methods. Also considered are the current limitations of laser velocimetry, the need for further improvements in the method, and potential future applications.

Description: The results of Reynolds-averaged time-dependent inviscid and turbulent compressible Navier-Stokes computations using the implicit finite-difference approach of Steger (1978), modified by incorporating a pressure boundary condition, (PBC) to account for wall interference are compared with experimental data on a NACA 64A010 airfoil (Johnson and Bachalo, 1980) in graphs and briefly characterized. The computational approach is the same as that used by King and Johnson (1980), but a 137 x 50 mesh is used instead of a 97 x 35 mesh, and special care is taken in resolving the nose, shock, and trailing-edge regions. Imposition of PBC is shown to improve significantly the accuracy of the computations for the flowfield on the upper surface of the airfoil, shifting the shock forward to its experimentally measured position in the case of turbulent flow. The failure of the method, even with PBC, to match the experimental shock location in the case of a flow with a separation bubble is attributed to inadequacies in the algebraic turbulence model employed (Baldwin and Lomax, 1978).

Description: Current progress in the computational analysis of rotary-wing flowfields is surveyed, and some typical results are presented in graphs. Topics examined include potential theory, rotating coordinate systems, lifting-surface theory (moving singularity, fixed wing, and rotary wing), panel methods (surface singularity representations, integral equations, and compressible flows), transonic theory (the small-disturbance equation), wake analysis (hovering rotor-wake models and transonic blade-vortex interaction), limitations on computational aerodynamics, and viscous-flow methods (dynamic-stall theories and lifting-line theory). It is suggested that the present algorithms and advanced computers make it possible to begin working toward the ultimate goal of turbulent Navier-Stokes calculations for an entire rotorcraft.

Description: Unsteady aerodynamic and aeroelastic stability calculations based upon transonic small disturbance (TSD) potential theory are presented. Results from the two-dimensional XTRAN2L code and the three-dimensional XTRAN3S code are compared with experiment to demonstrate the ability of TSD codes to treat transonic effects. The necessity of nonisentropic corrections to transonic potential theory is demonstrated. Dynamic computational effects resulting from the choice of grid and boundary conditions are illustrated. Unsteady airloads for a number of parameter variations including airfoil shape and thickness, Mach number, frequency, and amplitude are given. Finally, samples of transonic aeroelastic calculations are given. A key observation is the extent to which unsteady transonic airloads calculated by inviscid potential theory may be treated in a locally linear manner.

Description: An algebraic procedure for the generation of boundary-fitted grids about wing-fuselage configurations is presented. A wing-fuselage configuration is specified by cross sections and mathematically represented by Coons' patches. A configuration is divided into sections so that several grid blocks that either adjoin each other or partially overlap each other can be generated. Each grid has six exterior surfaces that map into a computational cube. Grids are first determined on the six boundary surfaces and then in the interior. Grid curves that are on the surface of the configuration are derived from the intersection of planes with the Coons' patch definition. Single-valued functions relating approximate arc lengths along the grid curves to a computational coordinate define the distribution of grid points. The two-boundary technique and transfinite interpolation are used to determine the boundary surface grids that are not on the configuration, and transfinite interpolation with linear blending functions is used to determine the interior grid.

Description: In the present treatment of the calculation of forces on a wing that is suddenly brought into motion at a constant speed, attention is given to the unsteady potential's contribution to the force balance. Total bound vorticity is produced at the initial impulse. The results obtained are independent of wing aspect ratio; as time increases, this effect on the drag force becomes smaller as the vortex emanating from the trailing edge is left behind. The second contributor to induced drag is the spanwise vorticity shedding that results from the spanwise load distribution of three-dimensional wings. This contribution grows with time as the length of the wake grows.

Description: Laminar heating distributions have been measured on a 1.9 percent scale model of a generic aeroassisted vehicle taking the shape of a spherically blunted, 13-deg/7-deg biconic whose forecone section is bent upward (by 7 deg) to furnish self-trim capability at a 20-deg angle-of-attack. The results thus obtained were compared with data gathered for a straight biconic. While no Reynolds number effect on heating was noted on the windward side of the forecone, the opposite was true of the leeward side, where a Reynolds number increase caused circumferential flow separation at lower angles of attack. Generally, windward heating was predicted to within 10 percent with a computer code solving the steady, three-dimensional parabolized Navier-Stokes equations.

Description: Flow and sound field data are presented for a 2.54 cm diameter air jet at a Mach number of 0.50 and a Reynolds number of 300,000. Distributions of mean velocity, turbulence intensities, Reynolds stress, spectral components of turbulence as well as of the near field pressure, together with essential characteristics of the far field sound are reported. This detailed set of data for one particular flow, erstwhile unavailable in the literature, is expected to help promoote and calibrate subsonic jet noise theories. 'Source locations' in terms of the turbulence maxima, coupling between the entrainment dynamics and the near pressure field, the sound radiation paths, and the balance in mass, momentum and sound energy fluxes are discussed. The results suggest that the large scale coherent structures of the jet govern the 'source locations' by controlling the turbulence and also strongly influence the near field pressure fluctuations.

Description: The CRAY multitasking system was developed in order to utilize all four processors and sharply reduce the wall clock run time. This paper describes the techniques used to modify the computational fluid dynamics code ARC3D for this run and analyzes the achieved speedup. The ARC3D code solves either the Euler or thin-layer N-S equations using an implicit approximate factorization scheme. Results indicate that multitask processing can be used to achieve wall clock speedup factors of over three times, depending on the nature of the program code being used. Multitasking appears to be particularly advantageous for large-memory problems running on multiple CPU computers.

Description: The NCOREL full-potential method with an entropy correction is presently applied to supersonic missile flowfield problems. After defining the salient characteristics of the method, a combination of linear theory with NCOREL and experimental data is used to isolate the nonlinear features of the supersonic flow so that the influence of geometry and flow conditions on the development of such flow nonlinearities can be appreciated. Comparisons of experimental longitudinal force and moment data with NCOREL and various linear theory predictions are presented for several generic missile airframe configurations of circular and elliptic cross section. The NCOREL code solves the nonconservative full potential equation in a spherical coordinate system; exact boundary conditions are defined on the missile surface.

Description: A comprehensive evaluation is made of experimental data compiled to date for the flowfields and aerodynamic forces that occur at high angles of attack for low aspect ratio wings with delta, rectangular, clipped delta, and strake/wing planform geometries. Attention is given to wing leading edge-generated vortex breakdown, aspect ratio and compressibility effects, and strake vortex effects on main wing areas. Although the nonlinear effects created by a wing-body combination significantly alter wing-alone aerodynamics, the wing-alone data presented are vital to the development of prediction methodologies for large angle of attack aerodynamics.

Description: A solution procedure is presented for the lifting transonic flow past modern rotor configurations in forward flight. In this procedure, the three-dimensional, unsteady Euler equations are solved in strong conservation form on a body-fitted moving coordinate system. A hybrid procedure of second order spatial accuracy and first order temporal accuracy is used to integrate the governing equations. In lifting flows, the effect of the elements of wake not captured by the computational procedure, and other aeroelastic effects are accounted for as local angle of attack corrections. Detailed comparisons with experimental data are presented for a 1/7 scale model of the Cobra OLS rotor, and for a three-bladed rotor tested in France. Some preliminary results are also presented for a three-dimensional blade vortex interaction problem.

Description: Rotary-wing computational fluid dynamics is reaching a point where many three-dimensional, unsteady, finite-difference codes are becoming available. This paper gives a brief review of five such codes, which treat the small disturbance, conservative and nonconservative full-potential, and Euler flow models. A discussion of the methods of applying these codes to the rotor environment (including wake and trim considerations) is followed by a comparison with various available data. These data include tests of advancing lifting and nonlifting, and hovering model rotors with significant supercritical flow regions. The codes are also compared for computational efficiency.

Description: The propagation characteristics of several helicopter airfoil profiles have been investigated using the transonic small disturbance equation. A test case was performed to generate a moving shock that propagated off the airfoil. Various grids were then examined to determine their ability to accurately capture these propagating shock waves. Finally, the case of airfoil-vortex interactions was thoroughly studied over a wide range of Mach numbers and airfoil shapes with particular emphasis on the transonic regime; this results in a highly conplicated fluctuation of lift, drag, and pitching moment. The calculated acoustic intensity levels, along with the details of the computational flow field, provide new insights into the understanding of transonic airfoil-vortex interactions.

Description: An experimental and theoretical investigation of rotor/wing aerodynamic interactions in hover is described. The experimental investigation consisted of both a large-scale and small-scale test. A 0.658-scale, V-22 rotor and wing was used in the large-scale test. Wind download, wing surface pressure, rotor performance, and rotor downwash data from the large-scale test are presented. A small-scale experiment was conducted to determine how changes in the rotor/wing geometry affected the aerodynamic interactions. These geometry variations included the distance between the rotor and wing, wing incidence angle, and configurations both with the rotor axis at the tip of the wing (tilt rotor configuration) and with the rotor axis at the center of the wing (compound helicopter oonfiguration). A wing with boundary-layer control was also tested to evaluate the effect of leading and trailing edge upper surface blowing on the wing download. A computationally efficient, semi-empirical theory was developed to predict the download on the wing. Finally, correlations between the theoretical predictions and test data are presented.

Description: It is proposed that the study of Rusak et al. (1985), which reports numerical modeling sensitivities on longitudinal force/moment properties for a vortex-lattice method incorporating free vortex filaments to represent the leading-edge vortex separation, employs a formula that is strongly affected by the particular points of analysis chosen. This results in a narrowly applicable curve fit, where numerical sensitivities of the theory are inappropriately traded off against physical effects that are not modeled in that theory. Attention is also given to questionable drag estimate computations.

Description: Steady, high speed, compressible separated flows modeled through numerical simulations resulting from solutions of the mass-averaged Navier-Stokes equations are reviewed. Emphasis is placed on benchmark flows that represent simplified (but realistic) aerodynamic phenomena. These include impinging shock waves, compression corners, glancing shock waves, trailing edge regions, and supersonic high angle of attack flows. A critical assessment of modeling capabilities is provided by comparing the numerical simulations with experiment. The importance of combining experiment, numerical algorithm, grid, and turbulence model to effectively develop this potentially powerful simulation technique is stressed.

Description: The analysis and the incorporation into a multigrid scheme of several vectorizable algorithms are discussed. von Neumann analyses of vertical-line, horizontal-line, and alternating-direction ZEBRA algorithms were performed; and the results were used to predict their multigrid damping rates. The algorithms were then successfully implemented in a transonic conservative full-potential computer program. The convergence acceleration effect of multiple grids is shown, and the convergence rates of the vectorizable algorithms are compared with those of standard successive-line overrelaxation (SLOR) algorithms.

Description: Numerical procedures for solving the thin-shear-layer Navier-Stokes equations and for the interaction of solutions to inviscid and boundary-layer equations are described and evaluated. To allow appraisal of the numerical and fluid dynamic abilities of the two schemes, they have been applied to one airfoil as a function of angle of attack at two slightly different Reynolds numbers. The NACA 0012 airfoil has been chosen because it allows comparison with measured lift, drag, and moment and with surface-pressure distributions. Calculations have been performed with algebraic eddy-viscosity formulations, and they include consideration of transition. The results are presented in a form that allows easy appraisal of the accuracy of both procedures and of the relative costs. The interactive procedure is computationally efficient but restrictive relative to the thin-layer Navier-Stokes procedure. The latter procedure does a better job of predicting drag than does the former. In both procedures, the location of transition is crucial for accurate or detailed computations, particularly at high angles of attack. When the upstream influence of pressure field through the shear layer is important, the thin-layer Navier-Stokes procedure has an edge over the interactive procedure.

Description: Recent progress in the development of finite element methodology for the prediction of aerothermal loads is described. Three dimensional, inviscid computations are presented, but emphasis is placed on development of an approach extendable to three dimensional viscous flows. Progress in key research areas is described. Initial 30 results from the computational procedure are described.

Description: The numerical method developed by Schiff and Sturek (1980) on the basis of the thin-layer parabolized Navier-Stokes equations of Schiff and Steger (1980) is extended to the case of turbulent supersonic flows on pointed bodies at high angles of attack. The governing equations, the numerical scheme, and modifications to the algebraic eddy-viscosity turbulence model are described; and results for three cones and one ogive-cylinder body (obtained using grids of 50 nonuniformly spaced points in the radial direction between the body and the outer boundary) are presented graphically and compared with published experimental data. The grids employed are found to provide sufficient spatial resolution of the leeward-side vortices; when combined with the modified turbulence model, they are shown to permit accurate treatment of flows with large regions of crossflow separation.

Description: The present use of a parabolized Navier-Stokes solver to accurately simulate the flowfield in a supersonic inlet yields good agreement between numerical analysis and experiment for a Mach 7.4 inlet under cruise conditions, with an internal compression ratio of 8. The significance of real gas effects on the performance calculation of a hypersonic inlet is demonstrated, with small changes in the ratio of specific heats resulting in a substantial change in the calculated pitot pressure ratio.

Description: The mechanism that locates a shock wave in a transonic flow in one and two dimensions is examined. It is found that in one dimension the shock is located by specifying the downstream pressure whereas in two dimensions the shock is located by the application of an entropy condition at the sonic line.

Description: Experimental data describing the transonic, turbulent, separated flow generated by an axisymmetric flow model are presented. The model consisted of a circular-arc bump affixed to a straight, circular cylinder aligned with the flow direction. Measurements of the mean velocity, turbulence intensity, and Reynolds shear-stress profiles were made in the separated flow. These data revealed dramatic changes in the shear-stress levels as the flow passed through the interaction to reattachment. Behavior of the turbulence reaction to the imposed pressure gradients was examined in terms of the mixing length and the excursions of the turbulence from equilibrium.

Description: An experimental investigation of rotor/wing aerodynamic interactions in hover is described. The investigation consisted of both a large-scale and a small-scale test. A 0.658-scale V-22 rotor and wing was used in the large-scale test. Wing download, wing surface pressure, rotor performance, and rotor downwash data from the large-scale test are presented. A small-scale experiment was conducted to determine how changes in the rotor/wing geometry affected the aerodynamic interactions. These geometry variations included the distance between the rotor and wing, wing incidence angle, wing flap angle, rotor rotation direction, and configurations both with the rotor axis at the tip of the wing (tilt rotor configuration) and with the rotor axis at the center of the wing (compound helicopter configuration).

Description: Detailed surface heat transfer and pressure distributions have been obtained in three-dimensional shock-wave boundary-layer interactions flow regions. The data described were obtained on fundamental shapes: planar wings with trailing edge flaps or spoilers and planar or cylindrical center bodies, representative of the aft portion of hypersonic aircraft. An overview of the work is presented; details of the projects are available in many reports in the open literature. Analytic, empiric methods are advanced for predicting the extent of separation and the increased heat transfer and pressure loads in three-dimensional separated flow regions.

Description: Experimental data are shown for the inception lengths of swept compression corner-generated and fin-generated shock/boundary layer interactions at Mach 2.95. These results are found to correlate on the basis of three different flow regimes. The inception lengths of these flows are dominated by a singularity at the cylindrical/conical boundary for swept corners and by an elongation due to shock wave sweepback for fin interactions. Similarity rules for both Re(delta)and shock generator geometry effects on inception lengths are demonstrated.

Description: A combined experimental and computational investigation of an axisymmetric turbulent shock-wave boundary-layer interaction flow is presented. Experimental measurements include both mean and fluctuating data obtained by LDV techniques and identify large scale unsteady motions associated with shock induced separation. Computations using the compressible Navier-Stokes equations, and a two-equation turbulence model are in relatively good agreement with experimental measurements. It is found that the large scale unsteady motions do not appear to have a critical impact on the ability to compute the mean properties of the flows investigated in this paper.

Description: A mathematically simple, turbulence closure model designed to treat transonic airfoil flows even with massive separation is described. Numerical solutions of the Reynolds-averaged, Navier-Stokes equations obtained with this closure model are shown to agree well with experiments over a broad range of test conditions.

Description: An account is given of methods for the estimation of a wing-body-tail missile configuration's aerodynamic performance by means of the 'component buildup' method, in which the overall aerodynamic loads for an airframe are built up from the assumed characteristics of isolated airframe components and then subjected to additional loads generated by component interference. Attention is given to the cases of missile airframes in steady flow at arbitrary angles of attack and bank; the unifying thread of the present treatment is slender body theory, together with its nonlinear extensions through the equivalent angle-of-attack concept. The estimation of the forces and moments acting on each of the fins is emphasized, so that control cross-coupling effects as well as longitudinal and lateral directional effects can be determined.

Description: This paper contributes to the understanding of the noise generation process of an airfoil encountering an unsteady upwash. By using a fast Fourier transform together with accurate airfoil response functions, the lift-time waveform for an airfoil encountering a delta function gust (the indicial function) is calculated for a flat plate airfoil in a compressible flow. This shows the interesting property that the lift is constant until the generated acoustic wave reaches the trailing edge. Expressions are given for the magnitude of this constant and for the pressure distribution on the airfoil during this time interval. The case of an airfoil cutting through a line vortex is also analyzed. The pressure-time waveform in the far field is closely related to the left-time waveform for the above problem of an airfoil entering a delta function gust. The effects of varying the relevant parameters in the problem are studied, including the observed position, the core diameter of the vortex, the vortex orientation and the airfoil span. The far field sound varies significantly with observer position, illustrating the importance of non-compactness effects. Increasing the viscous core diameter tends to smooth the pressure-time waveform. For small viscous core radius and infinite span, changing the vortex orientation changes only the amplitude of the pressure-time waveform, and not the shape.

Description: Far-field boundary conditions for the Euler equations are formulated and applied to transonic lifting flow over an airfoil in an unbounded domain. An expansion of the linearized small-disturbance equation in the far field is developed and the leading-order term, corresponding to a point vortex representation for the airfoil, is retained. A comprehensive evaluation across the Mach number range of the procedure's effectiveness in eliminating dependence of the numerical results on the boundary extent is presented. Extension of the method to three dimensions is also outlined.

Description: Flow through a two-dimensional duct with supersonic inflow is numerically investigated, from the viewpoint of the formation of Mach reflection, aerodynamic choking, and the possibility of constructing a curve similar to that for the quasi-one-dimensional flow in a converging-diverging duct. Such a curve can be used to determine whether a duct with a certain area ratio will or will not choke for a given inflow Mach number. Plots of pressure and mass flux contours are obtained for a given duct configuration. It is found that the two-dimensional flow always chokes at a higher Mach number than the corresponding quasi-one-dimensional flow for a given throat/inlet flow area ratio.

Description: The DFVLR laser Doppler anemometer is a CO2 continuous wave homodyne system designed for boundary layer wind measurements. During the last three years, it was mainly used in the wake-vortex program at Frankfurt airport for determination of vortex strength, transport, and lifetime. The strategy for that special type of measurement was previously reported in detail along with single experimental results. Therefore, herein is given a short summary of the data concerning questions of air traffic control. In addition to the experimental activities a computer model describing wake-vortex behavior was installed. It allows the comparison of the measured data with the hydrodynamically predicted quantities. On the other hand, it leads to an improved procedure for future wake-vortex measurements.

Description: A water tunnel flow visualization test on leading edge vortex flaps was conducted at the flow visualization facility of the NASA Ames Research Center's Dryden Flight Research Facility. The purpose of the test was to visually examine the vortex structures caused by various leading edge vortex flaps on the delta wing of an F-106 model. The vortex flaps tested were designed analytically and empirically at the NASA Langley Research Center. The three flap designs were designated as full-span gothic flap, full-span untapered flap, and part-span flap. The test was conducted at a Reynolds number of 76,000/m (25,000/ft). This low Reynolds number was used because of the 0.076-m/s (0.25-ft/s) test section flow speed necessary for high quality flow visualization. However, this low Reynolds number may have influenced the results. Of the three vortex flaps tested, the part-span flap produced what appeared to be the strongest vortex structure over the flap area. The full-span gothic flap provided the next best performance.

Description: For the past 3 years, a research program pertaining to the study of wing leading edge vortices at supersonic speeds has been conducted in the Fundamental Aerodynamics Branch of the High-Speed Aerodynamics Division at the Langley Research Center. The purpose of the research is to provide an understanding of the factors governing the formation and the control of wing leading-edge vortices and to evaluate the use of these vortices for improving supersonic aerodynamic performance. The studies include both experimental and theoretical investigations and focus primarily on planform, thickness and camber effects for delta wings. An overview of this research activity is presented.

Description: The Euler code FL057 was applied to a blunt nose smooth surface missile body shape. A range of angle of attacks was analyzed at Mach numbers of 0.55 and 2.0. A Mach number sweep from 0.55 to 2.0 was run for 12 degrees angle of attack. Experimental force and moment data were compared at Mach 2.0. The Euler code agreed with the experimental data over the linear portion of the Mach 0.55 data and over the entire angle-of-attack range at Mach 2.0.

Description: An experimental study was conducted to quantify the hysteresis associated with various vortex flow transition points and to determine the effect of planform geometry. The transition points observed consisted of the appearance (or disappearance) of trailing edge vortex burst and the transition to (or from) flat plate or totally separated flows. Flow visualization with smoke injected into the vortices was used to identify the transitions on a series of semi-span models tested in a low speed tunnel. The planforms tested included simple deltas (55 deg to 80 deg sweep), cranked wings with varying tip panel sweep and dihedral, and a straked wing. High speed movies at 1000 frames per second were made of the vortex flow visualization in order to better understand the dynamics of vortex flow, burst and transition.

Description: The vapor screen technique was successfully applied to an F-106B fighter aircraft during subsonic and transonic maneuvers. This system has allowed the viewing of multiple vortex systems on the wing upper surface at angles of attack less than 19 deg. In addition, similarities as well as differences were determined to exist between the vortex systems for a full scale semispan model and the flight vehicle at 20 deg incidence. Furthermore, variations in Reynolds number and Mach number were identified as to how they affect vortex system details at flight conditions.

Description: Some background information is provided for the Vortex Flow Aerodynamics Conference and that current slender wing airplanes do not use variable leading edge geometry to improve transonic drag polar is shown. Highlights of some of the initial studies combining wing camber, or flaps, with vortex flow are presented. Current vortex flap studies were reviewed to show that there is a large subsonic data base and that transonic and supersonic generic studies have begun. There is a need for validated flow field solvers to calculate vortex/shock interactions at transonic and supersonic speeds. Many important research opportunities exist for fundamental vortex flow investigations and for designing advanced fighter concepts.

Description: The need for experimentally determined 3-D velocity information is crucial to the understanding of highly 3-dimensional and vortical flow fields. In addition to gaining an understanding of the physics of flow fields, a correlation of velocity data is needed for advanced computational modelling. A double pass method for acquiring 3-D flow field information using a 2-D laser velocimeter (LV) is described. The design and implementation of a 3-D LV with expanded capabilities to acquire real-time 3-D flow field information are also described. Finally, the use of such an instrument in a wind tunnel study of a generic fighter configuration is described. The results of the wind tunnel study highlight the complexities of 3-D flow fields, particularly when the vortex behavior is examined over a range of angles of attack.

Description: The simulation of the leading edge vortex flow about a series of conical delta wings through solution of the Navier-Stokes and Euler equations is studied. The occurrence, the validity, and the usefulness of separated flow solutions to the Euler equations of particular interest. Central and upwind difference solutions to the governing equations are compared for a series of cross sectional shapes, including both rounded and sharp tip geometries. For the rounded leading edge and the flight condition considered, viscous solutions obtained with either central or upwind difference methods predict the classic structure of vortical flow over a highly swept delta wing. Predicted features include the primary vortex due to leading edge separation and the secondary vortex due to crossflow separation. Central difference solutions to the Euler equations show a marked sensitivity to grid refinement. On a coarse grid, the flow separates due to numerical error and a primary vortex which resembles that of the viscous solution is predicted. In contrast, the upwind difference solutions to the Euler equations predict attached flow even for first-order solutions on coarse grids. On a sufficiently fine grid, both methods agree closely and correctly predict a shock-curvature-induced inviscid separation near the leeward plane of symmetry. Upwind difference solutions to the Navier-Stokes and Euler equations are presented for two sharp leading edge geometries. The viscous solutions are quite similar to the rounded leading edge results with vortices of similar shape and size. The upwind Euler solutions predict attached flow with no separation for both geometries. However, with sufficient grid refinement near the tip or through the use of more accurate spatial differencing, leading edge separation results. Once the leading edge separation is established, the upwind solution agrees with recently published central difference solutions to the Euler equations.

Description: Two approaches to calculate turbulent vortical flows over delta wing configurations are illustrated. The first is for a simple delta wing at low speeds using the boundary layer approximation to treat the effects of the secondary separation. The second is for the supersonic case of a generic fighter using the NASA Ames parabolized Navier/Stokes method. Test/theory comparisons are given in both cases.

Description: The effectiveness of apex fences on a 60-deg delta wing at low speeds was experimentally investigated. Resembling highly swept spoilers in appearance, the fences are designed to fold out of the wing apex region upper surface near the leading edges, where they generate a powerful vortex pair. The intense suction of the fence vortices augments lift in the apex region, the resulting positive pitching moment being utilized to trim trailing edge flaps for lift augmentation during approach and landing at relatively low angles of attack. The fences reduce the apex lift at high angles of attack, leading to a desirable nose-down moment. The above projected functions of the apex fence device were validated and quantified through low speed tunnel tests, comprising upper surface pressure surveys on a semispan model and balance measurements on a geometrically similar fully span wing/body configuration. Fence parameters such as area, shape, hinge position and deflection angle were investigated. Typical results are presented indicating the apex fence potential in controlling the longitudinal characteristics of a tail-less delta.

Description: A new version of the free vortex sheet formulation is presented which has greatly improved convergence characteristics for a broad range of geometries. The enhanced convergence properties were achieved largely with extended modeling capabilities of the leading edge vortex and the near field trailing wake. Results from the new code, designated FVS-1, are presented for a variety of configurations and flow conditions with emphasis on vortex flap applications.

Description: Two different schemes are presented for including the effect of rotor wakes on the finie-difference prediction of rotor loads. The first formulation includes wake effects by means of a blade-surface inflow specification. This approach is sufficiently simple to permit coupling of a full-potential finite-difference rotor code to a comprehensive integral model for the rotor wake and blade motion. The coupling involves a transfer of appropriate loads and inflow data between the two computer codes. Results are compared with experimental data for two advancing rotor cases. The second rotor-wake modeling scheme is a split potential formulation for computing unsteady blade-vortex interactions. Discrete vortex fields are introduced into a three-dimensional, conservative, full-potential rotor code. Computer predictions are compared with two experimental blade-vortex interaction cases.

Description: Results are presented for the measured performance recently obtained on several airfoil concepts designed to achieve low drag by maintaining extensive regions of laminar flow without compromising high-lift performance. The wind tunnel results extend from subsonic to transonic speeds and include boundary-layer control through shaping and suction. The research was conducted in the NASA Langley 8-Ft Transonic Pressure Tunnel (TPT) and Low Turbulence Pressure Tunnel (LTPT) which have been developed for testing such low-drag airfoils. Emphasis is placed on identifying some of the major factors influencing the anticipated performance of low-drag airfoils.

Description: A review of the NACA and NASA low-drag airfoil research is presented with particular emphasis given to the development of mechanical high-lift flap systems and their application to general aviation aircraft. These flap systems include split, plain, single-slotted, and double-slotted trailing-edge flaps plus slat and Krueger leading-edge devices. The recently developed continuous variable-camber high-lift mechanism is also described. The state-of-the-art of theoretical methods for the design and analysis of multi-component airfoils in two-dimensional subsonic flow is discussed, and a detailed description of the Langley MCARF (Multi-Component Airfoil Analysis Program) computer code is presented. The results of a recent effort to design a single- and double-slotted flap system for the NASA high speed natural laminar flow (HSNLF) (1)-0213 airfoil using the MCARF code are presented to demonstrate the capabilities and limitations of the code.

Description: Two research studies are described which directly relate to the application of natural laminar flow (NLF) technology to transonic transport-type wing planforms. Each involved using state-of-the-art computational methods to design three-dimensional wing contours which generate significant runs of favorable pressure gradients. The first study supported the Variable Sweep Transition Flight Experiment and involves design of a full-span glove which extends from the leading edge to the spoiler hinge line on the upper surface of an F-14 outer wing panel. A wing was designed computationally for a corporate transport aircraft in the second study. The resulting wing design generated favorable pressure gradients from the leading edge aft to the mid-chord on both upper and lower surfaces at the cruise design point. Detailed descriptions of the computational design approach are presented along with the various constraints imposed on each of the designs.

Description: The Langley Research Center has a concentrated and directed effort under way to develop both conventional and non-intrusive diagnostic instrumentation. These instruments are being developed to operate over large Mach number, total temperature, and total pressure ranges. Efforts are being made to evaluate the measurements made by the various instruments to determine the most accurate and reliable instrument to be used under a given flow environment. Although only one flow visualization technique was described, there are many different types presently being used at Langley Research Center.

Description: A description of and results from a solution algorithm for the compressible Navier-Stokes equations are presented. The main features of the algorithm are second or third order accurate upwind discretization of the convection and pressure derivatives and a relaxation scheme for the unfactored implicit backward Euler time method, implemented in a finite-volume formulation. Upwind methods were successfully used to obtain solutions to the Euler equations for flows with strong shock waves. The particular upwind method being used is based on the flux vector splitting technique developed by Van Leer and both second and third order accurate discretizations were developed. Currently, the most widely used implicit solution technique for the Navier-Stokes equations use approximate factorization (AF) methods to treat multidimensional problems. The time integration scheme being used in the present algorithm corresponds to a line Gauss-Seidel relaxation method. This method produces good convergence rates for steady-state flows, and most of the algorithm was vectorized on the NASA Langley VPS 32 computer. The Navier-Stokes algorithm was tested for several two-dimensional flow problems. Solutions for the problems gave excellent results. The presented effort is directed toward the extension of the scheme to the full three-dimensional Navier-Stokes equations.

Description: An experimental verification of a high performance natural laminar flow (NLF) airfoil for low speed and high Reynolds number applications was completed in the Langley Low Turbulence Pressure Tunnel (LTPT). Theoretical development allowed for the achievement of 0.70 chord laminar flow on both surfaces by the use of accelerated flow as long as tunnel turbulence did not cause upstream movement of transition with increasing chord Reynolds number. With such a rearward pressure recovery, a concave type deceleration was implemented. Two-dimensional theoretical analysis indicated that a minimum profile drag coefficient of 0.0026 was possible with the desired laminar flow at the design condition. With the three-foot chord two-dimensional model constructed for the LTPT experiment, a minimum profile drag coefficient of 0.0027 was measured at c sub l = 0.41 and Re sub c = 10 x 10 to the 6th power. The low drag bucket was shifted over a considerably large c sub l range by the use of the 12.5 percent chord trailing edge flap. A two-dimensional lift to drag ratio (L/D) was 245. Surprisingly high c sub l max values were obtained for an airfoil of this type. A 0.20 chort split flap with 60 deg deflection was also implemented to verify the airfoil's lift capabilities. A maximum lift coefficient of 2.70 was attained at Reynolds numbers of 3 and 6 million.

Description: Several computational studies are currently being pursued that focus on various aspects of representing the entire lifetime of the viscous trailing vortex wakes generated by an aircraft. The formulation and subsequent near-wing development of the leading-edge vortices formed by a delta wing are being calculated at modest Reynolds numbers using a three-dimensional, time-dependent Navier-Stokes code. Another computational code was developed to focus on the roll-up, trajectory, and mutual interaction of trailing vortices further downstream from the wing using a two-dimensional, time-dependent, Navier-Stokes algorithm. To investigate the effect of a cross-wind ground shear flow on the drift and decay of the far-field trailing vortices, a code was developed that employs Euler equations along with matched asymptotic solutions for the decaying vortex filaments. And finally, to simulate the conditions far down stream after the onset of the Crow instability in the vortex wake, a full three-dimensional, time-dependent Navier-Stokes code was developed to study the behavior of interacting vortex rings.

Description: Separation-induced leading-edge vortices can dominate the flow about slender wings at moderate to high angles of attack, often with favorable aerodynamic effects. However, at the high angles of attack which are desirable for takeoff and landing as well as subsonic-transonic maneuver the vortices can breakdown or burst in the vicinity of the aircraft causing many adverse effects; these include lift loss, pitchup, and buffet. The flow in the core of leading-edge vortices is generally affiliated with the vortex breakdown phenomenon. A theory is presented for the flow in the core of separation-induced, leading-edge vortices at practical Reynolds numbers. The theory is based on matching inner and outer representations of the vortex. The inner representation models continuously distributed vorticity and includes an asymptotic viscous subcore. The outer representation models concentrated spiral sheets of vorticity and is fully three dimensional. A parameter is identified which closely tracks the vortex breakdown stability boundary for delta, arrow, and diamond wings.

Description: A finite-volume scheme for numerical integration of the Euler equations was extended to allow solution of the thin-layer Navier-Stokes equations in two and three dimensions. The extended algorithm, which is based on a class of four-stage Runge-Kutta time-stepping schemes, was made numerically efficient through the following convergence acceleration technique: (1) local time stepping, (2) enthalpy damping, and (3) residual smoothing. Also, the high degree of vectorization possible with the algorithm has yielded an efficient program for vector processors. The scheme was evaluated by solving laminar and turbulent flows. Numerical results have compared well with either theoretical or other numerical solutions and/or experimental data.

Description: The breakdown of the conservative potential approximation occurs for all the airfoils tested. It develops as soon as the shock waves appear in the flow field. Since shock waves are not properly represented by the potential approximation, it is conjectured that the breakdown is due to the isentropic shock jump condition of the potential approximation.

Description: A conservative finite-volume difference scheme is developed for the potential equation to solve transonic flow about airfoils and bodies in an arbitrarily shaped channel. The scheme employs a mesh which is a nearly conformal O mesh about the airfoil and nearly orthogonal at the channel walls. The mesh extends to infinity upstream and downstream, where the mapping is singular. Special procedures are required to treat the singularities at infinity, including computation of the metrics near those points. Channels with exit areas different from inlet areas are solved; a body with a sting mount is an example of such a case.

Description: Two research programs are described which directly relate to the application of natural laminar flow (NLF) technology to transonic transport-type wind planforms. Each involved using state-of-the-art computational methods to design three-dimensional wing contours which generate significant runs of favorable pressure gradients. The first program supported the Variable Sweep Transition Flight Experiment and involves design of a full-span glove which extends from the leading edge to the spoiler hinge line on the upper surface of an F-14 outer wing panel. Boundary-layer and static-pressure data will be measured on this design during the supporting wind-tunnel and flight tests. These data will then be analyzed and used to infer the relationship between crossflow and Tollmein-Schlichting disturbances on laminar boundary-layer transition. A wing was designed computationally for a corporate transport aircraft in the second program. The resulting wing design generated favorable pressure gradients from the leading edge aft to the mid-chord on both upper and lower surfaces at the cruise design point. Detailed descriptions of the computational design approach are presented along with the various constraints imposed on each of the designs. Wing surface pressure distributions, which support the design objective and were derived from transonic three-dimensional analyses codes, are also presented. Current status of each of the research programs is included in the summary.

Description: As a result of reductions in form drag and roughness drag, skin friction, drag, or viscous drag now represents a major contributor to the cruise drag of subsonic business and transport aircraft, and hence, is considered a barrier problem to further significant improvements in the aerodynamic efficiency of these aircraft. To meet the challenge, research in the areas of laminar-flow control and turbulence control/drag reduction was initiated at NASA Langley Research Center. The significance of this research is discussed.

Description: Vortex flows of interest to aerodynamicists cover a wide range of scales from a fraction of an inch in boundary layer flows to many feet in wake flows. In many applications these flows are poorly understood and, due to their complexity, present a challenge both analytically and experimentally. Four topics representing the spectrum of experimental and analytical vortex research are presented.

Description: The present effort represents a first attempt of numerical simulation of the flow field around a complete aircraft-like, lifting configuration utilizing the Reynolds averaged Navier-Stokes equations. The numerical solution generated for the experimental aircraft concept X24C-10D at a Mach number of 5.95 not only exhibited accurate prediction of detailed flow properties but also of the integrated aerodynamic coefficients. In addition, the present analysis demonstrated that a page structure of data collected into cyclic blocks is an efficient and viable means for processing the Navier-Stokes equations on the CRAY XMP-22 computer with external memory device.

Description: A fast diagonalized Beam-Warming algorithm is coupled with a zonal approach to solve the three dimensional Euler/Navier-Stokes equations. The computer code, called Transonic Navier-Stokes (TNS), uses a total of four zones for wing configurations (or can be extended to complete aircraft configurations by adding zones). In the inner blocks near the wing surface, the thin-layer Navier-Stokes equations are solved, while in the outer two blocks the Euler equations are solved. The diagonal algorithm yields a speedup of as much as a factor of 40 over the original algorithm/zonal method code. The TNS code, in addition, has the capability to model wind tunnel walls. Transonic viscous solutions are obtained on a 150,000-point mesh for a NACA 0012 wing. A three-order-of-magnitude drop in the L2-norm of the residual requires approximately 500 iterations, which takes about 45 min of CPU time on a Cray-XMP processor. Simulations are also conducted for a different geometrical wing called WING C. All cases show good agreement with experimental data.

Description: A review of past parabolized Navier-Stokes applications is presented. The equations, boundary conditions, the numerical method and the grid generation are all discussed. Results ranging from the low supersonic regime to the hypersonic regime are included.

Description: A computer code is under development whereby the thin-layer Reynolds-averaged Navier-Stokes equations are to be applied to realistic fighter aircraft configurations. This transonic Navier-Stokes code (TNS) utilizes a zonal approach in order to treat complex geometries and satisfy in-core computer memory constraints. The zonal approach was applied to isolated wing geometries in order to facilitate code development. The TNS finite difference algorithm, zonal methodology, and code validation with experimental data is addressed. Also addressed are some numerical issues such as code robustness, efficiency, and accuracy at high angles of attack. Special free-stream-preserving metrics proved an effective way to treat H-mesh singularities over a large range of severe flow conditions, including strong leading edge flow gradients, massive shock induced separation, and stall. Furthermore, lift and drag coefficients were computed for a wing up through CLmax. Numerical oil flow patterns and particle trajectories are presented both for subcritical and transonic flow. These flow simulations are rich with complex separated flow physics and demonstrate the efficiency and robustness of the zonal approach.

Description: Methods of analyzing and experimentally measuring the effect of ice accretion on airfoil sections are presented. Empirical and analytical methods for predicting airfoil performance degradation due to ice are discussed. Ice simulation techniques for aerodynamic testing are presented and compared to data with actual ice accretions. The results show that simulation techniques to imitate the effect of ice on airfoil performance work well in most cases. Comparisons between predicted and measured airfoil performance with ice accretions are presented. For rime ice cases, the predictions compared well with experiments; but for glaze ice, a need for improved methods are seen.

Description: An attempt is made to develop an efficient staggered cascade blade unsteady aerodynamics model for the neighborhood of March 1, representing the blade row by a rectilinear two-dimensional cascade of thin, flat plate airfoils. The equations of motion are derived on the basis of linearized transonic small perturbation theory, and an analytical solution is obtained by means of the Wiener-Hopf procedure. Making use of the transonic similarity law, the results obtained are compared with those of other linearized cascade analyses. A parametric study is conducted to find the effects of reduced frequency, stagger angle, solidity, and the location of the pitching axis on cascade stability.

Description: An overview of a new finite element method for the compressible Euler and Navier-Stokes equations is presented. The discretization is based on entropy variables. The method is developed within the framework of a Petrov-Galerkin formulation. Two perturbations are added to the weighting function; one is a generalization of the SUPG operator and the other is designed to enhance shock capturing capability. The treatment of boundary conditions and the consistent calculation of boundary fluxes are addressed. Results of numerical tests are presented which confirm the robustness and wide applicability of the method.