Design the WINGGRID - How to evaluate an airplane’s performance.

The wing grid changes the performance of an airplane by influencing its aerodynamic efficiency that is its Lift to Drag ratio L/D. With the wing grid the
span efficiency e can be set to over 4.

This opens several
new design options that can be exploited with the winggrid:

- Increased tip stall resistance
- Bigger wing volume
- Smaller span with same L/D
- Smoother ride due to higher wing loading
- Smaller control surfaces due a rectangular lift distribution and increased stall resistance
- Better climb ability
- Increased range
- Increased payload per span
- Smaller wing aspect ratio with same L/D
- Reduced ground effect
- Rectangular lift distribution leads to easy design due to near two dimensional flow over wing
- Reduced vortex energy increasing airport capacity

Step 1: performance  improvement analysis

For starting we want to find the optimum design of a winggrid in a given case, using the special features of a wing with winggrid:

The tool proposed is the L/D vs speed polar of the aircraft under study. Span efficiency e, span b, wetted surface Swet, reference surface Sref and average Cfavg (friction and interference) are the parameter frame entering the generalised equation L/D vs speed polar. Below we show a few examples of such a study with L/D polar

Required inputs for this step are:

Symbol

Definition

M

aircraft weigth

b

total wing span

Sref

total projected wing area (span * avg_chord)

Sfus_tail

wetted surface fuselage & stabilazers

Cf_avg

area weighted surface friction coefficient including interference drag

q

reference total head (or reference IAS speed)

e

spanefficiency for design baseline without winggrid

Result of this performance prediction is a list of alternatives for

Span

b

Span-efficiency

e

Reference area

Sref

Cl_cutoff

maximum Cl_max for limit of Cutoff up to which the winggrid will work reducing induced drag

plus the respective L/D polars, which are the basis for subsequent winggrid designs.

Step 2: winggrid design procedure
The span-efficiency asked for is satisfied using suitable combinations of:

The diagram below gives the span-efficiency as function of these two parameters, gven by the expression reflecting the airmass deflected with rectangular spanload, relative span winggrid L2/L, N blades:

e = (1 + L2/L * (N - 1)) * 4 / p

In general it is advisable to have 2 < Nb < 5 for reasons of structural simplicity and control of additional interference drag and friction on the blades with smaller chord than main wing.

With L2/L and Nb a further choice has to be made on the type of configuration, e.g.:

For the choosen configuration the stagger angle(s) have to be calculated for the defined Cl-cutoff value, see equation below

DA* = (alfa –alfa0) / (1 – HS*) * HS*; (alfa-alfa0) derived from Cl_cutoff/(2*p)

Checking this expression with the  tests gives for HS*, the critical value of the coupling parameter

0.6 < HS* <0.65

In order to calculate DA* we use the Cl-cutoff value for the wing with winggrid assuming rectangular spanload (note that DA* is the angle between line of grid and zero-lift direction of blade(s).
When applying the stagger calculus for single blades as in the equal angle configuration, blades are identified as part of the grid of blades and therefore again the respective Cl_cutoff values concerne average grid values (in contrast to the local values Cl-cutoff_local=Cl_cutoff/overlap).
Load distribution calculation needs definition of the overlap to be used. As overlap only influences this distribution (and not stagger angle necessary as assumed earlier), smaller overlaps giving flatter load distribution a value in the range 0.7 < overlap < 0.9 is advised, too small overlap results in high local Cl-values on the blades.
Different loads on the grid-elements result in different stagger angles for Cl_cutoff, for equal angle configuration see the typical "banana – look", whereas for equal blade load the mean stagger angle is constant for all blades, although in general the individual blade angles differ somewhat.

Step 3: blade profiles and Betz-correction
Based on the load distribution the blade profiles are optimized and the chord of the winggrid (LE first blade to TE last blade) is adjusted to obtain identical lift gradient for main profile and winggrid.
From the Cl-value for the grid-element of the blade considered the local Cl for the blade follows:

Cl-local = Cl_grid / overlap

This local value is the basis for optimizing the blade profile (Re,range Cl, CD).
We take note here, that the profile choosen for the blade in question should always have identical zero-lift orientation as the main wings profile section.
The aerodynamic loads of blades in a well designed grid reach higher values before separation (corner stall as with turbine grids) compared to the values at separation for single wings, typical values used for design would be: Clmax_blade < 2.5, Clmax_wing < 1.5.  A grid of blades does show (at least for deviation angles less than 0.4) nearly constant flow deviation for a quite important range of overlap, e.g. 1.0 > overlap > 0.6 is corresponding to a range of the factor k of 0.96 to 0.86 (Betz). In the final adjustement of winggrid total chord, this factor is compensated by magnifying the winggrid cross-section by 1/k. This stretching assures identical lift gradient of the winggrid and the main wing profile.

2.0 Stagger angle, high speed and low speed (approach to landing) performance
As borne out by our tests, in order to operate as induced drag reduction, the WINGGRID has to have a minimum stagger angle relative to chord of main wing it is attached to. This stagger angle should as a rule be twice the angle of attack at the lowest speed (speed-limit) for full reduction of induced drag. This speed-limit separates two distinct regimes of speed:

speed effect of the WINGGRID
above speed-limit full induced drag reduction
below speed-limit - fade out of drag reduction from 100% to 0% at 70% of speed-limit but fully resistent to wingtip stall because operating as multiple slit wing piece.
- for landing approach this will help by increasing sinking speed and improving wingtip stall resistance without any controls to be actuated

3.0 Lift distribution and angle of attack adjustements of winglets in WINGGRID
There are basically three different designs possible regarding lift distribution (for details pls ask):

DESIGN EFFECTS
all winglets individually adjusted for equal lift has to be adjusted over the whole speed range
all winglets set at equal angle of attack and fix operates over whole speed range without adjustement however first winglets have much higher lift load
only first winglet is adjusted this will very much relieve the higher lift load on the first winglets

4.0 Design for minimum interference drag of winggrid
A winggrid does cause additional drag, which if not controlled will offset partially the gain made on reduction of induced drag, especially so at the high speed end of an aircrafts L/D polar. In this contribution a few propositions are listed to aid a designer, how to acually control or even nullify parts of this additional drag components.
Background:
Quantitative confirmation of this additional drag was actually the result of the high precision L/D tests of idaflieg 99 made on the first full scale testbed with winggrid. In the tested designs configuration no consideration was given to optimize this drag components and around maximum L/D of the actual polar its amount of around 20% of Cd0 was also verified by detailed drag calculations.

Drag analysis of winggrid
The winggrid exhibits the following additional drag components excluding induced drag compared to drag of the wing considered without winggrid:

The next paragraph is meant to enumerate and explain a few current methods for this.

Minimizing interference drag for lifting profiles ending in a plate
Interference drag for such cases is caused by development of a horseshoe-vortex at the basis of the profile. This vortex carrying away kinetic energy causes the interference drag. Its size is directly a consequence of the boundary layer thickness on the plate near the leading edge of the profile. So we have two methods dealing with this drag components:
Minimizing the boundary layer thickness at profile leading edge on the plate.

This is conveniently achieved by
a)       having a bump  in the plate before the leading edge of the profile, see for information also the bump on modern transport aircraft housing the wheels at the   
          base of the wings.
b)       Placing a NACA ventilation entry port in front of the leading edge of the profile, in order to reset the thickness of the boundary layer to near zero at the    
          profile leading edge junction with the plate
c)       Having the plate length extending in front of the leading edge of the profile as short as possible (applies to the connection of main wing with inner endplate    
          winggrid also)

Avoiding the horseshoe-vortex by adding to the profile a suitable forebody, which helps the streamlines within the boundary layer upstream to surmount the difference in total pressure to the stagnation line on the profile outside the boundary layer on the wall. As the available literature on this method shows, such forebodies inside the boundary layer thickness have to be carefully adapted to the operational characteristics of the profile and the wall boundary layer however. Control of additional friction drag of endplates. To achieve this the main considerations are two:
- Use minimal area to cover blades and main profile sections
- Use careful rounding design at entry leading edge in order to achieve reasonable yaw angle tolerance without partial separations developping.

Summary
Additional drag components of the winggrid consist essentially of three contributions.

As a conservative estimate of the sum of these additional drags for design purpose we use at present a value around 8 % of Cd0 of the wing without winggrid having identical span and chord.

Back to home