Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “vortex methods”

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.

At least 1,063 records · Page 59

Downwash and Wake Behind Plain and Flapped Airfoils

Extensive experimental measurements have been made of the downwash angles and the wake characteristics behind airfoils with and without flaps and the data have been analyzed and correlated with the theory. A detailed study was made of the errors involved in applying lifting-line theory, such as the effects of a finite wing chord, the rolling-up of the trailing vortex sheet, and the wake. The downwash angles, as computed from the theoretical span load distribution by means of the Biot-Savart equation, were found to be in satisfactory agreement with the experimental results. The rolling-up of the trailing vortex sheet may be neglected, but the vertical displacement of the vortex sheet requires consideration. By the use of a theoretical treatment indicated by Prandtl, it has been possible to generalize the available experimental results so the predictions can be made of the important wake parameters in terms of the distance behind the airfoil trailing edge and the profile-drag coefficient. The method of application of the theory to design and the satisfactory agreement between predicted and experimental results when applied to an airplane are demonstrated.

Silverstein, Abe↗

MFLOP to GFLOP: The Impact on High Fidelity Based Computational Aeroelasticity

Aeroelasticity which involves strong coupling of fluids, structures and controls is an important element in designing an aircraft. Computational aeroelasticity using low fidelity methods such as the linear aerodynamic flow equations coupled with the modal structural equations are well advanced. Though these low fidelity approaches are computationally less intensive, they are not adequate for the analysis of modern aircraft which can experience complex flow/structure interactions. Even at moderate angles of attack supersonic aircraft can experience vortex induced aeroelastic oscillations. Near transonic speeds buffet associated structural oscillations are possible. Aircraft flying in transonic regime may experience a dip in the flutter speed. For accurate aeroelastic computations at these complex fluid/structure interaction situations, high fidelity equations such as the Navier-Stokes for fluids and the finite-elements for structures are needed. Computations using these high fidelity equations require large computational resources both in memory and speed. Current conventional supercomputers have reached their limitations both in memory and speed. As a result, parallel computers have evolved to overco me the limitations of conventional computers. This paper will address the transition that is taking place in computational aeroelasticity from conventional computers to parallel computers. The paper will address special techniques needed to take advantage of the architecture of new parallel computers. Results will be illustrated from computations made on iPSC/860 and IBM SP2 computer by using ENSAERO code that directly couples the Euler/Navier-Stokes flow equations with high resolution finite-element structural equations. Modifications required in both fluids and structural solvers in order to run efficiently on parallel computers will be discussed. Implementation of moving grids and fluid/structural interface on parallel computers will be discussed.

Guruswamy, Guru P.↗

Full scale visualization of the wing tip vortices generated by a typical agricultural aircraft

The trajectories of the wing tip vortices of a typical agricultural aircraft were experimentally determined by flight test. A flow visualization method, similar to the vapor screen method used in wind tunnels, was used to obtain trajectory data for a range of flight speeds, airplane configurations, and wing loadings. Detailed measurements of the spanwise surface pressure distribution were made for all test points. Further, a powered 1/8 scale model of the aircraft was designed, built, and used to obtain tip vortex trajectory data under conditions similar to that of the full scale test. The effects of light wind on the vortices were demonstrated, and the interaction of the flap vortex and the tip vortex was clearly shown in photographs and plotted trajectory data.

Cross, E. J., Jr.↗

Full scale visualization of the wing tip vortices generated by a typical agricultural aircraft

The trajectories of the wing tip vortices of a typical agricultural aircraft were experimentally determined by flight test. A flow visualization method, similar to the vapor screen method used in wind tunnels, was used to obtain trajectory data for a range of flight speeds, airplane configurations, and wing loadings. Detailed measurements of the spanwise surface pressure distribution were made for all test points. Further, a powered 1/8 scale model of the aircraft was designed, built, and used to obtain tip vortex trajectory data under conditions similar to that of the full-scale test. The effects of light wind on the vortices were demonstrated, and the interaction of the flap vortex and the tip vortex was clearly shown in photographs and plotted trajectory data.

Cross, E. J., Jr.↗

Interactions of airfoils with gusts and concentrated vortices in unsteady transonic flow

Unsteady interactions of concentrated vortices and distributed free-stream gusts with a stationary airfoil have been analyzed in two-dimensional transonic flow. A simple method of introducing such disturbances has been implemented numerically in the well-known transonic small-disturbance code LTRAN2, and calculations have been performed for two important classes of current aerodynamic problems. The first, which demonstrates many of the essential features of the interactions between helicopter rotor blades and their trailing-vortex wakes, is that of a discrete potential vortex convecting past an airfoil. The second is the response of a transonic airfoil to a transverse periodic gust, with and without the alleviation that can be achieved by the proper active control motion of a trailing-edge flap. In both cases, unsteady effects are found to play important roles in the shock-wave motion, in the overall flow-field development, and consequently, in the air loads on the airfoil.

Mccroskey, W. J.↗

Reduction of Blade-Vortex Interaction (BVI) noise through X-force control

Momentum theory and the longitudinal force balance equations of a single rotor helicopter are used to develop simple expressions to describe tip-path-plane tilt and uniform inflow to the rotor. The uniform inflow is adjusted to represent the inflow at certain azimuthal locations where strong Blade-Vortex Interaction (BVI) is likely to occur. This theoretical model is then used to describe the flight conditions where BVI is likely to occur and to explore those flight variables that can be used to minimize BVI noise radiation. A new X-force control is introduced to help minimize BVI noise. Several methods of generating the X-force are presented that can be used to alter the inflow to the rotor and thus increasing the likelihood of avoiding BVI during approaches to a landing.

Schmitz, Fredric H.↗

Recovery of the Antarctic Ozone Hole

The Antarctic ozone hole develops each year and culminates by early Spring. Antarctic ozone values have been monitored since 1979 using satellite observations from the TOMS and OMI instruments. The severity of the hole has been assessed using the minimum total ozone value from the October monthly mean (depth of the hole), the average size during the September-October period, and the ozone mass deficit. Ozone is mainly destroyed by halogen catalytic cycles, and these losses are modulated by temperature variations in the collar of the polar lower stratospheric vortex. In this presentation, we show the relationships of halogens and temperature to both the size and depth of the hole. Because atmospheric halogen levels are responding to international agreements that limit or phase out production, the amount of halogens in the stratosphere should decrease over the next few decades. We use two methods to estimate ozone hole recovery. First, we use projections of halogen levels combined with age-of-air estimates in a parametric model. Second, we use a coupled chemistry climate model to assess recovery. We find that the ozone hole is recovering at an extremely slow rate and that large ozone holes will regularly recur over the next 2 decades. Furthermore, full recovery to 1980 levels will not occur until approximately 2068. We will also show some error estimates of these dates and the impact of climate change on the recovery.

Newman, Paul A.↗

Design of Turbomachinery Blading in Three-Dimensional Flow By the Circulation Method: A Progress Report

In this paper, we report on the progress of the further development and application of an analytical theory for the design of turbomachinery blading in three-dimensional flow. After a section describing the recent development of the theory, we discuss the various applications including the investigations of the design of radial inflow turbines. Numerical examples are also presented to illustrate the application of the theory to the design of a free vortex fan, an axial turbine rotor, a free vortex compressor stage for low and high Mach numbers, a non-free vortex fan with trailing vortex sheets, and a counter-rotating propeller in the limit of a large number of blades. Further possible applications of the theory to the design of other aero- and hydrodynamic devices are discussed.

W R Hawthorne↗

Thermal reactor

A thermal reactor apparatus and method of pyrolyticaly decomposing silane gas into liquid silicon product and hydrogen by-product gas is disclosed. The thermal reactor has a reaction chamber which is heated well above the decomposition temperature of silane. An injector probe introduces the silane gas tangentially into the reaction chamber to form a first, outer, forwardly moving vortex containing the liquid silicon product and a second, inner, rewardly moving vortex containing the by-product hydrogen gas. The liquid silicon in the first outer vortex deposits onto the interior walls of the reaction chamber to form an equilibrium skull layer which flows to the forward or bottom end of the reaction chamber where it is removed. The by-product hydrogen gas in the second inner vortex is removed from the top or rear of the reaction chamber by a vortex finder. The injector probe which introduces the silane gas into the reaction chamber is continually cooled by a cooling jacket.

Levin, H.↗

Solution method for a hovering helicopter rotor using the Euler equations

A method for the calculation of the flow field of a helicopter rotor blade in hover is presented. The approach uses a finite volume solution of the three dimensional Euler equations for the blade near field. In the Euler solver the velocity field is decomposed into two parts. One is the induced velocity of the vortex wake extending below the blade, found from a free wake calculation procedure. The other part is the unknown additional velocity field of the rotor blade. This approach eliminates numerical diffusion of the rolled up wake vorticity due to truncation error and artificial viscosity. Also, the effects of the far wake are included in the limited computational domain. Solutions are presented for an isolated wing and a model helicopter rotor and compared to experiment.

Roberts, T. W.↗

Uniform high order spectral methods for one and two dimensional Euler equations

Uniform high order spectral methods to solve multi-dimensional Euler equations for gas dynamics are discussed. Uniform high order spectral approximations with spectral accuracy in smooth regions of solutions are constructed by introducing the idea of the Essentially Non-Oscillatory (ENO) polynomial interpolations into the spectral methods. The authors present numerical results for the inviscid Burgers' equation, and for the one dimensional Euler equations including the interactions between a shock wave and density disturbance, Sod's and Lax's shock tube problems, and the blast wave problem. The interaction between a Mach 3 two dimensional shock wave and a rotating vortex is simulated.

Cai, Wei↗

Transonic aeroelastic computations on wings using Navier-Stokes equations

Aeroelastic computations are made on fighter type wings in the transonic regime. The flow is modeled using the Navier-Stokes equations and is coupled with structural equations of motion. The flow equations are solved by a time accurate finite difference scheme with moving grids. The coupled aeroelastic equations of motion are solved using the linear acceleration method. The configuration adaptive dynamic grids are time-accurately generated using the aeroelastically deformed shape of the wing. Computations are made for oscillating rigid wings with moving shock waves in the presence of leading edge vortices. The computed results compare well with the experiment. Unsteady computations are made to demonstrate the shock-vortex interaction phenomenon on wings in ramp motion. Effects of flexibility and pitch rate are demonstrated for flows with vortices.

Guruswamy, Guru P.↗

Turbulence Measurements and Computations for the Predication of Broadband Noise in High Bypass Ratio Fans

Work was performed under this grant with a view to providing the experimental and computational results needed to improve the prediction of broadband stator noise in large bypass ratio aircraft engines. The central hypothesis of our study was that a large fraction of this noise was generated by the fan tip leakage vortices. More specifically, that these vortices are a significant component of the fan wake turbulence and they contain turbulent eddies of a type that can produce significant broadband noise. To test this hypothesis we originally proposed experimental work and computations with the following objectives: (1) to build a large scale two-dimensional cascade with a tip gap and a stationary endwall that, as far as possible, simulates the fan tip geometry, (2) to build a moving endwall for use with the large scale cascade, (3) to measure, in detail, the turbulence structure and spectrum generated by the blade wake and tip leakage vortex, for both endwall configurations, (4) to use the CFD to compute the flow and turbulence distributions for both the experimental configurations and the ADP fan, (5) to provide the experimental and CFD results for the cascades and the physical understanding gained from their study as a basis for improving the broadband noise prediction method. In large part these objectives have been achieved. The most important achievements and findings of our experimental and computational efforts are summarized below. The bibliography at the end of this report includes a list of all publications produced to date under this project. Note that this list is necessarily incomplete the task of publication (particularly in journal papers) continues.

Devenport, William J.↗

Entropy Splitting and Numerical Dissipation

A rigorous stability estimate for arbitrary order of accuracy of spatial central difference schemes for initial-boundary value problems of nonlinear symmetrizable systems of hyperbolic conservation laws was established recently by Olsson and Oliger (1994) and Olsson (1995) and was applied to the two-dimensional compressible Euler equations for a perfect gas by Gerritsen and Olsson (1996) and Gerritsen (1996). The basic building block in developing the stability estimate is a generalized energy approach based on a special splitting of the flux derivative via a convex entropy function and certain homogeneous properties. Due to some of the unique properties of the compressible Euler equations for a perfect gas, the splitting resulted in the sum of a conservative portion and a non-conservative portion of the flux derivative. hereafter referred to as the "Entropy Splitting." There are several potential desirable attributes and side benefits of the entropy splitting for the compressible Euler equations that were not fully explored in Gerritsen and Olsson. The paper has several objectives. The first is to investigate the choice of the arbitrary parameter that determines the amount of splitting and its dependence on the type of physics of current interest to computational fluid dynamics. The second is to investigate in what manner the splitting affects the nonlinear stability of the central schemes for long time integrations of unsteady flows such as in nonlinear aeroacoustics and turbulence dynamics. If numerical dissipation indeed is needed to stabilize the central scheme, can the splitting help minimize the numerical dissipation compared to its un-split cousin? Extensive numerical study on the vortex preservation capability of the splitting in conjunction with central schemes for long time integrations will be presented. The third is to study the effect of the non-conservative proportion of splitting in obtaining the correct shock location for high speed complex shock-turbulence interactions. The fourth is to determine if this method can be extended to other physical equations of state and other evolutionary equation sets. If numerical dissipation is needed, the Yee, Sandham, and Djomehri (1999) numerical dissipation is employed. The Yee et al. schemes fit in the Olsson and Oliger framework.

Yee, H. C.↗

Lattice Boltzmann Equation On a 2D Rectangular Grid

We construct a multi-relaxation lattice Boltzmann model on a two-dimensional rectangular grid. The model is partly inspired by a previous work of Koelman to construct a lattice BGK model on a two-dimensional rectangular grid. The linearized dispersion equation is analyzed to obtain the constraints on the isotropy of the transport coefficients and Galilean invariance for various wave propagations in the model. The linear stability of the model is also studied. The model is numerically tested for three cases: (a) a vortex moving with a constant velocity on a mesh periodic boundary conditions; (b) Poiseuille flow with an arbitras~y inclined angle with respect to the lattice orientation: and (c) a cylinder &symmetrically placed in a channel. The numerical results of these tests are compared with either analytic solutions or the results obtained by other methods. Satisfactory results are obtained for the numerical simulations.

Bouzidi, MHamed↗

Wind-Tunnel to Full-Scale Mapping of Winds and Loads for Launch-Vehicle Ground Wind Loads

A launch vehicle ground-wind-loads program was conducted at the NASA Langley Transonic Dynamics Tunnel. The objectives were to quantify key aerodynamic and structural characteristics that impact the occurrence of large wind-induced oscillations of a launch vehicle when exposed to ground winds prior to launch. Of particular interest is the dynamic response of a launch vehicle when a von Kármán vortex street forms in the wake of the vehicle resulting in quasiperiodic lift and drag forces. Vehicle response to these quasiperiodic forces can become quite large when the frequency of vortex shedding nears that of a lowly-damped structural mode thereby exciting a resonant response. The study of ground-wind-loads presents unique challenges to quantify significant characteristics of the approaching wind and to relate the wind-tunnel acquired static, dynamic, and gravitational loads to full-scale vehicles. This paper will explain the correlation process between model-scale and full-scale wind characteristics and resulting structural loads. Characterization of the wind at a launch site requires knowledge of the vehicle dynamics and the anemometer performance. Vehicle dynamics dictate the frequency range of interest in characterization of the atmospheric turbulence, and in defining the required data window-length to model dynamic oscillation build-up during wind-gust analysis. Knowledge of the anemometer performance is required to infer peak wind magnitudes since most anemometers cannot directly measure instantaneous peak values. Following proper wind-speed correlation, static and dynamic wind-tunnel loads must be converted to full-scale equivalent values to derive design and operational guidance. In addition to the direct computation of full-scale equivalent loads using relevant scaling parameters, the corrections applied to dynamic loads accounting for differences in structural damping will also be discussed. Furthermore, a phenomenon that manifests in many slender and flexible launch vehicles is a load magnification due to gravity resulting from mass offset under deflection. A method for scaling this contribution from model-scale to full-scale is presented. Finally, a comparison of wind-tunnel derived loads using the described methodologies to launch vehicle measured loads is presented.

Ground wind loads↗

A two-dimensional, iterative solution for the jet flap

A solution is presented for the jet-flapped wing in two dimensions. The main flow is assumed to be inviscid and incompressible. The flow inside the jet is considered irrotational and the upper and lower boundaries between the jet and free stream are assumed to behave as vortex sheets which allow no mixing. The solution is found to be in satisfactory agreement with two dimensional experimental results and other theoretical work for intermediate values of momentum coefficient, but the regions of agreement vary with jet exit angle. At small values of momentum coefficient, the trajectory for the jet, as computed by this method, has more penetration than that of other available data, while at high values of moment coefficient this solution results in less penetration of the jet into the main flow.

Herold, A. C.↗

Semi-Empirical Prediction of Aircraft Low-Speed Aerodynamic Characteristics

This paper lays out a comprehensive methodology for computing a low-speed, high-lift polar, without requiring additional details about the aircraft design beyond what is typically available at the conceptual design stage. Introducing low-order, physics-based aerodynamic analyses allows the methodology to be more applicable to unconventional aircraft concepts than traditional, fully-empirical methods. The methodology uses empirical relationships for flap lift effectiveness, chord extension, drag-coefficient increment and maximum lift coefficient of various types of flap systems as a function of flap deflection, and combines these increments with the characteristics of the unflapped airfoils. Once the aerodynamic characteristics of the flapped sections are known, a vortex-lattice analysis calculates the three-dimensional lift, drag and moment coefficients of the whole aircraft configuration. This paper details the results of two validation cases: a supercritical airfoil model with several types of flaps; and a 12-foot, full-span aircraft model with slats and double-slotted flaps.

Olson, Erik D.↗