Page 1 of 3 The Spitfire Wing – A Mathematical Model
Part 1
Spitfire Wing As Per Factory Drawings
WING PLANFORM
|
|
Figure 1 |
Figure 2 |
Figure 3 |
AIRFOIL
Ahead of the spar, the thick-skinned leading edge of the wing formed a strong and very rigid D-shaped box, which took most of the wing loads. At the time the wing was designed, this D-shaped leading edge was intended to house steam condensers for the evaporative cooling system intended for the PV XII engine.
In 1933 a family of new airfoil shapes were developed by N.A.C.A. based on a systematic study of aerodynamic characteristics with variations in airfoil thickness and mean-line form. These airfoils were ideal for lofting wings with variable thickness to chord ratios as they were easy to derive mathematically and/or graphically. The 2200 series airfoil gave the thickest nose for the D-shaped leading edge. I suspect the Canadian had a lot to do with the wing design, particularly with airfoil selection.
Drawing No. 33708 sheet 8 has tabulated values for the airfoil sections in standard NACA form (Figure 2), unlike the Mustang P 51-D were the airfoil ordinates readily delimit incidence.
WING TWIST
Wing twist, or washout is helical, unlike the Mustang wing which has a linear (lofted) twist.
WING GEOMETRY
Table 1 ‘Wing Geometry' and the accompanying notes completely define the Spitfire wing, resulting in the rear spar ordinates given under ‘Spar Geometry'. However, there is a problem with the rear spar. There is a discontinuity at the dihedral point resulting in a joggle as well as a twist in the spar web. More about this in Part 2. Table 1 is derived from Supermarine Drawing No. 33708 sheet 8 and is summarized in Figure 4.
|
Table 1 |
Figure 4 |
WING CHORD DISTRIBUTION
The wing plan-form of the Spitfire is not elliptical. Figure 5 gives a comparison between a true elliptical chord distribution and that from Table 1.
|
Figure 5 |
A good way to see how the chord distribution deviates from a true ellipse is to normalize the chord and semi-span axes to unity so that the ellipse becomes a circle (Figure 6).
We then convert the chord from Table 1 into polar coordinates (Table 2) and plot the deviation (Figure 7).
|
Figure 6 |


where

|
Table 2 |
Figure 7 |
We won't worry about the type of curve at this stage. The curve in Figure 7 has some noise and the first task is to smoothen it while keeping as close as possible
to the original ordinates. We use regression analysis (curve fitting) to achieve this.
The following formula gives a good fit.

where the coefficients in Table 3.
|
Table 3 |
Figure 8 and Table 4 show smoothened
.

We now derive a smooth chord distribution.
|
Figure 8 |
Table 4 |
Table 5 gives a smoothened chord distribution to three decimal places.
|
Table 5 |
Chord distribution is not elliptical. However, the slopes of the chord distribution curve in Figure 6 should be 0° at centre of aircraft and -90° at the tip.
In order to meet these conditions, the
Curve in Figure 7 should have 0 slope at
= 0 and
=
/2 .
We now manually (using CAD)draw a spline through
points in Table 2 such that the ends of the spline have zero slope. We then add new
points near the ends of the spline and read off the resultant
points (Table 6).
Table 7 and Figure 9 give the corrected
.
|
Table 6 |
Figure 9 |
Table 7 |
The following formula is used for curve fitting

with coefficients in Table 8.
|
Table 8 |
Now we derive a smooth chord distribution.
This is tabulated in Table 9 and plotted in Figure 10.
|
Table 9 |
This is a good fit shown by the residuals (Table 9) of one hundredths of an inch, that's one quarter of a millimeter. I don't think that the Spitfire wing was manufactured within these limits.
|
Figure 10 |
Table 10 gives a smoothened chord distribution to three decimal places.
|
Table 10 |
Brutal power is used to smoothen the curves. The type of curve-fitting equations used here are unimportant at this stage. The main object is to have smooth curves to work with in Part 2 to build mathematical models. Figure 10a shows the residuals of chord distribution and the noise is evident so Part 1 here deals with 'cleaning up' the curves.
|
Figure 10a |
LEADING EDGE CURVE
We now convert the leading edge offsets from Table 1 into polar coordinates (Table 11) and plot the deviation (Figure 11).

where

|
Table 11 |
The following formula is used for curve fitting the bell shape
with coefficients in Table 12
|
Table 12 |
We now derive a smooth chord distribution.
This is plotted on Table 13 and shown in Figure 11.
|
Figure 11 |
Table 13 |
Table 14 gives smoothened L.E. offsets to three decimal places.
|
Table 14 |
WING TWIST
Wing twist (or washout) is helical along the front spar datum-line: 2° at 31" semi-span to -½° at wing-tip and is calculated with the
following formula
This is tabulated in Table 15.
|
Table 15 |
FRONT SPAR TO REAR SPAR ORDINATES
Distance B along the chord is given by the formula
The horizontal distance from the front spar is
The vertical distance from the front spar is given by

Rear spar ordinates are given in Table 16 to 3 decimal places.
|
Table 16 |
AIRFOIL THICKNESS DISTRIBUTION
The airfoil thickness distribution in % chord is given by the formula
where
b = semi-span = 222.5
with coefficients in Table 17.
|
Table 17 |
This is tabulated in Table 18 and plotted in Figure 12.
|
Table 18 |
Figure 12 |
To summarize, Table 19 gives a smoothened wing to three decimals.
|
Table 19 |
Figure 13 shows the curvature comb of the wing planform as per Table 1.
|
Figure 13 |
Figure 14 shows the curvature as per Table 19.
|
Figure 14 |
Although we have a smooth curve in Figure 14, there are still inflections. We will get rid of these when we do mathematical modelling in Part 2.
End of Part 1


































