Aircraft Viscous Drag Reduction Using Riblets - Viswanath
Aircraft Viscous Drag Reduction Using Riblets - Viswanath
Aircraft Viscous Drag Reduction Using Riblets - Viswanath
Abstract
The last two decades have seen considerable research activity on the use of riblets for viscous drag reduction.
Experimental results concerning the performance of 3 M riblets on airfoils, wings and wing-body or aircraft
congurations at different speed regimes are reviewed; these applications bring in additional effects like pressure
gradients and three dimensionality. In addition to drag reduction, aspects of altered ow features due to riblets are
discussed based on detailed wind tunnel measurements at low speeds. The available results obtained from wind tunnels
as well as ight tests rmly establish the effectiveness of riblets from low speed to moderate supersonic Mach numbers.
With optimized riblets, skin friction drag reduction in the range of 58% have been measured on 2D airfoils at low
incidence and in mild adverse pressure gradients; strong evidence exist at low speeds to indicate that riblets are more
effective in adverse pressure gradients. On wings of moderate sweep relevant to transport aircraft, riblets remain
effective providing drag reduction comparable to 2D airfoils, as long as the local angle between the surface streamlines
and riblet orientation is relatively small (o101). Limited data available on wing-body congurations show that total
drag reduction of about 23% is likely. Certain suggestions for future research are outlined.
r 2002 Elsevier Science Ltd. All rights reserved.
Contents
1. Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 572
2. Some general remarks on the experimental evaluation of drag reduction due to riblets . . . 574
3. Airfoils . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 574
3.1. Airfoils at low speeds . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 574
3.2. NAL studies on airfoils at low speeds . . . . . . . . . . . . . . . . . . . . . . . . . . 576
3.3. Low-speed studies on NACA 0012 and GAW-2 airfoils . . . . . . . . . . . . . . . . . 577
3.3.1. Selection of riblets . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 577
3.3.2. Drag performance with incidence . . . . . . . . . . . . . . . . . . . . . . . . 577
3.3.3. Boundary layer development on the NACA 0012 airfoil . . . . . . . . . . . . 579
3.3.4. Turbulence intensity proles on the NACA 0012 airfoil . . . . . . . . . . . . 580
3.3.5. Reynolds shear stress proles on the NACA 0012 airfoil and quadrant analysis 581
3.4. Airfoils at transonic speeds . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 583
0376-0421/02/$ - see front matter r 2002 Elsevier Science Ltd. All rights reserved.
PII: S 0 3 7 6 - 0 4 2 1 ( 0 2 ) 0 0 0 4 8 - 9
572 P.R. Viswanath / Progress in Aerospace Sciences 38 (2002) 571600
Acknowledgements . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 598
References . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 598
at subsonic, transonic and low supersonic Mach there is enormous consistency in the available data
numbers are included. Effect of riblets on other design concerning the effectiveness of riblets in different speed
parameters like lift characteristics are discussed based on regimes and ight experiments.
limited data available. Finally, certain suggestions for
future research are outlined.
3. Airfoils
2. Some general remarks on the experimental evaluation It is relevant and appropriate to briey review early
of drag reduction due to riblets work on riblets in pressure gradients at low speeds
before moving over to airfoil applications. In most
Three different approaches have been employed in studies, the magnitude of the pressure gradient is often
literature for the determination of viscous drag reduc- described by the Clauser pressure gradient parameter b
tion arising from riblets. These are: (i) direct measure- [=(d*/tw)(dp/dx)]. While b can vary along the riblet
ment of wall shear stress using a skin friction balance; surface , it has been convenient to dene an average
(ii) use of 2D boundary layer momentum integral (MI) value of b over the riblet surface (the averaging
technique and (iii) use of an internal strain gauge procedure used may vary in different studies); the
balance (e.g. for bodies of revolution, wing-body average value of b so dened is to be understood as a
conguration). representative value for the given ow condition.
Direct measurement of wall shear stress using an Sometimes, the boundary layer momentum thickness,
accurate skin friction drag balance has been performed y; is used in place of d and the pressure gradient
in several investigations and error bands have been parameter is dened here as b (see nomenclature).
estimated. The use of 2D MI involves measurement of Choi [39] investigated the effects of longitudinal
velocity proles in the boundary layer, both ahead and pressure gradients on a at plate with machined riblets
immediately downstream of the test surface, for the (1.5 mm high and 2.5 mm pitch) for two values of
smooth as well as in the presence of riblets; this is pressure gradient parameter b of 3.1 and 0.16 at low
usually accomplished using a pitot rake. The mean skin speeds; the emphasis in the study was on the structure of
friction (over the test surface) is determined by calculat- near-wall turbulence but not on drag reduction. Based
ing the change in boundary layer momentum thickness on measurements of mean velocity, streamwise turbu-
y across the test surface. This method can provide lence intensity, wall shear stress uctuations, he sug-
satisfactory results provided the test surface is long gested that the effectiveness of riblets in reducing skin
enough to cause sufcient change in y: The application friction may remain under pressure gradients. The
of MI for ows with streamwise pressure gradients experiments carried out by Pulvin and Truong [25] in a
involves the important assumption that the pressure channel ow showed maximum viscous drag reduction
eld is virtually unaltered in the presence of riblets. of about 6% (h+B13) for mild pressure gradients
There is reasonable evidence (e.g. on airfoils at low jbjo0:1; while the drag reduction was much lower at
speeds) that changes in pressure distribution due to higher b (0.201.90). Nieuwstadt et al. [27] measured
riblets are negligible as long as the boundary layer is drag reduction using a skin friction balance over a range
attached and far from separating. Several studies have of b of 0.41.5 at low speeds ; they used machined v-
utilized the MI technique for the determination of the groove riblets (cut from PVC pipe) of size h=s=0.64
total drag of a 2D body (the classical wake survey and 0.36 mm. They observed skin friction drag reduction
method); this again involves measurement of velocity in the range of 47% (for h=0.36 mm) over the range of
proles in the wake (typically at 12 chord downstream b explored; the corresponding h+ range was 1013.
of the trailing-edge) for the smooth body as well as with These results provided support concerning the effective-
riblets. While the skin friction drag reduction can be ness of riblets in adverse pressure gradients; further-
obtained from total drag reduction (with the assumption more, an indication of a slight increase in drag reduction
that pressure drag does not change), estimate of at higher b was noted as well [27].
percentage skin friction drag reduction would require
knowledge of the pressure drag of the body. The use 3.1. Airfoils at low speeds
of an internal strain gauge balance for the measurement
of total drag of a body is relatively straight forward; of Relevant experimental details and test conditions of
course the balance output has to be properly corrected earlier investigations performed on airfoil, wings and
for internal forces, if any, to obtain the net drag force. wing-body congurations are summarized in Table 1.
Viscous drag reduction in percentage can be inferred Coustols and Cousteix [28] presented results of drag
from additional information on pressure drag. reduction on a LC100D airfoil using 3 M riblets at low
As we shall see in later sections, despite the use of the speeds. Riblet sheets of h=0.152 and 0.076 mm were
three different techniques for drag reduction assessment, tested. With riblets covering only the airfoil upper
Table 1
Experimental details
Airfoil
[28] LC 100D 0.40 5.37.95 105 0.152, 0.076 06 Skin friction drag reduction inferred from
measured total drag reduction (wake survey)
[40] NACA 0012 0.152 2.5 105 0.152, 0.076, 0.023 0 Total drag reduction measured (wake survey)
[11] CAST 7 0.20 3.3 106 (at 0.650.76 0.017, 0.023, 0.033, 0 Skin friction drag reduction inferred from
MN=0.76) 0.051 measured total drag reduction (wake survey)
[29] ADA-S1 0.15 3.0 106 0.600.76 0.018, 0.033 0.5 to 1.0 Skin friction drag reduction estimated from
measured total drag reduction (wake survey)
[32] NACA 0012 0.60 1 106 0.152, 0.076 06 Skin friction drag reduction estimated from
measured total drag reduction (wake survey)
[33] NACA 0012 0.60 1 106 0.152 6 12 Skin friction drag reduction estimated from
measured total drag reduction (wake survey)
[35] GAW-2 0.60 1 106 0.076 0 12 Total drag reduction measured (wake survey)
[54] GAW-2 0.60 1 106 0.152, 0.076 06 Skin friction drag reduction inferred from
measured total drag reduction (wake survey) ;
base drag reduction also measured
Swept wing
[46] ONERA D 0.20 2.654.25 105 0.152, 0.076, 0.051 0 22.5 Skin friction drag reduction inferred from
measured total drag reduction (wake survey)
[37] GAW-2 0.45 0.75 106 0.114 06 25 Skin friction drag reduction inferred from
measured total drag reduction (wake survey)
Wing-body
[11] AIRBUS 320 mac.0.381 3.9 106 (at 0.300.82 0.023 2 to 3 28 Total drag reduction measured (balance)
MN=0.70)
P.R. Viswanath / Progress in Aerospace Sciences 38 (2002) 571600
[30] Do228 0.508 (centre 1.373.09 106 0.076 5 to 20 Total drag reduction measured (balance)
section)
575
576 P.R. Viswanath / Progress in Aerospace Sciences 38 (2002) 571600
(or suction) surface, drag measurements were made generally lacking, even for a gross understanding of how
using wake survey over an incidence range of 061. They the ow is altered.
reported total drag reduction of about 2% at a 01 and Systematic investigations with 3 M riblets were under-
21, and no drag reduction at higher a; the corresponding taken in our laboratory, both at low speeds and
viscous drag reduction was estimated to be about 7%. transonic speeds, on several airfoils and wings which
The poor performance of riblets at higher a was are sketched in Fig. 2; certain broad experimental details
attributed partly to possible effects of boundary layer are given in Table 1. The experiments and analysis were
separation. Caram and Ahmed [40] studied the near and relatively detailed on the NACA 0012 airfoil [3133]
intermediate wake region of a NACA 0012 airfoil at covering an incidence range of 0121. The studies on the
zero incidence. They reported total drag reduction of GAW(2) airfoil [34,35] were undertaken primarily to
13.3% for h=0.152 mm and lower reduction for other conrm certain new ndings observed on the NACA
sizes (2.7% and 7.3% for h=0.076 and 0.023 mm, 0012 airfoil. Riblet effectiveness at transonic speeds was
respectively); the total drag reduction of 13.3% and assessed on a ADA-S1 supercritical airfoil [29].
7.3% which imply higher viscous drag reduction are In all the airfoil and wing congurations tested, the
much larger than the levels measured in zero or mild boundary layers on the upper and lower surfaces were
pressure gradient ows in literature; furthermore, the tripped near the leading edge and the total drag was
non-monotic variation of drag reduction with h+ raises measured (both with and without riblets) using the well-
some doubts on the accuracy of these measurements. known wake survey method; two-dimensionality of the
mean ow was rst established using the 2D momentum
3.2. NAL studies on airfoils at low speeds integral equation technique. In all the cases, the riblets
were applied typically over a streamwise distance of
The limited data base from early studies on airfoils at about 80% chord on both upper and lower surfaces. In
zero and low incidence and ows subjected to mild cases, where the surface pressure distributions were
adverse pressure gradients suggested that viscous drag measured, viscous drag reduction has been estimated
reduction comparable to those observed in zero pressure based on total drag reduction and knowledge of pressure
gradient ows was likely. However, the available drag at each a: The reference or baseline conguration
information on airfoils was inadequate to provide for drag reduction assessment was always the smooth
answers to many important issues: for example, effect airfoil (without riblets) in the presence of the same
of airfoil incidence on the performance of riblets, boundary layer trip used with riblets. Furthermore,
method for choosing optimum h+ for airfoil applica- unlike in some of the earlier studies in literature, the
tions, etc. Furthermore, information on boundary layer smooth airfoil was not covered with plain plastic lm
properties on airfoils in the presence of riblets was (about 0.1 mm thick) to compensate for the riblet
P.R. Viswanath / Progress in Aerospace Sciences 38 (2002) 571600 577
the symmetric airfoil section (NACA 0012) used and to adverse pressure gradients [42]; these results also suggest
study the behaviour of riblets on a cambered airfoil, that the increased effectiveness is not limited to v-
additional experiments [3435] were carried out on a grooves with h=s.
GAW-2 airfoil in the same wind tunnel, instrumentation The possible role of convex surface curvature on
and similar test conditions; only total drag measure- riblet effectiveness on the NACA 0012 results was
ments were made for the riblet size of 0.076 mm. The discussed by Sundaram et al. [32]. The surface curvature
results, included in Fig. 6, showed the same qualitative is identical on the upper and lower surfaces because
behaviour as the NACA 0012; maximum total reduction of symmetry; however, the curvature effects will in
of 10% occurred at a 61 with a value of b around 1.50. general be different due to differing boundary layer
It is interesting that both the airfoils tested reect peak growth. The combined effects of streamwise convex
drag reduction around a 61 which is essentially a surface curvature and mild pressure gradients on
feature associated with the test conditions of the the airfoil lower surface was found to be small from
experiments and Reynolds number in particular; at the point of view of drag reduction as seen from Fig. 7.
higher Reynolds numbers, the maximum drag reduction Furthermore, on the airfoil upper surface, the effective
may be expected to shift to a higher a: Signicant streamline curvature is likely to be smaller because
retardation of the near-wall ow arising from progres- of streamline divergence due to adverse pressure
sively increasing adverse pressure gradient with a may be gradients. Sundaram et al. [32] indicated that the
expected to limit the maximum effectiveness of riblets on adverse pressure gradient on the airfoil upper surface
an airfoil. was a dominant factor resulting in increased viscous
Additional support to the main nding of increased drag reduction on a riblet surface. It would be of
drag reduction in adverse pressure gradients is provided signicant value to assess riblet effectiveness on a convex
by the experimental results of Debisschop and Nieuw- wall in a zero pressure gradient ow in future experi-
stadt [42] on a at plate boundary layer with imposed ments.
adverse pressure gradients at low speeds. A specially
designed skin friction drag balance was used in the
experiments. For both v-groove (h=s) and trapezoidal 3.3.3. Boundary layer development on the NACA 0012
groove (h=0.5 s) riblets (manufactured from PVC airfoil
pipes), skin friction drag reduction in the range of 12 Sundaram et al. [32] examined the effects of riblets on
13% (corresponding h+ being 1016 for v-groove) was mean ow development in the boundary layer on the
observed at a value of b=2.20 (Fig. 8). They also airfoil. Velocity proles, measured at two x/c values of
concluded that the effectiveness of riblets increases in 0.25c and 0.50c on the airfoil upper surface at a 41; are
shown in Fig. 9 [32]. The velocities are higher (for
y+X20) with riblets and an increased intercept in the
log-law associated with the thickening of the sublayer is
seen; this is a feature well known in zero pressure
gradient ows [8,19,20]. Boundary layer proles mea-
sured just ahead of the trailing edge (x/c=0.964) on the
upper surface at a=01, 41 and 61 are displayed in
Fig. 10. The proles at a=01 and 41 showed features
similar to those discussed above except that the values of
A and B in the log-law were slightly different. At a=61
on the other hand, a decreased intercept with riblets
(from 4.93 to 2.10) was observed with the slope A
remaining the same: for this case, the lower intercept
was caused not by an increase in Cf or tw (as one may
guess), but from lower mean velocities for y+p200 as
may be seen from Fig. 11 (in contrast to higher velocities
observed generally); the lower mean velocities are
obviously a result of higher adverse pressure gradient.
The wall friction velocity is still lower with riblets at
a=61 (Fig. 12) and is associated with maximum viscous
drag reduction on the NACA 0012 airfoil. In the
Fig. 8. Evolution of relative skin friction vs pressure gradient; experiments of Debisschop and Nieuwstadt [42] an
open symbols, experiments by Nieuwstadt et al, 1993; closed increased intercept was observed at b=2.20 for which
symbol, Debisshop and Nieuwstadt 1996 (taken from a skin friction drag reduction of about 13% with
Debisshop and Nieuwstadt, 1996). v-grooves was recorded.
580 P.R. Viswanath / Progress in Aerospace Sciences 38 (2002) 571600
3.3.5. Reynolds shear stress profiles on the NACA 0012 the boundary layer thickness (B0.4d); the maximum
airfoil and quadrant analysis reduction was about 10% and occurred at y=dB0:10:
Normalized turbulent shear stress proles at x/ Similar level of reduction in /u0 v0 S in the wall region
c=0.964 on the airfoil upper surface at a 01 and 61 due to riblets has been measured by Walsh [22] and
are presented in Fig. 14; only a marginal reduction in / Suzuki and Kasagi [18] in zero pressure gradient ows.
u0 v0 S values was seen on the riblet surface at a 01: At Quadrant analysis of u0 v0 signals was performed [38]
a 61 (bB1.06), a noticeable reduction in /u0 v0 S on the with a view to gain some understanding of the relative
riblet surface was observed and over a large fraction of contribution of sweep and ejection events to the mean
Reynolds stress under the inuence of riblets at a 01
and 61; the measurements were made at selected y
locations in the boundary layer, the closest to the wall
being at y+B60. The results of percentage stress
(normalized by the mean turbulent shear stress) and
the time or duration associated with the four quadrants
at a 01 and 61 are presented in Figs. 15 and 16,
respectively. At a 01 (b=0.20), the duration of sweep
and ejection events are even closer in the wall region
under the inuence of riblets; the contribution to the
mean stress from sweep events seems slightly higher
while a small reduction from ejection events is notice-
able. The total stress from Q2 and Q4 motions is only
weakly reduced by riblets.
With an increased adverse pressure gradient at a 61
(bB1:06), the changes due to riblets both in percentage
duration and contributions Q2 and Q4 events are more
Fig. 12. Variation of friction velocity on upper surface at signicant. The equalization of the duration associated
x/c=0.964, NACA 0012 airfoil. with sweep and ejection motions are more spectacular
Fig. 13. Streamwise turbulent intensity proles and power spectral density of u0 on upper surface at x/c=0.964, NACA 0012 airfoil
(taken from [32]).
582 P.R. Viswanath / Progress in Aerospace Sciences 38 (2002) 571600
and extends to a greater distance from the wall contributions from Q2 and Q4 events on the riblet
(yB0:30d). A reversal in the contribution from sweep surface are normalized by the mean turbulent shear
and ejection events is clearly seen on the riblet surface; stress values on the smooth wall (Fig. 14) in order to
the contribution from sweep (Q4) events is enhanced by assess absolute changes, the following inference can be
about 10%, while a similar decrease is seen for the drawn: the contribution from sweep events will drop to
ejection (Q2) events for yo3 mm (corresponding y+ (nearly) the same level as the unribbed case, while the
being around 250). The total stress from Q2 and Q4 contribution of ejection events will show a further
(normalized by mean shear stress) remains about the decrease of about 510% (depending on the y position).
same (about 120%) as on the smooth wall. If the These results unambiguously suggest that, under high
drag reduction conditions, the ejection events are
appreciably inhibited by riblets and the reduction in
the mean Reynolds shear stress is caused essentially by
the lower contribution from ejection motion. The above
observations also indicate that a certain modication of
coherent motion (as discussed above) near the wall is a
key factor inuencing high skin friction drag reduction
from riblets. The Q1 and Q3 events are only altered to
a small extent by riblets.
In the light of the above observations, it is informative
to take stock of results obtained from quadrant analysis
on riblet surfaces in zero pressure gradient ows. Pulles
et al. [43] found that fractional contribution to Reynolds
stress from second quadrant (Q2) was reduced and
fourth quadrant (Q4) increased at y+=38; they
suggested that ejection was weaker and sweep stronger
on the grooved surface. The results from direct
numerical simulation of Choi et al. [44] have revealed
that riblets mitigate the positive Reynolds shear stress
producing events (Q2 and Q4) in drag reducing
congurations. The measurements of Benhalilou et al.
Fig. 14. Turbulent shear stress proles on upper surface at [45] made on v-grooves (s+=34, h+=17) have shown
a=01 and 61, x/c=0.964, NACA 0012 airfoil. that, close to ridge plane, contributions to the Reynolds
Fig. 15. Quadrant analysis of u0 v0 signals on upper surface at a=01, x/c=0.964, NACA 0012 airfoil.
P.R. Viswanath / Progress in Aerospace Sciences 38 (2002) 571600 583
Fig. 16. Quadrant analysis of u0 v0 signals on upper surface at a=61, x/c=0.964, NACA 0012 airfoil.
stress from ejection and sweep events are larger and number range of 0.300.815 and the Reynolds number
more frequent above the ridge (than over the smooth based on the riblet length (=550 mm) varied between
wall), while they are reduced above the valley. 3.818.7 106. 3 M riblets with h=0.023, 0.033, 0.051
In summary, available results from quadrant analysis and 0.076 mm were tested and they reported maximum
indicate that riblets modify the ejection (Q2) and sweep skin friction drag reduction of 78% (Fig. 17) and the
+
(Q4) events, while rst- and third-quadrant events are corresponding h w range was 1015 (hw =riblet height
nearly unchanged. These results imply that riblets do scaled w.r.t. inner variables of the turbulent boundary
+
alter the organized motion that is associated with layer at the wall; hw =average of hw over riblet length).
streamwise vortices near the wall; the changes seem The above results of Squire and Savill [26] and Coustols
more pronounced in adverse pressure gradients resulting and Schmitt [11] reveal that the effectiveness of riblets at
in relatively higher skin friction drag reduction. transonic Mach numbers is about as good as at low
speeds for zero and mild adverse pressure gradient
3.4. Airfoils at transonic speeds conditions.
Coustols and Schmitt [11] presented results of riblet
It is useful to review the results under zero and mild performance on a CAST 7 airfoil at zero incidence in the
pressure gradients at transonic speeds before addressing Mach number range of 0.650.76. 3 M riblets with
the performance of riblets on airfoils. Squire and Savill h=0.017, 0.023, 0.033 and 0.051 mm were applied
[26] were possibly the rst to study the effectiveness of between 15% and 100% chord length and pitot surveys
3 M riblets of varying sizes at a subsonic and transonic in the wake were utilized to measure the airfoil total
Mach number. They carried out tests in a wind tunnel at drag. Total drag reduction of about 3.3% was observed
relatively high Reynolds numbers and Mach numbers of for h=0.023 mm, while drag increases were found for
0.88 and 0.50. Skin friction reduction due to riblets was riblets with h=0.033 and 0.051 mm (Fig. 18); the results
determined by plotting the ratio of the change in y along for h=0.017 mm were about the same as at h=0.023 mm
the riblet surface to the corresponding change along the [11]. The above total drag reduction for h=0.023 mm
same length of smooth surface. Maximum skin friction implied a viscous drag reduction of about 7.58%
drag reduction of about 5% was observed for the zero occurring in the range of h
w of 1216 (hw is the mean
+
pressure gradient ow at both Mach numbers; the value of integrated hw parameters evaluated along the
corresponding mean h+ range was 1020. Riblet manipulated surface on pressure and suction sides).
performance in adverse pressure gradients at b 0:25 These results showed the effectiveness of riblets on a
was about the same as at b 0; while riblets were not transonic airfoil at zero incidence.
effective at b 0:50: Transonic evaluation of drag reduction on a ADA-S1
On a cylindrical body at zero incidence Coustols and supercritical airfoil over an incidence range of 0.51 to
Schmitt [11] made drag measurements in the Mach +1.01 was reported by Viswanath and Mukund [29].
584 P.R. Viswanath / Progress in Aerospace Sciences 38 (2002) 571600
Fig. 17. Synthesis of drag data, - - - - low-speed results (taken from [11]).
Fig. 18. Synthesis of drag measurements, CAST 7 aerofoil (taken from [11]).
Fig. 19 shows a sample of pressure distributions and reduction while a drag increase can be expected from
the corresponding h+ variations for h=0.018 and h=0.033 mm. Measurements of total drag from wake
0.033 mm; it would appear that riblets with survey method did conrm the above expectations
h=0.018 mm is a good choice for optimum drag [29]. Total drag measured for the optimized riblet
P.R. Viswanath / Progress in Aerospace Sciences 38 (2002) 571600 585
presence of crossow boundary layer, the riblet yaw The performance of riblets on an ONERA-D airfoil
angle effect could get pronounced with wing incidence. section, innitely swept at 22.51 (Fig. 22) and zero
We shall discuss some of these aspects based on the incidence was reported by Coustols [46]: 3 M riblets of
available data. h=0.152, 0.076 and 0.051 mm were used with grooves
Mclean et al. [12] reported estimates of (average) skin aligned parallel to the freestream direction and the
friction drag reduction from ight experiments on a T- manipulated area was about 85%. Maximum total drag
33 jet trainer in the Mach number range of 0.350.70; reduction of 33.5% was recorded (Fig. 23) for h w
the wing had a sweep of about 91 and 3 M riblets with (average value of h w over the manipulated length on
groove heights of 0.033 and 0.076 mm were applied only aerofoil upper surface) o10 in the chord Reynolds
on the wing upper surface covering a streamwise number range of 2.654.25 105 (the vertical lines in
distance of about 76% chord. Based on the estimates Fig. 23 refer to the experimental uncertainty); the
of momentum thickness (obtained from measured corresponding skin friction drag reduction was esti-
velocity proles), they inferred viscous drag reduction mated to be in the range of 56%. Coustols recognized
upto 6% in the h+(=s+) range of 1015. We shall that the above level of drag reduction was some what
discuss these results in detail in Section 5. lower than generally reported in 2D ows; his analysis
showed that the riblet yaw angle (j) was quite small (less
than 101 over 90% chord) suggesting that the grooves
must be effective.
Fig. 23. Variations of total drag coefcient, ONERA D aerofoil (taken from [46]).
P.R. Viswanath / Progress in Aerospace Sciences 38 (2002) 571600 587
made at x/c=0.960 on both upper and lower surfaces of lower surface (Fig. 30). Reduction in turbulent shear
the swept wing at zero incidence, which are displayed in stress values due to riblets (about 67%) in the wall
Figs. 29 and 30. As in 2D airfoil ows, riblets result in region of the upper surface can be observed whereas no
higher velocities in the boundary layer on the wing such reduction is seen on the wing lower surface
upper surface (Fig. 29) while the effects are much weaker (Fig. 30).
on the lower surface; these results indicate again that the A useful guideline that emerges is that, as long as j
wing upper (or suction) surface is the major contributor values are relatively small (o101), moderate wing sweep
to the observed drag reduction. The results of /u0 S may not adversely affect the performance of riblets. The
show a (max) reduction of about 8% due to riblets on results of /u0 S and /u0 v0 S proles under the inuence
the wing upper surface with negligible change on the of riblets (discussed above) show features qualitatively
similar to those on 2D airfoil ows at low speeds.
While certain broad conclusions on the effect of
moderate wing sweep on riblet performance emerge
from low-speed results, there is a denite need to
investigate in detail transonic effects (including aspects
of shock waveboundary layer interaction) on a trans-
port aircraft wing from the point of view of applications.
5. Flight investigations
Fig. 30. Streamwise turbulent intensity and Reynolds shear stress proles on swept wing at x/c=0.96, a=01, GAW-2 airfoil
(taken from [37]).
590 P.R. Viswanath / Progress in Aerospace Sciences 38 (2002) 571600
chord was 1.37 106 and 3.09 106, respectively. Riblet total drag reduction of a little less than 2% in the Mach
lms with h=0.076 mm were chosen which was a number range of 0.770.79, in conformity with their
compromise considering optimum h+ for the fuselage predictions [14]. There results are extremely signicant
and wings and the grooves were aligned parallel to the since they reect benets in actual ight applications.
freestream. Balance measurements (with tripping devices
on the fuselage and wings) showed drag reduction in the
range of 16%. The authors however indicated that drag 7. Riblets at supersonic speeds
reduction of 6% was rather unrealistic because of certain
inaccuracies in the balance measurements; a total drag It is reasonable to expect that effectiveness of riblets
reduction of 23% was considered more reliable. will prevail in the supersonic regime since the micro-
Recently limited tests at transonic speeds have been grooves will be well immersed in the low-speed/subsonic
performed in our laboratory on a wing-body model velocity eld atleast upto moderate supersonic Mach
whose sketch is included in Fig. 2. 3 M riblets with numbers. The potential of using riblets at supersonic
h=0.018 mm was applied on the wing upper surface and speeds has been discussed by Bushnell [6]. Limited wind
h=0.023 mm was used on the wing lower surface and on tunnel investigations including ight test results are
the cylindrical fuselage. Riblets were aligned parallel the available in supersonic ow.
freestream and their coverage on the wing-body was Robinson [47] noticed a reduction in the streamwise
about 90%. Six-component balance measurements were turbulence intensity near the wall and a thickening of the
made over an incidence range of 11 to 1.51. Circular end viscous sublayer due to 3 M riblets at a Mach number of
plates were attached at wing tips primarily to minimize 2.97. Gaudet [48] reported drag reduction results due to
the contribution of induced drag in the drag measure- 3 M riblets at a freestream Mach number of 1.25 and
ments. Preliminary examination of the data at a Mach in the range of Reynolds number of 2.411 106/m.
number of 0.76 suggests a total drag reduction of about Experiments were made on the turbulent boundary layer
34%; detailed analysis of the results is in progress. developing along the wall of the wind tunnel using
Flight evaluation of riblet performance on an Airbus riblets of h=0.051 mm. Maximum local skin friction
A320 aircraft has been reported by Szodruch [14]. About reduction upto 7% was measured using a drag balance.
70% of the aircraft surface was covered with 3 M riblets Coustols and Cousteix [49] assessed the performance
(Fig. 36) and total drag reduction was assessed based on of 3 M riblets on a cone-cylinder body at zero incidence
fuel burn saving in ight tests. These tests revealed a at freestream Mach numbers of 1.60, 2.0 and 2.50; the
test Reynolds number range was 4.9 to 22.3 106/m.
Riblet lms with height of 0.033, 0.051 and 0.076 mm
were investigated and accurate drag measurements were
made using a one-component balance. Typical drag
results for h=0.051 mm are shown in Fig. 37. At MN=
2.0 and 2.5, net drag reduction was observed for all three
values of h; at MN=1.60, some drag increase was
recorded for unit Reynolds numbers >15 106. These
results implied maximum viscous drag reduction upto
about 6% which occurred around a value of h w (average
value of h w along the manipulated length) of 10. It may,
therefore, be noted that the supersonic wind tunnel test
results of skin friction drag reduction in nominally zero
pressure gradient conditions are in good agreement with
the available data at low and transonic speeds.
The ight test results reported by Zuniga et al. [50]
provide additional support regarding the effectiveness of
3 M riblets at supersonic speeds. Riblets with height of
0.003 in (0.076 mm) and 0.0013 in (0.033 mm) were xed
on a F-104G ight test xture (FTF) (Fig. 38), and they
were tested in a design Reynolds number range selected
to provide maximum skin friction drag reduction. The
FTF experimental setup is presented in Fig. 39 and
riblets were applied to the right side test surface only.
Boundary layer rakes were employed to measure the
velocity proles for the smooth and riblet surfaces and
Fig. 36. A-320 no.1 covered with riblet lm (taken from [14]). the average wall skin friction over the test surface was
P.R. Viswanath / Progress in Aerospace Sciences 38 (2002) 571600 593
Fig. 39. Flight test xture experiment setup (taken from [50]).
(as much as 1020%) and Reynolds shear stress of h+ calculated for the measured pressure distributions
(about 1015%). It was therefore conjectured [54] on the airfoil upper surface showed that riblets with
that the combination of lower mean velocity gradient h=0.076 and 0.152 mm would be optimum at a 01
(@u=@y) close to the wall and reduced levels of /u0 S and and 61 (Fig. 44), respectively, considering viscous drag
/u0 v0 S in the wall region of the approaching reduction. They made measurements of base pressure,
boundary layer (ahead of the blunt base) may model surface pressures and total drag using wake
favourably affect the shear layer development because survey.
the mixing zone is relatively short (comparable to The base pressure coefcient (Fig. 45) for the baseline
trailing-edge thickness) on an airfoil; reduced mixing airfoil (without riblets) was positive at all a indicating a
in the free shear layer could lead to lower velocity base thrust. Interestingly, the base pressure showed an
along the dividing streamline resulting in higher base increase with h in the a range considered and it was as
pressure. high as 50% at a=61 for h=0.152 mm. Although the
With this in background, Channaraju and Viswanath increase in base thrust was large, the base drag
[54] made an assessment of riblets for base drag reduction, however, as a fraction of total drag of the
reduction on a GAW(2) airfoil (See Fig. 2) which has airfoil was only about 0.7% since the base drag
a trailing-edge thickness of 0.5% chord; on a model component of the baseline airfoil was itself small. Even
scale with a chord of 600 mm, the trailing-edge thickness this low level of drag reduction is of engineering value
was 3 mm. Riblet lms with a height of 0.076 and since it is an added benet due to riblets on a blunt
0.152 mm were applied between 0.12c and 0.96c on the trailing-edge airfoil. The increase in base pressure due to
airfoil top and bottom surfaces. Streamwise variations riblets was attributed to be a direct consequence of two
P.R. Viswanath / Progress in Aerospace Sciences 38 (2002) 571600 595
Fig. 40. Average skin-friction coefcient for 0.0030 in riblets compared with unit Reynolds number for Mach 1.21.6
(taken from [50]).
factors: a lower effective base height of the airfoil turbulent skin friction drag reduction on airfoils, wings
(including the boundary layer displacement thickness) and wing-body combinations in different speed regimes;
and altered shear layer development characteristics as these applications bring in issues of riblet performance
discussed above. It would be of signicant engineering in pressure gradients and in the presence of three-
interest to assess base drag reduction from riblets on a dimensionality. Based on the available experimental
supercritical airfoil with a blunt trailing edge at data, certain broad conclusions are drawn, which are
transonic speeds. informative both from the point of view of design
applications as well as ow features associated with
riblets.
10. Concluding remarks The data base generated in wind tunnels and ight
experiments at high Reynolds number rmly establish
Research investigating the performance of riblets for the effectiveness of riblets from low speed to moderate
viscous drag reduction has received considerable atten- supersonic Mach numbers. Taking into account (typical)
tion during the last two decades. Riblet lms with uncertainties in the measured data, the performance of
adhesive backing manufactured by 3 M company, USA riblets may be judged to be about the same across the
have been utilized very widely in riblet research both in Mach number range investigated. With optimized
wind tunnels and in ight tests. In this review, we have riblets, skin friction drag reduction in the range of
addressed primarily the effectiveness of 3 M riblets for 58% are achievable on 2D airfoils at low incidence and
596 P.R. Viswanath / Progress in Aerospace Sciences 38 (2002) 571600
[18] Suzuki Y, Kasagi N. Turbulent drag reduction mechanism [36] Sundaram S, Viswanath PR, Subashchandar N. Studies on
above a riblet surface. AIAA J 1994;32(9):178190. drag reduction using riblets on an innite swept wing.
[19] Choi KS. Near-wall structures of a turbulent boundary National Aerospace Laboratories Report PD-EA-9802,
layer with riblets. J Fluid Mech 1989;208:41758. India, 1998.
[20] Choi KS. Turbulence structure revisited: results and [37] Sundaram S, Viswanath PR, Subashchandar N. Viscous
implications from riblets research. In: So RMC, Speziale drag reduction using riblets on a swept wing. AIAA J
CG, Launder BE, editors. Near wall turbulent ows. 1999;37(7):8516.
Amsterdam: Elsevier Science Publishers, 1993. p. 699707. [38] Viswanath PR. Riblets on airfoils and wings: a review.
[21] Bechert DW, Bartenwerfer M. The viscous ows on AIAA Paper 99-3402, 1999.
surfaces with longitudinal ribs. J Fluid Mech 1989; [39] Choi KS. Effect of longitudinal pressure gradient on
206:10530. turbulent drag reduction with riblets. In: Coustols E,
[22] Walsh MJ. Drag characteristics of v-groove and transverse editor. Turbulence control by passive means. Dordrecht,
curvature riblets. In: Hough GR, editor. Progress in Netherlands: Kluwer Academic, 1990. p. 10921.
astronautics and aeronautics. New York: AIAA, 1980. p. [40] Caram JM, Ahmed A. Effects of riblets on turbulence in
16884. the wake of an airfoil. AIAA 1991;29(11):176970.
[23] Vukoslavcevic P, Wallace JM, Balint JL. Viscous drag [41] Desai SS, Kiske S. A computer program to calculate
reduction using streamwise aligned grooves. AIAA J turbulent boundary layer and wakes in compressible ow
1992;30(4):111922. with arbitrary pressure gradient based on Greens lag-
[24] Park SR, Wallace JM. Flow-eld alteration and viscous entrainment method. Ruhr University. Bericht No. 89/
drag reduction by riblets in a turbulent boundary layer. 1982, Bochum, Germany, 1982.
AIAA Paper 93-3256, 1993. [42] Debisschop JR, Nieuwstadt FTM. Turbulent boundary
[25] Pulvin Ph, Truong TV. Riblets in internal ows with layer in an adverse pressure gradient: effectiveness of
adverse pressure gradients. In: Gyr A, editor. Proceedings riblets. AIAA J 1996;34(5):9327.
of the IUTAM Symposium on Structure of Turbulence [43] Pulles CJA, Krishna Prasad K, Nieuwstadt FTM.
and Drag Reduction. Berlin: Springer, 1990. p. 56976. Turbulence measurements over longitudinal micro-
[26] Squire LC, Savill AM. Drag measurements on planar grooved surfaces. Appl Sci Res 1989;46:197208.
riblet surfaces at high subsonic speeds. Appl Sci Res [44] Choi H, Moin P, Kim J. Direct numerical simulation of
1989;46:22943. turbulent ow over riblets. J Fluid Mech 1993;255:50339.
[27] Nieuwstadt FTM, Wolthers W, Leijdens H, Krishna [45] Benhalilou M, Anseimet F, Fulachier L. Conditional
Prasad K, Schwarz-van Manen A. The reduction of skin Reynolds stress on a v-grooved surface. Phys Fluids
friction by riblets under the inuence of an adverse 1994;6(6):210117.
pressure gradient. Expt Fluids 1993;15:1726. [46] Coustols E. Performance of internal manipulators in
[28] Coustols E, Cousteix J. Experimental investigation of subsonic three-dimensional ows. In: Choi KS, editor.
turbulent boundary layers manipulated with internal Recent developments in turbulence management. Dor-
devices: riblets. In: Gyr A, editor. Proceedings of the drecht, Netherlands: Kluwer Academic, 1991. p. 4364.
IUTAM Symposium on Structure of Turbulence and Drag [47] Robinson SK. Effects of riblets on turbulence in a
Reduction. Berlin: Springer, 1990. p. 57784. supersonic boundary layer. AIAA Paper 88-2526, 1988.
[29] Viswanath PR, Mukund R. Turbulent drag reduction [48] Gaudet L. Properties of riblets at supersonic speed. Appl
using riblets on a supercritical airfoil at transonic speeds. Sci Res 1989;46:24554.
AIAA J 1995;33(5):9457. [49] Coustols E, Cousteix J. Performances of riblets in the
[30] Van Der Hoven JG, Bechert DW. Experiments with a 1:4.2 supersonic regime. AIAA J 1994;32(2):4313.
model of a commuter aircraft with riblets in a large wind [50] Zuniga FA, Anderson BT, Bertelrud A. Flight test results
tunnel. In: Choi KS, editor. Recent developments in of riblets at supersonic speeds. NASA Tech. Memo 4387,
turbulence management. Dordrecht, Netherlands: Kluwer 1992.
Academic Publishers, 1991. p. 324. [51] Nash JF. A discussion of two-dimensional turbulent base
[31] Sundaram S, Viswanath PR, Rudrakumar S. Studies on ows. Aeronautical Research Council Research Memor-
turbulent drag reduction on a NACA 0012 airfoil using anda 3468, UK, 1965.
riblets. National Aerospace Laboratories Report PD-EA- [52] Chapman DR. An analysis of base pressure at supersonic
9401, India, 1994. velocities and comparison with experiment. NACA Report
[32] Sundaram S, Viswanath PR, Rudrakumar S. Viscous drag 1051, 1951.
reduction using riblets on a NACA 0012 airfoil to [53] Korst HH. A theory of base pressure in transonic and
moderate incidence. AIAA J 1996;34(4):67682. supersonic ow. J Appl Mech 1956;23:593600.
[33] Subashchandar N, Rajeev K, Sundaram S. Drag reduction [54] Channaraju, Viswanath PR. Base drag reduction caused
due to riblets on NACA 0012 airfoil at higher angles of by riblets on a GAW-2 airfoil. J Aircr 1998;35(6):98892.
attack. National Aerospace Laboratories Report PD-EA- [55] Lynch FT, Klinge MD. Some practical aspects of viscous
9504, India, 1995. drag reduction concepts. SAE Technical Paper 912129, 1991.
[34] Subashchandar N, Rajeev K, Sundaram S. Drag reduction [56] Wilkinson SP, Anders JB, Lazos BS, Bushnell DM.
due to riblets on a GAW(2) airfoil. National Aerospace Turbulent drag reduction research at NASA Langley
Laboratories Report PD-EA-9601, India, 1996. progress and plans. Proceedings of International Con-
[35] Subashchandar N, Rajeev K, Sundaram S. Drag reduction ference on Turbulent Drag Reduction by Passive Means.
due to riblets on a GAW(2) airfoil. J Aircr 1999;36(5):8902. London: Royal Aeronautical Society, 1987.
600 P.R. Viswanath / Progress in Aerospace Sciences 38 (2002) 571600
[57] Bechert DW. Experiments on three-dimensional riblets. [60] Launder BE, Li SP. The prediction of riblet behaviour with
Proceedings of International Conference on Turbulent a Low-Reynolds number ke model. Aeronaut J
Drag Reduction by Passive Means. London: Royal 1992;96:3515.
Aeronautical Society, 1987. [61] Pollard A, Savill AM, Tullis S, Wang X. Simulating
[58] Pollard A. Passive and active control of near- turbulent ow over thin element and at valley v-shaped
wall turbulence. Prog Aerosp Sci. 1997;33(11/12): riblets. AIAA J 1996;341(11):22618.
689708. [62] Pollard A. Riblets and other methods of controlling near-
[59] Bechert DW, Bruse M, Hage W. Experiments with three- wall turbulence. In: Choi KS, Prasad KK, Truong TV,
dimensional riblets as an idealised model of shark skin. editors. Emerging techniques in drag reduction. UK:
Expt Fluids 2000;28:40312. Mechanical Engineering Publications, 1996, p. 4576.