Effect of Finite Span on the Airload Distributions for Oscillating Wings I : Aerodynamic Theory of Oscillating Wings of Finite Span
Explore the source record for details and available documents.
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.
Explore the source record for details and available documents.
The oscillating air forces on a two-dimensional wing oscillating in pitch about the midchord have been measured at various mean angles of attack and at Mach numbers of 0.35 and 0.7. The magnitudes of normal-force and pitching-moment coefficients were much higher at high angles of attack than at low angles of attack for some conditions. Large regions of negative damping in pitch were found, and it was shown that the effect of increasing the Mach number 0.35 to 0.7 was to decrease the initial angle of attack at which negative damping occurred. Measurements of the aerodynamic damping of a 10-percent-thick and of a 3-percent-thick finite-span wing oscillating in the first bending mode indicate no regions of negative damping for this type of motion over the range of variables covered. The damping measured at high angles of attack was generally larger than that at low angles of attack. (author)
Explore the source record for details and available documents.
Determination of the spontaneous oscillations of a wing or tail unit entail many difficulties, both the mathematical determination and the determination by static wing oscillation tests being far from successful and flight tests involving very great risks. The present paper gives a method developed at the Junkers Airplane Company by which the critical velocity with respect to spontaneous oscillations of increasing amplitude can be ascertained in flight tests without undue risks, the oscillation of the surface being obtained in the tests by the application of an external force.
Three wing models were oscillated in yaw about their vertical axes to determine the effects of systematic variations of frequency and amplitude of oscillation on the in-phase and out-of-phase combination lateral stability derivatives resulting from this motion. The tests were made at low speeds for a 60 degree delta wing, a 45 degree swept wing, and an unswept wing; the swept and unswept wings had aspect ratios of 4. The results indicate that large changes in the magnitude of the stability derivatives due to the variation of frequency occur at high angles of attack, particularly for the delta wing. The greatest variations of the derivatives with frequency take place for the lowest frequencies of oscillation; at the higher frequencies, the effects of frequency are smaller and the derivatives become more linear with angle of attack. Effects of amplitude of oscillation on the stability derivatives for delta wings were evident for certain high angles of attack and for the lowest frequencies of oscillation. As the frequency became high, the amplitude effects tended to disappear.
Unsteady flowfields of a two-dimensional oscillating wing are calculated using an implicit, finite-difference, Navier-Stokes numerical scheme using five widely used turbulence models. The objective of this study is to identify an appropriate turbulence model for accurate simulation of three-dimensional dynamic stall. Three unsteady flow conditions corresponding to attached flow, light-stall, and deep-stall of an oscillating wing experiment were chosen as test cases for computations. Results of unsteady airload hysteresis curves, harmonics of unsteady pressures, and instantaneous flow pictures are presented. Comparison of unsteady airloads with experiment show that all models are deficient in some sense and not a single model predicts all airloads consistently and in agreement with experiment for all flow conditions. For the attached flow condition, the Renormalization Group Theory (RNG), the Johnoson-King (J-K), and the Spalart-Allmaras (S-A) models have better performance. The Baldwin-Lomax (B-L) and the Baldwin-Barth (B-B) models fair poorly. At the light-stall condition, the results for the RNG, the J-K, and S-A models are in agreement with experiment for the upstroke but they all over predict the separation shown by the experiment and therefore have bigger hysteresis loops than experimental results. The B-B model results are also in good agreement for upstroke but have poor lift hysteresis for downstroke. It has superior drag and pitching-moment predictions. For deep-stall conditions, the airloads for the RNG, the B -B, and the S-A models have fair agreement with experiment, but the B-B model performed better at the extreme deep-stall condition. Overall, the RNG model provides significant improvement over the B-L model in all flow regimes with no additional computational cost. The Baldwin-Barth model is the most expensive of the models considered here, costing about 2.5 times that of the Baldwin-Lomax model. Finally, a brief discussion of the effects of grid density, time-step size, and numerical dissipation on the unsteady solutions are also presented.
A computer program was developed to account approximately for the effects of finite wing thickness in the transonic potential flow over an oscillating wing of finite span. The program is based on the original sonic-box program for planar wing which was previously extended to include the effects of the swept trailing edge and the thickness of the wing. Account for the nonuniform flow caused by finite thickness is made by application of the local linearization concept. The thickness effect, expressed in terms of the local Mach number, is included in the basic solution to replace the coordinate transformation method used in the earlier work. Calculations were made for a delta wing and a rectangular wing performing plunge and pitch oscillations, and the results were compared with those obtained from other methods. An input quide and a complete listing of the computer code are presented.
A computer program was developed to account approximately for the effects of finite wing thickness in transonic potential flow over an oscillation wing of finite span. The program is based on the original sonic box computer program for planar wing which was extended to account for the effect of wing thickness. Computational efficiency and accuracy were improved and swept trailing edges were accounted for. Account for the nonuniform flow caused by finite thickness was made by application of the local linearization concept with appropriate coordinate transformation. A brief description of each computer routine and the applications of cubic spline and spline surface data fitting techniques used in the program are given, and the method of input was shown in detail. Sample calculations as well as a complete listing of the computer program listing are presented.
The two-dimensional problem of the oscillating wing with aerodynamically balanced elevator is treated in the manner that the wing is replaced by a plate with bends and stages and the airfoil section by a mean line consisting of one or more straights. The computed formulas and tables permit, on these premises, the prediction of the pressure distribution and of the aerodynamic reactions of oscillating elevators and tabs with any position of elevator hinge in respect to elevator leading edge.
A finite difference method for solving the unsteady flow about harmonically oscillating wings is investigated. The procedure is based on separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady differential equation for small disturbances. Solutions are obtained using relaxation procedures. It is determined that there is a limit on reduced frequency, which is a function of Mach number and size of mesh region, above which the relaxation procedures will not converge. It is found that row line relaxation is more efficient than column relaxation and results are presented for a rectangular wing in harmonic pitch.
A computer program has been developed to account approximately for the effects of finite wing thickness in the transonic potential flow over an oscillating wing of finite span. The program is based on the original sonic-box program of Rodemich and Andrew, and accounts for the nonuniform flow caused by finite thickness by application of the local linearization concept. A brief description of each subroutine is given, and the method of input is shown in detail. A sample problem as well as a complete listing of the computer program are presented.
A finite difference method to solve the unsteady transonic flow about harmonically oscillating wings was investigated. The procedure is based on separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady differential equation for small disturbances. The differential equation for the unsteady velocity potential is linear with spatially varying coefficients and with the time variable eliminated by assuming harmonic motion. An alternating direction implicit procedure was investigated, and a pilot program was developed for both two and three dimensional wings. This program provides a relatively efficient relaxation solution without previously encountered solution instability problems. Pressure distributions for two rectangular wings are calculated. Conjugate gradient techniques were developed for the asymmetric, indefinite problem. The conjugate gradient procedure is evaluated for applications to the unsteady transonic problem. Different equations for the alternating direction procedure are derived using a coordinate transformation for swept and tapered wing planforms. Pressure distributions for swept, untaped wings of vanishing thickness are correlated with linear results for sweep angles up to 45 degrees.
A finite difference method for the solution of the transonic flow about a harmonically oscillating wing is presented. The partial differential equation for the unsteady transonic flow was linearized by dividing the flow into separate steady and unsteady perturbation velocity potentials and by assuming small amplitudes of harmonic oscillation. The resulting linear differential equation is of mixed type, being elliptic or hyperbolic whereever the steady flow equation is elliptic or hyperbolic. Central differences were used for all derivatives except at supersonic points where backward differencing was used for the streamwise direction. Detailed formulas and procedures are described in sufficient detail for programming on high speed computers. To test the method, the problem of the oscillating flap on a NACA 64A006 airfoil was programmed. The numerical procedure was found to be stable and convergent even in regions of local supersonic flow with shocks.
An analytical approach is presented to account for some of the nonlinear characteristics of the transonic flow equation for finite thickness wings undergoing harmonic oscillation at sonic flight speed in an inviscid, shock-free fluid. The thickness effect is accounted for in the analysis through use of the steady local Mach number distribution over the wing at its mean position by employing the local linearization concept and a coordinate transformation. Computed results are compared with that of the linearized theory and experiments. Application to a flutter problem is shown.
An analytical approach is presented to account for some of the nonlinear characteristics of the transonic flow equation for finite thickness wings undergoing harmonic oscillation at sonic flight speed in an inviscid, shock-free fluid. The thickness effect is accounted for in the analysis through use of the steady local Mach number distribution over the wing at its mean position by employing the local linearization concept and a coordinate transformation. Computed results are compared with that of the linearized theory and experiments. Based on the local linearization concept, an alternate formulation avoiding the limitations of the coordinate transformation method is presented.
Analytical and empirical studies of a finite difference method for the solution of the transonic flow about an harmonically oscillating wing are presented along with a discussion of the development of a pilot program for three-dimensional flow. In addition, some two- and three-dimensional examples are presented.
The design and usage of a pilot program using a finite difference method for calculating the pressure distributions over harmonically oscillating wings in transonic flow are discussed. The procedure used is based on separating the velocity potential into steady and unsteady parts and linearizing the resulting unsteady differential equation for small disturbances. The steady velocity potential which must be obtained from some other program, is required for input. The unsteady differential equation is linear, complex in form with spatially varying coefficients. Because sinusoidal motion is assumed, time is not a variable. The numerical solution is obtained through a finite difference formulation and a line relaxation solution method.
A method based on the Navier-Stokes equations was developed for determining analytically the three-dimensional unsteady flow patterns around oscillating wings. The Helmholz vorticity transport equations were discretized in three-dimensional finite element form from a variational formulation and integrated numerically. At each time step of the numerical integration the velocity field was calculated from the representation of the three-dimensional wing by a system of optimized distribution of vortices in space. During the numerical integration of the vorticity transport equations the time-dependent boundary conditions on the wing were specified as external constraint conditions. Examples of obtained results describing the three-dimensional unsteady flow around a wing were presented.