XFLR Airfoil Analysis
XFLR Airfoil Analysis
In a future issue...
and a contest announcement 26
Look forward to a Mark Nankivil walkaround of the
3 RC Soaring Digest Editorial Morelli M-200, and an announcement concerning the
Australia F3B Open International, 7-8 March 2008.
4 A Mid-summer trip to Volksrust
Model Analysis with XFLR5 27
A four-day excursion to Tamatie-Berg for some slope
flying, related by Piet Rheeders. XFLR5, a “virtual wind tunnel” allows a detailed
investigation of the aerodynamics of a specific
20 A photo by David Copple airframe. Francesco Meschia demonstrates all
Two Tangent ASH-26 models play in the sun. of the procedures necessary to evaluate the
X-Models 1.9 m Blade V-tail sloper.
21 Have Sailplane - Will Travel
Sloping at Armadillo World Headquarters
Tom Nagel (without a sailplane this time) finds out that
there actually is a slope site in Austin Texas - Mt. Bonnell.
24 At Last... The Troodons are finished! Back cover: Piet Rheeders' photo of the incoming
Not a small bird-like dinosaur, but a 3.65 m span F3J/F5J thunderstorm which brought an end to their flying at
sailplane. Simon Nelson gives a quick overview of his “Tamatie-Berg” Volksrust, South Africa, back in December.
most recent creation. (See Piet's article on this adventure starting on page 4.)
Panasonic DMC-LZ5, ISO 100, 1/250 ,sec, f`5.6, 37mm
February 2008 3
A Mid-summer trip to
Volksrust
Piet Rheeders, [email protected]
BERG club, South Africa
Mike May lets loose his Jart with a mighty heave. Note the whiplash on TX aerial.
February 2008 5
February 2008 5
6 R/C Soaring Digest
A selection of models that we took along.
February 2008 7
8 R/C Soaring Digest
Mike, Glen and myself
(from Durban) set out
in the early morning on
Wednesday the 2nd of
January 2008, and arrived
at Volksrust at around
10 AM. We had to wait,
however until 1 PM before
we could get to the top
of “Tamatie-Berg” due
to a locked farm gate.
This was no problem as
there was no wind from
any direction. This gave
us the time to go to our
overnight guest farm,
situated directly under the
northwesterly slope, to
unpack our baggage and
personal belongings.
The 7.5 meter span, 15 Kg, stand-off scale DG500 Elan of Mike May.
February 2008 9
Glen about to launch his weasel as high as possible in the week lift we had on Day 1.
February 2008 11
On Day 2 the morning broke with a partly This time round we had no problem As we arrived at the SW slope we found
cloudy sky and the SMS I received from with the locked farm gates. Once on top some more glider pilots joining us —
my friend Evan Shaw predicted SE wind we quickly assembled our gliders and Herman and Izak from the ETB club, and
later on in the day. then waited as the early morning mist Charl and Peter from the club BERG.
After breakfast, and just in case the wind disappeared and as the sun got higher in
did not work, I had two flights with my the sky.
IC-power glider before we departed for
the south slope.
February 2008 13
Glen’s pink Jart in a blue sky makes for a pretty sight.
February 2008 15
Pre-flight check on Mike May’s new monster 7.5 meter DG500 Elan.
Day 3 (Friday) we had similar conditions ships and managing some reasonable This happened to be so on Day 4
as Day 1 (Wednesday) and we ended flights. The indications, however, were (Saturday), with the wind direction from
up with light winds turning the full 360 that the next day (Saturday) was going to the northwest and picking up ever so
degrees before ending the day on the be the best day of our visit. gradually until it got strong enough for
northwest slope flying our light thermal Mike to maiden his big DG500.
February 2008 17
Mike had one flight with his 3M Swift
and was confidant that the air was now
good enough to fly the 15 Kg monster
DG 500.
The launch went well with Glen on the
left tip, Charl in the middle and Peter on
the right tip. Mike was quick to pick up
the left wing that dropped a bit because
of a late release by Glen.
The DG500 then climbed high above
the slope and soared for eight minutes
before Mike became aware of a
threatening thunderstorm on the lee-
side of the hill and decided to land
before the wind changed direction.
The landing was also good as the
powerful flaps slowed the DG500 down
and Mike had to tuck them away for a
moment to retain some speed before
touchdown.
The thunderstorm came closer and
we had to dismantle and pack up
our models in a hurry. [A photo of the
approaching thunderstorm serves as
the back cover for this issue.]
I don’t think that any one of us would
regret this outing, and as for myself I
have flown six of seven models that I
took with me, and so as we drove down
this magic mountain we are already The Berg team just before we departed for home on Sunday morning. From left to
planning the next trip to Volksrust and right: Piet, Glen, Jenny, Charl, Peter, Mike. Middle center: Blake
“Tamatie-Berg.” n
February 2008 21
the waterfront, built on a series of man-made lake front
parapets. You do not want to send a plane down there.
This would be a good place to slope something with an
electric motor.
The landing zone comes in two varieties: full of tourists
and full of trees. Luckily, the trees are of the Texas scrub
vegetation variety, and not particularly thorny. Landing will
require either dumping it in the trees on the downwind side
of the access trail, or enlisting a spotter to keep the tourists
at bay while you slide in on the gravel and stone access
path.
This would be a good place to slope something foamiferous
with an electric motor.
The very south end of the park affords a great view of
downtown Austin and the Gozer Building. Austin has only
one really tall building downtown, the proper name of which
February 2008 23
At last...
the Troodons are finished!
Simon Nelson, [email protected]
February 2008 25
In a future issue...
Virtual wind tunnels and simulators long time, and a numerical method for the aeromodeling realm, thanks to the
aircraft performance analysis, known increase in processing power of modern
Many aeromodelers are already familiar
as vortex lattice method, or VLM, was PCs and to a few suitable software
with XFoil, the famous virtual wind tunnel
developed since the 1930’s. I am not packages which are now available. Again
developed by Mark Drela and Harold
an aerodynamicist and so I cannot Mark Drela and its research group have
Youngren at MIT. XFoil, like many of its
describe the VLM in great detail, but released AVL, a powerful and complete
siblings, simulates with good accuracy
from the Internet I learned that VLM program with a “family feeling” with XFoil
the airflow around a “two dimensional”
involves dividing the lifting surfaces into in its command-line user interface. Much
airfoil, i.e. the purely ideal case, not likely
a fine mesh of panels. Each panel is like what happened to XFoil with Stefano
to be attained in practice, of an airfoil
surrounded by a horseshoe vortex, that Duranti’s “Profili 2,” VLM has become
section traveling through the air by itself,
extends chord wise to infinity; with a few more accessible to the non-specialists
without an associated wing, therefore
boundary conditions one can calculate thanks to XFLR5, a program developed
neglecting induced drag and any other
the lift and drag contribution of each by André Deperrois. Mr. Deperrois’s
planform induced effect.
vortex and so, by summing the individual work includes an interesting contribution:
This kind of analysis lends itself well to contributions, one eventually evaluates whereas the “classical” VLM analysis
the comparison of different airfoils but, the performance of the whole surface. assumes a purely inviscid flow around
on the other hand, requires additional the lifting bodies and is therefore a bit
Because of its numerical nature, the
work if one wishes to analyze the unrealistic for the Reynolds numbers
VLM didn’t really “take off” until enough
performance of a wing or a complete used by model aircraft, XFLR5
number crunching power became
aircraft, full-size or model, a task that postulates that the viscous and inviscid
available, with the advent of computers
involves many other variables and contributions to aerodynamic forces are
in the 1960’s; it proved since then to be
requires careful modeling of many other linearly independent, so that an inviscid
a very powerful tool and it was used to
factors. VLM output may be complemented
study and develop a large number of
Because of its extreme importance different aerodynamic configurations. by a viscous XFoil analysis to get a
for the whole aviation industry, the more realistic mathematical model.
In the last few years VLM has landed in Deperrois warns the users, though, that
latter problem was investigated for a
February 2008 27
the “independence hypothesis” is not
supported by a theoretical model, and
so XFLR5 results need to be considered
preliminary and experimental work and
model validation still needs to be done.
XFLR5 users need to keep this
philosophy in mind when approaching
the program. XFLR5 must be “fed”
with both a geometrical model of the
lifting surfaces to be analyzed and
a set of polars derived from viscous
analysis of the adopted airfoils, for
a range of Reynolds number and lift
coefficient broad enough to cover all
flying conditions. In other words, before
even starting XFLR5 the user must find
the coordinate files for all the airfoils
which will be used, and a 3-view of the
aircraft. Also we must remember that
XFLR5 is merely a simulator for a set
of physical behaviors, and can’t tell us
anything about an aircraft but what we
are prepared to ask; it’s a handy tool, but
teaching us how to design a successful
model is entirely out of its scope.
Much better, in my opinion, is to
approach XFLR5 with a practical touch,
to show what the program is capable of
and also what to ask the program for.
I have therefore chosen a well known
model, the 1.9m Blade from X-Models,
and I will use it as a test case to look for
X-Models 1.9m Blade answers to some common questions
such as lift distribution, stability analysis,
Photo by Stefano Bisio
and speed polar determination.
February 2008 29
Figure 4
Figure 6
and chord, and variable lift coefficient); routine will start working and in a little any plot we’re interested in (Cl vs. Cd, Cl
then we enter the desired Reynolds while will crank out the polars. When vs. alpha, Cm vs. alpha, and so on) by
range (I’d suggest to use a list of Re’s, the computation is over we may close choosing the appropriate item in “Polars
like the one shown in Figure 5) and the the window, then we open the “View” > View.”
angle of attack range (e.g. AoA ranging menu and choose “Polars.” The newly- The same procedure must be re-iterated
from -3° to +9°, in 0.25° steps). When computed polar family will be shown for the NACA 0007 tail section (Figure 7),
we click the “Analyze” button the XFoil (Figure 6), and we’ll be able to choose at first for the same range of Reynolds
February 2008 31
Figure 7
numbers used for the main wing section Now we have the 2D polar families for its planform and decalage. For this I will
and for AoA between -4 and +4 degrees; both wing and tail sections, it’s time to use actual geometric data I measured on
theory suggests that we should shift enter the realm of 3D simulation with the my own Blade.
the Reynolds range towards lower Re’s, vortex lattice method. First we must set We can enter the XFLR5 plane geometry
because the tail chord is shorter than the up in XFLR5 a simple three-dimensional routine first by selecting the “Wing
wing chord, but we’ll see later if this will model of the Blade, re-creating the wing Design” option from the “Application”
actually be necessary. planform, the wing dihedral, the tail with menu, then by choosing “Define a Plane”
February 2008 33
fit. Each panel is defined
by its airfoil, its root and
tip chords, its span,
its offset, its dihedral
relative to the horizontal
plane, and its twist;
XFLR5 uses metric units
by default, but it may
be set to use the U.S.
system via the “Units”
option in the “View”
menu. (I am from Europe
so I will use the metric
international system in
data and calculations.)
As the data are entered,
the program immediately
calculates the wing
surface, the mean
aerodynamic chord, and
the wing aspect ratio.
The wing of the Blade
has three degrees
central dihedral, that is
1.5 degrees per each
half-wing relative to the
horizontal plane. The
V-tail is modeled in
XFLR5 as an elevator
without fin: each half-
elevator has 35 degrees
dihedral to the horizontal
plane. I measured a
longitudinal dihedral
Figure 9 (decalage) of one
February 2008 35
of the whole model. The polar computation requires a set
of boundary conditions, and so we must decide whether
we want to either have constant airspeed, or constant lift,
or constant angle of attack. In the gliding task the wing
must generate, at any given moment, just enough lift to
balance the weight of the model, therefore we are mostly
interested in the second type of boundary condition. We
can tell XFLR5 about our choice by choosing “Define a
Polar Analysis” from the “Polars” menu, and by choosing
a “Type 2” analysis from the window that will pop up
(Figure 11). We must also tell XFLR5 the weight of our
model (my Blade weighs 1450 grams) and the position
of the center of gravity (“moment reference location” in
XFLR5 jargon) relative to the leading edge of the wing
root that is taken as the coordinates’ origin (79 mm for
my Blade). From the same window we notice that VLM is
the only applicable analysis method (the other one, the
lifting line theory method or LLT, can only be used for
an isolated wing, i.e. without tail) and we need to make
sure that the “Viscous” option is checked. That said,
we can dismiss the window with the OK button and we
are ready to run the computation task. In the rightmost
window pane we set up a sequence of angles of attack,
for instance from -2°to 7° in steps of 0.25°, we check the
“Store OpPoints” and “Store points outside the polar
mesh” boxes (see Figure 12) and we can start the job
by clicking the “Analyze” button. A window will pop up
with a processing log and maybe some final messages
warning us that some operating point fell outside the
Figure 11 polar mesh. We’ll see later what those log messages
mean, for now we want to jump directly to the results.
We choose (if not already done) the “3D” option in the
“View” menu and we are presented with a beautiful color
rendering of the distribution of the pressure coefficient
over the lifting surfaces (Figure 13). By altering the angle
of attack via the drop-down menu in the upper right part
of the window we can see how the pressure distribution
Figure 13
changes its shape accordingly. and maybe easier to understand if a list of variables we can assign to X and
Undoubtedly impressive, but somewhat little less appealing. For instance, we Y axes. If we assign vertical velocity (Vz)
hard to decipher because on the same may choose “Polars” from the “View” to Y axis and airspeed (Vinf) to X axis
plot we can find lift vectors, drag vectors menu, then we right-click on the graph we’ll get the speed polar for our 1.9m
and also downwash vectors. A number that will appear and last we choose Blade (Figure 14). This graph pictures
of different output modes for the same “Graph” and then “Variables” from the the relationship between airspeed and
data are available, equally interesting pop-up menu. We’ll be presented with a sinking speed, which is a true snapshot
February 2008 37
Figure 14
of the performance of a glider, be it airspeed of about 10 m/s), and another the pilot who wants to stay aloft in dead
full-size or a model. It is worth noticing point, different from the former, for which calm air for as long as possible must trim
that the plot has a minimum sink point the horizontal-speed-to-vertical-speed the sailplane to keep the minimum sink
which corresponds to an airspeed value ratio (also known as efficiency and glide airspeed, whereas the pilot who wants
slightly above stall speed (theoretically ratio) has a maximum (at 11 m/s airspeed to travel as far as possible must trim for
speaking, the 1.9m Blade could attain a for our Blade). These two points picture maximum glide ratio.
sinking speed as low as 40 cm/s at an two very important attitudes of a glider: Another, clearer, way to study the glide
ratio versus airspeed relationship is to which is theoretically above 25 at 11 trim speed is in normal conditions. The
ask XFLR5 to assign the “Glide Ratio Cl/ m/s airspeed. The peak is well defined, answer is readily found by asking XFLR5
Cd” to the Y axis while keeping airspeed that means that any drift from optimum to plot a pitching moment coefficient
(Vinf) on the X axis. The resulting chart airspeed translates into a marked loss of versus airspeed chart (Figure 16). We
(Figure 15) immediately displays not efficiency. see that the curve crosses the horizontal
only the maximum glide ration point, We should now ask ourselves how fast axis at an airspeed of about 18 m/s. This
but also the numerical value of the ratio, the Blade wants to fly, that is what its point is the only balance point for the
February 2008 39
stability later).
It is rather conspicuous that 18 m/s is
significantly faster than both 10 m/s
(minimum sink speed) and 11 m/s
(maximum glide ratio speed). This is
consistent with the kind of flight the
Blade is made for. It’s a fast glider
and not a thermal duration model, a
runner and not a floater. It wants to
“sting like a bee” rather than “float
like a butterfly.” It does not make
much sense to force it to fly at the
minimum sink airspeed, and its pilot
is probably willing to trade some
efficiency for more thrill. At 18 m/s
(40 mph) the Blade sinks in the air
at a rate of a 1.1 meter per second.
That means that a slope updraft with
a vertical velocity component of just
above 1 m/s (fairly common for most
slopes) will be enough to keep it
flying.
The lift distribution
We will now talk about another
way to look at the results of XFLR5
computations. So far we have
Figure 16 looked at charts that are a picture of
the behavior of the whole model as a
function of one of the flight variables.
whole flight envelope of our glider, since also suggests that any higher speed Each point in the curve is an “operating
here the model does not experience would cause the Blade to pitch up and point” of the model, that is a collection of
any pitching moment, neither upward decelerate, and any lower speed would data for lift, speed, efficiency, moment,
nor downward. The model simply make it pitch down and accelerate. That etc., computed for one particular value
flies undisturbed at 18 m/s: this is its means that we have chosen a center of of the independent variable we have
trim speed with neutral elevator trim gravity location that allows for a stable chosen for our analysis (i.e., in our
and 1 degree of decalage. The curve equilibrium (but we’ll see more about examples, the angle of attack). But
February 2008 41
the wing root) ought to generate a
lifting coefficient of 0.86 at Reynolds
23308. According to the XFoil polars
in its memory, XFLR5 doesn’t think
that our modified RG 15 airfoil can
accommodate such requirements.
In other words, XFLR5 foresees the
beginning of a tip stall at this angle
of attack. XFLR5 may not always
be right (in this particular case,
if we examine the family of airfoil
polars we can see there’s such a
wide a gap between the Re=20000
and Re=30000 polars that XFLR5
interpolation might not be able to
bridge reliably), but it is an interesting
clue nevertheless. It is also worth
mentioning that, despite XFLR5’s
opinion, the “real” Blade is not
particularly prone to tip stall problems
(or, at least, it hasn’t yet played
such bad tricks to me). This might
support our initial suggestion that, at
the wing tips, the RG 15 foil evolves
into a different section, one that can
generate more lift at low Reynolds
numbers.
Figure 18 Another variable we might want to
examine the distribution of is local
“View Log File”); at the end of the file we Cl = 0.86 could not be interpolated lift (not lift coefficient). This is the force
should see a similar message: Span pos = 937.49 mm, Re = 23 308, that is generated by every wing portion;
Cl = 0.86 could not be interpolated it is important because, according
...Alpha=5.75
Calculating elevator... to Prandtl’s theory, the induced drag
Calculating induced angles... contribution for a given wing area is
Calculating aerodynamic coefficients... What this message is telling us is that, minimized whenever the lift distribution
Calculating wing... at an AoA of 5.75°, the wing tip station shape is elliptical. We can ask XFLR5 for
Span pos = -937.49 mm, Re = 23 308, (937 mm, or about three feet, away from this distribution by right-clicking on the
February 2008 43
the CG at 100 mm from the leading
edge. The resulting CM vs. AoA curve
will be like the green curve in Figure
20: it is a monotonically increasing
curve, and therefore the model will
be unstable for any AoA in the flight
envelope.
We can take this even further by
tracing other Cm vs. AoA curves for
other CG locations. For instance,
in Figure 20 I have also plotted (in
blue) the curve that corresponds to
a Blade balanced at 95 mm from the
leading edge. This one is no longer
monotonic but has a maximum at
about two degrees. Balanced in this
way, the model is neutrally stable.
Strictly speaking, it will be neutral
only at the angle of attack for which
the dCM/dAoA derivative is null, but
in practice the CM shows such a
slight variation, and for such a broad
range of AoA’s, that the stabilizing
or de-stabilizing effect may be
negligible. The model is not willing to
react to any disturbance, but will fly
at and maintain any angle of attack
Figure 20 we can trim it to. As a consequence, it
will require a pilot’s constant attention
moment coefficient versus AoA curve The steeper the slope of the CM vs. AoA and input, and personally I would not
is decreasing from left to right is stable; curve, the greater the stabilizing or de- enjoy flying it at a slope with a brisk
on the other hand all configurations with stabilizing action will be. breeze.
an increasing CM vs. AoA curve are If we want to decrease a model stability, The use of camber flaps and
unstable, because every disturbance or even make it unstable, we just need variable camber
results in a pitch unbalance that tries to to move the CG aft; we can test that The Blade, like many other gliders, has
further increase the disturbance effect. by computing a new polar and placing control surfaces (inboard flaps and
February 2008 45
versus AoA plots (View > Polars, then
right-click > Graph > Variables) we
will probably see a tangled mess of
curves, but if we clean up things a bit
(Polars > Hide all polars) and then we
plot only a single Reynolds number
for each airfoil modification (by
choosing the airfoil and the desired
Re from the drop-down lists in the
toolbar, then checking the “Show
Curve” box) we should see a graph
like the one shown in Figure 22. This
graph shows that, in a limited AoA
range, different flap settings result in
a conspicuous variation in Cl, with a
limited influence on Cd.
We might now feel the temptation
to see what flap setting allows for
maximum Cl/Cd ratio, but we must
always keep in mind that these
are the 2D polars, and the results
may not be immediately translated
into performance prediction for the
complete model. My suggestion
is to spend some more time and
“persuade” XFLR5 to simulate the
complete model with differently
Figure 23 cambered airfoil sections, so that our
investigations will be based upon a
wing flaps) and of the airfoil thickness the polar family exactly like we did for 3D model of the plane behavior.
(100% in our case, because all wing the unmodified airfoil, and we might To set up a 3D model of the Blade
control surfaces in the Blade are top also want to repeat the whole process with flapped airfoil, first we shall enter
hinged); then we click the “OK” button for a broader range of flap deflection the 3D module of the program (via the
and we store the modified profile under angles (in addition to 0 and +2 degrees, “Application > Wing Design” menu).
a suitable name (for instance, “RG-15 I suggest +4, -2 and -4 degrees). Now, Then we recall the already-processed
7.8% +2”). Now we must pre-process if we look at the Cl versus AoA and Cd Blade model via the drop-down list
February 2008 47
tendency to oscillate around the
desired airspeed rather than keep it
throughout the flight.
This behavior can be investigated via
XFLR5 by creating a modified NACA
0007 airfoil with a leading edge flap at
75% chord, deflected by -1.5 degrees
(upward deflection, that is), and
instructing the program to simulate
a Blade endorsing this airfoil section
in its V-tail. We will get a CM vs. Vinf
graph like the one shown in Figure 25.
The model with “up trim” (magenta
curve) is now balanced at an airspeed
of less than 10 m/s, but the slope
of the curve near the balance point
is much steeper than the slope of
the curve that describes the “zero
trim” Blade (blue curve) at the 15
m/s balance point. This means that,
for instance, any drift from the trim
airspeed for the “trimmed up” Blade
will result in a pitch response three
times stronger than we would have
with the neutrally trimmed model;
therefore, we may expect a stronger
tendency to pitch oscillations.
Figure 25
Obviously there is a way to get the
model to balance at 10 m/s without
graph, Figure 24, the curve shows that as slow as required. But there’s a catch: requiring so much up trim; we just need
the zero trim airspeed is about 15 m/s, as we trim the model in this way, the to move the CG aft. If we re-compute the
slower than the trim speed for when pitching moment to airspeed relation model polar by setting the CG at 88 mm
no flaps were used, but faster than our will change, and in particular when from the leading edge, we see that CM
target speed for which more camber is we trim for “nose up” the model will goes to zero at the desired airspeed of
useful. This means that we should also become more pitch-sensitive to airspeed 10 m/s (Figure 26), and that the slope of
dial some up trim to make the Blade fly changes. This will result in a stronger the curve is not too step, as predicted.
February 2008 49
airspeed, especially at a turbulent
slope where it would be virtually
impossible to tell the former from
the latter. But, just for fiction, let’s
suppose we want to bring our Blade
to a thermal duration contest on a
flat field. In this case it’s easier to
distinguish the two flight conditions,
each in its own particular “task.”
We will seek the best glide ratio
when we want to explore as much
air mass as possible while hunting
for a thermal, and we will struggle to
fly at the minimum sink rate in dead
air, where there are no thermals to
be found, or while we are trying to
gain as much altitude as possible
in a weak thermal. In addition, if we
need to penetrate headwind or to
escape from a sinking air zone, it is
also important to fly faster than the
best glide ratio airspeed. In other
words, it becomes important for a
pilot to exploit the “extended polar”
of its model, that is the curve arising
from the envelope of the many
different polar curves associated with
Figure 27 variable camber. To be successful
in this, we need to study the curves
and match the theoretical capabilities
Trade-offs and conclusions need to deflect the leading edge flaps
of the model with our piloting skills. For
by 4 degrees downwards and keep an
By looking at Figure 26 we see that example, from the charts we see there
airspeed of 9 m/s (20 mph). We already
the 1.9-m Blade has its peak efficiency is hardly any gain, either in sink rate
mentioned that the Blade is a slope
(best glide ratio) when flying at about or in glide ratio, when we deflect the
model, designed to fly fast and that
11 m/s (25 mph) with “flat” wing (no leading edge upwards by more than two
probably nobody would want to fly it at
flaps deflected), and that when we want degrees; this may be translated in a radio
the minimum sink or at the best glide
to descend as slowly as possible we setup that will allow for a “speed” flight
February 2008 51