Aerodynamics Paper
Aerodynamics Paper
Aerodynamics Paper
2010
LIST OF CONTENTS
C.D. HARLEY 2
LIST OF CONTENTS
C.D. HARLEY 3
LIST OF FIGURES
LIST OF FIGURES
FIGURE 1.1: PHOTO OF AN AIRBUS A380 DURING LANDING PHASE WITH SPOILERS AND
HIGH LIFT DEVICES FULLY DEPLOYED. TAKEN FROM [6] ..................................... 18
FIGURE 2.4: TYPICAL REARWARD MOUNTED (APART FROM WHERE STATED) SPOILER
EFFECT ON LIFT, DRAG, PITCHING MOMENT AND CHANGE IN SURFACE PRESSURE
WITH DEFLECTION. A) AND B) ADAPTED FROM ESDU DATA SHEET [2]. C)
ADAPTED FROM MCLACHLAN ET AL [16] ............................................................ 27
FIGURE 3.1: ONE OF THE FIRST APPLICATIONS OF A FORWARD MOUNTED AND FORWARD
HINGED SPOILER TO AIRCRAFT WINGS TAKEN FROM WEICK AND WENZINGER [21]
............................................................................................................................ 30
FIGURE 3.2: COMPARISON OF ROLLING AND YAWING MOMENT COEFFICIENTS OBTAINED
WITH AILERONS AND SPOILERS DURING WIND TUNNEL TESTING, TAKEN FROM
WEICK ET AL [24]. .............................................................................................. 32
FIGURE 3.3: LATERAL CONTROL DEVICE CONFIGURATIONS ON THE FAIRCHILD 22
AIRPLANE DURING FLIGHT TESTS. TAKEN FROM WEICK ET AL [23]. .................... 33
FIGURE 3.4: TIME HISTORY CURVES SHOWING THE LAG CHARACTERISTICS OF VARIOUS
CONTROL SYSTEMS. INDICATED AIR SPEED, 22M/S, FULL CONTROL DEFLECTION
TAKEN FROM WEICK ET AL [23].......................................................................... 34
FIGURE 3.5: SUMMARY OF A SELECTION OF SPOILER AND DEFLECTOR COMBINATIONS
INVESTIGATED BY WENZIGNER AND ROGALLO. [25]........................................... 35
FIGURE 3.6: THREE QUARTER REAR-VIEW OF THE TEST AIRPLANE AS INSTRUMENTED FOR
FLIGHT SHOWING DEFLECTED FLAP, DROOPED AILERON, DEFLECTED SPOILER AND
OPEN SLOT. TAKEN FROM [26] ............................................................................ 36
FIGURE 3.7: EFFECT OF CHORDWISE LOCATION ON PITCHING MOMENTS FOR FLAPS AND
SPOILERS, TAKEN FROM PURSER AND MCKINNEY [28]. ...................................... 38
C.D. HARLEY 4
LIST OF FIGURES
FIGURE 3.9: PHOTO TAKEN FROM INSIDE A BOEING 737 DURING THE APPROACH PHASE . 40
FIGURE 3.10: SPOILER DEPLOYED DURING A) CRUISE CONDITIONS OF AN AIRBUS A333
AND B) LANDING PHASE OF AN AIRBUS A380. ..................................................... 41
FIGURE 3.11: MICRO DRAG GENERATOR (MDG) SYSTEM SHOWING MDG‟S ON THE UPPER
AND LOWER SURFACE ADAPTED FROM BAUER [37]. ............................................ 43
O
FIGURE 3.12: EFFECT OF SPOILER SWEEP ON THE LATERAL CONTROL POWER OF A 60
SWEPT TAILLESS AIRCRAFT CONFIGURATION. TAKEN FROM DORSET AND MEHL
[36] ..................................................................................................................... 44
FIGURE 3.13: LATERAL CONTROL POWER OF A SPOILER AND SSD APPLIED TO A 60O SWEPT
TAILLESS AIRCRAFT CONFIGURATION. TAKEN FROM DORSET AND MEHL [36] .... 45
FIGURE 3.14: NORMAL BLOWING FLUIDIC SPOILER CONCEPT TAKEN FROM LEOPOLD ET
AL [40]. ............................................................................................................... 46
FIGURE 3.16: DIAGRAM OF ,A) CROSS-SECTIONAL VIEW OF THE MODEL AT THE SLOT MID-
SPAN AND B) TOP VIEW OF THE SEMI-SPAN WIND TUNNEL MODEL USED BY
TAVELLA ET AL [42]. .......................................................................................... 49
FIGURE 4.1: OVERVIEW OF THE PROJECT TUNNEL TEST SECTION SHOWING MAJOR
DIMENSIONS AND MODEL ORIENTATION. ............................................................. 53
FIGURE 4.2: CROSS-SECTIONAL VIEWS OF THE A) MIGS, B) CFFS WIND TUNNEL MODEL
CONFIGURATIONS, AND C) WIND TUNNEL MODEL DETAILS. ................................. 54
FIGURE 4.3: WIND TUNNEL MODEL INTERNAL STEEL „SKELETON‟ LAYOUT ..................... 55
FIGURE 4.4: PHOTO OF THE WIND TUNNEL MODEL IN CFFS CONFIGURATION ................. 56
FIGURE 4.5: CFFS DESIGN USED IN THE WIND TUNNEL MODEL ....................................... 57
FIGURE 4.6: PHOTO OF THE WIND TUNNEL SETUP WITH VERTICALLY MOUNTED MODEL, 6-
COMPONENT FORCE BALANCE, SCANIVALVES FOR SURFACE PRESSURE AND WAKE
RAKE FOR DRAG MEASUREMENT. ........................................................................ 59
FIGURE 4.7: LAYOUT OF THE WIND TUNNEL SETUP. SOLID ARROWS SHOW THE DIRECTION
OF AIRFLOW FROM THE COMPRESSED SUPPLY, DASHED LINES INDICATE 0.35C
CFFS LOCATION AND PRESSURE SUPPLY PIPING. ................................................. 60
FIGURE 4.8: PHOTO OF THE WAKE RAKE USED DURING WIND TUNNEL EXPERIMENTS ...... 61
FIGURE 4.9: FLOW CHART OF EXPERIMENTAL PROCEDURE .............................................. 65
FIGURE 4.10: MIGS NORMAL TO A WALL PARTIALLY SUBMERGED BY A BOUNDARY LAYER
OF THICKNESS, 0.99. ............................................................................................ 71
C.D. HARLEY 5
LIST OF FIGURES
COEFFICIENT VALUES SHOWN ARE MEASURED FROM THE FORCE BALANCE AND
MOMENTUM LOSS IN THE WAKE. .......................................................................... 77
FIGURE 4.14: SPANWISE VARIATION OF STATIC PRESSURE IN THE CFFS PLENUM ........... 78
FIGURE 5.1: COMPUTATIONAL GRID USED FOR THE FLUIDIC SPOILER STUDY. CASE SHOWN
IS FOR SPOILER AT 0.65X/C .................................................................................. 85
FIGURE 5.2: EFFECT OF GRID REFINEMENT ON AEROFOIL LIFT COEFFICIENT, FOR UPPER
SURFACE CFFS AT 0.65C, WITH MAX. C, AND = 6 . ........................................ 86
O
FIGURE 5.3: CLOSE-UP OF THE MESH AROUND THE RA16SC AEROFOIL WITH INLAY
SPOILER ............................................................................................................... 88
FIGURE 5.5: CLOSE-UP OF THE MESH AROUND THE NACA 0018 AEROFOIL WITH NORMAL
BLOWING FLUIDIC SPOILER .................................................................................. 90
FIGURE 5.6: FLUIDIC SPOILER VALIDATION CASE FOR A NACA 0018 AIRFOIL WITH JET
O 5
ISSUING FROM THE LOWER SURFACE; = 0 , AND RE = 2X10 , EXPERIMENTAL C
= 0.48, SIMULATED C = 0.24.[40] ...................................................................... 92
FIGURE 6.1: COMPARISON OF THE BASELINE EXPERIMENTAL AND NUMERICAL LIFT, DRAG
AND SURFACE PRESSURE DISTRIBUTIONS. ............................................................ 94
FIGURE 6.2: EFFECT OF MICRO GEOMETRIC SPOILER (MIGS) LOCATION ON THE CHANGE
IN EXPERIMENTAL LIFT, DRAG AND PITCHING MOMENT WITH SPOILER HEIGHT.
ANGLE OF ATTACH RANGE: ○ = 0O, ∆ = 3O, □ = 6O. CORRESPONDING SURFACE
PRESSURE DISTRIBUTION PLOTS (IN FIG. 6.3) INDICATED IN VERTICAL AXIS (ACP,
BCP, CCP, DCP). .................................................................................................... 98
FIGURE 6.3: EFFECT 0.03C MIGS LOCATION ON CHANGE IN SURFACE PRESSURE
OF
DISTRIBUTIONS AT = 6 . ................................................................................... 99
O
FIGURE 6.4: GAIN (EFFICIENCY) OF MIGS FOR BOTH LOWER AND UPPER SURFACES, AND
BOTH CHORDWISE LOCATIONS ............................................................................. 99
FIGURE 6.7: GAIN (EFFICIENCY) OF CFFS FOR BOTH LOWER AND UPPER SURFACES, AND
BOTH CHORDWISE LOCATIONS ........................................................................... 103
FIGURE 6.8: EXPERIMENTAL LIFT-DRAG POLAR PLOTS COMPARED WITH DATA SHEET
(MIGS) AND CFD (CFFS) RESULTS .................................................................. 106
FIGURE 6.9: EXPERIMENTAL PITCHING MOMENT – LIFT POLAR PLOTS........................... 106
FIGURE 6.10: FORCE BALANCE DRAG COEFFICIENT FROM DUAL SURFACE BLOWING FOR
YAW CONTROL. UPPER SURFACE CFFS = 0.35X/C. LOWER SURFACE CFFS =
0.65X/C. ZERO PITCHING MOMENT (CM) AND LIFT (CL) LOCI INDICATED. .......... 109
C.D. HARLEY 6
LIST OF FIGURES
C.D. HARLEY 7
LIST OF TABLES
LIST OF TABLES
TABLE 4.1: WIND TUNNEL BALANCE FORCE AND TORQUE NOMINAL RANGES, IN FORCE
BALANCE AXES. ................................................................................................... 62
TABLE 4.2: UNCERTAINTY CALCULATION OF THE FORCE BALANCE LIFT COEFFICIENT. THE
EXPANDED UNCERTAINTY AT THE 95% CONFIDENCE LEVEL FOR CL = 0.6197 IS
9.15%. ................................................................................................................ 81
TABLE 5.1: OVERVIEW OF THE IMPORTANT MESH PARAMETERS USED IN THE MESH
REFINEMENT STUDY ............................................................................................ 86
TABLE 5.2: VALIDATION CASES CONDUCTED FOR THE CFD PROCESS ............................. 87
TABLE 9.1: PARAMETER RANGES FOR TEST DATA FOR PLAIN TRAILING-EDGE FLAPS .... 124
C.D. HARLEY 8
NOMENCLATURE
NOMENCLATURE
A Axial force N
c Chord length m
CL Lift coefficient -
CP Pressure coefficient -
C 𝜌 ∞ ℎ 𝑗 𝑉𝑗2 -
Blowing coefficient = 𝑞∞ 𝑐
D Drag force N
F Force N
h Slot height m
hs MiGS height m
k Coverage factor -
L Lift force N
M Pitching moment Nm
N Normal force N
C.D. HARLEY 9
NOMENCLATURE
P Static pressure Pa
PT Total pressure Pa
q Dynamic pressure Pa
Re Reynolds number -
Rx Reactionary force N
S Area m2
t Time s
V Velocity m/s
Y Side force N
C.D. HARLEY 10
NOMENCLATURE
Subscripts
Abbreviations
2D Two-Dimensional
3D Three-Dimensional
CC Circulation Control
C.D. HARLEY 11
NOMENCLATURE
LE Leading Edge
SSD Spoiler-Slot-Deflector
TE Trailing Edge
C.D. HARLEY 12
ABSTRACT
ABSTRACT
ABSTRACT OF THESIS submitted by Christopher Donald Harley for the Degree
of Doctor of Philosophy and entitled „Aerodynamic Performance of Low Form
Factor Spoilers‟. Submitted September 2010
The development of low form factor flight controls is driven by the benefits of reducing
the installed volume of the control device and/or minimising the change in external
geometry, with particular application to flight control of low observable aircraft. For
this work, the term „low form factor‟ does not refer to the aspect ratio of the control
device rather the overall installed volume. This thesis compares the use of low form
factor geometric and fluid devices on a NACA 0015 aerofoil section through two-
dimensional numerical analysis and low speed wind tunnel experiments. The geometric
spoiler is implemented as a small (boundary layer scale) variable height tab oriented
normal to the local surface, referred to as a Micro Geometric Spoiler (MiGS). The
fluidic spoiler is implemented as an air jet tangential to the local surface acting in the
forward direction, referred to as a Counter-Flow Fluidic Spoiler (CFFS). Two
chordwise spoiler locations were considered: 0.35c and 0.65c. Numerical analysis was
undertaken using a commercial CFD code using an unsteady solver and k-omega shear-
stress-transport turbulence model. Experimental forces and moments were measured
via an overhead force balance, integrated surface pressures and pressure wake survey.
Device performance is assessed against the magnitude of control achievable compared
to macro scale spoilers and trailing edge controls (effectiveness), the ratio of
aerodynamic output to control input (efficiency or gain), the shape of control response
curve (linearity), and the degree of control cross coupling.
Results show that the MiG and CFF spoilers work by a similar mechanism based on
inducing flow separation that increases the pressure ahead of the spoiler and reduces
the pressure downstream. Increasing control input increases drag and reduces lift,
however the change in pitching moment is dependent on chordwise location.
Chordwise location has a significant effect on effectiveness, efficiency, linearity and
separability. Forward MiGS location gives the largest drag gain however the control
response is strongly nonlinear with angle of attack and there is a significant undesirable
coupling of drag with pitching moment. Aft MiGS location significantly improves
control linearity and reduces pitching moment coupling however the drag gain is much
reduced. For the CFFS, the control linearity with respect to control input and angle of
attack is good for both forward and aft locations, with the aft location giving the largest
gain for lift and drag. The control response trends predicted from numerical analysis are
good, however a calibration factor of around ½ has to be applied to the control input
momentum to match the experimentally observed gains. Furthermore numerical control
drag polars under predict the change in lift with change in drag at low blowing rates.
Through the use of a CFFS device on both the upper and lower surfaces of a wing
section it is possible to generate control drag inputs fully decoupled from both lift and
pitching moment, thus potentially simplifying device control law implementation
within an integrated yaw control system.
C.D. HARLEY 13
DECLARATION AND COPYRIGHT STATEMENTS
DECLARATION
No portion of the work referred to in the thesis has been submitted in support of an
application for another degree or qualification of this or any other university or other
institute of learning.
COPYRIGHT STATEMENT
The author of this thesis (including any appendices and/or schedules to this thesis)
owns certain copyright or related rights in it (the “Copyright”) and s/he has given The
administrative purposes.
Copies of this thesis, either in full or in extracts and whether in hard or electronic copy,
may be made only in accordance with the Copyright, Designs and Patents Act 1988 (as
amended) and regulations issued under it or, where appropriate, in accordance with
licensing agreements which the University has from time to time. This page must form
The ownership of certain Copyright, patents, designs, trademarks and other intellectual
property (the “Intellectual Property”) and any reproductions of copyright works in the
thesis, for example graphs and tables (“Reproductions”), which may be described in
this thesis, may not be owned by the author and may be owned by third parties. Such
Intellectual Property and Reproductions cannot and must not be made available for use
without the prior written permission of the owner(s) of the relevant Intellectual
C.D. HARLEY 14
DECLARATION AND COPYRIGHT STATEMENTS
commercialisation of this thesis, the Copyright and any Intellectual Property and/or
C.D. HARLEY 15
ACKNOWLEDGMENTS
ACKNOWLEDGEMENTS
I would like to start by thanking Dr Bill Crowther, firstly for offering me the
opportunity to do this PhD, secondly for his supervision throughout and finally for his
brutally honest critique where necessary. I would like to thank Dr Alistair Revell for his
CFD knowledge and support, and Dr Rodger Edwards for his guidance and support. I
must also thank the BAE Systems/EPSRC sponsored FLAVIIR project for funding this
work.
I have been fortunate enough to spend part of the PhD career at the now closed offsite
facility, Barton. I would like to thank the long serving lab technicians, Mike, Ken and
Dave who gave a wealth of engineering knowledge and general advice. This was a true
research facility where I learnt the art of research with the help of my fellow PhD
students: Russell, Ken, Steve Michie, Chip, Phil, Matt, Steve Liddle, and Paul Wilde. A
The second part of my PhD has been spent at the George Begg Building, where my
University career began. I would like to thank the technicians Mike, Dave and Lee who
collectively helped me put the wind tunnel experiment together, and where there was
always an open ear willing to listen. I would also like to thank John Ashley, Ben,
Johnny Potts, Phil and Matt who were always willing to discuss research over a brew.
I would like to thank my closest family and friends and especially Kate, who has
supported me throughout this process, but who I have undoubtedly neglected over the
past few years, especially since I have been writing and working. Finally, to my Mum
and Dad, thank you being there every step of the way over the past 28 years.
C.D. HARLEY 16
CHAPTER 1 – INTRODUCTION
1 INTRODUCTION
The aim of this chapter is to introduce the “low form-factor” spoiler class
and set it in context within the broader field of aircraft spoilers for flight
control. The low form factor devices investigated during this study are
Finally, the aim and objectives for this thesis are defined.
Whilst much of aerodynamic design is concerned with maximising the efficiency with
which attached air flows around a solid surface, there are specific operational
incidences where the flow is required to depart from or 'separate' from the surface.
Typically a device used to cause flow separation is referred to as a 'spoiler', in the sense
that it 'spoils' the smooth flow around an aerodynamically contoured body. An aircraft
spoiler is a type of geometric flight control typically mounted on the upper surface of a
wing that when deflected causes the flow to separate. This results in an increase in
drag, a loss in lift and a change in pitching moment, with the sign of the pitching
Aircraft require control about three axis for flight control, which can be grouped under
longitudinal control (changes to rate of pitch), and lateral control (changes to rate of
roll and rate of yaw). Spoiler devices are suitable for lateral aircraft control due to their
ability to rapidly deploy, produce favourable yawing moments and typically cause
independently on either aircraft wing, asymmetric operation, can be used to provide the
majority of the lateral control authority required during flight [1]. Spoilers deployed
C.D. HARLEY 17
CHAPTER 1 – INTRODUCTION
dumping, and air braking, typically used during landing [2]. Spoilers can also provide
active control for flutter suppression [3], direct lift control [4] and gust load alleviation
[5]. Figure 1.1 shows a photo of an Airbus A380 during the landing phase with all eight
spoilers fully (per wing) deployed for lift dumping and air braking.
Spoilers
Aileron
Figure 1.1: Photo of an Airbus A380 during landing phase with spoilers and high lift
devices fully deployed. Taken from [6]
The motivation for the work in this thesis is based on the development of low 'form
factor' spoiler-like flight controls on lifting surfaces located in between of the leading
and trailing edges. Low form factor devices are defined by utilising minimum wing
volume for installation of the control devices and systems, and are compatible with best
practice for low observable design [69]. For the purposes of this work „low form factor‟
should not be confused with „low aspect ratio‟ controls. Two types of low form factor
spoiler are considered in the work: Micro Geometric Spoilers (MiGS) and Counter-
C.D. HARLEY 18
CHAPTER 1 – INTRODUCTION
Flow Fluidic Spoilers (CFFS). A schematic comparing the main flow field features of
each of the spoiler types compared with a macro geometric spoiler is shown in Figure
1.2.
Figure 1.2: Comparison of macro geometric spoiler (MaGS), micro geometric spoiler
(MiGS) and counter-flow fluidic spoiler (CFFS).
For the purposes of this work a micro geometric spoiler is defined as a device whose
deployed length scale is of a similar order to the local boundary layer thickness at the
point of operation. This contrasts with conventional 'macro geometric spoilers' where
the deployed length scale is much larger than the local boundary layer thickness. The
Micro Geometric Spoilers considered here are similar in function to what are referred
to as micro geometric tabs, which have been variously used for load alleviation on
aircraft wings [7], load alleviation on helicopter rotor blades [8], and load alleviation on
wind turbine blades [9]. The work in this thesis is distinct in that it considers the use of
MiGS and CFFS as flight control devices directly in comparison with macro geometric
A Counter-flow Fluidic Spoiler is a device that produces a thin tangential jet of air on
the surface of a wing in a direction opposing the local flow direction. Tangential
blowing in general for flow control has been widely studied with applications mainly
C.D. HARLEY 19
CHAPTER 1 – INTRODUCTION
separation location on a curved trailing edge [10] or, less commonly leading edge
concerned with increasing lift for high lift applications, however there have been recent
applications where CC has been used as control effectors [12], [13]. In this work a
trailing edge dual upper/lower surface co-flow slot device is used to provide lift
modulation in a positive and negative sense. The present work is distinct to the
foregoing in that a counter-flow tangential jet is used and the aim is primarily to
produce drag through flow separation, with change in lift a secondary (but still
important) consideration. In the same way that flow topology generated by suction is
structurally different to the topology of blowing, co-flow and counter-flow topology are
Whilst MiGS and CFFS are very different from a practical implementation point of
view, there are strong similarities in the nature of the flow control input they provide, in
that they both introduce a forcing in a counter-flow tangential direction. A MiGS does
instructive to compare these devices in the same study since this sheds light on the
An aim of this work is to show that the mechanism of flow control of MiGS and CFFS
MiGS and CFFS and two-dimensional numerical analyses of CFFS have been
C.D. HARLEY 20
CHAPTER 1 – INTRODUCTION
performed and are presented in later chapters of this thesis. The independent variables
used in the investigations are the actuation input (height for MiGS, and blowing
coefficient for CFFS), angle of attack, and chordwise geometric location of the device.
The dependent variables are lift, drag and pitching moment. Aerodynamic performance
Effectiveness considers the relative magnitude of the control output for a particular
configuration, whereas efficiency (or gain) describes the ratio of aerodynamic output to
control input (where the control input for MiGS is the spoiler height
nondimensionalised with wing chord, hs/c, and for CFFS is blowing coefficient, C).
Effectiveness is typically determined by the control gain at the test condition. Control
linearity refers to the degree to which control outputs are simply proportional to control
increases the complexity of the control system implementation and may limit the
ultimate authority of the control. Control cross coupling refers to the degree to which
changes in lift, drag and pitching moment due to control inputs are correlated. Ideally,
from an implementation point of view, the control outputs should be fully uncorrelated.
For a single control device this is not possible, however by the use of an upper and
lower surface device it is possible to remove the correlation between drag and lift, and
C.D. HARLEY 21
CHAPTER 1 – INTRODUCTION
Micro Geometric and Counter-Flow Fluidic Spoilers for flight control applications.
Set the work in context within the broader field of aircraft flight control devices
C.D. HARLEY 22
CHAPTER 2 – THEORY
2 THEORY
A very simple but useful model for predicting the qualitative performance of flight
controls on lifting surfaces can be derived by considering the effect of the control on
the camber of the local aerofoil section. For a two-dimensional symmetrical aerofoil
section at a fixed low angle of attack (below the stall angle), a positive increment in
camber produces a positive increment in lift, positive increment in drag and a negative
positive increment in drag and a positive increment in pitching moment. This behaviour
is summarised in Figure 2.1. Note that for camber change through rearward mounted
spoiler deflection, the model is qualitatively correct for lift and drag, however the sign
of the pitching moment will in general depend on the chordwise location of the spoiler
(forward spoiler locations will tend to produce a pitching moment change in the
opposite sense to aft located devices). The drag considered in this case is the profile
C.D. HARLEY 23
CHAPTER 2 – THEORY
drag, however one could also consider the pressure drag which does not take into
Figure 2.1: Qualitative illustration of the effect of aerofoil camber on lift drag and
pitching moment
spoilers
Figure 2.2 shows the steady time-averaged two-dimensional streamline topology for
fluidic spoiler (NBFS, a jet exhausts from a slot normal to the local surface) and b) ii) a
counter-flow fluidic spoiler (CFFS). These diagrams are based on streamline data from
CFD analyses; refer to Chapter 5 for details on CFD methodology used. It can be seen
in Figure 2.2 that the flow field from a geometric and fluidic spoiler is broadly similar,
but with a number of detailed differences regarding the number and location of the
C.D. HARLEY 24
CHAPTER 2 – THEORY
separation and attachment points. Upstream of all spoiler locations there is a separation
point, however the fluidic spoiler cases also have an attachment point due to the jet
spoiler locations a large recirculating region is formed, however unlike the geometric
and normal blowing spoilers the CFF spoiler does not have a small recirculating region.
From the above observations it is expected that the control response of micro geometric
and fluidic spoiler devices is similar. However, for macro geometric devices at large
deflection angles the control drag response is dominated by the increase in projected
area of the spoiler, so it is likely that there will be significant differences with fluidic
This section discusses the flow field around a two-dimensional aerofoil with deflected
geometric spoiler and presents the effect of spoiler deflection on the global
C.D. HARLEY 25
CHAPTER 2 – THEORY
ESDU data sheet item on spoiler aerodynamics [2] and with some additional material
A deflected spoiler will generally cause the flow over a lifting surface to separate. The
extent and detailed characteristics of the separation will in general depend on the
effective height of the spoiler (projected frontal height), the chordwise position, the
aerofoil section, the aerofoil incidence, and the free stream Mach number and Reynolds
number. Figure 2.3 shows the typical steady flow field features around an aerofoil with
deflected rearward mounted (macro) geometric spoiler at a low positive angle of attack,
adapted from Lee and Bodapati [15]. There are two main regions of separation
identified in the figure. One termed the hinge bubble is located just ahead of the spoiler
hinge position and encompasses the lower region of the spoiler. The second region of
separation originates from the spoiler tip and emanates aft of the spoiler as a free shear
layer. If this separation stays completely detached from the aerofoil surface an
forward located spoiler the separated flow from the spoiler tip can reattach ahead of the
trailing edge. This can cause a non-linear change in lift with spoiler deflection, as
Figure 2.3: Typical flowfield due to a rearward mounted macro geometric spoiler
Figure 2.4 shows the effect of a deflected spoiler on the global aerodynamic
consistent with the qualitative model of Figure 2.1 a), an increase in drag, consistent
with the qualitative model of Figure 2.1 b), and an increase in pitching moment,
C.D. HARLEY 26
CHAPTER 2 – THEORY
consistent with the qualitative model of Figure 2.1 c). Chordwise location is an
interesting variable, such that a forward mounted upper surface spoiler at low
deflections can cause an increase in lift, due to flow reattaching downstream of the
Figure 2.4: Typical rearward mounted (apart from where stated) spoiler effect on
lift, drag, pitching moment and change in surface pressure with deflection. a) and
b) adapted from ESDU data sheet [2]. c) adapted from McLachlan et al [16]
Separation of the flow over a lifting surface causes a modification of the surface
pressure distribution and therefore changes the overall forces and moments. A
qualitative model for the change in surface pressure distribution for an aerofoil with
deflected upper surface rearward mounted spoiler is shown in Figure 2.5. There is an
increase in pressure ahead of the spoiler and a decrease in pressure downstream of the
spoiler. The pressure on the lower surface decreases with spoiler deflection. This
decrease is relatively small compared to the changes on the upper surface and is
associated with the reduction in circulation around the aerofoil, due to a change in the
C.D. HARLEY 27
CHAPTER 2 – THEORY
Figure 2.5: Effect of increasing spoiler deflection/height on the wing section surface
pressure distribution
C.D. HARLEY 28
CHAPTER 3 – LITERATURE REVIEW
3 LITERATURE REVIEW
The first part of this section will review the early development of geometric
spoilers for lateral aircraft control. This will be followed by the application
of spoilers for lateral control of tailless aircraft. The section will then
conclude with a review of fluidic spoiler type devices for lateral aircraft
control.
Hinged flight controls on the trailing edges of aerodynamic surfaces are relatively
simple to implement, are effective in producing the control moment magnitudes needed
for flight control, and the control response characteristics are generally linear over a
useful range of angle of attack and control surface deflection. As such, most aircraft in
the early part of the 20th Century used combinations of aileron, elevator and rudder for
flight control. One issue with ailerons for roll control is the drag on the down going
aileron generates a yawing moment in the adverse sense, in that it generates a yaw rate
of opposite sense to that required for a coordinated turn [17]. A further issue with use
of ailerons is that control deflection generates significant pitching moment, which tends
to twist the wing in the opposite sense to the deflection of the control. In the first
instance this leads to a reduction in aileron effectiveness, moreover with slender wings
and high speed flight this may lead to control reversal [18]. In light of these issues, a
number of studies were initiated in the 1920‟s looking at the use of alternative flight
controls for lateral control, in particular, spoiler type controls for which the yawing
C.D. HARLEY 29
CHAPTER 3 – LITERATURE REVIEW
moment is pro turn and there is significantly less pitching moment for a given change
Some of the earliest spoiler type devices included “Baffle flaps” [19] and “Projecting
flaps” [20]. Baffle flaps were leading edge spoilers on the upper surface of a wing,
aimed at decreasing lift, increasing drag and thereby causing a rolling moment and a
pro turn yawing moment. Projecting flaps were very similar to baffle flaps, however,
were not limited to only forward chordwise location, but various locations across the
wing surface. Figure 3.1 shows a typical forward mounted spoiler from the early
1930‟s. References for both forms of spoiler terminology can be traced back to the
early 1920‟s. Around the early 1930‟s these terms were quickly dropped for the more
common “spoiler” terminology, but it is unclear who first used the term.
Figure 3.1: One of the first applications of a forward mounted and forward hinged
spoiler to aircraft wings taken from Weick and Wenzinger [21]
Eliminating adverse yaw was an important driving factor in the first spoiler
NACA during this period aiming to gain a full understanding of spoiler devices for
lateral aircraft control[17], [22], [23], [19-21], [24-33], the main findings of which were
reported to the state of congress in NACA‟s 18th annual report [22]. The report states
C.D. HARLEY 30
CHAPTER 3 – LITERATURE REVIEW
An actual definition of satisfactory stability and control was not provided in the report
and was not attempted until flight investigations a few years later [23]. The flight
investigations tested control devices for their lateral aircraft control potential, in
particular ailerons and spoilers. The required flight characteristics for lateral aircraft
The maximum rolling rate is aircraft and pilot sensitive but was not given a
definitive limit.
The response of the aircraft to any movement of the lateral control surface
objectionable.
Finally, the action should be graduated so that the acceleration and maximum
rate of roll increase smoothly and regularly as the stick deflection is increased.
Wind tunnel tests performed by Weick et al [24] showed forward mounted spoilers to
reduces significantly. Figure 3.2 compares the rolling and yawing moments of a spoiler
deflected on the upper surface of a starboard wing, with a positively deflected aileron
on the starboard wing, as the angle of attack is increased. The aileron rolling moment is
increases with angle of attack up to 16 degrees. The opposite is observed in the yawing
moment response, where the spoiler yawing moment is independent of angle of attack
up to 12 degrees and the aileron yawing moment decreases with angle of attack. The
C.D. HARLEY 31
CHAPTER 3 – LITERATURE REVIEW
constant positive yawing moment of the spoiler is a favourable (pro turn) yawing
moment. However, this report lacks any comparison of the effect of pitching moment
Figure 3.2: Comparison of rolling and yawing moment coefficients obtained with
ailerons and spoilers during wind tunnel testing, taken from Weick et al [24].
by Weick et al [23]. The Fairchild 22 aircraft was used with the configurations shown
in Figure 3.3, consisting of three types of forward mounted spoiler; a rearward hinged
C.D. HARLEY 32
CHAPTER 3 – LITERATURE REVIEW
Figure 3.3: Lateral control device configurations on the Fairchild 22 airplane during
flight tests. Taken from Weick et al [23].
The aircraft time response due to inputs from these spoiler configurations and an
aileron control input is shown in Figure 3.4. Pilots reported a delay in aircraft response
to a control input from these forward mounted spoilers. The figure shows the aircraft
control response due to an aileron deflection has a lag of 0.1s and a spoiler deflection
has a lag of between 0.4s and 0.6s, based on a rate of roll of 0.1rad/s. The time lag can
be transformed into a time constant dependent on the aircraft velocity and chord length
to be, 𝑡𝑉 𝑐 = 5.9 𝑡𝑜 8.8. The spoiler lag was only observed in the rate of roll, whereas
C.D. HARLEY 33
CHAPTER 3 – LITERATURE REVIEW
no lag was observed in the rate of yaw. Whilst the authors don‟t suggest a reason for
this, it could be due to the change of local wing section surface pressure caused by the
separated flow downstream of the spoiler, whereas the main component of the drag is
probably caused by the pressure drag on the front face of the spoiler which would
adjust more quickly than the local wing section surface pressure. Subsequent research
location.
Figure 3.4: Time history curves showing the lag characteristics of various control
systems. Indicated air speed, 22m/s, full control deflection Taken from Weick et al
[23].
Wenzigner and Rogallo [25] showed the lag time of the rolling moment coefficient of a
spoiler deflection reduced by a half when a spoiler was moved from 0.29c to 0.55c, and
the maximum rolling moment reduced by ~20%. Wenzigner and Rogallo also
C.D. HARLEY 34
CHAPTER 3 – LITERATURE REVIEW
surface spoilers or deflectors, upper surface spoilers and deflector combinations, and
upper surface spoilers and deflector combinations with a slot allowing airflow in
some of the investigated spoilers in combination with ailerons and split flaps. The
i) Forward hinged spoiler, in ii) Forward hinged spoiler, in iii) Retractable spoiler, in
combination with spilt flap combination with aileron combination with split flap
i) Spoiler and Deflector, in ii) Spoiler and Deflector, in iii) Retractable spoiler and
combination with split flap combination with split flap Deflector, in combination
with split flap
c) Spoiler Deflector combinations
In terms of effect on lift, it was observed that spoilers or deflectors alone were not as
effective as the combination of a spoiler and deflector. This effectiveness was increased
C.D. HARLEY 35
CHAPTER 3 – LITERATURE REVIEW
further by adding a slot in between the spoiler and deflector such that airflow from the
lower surface is fed through the aerofoil to the upper surface behind the spoiler. This
causes a loss of pressure on the lower surface of the aerofoil, and an increase in
Following the improvement in lag time for a rearward mounted spoiler, a series of
configurations. One such lateral control configuration, shown in Figure 3.6, is the use
of a spoiler and aileron at the same spanwise station. This configuration in the form of
high lift flap and spoiler went on to become the standard configuration for transport
aircraft.
Figure 3.6: Three quarter rear-view of the test airplane as instrumented for flight
showing deflected flap, drooped aileron, deflected spoiler and open slot. Taken from
[26]
With the advent of the Second World War, aircraft technology advanced dramatically
and with it the requirements for lateral control changed. Investigations of lateral control
C.D. HARLEY 36
CHAPTER 3 – LITERATURE REVIEW
effectiveness with speed. The overall performance of which was stated “suitable to
replace ailerons”, however, with the caveat that a time response investigation were
As aircraft speeds increased, wing twist due to aileron deflection became large enough
high speed effect once again made spoilers attractive as lateral control devices due to
their lower associated pitching moments. An investigation on the effect of wing twist
caused by a spoiler deflection was performed by Fitzpatrick and Furlong [27]. Figure
3.7 shows a comparison of the effect of flaps and spoilers on the pitching moment with
the hinge point of the control surface at various chord-wise locations. The change in
pitching moment from the flap and spoiler deflection is of the same orientation; in this
case a negative flap deflection and upper surface spoiler deflection cause a nose up or
positive pitching moment. The figure shows that the negative flap deflection causes a
positive pitching moment, nose up, at all hinge locations, whereas a spoiler deflected
and located ahead of 0.4c causes a negative or nose down pitching moment. The overall
effect of a spoiler deflection on the aircraft pitching moment will depend on the spoiler
longitudinal location with respect to the aircraft centre of gravity. However, a spoiler
C.D. HARLEY 37
CHAPTER 3 – LITERATURE REVIEW
Figure 3.7: Effect of chordwise location on pitching moments for flaps and spoilers,
taken from Purser and McKinney [28].
This thesis reports and discusses the two dimensional aerodynamic performance of a
MiGS and CFFS device. For future research directions it is important to highlight the
Figure 3.8 shows the effect of wing sweep on the lateral control power of spoilers and
flaps taken from Letko [30]. The spoilers and flaps span the same amount of semi-span,
with the spoilers hinged at a constant 0.7c, and the flaps hinged at a constant 0.75c.
There is a loss in rolling-moment with angle of sweepback at a fixed Mach number and
angle of attack for both the flaps and spoilers. The rate of loss in rolling moment with
increasing sweep angle is similar for both flap and spoiler. Both the flap and spoiler
have an inherent directional stability mechanism, whereby the windward wing is more
C.D. HARLEY 38
CHAPTER 3 – LITERATURE REVIEW
Conventional wing/tail transport aircraft design has converged to the use of upper
surface spoilers typically mounted on the wing rear spar. Civil transport aircraft today,
such as the Boeing B737, use asymmetric aileron deflection and or spoiler deflection to
minimise any adverse yawing moment for roll control. Figure 3.9 shows a control
C.D. HARLEY 39
CHAPTER 3 – LITERATURE REVIEW
approach. Conventional spoilers are suitable lateral aircraft control devices particularly
during the approach phase where a reduction in lift and an increase in drag are required.
Figure 3.9 b) shows the spoilers deflected for the production of roll rate modification,
while Figure 3.9 c) shows all spoilers fully deployed for lift dumping and drag increase
on landing. What is not apparent from these pictures is that spoilers are deflected
asymmetrically for a roll manoeuvre, and deflected symmetrically for lift dumping and
air braking.
Figure 3.9: Photo taken from inside a Boeing 737 during the approach phase
Figure 3.10 shows spoilers deployed during a) cruise conditions for roll manoeuvre and
b) the landing phase for lift dumping and air braking. The Airbus A380 deploys a total
of 8 spoiler control surfaces on each wing covering 60% of the wing span during the
landing phase.
C.D. HARLEY 40
CHAPTER 3 – LITERATURE REVIEW
Figure 3.10: Spoiler deployed during a) cruise conditions of an Airbus A333 and b)
landing phase of an Airbus A380.
For tailless/finless aircraft where spoilers may be the only means of yaw control,
implementation is more varied compared to aircraft with vertical lifting surfaces for
yaw control. Tailless/finless aircraft such as the B2-spirit and diamond planform
Pegasus UCAV use split flaps, and upper and lower surface spoiler devices respectively
There have been a number of studies investigating innovative control effectors for
tailless aircraft [35], [36]. These studies aimed to implement control effectors that could
C.D. HARLEY 41
CHAPTER 3 – LITERATURE REVIEW
efficiency and improve aircraft manoeuvrability. These control aims have moved on
since the requirement for “satisfactory stability and controllability” in the 1930‟s.
that focused on “micro flow management” for aircraft control [35]. Micro flow
management refers to the use of micro geometric/fluidic devices for flow control that
can provide sufficient changes to the global aerodynamic coefficients for aircraft flight
control. Micro geometric flow devices include micro drag bumps or micro drag
generators (MDG), spoilers and splitter plates, which due to having heights of the order
of the local boundary layer, are attractive for low observable military applications.
Alone, micro flow devices may not provide comparable effectiveness compared to
conventional controls, however, surface contouring technologies that alter the pressure
field over a wing such that only a small disturbance is required to cause a control
controls devices .
The micro drag generator (MDG) concept consists of a number of micro tabs with
heights similar to the local boundary layer, distributed at regular chordwise intervals
across a spanwise station of a wing [37]. These MDG‟s extrude from the wing surface
to cause a local surface flow separation, which in turn causes an increase in wing drag.
A high profile drag modulating system such as the MDG‟s would also lend itself to a
steep-descent manoeuvre for transport aircraft where high profile drag is very
beneficial as discussed by Filippone [38]. The MDG concept is shown in Figure 3.11.
C.D. HARLEY 42
CHAPTER 3 – LITERATURE REVIEW
Figure 3.11: Micro drag generator (MDG) system showing MDGs on the upper and
lower surface adapted from Bauer [37].
Wood and Bauer suggest the performance of a MDG system would be equivalent to a
single spoiler device with the same overall projected area. A simple analysis is
performed in this thesis that calculates the effect or the drag due to the projected face of
the MiGS (similar to the MDG concept). It was found that 20% to 40% of the total drag
was due to the frontal area of the MiGS, but with high dependence on chord-wise
location. With a system of MiGS devices any amplification similar as that observed for
with increasing numbers of MiGS devices would be worthwhile, assuming the cost of
An MDG system may be attractive for its low observable characteristics; however the
wing skin is often part of the load bearing structure which may be compromised by the
large number of discontinuities in the wing. The need for additional strengthening may
lead to an increase in aircraft weight. A single MiGS or CFFS device would therefore
The Innovative Control Effectors (ICE) investigations from NASA, reported in 1996 on
the lateral control effectiveness of control devices for tailless aircraft [36]. Spoiler type
C.D. HARLEY 43
CHAPTER 3 – LITERATURE REVIEW
The effect of spoiler sweep angle on a diamond wing planform is shown in Figure 3.12.
At high angles of attack (>20o), the figure shows that aircraft response from a spoiler
deflection is dependent on spoiler sweep angle. It is not clear whether this result is
planform specific, however due to the lack of literature on effectiveness due to spoiler
Figure 3.12: Effect of spoiler sweep on the lateral control power of a 60o swept
tailless aircraft configuration. Taken from Dorset and Mehl [36]
surface spoiler configuration for tailless aircraft lateral control. The lateral effectiveness
of a spoiler and SSD applied to the same 60o swept tailless aircraft configuration is
shown in Figure 3.13. The figure shows the SSD provides nearly 50% improvement in
lateral control power compared to conventional spoilers at 20o angle of attack. SSD‟s
were also shown to remove the non-linear control response of an upper surface spoiler
at small deflections.
C.D. HARLEY 44
CHAPTER 3 – LITERATURE REVIEW
The SSD device benefits from greater effectiveness than a conventional spoiler,
however the implications with airflow movement between the upper and lower wing
surfaces provides other engineering problems that have hindered further development
of SSD‟s.
Figure 3.13: Lateral control power of a spoiler and SSD applied to a 60o swept
tailless aircraft configuration. Taken from Dorset and Mehl [36]
3.2.1 Introduction
This section reviews a number of studies investigating fluidic spoilers for lateral
aircraft control. Jets of air issuing from the lower wing surfaces were initially
investigated for VSTOL applications during the 1970‟s [39]. Particular focus was on
hover and the transition phases (from forward flight to hover, or from hover to forward
flight). Both favourable and unfavourable effects can be identified during these phases
C.D. HARLEY 45
CHAPTER 3 – LITERATURE REVIEW
of the VSTOL process. For example, during hover, a loss of lift due to jet thrust can be
caused by the jet entraining the freestream air, and during the transition phase an
increase in lift can be caused by effects of the jet-cross flow interaction on the
aerodynamics of the wing. The favourable effect during the transition phase has
motivated a number of studies investigating the potential use of jets issuing from the
Leopold et al [40] presents the data from a two-dimensional wind tunnel investigation
theory section when comparing the flow topology of geometric and fluidic spoilers.
Figure 3.14 shows a jet issuing from the lower surface of the wing with the aim of
Figure 3.14: Normal blowing fluidic spoiler concept taken from Leopold et al [40].
The model consisted of a NACA 0018 aerofoil with a slot running 90% of the span to
minimise three-dimensional effects, and located at a 0.5c. A uniform flow along the slot
was obtained by the use of a plenum and internal vanes, for which there are no
schematics.
The results show that at zero angle of attack and for a blowing coefficient, C = 0.48,
the sectional lift coefficient increases by ~1.0 (this does not include the lift due to the
C.D. HARLEY 46
CHAPTER 3 – LITERATURE REVIEW
jet thrust). A sectional lift coefficient of 1.0 is approximately equal to an angle of attack
of 10o for a NACA 0018 aerofoil. The lift was calculated from the aerofoil surface
pressure distribution, which with a low density of pressure tapping‟s can be a source of
error, however there appears to be sufficient data to capture the changes in surface
pressure. Figure 3.15 shows the model pressure distribution due to a blowing
coefficient of 0.48. The effect of a NBFS is very similar to that of a geometric spoiler,
with an increase in pressure ahead of the jet, a reduction in pressure downstream of the
jet, and a reduction in pressure on the aerofoil surface opposite to the jet.
but does not calculate the associated drag from these results. There are also no force
balance measurements and therefore no reference to jet blowing effect on the model
pitching moment coefficient. Although achieving impressive lift gains with the use of a
jet, without understanding the effect of normal blowing on drag and pitching moment
there cannot be a valid conclusion on the potential of normal blowing for lateral aircraft
CP
-2.5
-2
-1.5
-1
-0.5
0
0.5 Spoiler located at
1 mid chord, 0.5c
1.5
0 0.2 0.4 0.6 0.8 1
X/C
Figure 3.15: Pressure distribution of the normal blowing fluidic spoiler taken from
Leopold et al [40]. Blowing coefficient, C = 0.48, angle of attack, = 0o, Reynolds
number, Re = 2x105.
C.D. HARLEY 47
CHAPTER 3 – LITERATURE REVIEW
Walchli and Langan [41] applied the NBFS concept to a highly swept semi-span high
aspect ratio wing in high subsonic and transonic regimes. The normal blowing spoiler
configuration consisted of a row of holes just ahead of the flap hinge line. In all cases
tested, increasing the blowing coefficient reduced the lift with the exception at a Mach
number of 0.9, where the lift increased. This change in lift is due to influencing the
shock location. The effectiveness of blowing appears to plateaux after a small amount
applied to a realistic wing planform in the transonic regime, known to the author.
Yaw control using a NBFS was investigated by Tavella et al [42]. The concept is very
similar to that of Leopold et al and appears to be from the same research group.
Leopold et al showed that a lower surface NBFS can cause a substantial increase in
aerofoil lift. Tavella et al targets an increase in drag to cause a yawing moment, with
minimal changes in lift. A NBFS was applied to the outboard region of a rectangular
semi-span wing of constant NACA 0018 cross-section. The blowing concept is shown
C.D. HARLEY 48
CHAPTER 3 – LITERATURE REVIEW
V∞
Results show that the overall effect on lift is negligible over the angle of attack range of
-4o to +4o up to a blowing coefficient of 0.0375. The change in drag due to a blowing
coefficient was constant over the angle of attack range of -4o to +4o. It appears that the
NBFS is sensitive to wing tip effects, however this cannot be confirmed due to the lack
of lift and drag data from a comparable NBFS study. This study also does not present
the effect normal blowing on the pitching moment of this model, therefore a conclusion
on the relative effectiveness of a NBFS for aircraft lateral control cannot be made.
C.D. HARLEY 49
CHAPTER 3 – LITERATURE REVIEW
For aircraft control systems it is important to understand the effect on the lift, drag and
pitching moment due to actuation. Currently there is no literature known to the author
that combines this information for a fluidic spoiler device, and also compares the
turn) yawing moment with deflection, and relatively higher control reversal
Early wind tunnel testing showed that forward mounted spoilers were capable
ailerons. Moving the spoiler further aft reduced the aircraft response time, but
flight control may require further understanding on the potential cost to the
the lower surface of a NACA0018 aerofoil section can increase the lift by a CL
A NBFS located in near the wing tip of a semi span wing can provide increases
C.D. HARLEY 50
CHAPTER 3 – LITERATURE REVIEW
0.0375.
A general observation of most of the papers reviewed here is the lack of lift,
from a single control input. Therefore a full picture of the control characteristics
cannot be performed.
C.D. HARLEY 51
CHAPTER 4 – RESEARCH METHODS
4 EXPERIMENTAL METHODS
uncertainties.
Wind tunnel tests were performed in the open circuit Project Wind Tunnel at the
University of Manchester. The tunnel has an octagonal test section with maximum
dimensions of 1m x 1.1m and a length of 2m, Figure 4.1. The maximum test section
velocity is ~50m/s. The test section velocity was measured using calibrated static
pressure tappings in the tunnel settling chamber ahead of the contraction cone and just
ahead of the test section. The tunnel overhead six-component force/torque balance was
used to measure forces and moments. The tests were conducted at a velocity of 24.5m/s
corresponding to a Reynolds number based on model chord length of 6x105. Whilst this
Reynolds number is clearly lower than that expected for full scale application, it is
sufficiently high to avoid gross low Reynolds number effects within the angle of attack
C.D. HARLEY 52
CHAPTER 4 – RESEARCH METHODS
6 component
force/torque
balance mount
Vertically
mounted
model
1.1
m 2m
1m
End plates
V∞ grounded to
test section
Figure 4.1: Overview of the project tunnel test section showing major dimensions and
model orientation.
Details of the wind tunnel model geometry are shown in Figure 4.2. The wind tunnel
model has a NACA 0015 aerofoil section, constant chord of 350mm and a span of
730mm. Use of a symmetrical aerofoil section and upper and lower surface spoilers
locations removed the need to test at negative angles of attack. A relatively thick (15%)
section was used to increase the likely changes in drag obtained from the spoiler and
The model core was manufactured from 'blue' polyurethane foam using an outsourced
CNC hot wire cutting service. Recesses were cut out of the surface of the foam for 28
thin tube pipes for surface pressure measurement. The model skin was made from
obechi veneer skins bonded with epoxy to the foam core and finished with carbon fibre
C.D. HARLEY 53
CHAPTER 4 – RESEARCH METHODS
Figure 4.2: Cross-sectional views of the a) MiGS, b) CFFS wind tunnel model
configurations, and c) wind tunnel model details.
The model was mounted using the steel internal structure or „skeleton‟ shown in Figure
4.3. Two 10mm steel rods running the length of the model allowed 5mm thick steel
plates to be clamped to the model at either end. These plates were recessed into the
model profile. This steel structure created extra stiffness and was the mounting point
for a 25mm steel rod that connected to the force balance via a 5mm steel plate. Two
CFFS plenums were recessed on opposite sides of the model so that the slots were
located at 0.35c and 0.65c. This allowed a combined upper and lower surface CFFS
only spanned the middle 55% of the model, hence the model is only partially two-
dimensional. A photo of the model in its CFFS configuration is shown in Figure 4.4.
Due to the partial span length of the CFFS slot, appropriate scaling of the balance data
for the CFFS results is required to obtain equivalent full span (two-dimensional) data
suitable for comparison with the surface pressure measurements, momentum loss in the
wake and numerical results. The required scaling factor is the ratio of the slot length to
C.D. HARLEY 54
CHAPTER 4 – RESEARCH METHODS
span length and has been applied to all CFFS force balance measurements presented in
this thesis. This is explained in section Error! Reference source not found..
300x50x5mm
Mounting cylinder, Force balance
25mm diameter mounting plate
Leading edge of
vertically
mounted model
Steel rods
running the
Steel plate recessed in
model span,
model tip, nuts on end
725mm long,
used to clamp the model
10mm diameter
between the steel plates,
5x150x20mm
The MiGS test articles are implemented as variable height tabs oriented normal to the
local surface. Tabs were machined from aluminium box section to give heights of
0.01c, 0.02c and 0.03c. The MiGS span was set as the same as that for the CFFS, i.e.
55% of the model span so that the same correction factor was applied to both sets of
results.
C.D. HARLEY 55
CHAPTER 4 – RESEARCH METHODS
Pressure
tapping
Clamping rod piping
Pneumatic
supply pipe for Pneumatic
0.35c CFFS supply pipe for
0.65c CFFS
The design objective for the CFFS system is to provide maximum efficiency from a
given pneumatic supply with minimum installed volume. Crowther et al [44] provides a
useful discussion of the design variables for two fluidic control devices. The paper
identifies the important design features of fluidic devices, which provided the starting
C.D. HARLEY 56
CHAPTER 4 – RESEARCH METHODS
The CFFS device produces a thin uniform velocity jet sheet from a slot within a
pressurised plenum. The slot is recessed within the aerodynamic mould line of the
model in keeping with the practical aim of minimising the installed drag penalty when
the device is not operating. Once leaving the slot the jet sheet remains attached to the
curved surface that starts as tangent to the slot exit and ends as tangent to the wing
surface Figure 4.5 c). The jet sheet is encouraged to stay attached to the wing surface
downstream of the slot exit through the so called Coanda effect [45-47]. A number of
flight control applications make use of the Coanda effect including circulation control
[48], [49] and fluidic thrust vectoring [50]. The reader may refer to these references for
C.D. HARLEY 57
CHAPTER 4 – RESEARCH METHODS
Figure 4.5 shows the design of the CFFS plenum for the wind tunnel test. The key
1. The top plate is located in position using the housing knuckle and top plate
groove provided an airtight seal. The top plate was secured by two rows of
2. The support pillars were located just behind the point at which the top plate
begins to taper from 2mm to 0.5mm. This helps the formation of a uniform
3. Good practice suggests that a contraction ratio defined by the slot height to
internal plenum height of at least 10 should be used to obtain good slot flow
quality. At the location where the pillars meet the top plate the contraction
ratio is ~25. A contraction ratio of at least 10 also implies that the plenum
static pressure will be 99% of the total pressure in the flow, which allows a
method of jet velocity calculation from the plenum static pressure [44].
4. The slot lip height was defined by the limit of manufacturability as 0.2mm ±
the plenum top plate was manufactured from steel sheet, while the plenum
5. The angular range over which a coanda jet will remain attached to a curved
surface is strongly dependent on the slot curvature (h/R), with smaller values
slot curvature of 5 per cent for the CFFS slots in this study.
C.D. HARLEY 58
CHAPTER 4 – RESEARCH METHODS
experiments, both force balance and wake pressure techniques were used to improve
of the wind tunnel set-up is shown in Figure 4.6 and a diagram of the model installation
6-component Pneumatic
force/torque supply valve
balance
Vertically Wake
mounted Rake and
model mounting
Flexible
Scanivalves and pneumatic
traverse for the supply
wake survey piping
Figure 4.6: Photo of the wind tunnel setup with vertically mounted model, 6-
component force balance, scanivalves for surface pressure and wake rake for drag
measurement.
C.D. HARLEY 59
CHAPTER 4 – RESEARCH METHODS
Figure 4.7: Layout of the wind tunnel setup. Solid arrows show the direction of airflow
from the compressed supply, dashed lines indicate 0.35c CFFS location and pressure
supply piping.
The drag coefficient is calculated from the momentum loss in the wake caused by the
model drag. The wake rake technique uses integrated total pressure measurements in
the wake of a body to determine the momentum loss in the wake. For the present
experiments, the survey location was one chord length downstream of the trailing edge
of the model. This is far enough downstream that the static pressure variation across the
wake can be assumed to be negligible [52]. The wake rake consisted of 40 total
pressure tubes (single tube shown on fig 21) with a centre distance of 3mm. Two static
pressure tubes were located at the extremities of the total pressure tube array. The static
pressure holes were located eight tube diameters downstream of the tube tip as
C.D. HARLEY 60
CHAPTER 4 – RESEARCH METHODS
closely adjacent tubes the tips of the total pressure tubes were flattened perpendicular to
The wake rake was mounted to a vertical rod that was attached to a manual traverse
underneath the test section. Once the Scanivalve had stepped through each total and
static pressure probe, the manual traverse was used to move the mounting rod to the
next location, such that each wake survey overlapped the previous by at least 20% of
the wake rake measurement length (120mm). The number of movements of the wake
rake was defined by the size of the wake, a total of 2-3 movements was sufficient for
the baseline aerofoil at all angles of attack measured. Figure 4.8 shows a close up of the
wake rake Pitot and static tube array mounted in the wind tunnel.
Static
pressure
tubes
Total
pressure
tubes
Flexible
pressure
tubing
Figure 4.8: Photo of the wake rake used during wind tunnel experiments
Direct force and moment measurements were taken using the overhead 6-component
force/torque balance. The load ranges of the force balance are shown in Table 4.1.
C.D. HARLEY 61
CHAPTER 4 – RESEARCH METHODS
Assuming a 2π lift curve slope and a maximum angle of attack of 6o, the maximum lift
produced by the baseline model is ~61N, which is in the range of the balance side force
load limit. For all components the nominal force/moment range is equivalent to 1V or
+\- 1V signal from the displacement transducer. A balance calibration matrix was
obtained by loading each displacement transducer individually through the balance load
centre, located at the centre of the tunnel test section, and correlating the voltage output
to the load input. This balance calibration matrix was used to transform all balance raw
Table 4.1: Wind tunnel balance force and torque nominal ranges, in force balance
axes.
Static and total pressure data was recorded from the following sources:
All pressure measurements except for the CFFS plenums and wing surface pressure
of 5mbar. The CFFS plenum pressure was measured using the Sensortechnics
C.D. HARLEY 62
CHAPTER 4 – RESEARCH METHODS
surface static pressure was measured using the HCXM020D6 pressure transducer limited to
regulated 5V supply.
Data was sampled using the National Instruments PCI-6229 card which is capable of
up to 4 analogue output channels. The PCI-6229 card was controlled through the
the same user interface read all channels required. Scanivalves were driven by a 5V
signal from the PCI card that switched the Scanivalve to the next pressure port after
which a pressure measurement was taken before switching to the next pressure port. A
delay of 1 second was left between switching ports and taking readings to allow for a
settling time.
The following tests were performed with the aim of providing measurement calibration
and validation:
C.D. HARLEY 63
CHAPTER 4 – RESEARCH METHODS
incidence. The outcome of this test was that the model surface pressure
distribution was not symmetric at the model zero incidence. This was caused
results comparisons it was decided to use a zero incidence as the location that
gave a symmetric surface pressure distribution. The offset angle in this case
6. A trip strip analysis was performed to find the trip strip height that would be
just sufficient to cause a turbulent boundary layer over the model surface. This
These phases were split into two due to time requirements for initial validation of test
measurements. The first measurements included the model at zero lift incidence in
MiGS and CFFS configurations. This allowed measurement validation before the bulk
A flow chart of the experimental procedure (shown in Figure 4.9) is described below.
1. The model was set to the required configuration and incidence. The first zero
reading (1) was taken. For the CFFS cases the pneumatic supply was turned to
2. Once the required pressure had settled, a second zero measurement was taken.
This zero was taken so that the thrust of the jet could be calculated if required.
(For the baseline and MiGS tests no extra zero is required was taken). The
3. Once the velocity was stable, the scanivalves were started at the current wake
C.D. HARLEY 64
CHAPTER 4 – RESEARCH METHODS
rake location.
4. Once the Scanivalve results were obtained the wake rake was moved to the
next position, and the wake rake Scanivalve was started again (3). This loop
was performed until the entire model wake was recorded. The wind tunnel
5. A tunnel off zero reading was taken. For the CFFS cases then pneumatic
For the results, the zero readings of 2) and 5) were average and subtracted from the
baseline and MiGS readings to allow for bias in measurements over the testing time.
For the CFFS cases the zero readings of 1) and 6) were used.
The flow conditions in a wind tunnel are not the same as real flight conditions because
the air is bounded by walls. However, a wind tunnel aims to simulate actual flight
conditions and therefore requires boundary corrections that take the effect of the walls
into account. Boundary corrections for two-dimensional cases can be split into
C.D. HARLEY 65
CHAPTER 4 – RESEARCH METHODS
Horizontal buoyancy is the static pressure variation in the wind tunnel test section in
the streamwise direction due to the thickening of the boundary layer at the walls. The
horizontal buoyancy is usually insignificant for wing models, however, this variation is
minimised even further in the Project tunnel due to a test-section that expands in the
longitudinal axis.
The solid blockage is the ratio of the frontal area of the model to the test section area.
Using the maximum model angle of attack and assuming a constant test section area the
solid blockage comes out to 0.1. This is at the higher end of a typical solid blockage but
The wake blockage is similar to the solid blockage but in terms of the effect of the
wake of a body in the test section. Pope et al [54] provide a two-dimensional wind
tunnel test example to show the impact of the solid and wake blockages. Applying the
blockage corrections to a representative case was found to cause <1% change in the
drag measured. Since this is small compared to other measurement uncertainties basis,
corrections due to blockage effects have not been applied to wind tunnel results.
The streamline curvature refers to the alteration of the streamline curvature due to the
flow around a body in the wind tunnel. If the wing chord is less than 0.7 times the
tunnel height this effect on the distribution of lift may be neglected, the value in this
The static pressure was measured at 28 locations around the surface of the wind tunnel
model at a constant spanwise location (distribution shown in Figure 4.2 c). The
spanwise location was offset 1/8c from the model centre line to avoid potential
C.D. HARLEY 66
CHAPTER 4 – RESEARCH METHODS
pressure coefficient at each location on the wing surface was calculated [14] using:
𝑃 − 𝑃∞
𝐶𝑝 = ( 4.1 )
𝑞∞
Where: q∞, is the dynamic pressure in the test section, P∞, is the freestream static
The wind tunnel force/torque balance default set up is to measure forces and moments
vertically mounted wing which means an initial reordering of forces and moments is
required:
𝐴 𝐹𝑥
𝑌 𝐹𝑧
𝑁 𝐹𝑦
𝑇𝑥 =
𝑇𝑥
𝑇𝑦 𝑇𝑧
𝑇𝑧 𝑚𝑜𝑑𝑒𝑙 𝑎𝑥𝑒𝑠 𝑇𝑦 𝑏𝑎𝑙𝑎𝑛𝑐𝑒 𝑎𝑥𝑒𝑠
The drag, lift and pitching moment can then be transformed from the model axes to the
𝐷 𝑁 sin 𝛼 + 𝐴 cos 𝛼
( 4.2 )
𝐿 = 𝑁 cos 𝛼 − 𝐴 sin 𝛼
𝑀 𝑤𝑖𝑛𝑑 𝑎𝑥𝑒𝑠 𝑇𝑦
The balance centre is located at c/2 on the model chord line. Moments at c/4 are thus
obtained by:
C.D. HARLEY 67
CHAPTER 4 – RESEARCH METHODS
𝑐
𝑀𝑐 =𝑀−𝐿 ( 4.3 )
4 2
The forces and moments in wind axes calculated from equations ( 4.2 ) and ( 4.3 ) are
nondimensionalised using the wind tunnel dynamic pressure and appropriate reference
length/area:
𝐿
𝐶𝐿 =
𝑞∞ 𝑐
𝐷
𝐶𝐷 = ( 4.4 )
𝑞∞ 𝑐
𝑀𝑐 4
𝐶𝑀 𝑐 =
4 𝑞∞ 𝑐 2
The following equation is used to calculate the lift due to model surface pressure
distribution [14]:
𝑇𝐸 𝑇𝐸
𝑥 𝑥
𝐶𝑁𝐶𝑃 = 𝐶𝑃𝑙𝑜𝑤𝑒𝑟 𝑑 − 𝐶𝑃𝑢𝑝𝑝𝑒𝑟 𝑑
𝑐 𝑐 ( 4.5 )
𝐿𝐸 𝐿𝐸
Where, lower and upper refer to the lower and upper surface of the model, and LE and
Note that calculation of the lift coefficient from the surface pressure distribution does
not take into account any lift due to shear stress (which is very small compared to the
C.D. HARLEY 68
CHAPTER 4 – RESEARCH METHODS
For a configuration without any blowing from the model the wake drag is calculated
using [54]:
𝑃𝑇 − 𝑃 𝑃𝑇 − 𝑃 𝑦
𝐶𝐷𝑊𝑎𝑘𝑒 = 2 − 𝑑
𝑃𝑇∞ − 𝑃∞ 𝑃𝑇∞ − 𝑃∞ 𝑐 ( 4.6 )
The CFFS jet velocity is calculated from the plenum pressure using isentropic flow
relations:
𝛾−1
𝑃∞ −𝛾 2
𝑉𝑗 = 𝑎 −1
𝑃𝐶𝐹𝐹𝑆 𝛾−1 ( 4.7 )
Where a is the speed of sound, P∞ is the free stream static pressure, PCFFS is the mean
pressure in the CFFS plenum and γ is the ratio of specific heats of air (=1.4). The two-
𝜌∞ ℎ𝑗 𝑉𝑗2 ( 4.8 )
𝐶𝜇 =
𝑞∞ 𝑐
Where ρ∞ is the freestream air density, hj is the slot height, Vj is the velocity of the jet
calculated using equation ( 4.7 ), and c is the wing local chord length. This calculation
does not allow for total pressure losses in the contraction and assumes a top hat velocity
profile at the slot exit. In practice there will be some pressure losses and the exit profile
will include a boundary layer at each side. Both these effects will mean that the
calculated blowing coefficient is slightly higher than the actual delivered blowing
coefficient, however the effect is likely to be small and the level of uncertainty from
C.D. HARLEY 69
CHAPTER 4 – RESEARCH METHODS
the impact the wing surface pressure has on the jet velocity and therefore the blowing
coefficient. For an upper surface CFFS, as the wing angle of attack increases, the
pressure on the surface would decrease, causing the jet expansion and therefore an
increase in jet velocity. This complex correlation between plenum pressure, surface
pressure and jet velocity has been studied in a number of articles in [71]. Due to the
small angle of attack range over which measurements were performed in this thesis, the
impact of the surface pressure modification on the jet velocity has been assumed
negligible.
As part of the experimental analysis it is of interest to separate the aerofoil profile drag
from the drag directly due to the projected area of the spoiler. To achieve this, the MiG
spoiler drag was estimated using a combination of empirical [55] and theoretical
methods [56]. The empirical method is used to calculate the drag on a plate normal to a
calculated from laminar boundary layer theory. This approach assumes a laminar
boundary layer, however the flow is turbulent. Therefore this approach is not accurate,
but does allow one to estimate the impact of the MiGS drag due to the frontal area, and
therefore provide a method of calculating the aerofoil profile drag. The total drag of the
MiGS is equal to the drag on the MiGS within the boundary layer and the drag on the
MiGS within the freestream. Figure 4.10 shows a partially submerged MiGS by a shear
C.D. HARLEY 70
CHAPTER 4 – RESEARCH METHODS
From empirical methods the effective velocity in the boundary layer is:
2
𝑛 ( 4.9 )
𝑉𝐵𝐿 = 𝑉2
𝑛+2 ∞
Where n = 6, for a flat plate velocity profile constant, V∞ is the velocity in the
freestream.
The force on a plate normal to a wall is equal to the force due to the boundary layer,
1 1 2
𝐹𝑝𝑙𝑎𝑡𝑒 = 𝐶𝐷𝑝𝑙𝑎𝑡𝑒 2𝜌𝑉∞2 ℎ𝑒𝑓𝑓 + 𝐶𝐷𝑝𝑙𝑎𝑡𝑒 2𝜌𝑉𝐵𝐿 𝛿0.99 ( 4.10 )
Where CDplate is the drag on the flat plate, taken from boundary layer theory, and heff is
the effective height of the MiGS outside of the boundary layer thickness, shown in
layer of thickness, 0.99.. Rearranging in terms of drag coefficient and dynamic pressure
1 𝑛
𝐹𝑝𝑙𝑎𝑡𝑒 = 𝐶𝐷𝑝𝑙𝑎𝑡𝑒 2𝜌𝑉∞2 ℎ𝑒𝑓𝑓 + 𝛿0.99 𝑛+2 ( 4.11 )
C.D. HARLEY 71
CHAPTER 4 – RESEARCH METHODS
Using the force on the plate, the drag coefficient due to the plate can be calculated:
𝑛
𝐶𝐷𝑝𝑙𝑎𝑡𝑒 ℎ𝑒𝑓𝑓 + 𝛿0.99 𝑛+2 ( 4.13 )
𝐶𝐷𝑀𝑖𝐺𝑆 =
𝑐
The drag calculated, CDMiGS, is then subtracted from the measured drag coefficient to
4.3.9 Control volume analysis for CFFS wind tunnel model configuration
A control volume analysis is required to define the impact of the added momentum and
mass flow from the jet of the CFFS on the total drag. The method is limited by the
following assumptions:
The freestream flow enters the control volume from the inlet on the left and
The freestream flow at the inlet and outlet is perpendicular to the inlet and
The upper and lower boundaries are parallel to the freestream flow at the
The pressures at control volume boundaries are equal to the pressures in the
Following a method presented in [57], [58], Figure 4.11 depicts a control volume
surrounding the aerofoil with CFFS device shown by the dotted line from a to i. The
C.D. HARLEY 72
CHAPTER 4 – RESEARCH METHODS
resultant force acting on the control volume is equal to the momentum variation across
𝐹= 𝜌𝑉 ∙ 𝒅𝑺 ∙ 𝑽 ( 4.14 )
𝑠
𝑎 𝑏
−𝑝1 𝑆1 + 𝑝2 𝑆2 + (𝑝𝑗 𝑆𝑗 )𝑥 + 𝑅𝑥 = − 𝜌𝑉1 𝑑𝑦𝑉1 + 𝜌𝑉2 𝑑𝑦𝑉2 − 𝑚𝑗 𝑉𝑥𝑗 ( 4.15 )
𝑖 ℎ
The total drag of the aerofoil can be defined as the x-component of the total force acting on
the aerofoil. The total force is made up of the pressure and shear stress acting on the aerofoil
surface, the thrust of the jet acting on the aerofoil, and the pressure at the slot acting on the
𝑎 𝑏
𝐷 = − 𝑅𝑥 + 𝑚𝑗 𝑉𝑥𝑗 + (𝑝𝑗 𝑆𝑗 )𝑥 = 𝜌𝑉1 𝑑𝑦𝑉1 − 𝜌𝑉2 𝑑𝑦𝑉2 − 𝑝1 𝑆1 + 𝑝2 𝑆2
𝑖 ℎ
( 4.16 )
𝑎 𝑏 ( 4.17 )
𝜌𝑉1 𝑑𝑦 = 𝜌𝑉2 𝑑𝑦 − 𝑚𝑗
𝑖 ℎ
Assuming 𝑝1 = 𝑝2 , and using equations ( 4.16 ) and ( 4.17 ), the drag on the
aerofoil is:
𝑏
𝐷= 𝜌𝑉2 𝑉1 − 𝑉2 𝑑𝑦 − 𝑚𝑗 𝑉1 ( 4.18 )
ℎ
From the general assumptions of control volume analysis, the velocity at station 1 is equal to
𝑏
𝐷= 𝜌𝑉2 𝑉∞ − 𝑉2 𝑑𝑦 − 𝑚𝑗 𝑉𝑗 ( 4.19 )
ℎ
C.D. HARLEY 73
CHAPTER 4 – RESEARCH METHODS
Converting to coefficient form by dividing through by dynamic pressure and wing area gives:
𝑉∞
𝐶𝐷 = 𝐶𝐷𝑊𝑎𝑘𝑒 − 𝐶𝜇 ( 4.20 )
𝑉𝑗
Therefore the total drag is the drag from the wake survey minus a correction factor due
to conservation of mass. This means the thrust of the jet is a part of the wake survey
measurement. The same result is reached for circulation control aerofoils [57], [58] and
injection only aerofoils [59] whose jet is in the positive drag direction.
The slot length was varied from 100mm to 400mm in 100mm intervals and
measurements were taken using the force balance. The effect of slot length with change
in lift is shown in Figure 4.12. For the purposes of this work the gradients of the results
have been assumed within a usable tolerance. The plot shows that there are negligible
slot end effects. Therefore extending the result to a full span slot requires a scaling
C.D. HARLEY 74
CHAPTER 4 – RESEARCH METHODS
Figure 4.12: Effect of slot length on change in lift coefficient (force balance
measurement). ○ = 100mm,□ = 200mm, ∆ = 300mm, ◊ = 400mm.
An investigation into the effect of the sampling period on the drag coefficient measured
from the wake survey method was performed. The time period per sample taken at each
total pressure port in the wake rake was increased from 3 seconds upwards until the
drag calculated converged to a constant result. This investigation helps to reduce the
errors in the wake survey method from potential fluctuations of the tunnel speed. The
sampling period investigation was performed with the wind tunnel model in its baseline
Figure 4.13 shows the change in drag coefficient with increasing sampling period from
varies by 10%. Only above 6 seconds does the variation of the drag measurement begin
to settle. Based on these results a sampling period of 10 seconds was chosen to be used
for the rest of the wind tunnel testing. The same sampling period was also used for all
C.D. HARLEY 75
CHAPTER 4 – RESEARCH METHODS
Figure 4.13: The effect of sampling period on wake survey drag coefficient
measurements at a sampling rate of 1kHz.
“Scale effects” refer to differences between the experimental and actual flight operation
flow conditions. A trip strip analysis is one method of reducing the differences between
these flow conditions, by tripping the surface flow from laminar to turbulent at 0.1c for
all angles of attack tested. The location of transition is not fixed, but moves with angle
of attack. It is assumed that a transition fixed at 0.1c is suitable for this investigation
[54]. Figure 4.14 shows the results of a transition study on the model using two-
At trip strip heights of less than 0.3mm the drag coefficient increases rapidly. Beyond a
height of 0.3mm any added height leads to a constant increase in drag. This indicates
that boundary layer transition has been reached with a pinked tape of thickness 0.3mm.
The trip strip drag correction is the delta between the (fully established transition or)
chosen trip strip height to that extrapolated back to zero trip strip height. This delta
corrects for the pressure drag caused by the trip strip, while the trip strip ensures
C.D. HARLEY 76
CHAPTER 4 – RESEARCH METHODS
laminar to turbulent transition occurs at the proper location. Due to the discrepancy
with the force balance data, a trip strip height of 0.5mm was used in all following tests.
Figure 4.14: Results of a transition study using a series of two-dimensional tape (or
pinked tape) of increasing thickness located at 0.1c from the leading edge on the upper
and lower surfaces of the model. The drag coefficient values shown are measured from
the force balance and momentum loss in the wake.
Based on [44], the design objectives for the CFFS plenum were:
2. To minimise the static pressure drop from plenum entry to plenum exit
The design was tested by measuring the static pressure distribution across the span on
the internal lower surface of the plenum, shown by Figure 4.15. The four pressure
tappings are equally spaced along the span of the plenum. All pressure readings are
C.D. HARLEY 77
CHAPTER 4 – RESEARCH METHODS
30
25
20
15
10
0
1 2 3 4
Pressure tapping location
Figure 4.15: Spanwise variation of static pressure in the CFFS plenum
terms of three parameters: the mean value, the expanded uncertainty and the confidence
The method of uncertainty calculation in this thesis is based on that defined by the
Institute [60]. The following process was used to calculate the measurement uncertainty
of the lift coefficient from the force balance in the following sections:
𝑠 𝑥
𝑢 𝑥 =
𝑛 ( 4.21 )
C.D. HARLEY 78
CHAPTER 4 – RESEARCH METHODS
Where, s(x), is the standard deviation of the sample and, n, is the number of
samples taken.
in a calibration file for the device) used to then calculate the calibration standard
uncertainty.
𝑈 𝑥
𝑢𝑐𝑎𝑙 𝑥 =
𝑘 ( 4.22 )
Where, U(x), is the expanded uncertainty and, k, is the coverage factor (1.96 for
3. The sample and calibration standard uncertainties are then combined to give the
2 2
𝑢= 𝑢 𝑥 + 𝑢𝑐𝑎𝑙 𝑥 ( 4.23 )
4. The effect of the measurement on the output quantity (e.g. lift coefficient) is
𝜕𝑓
𝑐𝑖 = ( 4.24 )
𝜕𝑥𝑖
2 2 2 2
𝑢𝑐2 𝑌 = 𝑐𝑥1 𝑢 𝑥1 + 𝑐𝑥2 𝑢 𝑥2 + ⋯ ( 4.25 )
C.D. HARLEY 79
CHAPTER 4 – RESEARCH METHODS
The lift coefficient from equation ( 4.4 ) can be written in terms of its input quantities
𝑁 cos 𝛼 − 𝐴 sin 𝛼
𝐶𝐿 =
𝑞∞ 𝑐 ( 4.26 )
The input quantities are the axial force, A, the normal force, N, and the dynamic
pressure, q∞. The following table sets out the calculation of the uncertainty associated
with the force balance lift coefficient of the baseline model configuration at an angle of
attack of 6 degrees. The nominal lift coefficient is 0.6197, and the calculated expanded
uncertainty at the 95% confidence level is 9.15%. This is an acceptable uncertainty for
C.D. HARLEY 80
Expanded Standard 𝝏𝑪𝑳
xi uncertainty, K uncertainty, 𝒄𝒊 = c cu (cu)2
𝝏𝒙𝒊
U u
C.D. HARLEY
𝜕𝐶𝐿 𝑁 cos 𝛼
A 1.3133 1.96 0.67 = -0.00113 -0.00075 5.689e-7
𝜕𝐴 𝑞∞ 𝑐
𝜕𝐶𝐿 −𝐴 sin 𝛼
N 5.2921 1.96 2.7 = 0.01071 0.028919 0.000836
𝜕𝑁 𝑞∞ 𝑐
Nominal
0.6197
value
81
Table 4.2: Uncertainty calculation of the force balance lift coefficient. The expanded
CHAPTER 4 – RESEARCH METHODS
CHAPTER 5 – COMPUTATIONAL METHODS
5 COMPUTATIONAL METHODS
flow fluidic spoilers. The numerical solver is described along with the
mesh, the case set-up and convergence metrics. Results from two
validation cases are then presented, one for a macro geometric spoiler [61]
and one for a normal blowing fluidic spoiler [40], where the aim is to
The computational fluid dynamics (CFD) code Fluent (version 6.1) was used for the
was employed using the fully implicit unsteady solver and a second-order-discretisation
scheme in time and space for all variables. The Semi-Implicit Method for Pressure-
Linked Equations (SIMPLE) [62] algorithm is suitable for incompressible flows [70]
of locations in the unsteady flow field region, as well as the global aerodynamic force
coefficients. The maximum number of inner-loop iterations (per time step) was set to
20, which was sufficient to ensure inner-loop convergence (inner step convergence
criteria was set to 10-5); the computation required approximately 6,000 to 8,000 outer-
C.D. HARLEY 82
CHAPTER 5 – COMPUTATIONAL METHODS
To ensure efficient convergence, computations were initially obtained for steady state
flow up to 1000 iterations, before then being restarted using the unsteady solver. A
constant time step size of t = 1x10-4 was used, which ensures that the Courant number
is less than 1 in the majority of the domain, and a maximum Courant number
The k- Shear Stress Transport (SST) turbulence model of Menter [63] was used in all
cases presented. This selection was based on numerical studies performed by Choi et al
[64] who tested a range of turbulence models for flow past an inlay spoiler; the k-
The SST model was developed based on observations that the k-e model [65] had
difficulties at the surface of a wall due to the non-zero value of epsilon (the rate of
turbulence dissipation), whilst the standard k-omega model [66] was known to have an
unstable dependence upon freestream turbulence levels. The SST effectively employs
the k-omega model in the near-wall region, whilst switching to the k-epsilon equation
away from the wall. As such, the SST model is able to be used in the near-wall viscous
region without the need for additional modelling and without adversely influencing the
The numerical meshes used have a near wall resolution sufficient to ensure that the
non-dimensional wall distance, y+, attains a value of less than unity at all times.
C.D. HARLEY 83
CHAPTER 5 – COMPUTATIONAL METHODS
The lift and drag forces for all cases were calculated from the surface pressure and
shear-stress distribution around the aerofoil. These forces are then summed and
converted to coefficient form by dividing by the freestream dynamic pressure and the
chord length, as is shown in equation ( 4.4 ). For the CFFS case the momentum of the
jet is taken into account when calculating the lift and drag coefficients.
was used in all cases in order to minimise the number of cells and therefore the
computational time per iteration. A structured mesh block of regular rectangular cells
was used in the immediate vicinity of the aerofoil, from the surface up to a distance of
0.3c away from the aerofoil surface, so as to provide an accurate resolution of the
boundary layer, and to ensure a y+ value of less than unity at the first node away from
the wall. The boundary layer refinement region was projected downstream by 1 chord
in order to provide detailed resolution of the wake profile. Beyond this region, at a
suitable distance downstream, an unstructured mesh was employed until the end of the
domain, approximately 10 chord lengths away from the aerofoil, where the mesh ended
as a circular far-field boundary split into an inlet on the upstream side and an outlet on
the downstream side. These details are shown in Figure 5.1. Mesh refinement is also
concentrated around the jet slot, to capture the interaction of the jet with the oncoming
freestream flow.
C.D. HARLEY 84
CHAPTER 5 – COMPUTATIONAL METHODS
Figure 5.1: Computational grid used for the fluidic spoiler study. Case shown is for
spoiler at 0.65x/c
In order to ensure mesh independence of the results, three different mesh resolutions
The increase in mesh density at each state was split in the perpendicular and parallel
directions from the aerofoil surface in the structured mesh region. The increase in cell
Figure 5.2 displays the evolution of predicted lift coefficient with numerical time, and
highlights the impact of mesh refinement on the results. Both the coarse and medium
mesh sizes are completely unable to capture any periodic oscillation of the separation in
the wake. However, this was captured by the fine mesh, which was therefore selected to
C.D. HARLEY 85
CHAPTER 5 – COMPUTATIONAL METHODS
Average Courant
Mesh Number of Cells Max y+
number
Table 5.1: Overview of the important mesh parameters used in the mesh refinement
study
Figure 5.2: Effect of grid refinement on aerofoil lift coefficient, for upper surface
CFFS at 0.65c, with max. C, and = 6o.
The computational procedure described above was validated against experimental data
for two cases; a geometric and fluidic spoiler case as summarised in Table 5.2. Results
are presented in the form of the surface pressure distribution comparison between
experiment and computation. The flow was deemed to have converged when both the
velocity probes downstream of the spoiler and the total lift coefficient reached periodic
solutions.
C.D. HARLEY 86
CHAPTER 5 – COMPUTATIONAL METHODS
Reynolds
Validation Case # of Cells t Average y+
#
The flow around an aerofoil with an inlay spoiler was studied experimentally by
Consigny et al [61]. This was selected as a validation case since another relevant
computational study also investigated this case, and so it was possible to directly
compare the computational prediction [67]. The aerofoil is at zero incidence, the spoiler
is deflected to 20o and the Reynolds number with respect to chord length is 1.9x106.
A close up of the mesh used is shown in Figure 5.3. Note that the mesh refinement
increases towards the aerofoil and spoiler surfaces in order to accurately resolve the
boundary layer and flow structures in these regions. The mesh refinement downstream
of the spoiler is designed to capture the unsteady oscillation of the wake resulting from
the separation. The Courant number is less than unity for the majority of the flow field.
C.D. HARLEY 87
CHAPTER 5 – COMPUTATIONAL METHODS
Figure 5.3: Close-up of the mesh around the RA16SC aerofoil with inlay spoiler
The aerofoil profile used by Consigny et al was the experimental aerofoil RA16SC,
which is a rare geometry and was not obtainable during this work. As such the profile
was extrapolated digitally from an electronic journal paper; for this reason it was not
possible to capture the geometry of the leading edge to a satisfactorily high level of
accuracy and so minor discrepancies were expected in this region (see Figure 5.4).
However, since the primary aim of this validation case was to evaluate the ability of the
present computational method in the prediction the separated flow resulting from the
spoiler, it was deemed that this loss of accuracy in the vicinity of the leading edge was
acceptable.
As shown in Figure 5.4, the computation provides a good approximation of the surface
pressure distribution in the region of separated flow, a feature of the flow which has
discrepancy at the spoiler hinge point (0.5x/c) on the upper surface, which was also
C.D. HARLEY 88
CHAPTER 5 – COMPUTATIONAL METHODS
observed by Filippone [67]. This is most likely due to geometrical differences between
CP
-1
Lower
-0.8 surface
-0.6
-0.4
-0.2
0
0.2
Upper
0.4 surface
0.6 Present Simulation
0.8
Experimental data
1
0 0.2 0.4 0.6 0.8 X/C 1
Figure 5.4: Geometric spoiler validation case for a 16% supercritical wing with
inlay spoiler, = 20o, = 0o, and Re = 1.9x106 [61].
The flow around a normal blowing fluidic spoiler case at Re = 2x105 was reported by
Leopold et al [40]. With respect to the spoiler jet, normal blowing refers to blowing
from an orifice/slot such that the jet is normal to the local aerofoil surface. The
useful validation of the application of CFD to a case with a fluidic spoiler. The aerofoil
section is a NACA 0018, the jet slot width is 0.0067c and located at the mid chord with
The mesh used in this validation case is shown in Figure 5.5. Notice the refinement
around the aerofoil surface and particularly in the region of the jet plume. The Courant
number is less than unity for the majority of the flow field except for the refinement
C.D. HARLEY 89
CHAPTER 5 – COMPUTATIONAL METHODS
region around the jet orifice, where some larger values could not be efficiently avoided,
Figure 5.5: Close-up of the mesh around the NACA 0018 aerofoil with normal
blowing fluidic spoiler
Figure 5.6 shows a comparison of the predicted CP and that reported from the
experiment for the normal blowing case. Initial results indicated that when the jet
momentum coefficient was the same as that reported in the experiment, the numerical
prediction was poor in the vicinity of the spoiler. This could be due to differences in
the dimensionality between the experimental and numerical flow; the current
conditions. Given that the aim of this validation is to evaluate the accuracy of the
current numerical approach in the application to the flow around a fluidic spoiler, it
was decided to adjust the jet momentum coefficient to a value that was more
simulated jet momentum coefficient was found to provide the best agreement, and was
C.D. HARLEY 90
CHAPTER 5 – COMPUTATIONAL METHODS
therefore used to obtain the predictions plotted in Figure 5.6. The figure shows good
agreement upstream of the jet location; though the agreement is not quite so good
downstream of the jet. In particular, the predicted flow on the upper surface of the
aerofoil, downstream of the jet, indicates regions of pressure which are lower than
those in the experiment. This is a fairly common occurrence in the unsteady modelling
generated in the spanwise direction of a 3D flow are not able to be captured, and so a
possible sources of error, a large number of iterations were performed to ensure that
these results are sufficiently time-averaged, and grid resolution was improved as much
distributions for the CFFS cases, in the results section of this thesis, the numerical
blowing coefficient will be half of that of the experimental value (and stated in the
While it has not been possible to achieve a perfect agreement with the case of a fluidic
spoiler, a sufficiently accurate prediction is obtained, and the jet momentum coefficient
has been calibrated for this type of flow. If one were aiming to improve the accuracy of
this prediction further, one would most likely need to consider the simulation of a 3D
domain, together with a more advanced turbulence simulation technique such as Large
Eddy Simulation (LES). Both these measures increase computational cost considerably
C.D. HARLEY 91
CHAPTER 5 – COMPUTATIONAL METHODS
CP
-2.5
Lower surface
-2
-1.5
-1
-0.5
Upper
surface
0
Present Simulation
0.5
Experimental data
1
0 0.2 0.4 0.6 0.8 X/C
1
Figure 5.6: Fluidic spoiler validation case for a NACA 0018 airfoil with jet issuing
from the lower surface; = 0o, and Re = 2x105, experimental C = 0.48, simulated C =
0.24 [40].
C.D. HARLEY 92
CHAPTER 6 – RESULTS & DISCUSSION
The aim of this chapter is to present and discuss the results of the Micro
A comparison of baseline lift, drag and surface pressure between different measurement
techniques, CFD and published data for the baseline model is shown in Figure 6.1. The
experimental lift curve slope, Figure 6.1 a), is consistent between force balance data
and integrated surface pressure data, and these measurements are also consistent with
CFD data and data from the literature, providing evidence that the experimental and
computational methods are suitable for the evaluation of lift. A similar comparison for
drag is shown in Figure 6.1 b). In this case there is significant discrepancy between
drag measurement from the force balance and the other sources. It is believe that this is
due to inaccuracies in accounting for support interference and tare effects, and lack of
two-dimensionality flow towards the wing tips. As a result of this discrepancy all
measurements, apart from for the control coupling tests presented at the end of the
results section where balance data had to be used for experimental expediency. The
baseline surface pressure data, Figure 6.1 c) and d), show reasonably good agreement
computational tools used are fit for purpose at least for the benign baseline cases.
C.D. HARLEY 93
CHAPTER 6 – RESULTS & DISCUSSION
a) Lift b) Drag
Figure 6.1: Comparison of the baseline experimental and numerical lift, drag and surface
pressure distributions. Literature results from [43].
force, moment and pressure results for the Micro Geometric Spoilers (MiGS). Lift,
Drag and Moment data for varying spoiler deflection and chordwise location is shown
in Figure 6.2. The broad conclusion from these data is that the MiGS is behaving in a
similar manner to that expected from a macro scale spoiler from simple base area and
camber considerations, i.e. spoiler deflection on either the upper or lower surface
C.D. HARLEY 94
CHAPTER 6 – RESULTS & DISCUSSION
increases drag, whereas upper surface spoiler deflection reduces lift and lower surface
Focussing on the lift data, Figure 6.2 a) and b), the spoiler chordwise location makes a
relatively significant difference to the control response. For the forward location the lift
reversal at low spoiler deflections. This is similar to the effect by small deflection of
forward located macro geometric spoilers (as seen in Figure 2.4 a). At the aft location
the response is monotonic. Considering now the effect of upper or lower surface spoiler
the symmetric geometry. For positive angle of attack the change in lift from the upper
surface is always greater than the change in lift from the lower surface. However, for
the forward spoiler location on the upper surface there is significant change in control
response with incidence whereas for the aft position there is very little change. In
summary, for use of a MiGS for lift control, the aft position has better linearity
compared to forward, however the forward location is more efficient (larger gain).
Considering now the drag data, Figure 6.2 c) and d), it can be seen that deflection of a
spoiler at either the upper or lower, or fore or aft location generates an increase in drag,
(consistent with increasing base area and reduced downstream pressure recovery),
however the magnitude of the response is much larger for the forward location.
Furthermore, the magnitude of control response is proportional to alpha for the forward
alpha. In terms of upper/lower surface location, for the forward station the control drag
gain is significantly higher for the upper surface compared to the lower surface. On the
other hand, for the aft station the spoiler drag response is symmetric upper/lower. The
C.D. HARLEY 95
CHAPTER 6 – RESULTS & DISCUSSION
dotted lines in Figure 6.2 b) and c) show the total drag minus the calculated base drag
of the spoilers (as detailed in the method section) as a way of illustrating the relative
contribution of the aerofoil section and the spoiler itself to the production of drag. For
the forward located spoiler the majority of the drag comes from the aerofoil component
whereas for the aft location the majority of drag comes from the calculated base drag
component. As mentioned in the research method section, the base drag estimate will
be an over estimate since the reference velocity was based on the free stream rather
than the local boundary layer velocity, however even with this uncertainty it can be
seen that if drag is required then it may be advantageous to place MiGSs in a forward
location on an aerofoil, however tests of the time response due to a control input are
required to make sure the lag is below a minimum limit. In summary, for use of a
MiGS for drag control, the aft location has relatively poor efficiency with the drag
increment mainly due to the spoiler base drag increment. The forward location has
much higher efficiency, however the magnitude of control response is strongly coupled
The MiGS generated pitching moments about the quarter chord are shown in Figure 6.2
e) and f). Of particular note is the change in sign of the pitching moment response with
spoiler deflection between the fore and aft spoiler locations, which is observed in macro
geometric spoilers in Figure 3.7. For the fore location, upper surface spoiler deflection
produces a nose down (negative) pitching moment in association with a decrease in lift
(and an increase in drag). This implies that the centre of pressure of the loading
increment is aft of the quarter chord. For the aft spoiler location at low incidence the
pitching moment trend with spoiler deflection is opposite to the fore location, even
though the drag and lift response has the same sign as for the fore location.
C.D. HARLEY 96
CHAPTER 6 – RESULTS & DISCUSSION
Change in surface pressure distribution generated from MiGS deflection relative to the
baseline is shown in Figure 6.3 at 6o angle of attack. These results demonstrate that the
pressure behind the spoiler, consistent with the established mechanism for macro
geometric spoilers. Furthermore it can be seen that the magnitude of pressure change is
largest for the forward located spoiler on the upper surface, consistent with observed
The efficiency of the MiGS defined in section 1.1 as the change in aerodynamic
is shown in Figure 6.4. Due to the nonlinear nature of a number of the control
responses, the gain is defined by the gradient of the line at the maximum spoiler height.
Comparing the efficiency of the forward and aft spoiler locations, it can be seen that the
upper surface forward location provides the maximum efficiency for all aerodynamic
coefficients show. The aft spoiler location provides greater consistency between the
C.D. HARLEY 97
CHAPTER 6 – RESULTS & DISCUSSION
a) Lift, MiGS at 0.35x/c. Showing CL of b) Lift, MiGS at 0.65x/c. Showing CL of
20% chord flap deflected to 14o. 20% chord flap deflected to 14o.
Figure 6.2: Effect of Micro Geometric Spoiler (MiGS) location on the change in
experimental lift, drag and pitching moment with spoiler height. Angle of attach range: ○ =
0o, ∆ = 3o, □ = 6o. Corresponding surface pressure distribution plots (in Fig. 6.3) indicated
in vertical axis (ACP, BCP, CCP, DCP).
C.D. HARLEY 98
CHAPTER 6 – RESULTS & DISCUSSION
-4
Figure 6.4: Gain (efficiency) of MiGS for both lower and upper surfaces, and both
chordwise locations
C.D. HARLEY 99
CHAPTER 6 – RESULTS & DISCUSSION
Force and moment data for CFFS plotted in the same format as for the MiGS is shown
in Figure 6.5. The first observation is that the overall form of the lift and drag plots are
broadly similar for the MiGS and the CFFS, confirming that whilst the implementation
of these two types of control is very different, the fundamental control response for the
different control locations is similar. In terms of differences, the forward upper surface
CFFS has a monotonic lift control response with blowing coefficient, unlike the MiGS
in the same location. Also, for drag, the control gain for the aft CFFS is greater than for
the fore location, which is the opposite way round to the MiGS case.
In order to understand the level of „amplification‟ generated by the CFFS the overall
measured drag with the calculated jet thrust component subtracted is shown in Figure
6.5 c) and d) as dotted lines. For the aft location at higher blowing rates the fluidic gain
of the system is around 2 (one unit of momentum gives two units of drag). At the
forward location the fluidic gain is reduced, with a gain of less than two at higher
The moment data for the CFFS is shown in Figure 6.5 e) and f). The change in moment
magnitude for a given change in lift is roughly similar between the MiGS and CFFS,
however for drag, the moment coupling is significantly reduced for the aft CFFS
position compared to the best case moment coupling for the MiGS. The effect of
amplification is carried over from the drag plots into the pitching moment plots by the
dotted line, which shows the pitching moment caused by the jet. It can be seen that the
effect of the jet on the pitching moment is almost negligible compared to the pitching
Pressure distributions for the CFFS are shown in Figure 6.6. Comparison with the
pressure distributions for the MiGS shows that the basic pressure signature of the CFFS
is similar to that of the MiGS, i.e. actuation generates an increase in pressure ahead of
the device and a decrease behind. Notice, however, that the device pressure signature is
projected further ahead of the device location for the CFFS compared to the MiGS.
The CFD data is plotted against the upper surface blowing CFFS in plots Figure 6.6 a)
and b). The simulated C value is half of the experiment as discussed in the
experimental and simulated CFFS case than observed in the normal blowing case study.
The main difference is observed ahead of the spoiler, where surface flow separation
occurs. The experimental results show a further forward projected pressure signature
in the spanwise direction of a 3D flow are not able to be captured in the simulated case.
The efficiency of the CFFS is shown in Figure 6.7. As with the MiGS efficiency plot
the gain is defined by the gradient of the line at the maximum control input. Unlike for
the MiGS, the forward located CFFS provides less overall efficiency across most of the
aerodynamic coefficients than the aft location, on both the upper and lower surfaces.
This difference is expected to be caused by the lack of increase in base observed in the
MiGS case, such that the change aerodynamic coefficients are mainly caused by change
a) Lift, CFFS at 0.35x/c. Showing CL of b) Lift, CFFS at 0.65x/c. Showing CL of
20% chord flap deflected to 14o. 20% chord flap deflected to 14o.
Figure 6.5: Effect of Counter-Flow Fluidic Spoiler (CFFS) location on the change in
experimental lift, drag and pitching moment with blowing coefficient. Angle of attach
range: ○ = 0o, ∆ = 3o, □ = 6o. Corresponding surface pressure distribution plots (in Fig.
6.6) indicated in vertical axis (ACP, BCP, CCP, DCP).(z is the moment arm for the jet)
-4
In order to aid visualisation of the coupling between lift, drag and pitching moments
generated by the MiGS and CFFS, Figure 6.8 and Figure 6.9 show the lift, drag and
pitching moment data from Figure 6.2 and Figure 6.5 plotted as control response polars.
The data presented in Figure 6.8 and Figure 6.9 is interpreted as follows. In each plot
the three open circles joined by solid lines is the baseline data for zero control input at
zero, three and six degrees angle of attack. The other open symbols joined by a dashed
line are for increasing values of control deflection, with values as identified in the
legend for each figure. The angle of attack for each data point is identified by a dotted
line joining data points at constant angle of attack to a baseline angle of attack. The
Lift/Drag polars in Figure 6.8 also show comparative data from data sheets for the
MiGS and from CFD for the CFFS. These data are identified by filled symbols, linked
Consider first the Drag/Lift polars, Figure 6.8. Linking with the discussion from the
results presented on control response as a function of control input, it can be seen that
for both MiGS and CFFS control input always produces a positive increment in drag
(point on lines of constant alpha always move upwards with increasing control
deflection), however the sign of the lift change depends on control upper/lower
location, with upper surface controls reducing lift (points on lines of constant alpha
move left) and lower surface controls increasing lift (points on lines of constant alpha
move right). Comparison with extrapolated data sheet [2] values for small deflections
of a 'macro' geometric spoiler at 0.35c and 0.65c with the MiGS are shown in Figure
6.8 a ii) and Figure 6.8 a iv). For the forward case (0.35c) the agreement between the
present experimental data and the data sheet is reasonably good. At the aft location
(0.65c) the agreement is less good, with the data sheet values considerably over
predicting the drag. This is not unexpected since the data sheet methods do not take into
account the loss in spoiler effectiveness due to boundary layer emersion at small
deflections. Comparison of upper surface CFFS CFD results with experimental data are
given in Figure 6.8 b ii) and Figure 6.8 b iv). The CFD trends with increasing blowing
are consistent with the experimental data, i.e. blowing reduces lift and increases drag.
However, for the forward spoiler locations CFD over predicts both the reduction in lift
and the increase in drag for a given blowing coefficient, whereas for the aft spoiler the
CFD over predicts lift and under predicts drag. Comparing the gradients of the L/D
CFD and experimental data, shows very similar results for the aft location but not for
the fore location. This discrepancy in CFD and experimental data at the fore location
was investigated and initially thought due to improper resolution of the separation
location ahead of the jet. A number of increasingly dense meshes were used to capture
this, during which the monitoring of velocity at a number of points showed no period
excitation was captured. The forward mounted CFFS device causes the leading edge
stagnation and stagnation point ahead of the CFFS to be much closer than in the aft
case. This is thought to cause a higher frequency separation relative to the aft CFFS and
Figure 6.8: Experimental lift-drag polar plots compared with data sheet (MiGS) and CFD
(CFFS) results
Consider now the moment/Lift polars shown in Figure 6.9. The data presented in this
form confirms that the pitching moment response of a forward mounted MiGS and
CFFS are similar, however there are a noticeable few differences. A lower surface
mounted CFFS produces very little pitching moment for a given lift compared to the
MiGS. Also, unlike the MiGS the upper surface forward mounted CFFS of Figure 6.9 b
ii) produces a negative pitching moment of a similar order to the aft location, at all lift
coefficients tested. The pitching moment response from both lower surface locations
are equivalent, implying that the change of pitching moment due to CFFS actuation is
The final set of data explores the simultaneous use of both an upper surface and lower
surface CFFS, with the aim of understanding the potential for use of dual CFFS devices
for control of drag at constant lift and pitching moment. It was decided to use an upper
surface CFFS at the forward location and a lower surface CFFS at the aft station based
on the control authority demonstrated in these positions when used exclusively (only
one control operating at a time). Since there is a fair degree of decoupling between the
pressure distribution on the controlled side and the pressure distribution on the opposite
(uncontrolled) side, it is hypothesized that the control effect from non exclusive
actuation (both controls working at the same time) will be similar to the sum of the
exclusive effects. Measurement of drag, lift, and moment were obtained from the wind
tunnel force balance for varying blowing through both the upper and lower CFFS
Figure 6.10 shows contour maps of the drag data plotted using the lower surface
blowing as the x axis and upper surface blowing as the y axis. Overlaid on top of this
are loci of upper and lower surface blowing for zero change in lift and zero change in
pitching moment, labelled CM=0 and CL=0, respectively. Thus by choosing blowing
control input pairs that correlate with either of the two loci it is possible to generate a
finite drag with either zero change in pitching moment or zero change in lift, which is
the desired result from a control independence point of view. As it happens, the pitch
and lift loci are approximately overlaid for both zero and six degrees alpha, so
following either loci means that drag can be obtained with zero change in pitching
moment and zero change in lift. Whilst this is a fortuitous result, it arises because for
the present configuration there is proportionality between control lift and control
pitching moment; hence if either is driven to zero then the other will be zero also.
Finally, in comparing the achievable drag control at zero and 6 degrees alpha, it can be
seen that the magnitude of drag obtainable is significantly greater at the higher angle of
attack; however, these high values of drag cannot be obtained along the zero change in
pitching moment and lift loci. Indeed, the maximum achievable drag with no pitch and
lift coupling is less for the six degrees angle of attack case compared to the zero angle
of attack.
a) 0 degrees alpha
b) 6 degrees alpha.
Figure 6.10: Force balance drag coefficient from dual surface blowing for yaw control.
Upper surface CFFS = 0.35x/c. Lower surface CFFS = 0.65x/c. Zero pitching moment
(CM) and lift (CL) loci indicated.
OPPORTUNITIES
7.1 Conclusions
Micro geometric spoilers (MiGS) and Counter-flow fluidic spoilers (CFFS) can be
placed in a similar class of spoiler, ‘low form factor spoiler’, due to their similar
to the local boundary layer thickness at the point of operation. This is contrast to
much larger than the local boundary layer thickness. The MiGS acts as a barrier
to the smooth near wall flow effectively reducing the momentum in the co-flow
direction.
edge and in opposition to the local freestream direction and that the control is
placed on the upper or lower surface away from the leading and trailing edges.
This is significant in that leading and trailing edges are typically already highly
systems placement and hence the capability to implement new flight controls
outside these areas is advantageous. The CFFS acts as a barrier to the smooth
modifications. Geometric and fluidic spoilers have also been shown to have
differences.
Both geometric and fluidic spoilers cause a large recirculating region that
configuration and flow field conditions. Both geometric and fluidic spoilers
have small recirculating regions similar to those termed “hinge bubbles” just
ahead and behind the spoiler base or slot, apart from the CFFS, whose large
recirculating region is entrained by the jet such that there is no minor aft
recirculation.
Relatively small scale, low speed two dimensional wind tunnel experiment has
pressure tapping‟s and wake survey. The force balance drag measurements
configuration therefore the wake survey drag was used for subsequent drag data
A control volume analysis of the wind tunnel test section containing CFFS
model configuration has shown that a correction related to the additional mass
flow of the jet is required to obtain the actual drag coefficient, similar to that
CFD can provide realistic results for aerofoils at low angles of attack with
model and macro geometric spoiler configurations from literature have shown
coefficients.
Accurate correlation between CFD and experimental results for the CFFS cases
was only found when the calibration factor used in the fluidic spoiler validation
case was used. The use of a calibration factor is considered acceptable practice
on the basis that the effective discharge coefficient of the blowing slot is
unknown and hence there needs to be some correction applied to the blowing
coefficient used in the CFD and the actual blowing coefficient delivered by the
It was also found that for the forward CFFS location using the methods
described; the CFD methodology used could not qualitatively predict the global
results, it is clear that the good qualitative agreement is found for the aft CFFS
case; however there are differences in the forward location results. A mesh
refinement study failed to capture the periodic oscillation in the flow field,
URANS in capturing the high frequency oscillations that may dominate the
flow field.
CFD shows that the fundamental fluid mechanism for both MiGS and CFFS is
Both the MiGS and CFFS generate an adverse pressure gradient ahead of the
this separation is approximately proportional to the control input, that is, larger
MiGS deflection or increased blowing from the CFFS shifts the separation point
(becomes more negative) with increasing spoiler input, however, the magnitude
of the change is generally smaller than pressure change ahead of the device.
This overall mechanism is similar to the way in which macro spoiler devices
For an aerofoil at a positive lift coefficient, spoiler input on the upper surface
reduces the lift and increases the drag, whereas a spoiler on the lower surface
The fore/aft location of the spoiler devices has important effects on the device
control response and these effects are different for the MiGS and CFFS. For the
MiGS, the greatest lift and drag control gain (efficiency) is for a forward located
device on the upper surface, however the control response varies considerably
with angle of attack. This coupled with the potential lag in control response
problematic. For the CFFS, the greatest lift and drag gain is for the aft location
pitching moment generated by the MiGS is nose down for the forward location
and nose up for the aft location on both upper surfaces. For the CFFS, the
By placement of a CFFS device on the upper and lower surfaces and use of
in drag with zero change in both lift and pitching moment. This potentially
The work presented in this thesis combined the micro geometric spoiler and counter-
flow fluidic spoiler in to the low form factor spoiler class. The two-dimensional
aerodynamic performance of these two devices has been evaluated and shown to
trailing edge devices. Following on from this thesis, a number of areas may benefit
The URANS methodology applied in this thesis has shown good agreement
with experiment for aerofoils with deflected geometric spoiler, however, the
fluidic spoiler has proven more difficult to simulate. In particular two areas
the ability to capture the period excitation within the flowfield due to the fluidic
spoiler in the fore (0.35x/c) location separating the flow near the leading edge.
This thesis has evaluated low form factor devices in two-dimensions, however
applying the devices to a swept wing with and without end plates, to distinguish
Both the MiGS and CFFS are susceptible to interaction with a wide range of
radar frequencies due to the discontinuity in the aerodynamic mould line of the
wing. The MiGS discontinuity consists of a high aspect ratio physical spoiler,
whereas the CFFS discontinuity consists of a high aspect ratio slot in the wing
increase effectiveness, which in turn will increase radar signature. However, the
counter-flow fluidic spoiler has a fixed discontinuity and therefore a fixed radar
signature at all actuation conditions. Therefore the CFFS may lend itself to low
efficiency or effectiveness.
solutions are a mechanical slot lip actuator and a flexible slot lip. A
therefore minimal lag in the device activation. A flexible slot lip could
The macro geometric spoiler located towards the leading edge has been shown
to increase the lag of the aircraft response due to spoiler deployment. Moving
the spoiler further towards the trailing edge reduced the lag but also reduced the
investigation of the lag due to a MiGS and CFFS is suggested through the use of
Flapless flight refers to the ability of a aircraft to perform a full flight operation
without the use of geometric control surfaces. This can be performed by the use
of fluidic devices for aircraft control such as circulation control and fluidic
engine bleed, high control authority but low throttle phases of flight, such as the
approach and landing phase, may cause issues. A CFFS system could provide a
method which allows maximum throttle and therefore high control authority at
industry. This is due to academic studies residing in a much lower technology readiness
level (TRL) than that of industry. To bridge this gap a number of institutions and
private companies are creating collaborative projects, such as the FLAVIIR program,
driven by BAE Systems. With this in mind, and to provide an additional perspective to
the further research opportunities section, the following table presents the relevant steps
required to take low form factor technologies, in particular the CFFS device, from its
current TRL of 2/3 to a TRL of 6. Each level has been targeted separately, presenting
8 REFERENCES
[1] Mack, M.D, Seetharam, H. C., Kuhn, W. G., Bright, J. T., “Aerodynamics of
spoiler control devices,” presented at the AIAA Aircraft Systems and Technology
meeting, New York, 1979.
[2] Flower, J.W., “Lift and rolling moment due to spoilers on wings with flaps
undeflected at subsonic speeds,” ESDU, 90030b, pp. 1-30, 1999.
[3] Waszak, M. R., “Robust Multivariable Flutter suppression for benchmark active
control technology wind tunnel model,” AIAA, vol. 24, no. 1, pp. 147-153, Feb.
2001.
[4] Il, D. A., Gentry, G. L. and Stickle, J. W., Wind tunnel study of slot spoilers for
direct lift control. NASA, Langley Research Centre, NASA-TN-D-6672, 1972.
[5] Croom, D. R., Dunham, R. E., Low-speed wind-tunnel investigation of span load
alteration, forward-located spoilers, and splines as trailing-vortex hazard
alleviation devices on a transport aircraft model. NASA, Langley Research
Centre, NASA-TN-D-8133, 1975.
[6] “www.airliner.net.”
[7] Standish, K. J., Baker, J.P. and van Dam, C. P., “Two-Dimensional Wind Tunnel
and Computational Investigation of a Microtab Modified Airfoil,” Journal of
Aircraft, vol. 44, no. 2, pp. 563-572, Mar. 2007.
[8] Standish, K. J., Van Dam, C. P., “Computational analysis of a microtab-based
aerodynamic load control system for rotor blades,” Journal of the American
Helicopter Society, vol. 50, no. 3, pp. 249-258, Jul. 2005.
[9] Nakafuji, D. Y., Active Load Control & Mitigation Using Microtabs: A Wind
Energy Application. U.S. Department of Energy, 2003.
[10] Englar, R. J., “Circulation control for high lift and drag generation on STOL
aircraft,” Journal of Aircraft, vol. 12, no. 5, pp. 457-463, 1975.
[11] Bradley, R. D., Wray, W. O., “A conceptual study of leading-edge-vortex
enhancement by blowing,” presented at the 6th AIAA Fluid and Plasma Dynamics
Conference, Palm Springs, California., pp. 33-48, 1973.
[12] Crowther, W. J. Wilde, P. I. A and Harley, C. D., “Application of circulation
control for three axis control of a tailless flight vehicle,” Proceedings of the
Institute of Mechancial Engineers, vol. 224, pp. 373-386, 2009.
[13] Crowther, W. J., Wilde, P. I. A and Buonanno, A., “Integrated design of fluidic
flight controls for a gas turbine powered aircraft,” presented at the 26th AIAA
Applied Aerodynamics Conference, Honolulu, Hawaii, 2008.
[14] Anderson, J. D., Fundamentals of Aerodynamics, Third. McGraw-Hill, 2001.
[15] Lee, C. S., Bodapati, S., “Experimental investigations of the flowfield of an airfoil
with spoiler,” AIAA, vol. 25, no. 11, pp. 1411-1416, Nov. 1987.
[16] McLachlan, B. G., Karamcheti, K. and Ayoub, A., Hadjidakis, G., “A study of the
unsteady flow field of an airfoil with deflected spoiler,” presented at the AIAA
Atmospheric flight mechanics conference, Gatlinburg, Tennessee, 1983.
[17] Weick, F. E., Wenzinger, C., J., Wind-tunnel research comparing lateral control
devices, particularly at high angels of attack. I - Ordinary Ailerons on rectangular
wings. NACA, Langley Memorial Aeronautical Laboratory, NACA-TR-419,
1933.
[18] Toll, T. A., Summary of lateral-control research. NACA, Langley Research
Centre, NACA-TN-1245, 1947.
Cambridge, 1969.
[58] Williams, R. M., Harvey, J. H., Two dimensional subsonic wind tunnel tests on a
20 percent thick, 5 percent cambered circulation control airfoil. Naval ship
research and development centre, NACA-TN-AL-176, 1970.
[59] Gao, W. Zha, G. and Paxton, C. D, “Jet effects on coflow jet airfoil performance,”
Journal of Aircraft, vol. 45, no. 6, pp. 1222-1231, Jun. 2007.
[60] General metrology - Part 4: Practical guide to measurement uncertainty. British
Standards Institute, 2004.
[61] Consigny, H., Gravelle, A. and Molinaro, R., “Aerodynamic characteristic of a
two-dimensional moving spoiler in subsonic and transonic flow,” Journal of
Aircraft, vol. 21, no. 9, pp. 687-694, Sep. 1984.
[62] Patankar, S.V., Spalding, D. B., “A calculation procedure for heat mass and
momentum transfer in three-dimensional parabolic flows,” In. J. Heat Mass
Transfer, vol. 15, pp. 1787-1806, 1972.
[63] Menter, F. R., “Two-equation eddy-viscosity turbulence models for engineering
applications,” AIAA, vol. 32, no. 8, pp. 1598-1605, 1994.
[64] Chang, K. S., Choi, S. W., “Parametric Study of Transient Spoiler Aerodynamics
with Two-Equation Turbulence Models,” AIAA, vol. 38, no. 5, pp. 888-894, Sep.
2001.
[65] Launder, B. E., Sharma, B. I., “Applications of the energy dissipation model of
turbulence to the calculation of flow near a spinning disc,” Letters in Heat and
Mass Transfer, vol. 1, no. 2, pp. 131-138, 1974.
[66] Wilcox, D. C., “Re-assessment of the scale-determining equation for advanced
turbulence models,” AIAA, vol. 26, no. 11, pp. 1299-1310, 1988.
[67] Filippone, A., “Inverted jet spoilers for aerodynamic control,” Journal of Aircraft,
vol. 46, no. 4, pp. 1240-1252, 2009.
[68] Increments in aerofoil lift coefficient at zero angle of attack and in maximum lift
coefficient due to deployment of a plain trailing-edge flap, with or without a
leading-edge high-lift device, at low speeds. ESDU, data sheet 94028, 1994.
[69] Jenn, D, Radar and Laser Cross Section Enginering, Second Edition, AIAA, 2005
[70] FLUENT 6.3. users guide.
[71] Jones, G.S. Pneumatic Flap Performance for a 2D Circulation Control Airfoil, Steady
and Pulsed. 2004 NASA/ONR Circulation Control Workshop, NASA Langley
Research Center, Hampton, VA, United States, 2005.
9 APPENDIX
The theoretical lift due to the deflection of a plain flap has been compared with the
change in lift caused by the micro geometric spoiler and counter-flow fluidic spoiler in
the results section. The theoretical and empirical working used to obtain the change in
lift due to a flap deflection presented here is based on that presented in ESDU [68].
Figure 9.1 shows a schematic of a two-dimensional aerofoil and plain trailing edge flap
to which the theory refers. This theory extends into three-dimensional applications,
however, only the two-dimensional case is required in this study. The increment in lift
Where Jp is an empirical constant based on aerofoil geometry, flap is the flap deflection
1
2 2
2𝑐𝑓𝑙𝑎𝑝 2𝑐𝑓𝑙𝑎𝑝
𝑎𝑓𝑙𝑎𝑝 = 𝜋 − cos−1 −1 + 1− −1
𝑐′ 𝑐′ ( 28 )
Where cflap is the flap chord length, and c’ is the mean aerodynamic chord length.
The limiting parameters of the theory are shown in the following table. The predicted
comparison with micro geometric and fluidic spoilers is justified for this study.
Parameter Range
Table 9.1: Parameter ranges for test data for plain trailing-edge flaps
Figure 9.1: Geometric definition of an aerofoil with a plain flap control surface
for the determination of the theoretical lift coefficient