SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
This report deals with the flow past an infinite screen of thin airfoil (two-dimensional problem). The vortex distribution across the profile is established with appropriate expansion in series and the velocity distribution lift, moment, and profile shape deduced. Inversely, the distribution is deduced from the vorticity. The method is the extension of the Birnbaum-Glauert method for the isolated wing.
Attached as appendices to this report are documents describing work performed on the simulation of a landing powered-lift delta wing, the tracking of flow features using overset grids, and the simulation of flaps on the Wright Patterson Lab's fighter-lift-and-control (FLAC) wing. Numerical simulation of a powered-lift landing includes the computation of flow about a delta wing at four fixed heights as well as a simulated landing, in which the delta wing descends toward the ground. Comparison of computed and experimental lift coefficients indicates that the simulations capture the qualitative trends in lift-loss encountered by thrust-vectoring aircraft operating in ground effect. Power spectra of temporal variations of pressure indicate computed vortex shedding frequencies close to the jet exit are in the experimentally observed frequency range; the power spectra of pressure also provide insights into the mechanisms of lift oscillations. Also, a method for using overset grids to track dynamic flow features is described and the method is validated by tracking a moving shock and vortices shed behind a circular cylinder. Finally, Chimera gridding strategies were used to develop pressure coefficient contours for the FLAC wing for a Mach no. of 0.18 and Reynolds no. of 2.5 million.
Polar statospheric clouds (PSC) and stratospheric aerosol can be observed by ground-based UV-visible spectrometry by looking at the variation of the color of the sky during twilight. A radiative transfer model shows that reddenings are caused by high altitude (22-28 km) thin layers of scatterers, while low altitude (12-20 km) thick ones result in blueings. The color index method applied on 4 years of observations at Dumont d'Urville (67 deg S), from 1988 to 1991, shows that probably because the station is located at the edge of the vortex, dense PSC are uncommon. More unexpected is the existence of a systematic seasonal variation of the color of the twilight sky - bluer at spring - which reveals the formation of a dense scattering layer at or just above the tropopause at the end of the winter. Large scattering layers are reported above the station in 1991, first in August around 12-14 km, later in September at 22-24 km. They are attributed to volcanic aerosol from Mt Hudson and Mt Pinatubo respectively, which erupted in 1991. Inspection of the data shows that the lowest entered rapidly into the polar vortex but not the highest which remained outside, demonstrating that the vortex was isolated at 22-26 km.
Strong interactions can occur between the flow about an aerospace vehicle and its structural components resulting in several important aeroelastic phenomena. These aeroelastic phenomena can significantly influence the performance of the vehicle. At present, closed-form solutions are available for aeroelastic computations when flows are in either the linear subsonic or supersonic range. However, for aeroelasticity involving complex nonlinear flows with shock waves, vortices, flow separations, and aerodynamic heating, computational methods are still under development. These complex aeroelastic interactions can be dangerous and limit the performance of aircraft. Examples of these detrimental effects are aircraft with highly swept wings experiencing vortex-induced aeroelastic oscillations, transonic regime at which the flutter speed is low, aerothermoelastic loads that play a critical role in the design of high-speed vehicles, and flow separations that often lead to buffeting with undesirable structural oscillations. The simulation of these complex aeroelastic phenomena requires an integrated analysis of fluids and structures. This report presents a summary of the development, applications, and procedures to use the multidisciplinary computer code ENSAERO. This code is based on the Euler/Navier-Stokes flow equations and modal/finite-element structural equations.
The key to accurate predictions of rotorcraft aerodynamics, acoustics, and dynamics lies in the accurate representation of the rotor wake. The vortical wake computed by rotorcraft CFD analyses typically suffer from numerical dissipation before the first blade passage. With some a priori knowledge of the wake trajectory, grid points can be concentrated along the trajectory to minimize the dissipation. Comprehensive rotorcraft analyses based on lifting-line theory rely on classical vortex models and/or semi-empirical information about the tip vortex structure. Until the location, size, and strength of the trailed tip vortex can be measured over a range of wake ages, the analyses will continue to be adjusted on a trial and error basis in order to correctly predict blade airloads, acoustics, dynamics, and performance. Using the laser light sheet technique, tip vortex location can be acquired in a straightforward manner. Measuring wake velocities and vortex core size, however, has been difficult and tedious using point-measurement techniques such as laser velocimetry. Recently, the Particle Image Velocimetry (PIV) technique has proven to be an efficient method for acquiring velocity measurements over relatively large areas and volumes of a rotor wake. The work reported to date, however, has been restricted to 2-component velocity measurements of the rotor wake. Three-component velocity measurements of a hovering rotor wake were acquired at NASA Ames Research Center in May 1999. This experiment represents a major step toward understanding the detailed structure of a rotor wake. This paper will focus primarily on the experimental technique used in acquiring this data. The accuracy and limitations of the current technique will also be discussed. Representative velocity field measurements will be included.
A proposed flight test program to measure the characteristics of wake vortices behind a T-33 aircraft was investigated. A number of facets of the flight tests were examined to define the parameters to be measured, the anticipated vortex characteristics, the mutual interference between the probe aircraft and the wake, the response of certain instruments to be used in obtaining measurements, the effect of condensation on the wake vortices, and methods of data reduction. Recommendations made as a result of the investigation are presented.
An attempt is made to summarize the state of rotary wing flow fields. The theory is traced from its origin as a momentum theory estimate of average interference, through simple vortex theory, to its present status where it is indispensible in calculating blade loads. A comparison of the theory with flow measurements are presented. The modern efforts toward using more detailed digital methods to obtain blade load distribution are described.
A numerical simulation of the final stages of transition to turbulence in plane channel flow is reported. Three dimensional, incompressible Navier-Stokes equations are numerically integrated to obtain the time-evolution of two and three dimensional finite amplitude disturbances. Computations are performed on the CYBER-203 vector processor for a 32x51x32 grid. Results are presented for no-slip boundary conditions at the solid walls as well as for periodic suction-blowing to simulate active control of transition by mass transfer. Solutions indicate that the method is capable of simulating the complex character of vorticity dynamics during the various stages of transition and final breakdown. In particular, evidence points to the formation of a lambda-shape vortex and the subsequent system of horseshoe vortices inclined to the main flow direction as the main elements of transition. Calculations involving suction-blowing indicate that interference with a wave of suitable phase and amplitude reduces the disturbance growth rates.
Results obtained from a numerical simulation of the final stages of transition to turbulence in plane channel flow are described. Three dimensional, incompressible Navier-Stokes equations are numerically integrated to obtain the time evolution of two and three dimensional finite amplitude disturbances. Computations are performed on CYBER-203 vector processor for a 32x51x32 grid. Results are presented for no-slip boundary conditions at the solid walls as well as for periodic suction blowing to simulate active control of transition by mass transfer. Solutions indicate that the method is capable of simulating the complex character of vorticity dynamics during the various stages of transition and final breakdown. In particular, evidence points to the formation of a lambda-shape vortex and the subsequent system of horseshoe vortices inclined to the main flow direction as the main elements of transition. Calculations involving periodic suction-blowing indicate that interference with a wave of suitable phase and amplitude reduces the disturbance growth rates.
Theoretical and experimental research on unsteady aerodynamic flows is discussed. Complex flow fields that involve separations, vortex interactions, and transonic flow effects were investigated. Flow visualization techniques are used to obtain a global picture of the flow phenomena before detailed quantitative studies are undertaken. A wide variety of methods are used to visualize fluid flow and a sampling of these methods is presented. It is emphasized that the visualization technique is a thorough quantitative analysis and subsequent physical understanding of these flow fields.
A three-dimensional, quasi-steady, full potential flow solver was adapted to include realistic rotor wake influence for the aerodynamic analysis of helicopter rotors. The method uses an inner/outer domain technique to accommodate wake effects. Nonlinear flow is computed in the inner domain using a finite difference solution method. The wake is modeled using prescribed wake techniques to allow for the inclusion of realistic wake geometries. Portions of the wake passing inside the inner domain are treated using an embedded vortex technique. The procedure couples the wake influence with the inner domain solution in a consistent and efficient solution process. Correlation with measured lifting transonic data in hover and forward flight is shown which demonstrates the merits of the approach.
The properties of vortical flows are studied using both theoretical analysis and computational flow fields. The consequences of two definitions of the vortex core, a minimum in the streamline curvature and a maximum in the normalized helicity, are examined. Analysis indicates that several criteria must be met if the cores defined by these two methods are to coincide. In certain regions of the flow, computational flow fields indicate that these two definitions are coincident, that the velocity and vorticity fields are aligned at the core, and that extrema in the velocity magnitude, vorticity magnitude, pressure, and density occur at the core.
A pressure-sensitive paint (PSP) technique was applied in a wind tunnel experiment in the NASA Langley Research Center 8-Foot Transonic Pressure Tunnel to quantify the effect of wing fillets on the global vortex-induced surface static pressure field about a sharp leading-edge 76o/40o double delta wing, or strake-wing, model at subsonic and transonic speeds. Global calibrations of the PSP were obtained at M = 0.50, 0.70, 0.85, 0.95, and 1.20, a Reynolds number per unit length of 2.0 million, and angles of attack from 10 degrees to 20 degrees using an in-situ method featuring the simultaneous acquisition of electronically-scanned pressures (ESP) at discrete locations on the model. The mean error in the PSP measurements relative to the ESP data was approximately 2 percent or less at M = 0.50 to 0.85 but increased to several percent at M =0.95 and 1.20. The PSP pressure distributions and pseudo-colored planform view pressure maps clearly revealed the vortex-induced pressure signatures at all Mach numbers and angles of attack. Small fillets having a parabolic or diamond planform situated at the strake-wing intersection were designed to manipulate the vortical flows by, respectively, removing the leading-edge discontinuity or introducing additional discontinuities. The fillets caused global changes in the vortex-dominated surface pressure field that were effectively captured in the PSP measurements. The vortex surface pressure signatures were compared to available off-surface vortex cross-flow structures obtained using a laser vapor screen (LVS) flow visualization technique. The fillet effects on the PSP pressure distributions and the observed leading-edge vortex flow characteristics were consistent with the trends in the measured lift, drag, and pitching moment coefficients.
A computational system for estimation of nonlinear aerodynamic characteristics of wings at supersonic speeds was developed and was incorporated in a computer program. This corrected linearized theory method accounts for nonlinearities in the variation of basic pressure loadings with local surface slopes, predicts the degree of attainment of theoretical leading edge thrust, and provides an estimate of detached leading edge vortex loadings that result when the theoretical thrust forces are not fully realized.
The capability of two numerical methods, one for transonic and one for supersonic flows, to predict the flow fields about representative fighter aircraft forebodies in the vicinity of the engine inlets was examined. The Mach number range covered was 0.9 to 2.5 and the angle-of-attack range was 0 deg to 25 deg. The computer progams that implement each of the numerical methods are described as to their features and usage, and results are compared with comprehensive wind tunnel data. Although both prediction methods were inviscid, results show that the aerodynamic effects of the forebody, with and without a wing, can be simulated fairly well. Futher work is needed to include the effects of viscosity, including vortex shedding.
The ability of a low-order panel method to calculate the aerodynamic loads on wings caused by interaction with wake vortices was studied. The loads were calculated for various positions of a downstream following wing relative to an upstream vortex-generating wing. Calculated vortex-induced span loads and rolling-moment coefficients on the following wing were compared with experimental data. A good agreement with experiment was obtained when the following wing was located more than one following-wing chord length from the tip vortex. The predictions deteriorated as the following wing was placed closer to the vortex. At large downstream distances (approximately 10 generating-wing chord lengths), induced rolling-moment coefficients on the following wing were consistently overestimated. Despite the strong interaction between the wake-vortex filaments and surface doublet panels, the accuracy of the calculations was in most cases independent of the panel distribution and density. A good agreement between theoretical and experimental loads was obtained with a minimum of experimentation with panel arrangements.
A computer program, VSAERO, has been applied to a number of V/STOL configurations with a view to advancing prediction techniques for the low-speed aerodynamic characteristics. The program couples a low-order panel method with surface streamline calculation and integral boundary layer procedures. The panel method--which uses piecewise constant source and doublet panels-includes an iterative procedure for wake shape and models boundary layer displacement effect using the source transpiration technique. Certain improvements to a basic vortex tube jet model were installed in the code prior to evaluation. Very promising results were obtained for surface pressures near a jet issuing at 90 deg from a flat plate. A solid core model was used in the initial part of the jet with a simple entrainment model. Preliminary representation of the downstream separation zone significantly improve the correlation. The program accurately predicted the pressure distribution inside the inlet on the Grumman 698-411 design at a range of flight conditions. Furthermore, coupled viscous/potential flow calculations gave very close correlation with experimentally determined operational boundaries dictated by the onset of separation inside the inlet. Experimentally observed degradation of these operational boundaries between nacelle-alone tests and tests on the full configuration were also indicated by the calculation. Application of the program to the General Dynamics STOL fighter design were equally encouraging. Very close agreement was observed between experiment and calculation for the effects of power on pressure distribution, lift and lift curve slope.
The application of an upwind-biased implicit approximate factorization Navier-Stokes algorithm to the unsteady impulsive start-up flow over a circular cylinder at Reynolds number 1200 is described. The complete form of the compressible Navier-Stokes equations is used, and the algorithm is second-order accurate in both space and time. The development with time of the shape and size of the separated vortical flow region is computed, as well as the time-variation of several boundary layer parameters and profile shapes. Computations, in general, show excellent agreement with experiment, although the present method predicts a more rapid onset of reversed flow on the cylinder than evidenced in experiment. The changes that the vortical region behind the cylinder undergoes as the symmetric flow transitions to periodic vortex shedding are discussed. The flow becomes periodic with a Strouhal frequency of 0.222, which compares well with the experimental value of approximately 0.21. The effect of grid density on the development of the unsteady flow is also shown.