A simple solution for unswept three-dimensional wings can be obtained by using Prandtl's lifting line model. For incompressible, inviscid flow, the wing is modelled as a single bound vortex line located at the 1/4 chord position and an associated shed vortex sheet.
where Γ is the circulation (vortex strength), defined a function across the wing span, V∞ is free-stream velocity, b is wing span, s is wing semi-span, c is wing chord and y is the distance across the span measured from the wing root.
The mapping of angle (θ) to semi-span (s) position is done using Fourier series and allows variation of the model to suit different geometries.
The span-wise lift distribution is assumed to be approximately elliptical with a small modifications due to wing plan-form geometry. The vortex line strength can thus be modelled using the following Fourier series approximation,
The required strength of the distribution coefficients (An) for a given geometry and set of free-stream conditions can be calculated by applying a surface flow boundary condition. The equation used is based on the usual condition of zero flow normal to the surface.
For 3-D wings this condition is applied at several span-wise locations. By matching flow and surface angles, the normal velocity boundary condition can be restated as a requirement for the flow angles at the section to be in balance. Unlike 2-D section flow where flow angles are set by freestream direction and surface angles only, in 3-D wing flow an extra flow angle component is introduced by the shed vortices that are produced and trail behind the wing. This trailing vortex sheet produces a downwash.
For a 3-D wing, the local section flow angle must be equal to the sum of the wing's angle of attack, the section twist and the downwash induced flow angle.
where
and where α is the 3-D wing angle of attack, θt is the wing twist angle and wi is the velocity induced by trailing vortex sheet.
The downwash velocity is caused by the shed vorticies trailing behind the wing. The vortex strength in the trailing sheet will be a function of the changes in vortex strength along the wing span. The mathematical function describing the vortex sheet strength is thus obtained by differentiating the bound vortex distribution.
Thus the downwash at any span position on the wing can be found by integrating the influence of individual elements of the trailing sheet. Each sheet element behaves like 1/2 of an infinite vortex line so
where r is the distance across the span between the vortex element and the point at which downwash is being calculated. r = ( y – yi). The full downwash will be,
On substitution of the mapping variables and integration, the result is,
The 2-D section lift coefficient is a function of the local flow incidence and the bound vortex strength at this span location. So that,
where a0 is the section lift slope (dCL/dα2D), α0 is the section zero lift angle and c is the chord length.
This can be rearranged in terms of vortex strength,
and substituting for vortex strength and induced angle produces the following boundary condition equation,
This final boundary equation contains all the unknown coefficients of the wing model's vortex distribution, along with the wing's geometry and the stream conditions. It can be used to find coefficients A1,A2, A3, ...
Assuming that higher order coefficients become increasingly small and make negligible contribution to the result, one method of solution is to truncate the series at term AN. By applying the boundary condition at N span locations a set of simultaneous linear equations can be constructed. This set can be solved for coefficients A1 to AN.
A cosine distribution of span-wise locations should be used for the boundary conditions to match the assumed wing loading distribution. Clearly the number of coefficients used will determine the accuracy of the solution. If the wing loading is highly non-elliptical then a larger number of coefficients should be included. This occurs when analysing wings with part span flaps. This type of geometry causes a discontinuity in the spanwise loading and hence requires a much larger number of coefficients to accurately describe the distribution. Where the the wing loading is symmetric about the wing root, the contribution of even functions will become zero. Coefficients A2,A4,A6, ... are all zero and can be dropped from the analysis.
Once the coefficients of the load distribution are known the total lift of the wing can be found by integration.
which reduces to a lift coefficient of
where AR is the aspect ratio of the wing
A consequence of the downwash flow is that the direction of action of each section's lift vector is rotated relative to the free-stream direction. The local lift vectors are rotated backward and hence give rise to a lift induced drag. While the overall governing equations are potential flow and hence do not give rise to friction or pressure drag, this lift induced drag will be a significant component of the overall drag of the wing. The downwash velocity induced at any span location can be calculated once the strength of the wing loading is known. The variation in local flow angles can then be found. By integrating the component of section lift coefficient that acts parallel to the free-stream across the span, the induced drag coefficient can be found.
which produces the following induced drag coefficient ,
No information about pitching moment coefficient can be deduced from lifting line theory since the lift distribution is collapsed to a single line along the 1/4 chord. All forces will be centered at the 1/4 chord in this model irrespective of aerofoil section geometry.
If the wing planform is elliptical,c=c0cos(θ) then it can be shown that the wing load distribution is also a purely elliptical function and all coefficients except A1 will vanish.
Also for simple straight, untwisted wings, this solution may be used as a first approximation. For these geometries all other coefficients will be much less than A1 and the loading will be approximately elliptical. In this case, a single general boundary condition equation results containing only one unknown, the vortex line strength at the wing root.
The solution of the boundary condition in this case leads to a constant downwash across the span and the following simple answers for lift coefficient and induced drag coefficient.
Houghton & Carruthers "Aerodynamics for Engineering Students", Arnold Ed 3 1982.
The following computer program allows the user to define wing plan-forms (without sweep) and to define wing root and wing tip section properties. The program assumes a linear variation of section properties between wing root and tip. For the MS Windows downloadable version, and that the wing loading is assumed to be symmetric about the wing root.
Section lift data in terms of a0, α0, ie. CL(2D) = a0(α2D - α0) are required for the wing root and the wing tip sections. This information can be obtained from published 2-D experimental data or theoretical techniques such as thin aerofoil theory.
The program uses the above lifting line equations to get solutions for lift coefficient versus angle of attack and induced drag coefficient versus lift coefficient2. For a given angle of attack the program will display the resulting distribution of section lift coefficient across the span, the distribution of downwash at the wing and a listing of the solution Fourier coefficients for this angle.
A flapped section can also be input. The percentage of wing span with flap must be input to create a flapped wing section. The flap section properties are assumed to be those entered for the wing root section. These section properties will be kept constant across the flapped portion of the span. The section properties used outboard of the flap will also be constant and assumed to be equal to those of the wing tip.
Prandtl Lifting Line Solver for 3D unswept wings (MS Windows downloadable)