ALBERT

Description: The model is a replica of a modified MAX-103 kit aircraft that our Parks College of St. Louis University Student Design Group built and modified from a tail wheel to a tricycle configuration. A model was tested in the Parks College low-speed wind tunnel. I hope to initiate flight-testing upon my second return to. St. Louis. The combined data using wind tunnel, water tunnel, RC, flight-testing and analytical results will be very valuable for assessing the correlation between the different methods of analyses, since at present it is almost impossible to accurately predict flight characteristics from anything but in-situ tests. Unfortunately, political/financial reasons dictated using a generic wing rather than a specific model in the NASA-DFRC water tunnel.

Description: A hybrid boundary element finite volume method for unsteady transonic flow computation has been developed. In this method, the unsteady Euler equations in a moving frame of reference are solved in a small embedded domain (inner domain) around the airfoil using an implicit finite volume scheme. The unsteady full-potential equation, written in the same frame of reference and in the form of the Poisson equation. is solved in the outer domain using the integral equation boundary element method to provide the boundary conditions for the inner Euler domain. The solution procedure is a time-accurate stepping procedure, where the outer boundary conditions for the inner domain are updated using the integral equation -- boundary element solution over the outer domain. The method is applied to unsteady transonic flows around the NACA0012 airfoil undergoing pitching oscillation and ramp motion. The results are compared with those of an implicit Euler equation solver, which is used throughout a large computational domain, and experimental data.

Description: Large-eddy simulation (LES) has matured to the point where application to complex flows is described. The extension to higher Reynolds numbers leads to an impractical number of grid points with existing structured-grid methods. Furthermore, most real world flows are rather difficult to represent geometrically with structured grids. Unstructured-grid methods offer a release from both of these constraints. However, just as it took many years for structured-grid methods to be well understood and reliable tools for LES, unstructured-grid methods must be carefully studied before we can expect them to attain their full potential.

Description: The concept studied during the summer NASA/ASEE Fellowship provides a means of lowering drag and a means for directional control of supersonic and hypersonic vehicles. Low drag and efficient directional control are essential for the success of aircraft, atmospheric entry vehicles, missiles, and other vehicles in supersonic and hypersonic flight. Drag reduction can result in increased vehicle range, increased speed, improved fuel efficiency, increased lift/drag ratio, and increased climb rate. For high supersonic and hypersonic vehicles heat transfer considerations dictate the design of the nose and leading edge. The heat transfer to such vehicles is most severe at stagnation points which occur on the leading edges and nose of the vehicle. Theoretical formulations, experimental data, and semi-empirical formulas all agree in the fact that stagnation point heat transfer is inversely proportional to the square root of the nose or leading edge radius. Thus, the noses and leading edges of supersonic and hypersonic vehicles are typically blunted so that the heat transfer and structural loads will be manageable. However, much of the wave drag experienced by these vehicles is due to nose blunting.

Description: A theoretical and experimental program is underway at NASA Ames Research Center to first obtain a better understanding of the hazard posed by the vortex wakes of subsonic transports, and then to develop methods on how to modify the wake-generating aircraft in order to make the vortices less hazardous. This paper summarizes results obtained in the 80- by 120-Foot Wind Tunnel at NASA Ames Research Center on the characteristics of the vortex wakes that trail from 0.03 scale models of a B-747 and of a DC-10. Measurements are first described that were taken in the wakes with a hot-film anemometer probe, and with wings that range in size from 0.2 to 1.0 times the span of the wake generating models at downstream distances of 81 ft and 162 ft. behind the wake-generating model; i.e., at scale distances of 0.5 and 1.0 mile. The data are then used to evaluate the accuracy of a vortex-lattice method for prediction of the loads induced on following wings by vortex wakes.

Description: Results from two laboratory studies and two numerical studies are presented. In the first laboratory study, measurements of the strength of vortices from a three-dimensional (3-D) model wing are presented. The measurements follow the vortices as they evolve in time from a two-dimensional (2-D) line vortex pair to the development and migration of 3-D vortex rings. It is shown that the resulting vortex rings can contain up to 40 percent of the initial vortex circulation. Thus, the formation of vortex rings may not necessarily signal the end of the wake hazard to following aircraft. In the second laboratory study, we present the results of an experiment which shows how the spanwise drag distribution affects wake-vortex evolution. In this experiment, we modified the spanwise drag distribution on a model wing while keeping the total lift and drag constant. The results show that adding drag on or near the centerline of the wing has a larger effect than adding drag at or near the wingtips. These measurements complement the results of NASA studies in the 1970s. In the first numerical study, results of 3-D numerical calculations are presented which show that the vortex Reynolds number has a significant influence on the evolution and migration of wake vortices. When the Reynolds number is large, 3-D vortex rings evolve from the initially 2-D line vortex pairs. These vortex rings then migrate vertically. When the Reynolds number is lower, the transition of vorticity from 2-D to 3-D is delayed. When the Reynolds number is very low, the vortices never transition to 3-D, and the vertical migration is significantly reduced. It is suggested that this effect may have been important in previous laboratory wake-evolution studies. A second numerical study shows the influence that vertical wind shear can have on trailing vortex evolution.

Description: The lifting body/linear aerospike is one of three configurations being studied for a single stage to orbit (SSTO) vehicle. A preliminary aerodynamics database existed for then current lifting body configurations, however, this database was developed without considering plume effects. A combined effort by the Computational Fluid Dynamics (CFD) and the Experimental Fluids Dynamics Branches was undertaken to determine first order effects of plume/external flow interactions on vehicle aerodynamics of this lifting body/linear aerospike configuration. Of interest were plume pumping/entrainment at low Mach numbers and plume induced separation of flow over the vehicle at higher altitudes. The CFD analysis included combinations of four Mach numbers, two angles of attack, and four throttle settings. The majority of the CFD was two dimensional centerline analysis of the lifting body/aerospike. Incremental plume effects were derived by comparing the power-on, power-off, and throttled cases and were extrapolated to the preliminary aerodynamic database. The plume had little effect on the vehicle aerodynamics for supersonic freestream velocities. At subsonic freestream velocities, the plume affected the vehicle aerodynamics through both jet pumping/entrainment and the jet flap effect.

Description: Flow computations around two-element and three-element configurations are presented and compared to detailed experimental measurements. The k-epsilon-v(exp 2)(bar) model has been applied and the ability of the model to capture streamline curvature effects, wake-boundary layer confluence, and laminar/turbulent transition is discussed. The numerical results are compared to experimental datasets that include mean quantities (velocity and pressure coefficient) and turbulent quantities (Reynolds normal and shear stresses).

Description: The k-epsilon-v(sup 2)(bar) model has been investigated to quantify its predictive performance on two high-lift configurations: 2D flow over a single-element aerofoil, involving closed-type separation; 3D flow over a prolate spheroid, involving open-type separation. A 'code-friendly' modification has been proposed which enhances the numerical stability, in particular, for explicit and uncoupled flow solvers. As a result of introducing Reynolds-number dependence into a coefficient of the s-equation, the skin-friction distribution for the by-pass transitional flow over a flat plate is better predicted. In order to improve deficiencies arising from the Boussinesq approximation, a nonlinear stress-strain constitutive relation was adopted, in which the only one free constant is calibrated on the basis of DNS data, and the Reynolds-stress anisotropy near the wall is fairly well represented.

Description: Jet plume impingement forces acting on large flexible space structures may precipitate dynamically unstable behavior during space flights. Typical operating conditions in space involve rarefied gas flow regimes which are intrinsically distinct from continuum gas flow and are normally modeled using the kinetic theory of gas flow. Docking and undocking operations of the Space Shuttle with the Russian Mir space laboratory represent a scenario in which the stability boundaries of solar panels may be of interest. Extensive literature review of research work on the dynamic stability of rectangular panels in rarefied gas flow conditions indicated the lack of published reports dealing with this phenomenon. A recently completed preliminary study for NASA JSC dealing with the mathematical analysis of the stability of two-degree-of-freedom elastically supported rigid panels under the effect of rarefied gas flow was reviewed. A test plan outline is prepared for the purpose of conducting a series of experiments on four rectangular rigid test articles in a vacuum chamber under the effect of continuous and pulsating Nitrogen jet plumes. The purpose of the test plan is to gather enough data related to a number of key parameters to allow the validation of the two-degree-of-freedom mathematical model. The hardware required careful design to select a very lightweight material while satisfying rigidity and frequency requirements within the constraints of the test environment. The data to be obtained from the vacuum chamber tests can be compared with the predicted behavior of the theoretical two-degree-of-freedom model. Using the data obtained in this study, further research can identify the limitations of the mathematical model. In addition modifications to the mathematical model can be made, if warranted, to accurately predict the behavior of rigid panels under rarefied gas flow regimes.

Description: A proposed wing body reusable launch vehicle was tested in the Marshall Space Flight Center (MSFC) 14 x 14 inch trisonic wind tunnel during the winter of 1994. this test resulted in the vehicle's subsonic and transonic (Mach 03 to 1.96) longitudinal and lateral aerodynamic characteristics. The effects of control surface deflections on the basic vehicle aerodynamics including a body flap, elevons, ailerons, and tip fins are presented.

Description: To support the reusable launch vehicle concept study, the aerodynamic data and surface pressure for WB001 were predicted using three computational fluid dynamic (CFD) codes at several flow conditions between code to code and code to aerodynamic database as well as available experimental data. A set of particular solutions have been selected and recommended for use in preliminary conceptual designs. These computational fluid dynamic (CFD) results have also been provided to the structure group for wing loading analysis.

Description: This paper deals with nonplanar wing concepts -- their advantages and possible applications in a variety of aircraft designs. A brief review and assessment of several concepts from winglets to ring wings is followed by a more detailed look at two recent ideas: exploiting nonplanar wakes to reduce induced drag, and applying a 'C-wing' design to large commercial transports. Results suggest that potential efficiency gains may be significant, while several nonaerodynamic characteristics are particularly interesting.

Description: North American Marine Jet (NAMJ), Inc. received assistance from Marshall Space Flight Center engineers in the Computational Fluid Dynamics (CFD) branch of the Structure and Dynamics Laboratory in improving the proposed design of a new impeller for their jet-propulsion systems. Marshall used advanced CFD techniques, which included creating a three-dimensional computer model of the impeller for analysis. With Marshall input, the company modified the design, then Marshall used a computer model to make a solid polycarbonate model. The rapid prototyping allowed the company to avoid many time- consuming and costly steps in creating the impeller model. NAMJ is now able to compete with Pacific-area and European manufacturers who have traditionally dominated the market.

Description: In this paper we analyze inviscid aerodynamic shape optimization problems governed by the full potential and the Euler equations in two and three dimensions. The analysis indicates that minimization of pressure dependent cost functions results in Hessians whose eigenvalue distributions are identical for the full potential and the Euler equations. However the optimization problems in two and three dimensions are inherently different. While the two dimensional optimization problems are well-posed the three dimensional ones are ill-posed. Oscillations in the shape up to the smallest scale allowed by the design space can develop in the direction perpendicular to the flow, implying that a regularization is required. A natural choice of such a regularization is derived. The analysis also gives an estimate of the Hessian's condition number which implies that the problems at hand are ill-conditioned. Infinite dimensional approximations for the Hessians are constructed and preconditioners for gradient based methods are derived from these approximate Hessians.

Description: A numerical investigation was conducted to assess the accuracy of two turbulence models when computing non-axisymmetric nozzle-afterbody flows with propulsive jets. Navier-Stokes solutions were obtained for a Convergent-divergent non-axisymmetric nozzle-afterbody and its associated jet exhaust plume at free-stream Mach numbers of 0.600 and 0.938 at an angle of attack of 0 deg. The Reynolds number based on model length was approximately 20 x 10(exp 6). Turbulent dissipation was modeled by the algebraic Baldwin-Lomax turbulence model with the Degani-Schiff modification and by the standard Jones-Launder kappa-epsilon turbulence model. At flow conditions without strong shocks and with little or no separation, both turbulence models predicted the pressures on the surfaces of the nozzle very well. When strong shocks and massive separation existed, both turbulence models were unable to predict the flow accurately. Mixing of the jet exhaust plume and the external flow was underpredicted. The differences in drag coefficients for the two turbulence models illustrate that substantial development is still required for computing very complex flows before nozzle performance can be predicted accurately for all external flow conditions.

Description: Extensive pressure measurements and off-surface flow visualization were obtained on the forebody and strakes of the NASA F-18 High Alpha Research Vehicle (HARV) equipped with actuated forebody strakes. Forebody yawing moments were obtained by integrating the circumferential pressures on the forebody and strakes. Results show that large yawing moments can be generated with forebody strakes. At angles of attack greater than 40 deg., deflecting one strake at a time resulted in a forebody yawing moment control reversal for small strake deflection angles. At alpha = 40 deg. and 50 deg., deflecting the strakes differentially about a 20 deg. symmetric strake deployment eliminated the control reversal and produced a near linear variation of forebody yawing moment with differential strake deflection. At alpha = 50 deg. and for 0 deg. and 20 deg. symmetric strake deployments, a larger forebody yawing moment was generated by the forward fuselage (between the radome and the apex of the leading-edge extensions), than on the radome where the actuated forebody strakes were located. Cutouts on the flight vehicle strakes that were not on the wind tunnel models are believed to be responsible for deficits in the suction peaks on the flight radome pressure distributions and differences in the forebody yawing moments.

Description: An investigation was conducted in the static test facility of the Langley 16-Foot Transonic Tunnel to determine the internal performance of a two-dimensional convergent-divergent nozzle. The nozzle design was tested with dry and afterburning throat areas, which represent different power settings and three expansion ratios. For each of these configurations, three trailing-edge geometries were tested. The baseline geometry had a straight trailing edge. Two different shaping techniques were applied to the baseline nozzle design to reduce radar observables: the scarfed design and the sawtooth design. A flat plate extended downstream of the lower divergent flap trailing edge parallel to the model centerline to form a shelf-like expansion surface. This shelf was designed to shield the plume from ground observation (infrared radiation (IR) signature suppression). The shelf represents the part of the aircraft structure that might be present in an installed configuration. These configurations were tested at nozzle pressure ratios from 2.0 to 12.0.

Description: A detailed experimental investigation to understand and quantify the development of loss and blockage in the flow field of a transonic, axial flow compressor rotor has been undertaken. Detailed laser anemometer measurements were acquired upstream, within, and downstream of a transonic, axial compressor rotor operating at design and off-design conditions. The rotor was operated at 100%, 85%, 80%, and 60% of design speed which provided inlet relative Mach numbers at the blade tip of 1.48, 1.26, 1.18, and 0.89 respectively. At design speed the blockage is evaluated ahead of the rotor passage shock, downstream of the rotor passage shock, and near the trailing edge of the blade row. The blockage is evaluated in the core flow area as well as in the casing endwall region. Similarly at pm speed conditions for the cases of (1) where the rotor passage shock is much weaker than that at design speed and (2) where there is no rotor passage shock, the blockage and loss are evaluated and compared to the results at design speed. Specifically, the impact of the rotor passage shock on the blockage and loss development, pertaining to both the shock/boundary layer interactions and the shock/tip clearance flow interactions, is discussed. In addition, the blockage evaluated from the experimental data is compared to (1) an existing correlation of blockage development which was based on computational results, and (2) computational results on a limited basis. The results indicate that for this rotor the blockage in the endwall region is 2-3 times that of the core flow region and the blockage in the core flow region more than doubles when the shock strength is sufficient to separate the suction surface boundary layer. The distribution of losses in the care flow region indicate that the total loss is primarily comprised of the shock loss when the shock strength is not sufficient to separate the suction surface boundary layer. However, when the shock strength is sufficient to separate the suction surface boundary layer, the profile loss is comparable to the shock loss and can exceed the shock loss.

Description: The Chimera overset grid method is reviewed and discussed in the context of a method of solution and analysis of unsteady three-dimensional viscous flows. The state of maturity of the various pieces of support software required to use the approach is discussed. A variety of recent applications of the method is presented. Current limitations of the approach are defined.

Description: This final report includes five publications that resulted from the studies of the global flowfield about the V-22 Tiltrotor Aircraft. The first of the five is 'The Chimera Method of Simulation for Unsteady Three-Dimensional Viscous Flow', as presented in 'Computational Fluid Dynamics Review 1995.' The remaining papers, all presented at AIAA conferences, are 'Unsteady Simulation of the Viscous Flow About a V-22 Rotor and Wing in Hover', 'An Efficient Means of Adaptive Refinement Within Systems of Overset Grids', 'On the Spatial and Temporal Accuracy of Overset Grid Methods for MOving Body Problems', and 'Moving Body Overset Grid Methods for Complete Aircraft Tiltrotor Simulations.'

Description: A set of compressible flow relations describing flow properties across oblique shock waves, derived for a thermally perfect, calorically imperfect gas, is applied within the existing thermally perfect gas (TPG) computer code. The relations are based upon a value of cp expressed as a polynomial function of temperature. The updated code produces tables of compressible flow properties of oblique shock waves, as well as the original properties of normal shock waves and basic isentropic flow, in a format similar to the tables for normal shock waves found in NACA Rep. 1135. The code results are validated in both the calorically perfect and the calorically imperfect, thermally perfect temperature regimes through comparisons with the theoretical methods of NACA Rep. 1135, and with a state-of-the-art computational fluid dynamics code. The advantages of the TPG code for oblique shock wave calculations, as well as for the properties of isentropic flow and normal shock waves, are its ease of use, and its applicability to any type of gas (monatomic, diatomic, triatomic, polyatomic, or any specified mixture thereof).

Description: A modification of airfoil section geometry is examined for improvement of the leading edge pressures predicted by the Computational Aeroelasticity Program - Transonic Small Disturbance (CAP-TSD). Results are compared with Eppler solutions to assess improvement. Preliminary results indicate that a fading function modification of section slopes is capable of significant improvements in the pressures near the leading edge computed by CAP-TSD. Application of this modification to airfoil geometry before use in CAP-TSD is shown to reduce the nonphysical pressure peak predicted by the transonic small disturbance solver. A second advantage of the slope modification is the substantial reduction in sensitivity of CAP-TSD steady pressure solutions to the computational mesh.

Description: High quality airloads data have been obtained on an instrumented UH-60A in flight and these data provide insight into the aerodynamic limiting behavior of the rotor. At moderate weight coefficients and high advance ratio limiting performance is largely caused by high drag near the blade tip on the advancing side of the rotor as supercritical flow develops on the rotor with moderate to strong, shocks on both surfaces of the blade. Drag divergence data from two-dimensional airfoil tests show good agreement with the development of the supercritical flow regions. Large aerodynamic pitching moments are observed at high advance ratio, as well, and these pitching moments are the source of high torsional moments on the blade and control system loads. These loads occur on the advancing side of the disk and are not related to blade stall which does not occur for these weight coefficients. At high weight coefficients aerodynamic and structural limits are related to dynamic stall cycles that begin on the retreating side of the blade and, for the most severe conditions, carry around to the advancing side of the blade at the presumed first frequency of the blade/control system.

Description: This report presents a unified method for subsonic and supersonic jet analysis using the three-dimensional Navier-Stokes code PAB3D. The Navier-Stokes code was used to obtain solutions for axisymmetric jets with on-design operating conditions at Mach numbers ranging from 0.6 to 3.0, supersonic jets containing weak shocks and Mach disks, and supersonic jets with nonaxisymmetric nozzle exit geometries. This report discusses computational methods, code implementation, computed results, and comparisons with available experimental data. Very good agreement is shown between the numerical solutions and available experimental data over a wide range of operating conditions. The Navier-Stokes method using the standard Jones-Launder two-equation kappa-epsilon turbulence model can accurately predict jet flow, and such predictions are made without any modification to the published constants for the turbulence model.

Description: An experimental wind tunnel test of a 65 deg. delta wing model with interchangeable leading edges was conducted in the Langley National Transonic Facility (NTF). The objective was to investigate the effects of Reynolds and Mach numbers on slender-wing leading-edge vortex flows with four values of wing leading-edge bluntness. Experimentally obtained pressure data are presented without analysis in tabulated and graphical formats across a Reynolds number range of 6 x 10(exp 6) to 84 x 10(exp 6) at a Mach number of 0.85 and across a Mach number range of 0.4 to 0.9 at Reynolds numbers of 6 x 10(exp 6) and 60 x 10(exp 6). Normal-force and pitching-moment coefficient plots for these Reynolds number and Mach number ranges are also presented.

Description: In this paper, we propose a new technique for the numerical treatment of external flow problems with oscillatory behavior of the solution in time. Specifically, we consider the case of unbounded compressible viscous plane flow past a finite body (airfoil). Oscillations of the flow in time may be caused by the time-periodic injection of fluid into the boundary layer, which in accordance with experimental data, may essentially increase the performance of the airfoil. To conduct the actual computations, we have to somehow restrict the original unbounded domain, that is, to introduce an artificial (external) boundary and to further consider only a finite computational domain. Consequently, we will need to formulate some artificial boundary conditions (ABC's) at the introduced external boundary. The ABC's we are aiming to obtain must meet a fundamental requirement. One should be able to uniquely complement the solution calculated inside the finite computational domain to its infinite exterior so that the original problem is solved within the desired accuracy. Our construction of such ABC's for oscillating flows is based on an essential assumption: the Navier-Stokes equations can be linearized in the far field against the free-stream back- ground. To actually compute the ABC's, we represent the far-field solution as a Fourier series in time and then apply the Difference Potentials Method (DPM) of V. S. Ryaben'kii. This paper contains a general theoretical description of the algorithm for setting the DPM-based ABC's for time-periodic external flows. Based on our experience in implementing analogous ABC's for steady-state problems (a simpler case), we expect that these boundary conditions will become an effective tool for constructing robust numerical methods to calculate oscillatory flows.

Description: In the numerical computations of flows around complex configurations, accurate calculations of force and moment coefficients for aerodynamic surfaces are required. When overset grid methods are used, the surfaces on which force and moment coefficients are sought typically consist of a collection of overlapping surface grids. Direct integration of flow quantities on the overlapping grids would result in the overlapped regions being counted more than once. The FOMOCO Utilities is a software package for computing flow coefficients (force, moment, and mass flow rate) on a collection of overset surfaces with accurate accounting of the overlapped zones. FOMOCO Utilities can be used in stand-alone mode or in conjunction with the Chimera overset grid compressible Navier-Stokes flow solver OVERFLOW. The software package consists of two modules corresponding to a two-step procedure: (1) hybrid surface grid generation (MIXSUR module), and (2) flow quantities integration (OVERINT module). Instructions on how to use this software package are described in this user's manual. Equations used in the flow coefficients calculation are given in Appendix A.

Description: A variety of empirical and computational fluid dynamics two-dimensional (2-D) dynamic stall models were compared to recently obtained three-dimensional (3-D) dynamic stall data in a workshop on modeling of 3-D dynamic stall of an unswept, rectangular wing, of aspect ratio 10. Dynamic stall test data both below and above the static stall angle-of-attack were supplied to the participants, along with a 'blind' case where only the test conditions were supplied in advance, with results being compared to experimental data at the workshop itself. Detailed graphical comparisons are presented in the report, which also includes discussion of the methods and the results. The primary conclusion of the workshop was that the 3-D effects of dynamic stall on the oscillating wing studied in the workshop can be reasonably reproduced by existing semi-empirical models once 2-D dynamic stall data have been obtained. The participants also emphasized the need for improved quantification of 2-D dynamic stall.

Description: The use of canards in advanced aircraft for control and improved aerodynamic performance is a topic of continued interest and research. In addition to providing maneuver control and trim, the influence of canards on wing aerodynamics can often result in increased maximum lift and decreased trim drag. In many canard-configured aircraft, the main benefits of canards are realized during maneuver or other dynamic conditions. Therefore, the detailed study and understanding of canards requires the accurate prediction of the non-linear unsteady aerodynamics of such configurations. For close-coupled canards, the unsteady aerodynamic performance associated with the canard-wing interaction is of particular interest. The presence of a canard in close proximity to the wing results in a highly coupled canard-wing aerodynamic flowfield which can include downwash/upwash effects, vortex-vortex interactions and vortex-surface interactions. For unsteady conditions, these complexities of the canard-wing flowfield are further increased. The development and integration of advanced computational technologies provide for the time-accurate Navier-Stokes simulations of the steady and unsteady canard-wing-body flox,fields. Simulation, are performed for non-linear flight regimes at transonic Mach numbers and for a wide range of angles of attack. For the static configurations, the effects of canard positioning and fixed deflection angles on aerodynamic performance and canard-wing vortex interaction are considered. For non-static configurations, the analyses of the canard-wing body flowfield includes the unsteady aerodynamics associated with pitch-up ramp and pitch oscillatory motions of the entire geometry. The unsteady flowfield associated with moving canards which are typically used as primary control surfaces are considered as well. The steady and unsteady effects of the canard on surface pressure integrated forces and moments, and canard-wing vortex interaction are presented in detail including the effects of the canard on the static and dynamic stability characteristics. The current study provides an understanding of the steady and unsteady canard-wing-body flowfield. Emphasis is placed on the effects of the canard on aerodynamic performance as well as the detailed flow physics of the canard-wing flowfield interactions. The computational tools developed to accurately predict the time-accurate flowfield of moving canards provides for the capability of coupled fluids-controls simulations desired in the detailed design and analysis of advanced aircraft.

Description: An experimental and computational investigation of the effect of lift enhancing tabs on a two-element airfoil was conducted. The objective of the study was to develop an understanding of the flow physics associated with lift enhancing tabs on a multi-element airfoil. A NACA 63(sub 2)-215 ModB airfoil with a 30 percent chord Fowler flap was tested in the NASA Ames 7 by 10 foot wind tunnel. Lift enhancing tabs of various heights were tested on both the main element and the flap for a variety of flap riggings. Computations of the flow over the two-element airfoil were performed using the two-dimensional incompressible Navier-Stokes code INS2D-UP. The computer results predict all of the trends in the experimental data quite well. When the flow over the flap upper surface is attached, tabs mounted at the main element trailing edge (cove tabs) produce very little change in lift. At high flap deflections. however, the flow over the flap is separated and cove tabs produce large increases in lift and corresponding reductions in drag by eliminating the separated flow. Cove tabs permit high flap deflection angles to be achieved and reduce the sensitivity of the airfoil lift to the size of the flap gap. Tabs attached to the flap training edge (flap tabs) are effective at increasing lift without significantly increasing drag. A combination of a cove tab and a flap tab increased the airfoil lift coefficient by 11 percent relative to the highest lift tab coefficient achieved by any baseline configuration at an angle of attack of zero percent and the maximum lift coefficient was increased by more than 3 percent. A simple analytic model based on potential flow was developed to provide a more detailed understanding of how lift enhancing tabs work. The tabs were modeled by a point vortex at the training edge. Sensitivity relationships were derived which provide a mathematical basis for explaining the effects of lift enhancing tabs on a multi-element airfoil. Results of the modeling effort indicate that the dominant effects of the tabs on the pressure distribution of each element of the airfoil can be captured with a potential flow model for cases with no flow separation.

Description: While numerically solving a problem initially formulated on an unbounded domain, one typically truncates this domain, which necessitates setting the artificial boundary conditions (ABC's) at the newly formed external boundary. The issue of setting the ABC's appears to be most significant in many areas of scientific computing, for example, in problems originating from acoustics, electrodynamics, solid mechanics, and fluid dynamics. In particular, in computational fluid dynamics (where external problems present a wide class of practically important formulations) the proper treatment of external boundaries may have a profound impact on the overall quality and performance of numerical algorithms. Most of the currently used techniques for setting the ABC's can basically be classified into two groups. The methods from the first group (global ABC's) usually provide high accuracy and robustness of the numerical procedure but often appear to be fairly cumbersome and (computationally) expensive. The methods from the second group (local ABC's) are, as a rule, algorithmically simple, numerically cheap, and geometrically universal; however, they usually lack accuracy of computations. In this paper we first present a survey and provide a comparative assessment of different existing methods for constructing the ABC's. Then, we describe a relatively new ABC's technique of ours and review the corresponding results. This new technique, in our opinion, is currently one of the most promising in the field. It enables one to construct such ABC's that combine the advantages relevant to the two aforementioned classes of existing methods. Our approach is based on application of the difference potentials method attributable to V. S. Ryaben'kii. This approach allows us to obtain highly accurate ABC's in the form of certain (nonlocal) boundary operator equations. The operators involved are analogous to the pseudodifferential boundary projections first introduced by A. P. Calderon and then also studied by R. T. Seeley. The apparatus of the boundary pseudodifferential equations, which has formerly been used mostly in the qualitative theory of integral equations and PDE'S, is now effectively employed for developing numerical methods in the different fields of scientific computing.

Description: Aerodynamic shape design has long persisted as a difficult scientific challenge due its highly nonlinear flow physics and daunting geometric complexity. However, with the emergence of Computational Fluid Dynamics (CFD) it has become possible to make accurate predictions of flows which are not dominated by viscous effects. It is thus worthwhile to explore the extension of CFD methods for flow analysis to the treatment of aerodynamic shape design. Two new aerodynamic shape design methods are developed which combine existing CFD technology, optimal control theory, and numerical optimization techniques. Flow analysis methods for the potential flow equation and the Euler equations form the basis of the two respective design methods. In each case, optimal control theory is used to derive the adjoint differential equations, the solution of which provides the necessary gradient information to a numerical optimization method much more efficiently then by conventional finite differencing. Each technique uses a quasi-Newton numerical optimization algorithm to drive an aerodynamic objective function toward a minimum. An analytic grid perturbation method is developed to modify body fitted meshes to accommodate shape changes during the design process. Both Hicks-Henne perturbation functions and B-spline control points are explored as suitable design variables. The new methods prove to be computationally efficient and robust, and can be used for practical airfoil design including geometric and aerodynamic constraints. Objective functions are chosen to allow both inverse design to a target pressure distribution and wave drag minimization. Several design cases are presented for each method illustrating its practicality and efficiency. These include non-lifting and lifting airfoils operating at both subsonic and transonic conditions.

Description: A 1/8-scale model of a fan-in-wing concept considered for development by Grumman Aerospace Corporation for the U.S. Army was tested in the Langley 14- by 22-Foot Subsonic Tunnel. Hover testing, which included height above a pressure-instrumented ground plane, angle of pitch, and angle of roll for a range of fan thrust, was conducted in a model preparation area near the tunnel. The air loads and surface pressures on the model were measured for several configurations in the model preparation area and in the tunnel. The major hover configuration change was varying the angles of the vanes attached to the exit of the fans for producing propulsive force. As the model height above the ground was decreased, there was a significant variation of thrust-removed normal force with constant fan speed. The greatest variation was generally for the height-to-fan exit diameter ratio of less than 2.5; the variation was reduced by deflecting fan exit flow outboard with the vanes. In the tunnel angles of pitch and sideslip, height above the tunnel floor, and wind speed were varied for a range of fan thrust and different vane angle configurations. Other configuration features such as flap deflections and tail incidence were evaluated as well. Though the V-tail empennage provided an increase in static longitudinal stability, the total model configuration remained unstable.

Description: A Mathematical theory is developed to perform the calculations necessary to determine the wave drag for slender bodies of non-circular cross section. The derivations presented in this report are based on extensions to supersonic linearized small perturbation theory. A numerical scheme is presented utilizing Fourier decomposition to compute the pressure coefficient on and about a slender body of arbitrary cross section.

Description: A full-scale helicopter rotor test was conducted in the NASA Ames 80- by 120-Foot Wind Tunnel with a four-bladed S-76 rotor system. Rotor performance and loads data were obtained over a wide range of rotor shaft angles-of-attack and thrust conditions at tunnel speeds ranging from 0 to 100 kt. The primary objectives of this test were (1) to acquire forward flight rotor performance and loads data for comparison with analytical results; (2) to acquire S-76 forward flight rotor performance data in the 80- by 120-Foot Wind Tunnel to compare with existing full-scale 40- by 80-Foot Wind Tunnel test data that were acquired in 1977; (3) to evaluate the acoustic capability of the 80- by 120- Foot Wind Tunnel for acquiring blade vortex interaction (BVI) noise in the low speed range and compare BVI noise with in-flight test data; and (4) to evaluate the capability of the 80- by 120-Foot Wind Tunnel test section as a hover facility. The secondary objectives were (1) to evaluate rotor inflow and wake effects (variations in tunnel speed, shaft angle, and thrust condition) on wind tunnel test section wall and floor pressures; (2) to establish the criteria for the definition of flow breakdown (condition where wall corrections are no longer valid) for this size rotor and wind tunnel cross-sectional area; and (3) to evaluate the wide-field shadowgraph technique for visualizing full-scale rotor wakes. This data base of rotor performance and loads can be used for analytical and experimental comparison studies for full-scale, four-bladed, fully articulated rotor systems. Rotor performance and structural loads data are presented in this report.

Description: Prediction of the loss of wing leading-edge thrust and the accompanying increase in drag due to lift, when flow is not completely attached, presents a difficult but commonly encountered problem. A method (called the previous method) for the prediction of attainable leading-edge thrust and the resultant effect on airplane aerodynamic performance has been in use for more than a decade. Recently, the method has been revised to enhance its applicability to current airplane design and evaluation problems. The improved method (called the present method) provides for a greater range of airfoil shapes from very sharp to very blunt leading edges. It is also based on a wider range of Reynolds numbers than was available for the previous method. The present method, when employed in computer codes for aerodynamic analysis, generally results in improved correlation with experimental wing-body axial-force data and provides reasonable estimates of the measured drag.

Description: The numerical study of aeroacoustic problems places stringent demands on the choice of a computational algorithm, because it requires the ability to propagate disturbances of small amplitude and short wavelength. The demands are particularly high when shock waves are involved, because the chosen algorithm must also resolve discontinuities in the solution. The extent to which a high-order-accurate shock-capturing method can be relied upon for aeroacoustics applications that involve the interaction of shocks with other waves has not been previously quantified. Such a study is initiated in this work. A fourth-order-accurate essentially nonoscillatory (ENO) method is used to investigate the solutions of inviscid, compressible flows with shocks in a quasi-one-dimensional nozzle flow. The design order of accuracy is achieved in the smooth regions of a steady-state test case. However, in an unsteady test case, only first-order results are obtained downstream of a sound-shock interaction. The difficulty in obtaining a globally high-order-accurate solution in such a case with a shock-capturing method is demonstrated through the study of a simplified, linear model problem. Some of the difficult issues and ramifications for aeroacoustics simulations of flows with shocks that are raised by these results are discussed.

Description: On a swept wing, contamination along the leading edge, Tollmien-Schlichting waves, stationary or traveling crossflow vortices, and/or Taylor-Gortler vortices can cause the catastrophic breakdown of laminar to turbulent flow, which leads to increased skin-friction drag for the aircraft. The discussion in this Note will be limited to disturbances which evolve along the attachment line (leading edge of swept wing). If the Reynolds number of the attachment-line boundary layer is greater than some critical value, then the complete wing is inevitably engulfed in turbulent flow. Essentially, there are two critical Reynolds number points that must be considered. The first is for small-amplitude disturbances, and the second is for bypass transition. The present study will use direct numerical simulations to validate a linear 2D-eigenvalue prediction method based on parabolized stability equations by Lin and Malik. This method is considered because it suggests that a number of symmetric and asymmetric modes exist and are stable or unstable on the attachment line depending on the Reynolds number. If validated, the approach would predict a number of modes which are linearly damped in the Reynolds number regime 100 to 245; however, these modes may grow nonlinearly and provide an explanation to this region.

Description: The linear and the nonlinear stability of disturbances that propagate along the attachment line of a three-dimensional boundary layer is considered. The spatially evolving disturbances in the boundary layer are computed by direct numerical simulation (DNS) of the unsteady, incompressible Navier-Stokes equations. Disturbances are introduced either by forcing at the in ow or by applying suction and blowing at the wall. Quasi-parallel linear stability theory and a nonparallel theory yield notably different stability characteristics for disturbances near the critical Reynolds number; the DNS results con rm the latter theory. Previously, a weakly nonlinear theory and computations revealed a high wave-number region of subcritical disturbance growth. More recent computations have failed to achieve this subcritical growth. The present computational results indicate the presence of subcritically growing disturbances; the results support the weakly nonlinear theory. Furthermore, an explanation is provided for the previous theoretical and computational discrepancy. In addition, the present results demonstrate that steady suction can be used to stabilize disturbances that otherwise grow subcritically along the attachment line.

Description: The sonic boom propagation codes reviewed in this study, SHOCKN and ZEPHYRUS, implement current theory on air absorption using different computational concepts. Review of the codes with a realistic atmosphere model confirm the agreement of propagation results reported by others for idealized propagation conditions. ZEPHYRUS offers greater flexibility in propagation conditions and is thus preferred for practical aircraft analysis. The ZEPHYRUS code was used to propagate sonic boom waveforms measured approximately 1000 feet away from an SR-71 aircraft flying at Mach 1.25 to 5000 feet away. These extrapolated signatures were compared to measurements at 5000 feet. Pressure values of the significant shocks (bow, canopy, inlet and tail) in the waveforms are consistent between extrapolation and measurement. Of particular interest is that four (independent) measurements taken under the aircraft centerline converge to the same extrapolated result despite differences in measurement conditions. Agreement between extrapolated and measured signature duration is prevented by measured duration of the 5000 foot signatures either much longer or shorter than would be expected. The duration anomalies may be due to signature probing not sufficiently parallel to the aircraft flight direction.

Description: The effect of a surface wave on transition in three-dimensional boundary-layer flow over an infinite swept wing was studied. The mean flow computed using interacting boundary-layer theory, and transition was predicted using linear stability theory coupled with the empirical eN method. It was found that decreasing the wave height, sweep angle, or freestream unit Reynolds number, and increasing the freestream Mach number or suction level all stabilized the flow and moved transition onset to downstream locations.

Description: Pressure distributions on three NACA 1-series inlets have been obtained in the Langley 16-Foot Transonic Tunnel. The cowl diameter ratio (ratio of cowl highlight diameter to cowl maximum diameter) was 0.85 for all three inlets. The cowl length ratio (ratio of cowl length to cowl maximum diameter) was 1.0 for two of the inlets (NACA 1-85-100) and 0.439 for the other (NACA 1-85-43.9) inlet. One of the inlets with a cowl length ratio of 1.0 had an internal contraction ratio (ratio of highlight area to throat area) of 1.009 and the other had a contraction ratio of 1.250. The inlet with a cowl length ratio of 0.439 also had an internal contraction ratio of 1.250. All three inlets had longitudinal rows of static pressure orifices on the top and bottom external cowl surfaces. The inlet with a contraction ratio of 1.009 also had a row of static pressure orifices on the side of the cowl (external surface). The two inlets with a contraction ratio of 1.250 had a longitudinal row of static pressure orifices on the diffuser surface.

Description: The initial structure and axial decay of an array of streamwise vortices embedded in a turbulent pipe boundary layer is experimentally investigated. The vortices are shed in counter-rotating fashion from an array of equally-spaced symmetric airfoil vortex generators. Vortex structure is quantified in terms of crossplane circulation and peak streamwise vorticity. Flow conditions are subsonic and incompressible. The focus of this study is on the effect of the initial spacing between the parent vortex generators. Arrays with vortex generators spaced at 15 and 30 degrees apart are considered. When the spacing between vortex generators is decreased the circulation and peak vorticity of the shed vortices increases. Analysis indicates this strengthening results from regions of fluid acceleration in the vicinity of the vortex generator array. Decreased spacing between the constituent vortices also produces increased rates of circulation and peak vorticity decay.

Description: Computer programs to produce the ordinates for airfoils of any thickness, thickness distribution, or camber in the NACA airfoil series were developed in the early 1970's and are published as NASA TM X-3069 and TM X-3284. For analytic airfoils, the ordinates are exact. For the 6-series and all but the leading edge of the 6A-series airfoils, agreement between the ordinates obtained from the program and previously published ordinates is generally within 5 x 10(exp -5) chord. Since the publication of these programs, the use of personal computers and individual workstations has proliferated. This report describes a computer program that combines the capabilities of the previously published versions. This program is written in ANSI FORTRAN 77 and can be compiled to run on DOS, UNIX, and VMS based personal computers and workstations as well as mainframes. An effort was made to make all inputs to the program as simple as possible to use and to lead the user through the process by means of a menu.

Description: The use of gradient based optimization algorithms in inverse design is well established as a practical approach to aerodynamic design. A typical procedure uses a simulation scheme to evaluate the objective function (from the approximate states) and its gradient, then passes this information to an optimization algorithm. Once the simulation scheme (CFD flow solver) has been selected and used to provide approximate function evaluations, there are several possible approaches to the problem of computing gradients. One popular method is to differentiate the simulation scheme and compute design sensitivities that are then used to obtain gradients. Although this black-box approach has many advantages in shape optimization problems, one must compute mesh sensitivities in order to compute the design sensitivity. In this paper, we present an alternative approach using the PDE sensitivity equation to develop algorithms for computing gradients. This approach has the advantage that mesh sensitivities need not be computed. Moreover, when it is possible to use the CFD scheme for both the forward problem and the sensitivity equation, then there are computational advantages. An apparent disadvantage of this approach is that it does not always produce consistent derivatives. However, for a proper combination of discretization schemes, one can show asymptotic consistency under mesh refinement, which is often sufficient to guarantee convergence of the optimal design algorithm. In particular, we show that when asymptotically consistent schemes are combined with a trust-region optimization algorithm, the resulting optimal design method converges. We denote this approach as the sensitivity equation method. The sensitivity equation method is presented, convergence results are given and the approach is illustrated on two optimal design problems involving shocks.

Description: This paper presents the results of a computational flow analysis of the McDonnell Douglas single-stage-to-orbit vehicle concept designated as the 24U. This study was made to determine the aerodynamic characteristics of the vehicle with and without body flaps over an angle of attack range of 20-40 deg. Computations were made at a flight Mach number of 20 at 200,000 ft. altitude with equilibrium air, and a Mach number of 6 with CF4 gas. The software package FELISA (Finite Element Langley imperial College Sawansea Ames) was used for all the computations. The FELISA software consists of unstructured surface and volume grid generators, and inviscid flow solvers with (1) perfect gas option for subsonic, transonic, and low supersonic speeds, and (2) perfect gas, equilibrium air, and CF4 options for hypersonic speeds. The hypersonic flow solvers with equilibrium air and CF4 options were used in the present studies. Results are compared with other computational results and hypersonic CF4 tunnel test data.

Description: An evaluation was made on the effects of integrating the required aircraft components with hypersonic high-lift configurations known as waveriders to create hypersonic cruise vehicles. Previous studies suggest that waveriders offer advantages in aerodynamic performance and propulsion/airframe integration (PAI) characteristics over conventional non-waverider hypersonic shapes. A wind-tunnel model was developed that integrates vehicle components, including canopies, engine components, and control surfaces, with two pure waverider shapes, both conical-flow-derived waveriders for a design Mach number of 4.0. Experimental data and limited computational fluid dynamics (CFD) solutions were obtained over a Mach number range of 1.6 to 4.63. The experimental data show the component build-up effects and the aerodynamic characteristics of the fully integrated configurations, including control surface effectiveness. The aerodynamic performance of the fully integrated configurations is not comparable to that of the pure waverider shapes, but is comparable to previously tested hypersonic models. Both configurations exhibit good lateral-directional stability characteristics.

Description: Hypersonic boundary layer measurements over a flared cone were conducted in a Mach 6 quiet wind tunnel at a freestream unit Reynolds number of 2.82 million/ft. This Reynolds number provided laminar-to-transitional flow over the cone model in a low-disturbance environment. Four interchangeable nose-tips, including a sharp-tip, were tested. Point measurements with a single hot-wire using a novel constant voltage anemometer were used to measure the boundary layer disturbances. Surface temperature and schlieren measurements were also conducted to characterize the transitional state of the boundary layer and to identify instability modes. Results suggest that second mode disturbances were the most unstable and scaled with the boundary layer thickness. The second mode integrated growth rates compared well with linear stability theory in the linear stability regime. The second mode is responsible for transition onset despite the existence of a second mode subharmonic. The subharmonic disturbance wavelength also scales with the boundary layer thickness. Furthermore, the existence of higher harmonics of the fundamental suggests that nonlinear disturbances are not associated with 'high' free stream disturbance levels. Nose-tip radii greater than 2.7% of the base radius completely stabilized the second mode.

Description: Flight experiments were conducted on Ames Research Center's V/STOL Systems Research Aircraft (VSRA) to assess the influence of advanced control modes and head-up displays (HUD's) on flying qualities for precision approach and landing operations. Evaluations were made for decelerating approaches to hover followed by a vertical landing and for slow landings for four control/display mode combinations: the basic YAV-8B stability augmentation system; attitude command for pitch, roll, and yaw; flightpath/acceleration command with translational rate command in the hover; and height-rate damping with translational-rate command. Head-up displays used in conjunction with these control modes provided flightpath tracking/pursuit guidance and deceleration commands for the decelerating approach and a mixed horizontal and vertical presentation for precision hover and landing. Flying qualities were established and control usage and bandwidth were documented for candidate control modes and displays for the approach and vertical landing. Minimally satisfactory bandwidths were determined for the translational-rate command system. Test pilot and engineer teams from the Naval Air Warfare Center, the Boeing Military Airplane Group, Lockheed Martin, McDonnell Douglas Aerospace, Northrop Grumman, Rolls-Royce, and the British Defense Research Agency participated in the program along with NASA research pilots from the Ames and Lewis Research Centers. The results, in conjunction with related ground-based simulation data, indicate that the flightpath/longitudinal acceleration command response type in conjunction with pursuit tracking and deceleration guidance on the HUD would be essential for operation to instrument minimums significantly lower than the minimums for the AV-8B. It would also be a superior mode for performing slow landings where precise control to an austere landing area such as a narrow road is demanded. The translational-rate command system would reduce pilot workload for demanding vertical landing tasks aboard ship and in confined land-based sites.

Description: This paper presents aerodynamic characteristics and pressure distributions for an executive-jet modified airfoil and discusses drag reduction relative to a baseline airfoil for two cruise design points. A modified airfoil was tested in the adaptive-wall test section of the NASA Langley 0.3-Meter Transonic Cryogenic Tunnel (0.3-m TCT) for Mach numbers ranging from 0.250 to 0.780 and chord Reynolds numbers ranging from 3.0 x 10(exp 6) to 18.0 x 10(exp 6). The angle of attack was varied from minus 2 degrees to almost 10 degrees. Boundary-layer transition was fixed at 5 percent of chord on both the upper and lower surfaces of the model for most of the test. The two design Mach numbers were 0.654 and 0.735, chord Reynolds numbers were 4.5 x 10(exp 6) and 8.9 x 10(exp 6), and normal-force coefficients were 0.98 and 0.51. Test data are presented graphically as integrated force and moment coefficients and chordwise pressure distributions. The maximum normal-force coefficient decreases with increasing Mach number. At a constant normal-force coefficient in the linear region, as Mach number increases an increase occurs in the slope of normal-force coefficient versus angle of attack, negative pitching-moment coefficient, and drag coefficient. With increasing Reynolds number at a constant normal-force coefficient, the pitching-moment coefficient becomes more negative and the drag coefficient decreases. The pressure distributions reveal that when present, separation begins at the trailing edge as angle of attack is increased. The modified airfoil, which is designed with pitching moment and geometric constraints relative to the baseline airfoil, achieved drag reductions for both design points (12 and 22 counts). The drag reductions are associated with stronger suction pressures in the first 10 percent of the upper surface and weakened shock waves.

Description: This paper compares and evaluates numerical and experimental flowfields of the RAH-66 Comanche helicopter. The numerical predictions were obtained by solving the Thin-Layer Navier-Stokes equations. The computations use actuator disks to investigate the main and tail rotor effects upon the fuselage flowfield. The wind tunnel experiment was performed in the 14 x 22 foot facility located at NASA Langley. A suite of flow conditions, rotor thrusts and fuselage-rotor-tail configurations were tested. In addition, the tunnel model and the computational geometry were based upon the same CAD definition. Computations were performed for an isolated fuselage configuration and for a rotor on configuration. Comparisons between the measured and computed surface pressures show areas of correlation and some discrepancies. Local areas of poor computational grid-quality and local areas of geometry differences account for the differences. These calculations demonstrate the use of advanced computational fluid dynamic methodologies towards a flight vehicle currently under development. It serves as an important verification for future computed results.

Description: As a rotor blade moves through the air, it sheds vortices. These vortices shed along the length of the blade over time form the wake. The strongest vortices of the wake are those trailing from the tip of the blade. When a rotating blade system moves under certain operating conditions, each blade will impinge on the tip vortices shed by itself or other blades. This impingement is called a blade-vortex interaction, or BVI. Although the blade and trailing tip vortices interact with many different orientations, one of the two extremes, either parallel or perpendicular interaction, is usually modelled. In a perpendicular interaction, the portion of the blade that is actually interacting with the travelling vortex at any given time is very small. A parallel interaction, however, has the largest concurrent interaction with the blade, as a result this case is given the most attention. One of the most commonly studied occurrences of blade-vortex interactions is associated with low-speed descending rotorcraft flight. BVI occur when the tip vortices shed by the blades intersect the plane of the rotor. BVI cause local pressure changes over the blades which are responsible, in part, for the acoustic signature of the rotorcraft. The local pressure changes also cause vibrations which lead to fatigue of both the blades and the mechanical components driving the blades.

Description: An experimental wind tunnel test of a 65 deg delta wing model with interchangeable leading edges was conducted in the Langley National Transonic Facility (NTF). The objective was to investigate the effects of Reynolds and Mach numbers on slender-wing leading-edge vortex flows with four values of wing leading-edge bluntness. Experimentally obtained pressure data are presented without analysis in tabulated and graphical formats across a Reynolds number range of 6 x 10(exp 6) to 120 x 10(exp 6) at a Mach number of 0.85 and across a Mach number range of 0.4 to 0.9 at Reynolds numbers of 6 x 10(exp 6) and 60 x 10(exp 6). Normal-force and pitching-moment coefficient plots for these Reynolds number and Mach number ranges are also presented.

Description: A test was conducted on a model of the NACA 0012 airfoil section with a solid upper surface or a porous upper surface with a cavity beneath for passive venting. The purposes of the test were to investigate the aerodynamic characteristics of an airfoil with full-chord porosity and to assess the ability of porosity to provide a multipoint or self-adaptive design. The tests were conducted in the Langley 8-Foot Transonic Pressure Tunnel over a Mach number range from 0.50 to 0.82 at chord Reynolds numbers of 2 x 10(exp 6), 4 x 10(exp 6), and 6 x 10(exp 6). The angle of attack was varied from -1 deg to 6 deg. At the lower Mach numbers, porosity leads to a dependence of the drag on the normal force. At subcritical conditions, porosity tends to flatten the pressure distribution, which reduces the suction peak near the leading edge and increases the suction over the middle of the chord. At supercritical conditions, the compression region on the porous upper surface is spread over a longer portion of the chord. In all cases, the pressure coefficient in the cavity beneath the porous surface is fairly constant with a very small increase over the rear portion. For the porous upper surface, the trailing edge pressure coefficients exhibit a creep at the lower section normal force coefficients, which suggests that the boundary layer on the rear portion of the airfoil is significantly thickening with increasing normal force coefficient.

Description: An experimental wind tunnel test of a 65 deg delta wing model with interchangeable leading edges was conducted in the Langley National Transonic Facility (NTF). The objective was to investigate the effects of Reynolds and Mach numbers on slender-wing leading-edge vortex flows with four values of wing leading-edge bluntness. Experimentally obtained pressure data are presented without analysis in tabulated and graphical formats across a Reynolds number range of 6 x 10(exp 6) to 120 x 10(exp 6) at a Mach number of 0.85 and across a Mach number range of 0.4 to 0.9 at Reynolds numbers of 6 x 10(exp 6), 60 x 10(exp 6), and 120 x 10(exp 6). Normal-force and pitching-moment coefficient plots for these Reynolds number and Mach number ranges are also presented.

Description: The purpose of this report is to document the results of experiments performed at the University of Kansas and at the NASA Langley Research Center (LaRC) into the use of shaped cavities to generate vortices in supersonic flow, as well as the progress made in simulating the observed flow using the PAB3D flow solver. The investigation was performed on 18 different cavity configurations installed in a convergent-divergent nozzle at the Jet Exit Facility at the LaRC. Pressure sensitive paint, static-pressure ports, focusing Schliern, and water tunnel flow visualization techniques were used to study the nature of the flow created by these cavities. The results of these investigations revealed that a shaped cavity can generate a pair of counter-rotating streamwise vortices in supersonic flow by creating weak, compression Mach waves and weak shocks. The PAB3D computer program, developed at the LaRC, was used to attempt to reproduce the experimental results. Unfortunately, due to problems with matching the grid blocks, no converged results were obtained. However, intermediate results, as well as a complete definition of the grid matching problems and suggested courses of actions are presented.

Description: This paper develops a dynamic model for pressure sensors in continuum and rarefied flows with longitudinal temperature gradients. The model was developed from the unsteady Navier-Stokes momentum, energy, and continuity equations and was linearized using small perturbations. The energy equation was decoupled from momentum and continuity assuming a polytropic flow process. Rarefied flow conditions were accounted for using a slip flow boundary condition at the tubing wall. The equations were radially averaged and solved assuming gas properties remain constant along a small tubing element. This fundamental solution was used as a building block for arbitrary geometries where fluid properties may also vary longitudinally in the tube. The problem was solved recursively starting at the transducer and working upstream in the tube. Dynamic frequency response tests were performed for continuum flow conditions in the presence of temperature gradients. These tests validated the recursive formulation of the model. Model steady-state behavior was analyzed using the final value theorem. Tests were performed for rarefied flow conditions and compared to the model steady-state response to evaluate the regime of applicability. Model comparisons were excellent for Knudsen numbers up to 0.6. Beyond this point, molecular affects caused model analyses to become inaccurate.

Description: The Boeing-Vertol VR-7 airfoil was experimentally studied with steady and pulsed upper-surface blowing for sinusoidal pitching oscillations described by alpha = alpha(sub m) + 10 deg sin(omega t). The tests were conducted in the U.S. Army Aeroflightdynamics Directorate's Water Tunnel at NASA Ames Research Center. The experiment was performed at a Reynolds number of 100,000. Pitch oscillations with alpha(sub m) = 10 deg and 15 deg and with reduced frequencies ranging from k = 0.005 to 0.15 were examined. Blowing conditions ranged from C(sub mu) = 0.03 to 0.66 and F(+) = 0 to 3. Unsteady lift, drag, and pitching-moment loads were measured, and fluorescent-dye flow visualizations were obtained. Steady, upper-surface blowing was found to be capable of trapping a separation bubble near the leading edge during a portion of the airfoil's upward rotation. When this occurred, the lift was increased significantly and stall was averted. In all cases, steady blowing reduced the hysteresis amplitudes present in the loads and produced a large thrust force. The benefits of steady blowing diminished as the reduced frequency and mean angle of oscillation increased. Pulsed blowing showed only marginal benefits for the conditions tested. The greatest gains from pulsed blowing were achieved at F(+) = 0.9.

Description: Extensive measurements are made in a transitioning swept-wing boundary layer using hot-film, hot-wire and cross-wire anemometry. The crossflow-dominated flow contains stationary vortices that breakdown near mid-chord. The most amplified vortex wavelength is forced by the use of artificial roughness elements near the leading edge. Two-component velocity and spanwise surface shear-stress correlation measurements are made at two constant chord locations, before and after transition. Streamwise surface shear stresses are also measured through the entire transition region. Correlation techniques are used to identify stationary structures in the laminar regime and coherent structures in the turbulent regime. Basic techniques include observation of the spatial correlations and the spatially distributed auto-spectra. The primary and secondary instability mechanisms are identified in the spectra in all measured fields. The primary mechanism is seen to grow, cause transition and produce large-scale turbulence. The secondary mechanism grows through the entire transition region and produces the small-scale turbulence. Advanced techniques use Linear Stochastic Estimation (LSE) and Proper Orthogonal Decomposition (POD) to identify the spatio-temporal evolutions of structures in the boundary layer. LSE is used to estimate the instantaneous velocity fields using temporal data from just two spatial locations and the spatial correlations. Reference locations are selected using maximum RMS values to provide the best available estimates. POD is used to objectively determine modes characteristic of the measured flow based on energy. The stationary vortices are identified in the first laminar modes of each velocity component and shear component. Experimental evidence suggests that neighboring vortices interact and produce large coherent structures with spanwise periodicity at double the stationary vortex wavelength. An objective transition region detection method is developed using streamwise spatial POD solutions which isolate the growth of the primary and secondary instability mechanisms in the first and second modes, respectively. Temporal evolutions of dominant POD modes in all measured fields are calculated. These scalar POD coefficients contain the integrated characteristics of the entire field, greatly reducing the amount of data to characterize the instantaneous field. These modes may then be used to train future flow control algorithms based on neural networks.

Description: An experimental and computational investigation of the effect of lift-enhancing tabs on a two-element airfoil has been conducted. The objective of the study was to develop an understanding of the flow physics associated with lift-enhancing tabs on a multi-element airfoil. An NACA 63(2)-215 ModB airfoil with a 30% chord fowler flap was tested in the NASA Ames 7- by 10-Foot Wind Tunnel. Lift-enhancing tabs of various heights were tested on both the main element and the flap for a variety of flap riggings. A combination of tabs located at the main element and flap trailing edges increased the airfoil lift coefficient by 11% relative to the highest lift coefficient achieved by any baseline configuration at an angle of attack of 0 deg, and C(sub 1max) was increased by 3%. Computations of the flow over the two-element airfoil were performed using the two-dimensional incompressible Navier-Stokes code INS2D-UP. The computed results predicted all of the trends observed in the experimental data quite well. In addition, a simple analytic model based on potential flow was developed to provide a more detailed understanding of how lift-enhancing tabs work. The tabs were modeled by a point vortex at the air-foil or flap trailing edge. Sensitivity relationships were derived which provide a mathematical basis for explaining the effects of lift-enhancing tabs on a multi-element airfoil. Results of the modeling effort indicate that the dominant effects of the tabs on the pressure distribution of each element of the airfoil can be captured with a potential flow model for cases with no flow separation.

Description: The unsteady, three-dimensional Navier-Stokes equations coupled with the Euler equations of rigid-body dynamics are sequentially solved to simulate and analyze the aerodynamic response of a high angle of attack delta wing undergoing oscillatory motion. The governing equations of fluid flow and dynamics of the multidisciplinary problem are solved using a time-accurate solution of the laminar, unsteady, compressible, full Navier- Stokes equations with the implicit, upwind, Roe flux-difference splitting, finite-volume scheme and a four-stage Runge-Kutta scheme, respectively. The primary model under consideration consists of a 65 deg swept, sharp-edged, cropped delta wing of zero thickness at 20 deg angle of attack. In a freestream of Mach 0.85 and Reynolds number of 3.23 x 10(exp 6), the flow over the upper surface of the wing develops a complex shock system which interacts with the leading-edge primary vortices producing vortex breakdown. The effect of the oscillatory motion of the wing on the vortex breakdown and overall aerodynamic response is detailed to provide insight to the complicated physics associated with unsteady flows and the phenomenon of wing rock. Forced sinusoidal single and coupled mode rolling and pitching motion is presented for the wing in a transonic freestream. The Reynolds number, frequency of oscillation, and the phase angle are varied. Comparison between the single and coupled mode forced rolling and pitching oscillation cases illustrate the effects of coupling the motion. This investigation shows that even when coupled, forced rolling oscillation at a reduced frequency of 2(pi) eliminates the vortex breakdown which results in an increase in lift. The coupling effect for in phase forced oscillations show that the lift coefficient of the pitching-alone case and the rolling-moment coefficient of the rolling-alone case dominate the resulting response. However, with a phase lead in the pitching motion, the coupled motion results in a non-periodic response of the rolling moment. The second class of problems involve releasing the wing in roll to respond to the flowfield. Two models of sharp-edged delta wings, the previous 65 deg swept model and an 80 deg swept, sharp-edged delta wing, are used to observe the aerodynamic response of a wing free to roll in a transonic and subsonic freestream, respectively. These cases demonstrate damped oscillations, self-sustained limit cycle oscillations, and divergent rolling oscillations. Ultimately, an active control model using a mass injection system was applied on the surface of the wing to suppress the self-sustained limit cycle oscillation known as wing rock. Comparisons with experimental investigations complete this study, validating the analysis and illustrating the complex details afforded by computational investigations.

Description: NASA Dryden Flight Research Center has developed a second-generation flight test fixture for use as a generic test bed for aerodynamic and fluid mechanics research. The Flight Test Fixture 2 (FTF-2) is a low-aspect-ratio vertical fin-like shape that is mounted on the centerline of the F-I5B lower fuselage. The fixture is designed for flight research at Mach numbers to a maximum of 2.0. The FTF-2 is a composite structure with a modular configuration and removable components for functional flexibility. This report documents the flow environment of the fixture, such as surface pressure distributions and boundary-layer profiles, throughout a matrix of conditions within the F-15B/FTF-2 flight envelope. Environmental conditions within the fixture are presented to assist in the design and testing of future avionics and instrumentation. The intent of this document is to serve as a user's guide and assist in the development of future flight experiments that use the FTF-2 as a test bed. Additional information enclosed in the appendices has been included to assist with more detailed analyses, if required.

Description: An optimization approach to rocket nozzle design, based on computational fluid dynamics (CFD) methodology, is investigated for low-Reynolds-number cases. This study is undertaken to determine the benefits of this approach over those of classical design processes such as Rao's method. A CFD-based optimization procedure, using the parabolized Navier-Stokes (PNS) equations, is used to design conical and contoured axisymmetric nozzles. The advantage of this procedure is that it accounts for viscosity during the design process; other processes make an approximated boundary-layer correction after an inviscid design is created. Results showed significant improvement in the nozzle thrust coefficient over that of the baseline case; however, the unusual nozzle design necessitates further investigation of the accuracy of the PNS equations for modeling expanding flows with thick laminar boundary layers.

Description: Two-dimensional incompressible viscous driven-cavity flows are computed for Reynolds numbers on the range 100-20,000 using a loosely coupled, implicit, second-order centrally-different scheme. Mesh sequencing and three-level V-cycle multigrid error smoothing are incorporated into the symmetric Gauss-Seidel time-integration algorithm. Parametrics on the numerical parameters are performed, achieving reductions in solution times by more than 60 percent with the full multigrid approach. Details of the circulation patterns are investigated in cavities of 2-to-1, 1-to-1, and 1-to-2 depth to width ratios.

Description: The behavior of the velocity gradient tensor, A(ij)=delta u(i)/delta x(j), was studied using three turbulent flows obtained from direct numerical simulation The flows studies were: an inviscid calculation of the interaction between two vortex tubes, a homogeneous isotropic flow, and a temporally evolving planar wake. Self-similar behavior for each flow was obtained when A(ij) was normalized with the mean strain rate. The case of the interaction between two vortex tubes revealed a finite sized coherent structure with topological characteristics predictable by a restricted Euler model. This structure was found to evolve with the peak vorticity as the flow approached singularity. Invariants of A(ij) within this structure followed a straight line relationship of the form: gamma(sup 3)+gammaQ+R=0, where Q and R are the second and third invariants of A(ij), and the eigenvalue gamma is nearly constant over the volume of this structure. Data within this structure have local strain topology of unstable-node/saddle/saddle. The characteristics of the velocity gradient tensor and the anisotropic part of a related acceleration gradient tensor H(ij) were also studied for a homogeneous isotropic flow and a temporally evolving planar wake. It was found that the intermediate principal eigenvalue of the rate-of-strain tensor of H(ij) tended to be negative, with local strain topology of the type stable-node/saddle/saddle. There was also a preferential eigenvalue direction. The magnitude of H(ij) in the wake flow was found to be very small when data were conditioned at high local dissipation regions. This result was not observed in the relatively low Reynolds number simulation of homogeneous isotropic flow. A restricted Euler model of the evolution of A(ij) was found to reproduce many of the topological features identified in the simulations.

Description: A three-dimensional, linearized, Euler analysis is being developed to provide an efficient unsteady aerodynamic analysis that can be used to predict the aeroelastic and aeroacoustic response characteristics of axial-flow turbomachinery blading. The field equations and boundary conditions needed to describe nonlinear and linearized inviscid unsteady flows through a blade row operating within a cylindrical annular duct are presented. In addition, a numerical model for linearized inviscid unsteady flow, which is based upon an existing nonlinear, implicit, wave-split, finite volume analysis, is described. These aerodynamic and numerical models have been implemented into an unsteady flow code, called LINFLUX. A preliminary version of the LINFLUX code is applied herein to selected, benchmark three-dimensional, subsonic, unsteady flows, to illustrate its current capabilities and to uncover existing problems and deficiencies. The numerical results indicate that good progress has been made toward developing a reliable and useful three-dimensional prediction capability. However, some problems, associated with the implementation of an unsteady displacement field and numerical errors near solid boundaries, still exist. Also, accurate far-field conditions must be incorporated into the FINFLUX analysis, so that this analysis can be applied to unsteady flows driven be external aerodynamic excitations.

Description: We examine technology trade-offs related to grayscale resolution, spatial resolution, and error diffusion for tessellated display systems. We present new empirical results from our psychophysical study of these trade-offs and compare them to the predictions of a model of human vision.

Description: This study focuses on the prediction of airfoil characteristics, including lift and drag over a range of Reynolds numbers. Two different turbulence models, which represent two different types of models, are tested. The first is a standard isotropic eddy-viscosity two-equation model, and the second is an explicit algebraic stress model (EASM). The turbulent flow field over a general-aviation airfoil (GA(W)-2) at three Reynolds numbers is studied. At each Reynolds number, predicted lift and drag values at different angles of attack are compared with experimental results, and predicted variations of stall locations with Reynolds number are compared with experimental data. Finally, the size of the separation zone predicted by each model is analyzed, and correlated with the behavior of the lift coefficient near stall. In summary, the EASM model is able to predict the lift and drag coefficients over a wider range of angles of attack than the two-equation model for the three Reynolds numbers studied. However, both models are unable to predict the correct lift and drag behavior near the stall angle, and for the lowest Reynolds number case, the two-equation model did not predict separation on the airfoil near stall.

Description: Large-eddy simulations of a turbulent boundary layer with Reynolds number based on displacement thickness equal to 3500 were performed with two grid resolutions. The computations were continued for sufficient time to obtain frequency spectra with resolved frequencies that correspond to the most important structural frequencies on an aircraft fuselage. The turbulent stresses were adequately resolved with both resolutions. Detailed quantitative analysis of a variety of statistical quantities associated with the wall-pressure fluctuations revealed similar behavior for both simulations. The primary differences were associated with the lack of resolution of the high-frequency data in the coarse-grid calculation and the increased jitter (due to the lack of multiple realizations for averaging purposes) in the fine-grid calculation. A new curve fit was introduced to represent the spanwise coherence of the cross-spectral density.

Description: This progress report presents the results of an investigation focused on parameter identification for the NASA F/A-18 HARV. This aircraft was used in the high alpha research program at the NASA Dryden Flight Research Center. In this study the longitudinal and lateral-directional stability derivatives are estimated from flight data using the Maximum Likelihood method coupled with a Newton-Raphson minimization technique. The objective is to estimate an aerodynamic model describing the aircraft dynamics over a range of angle of attack from 5 deg to 60 deg. The mathematical model is built using the traditional static and dynamic derivative buildup. Flight data used in this analysis were from a variety of maneuvers. The longitudinal maneuvers included large amplitude multiple doublets, optimal inputs, frequency sweeps, and pilot pitch stick inputs. The lateral-directional maneuvers consisted of large amplitude multiple doublets, optimal inputs and pilot stick and rudder inputs. The parameter estimation code pEst, developed at NASA Dryden, was used in this investigation. Results of the estimation process from alpha = 5 deg to alpha = 60 deg are presented and discussed.

Description: In an effort to increase airport productivity, several wind-tunnel and flight-test programs are currently underway to determine safe reductions in separation standards between aircraft. These programs are designed to study numerous concepts from the characteristics and detection of wake vortices to the wake-vortex encounter phenomenon. As part of this latter effort, computational tools are being developed and utilized as a means of modeling and verifying wake-vortex hazard encounters. The objective of this study is to assess the ability of PMARC, a low-order potential-flow panel method, to predict the forces and moments imposed on a following business-jet configuration by a vortex interaction. Other issues addressed include the investigation of several wake models and their ability to predict wake shape and trajectory, the validity of the velocity field imposed on the following configuration, modeling techniques and the effect of the high-lift system and the empennage. Comparisons with wind-tunnel data reveal that PMARC predicts the characteristics for the clean wing-body following configuration fairly well. Non-linear effects produced by the addition of the high-lift system and empennage, however, are not so well predicted.

Description: In generating the lift on a wing, the static stall is a severe barrier. As the angle of attack, alpha, increases to the stall angle, alpha(sub stall) the flow separation point on the upper surface of the wing moves to the leading edge, so that on a two-dimensional airfoil or a large-aspect-ratio wing, the lift abruptly drops to a very low level. Therefore, the first generation of aeronautical flow type, i.e., the attached steady flow, has been limited to alpha less than alpha(sub stall). Owing to the obvious importance in applications, therefore, a great effort has been made in the past two decades to enlarge the range of usable angles of attack by various flow controls for a large-aspect-ratio wing. Basically, relevant works fall into two categories. The first category is usually refereed to as separation control, which concentrates on partially separated flow at alpha less than alpha(sub stall). Since the first experimental study of Collins and Zelenevitz, there has been ample literature showing that a partially separated flow can be turned to almost fully attached by flow controls, so that the lift is recovered and the stall is delayed (for a recent work see Seifert et al.). It has been well established that, in this category, unsteady controls are much more effective than steady ones and can be realized at a very low power-input level (Wu et al.; Seifert et al.). The second and more ambitious category of relevant efforts is the post-stall lift enhancement. Its possibility roots at the existence of a second lift peak at a very high angle of attack. In fact, As alpha further increases from alpha(sub stall), the completely separated flow develops and gradually becomes a bluff-body flow. This flow gives a normal force to the airfoil with a lift component, which reaches a peak at a maximum utilizable angle of attack, alpha(sub m) approx.= 40 deg. This second peak is of the same level as the first lift peak at alpha(sub stall). Meanwhile, the drag is also quickly increased (e.g., Fage and Johansen ; Critzos et al.). Figure 1 shows a typical experimental lift and drag coefficients of NACA-0012 airfoil in this whole range of angle of attack. Obviously, without overcoming the lift crisis at alpha(sub stall) the second lift peak is completely useless. Thus, the ultimate goal of post-stall lift enhancement is to fill the lift valley after stall by flow controls, so that a wing and/or flap can work at the whole range of 0 deg less than alpha less than alpha(sub m). Relevant early experimental studies have been extensively reviewed by Wu et al., who concluded that, first, similar to the leading-edge vortex on a slender wing, the lift enhancement on a large-aspect-ratio wing should be the result of capturing a vortex on the upper surface of the wing; and, second, using steady controls cannot reach the goal, and one must rely on unsteady controls with low-level power input as well. Wu et al. also conjectured that the underlying physics of post-stall lift enhancement by unsteady controls consists of a chain of mechanisms: vortex layer instability - receptivity resonance - nonlinear streaming.

Description: A process is presented by which aeroelastic analysis is performed by using an advanced computational fluid dynamics (CFD) code coupled with an advanced computational structural dynamics (CSD) code. The process is demonstrated on an F/A-18 Stabilator using NASTD (an in-house McDonnell Douglas Aerospace East CFD code) coupled with NASTRAN. The process is also demonstrated on an aeroelastic research wing (ARW-2) using ENSAERO (an in-house NASA Ames Research Center CFD code) coupled with a finite element wing-box structures code. Good results have been obtained for the F/A-18 Stabilator while results for the ARW-2 supercritical wing are still being obtained.

Description: An experimental investigation of the effect of leading-edge radius, camber, Reynolds number, and boundary-layer state on the incipient separation of a delta wing at supersonic speeds was conducted at the Langley Unitary Plan Wind Tunnel at Mach number of 1.60 over a free-stream Reynolds number range of 1 x 106 to 5 x 106 ft-1. The three delta wing models examined had a 65 deg swept leading edge and varied in cross-sectional shape: a sharp wedge, a 20:1 ellipse, and a 20:1 ellipse with a -9.750 circular camber imposed across the span. The wings were tested with and without transition grit applied. Surface-pressure coefficient data and flow-visualization data indicated that by rounding the wing leading edge or cambering the wing in the spanwise direction, the onset of leading-edge separation on a delta wing can be raised to a higher angle of attack than that observed on a sharp-edged delta wing. The data also showed that the onset of leading-edge separation can be raised to a higher angle of attack by forcing boundary-layer transition to occur closer to the wing leading edge by the application of grit or the increase in free-stream Reynolds number.

Description: This paper presents a survey of the effects of Reynolds number on the low- speed lift characteristics of wings encountering separated flows at their leading and side edges, with emphasis on the region near the stall. The influence of leading-edge profile and Reynolds number on the stall characteristics of two- dimensional airfoils are reviewed first to provide a basis for evaluating three- dimensional effects associated with various wing planforms. This is followed by examples of the effects of Reynolds number and geometry on the lift characteristics near the stall for a series of three-dimensional wings typical of those suitable for high-speed aircraft and missiles. Included are examples of the effects of wing geometry on the onset and spanwise progression of turbulent reseparation near the leading edge and illustrations of the degree to which simplified theoretical approaches can be useful in defining the influence of the various geometric parameters. Also illustrated is the manner in which the Reynolds number and wing geometry parameters influence whether the turbulent reseparation near the leading edge results in a sudden loss of lift, as in the two-dimensional case, or the formation of a leading-edge vortex with Rs increase in lift followed by a gentle stall as in the highly swept wing case. Particular emphasis is placed on the strong influence of 'induced camber' on the development of turbulent reseparation. R is believed that the examples selected for this report may be useful in evaluating viscous flow solutions by the new computational methods based on the Navier-Stokes equations as well as defining fruitful research areas for the high-Reynolds-number wind tunnels.

Description: In the work reported herein, a simplified, uncoupled, zonal procedure is utilized to assess the capability of numerically simulating icing effects on a Boeing 727-200 aircraft. The computational approach combines potential flow plus boundary layer simulations by VSAERO for the un-iced aircraft forces and moments with Navier-Stokes simulations by NPARC for the incremental forces and moments due to iced components. These are compared with wind tunnel force and moment data, supplied by the Boeing Company, examining longitudinal flight characteristics. Grid refinement improved the local flow features over previously reported work with no appreciable difference in the incremental ice effect. The computed lift curve slope with and without empennage ice matches the experimental value to within 1%, and the zero lift angle agrees to within 0.2 of a degree. The computed slope of the un-iced and iced aircraft longitudinal stability curve is within about 2% of the test data. This work demonstrates the feasibility of a zonal method for the icing analysis of complete aircraft or isolated components within the linear angle of attack range. In fact, this zonal technique has allowed for the viscous analysis of a complete aircraft with ice which is currently not otherwise considered tractable.

Description: In this video a wind tunnel investigation of aerodynamic and structural deflection characteristics of an inflatable airplane is shown. The film includes scenarios during wind tunnel tests of an inflatable airplane in the Langley Full Scale Tunnel with the main objective of obtaining load factors prior to wing buckle of 4.5 to 5.0 g. The inflation pressure during the test was indicated to be 7.0 psi.

Description: To determine the flow field characteristics of 12 planform geometries, a flow visualization investigation was conducted in the Langley 16- by 24-Inch Water Tunnel. Concepts studied included flat plate representations of diamond wings, twin bodies, double wings, cutout wing configurations, and serrated forebodies. The off-surface flow patterns were identified by injecting colored dyes from the model surface into the free-stream flow. These dyes generally were injected so that the localized vortical flow patterns were visualized. Photographs were obtained for angles of attack ranging from 10' to 50', and all investigations were conducted at a test section speed of 0.25 ft per sec. Results from the investigation indicate that the formation of strong vortices on highly swept forebodies can improve poststall lift characteristics; however, the asymmetric bursting of these vortices could produce substantial control problems. A wing cutout was found to significantly alter the position of the forebody vortex on the wing by shifting the vortex inboard. Serrated forebodies were found to effectively generate multiple vortices over the configuration. Vortices from 65' swept forebody serrations tended to roll together, while vortices from 40' swept serrations were more effective in generating additional lift caused by their more independent nature.

Description: An experimental investigation of the effects of angle of attack on hypersonic boundary-layer stability on a flared-cone model was conducted in the low-disturbance Mach-6 Nozzle-Test Chamber Facility at NASA Langley Research Center. This unique facility provided a 'quiet' flow test environment which is well suited for stability experiments because the low levels of freestream 'noise' minimize artificial stimulation of flow-disturbance growth. Surface pressure and temperature measurements documented the adverse-pressure gradient and transition-onset location. Hot-wire anemometry diagnostics were applied to identify the instability mechanisms which lead to transition. In addition, the mean flow over the flared-cone geometry was modeled by laminar Navier-Stokes computations. Results show that the boundary layer becomes more stable on the windward ray and less stable on the leeward ray relative to the zero-degree angle-of-attack case. The second-mode instability dominates the transition process at a zero-degree angle of attack, however, on the windward ray at an angle of attack this mode was completely stabilized. The less-dominant first-mode instability was slightly destabilized on the windward ray. Non-linear mechanisms such as saturation and harmonic generation are identified from the flow-disturbance bispectra.

Description: Solutions of aerodynamic characteristics are presented for the Galileo Probe entering Jupiter's hydrogen-helium atmosphere at a nominal relative velocity of 47.4 km/s. Focus is on predicting the aerodynamic drag coefficient during the transitional flow regime using the direct simulation Monte Carlo (DSMC) method. Accuracy of the probe's drag coefficient directly impacts the inferred atmospheric properties that are being extracted from the deceleration measurements made by onboard accelerometers as part of the Atmospheric Structure Experiment. The range of rarefaction considered in the present study extends from the free molecular limit to continuum conditions. Comparisons made with previous calculations and experimental measurements show the present results for drag to merge well with Navier-Stokes and experimental results for the least rarefied conditions considered.

Description: A computational investigation was conducted to support the development of a semi-span model test capability in the NASA LaRC's National Transonic Facility. This capability is required for the testing of high-lift systems at flight Reynolds numbers. A three-dimensional Navier-Stokes solver was used to compute the low-speed flow over both a full-span configuration and a semi-span configuration. The computational results were found to be in good agreement with the experimental data. The computational results indicate that the stand-off height has a strong influence on the flow over a semi-span model. The semi-span model adequately replicates the aerodynamic characteristics of the full-span configuration when a small stand-off height, approximately twice the tunnel empty sidewall boundary layer displacement thickness, is used. Several active sidewall boundary layer control techniques were examined including: upstream blowing, local jet blowing, and sidewall suction. Both upstream tangential blowing, and sidewall suction were found to minimize the separation of the sidewall boundary layer ahead of the semi-span model. The required mass flow rates are found to be practicable for testing in the NTF. For the configuration examined, the active sidewall boundary layer control techniques were found to be necessary only near the maximum lift conditions.

Description: An experimental wind tunnel test of a 65 deg delta wing model with interchangeable leading edges was conducted in the Langley National Transonic Facility (NTF). The objective was to investigate the effects of Reynolds and Mach numbers on slender-wing leading-edge vortex flows with four values of wing leading-edge bluntness. Experimentally obtained pressure data are presented without analysis in tabulated and graphical formats across a Reynolds number range of 6 x 10(exp 6) to 36 x 10(exp 6) at a Mach number of 0.85 and across a Mach number range of 0.4 to 0.9 at a Reynolds number of 6 x 10(exp 6). Normal-force and pitching-moment coefficient plots for these Reynolds number and Mach number ranges are also presented.

Description: One of the primary reasons for developing quiet tunnels is for the investigation of high-speed boundary-layer stability and transition phenomena without the transition-promoting effects of acoustic radiation from tunnel walls. In this experiment, a flared-cone model under adiabatic- and cooled-wall conditions was placed in a calibrated, 'quiet' Mach 6 flow and the stability of the boundary layer was investigated using a prototype constant-voltage anemometer. The results were compared with linear-stability theory predictions and good agreement was found in the prediction of second-mode frequencies and growth. In addition, the same 'N=10' criterion used to predict boundary-layer transition in subsonic, transonic, and supersonic flows was found to be applicable for the hypersonic flow regime as well. Under cooled-wall conditions, a unique set of continuous spectra data was acquired that documents the linear, nonlinear, and breakdown regions associated with the transition of hypersonic flow under low-noise conditions.

Description: A wind-tunnel investigation has been performed to determine the effect of winglets on the induced drag of a low-aspect-ratio wing configuration at Mach numbers between 0.30 and 0.85 and a nominal angle-of-attack range from -2 deg to 20 deg. Results of the tests at the cruise lift coefficient showed significant increases in lift-drag ratio for the winglet configuration relative to a wing-alone configuration designed for the same lift coefficient and Mach number. Further, even larger increases in lift-drag ratio were observed at lift coefficients above the design value at all Mach numbers tested. The addition of these winglets had a negligible effect on the static lateral-directional stability characteristics of the configuration. No tests were made to determine the effect of these winglets at supersonic Mach numbers, where increases in drag caused by winglets might be more significant. Computational analyses were also performed for the two configurations studied. Linear and small-disturbance formulations were used. The codes were found to give reasonable performance estimates sufficient for predicting changes of this magnitude.

Description: Despite the 80-year history of classical wing theory, considerable research has recently been directed toward planform and wake effects on induced drag. Nonlinear interactions between the trailing wake and the wing offer the possibility of reducing drag. The nonlinear effect of compressibility on induced drag characteristics may also influence wing design. This thesis deals with the prediction of these nonlinear aspects of induced drag and ways to exploit them. A potential benefit of only a few percent of the drag represents a large fuel savings for the world's commercial transport fleet. Computational methods must be applied carefully to obtain accurate induced drag predictions. Trefftz-plane drag integration is far more reliable than surface pressure integration, but is very sensitive to the accuracy of the force-free wake model. The practical use of Trefftz plane drag integration was extended to transonic flow with the Tranair full-potential code. The induced drag characteristics of a typical transport wing were studied with Tranair, a full-potential method, and A502, a high-order linear panel method to investigate changes in lift distribution and span efficiency due to compressibility. Modeling the force-free wake is a nonlinear problem, even when the flow governing equation is linear. A novel method was developed for computing the force-free wake shape. This hybrid wake-relaxation scheme couples the well-behaved nature of the discrete vortex wake with viscous-core modeling and the high-accuracy velocity prediction of the high-order panel method. The hybrid scheme produced converged wake shapes that allowed accurate Trefftz-plane integration. An unusual split-tip wing concept was studied for exploiting nonlinear wake interaction to reduced induced drag. This design exhibits significant nonlinear interactions between the wing and wake that produced a 12% reduction in induced drag compared to an equivalent elliptical wing at a lift coefficient of 0.7. The performance of the split-tip wing was also investigated by wing tunnel experiments. Induced drag was determined from force measurements by subtracting the estimated viscous drag, and from an analytical drag-decomposition method using a wake survey. The experimental results confirm the computational prediction.

Description: Under the High Speed Research (HSR) Program, the NASA Lewis Research Center, the NASA Langley Research Center, and the aerospace industry have been developing exhaust nozzle concepts for a future High Speed Civil Transport (HSCT). For such an aircraft to be both economically viable and environmentally acceptable, the exhaust nozzles must combine highly efficient operation throughout the flight envelope with low noise levels at takeoff. A team composed of researchers from Lewis, Langley, and McDonnell Douglas created a two-dimensional data base using computational fluid dynamics. First, each group validated their computational fluid dynamics code against existing experimental data to gain confidence in the code's results. Next, the organizations analyzed several configurations, and the results were used to create the new data base. Then, the analysis methods were updated to account for the effects of two-dimensional nozzles. This will result in a more accurate analysis of the propulsion system for the High Speed Civil Transport.

Description: The NASA Lewis Research Center is investigating the wave rotor for use as a core gas generator in future gas turbine engines. The device, which uses gas-dynamic waves to transfer energy directly to and from the working fluid through which the waves travel, consists of a series of constant-area passages that rotate about an axis. Through rotation, the ends of the passages are periodically exposed to various circumferentially arranged ports that initiate the traveling waves within the passages.

Description: The design and implementation of Doppler Global Velocimetry (DGV) for testing in the Langley Unitary Plan Wind Tunnel is presented. The discussion begins by outlining the characteristics of the tunnel and the test environment, with potential problem areas highlighted. Modifications to the optical system design to implement solutions for these problems are described. Since this tunnel entry was the first ever use of DGV in a supersonic wind tunnel, the test series was divided into three phases, each with its own goal. Phase I determined if condensation provided sufficient scattered light for DGV applications. Phase II studied particle lag by measuring the flow about an oblique shock above an inclined flat plate. Phase III investigated the supersonic vortical flow field above a 75-degree delta wing at 24-degrees angle of attack. Example results from these tests are presented.

Description: As opposed to previous explanations based on the effects of anharmonicity of simple diatomic molecules, traces of water vapor are suggested to be the most likely cause of the anomalously fast vibrational relaxation of such gases observed in supersonic and hypersonic nozzles. The mechanism is the strong V-VR coupling with H2O molecules that dramatically facilitates the collisional transfer of vibrational energy. Slight moisture content is thus a real world aspect of gas dynamics that must be considered in characterizations of shock tubes, reflected shock tunnels, and expansion tubes.

Description: The development of advanced-design ultrahigh bypass ratio engines has led to renewed interest in the study of the flutter of bladed disks. Previously, two fundamental approaches were used in flutter calculations: frequency domain analysis and time-domain analysis. With the development of time-marching computational fluid dynamics (CFD) flow solvers, both approaches have been used with equal ease. In the present work at the NASA Lewis Research Center, substantial computational savings have been achieved by applying a numerical eigensolver to a nonlinear, time-marching fluid-structure interaction system solver for flutter prediction.

Description: The National PARC (NPARC) Alliance was established by the NASA Lewis Research Center and the Air Force Arnold Engineering Development Center to provide the U.S. aeropropulsion community with a reliable Navier-Stokes code for simulating the nonrotating components of propulsion systems. Recent improvements to the turbulence model capabilities of the NPARC code have significantly improved its capability to simulate turbulent flows. Specifically, the Chien k-epsilon and Wilcox k-omega turbulence models were implemented at Lewis. Lewis researchers installed the Chien k-epsilon model into NPARC to improve the code's ability to calculate turbulent flows with attached wall boundary layers and free shear layers. Calculations with NPARC have demonstrated that the Chien k-epsilon model provides more accurate calculations than those obtained with algebraic models previously available in the code. Grid sensitivity investigations have shown that computational grids must be packed against the solid walls such that the first point off of the wall is placed in the laminar sublayer. In addition, matching the boundary layer and momentum thicknesses entering mixing regions is necessary for an accurate prediction of the free shear-layer growth.

Description: In a continuing in-house program at the NASA Lewis Research Center, analytical and numerical methods are being developed to apply radiative analysis to predict transient temperature distributions and heat flows in partially transmitting materials. Results have been obtained for a single plane layer, and a transient analysis is being developed for a two-layer composite where each layer has a different refractive index. Because the ceramic refractive indices are larger than one, internal reflections are produced at the surfaces and at the internal interface. Reflections tend to distribute energy within a layer, and this affects the transient temperature distributions.

Description: A gas-jet diffusion flame is similar to the flame on a Bunsen burner, where a gaseous fuel (e.g., propane) flows from a nozzle into an oxygen-containing atmosphere (e.g., air). The difference is that a Bunsen burner allows for (partial) premixing of the fuel and the air, whereas a diffusion flame is not premixed and gets its oxygen (principally) by diffusion from the atmosphere around the flame. Simple gas-jet diffusion flames are often used for combustion studies because they embody the mechanisms operating in accidental fires and in practical combustion systems. However, most practical combustion is turbulent (i.e., with random flow vortices), which enhances the fuel/air mixing. These turbulent flames are not well understood because their random and transient nature complicates analysis. Normal gravity studies of turbulence in gas-jet diffusion flames can be impeded by buoyancy-induced instabilities. These gravitycaused instabilities, which are evident in the flickering of a candle flame in normal gravity, interfere with the study of turbulent gas-jet diffusion flames. By conducting experiments in microgravity, where buoyant instabilities are avoided, we at the NASA Lewis Research Center hope to improve our understanding of turbulent combustion. Ultimately, this could lead to improvements in combustor design, yielding higher efficiency and lower pollutant emissions. Gas-jet diffusion flames are often researched as model flames, because they embody mechanisms operating in both accidental fires and practical combustion systems (see the first figure). In normal gravity laboratory research, buoyant air flows, which are often negligible in practical situations, dominate the heat and mass transfer processes. Microgravity research studies, however, are not constrained by buoyant air flows, and new, unique information on the behavior of gas-jet diffusion flames has been obtained.

Description: The Interface Configuration Experiment (ICE) is actually a series of experiments that explore the striking behavior of liquid-vapor interfaces (i.e., fluid surfaces) in a low gravity environment under which major shifts in liquid position can arise from small changes in container shape or contact angle. Although these experiments are designed to test current mathematical theory, there are numerous practical applications that could result from these studies. When designing fluid management systems for space-based operations, it is important to be able to predict the locations and configurations that fluids will assume in containers under low-gravity conditions. The increased ability to predict, and hence control, fluid interfaces is vital to systems and/or processes where capillary forces play a significant role both in space and on the Earth. Some of these applications are in general coating processes (paints, pesticides, printing, etc.), fluid transport in porous media (ground water flows, oil recovery, etc.), liquid propellant systems in space (liquid fuel and oxygen), capillary-pumped loops and heat pipes, and space-based life-support systems. In space, almost every fluid system is affected, if not dominated, by capillarity. Knowledge of the liquid-vapor interface behavior, and in particular the interface shape from which any analysis must begin, is required as a foundation to predict how these fluids will react in microgravity and on Earth. With such knowledge, system designs can be optimized, thereby decreasing costs and complexity, while increasing performance and reliability. ICE has increased, and will continue to increase this knowledge, as it probes the specific peculiarities of current theory upon which our current understanding of these effects is based. Several versions of ICE were conducted in NASA Lewis Research Center's drop towers and on the space shuttle during the first and second United States Microgravity Laboratory missions (USML-1 and USML-2). Additional tests are planned for the space shuttle and for the Russian Mir space station. These studies will focus on interfacial problems concerning surface existence, uniqueness, configuration, stability, and flow characteristics.

Description: The development of thermodynamics and statistical mechanics is very important in the history of physics, and it underlines the difficulty in dealing with systems involving many bodies, even if those bodies are identical. Macroscopic systems of atoms typically contain so many particles that it would be virtually impossible to follow the behavior of all of the particles involved. Therefore, the behavior of a complete system can only be described or predicted in statistical ways. Under a grant to the NASA Lewis Research Center, scientists at the Case Western Reserve University have been examining the use of modern computing techniques that may be able to investigate and find the behavior of complete systems that have a large number of particles by tracking each particle individually. This is the study of molecular dynamics. In contrast to Monte Carlo techniques, which incorporate uncertainty from the outset, molecular dynamics calculations are fully deterministic. Although it is still impossible to track, even on high-speed computers, each particle in a system of a trillion trillion particles, it has been found that such systems can be well simulated by calculating the trajectories of a few thousand particles. Modern computers and efficient computing strategies have been used to calculate the behavior of a few physical systems and are now being employed to study important problems such as supersonic flows in the laboratory and in space. In particular, an animated video (available in mpeg format--4.4 MB) was produced by Dr. M.J. Woo, now a National Research Council fellow at Lewis, and the G-VIS laboratory at Lewis. This video shows the behavior of supersonic shocks produced by pistons in enclosed cylinders by following exactly the behavior of thousands of particles. The major assumptions made were that the particles involved were hard spheres and that all collisions with the walls and with other particles were fully elastic. The animated video was voted one of two winning videos in a competition held at the meeting of the American Physical Society's Division of Fluid Dynamics, held in Atlanta, Georgia, in November 1994. Of great interest was the result that in every shock there were a few high-speed precursor particles racing ahead of the shock, carrying information about its impending arrival. Most recently, Dr. Woo has been applying molecular dynamics techniques to the problem of determining the drag produced by the space station truss structure as it flies through the thin residual atmosphere of low-Earth orbit. This problem is made difficult by the complex structure of the truss and by the extreme supersonic nature of the flow. A fully filled section of the truss has already been examined, and drag predictions have been made. Molecular dynamics techniques promise to make realistic drag calculations possible even for very complex partially filled truss segments flying at arbitrary angles.

Description: We have developed a calculational model that treats all the components of an orifice pulse tube cooler. We base our analysis on 1-dimensional thermodynamic equations for the regenerator and we assume that all mass flows, pressure oscillations and temperature oscillations are small and sinusoidal. Non-linear pressure drop effects are included in the regenerator to account for finite pressure amplitude effects. The resulting mass flows and pressures are matched at the boundaries with the other components of the cooler: compressor, aftercooler, cold heat exchanger, pulse tube, hot heat exchanger, orifice and reservoir. The results of the calculation are oscillating pressures, mass flows and enthalpy flows in the main components of the cooler. By comparing with the calculations of other available models, we show that our model is very similar to REGEN 3 from NIST and DeltaE from Los Alamos National Lab. Our model is much easier to use than other available models because of its simple graphical interface and the fact that no guesses are required for the operating pressures or mass flows. In addition, the model only requires a few minutes of running time allowing many parameters to be optimized in a reasonable time. A version of the model is available for use over the World Wide Web at http://irtek.arc.nasa.gov. Future enhancements include adding a bypass orifice and including second order terms in steady mass streaming and steady heat transfer. A two-dimensional anelastic approximation of the fluid equations will be used as the basis for the latter analysis. Preliminary results are given in dimensionless numbers appropriate for oscillating compressible flows. The model shows how transverse heat transfer reduces enthalpy flow, particularly for small pulse tubes. The model also clearly shows mass recirculation in the open tube on the order of the tube length. They result from the higher order Reynolds stresses. An interesting result of the linearized approach is that the steady mass streaming does not affect the enthalpy flow at second order. The major effect of recirculating mass streaming is to increase transverse temperature gradients, which leads to higher entropy production and reduced efficiency.

Description: The Harmonic Point Source (HPS) is the simplest possible form of 3D disturbance yet even this case is not fully understood. One of the characteristics of HPS generated Tolmien Schlichting (TS) waves is that the maximum rms amplitude downstream of the source occurs away from the centerline. As the waves fan out in the spanwise direction they can encounter Klebanoff modes (weak streamwise vortices) which may be present as background disturbances within the layer. The vortices originate near the leading edge and they appear to be caused by amplification of almost immeasurably small nonuniformities in the free stream. The vortices can be steady but they often appear as very low frequency background unsteadiness (as observed by Klebanoff). The vortices locally distort the mean flow and therefore (by definition) they are nonlinear phenomena e.g. the growth rate of the TS waves is altered. Kendall (private communication) has demonstrated total suppression of TS waves by deliberately introducing strong vortices generated by a delta wing. For the very weak vortices that are often encountered in experiments, their effect is characterized by a local reduction in the rms TS magnitude. The laminar wake behind a fine wire is used to generate Klebanoff modes and the interactions with HPS generated TS waves are explored in a controlled manner. It is shown that a proportion of the reduction in the observed TS magnitude can be attributed to washout of the phase-averaged signals owing to phase jitter.

Description: Calculations have been carried out on the adiabatic laminar boundary layer developing on the surface of a two-dimensional supersonic nozzle, consisting of contoured nozzle blocks and flat sidewalls. Two- and three-dimensional Navier-Stokes codes, as well as two-dimensional boundary-layer codes have been employed. Thee codes have been adapted to the characteristics of a specific wind tunnel nozzle, so that their numerical results could be directly compared with experimental data obtained in the same nozzle. Such comparisons have been made for the boundary-layer growth on the contoured nozzle, and for the boundary-layer growth, surface streamlines and surface shear on the sidewalls. The three-dimensional Navier-Stokes code was found to be the only one to correctly predict the mean boundary-layer flow on both the sidewalls and the contoured nozzle. Theory and experiment both indicate that the sidewall flow is highly three-dimensional, with non-uniform shear, comer vortices and a boundary layer strongly distorted by cross flows induced by lateral pressure gradients.

Description: RANS-MP, a new implementation of a single-grid Navier-Stokes solver using the diagonalized Beam-Warming approximate-factorization scheme, is presented. This first release of the completely rewritten solver employs the following optimizations: (1) Bi-directional multi-partition method for the ADI solver part; this improves granularity and load balance; (2) Improved cache usage through elimination of non-unit-stride array access (possible in part due to multi-partitioning); (3) Preprocessing of communicating boundary conditions to streamline logic during time stepping; (4) Truly parallel, high-performance I/O using the newly-developed MPI-IO library; (5) Elimination of large amounts of redundant operations through efficient use of workspace. Results of some realistic wing computations on the IBM SP2 computer will be presented. We will demonstrate that excellent absolute performance and scalability are obtained with RANS-MP, even for relatively small grid sizes. Besides high performance, an outstanding feature of RANS-MP is its true portability, due to the use of the portable message passing and I/O libraries MPI and MPI-IO.