1. Introduction
Turbulent boundary layers (TBL) at high Reynolds numbers over smooth and rough walls are common in a variety of engineering applications and natural environments. The variation of the mean velocity (
$U$
) with wall-normal position (
$y$
) of a TBL with thickness
$\delta$
in the outer region of smooth/rough walls is usually represented with a log-wake composite profile
\begin{align} U^+_{smooth} & = \frac {1}{\kappa } ln(y^+) + B + \frac {\Pi }{\kappa }W(\eta ), \ \textrm {with}\, \eta = \frac {y}{\delta } ,\end{align}
\begin{align} U^+_{rough} &= \frac {1}{\kappa } ln\left (\frac {y-d}{y_0}\right ) + \frac {\Pi }{\kappa } W(\eta ) \ \textrm {with} \, \eta = \frac {y-d}{\delta -d} ,\end{align}
\begin{align} \textrm {where}\,\ W(\eta )& = 2 \eta ^2 (3- 2 \eta ) - \frac {1}{\pi } \eta ^2 (1- \eta )(1 - 2 \eta ) ,\end{align}
where
$U^+ = U/U_\tau$
is the mean velocity scaled with skin friction velocity (
$U_\tau = \sqrt {\tau _w/\rho }$
,
$\tau _w$
is the wall-shear stress and
$\rho$
is the fluid density),
$y^+ = yU_\tau /\nu$
is the inner scaled wall-normal position (
$\nu$
is the kinematic viscosity) and
$\kappa$
and
$B$
are the von Kármán constant and smooth wall intercept, typically set at 0.39 and 4.3, respectively (Marusic et al. Reference Marusic, Monty, Hultmark and Smits2013). For a rough wall,
$y_0$
is the roughness length, and
$d$
is the zero-plane displacement (or a virtual origin for the log region), and both are flow/surface-specific quantities. The outer region for a rough wall starts typically
$3-5$
representative roughness heights above the surface (Jiménez Reference Jiménez2004; Schultz & Flack Reference Schultz and Flack2007; Chung et al. Reference Chung, Hutchins, Schultz and Flack2021). The roughness length scale (
$y_0$
) is analogous to Nikuradse’s equivalent sandgrain roughness (
$k_s$
) and they are trivially related to each other (Chung et al. Reference Chung, Hutchins, Schultz and Flack2021). Finally,
$\Pi$
is Cole’s wake strength and
$W$
is the functional form for the wake. There are several implementations for the wake function (
$W$
), and in this work, we only consider the form in (1.3) provided by Lewkowicz (Reference Lewkowicz1982).
Townsend (Reference Townsend1956) first proposed outer-layer similarity between smooth and rough walls (at least for zero-pressure-gradient – ZPG flows) where the flow in the outer region is not different between different surface conditions and the influence of roughness is limited to the near wall region (which is below the outer region). This implies that turbulent motions (mean flow, turbulence statistics and even structures) may behave similarly regardless of surface conditions in the region outside of the immediate roughness layer, at sufficiently high Reynolds numbers (Schultz & Flack Reference Schultz and Flack2007). For the mean flow, this implies that the outer-wake region (i.e. the value of
$\Pi$
as well as the function
$W$
) is similar between smooth and rough walls. This is typically assessed with the mean velocity in deficit form as given by (1.4)
\begin{equation} \frac {U_{99} - U}{U_\tau } = -\frac {1}{\kappa } ln\left (\frac {y-d}{\delta - d}\right ) + \frac {\Pi }{\kappa } \left [2 - W\left (\frac {y-d}{\delta -d}\right )\right ] ,\end{equation}
where
$U_{99}$
is the boundary layer edge velocity.
Previous works have shown that mean profiles in deficit form do indeed collapse between smooth and rough walls for ZPG flows provided the representative roughness height (
$k$
) is small compared with the boundary layer thickness (
$\delta$
). Different studies have reported different thresholds for this ratio ranging from
$k/\delta$
= 0.02 to 0.1 (Jiménez Reference Jiménez2004; Castro Reference Castro2007). This threshold appears to depend on the type of roughness as well as the scale separation achieved in the flow (i.e. Reynolds number). The presence of outer-layer similarity together with the knowledge of
$y_0$
(roughness length of a given surface) allows us to develop models that can be used to calculate skin friction and other boundary layer parameters at higher Reynolds numbers (Castro Reference Castro2007; Monty et al. Reference Monty, Dogan, Hanson, Scardino, Ganapathisubramani and Hutchins2016).
Most realistic systems do not operate under ZPG conditions due to surface curvature or external flow effects. When studying pressure gradient (PG) flows, it is essential to distinguish between equilibrium and non-equilibrium conditions since that could affect the nature of outer-layer similarity as well as other characteristics of the boundary layer flow. Equilibrium flows are those in which the mean velocity profiles and flow statistics are invariant with the streamwise position. However, the only true equilibrium flow is that of a smooth wall sink flow (Townsend Reference Townsend1956; Rotta Reference Rotta1962). A near equilibrium boundary layer is defined by Marusic et al. (Reference Marusic, McKeon, Monkewitz, Nagib, Smits and Sreenivasan2010) as one where the mean velocity deficit exhibits self-similarity at a high enough Reynolds number. For a near equilibrium flow, the Clauser parameter
$\beta$
, defined as in (1.5), must be constant (Clauser Reference Clauser1954)
\begin{equation} \beta = \frac {\delta ^*}{\tau _w}\frac {dP}{dx} ,\end{equation}
where
$\beta$
is the PG parameter,
$\delta ^*$
is the displacement thickness,
$\tau _w$
is the wall-shear stress and
${\rm d}P/{\rm d}x$
is the streamwise PG. It remains unclear how non-equilibrium conditions, i.e. streamwise variation in
$\beta$
generated by PGs, affect both the flow over a rough surface and the wall similarity hypothesis.
Extensive studies have been performed on smooth wall flows under various PGs. However, work on rough walls (with comparable PGs) is more limited. Adverse pressure gradients (APG) have received more attention due to their association with flow separation and the resulting increase in drag. The most prominent effect of an APG is its impact on the wake region. Hot-wire measurements of Monty, Harun & Marusic (Reference Monty, Harun and Marusic2011) and laser doppler velocimetry measurements in the work of Aubertine & Eaton (Reference Aubertine and Eaton2005) over smooth walls show larger wake strengths (i.e. larger values of
$\Pi$
) with APG. Smooth wall Direct Numerical Simulation (DNS) results provided by Lee & Sung (Reference Lee and Sung2009) and validated by Monty et al. (Reference Monty, Harun and Marusic2011) showed that the wall-normal extent of the log region is limited under APG conditions. The skin friction coefficient has been found to decrease with APG in the works of Shin & Song (Reference Shin and Song2015b
) and Volino (Reference Volino2020); however, it is difficult to quantify these effects as these results were derived from the mean velocity profile. The work of Monty et al. (Reference Monty, Harun and Marusic2011) also showed this decrease using independent skin friction measurements but only for a smooth wall case.
The earliest experiments on rough wall boundary layers with APG were carried out by Perry & Joubert (Reference Perry and Joubert1963), who reported that the roughness function was independent of the PGthe flow experienced. However, this conclusion has been challenged by more recent works. The experiments of Pailhas, Touvet & Aupoix (Reference Pailhas, Touvet and Aupoix2008) found that an APG affected the value of
$k_s$
(or
$y_0$
). Particle image velocimetry of flows over ribs, carried out by Tsikata & Tachie (Reference Tsikata and Tachie2013), concluded that the combined effect increases
$k_s$
, while amplifying the wake and reducing the length of the log region. Tay, Kuhn & Tachie (Reference Tay, Kuhn and Tachie2009a
) found that
$k_s$
also increases with APG and that the effects on a TBL of both roughness and APG augment each other. Likewise, Tachie (Reference Tachie2007) concluded that the combined effect of roughness and APG was greater than that of roughness on its own, resulting in a larger roughness sublayer. Hot-wire measurements by Shin & Song (Reference Shin and Song2015b
) concluded differently that APGs reduce the effect of rough walls and reduce the skin friction compared with zero PG. Song & Eaton (Reference Song and Eaton2002), by means of laser Doppler anemometry, showed an earlier separation on a rough wall with APG compared with a smooth wall, supported by the work of Aubertine, Eaton & Song (Reference Aubertine, Eaton and Song2004). Turbulent boundary layers with favourable pressure gradients (FPG) have been studied with less attention in the literature. Over smooth walls, FPG boundary layers have been shown to increase the log layer length due to relaminarisation effects induced by the flow acceleration (Piomelli, Balaras & Pascarelli Reference Piomelli, Balaras and Pascarelli2000). As one may expect, the DNS simulations of Yuan & Piomelli (Reference Yuan and Piomelli2015) showed that TBLs under FPGs do not relaminarise in the case of rough wall flows. The work of Tay et al. (Reference Tay, Kuhn and Tachie2009b
) and Ghanadi & Djenidi (Reference Ghanadi and Djenidi2022) showed FPGs also result in thinner boundary layers with a smaller wake strength (i.e. smaller
$\Pi$
), if compared with a ZPG case. Skin friction has been shown to increase over smooth and rough walls by the works of Tay et al. (Reference Tay, Kuhn and Tachie2009
Reference Tay, Kuhn and Tachie
b
) and Shin & Song (Reference Shin and Song2015a
).
The aforementioned studies have investigated the effects of PGs over rough walls, however, often for a single PG type. Very few studies have examined the combined effect of FPGs and APGs. Limited work has been carried out on this subject over smooth walls. The work of Bobke et al. (Reference Bobke, Vinuesa, Örlü and Schlatter2017) looks at different non-equilibrium APG histories over a flat plate and an aerofoil surface. They found that different streamwise developments of
$\beta$
imply that even for equal flow states (
$Re_{\tau }$
and
$\beta$
) at the measurement location result in different velocity profiles and turbulence statistics. This phenomenon is due to the historical effects of the PGs on TBL development, which is also shown by Sanmiguel Vila et al. (Reference Sanmiguel Vila, Örlü, Vinuesa, Schlatter, Ianiro and Discetti2017). The effects of PG histories on flows over rough walls are poorly understood, yet it is crucial for predicting drag over rough surfaces and enhancing system efficiencies. The work of Fritsch et al. (Reference Fritsch, Vishwanathan, Roy, Lowe and Devenport2022) and Vishwanathan et al. (Reference Vishwanathan, Fritsch, Lowe and Devenport2023) as well as Volino & Schultz (Reference Volino and Schultz2023) demonstrated the variation in
$k_s$
under different PGs. These studies had contrasting conclusions with one suggesting that
$k_s$
(or
$y_0$
) is independent of PG histories (although the values of
$\beta$
explored were not very strong and the variance in
$k_s$
was much as 50 % across cases, but, without any specific trends with PG) while the other showed that
$k_s$
increases with FPG and decreases with APG. Vishwanathan et al. (Reference Vishwanathan, Fritsch, Lowe and Devenport2023) suggested that the variation in
$k_s$
is due to the choice of extent of log region during the fitting process, which is necessary to determine
$k_s$
(or
$\Delta U^+$
). Volino & Schultz (Reference Volino and Schultz2023) also indicated a dependence of
$k_s$
on
$k/\delta$
, suggesting that these results may reflect a lack of scale separation. In should be noted that all these studies have significant uncertainty in their results also due to indirect wall-shear stress measurements. Therefore, any fitting process and determination of parameters depend on the value of skin friction.
There is still a need for high-fidelity experimental data of boundary layers over smooth and rough walls experiencing different PG histories. This type of data can then be used for developing new predictive models for skin friction where the non-equilibrium effects can be captured. With access to high-fidelity data with sufficient scale separation, it may be possible to develop these models following the approaches of Perry, Marusic & Jones (Reference Perry, Marusic and Jones2002) or Castro (Reference Castro2007), where bulk boundary layer characteristics can be determined using momentum integral approaches. Alternately, an empirical relationship for skin friction can be developed (for example, Vinuesa et al. Reference Vinuesa, Örlü, Sanmiguel Vila, Ianiro, Discetti and Schlatter2017). Overall, the above review points to several open questions that need to further explored. These include:
$(i)$
Does the value of
$y_0$
(or
$k_s$
) depend on scale separation or PG history?
$(ii)$
Can we translate information from smooth wall flows with a given PG history to rough walls with similar histories?,
$(iii)$
Is it possible develop new prediction/correlation models that can infer boundary layer properties with limited measurements? and
$(iv)$
Would the history or strength of PGs applied to smooth/rough walls influence the applicability of these models?
In this study, we address some of these questions through detailed measurements of smooth and rough wall boundary layers in a region where the flow locally has zero PG, but, has experienced very different PG histories. Hot-wire, oil film interferometry and floating-element drag-balance measurements are carried out to gain new insights on the mean flow. Based on the observations from the data, we develop a correlation model that can be used to predict the local skin friction coefficient that includes history effects for smooth and rough walls. The paper is organised in the following sections. Section 2 discusses the experimental methods used and we present the mean velocity and skin friction data in § 3. The data are further reduced and analysed in the context of predictive models in § 4 with final conclusions and further recommendations in § 5.
2. Methodology
This section describes the experiments carried out and the PG histories imposed on smooth and rough wall TBLs.
2.1. Facility
(a) Experimental set-up used for hot-wire anemometry (HWA) measurements over both smooth and rough walls. Here,
$x/c=-1$
is located 5.28 m from the start of the test section. ① is the upstream Pitot tube from which
$U_0$
is set, ② shows the 16 pressure taps of the rough wall, ③ the NACA 0012 aerofoil of 1.25 m chord, ④ is the location of the drag balance used for skin friction measurements on the rough wall, ⑤ is the traverse to which ⑥ the HWA probe is mounted. (b) A 0. 25 m × 0.25 m section of the smooth wall constructed from aluminium sandwich panels. (c) A 0. 25 m × 0.25 m section of the rough wall constructed from plywood topped with 3 mm acrylic with 3 mm roughness mounted on top.
Experiments were carried out in the boundary layer wind tunnel at the University of Southampton, consisting of 5 sections of 2.4 m long, 1 m high and 1.2 m wide. A shallow ramp (
$\approx 5^\circ$
) is fitted at the start of the test section to remove the step up to the test surface. A three-dimensional turbulator strip of height 0.5 mm is placed at the top of the ramp to trip the boundary layer. The turbulence intensity in the free-stream region of the wind tunnel is 0.6 %. A simplified diagram of the tunnel and the set-up from 5 to 9.8 m is shown in figure 1, showing the main elements of the experimental set-up. The experiment uses a NACA0012 aerofoil of 1.25 m chord and 1.2 m width mounted on four actuators to adjust the aerofoil’s angle of attack. The aerofoil is mounted such that the leading edge is 6.53 m far from the inlet of the test section. The angle of attack and the wing’s distance from the wall
$(h)$
are adjusted for the different experimental campaigns. The height,
$h$
, is defined as the wall-normal distance from the wall to the quarter chord point. This set-up is similar to that of Fritsch et al. (Reference Fritsch, Vishwanathan, Roy, Lowe and Devenport2022) and Vishwanathan et al. (Reference Vishwanathan, Fritsch, Lowe and Devenport2023), however, we achieve stronger/longer PG histories due to the length of this wing as well as its location in the free stream. The position of the wing is set by measuring the height of all four corners of the wing above the tunnel floor. This limits the error of setting the wing using only the actuator encoders.
A Pitot tube is mounted one chord upstream of the leading edge of the aerofoil when set to
$0^\circ$
. The pressure difference is measured using a Furness FCO560 micromanometer. This Pitot sets
$U_\infty$
throughout this experiment. Temperature and pressure inside the tunnel at the at the exit of the contraction are measured using an RTD TST414 thermometer and Setra 278 barometric pressure transducer, respectively. These are sampled via the tunnel control system after every hot-wire point. The tunnel is kept at a constant temperature throughout an experimental run via the tunnel heat exchanger.
Smooth wall measurements were carried out using aluminium honeycomb sandwich panels, a section of which can be seen in figure 1(b). These are sized for each wind tunnel section to reduce the number of joints and have a thickness of 26.5 mm. For the section where measurements are being taken, the middle half of the tunnel is replaced with 10 mm thick safety glass. Firstly, this reduces the conduction effects of the hot-wire close to the tunnel floor. Secondly it allows optical access for wall-shear-stress measurements. The exception was for the upstream smooth wall, which used the aluminium wall. For rough wall measurements, the floor consists of 15 mm plywood topped with 3 mm PVC onto which an expanded metal mesh is mounted. The roughness runs from the start of the test section for approximately 10.5 m downstream. The metal mesh used has dimensions of 62 mm
$\,\times\,$
30 mm. The longest dimension is in the spanwise direction with a 3 mm height, resulting in an open area of 73 %.
2.2. Parameters
The main part of this study looks at five different angles of attack:
$-8^\circ$
,
$-4^\circ$
,
$0^\circ$
,
$4^\circ$
and
$8^\circ$
. To examine the effects of PGs, measurements are taken at one chord downstream of the trailing edge. For the smooth wall case, this distance is
$116.1\delta _0$
from the test section inlet, meanwhile, for the rough wall case,
$55.4\delta _0$
, where
$\delta _0$
is defined as the boundary layer thickness one chord upstream of the aerofoil. For the smooth wall, measurements are taken at free-stream velocities of 10, 20 and 30 m s−1. Rough wall measurements were taken at 10 m–30 m s−1 (in steps of 5 m s−1). These speeds corresponds to
$6.0\times 10^6 \lt Re_L \lt 19.6\times 10^6$
. The quarter chord is kept at 0.5 m from the wind tunnel floor for all these cases. Rough wall measurements were taken with the quarter chord at 0.4 m for
$-10^\circ$
,
$-8^\circ$
and
$-4^\circ$
at 20, 25 and 30 m s−1. Smooth wall data were only taken at this height for the
$-10^\circ$
and
$-8^\circ$
cases. The reason for this was to provide a greater range of PG histories and strengths. Measurements were also taken for a height of 0.5 m for
$-8^\circ$
and
$8^\circ$
one chord upstream of the aerofoil for various speeds mentioned above. The measurement station is
$67.9\delta _0$
from the inlet of the test section for the smooth wall case and
$32.4\delta _0$
for the rough wall case, with
${3.3\times 10^6} \lt Re_L \lt {10.4\times 10^6}$
. The ZPG measurements were taken for the rough wall from 15 m s−1 to 35 m s−1 (
${8.9\times 10^6} \lt Re_L \lt {21.0\times 10^6}$
) in steps of 5 m s−1. Meanwhile, for the smooth wall, the skin friction and velocity profiles data from Wangsawijaya, Jaiswal & Ganapathisubramani (Reference Wangsawijaya, Jaiswal and Ganapathisubramani2023) and Aguiar Ferreira, Costa & Ganapathisubramani (Reference Aguiar Ferreira, Costa and Ganapathisubramani2024) are used. This data are taken at the same measurement station as the other data in this work. Aguiar Ferreira et al. (Reference Aguiar Ferreira, Costa and Ganapathisubramani2024) also has direct skin friction measurements using oil film interferometry (OFI). When plotting different velocities the transparency is reduced when plotting different free-stream velocities (i.e. increasing Reynolds number). We note that data here capture an extended range of Reynolds numbers, ensuring sufficient scale separation.
2.3. Pressure distribution
Pressure taps are fitted to the wind tunnel floor to measure the wall pressure. For the smooth wall, twenty tubes with an inner diameter of 0.6 mm are fitted to the floor space 0.24 m apart. For rough wall measurements, sixteen pressure taps of 0.5 mm inner diameter were used, spaced approximately 0.265 m apart. Panel method simulations were carried out, suggesting that the upstream and downstream influence of the aerofoil extended one chord. As a result, the taps were placed in this region. This configuration of pressure taps is shown in figure 1(a) by the small vertical lines on the underside of the tunnel floor for the rough wall. The mean pressure distribution was recorded using a 64-channel ZOC33/64 Px pressure scanner. The pressure difference was taken referenced to the atmospheric pressure. The pressure data were sampled at multiple points during the hot-wire sweep for each PG case, temperature and pressure data are taken simultaneously.
Mean PG for both smooth and rough walls with respect to
$x/c$
. For the rough wall cases at
$h$
= 0.5 m the following symbols are used;
$-8^\circ$
: 
,
$-4^\circ$
: 
,
$0^\circ$
: 
,
$4^\circ$
: 
and
$8^\circ$
: 
. For
$h$
= 0.4 m the rough wall cases are denoted by;
$-10^\circ$
: 
,
$-8^\circ$
: 
and
$-4^\circ$
: 
. The ZPG case is given by 
. The smooth wall is shown by the same symbols and colours, however, they are left unfilled.
Figure 2 shows the PG history,
$dC_p/d(x/c)$
for the different cases tested during this experiment. The pressure coefficient here is given by
$C_p = (P-P_s) / (0.5 \rho U_\infty ^2)$
; here,
$P_s$
is taken to be the static pressure at
$x/c = -1$
. Further details on the pressure coefficient distributions can be found in Appendix A. For the 0.5 m cases, there are two distinct history types. The first are those that have an FPG followed by an APG (
$-8^\circ$
,
$-4^\circ$
,
$0^\circ$
). Secondly, those with an APG followed by an FPG (
$4^\circ$
$8^\circ$
). All 0.4 m cases fall into the first category. They have greater strength than the 0.5 m cases due to the aerofoil being closer to the wall. The first group of cases will be called APG cases throughout this work. While the second group will be the FPG cases. This is because it is assumed that the PG type experienced second will be more dominant in the resulting boundary layer one chord downstream.
The PG histories show good agreement between there smooth and rough wall cases. There is some variation in the PG around the hot-wire measurement station. In some cases, the slight FPG comes from the acceleration around the hot-wire traverse. In other cases, pressure distribution was taken when the traverse was removed. The effect of the traverse on the different boundary layers is assumed to be minimal and equal in all cases. Measurements were also taken in a nominal ZPG case, with the wing removed from the wind tunnel. This data were seen to have a slight FPG. The flow accelerates due to the boundary layer growing and the tunnel being a fixed cross-section.
The cases presented are non-equilibrium pressure distributions since
$\beta = (\delta ^* / \tau _0) \cdot ({\rm d}P/{\rm d}x)$
is not constant. The PG histories confirm the prediction of the panel method that the influence of the aerofoil extends one chord upstream and downstream of the aerofoil. The boundary layer one chord upstream should be the same since they have had the same PG history. The hot-wire is, therefore, placed one chord downstream of the aerofoil so that the local PG is the same for all cases. Thus, any difference can be said to be due to the upstream PG history.
2.4. Wallshear stress
When studying TBLs, an important quantity is the skin friction; however, the magnitude of these forces is of the order of grams, making it difficult to take direct measurements. The OFI is used to measure skin friction directly for smooth wall measurements. Silicon-based oil was used to generate the fringe patterns and imaged using a Lavision ImagerProLX 16MP camera, onto which a Sigma 105 mm F2.8 lens is mounted. A Phillips 35W SOX-E bulb provides a single wavelength of light. The wall-shear stress is calculated based on the thinning rate of the oil, calculated from the rate of change of the fringe pattern spacing.
Oil film interferometry is not possible for rough wall cases as it does not provide the true measure of skin friction (in addition to being impossible to implement in the case of current roughness). Therefore, a 0.20
$\,\times\,$
0.20 m drag balance is located in the tunnel floor with the centre of the balance located 8.8m downstream from the start of the test section along the centreline. The balance contains a floating element consisting of a flat plate mounted on a floating element hung from the outer casing with four thin flexures in the corners. The smooth, flat metal plate sits level with the outer casing and wind tunnel floor. Upon this flat plate, the roughness element patch is attached. An electromagnet and distance sensor with a range of
$\pm 250\,\unicode{x03BC} {\mathrm m}$
keeps the displacement at zero when a force is applied. The voltage applied to the electromagnet to maintain zero displacements; for more details, see Aguiar Ferreira et al. (Reference Aguiar Ferreira, Costa and Ganapathisubramani2024). Calibration is performed using a known load from 0 to 45 g, providing a calibration that can be used to convert the measured voltage during experiments into a drag force.
2.5. Hot-wire anemometry
Hot-wire anemometry is used to acquire velocity profiles from measurements at a single location. A vertical traverse is mounted from the roof as shown in figure 1(a) onto which a Dantec 55H21 probe holder is mounted. The traverse is located so that the hot-wire sits 9.0 m downstream of the start of the test section. The offset between the hot-wire and drag balance is 0.1 m to prevent any minor interference from the balance. A single in-house wire probe similar to the Dantec 55P05 probe is used. Consisting of which a
$5\,\unicode{x03BC}$
m tungsten wire is soldered; this is then coated with copper, leaving a sensing length of 1 mm. This results in a length-to-diameter ratio of 200, meeting the requirement in Ligrani & Bradshaw (Reference Ligrani and Bradshaw1987) that 1/d should not be less than 200 to prevent the conduction of the supports affecting the result. The dimensionless wire length
$l^+$
given by
$(l U_\tau / \nu )$
varies in the range
$21 \lt l^+ \lt 74$
for the smooth wall and in the range
$33 \lt l^+ \lt 140$
for the rough wall cases. In order to measure the upstream boundary layers, the traverse is moved so that the hot-wire probe is mounted at 5.3 m from the test section start. The rest of the set-up remains the same as for the downstream cases.
The overheat ratio was set to 0.8 throughout the experiment. The sampling time (T) varies from case to case, ensuring at least 19,000 boundary layer turnover times (
$TU/\delta$
). The signal is read using an NI USB-6 212 16-bit DAQ and is also used to read the Pitot tube micromanometers. The wire’s initial position is calibrated using a microscope camera. The wire is then moved towards the wall to set the initial position. The voltages are adjusted based on deviations from the initial calibration temperature to account for any temperature drift throughout a run. Calibration is carried out with the wing at
$0^\circ$
at a height of around 0.55 m from the floor to ensure the wing does not affect the calibration process. The probe is mounted half-way between the boundary layer and the wing. A second Pitot tube is mounted on the traverse at the same height as the hot-wire probe for calibration purposes. Calibration fitting uses a fourth-order polynomial to convert voltages to velocities.
3. Results
This section examines the mean boundary layer velocity profiles from both upstream and downstream of the aerofoil. Direct wall-shear stress measurements are also used. The section starts by examining the upstream boundary layer before it experiences the PG history and, subsequently, the downstream location after the flow has experienced different histories of PGs.
3.1. Incoming flow
The pressure distributions in § 2.3 showed a ZPG region approximately one chord upstream of the leading edge. Therefore, it is expected that the incoming boundary layer, one chord upstream of the leading edge, is invariant to the angle of attack.
Mean velocity profiles at similar
$Re_\theta$
taken at a hot-wire location of 5.3 m from the test section start for
$-8^\circ$
and
$8^\circ$
with the quarter chord at a height of 0.5 m for (a) smooth wall at 30 m s−1, (b) rough wall at 10 m s−1. Symbols and colours are as per figure 2.
The TBL for a smooth wall case, shown in figure 3(a), exhibits invariance to the angle of attack. Table 1 shows that the value of the boundary layer thickness (
$\delta$
) is invariant across all the tested free-stream speeds. The flow over the tested rough wall, as shown in figure 3(b), exhibits similar results, with the profiles collapsing between two angles of attack (−8
$^\circ$
and 8
$^\circ$
). Table 1 shows that the boundary layer thickness remains invariant even in the rough wall case. There is a slight variation between different cases; however,
$\delta$
is approximately 0.16 m in both cases tested. Upstream testing was conducted at two angles of attack: −8
$^\circ$
and 8
$^\circ$
. Since no dependence was observed at these extreme angles, it can be reasonably assumed that intermediate angles will show no dependence. This result validates the assumption made during the planning stage that the relevant PG history extends from one chord length upstream to one chord length downstream.
Summary of key boundary layer properties for two angles of attack one chord upstream of the aerofoil. Surface given as SW for smooth wall and RW for rough wall.
3.2. Mean flow – after experiencing pressure-gradient history
(a) Cut down version of figure 2 repeated to aid interpretation showing
${\rm d}Cp/{\rm d}(x/c)$
variation for rough wall cases, (b) Inner scaled velocity profiles at
$Re_\tau \approx 6800-8300$
for both smooth and rough wall cases at 0.5 m. The dashed black line shows the log region from (1.1). (c) Rough wall velocity profiles for 20 m s−1 for the 0.4 m, 0.5 m and ZPG cases. In both plots,
$d$
is the zero-plane displacement, which for a smooth wall is zero. The
$x$
axis is scaled using
$y_0$
, this results in the collapse of the log region of the profiles. Symbols and colours are as per figure 2.
The velocity profiles of the rough and smooth wall TBLs are inner scaled using the directly measured friction velocity. Figure 4(b) shows the inner scaled profiles at
$Re_\tau \approx 8000$
, i.e. free-stream speeds of
$30$
m s−1 for the smooth wall and
$10$
m s−1 for the rough wall cases. Table 2 gives the variation of the main boundary layer parameters used throughout the investigation. For the smooth wall case,
$Re_\tau$
ranges from 8310 for the
$-8^\circ$
to 6900 for the
$8^\circ$
. The variation in
$Re_\tau$
is much lower for the rough wall, with values varying from 7830to 8330 for the
$-8^\circ$
and
$4^\circ$
, respectively.
Summary of hot-wire data taken 9.03 m from the inlet of the wind tunnel for different PG histories.
The boundary layer profiles of the flow over the smooth wall collapse into the log region. It is seen that any case with APG just upstream of the measurement location results in a larger wake region (i.e. larger
$\Pi$
) and earlier deviation from the log law region. As the angle of attack increases, the resulting PG just upstream of the measurement location becomes more favourable, the wake becomes smaller (i.e. smaller
$\Pi$
) and the log region extends further away from the wall. The rough wall cases show the same wake and log region trends. In figure 4(b), there is a clear downward shift of the profiles from the smooth wall cases due to the roughness effects. As explained in § 1, this is because of the extra momentum loss and increased drag which depends on the type of roughened surface.
The roughness length scale chosen throughout this work is
$y_0$
. It was chosen since all the flow measurements were taken within the fully rough regime, as was shown by the skin friction measurements. In (1.2), the two unknowns in the log region are
$d$
(the zero-plane displacement) and
$y_0$
(the roughness length scale). Using the measurement value of skin friction, we first fit the zero-plane displacement while ensuring
$d$
should be less than the roughness height,
$k$
(Castro & Vanderwel (Reference Castro and Vanderwel2021)). The value of
$d$
is fitted using the diagnostic function,
$\Xi = ({1}/{U_\tau }) ({{\rm d}U}/{{\rm d}y}) (y - d)$
, which is equal to
$1/\kappa$
in the log region. The value of d is chosen to give the longest log region possible within the acceptable error range. The acceptable error range is defined such that the average deviation of the points chosen to be in the log region is less than
$\pm 5\, \%$
from
$1/\kappa$
. The resultant values of zero-plane displacement are shown in table 3. The results show APG just upstream of the measurement location reduces the value of the zero-plane displacement, and FPG increases it. For the strong APG cases, the value of zero-plane displacement is approximately zero. This suggests that the range of
$d$
values chosen as bounds may be limiting and that
$d$
could well be negative; however, this would be inconsistent with the recommendations from previous work. Further work may be necessary to examine the influence of PGs on this quantity. Since most of our work involves high Reynolds numbers, the exact choice of
$d$
has minimal impact on the results in the following sections, so we opt to leave it unchanged. The value of
$d_{ZPG}$
was found to be 0.46
$k$
and this is consistent with the work of Squire et al. (Reference Squire, Morrill-Winter, Hutchins, Schultz, Klewicki and Marusic2016), who, for ZPG flows, suggested choosing
$d$
as
$k/2$
.
Values of
$d / d_{ZPG}$
and
$y_0 / y_{0_{ZPG}}$
for different PG histories with values of
$d_{ZPG} = 0.00137m$
and
$y_{0_{ZPG}} = 0.000462m$
.
The method for finding
$d$
defines the bounds of the log region where the error remains within an acceptable range. Within this region, we obtain
$y_0$
as an offset from the smooth wall. The results are shown in table 3. While for the zero-plane displacement, there is a trend shown with PG history, for the roughness length scale, there is no clear trend. The
$-8^\circ$
case shows a small increase in the roughness length scale compared with the ZPG value. However, the cases at 0.4 m have roughness length scale values lower than the ZPG case. The absence of a clear trend with the PG history, along with the minimal variation in values, suggests that
$y_0$
is unaffected by the flow history. Any differences are attributed to the fitting process and the selected boundary layer region. Especially for high APG cases just upstream of the measurement location where the log region is small. This means that some of the wake region will likely be fitted to the log region, thus affecting the value of
$y_0$
. Regardless, the maximum deviation in
$y_0$
across the different cases is less than 20 % and, in fact, less than 10 % for the majority of cases examined here. This is much smaller than the deviations reported in Vishwanathan et al. (Reference Vishwanathan, Fritsch, Lowe and Devenport2023) and is presumably because of the scale separation that was achieved in this study where a considerable log region can be identified across all profiles. Moreover, an independent measure of
$U_\tau$
limits the uncertainty in fitting leading to better estimates of
$y_0$
.
The inner scale velocity profiles as a function
$(y-d)/y_0$
is shown in figure 4(c). This scaling provides a good collapse for all velocity profiles in the log region. The scatter between the different cases in the log region is reduced compared with the rough wall cases in figure 4(b). All profiles also exhibit a clear wake region that changes with the nature of the PG just upstream of the measurement location. However, it should be noted that this local wake is an integral effect of the entire PG history experienced by the flow. It is shown that the profile of the ZPG case has a wake between those at
$0^\circ$
and
$4^\circ$
at 0.5 m. This is expected because the type of PGs reverses order between these two cases. It can be seen that the
$-4^\circ$
at 0.4 m causes a larger wake than at 0.5 m but it is smaller than that is seen at
$-8^\circ$
case at 0.5 m. This pattern in the velocity profiles follows the same order seen in figure 2 between different cases. There is clearly a complex relationship between the wake profile and the imposed PG history. Be that as it may, the correlation between
$C_f$
and wake is consistent regardless of the history for both smooth and rough walls.
(a) Cut down version of figure 2 repeated to aid interpretation showing
${\rm d}Cp/{\rm d}(x/c)$
for
$-8^\circ$
,
$0^\circ$
and
$8^\circ$
at 0.5 m for the rough wall, (b) comparison of the velocity deficit profiles for
$-8^\circ$
,
$0^\circ$
and
$8^\circ$
at 0.5 m for both smooth and rough walls. (c) Shows the variation in velocity deficit profile for rough wall with Reynolds number for
$-8^\circ$
at a height of 0.5 m for 10, 20 and 30 m s−1. In both plots,
$d$
is the zero-plane displacement, which is zero for a smooth wall. Symbols and colours are as per figure 2.
The velocity deficit profiles enable the examination of the outer wake in more detail. The results for three angles of attack are shown in figure 5(b). For the
$-8^\circ$
case, there is a good collapse of the profiles between the smooth and rough wall cases. This would suggest that the integral effect of the PG and roughness on the outer region is similar to that of the smooth wall for this combination of PG history. However, no collapse occurs as the PG becomes more favourable immediately upstream of the measurement location. There is no outer-layer similarity because the boundary layer growth of the rough wall is larger than that of the smooth wall. Therefore, the integral effects of roughness and PG between the smooth and rough walls are not consistent. As a result, the rough wall in the presence of an FPG just upstream of the measurement location has smaller wake strengths (lower
$\Pi$
) than that of a smooth wall. Therefore, it may not be trivial to have information for a smooth wall with a given PG history (even at similar
$Re_\tau$
and identical local
$\beta$
) and use that to infer properties of a rough wall. As suggested by Volino & Schultz (Reference Volino and Schultz2023) and Vishwanathan et al. (Reference Vishwanathan, Fritsch, Lowe and Devenport2023), it appears important to match the
$\beta$
history to obtain complete similarity but that is almost impossible to devise in experiments (since
$\beta$
is an output while
${\mathrm d}P/{\mathrm d}x$
is the only input). Therefore, we need new relationships that will allow us to infer information about these flows based on local measurements.
Figure 5(c) shows that deficit profiles (for
$-8^\circ$
) collapse across different Reynolds numbers, and similar trends are observed for the other angles of attacks across smooth and rough wall cases. Based on this observation, the wake parameter,
$\Pi$
for each case is calculated using all available velocities for a given angle of attack. The fitting is carried out using (1.4), which only depends on directly measured values. The results of this fit are seen in table 2 in the column labelled
$\Pi$
. The values of
$\Pi$
obtained with the fitting process confirm that TBLs under APG just upstream of the measurement location have larger wake strengths compared with ZPG flows. In contrast, FPGs reduce the wake strength. As shown in the deficit profiles, the wake values of the APG cases are similar. The variation in
$\Pi$
shows some interesting trends. Firstly, for 0.5 m cases, the smooth wall wake strength is always greater than the rough wall. This is despite the similar PG histories shown in figure 2. As explained previously, this is due to the difference in boundary layer thicknesses and the resulting acceleration of the flow. Further evidence of this is that for the ZPG TBLs, the wake strength is much higher for the smooth wall, suggesting an FPG effect.
The trends observed for the case where the wing is mounted at a distance of 0.5 m far from the wall do not occur at 0.4 m. With the aerofoil mounted closer to the wall, and so with stronger APG conditions, the smooth wall TBL exhibits a lower wake strength if compared with the rough wall velocity profiles. As shown in figure 2, the match between smooth and rough wall cases worsens as the PG strength increases. However, the difference suggests that the smooth wall has stronger peak PGs. Therefore, one might expect the smooth wall to have a larger wake due to the APG; however, this is not what is observed. One possible explanation is that a thicker boundary layer is more susceptible to APG rather than FPG and thus results in stronger wake strength.
3.3. Skin friction – after experiencing pressure-gradient history
Most of the previous experimental studies on PG effects on TBLs over both smooth and rough walls had inferred skin friction from the velocity profiles. These methods introduce uncertainty into the measurements. Therefore, in this presented work, we aim to improve the experimental investigation by directly measuring the wall-shear stress, as outlined in § 2.4. The skin friction coefficient for the smooth wall is shown in figure 6(a). As expected for the ZPG smooth wall, the skin friction coefficient reduces as the Reynolds number increases (Schultz & Flack (Reference Schultz and Flack2013)). The
$8^\circ$
angle exhibits the strongest FPG just upstream of the measurement location and therefore exhibits the highest skin friction. Conversely, the
$-8^\circ$
angle is characterised by the strongest APG just upstream of the measurement location and hence the lowest skin friction. The other angles are arranged in order of increasing angle of attack between these two cases. The ZPG case lies between the
$0^\circ$
and
$4^\circ$
cases. From figure 2 this might be expected since these are the mildest two cases.
$0^\circ$
case experiences a mild APG upstream of the hot-wire while the
$4^\circ$
experiences an FPG region. Therefore, it makes sense that the ZPG cases fit between these two cases. These results indicate that the immediate upstream PG is more critical than those further upstream. This would indicate that any model that includes history effects should account for this variation in importance.
Skin friction coefficient one chord downstream of the trailing edge of the aerofoil for both 0.4 m and 0.5 m cases. (a) Skin friction coefficient for a smooth wall and (b) skin friction coefficient for a rough wall. Symbols and colours are as per figure 2.
The rough wall skin friction coefficients, in figure 6(b), do not depend on
$Re_x$
, meaning that the flow is in the fully rough regime. For the smooth wall, the average range, defined as the average of the range-to-mean ratio was 12 %, while for the rough wall, it was nearly constant at 2.4 %. The order of cases observed with the smooth wall is replicated with the rough wall at the 0.5 m height for the angles tested. The ZPG case follows this trend, positioning between the
$0^\circ$
and
$4^\circ$
cases. The rough wall cases at 0.4 m show similar trends with decreasing skin friction as the angle of attack becomes increasingly negative. As expected, the
$-4^\circ$
and
$-8^\circ$
cases at 0.4 m have lower skin friction than the equivalents at 0.5 m since the aerofoil is closer to the surface and thus has a stronger PG history. The local skin friction measurement clearly retains the history of the PG type and strength. Furthermore, the PG type in the second half of the domain is more dominant in overall local skin friction than the PG type in the first half of the domain.
3.4. Skin friction estimation from mean flow
Direct skin friction measurements are rare, with a majority of previous studies relying on the mean velocity profile to estimate the friction velocity and, therefore, the skin friction coefficient. With the data presented in the previous sections, it is possible to discern the difference between skin friction estimation using mean flow and direct skin friction measurements. The log law fitting method for the smooth wall is a simple problem using log region in (1.1), where the only unknown from the raw data are
$U_\tau$
. For the rough wall cases, we will use the log region in (1.2) assuming
$y_{0_{ZPG}}$
to be the
$y_0$
for all cases due to the very small variation. Therefore, when fitting the log law to this section, the zero-plane displacement was fixed at half the roughness height (0.5
$k$
= 1.5 mm), following the method used by Squire et al. (Reference Squire, Morrill-Winter, Hutchins, Schultz, Klewicki and Marusic2016). This approach leaves only one unknown parameter,
$U_\tau$
, which can be determined through curve fitting. Figure 7(a) shows the diagnostic function, defined as
$\Xi = (y-d) \cdot ({\rm d}U^+/{\rm d}y)$
. The black dashed line represents
$1/\kappa$
, as
$\Xi$
is equal to this value in the log region, as indicated by (1.1) and (1.2). This allows the extent of the log law to be assessed by examining deviations from
$1/\kappa$
. It can be seen that the smooth wall has an earlier departure from the log than the rough wall cases. Therefore, for the smooth wall cases, the log law is fitted with a minimum of seven points up to
$0.15\delta$
, while for the rough wall cases, this is extended up to
$0.2\delta$
.
(a) Diagnostic plot showing
$\Xi = (y-d) \cdot ({\mathrm d}U^+/{\mathrm d}y)$
for
$-8^\circ$
,
$0^\circ$
and
$8^\circ$
, both smooth (30 m s−1) and rough wall (10 m s−1) are shown at matched
$Re_\tau \approx 6800-8300$
. The black dashed line shows
$1/\kappa$
. (b) Comparison of
$U_\tau /U_{99}$
from log law fitting vs
$U_\tau /U_{99}$
from direct measurement techniques for both smooth and rough walls. The black dashed line is that of
$y = x$
, which would be true for a perfect prediction from in-direct methods. Symbols and colours are as per figure 2.
Overall, figure 7(b) shows good agreement between the direct measurement techniques and the predicted value from the mean flow profiles. For the smooth wall cases, the error varies between 0.2 % and 4.2 %, while for rough wall cases, the error is between 0.6 % and 4.8 %. The largest percentage error for the smooth wall cases is seen for the 10 m s−1 cases. This was expected since the values of
$U_\tau$
are the smallest. The largest error for the rough wall cases is seen in the FPG cases where the log law fitting method underpredicts the skin friction. These boundary layers have the longest log region and highest skin friction. Despite having a longer log region the error can be attributed to the fitting process. For consistency, the log region is assumed to end at
$0.2\delta$
; however, as seen in figure 7(a), the log region can be said to extend well into the outer region of the flow. Smooth wall data are expected to show good agreement since there is only one unknown in the fitting problem. For rough walls, the scatter is minimal since the
$y_0$
is assumed a priori for a given surface. Otherwise, the scatter would be significantly greater as both
$y_0$
and
$U_\tau$
would need to be fitted simultaneously and are interdependent. We note that the work here shows that fitting methods do indeed work where the flow has locally reached a zero PG (i.e. relaxing flows). However, it is unclear if this is still the case when the flow is locally subjected to a PG and that requires further work.
4. Development of a correlation model for skin friction
Examination of mean flow characteristics indicated that the PG just upstream of the measurement location was important. The data showed the wake strength of smooth and rough walls with the same input PGs are not similar. However, the roughness length in the measurement location did not exhibit any history effects and that suggests that any changes previously reported would have been due to a lack of scale separation or direct skin friction measurements. Finally, the data also showed that there is a clear correlation between local skin friction and local wake strength regardless of PG history. Therefore, it may be possible to develop a correlation model for skin friction and wake strength with some further modelling assumptions. This is explored further in this section.
We take inspiration from work of Vinuesa et al. (Reference Vinuesa, Örlü, Sanmiguel Vila, Ianiro, Discetti and Schlatter2017) who showed that for APG flows, the local skin friction can be predicted based on the skin friction (
$C_f$
) and shape factor (
$H = \delta ^*/\theta$
) of ZPG flows, and the streamwise-averaged PG parameter (
$\overline {\beta }$
). They showed that for weak distribution of
$\beta$
(i.e. ranging from
$0\sim 2$
)
\begin{equation} C_f^{APG} = \frac {C_f^{ZPG}}{(H^{ZPG})^{\overline {\beta }/2}}\ \ \textrm {where}\ \ \overline {\beta } = \frac {1}{x_{DS} - x_{US}}\int _{x_{US}}^{x_{DS}} \beta (x) {\rm d}x .\end{equation}
Here, the
$\overline {\beta }$
is the streamwise-averaged PG parameter, averaged from an upstream streamwise location,
$x_{US}$
, up to a downstream streamwise location,
$x_{DS}$
(from
${Re_\theta }_{US}$
to
${Re_\theta }_{DS}$
). It is unclear if this
$\overline {\beta }$
will be feasible for flows that experience both APG and FPG (or vice versa) in succession as well as for flows with surface roughness. For experiments, it is often impracticable to obtain complete
$\beta$
history due to the need for many streamwise measurement stations (either for smooth or rough walls). Moreover, it is unclear if the streamwise average as proposed in (4.1) is sufficient to capture history effects. The previous section showed that the locations just upstream of a given point are more important than locations further upstream. Therefore, we need to revise the approach to get better skin friction models for arbitrary PG histories. To tackle this challenge, we evaluate the mean velocity profile relationships in (1.1) and (1.2) at
$y = \delta$
, to get
$U^+_{99}$
and this is directly related to
$C_f$
(
$=2/{U^+_{99}}^2$
)
\begin{equation} \sqrt {\frac {2}{C_f^{PG}}} - \sqrt {\frac {2}{C_f^{ZPG}}} = \frac {1}{\kappa } ln \left (\frac {y_0^{ZPG}}{y_0^{PG}}\right ) + \frac {2}{\kappa } (\Pi ^{PG} - \Pi ^{ZPG}) + \frac {1}{\kappa } ln\left (\frac {\delta _{PG}^+}{\delta _{ZPG}^+}\right ). \end{equation}
Here, the superscript
$PG$
refers to an arbitrary PG case and
$ZPG$
is the zero-PG case. Here,
$\delta ^+ = Re_\tau = U_\tau \delta /\nu$
is the friction-velocity-based Reynolds number, and the final term will be zero if we match the Reynolds numbers between PG and ZPG cases. This equation only depends on
$y_0$
and
$\Pi$
for matched
$Re_\tau$
cases. If
$y_0$
does not change with PG as established in § 3, then, the change in
$C_f$
is entirely due to changes in
$\Pi$
. It should be noted that the work of Castro (Reference Castro2007) developed similar correlations to obtain the skin friction of rough walls flows at ZPG conditions where he showed that
$C_f = f(\theta /y_0, H, \Pi )$
. In fact, the variation in local
$\Pi$
accounts for flows that do not satisfy outer-layer similarity. However, it is possible to interpret the relationship to be an effect of external PG history. In that case, the correlations in Castro (Reference Castro2007) are analogous to the relationship in (4.2). Both require knowledge of the local value of
$\Pi$
(in addition to the value of
$H$
and
$y_0$
) to determine the local skin friction.
The difference in skin friction from (4.2) is plotted against the true difference in skin friction obtained from direct measurements in figure 8(a) for all cases. It can be seen that there is an excellent agreement with all of the points lying along the diagonal line. Figure 8(b) shows the relative contribution of each term of (4.2) to the skin friction difference. As to be expected from the data presented in table 2, for a given free-stream speed, the
$\delta ^+$
term is negligible and is seen to only have a marginal contribution to the overall skin friction. The contribution of the
$y_0$
difference is also very small since any difference between them is negligible. The dominant contribution is from the
$\Pi$
term as shown in figure 8(b). Equation (4.2) therefore provides a solution to predict the skin friction increase due to an unknown PG history, if the wake strength (
$\Pi$
) is known.
(a) Predicted difference in skin friction from (4.2) against the measured skin friction difference for the rough wall from the drag balance for 15, 20, 25 and 30 m s−1 for all PGs histories. The value of
$\Pi$
is taken from fitting the velocity profile to (1.1) and (1.2). The black dashed line is that of
$y = x$
, which would be true for a perfect prediction. (b) Relative contribution of each term in (4.2) to the overall drag of the surface at 20 m s−1. Symbols and colours are as per figure 2.
We need to be able to predict
$\Pi$
from a known flow history of PGs in order to determine the skin friction. Following Perry et al. (Reference Perry, Marusic and Jones2002), the obvious parameter that can be used in this correlation is the PG parameter
$\beta$
. For our current problem, we cannot use
$\beta$
since we only have data at a single streamwise location and
$\beta$
at this location is zero (since local
${\rm d}P/{\rm d}x \approx 0$
). Moreover, it will not be possible for us to evaluate
$\overline {\beta }$
as proposed by Vinuesa et al. (Reference Vinuesa, Örlü, Sanmiguel Vila, Ianiro, Discetti and Schlatter2017) as the streamwise distribution of
$\delta ^*$
and
$\tau _w$
is not available either (and it will not be for most studies as they are both outputs for a given
${\rm d}P/{\rm d}x$
history). Therefore, we introduce a new parameter,
$\Delta \beta$
that can account for the PG history. This parameter is defined as
\begin{equation} \Delta \beta = \left (\frac {\delta ^*}{\tau _w}\right )_{DS} \left [\frac {1}{x_{DS}-x_{US}}\int _{x_{US}}^{x_{DS}}\left (\frac {{\rm d}P}{{\rm d}x}\right ) w(x) {\rm d}x\right ]\ \ \textrm {where}\ w(x) = \frac {x-x_{US}}{x_{DS}-x_{US}} ,\end{equation}
where
$(\delta ^*/\tau _w)_{DS}$
is the ratio of displacement thickness to the wall-shear stress at the downstream measurement location, the integral term within
$[\cdot ]$
that includes
${\rm d}P/{\rm d}x$
distribution and a weighting function
$w(x)$
is the weighted integral of streamwise PG history between the two streamwise locations. We hypothesise that local values of
$\delta ^*$
and
$\tau _w$
at the measurement location already have history effects incorporated in them, and it may not be necessary to include them in the integral. Since the PG history closer to the measurement station has a greater influence compared with the upstream regions, some weighting should be applied to obtain a weighted PG history. For simplicity, a linear weighting as shown in (4.3) is applied.
Difference between
$\Pi ^{PG}$
and
$\Pi ^{ZPG}$
for smooth wall as a function of
$\Delta \beta$
. Only the 20 m s−1 data are shown. Here, the
$\delta ^*$
is calculated from the estimated profile with the near wall based on the Musker profile (Musker Reference Musker1979) and outer wake as given in (1.1). The black dashed line is the best fit to the data. Symbols and colours are as per figure 2. https://cocalc.com/Cambridge/S0022112025003209/JFM-Notebooks/files/Figure9.
Using our experimental data for smooth walls, we develop a correlation between
$\Delta \beta$
and
$\Pi$
for the measurement location. Figure 9 shows
$\Pi ^{PG} - \Pi ^{ZPG}$
as a function of
$\Delta \beta$
just for the smooth wall case. Only the 20 m s−1 case is plotted here while the trend appears to hold for other free-stream speeds as well. The blacked dash line gives the best fit linear relationship given by (4.4)
\begin{equation} \Pi ^{PG} - \Pi ^{ZPG} = 0.94\Delta \beta .\end{equation}
This suggests that once the weighted integral of the PG history (which can be provided as an input) is known, we can infer the local values of
$\delta ^*$
,
$\tau _w$
and
$\Pi ^{PG}$
. This form is different to what previous works have found, such as Das (Reference Das1987), who showed that local
$\beta$
varies as
$\Pi ^2$
. However, this would not appear to be case for the flows considered in the current study as local
$\beta \approx 0$
. The work of Perry et al. (Reference Perry, Marusic and Jones2002) gave a different functional form using a theoretical relationship,
$\beta = 0.5 + A\Pi ^{4/3}$
based on the attached eddy hypothesis (where the
$A$
is a constant derived from data). This relationship under predicts the value of
$\Pi$
as it does not fully capture the non-equilibrium effects. In fact, Perry et al. (Reference Perry, Marusic and Jones2002) included history effects, especially strong streamwise changes in
$\beta$
, through a gradient parameter
$\zeta$
that captures
${\rm d}\Pi /{\rm d}x$
. This gradient parameter together with the evolution equations (momentum integral) can be used to predict the streamwise evolution, which can be further calibrated using experimental data. In the current work, the history effects are captured with
$\Delta \beta$
through integration of
${\rm d}P/{\rm d}x$
weighted by a function,
$w(x)$
, over a fixed streamwise distance where the PG effects are present. Changing this weighting function will result in an altered relationship between
$\Delta \beta$
and
$\Pi$
. For example, an error function (rather than linear) that goes from
$0$
at
$x_{US}$
to
$1$
at
$x_{DS}$
(with 0.5 at the midpoint between
$x_{US}$
and
$x_{DS}$
) lead to approximately half the slope (i.e.
$\Pi ^{PG} - \Pi ^{ZPG} = 0.45 \Delta \beta$
). However, it does not take away from the nature of the relationship between the two quantities. Comparing the current approach with that of Perry et al. (Reference Perry, Marusic and Jones2002) is more involved, but the first steps towards this comparison is presented in Appendix B. Further work is required to reconcile the similarities/differences between (4.4) (and various different weighting functions) and the work of Perry et al. (Reference Perry, Marusic and Jones2002) which is beyond the scope of this study.
(a) Predicted value of
$C_f$
using minimisation function of equation 4.4 and the predicted fit of the velocity profile.Here,
$\delta ^*$
is provided as calculated from the hot-wire velocity profile. This is compared with the measured value of
$C_f$
with the black dashed line showing
$y=x$
, a perfect prediction. Data shown for the 20 m s−1 cases. (b) The predicted value of
$C_f$
using the minimisation function of (4.4) and the predicted fit of the velocity profile.Here,
$\delta ^*$
is calculated using the velocity profile in (1.2) where the value of
$\Pi$
is implicitly included. The black dashed line shows
$y=x$
, which would be true for a perfect prediction. Data shown for the 20 m s−1 cases. Symbols and colours are as per figure 2. https://cocalc.com/Cambridge/S0022112025003209/JFM-Notebooks/files/Figure10.
We can test the predictive capability of the above relationship on flow over rough walls experiencing arbitrary PGs as it is clear from figure 8(b) that most of the skin friction comes from the changes in the wake parameter, especially if roughness length is not altered by the PG history. Given the weighted integral of the upstream PG history (as depicted in (4.3)), we solve (4.2) and (4.4) simultaneously to obtain
$C_f^{PG}$
and
$\Pi ^{PG}$
provided
$C_f^{ZPG}$
,
$y_0$
and
$\Pi ^{ZPG}$
are known. Solving these equations also requires an input of
$\delta ^*$
or the velocity profile data at the location of prediction. Figure 10(a) shows the prediction of skin friction of rough walls experiencing different PGs from this minimisation. The figure shows good agreement between prediction and data, demonstrating the suitability of the derived correlations for both smooth and rough walls. In fact, we can go a step further and include the calculation of
$\delta ^*$
as part of the minimisation process. In this case, the value of
$\delta ^*$
can be calculated from casting the mean velocity profile in (1.2) in the appropriate form for displacement thickness, and therefore it depends only on input values for ZPG case (at matched
$Re_\tau$
). Consequently,
$\delta ^*$
can be implicitly determined as a part of the minimisation. Figure 10(b) shows that the agreement between measured and predicted skin friction using this approach and it appears to be just as good as the former (where
$\delta ^*$
was given as an input). Overall, this shows the merit of the derived correlation-based method to predict skin friction of smooth and rough wall flows experiencing arbitrary PG histories. As such this method does have some limitations. Although the local values of
$\beta$
over the PG history can attain large values (preliminary estimates show
$\beta$
up to 5 in our flows), the range of integrated PG strengths is limited (i.e.
$0\lt \Delta \beta \lt 1.5$
). We note that this integrated value is still higher than local
$\beta$
achieved in some previous studies (Vishwanathan et al. Reference Vishwanathan, Fritsch, Lowe and Devenport2023). It is also unclear if the shape of the PG history is critical for this correlation. As such we are unable to validate this as most previous work do not have direct skin friction measurement or the scale separation required to use the proposed correlation. Exploring these effects should be considered in future works.
5. Conclusions
Hot-wire and skin friction measurements have been presented for smooth and rough wall TBLs, showing the effect of non-equilibrium PGs. Using a NACA0012 aerofoil, strong PG histories were created. The effects were investigated by measuring the velocity profiles one chord downstream of the trailing edge. Firstly, it was demonstrated that the PG is approximately zero one chord length upstream of the aerofoil. Consequently, the boundary layer one chord upstream remains invariant to the angle of attack. Therefore, any changes downstream of the wing are due to the different PG history. Several angles of attack of the wing produce different PG distributions along the streamwise direction. It was assumed that the PG type experienced closer to the measurement station was more dominant compared with the flow history further upstream. The direct skin friction measurements supported this conclusion. They showed that cases experiencing a FPG followed by an APG have lower skin friction compared with cases experiencing them in the reverse order.
The velocity profiles taken one chord downstream of the aerofoil for both smooth and rough walls are at approximately matched
$Re_\tau \approx 6800-8300$
. There is a clear downward shift in the profiles for the rough wall compared with the smooth wall, presumably due to the additional drag. There is also a larger wake for APG immediately upstream, while an upstream FPG is found to suppress the wake. A key result was that the variation in the roughness length scale
$y_0$
is not significant. The observed difference is assumed to be attributable to the fitting process and is thus considered invariant. Furthermore, it was seen that if
$y_0$
is known from ZPG measurements, it is possible to predict the skin friction for the rough wall within 5%. This error is comparable to the error observed in smooth wall log law fitting.
The mean flow results enabled an investigation into how the effect of PGs can be predicted in TBLs. We examined modelling the difference in skin friction between a PG case and a ZPG case for a given surface. Inspired by the works of Perry et al. (Reference Perry, Marusic and Jones2002); Castro (Reference Castro2007); and Vinuesa et al. (Reference Vinuesa, Örlü, Sanmiguel Vila, Ianiro, Discetti and Schlatter2017) we developed a new correlation between the local skin friction and wake strength parameter, which in turn depends on a weighted-averaged PG parameter,
$\Delta \beta$
. The correlation was first trained on smooth wall data and then applied to flow over rough walls for predictions. We showed that this correlation can predict the skin friction of flows over smooth and rough wall boundary layers with arbitrary PG history effects (given some input parameters from ZPG flows). Further work is necessary to understand the meaning of this new quantity and its ability to capture history effects. Work is also needed to capture the physical mechanisms (i.e. flow structure) that may be responsible for these observations.
Supplementary material
Computational Notebook files are available as supplementary material at https://doi.org/10.1017/jfm.2025.320 and online at https://cocalc.com/Cambridge/S0022112025003209/JFM-Notebooks.
Acknowledgements
We gratefully acknowledge the help of N. Agathangelou for taking the upstream smooth wall data as part of his undergraduate project.
Funding
We gratefully acknowledge the financial support from EPSRC (Grant Ref no: EP/W026090/1) and European Office for Airforce Research and Development (Grant No: FA8655-22-1-7163, Programme Manager: Dr D. Smith).
Declaration of interests
The authors report no conflict of interest.
Data availability statement
Data published in this article are available on the University of Southampton repository (DOI: 10.5258/SOTON/D3363).
Author contributions
T.P. designed the set-up, conducted experiments, post-processed HWA and balance data and wrote the first draft. M.V. conducted experiments, post-processed OFI data and edited drafts. P.J. conducted experiments and edited drafts. B.G. provided the ideas for the experiment, acquired funding, edited drafts and managed the project.
Appendix A. Further details on the pressure coefficient and its gradient
(a) Mean
$C_p$
pressure distribution history from one chord in front of the aerofoil to one chord behind the aerofoil for the five cases at
$h$
= 0.5 m for both smooth and rough walls. For the
$h$
= 0.4 m cases, the three rough and two smooth wall cases are shown. Solid lines are used for smooth wall datasets, while dashed lines are used for rough wall data. (b) Smooth wall
${{\rm d}C_p}/{{\rm d}(x/c)}$
for the
$-8^\circ$
and
$8^\circ$
cases for all Reynolds numbers. Symbols and colours are as per figure 2. As the Reynolds number increases, the opacity of the marker is increased.
The mean pressure distribution from all speeds is presented for smooth and rough walls at the five angles of attack for 0.5 m height in figure 11(a). Also presented are the cases with a quarter chord height of 0.4 m. It can be seen that both the smooth and rough wall datasets have very similar shapes. At 20 m s−1 the typical uncertainty in the
$C_p$
value is
$\pm 0.025$
. Furthermore, the maximum pressure coefficient is located at the same streamwise location. There are slight differences in the pressure coefficient between the rough and smooth walls. The smooth wall for strong negative angles of attack has higher magnitude peak values compared with the rough wall. The reason for this is likely due to the proximity of the taps to the roughness elements. This results in a lower pressure coefficient than for taps, which are further from roughness elements.
Figure 11(b) shows the PG for the
$-8^\circ$
and
$8^\circ$
cases for the smooth wall for different velocities. This shows an important result that the PG history is invariant to the Reynolds numbers, shown here for the most extreme cases. The small variations across the different speeds are due to minor differences in the boundary layer thicknesses, changing the tunnel’s effective cross-section. The same result is seen for the rough wall data with Reynolds number.
Appendix B. First steps towards reconciling model presented in Perry et al. (Reference Perry, Marusic and Jones2002) with the current work
The work of Perry et al. (Reference Perry, Marusic and Jones2002) presented a model which incorporates PG history effects into calculating
$\Pi$
. The key parameter for their work is
$\zeta$
, which is given by (B1). Using this parameter, they define a relationship for
$\beta$
, which is shown in (B2). This gives
$\beta = f(\Pi , {\rm d}\Pi /{\rm d}x)$
that, for many flows, presents a problem since knowing
${\rm d}\Pi /{\rm d}x$
is often not practical or possible to predict. It requires measurements at regular intervals with skin friction and boundary layer measurements at each station
\begin{equation} \zeta = \frac {U_{99}}{U_\tau } \delta \frac {{\rm d}\Pi }{{\rm d}x} ,\end{equation}
\begin{equation} \beta = \begin{cases} -0.5 + 1.2\Pi ^{\frac {4}{3}} + \zeta ^2 (1.10/\Pi ^2) & \text{if $\zeta \geqslant 0$},\\ -0.5 + 1.2\Pi ^{\frac {4}{3}} + \zeta (0.62 + 0.25 \Pi ) & \text{if $\zeta \lt 0$}.\\ \end{cases} \end{equation}
Using the definition of
$\beta$
it is possible to rewrite (B2) in the form of (4.3) as given in (B3). The left-hand side is the weighted integral of the PG and the right-hand side comes from (B2)
\begin{equation} \frac {1}{x_{DS}-x_{US}} \int ^{x_{DS}}_{x_{US}} \frac {{\rm d}P}{{\rm d}x} w(x) {\rm d}x = \frac {1}{x_{DS}-x_{US}} \int ^{x_{DS}}_{x_{US}} \frac {\tau _w}{\delta ^*}f\left (\Pi , \frac {{\rm d}\Pi }{{\rm d}x}\right ) w(x) {\rm d}x ,\end{equation}
where
$f(\Pi , {\rm d}\Pi /{\rm d}x)$
is given by (B2).
Comparison of variation
$\Pi ^{PG} - \Pi ^{ZPG}$
with
$\Delta \beta$
as given by (4.4) and variation in
$\Pi ^{PG} - \Pi ^{ZPG}$
with
$\beta$
from the model of Perry et al. (Reference Perry, Marusic and Jones2002). The black line is showing the fit of (4.4) and the labelled lines show variation of
$\Pi _{PG} - \Pi _{ZPG}$
for different
$\zeta$
values from (B2).
Figure 12 shows contours of
$\Pi$
and
$\beta$
from Perry et al. (Reference Perry, Marusic and Jones2002) for various values of
$\zeta$
. A flow with a given PG history will trace out a path in this space. However, the
$\Delta \beta$
term in the current study is a weighted integral of the PG term, and this will take on a value that can only be obtained through the integration of the function in this space. Figure 12 also shows the curve that relates
$\Delta \beta$
against
$\Pi$
from (4.4). It is clear that the relationship in (4.4) jumps from across different
$\zeta$
curves for different PG histories. This can also be seen from (B3) where
$\Delta \beta$
is indeed an integral across different curves. The effect of this space is, in fact, captured through the weighting function, and at this stage, linear weighting appears to capture the trends reasonably well. However, further work is required to tune the weighting function for different PG histories and to reconcile the similarities/differences with Perry et al. (Reference Perry, Marusic and Jones2002).
Appendix C. Comparison of current work with Castro (Reference Castro2007)
Castro (Reference Castro2007) provided a three parameter model for skin friction
$C_f = f( {\theta }/{y_0}, H, \Pi )$
over rough walls in ZPG conditions (where
$H$
is the shape factor of the boundary layer and
$\theta$
is the momentum thickness). This work showed that for a given value of
$\Pi$
and H the variation in
$C_f$
with
$\theta / y_0$
can be calculated using (C1). Alternately, it is also possible to obtain a variation of
$C_f$
with
$H$
following for different values of
$\Pi$
as shown below
\begin{eqnarray} \sqrt {\dfrac {2}{C_f}} & = -\dfrac {1}{\kappa } ln \left ( \dfrac {1}{H} \sqrt {\dfrac {C_f}{2}} \right ) + \dfrac {1}{\kappa } ln \left (\dfrac {\theta }{y_0}\right ) + \dfrac {2 \Pi }{\kappa } - \dfrac {1}{\kappa } ln\left (\dfrac {1 +\Pi }{\kappa }\right ), \end{eqnarray}
\begin{eqnarray} \sqrt {\dfrac {2}{C_f}} & = \left (\dfrac {H}{H-1}\right ) \left [\dfrac {2.009 + 3.018 \Pi + 1.486 \Pi ^2 }{\kappa (0.983 + \Pi )}\right ]. \end{eqnarray}
Three parameter model based on Castro (Reference Castro2007) where it is assumed
$C_f = f(\theta /y_0, H, \Pi )$
. (a) Variation in
$C_f$
with
$\theta /y_0$
curves show predicted
$C_f$
variation for each cases
$\Pi$
from (C1). (b) Variation in
$C_f$
with
$H$
curves show predicted
$C_f$
variation for each case
$\Pi$
from (C2). Symbols and colours are as per figure 2.
Note that the constants in (C2) are taken directly from Castro (Reference Castro2007) and they depend on the type of wake profile fitted to the velocity data. This is consistent with the use of Lewkowicz (Reference Lewkowicz1982) polynomial wake profile, which is also used in the current study.
This analysis can be extended for flows with PGs. Figure 13(a) and 13(b) shows the predictions of
$C_f$
as a function of
$\theta /y_0$
as well as
$C_f$
versus
$H$
for the different rough wall flows examined in this study. Lines of constant
$\Pi$
are shown in both figures and these lines are based on the value of
$\Pi$
from the new correlation developed in § 4 for different PG histories.
Both figures show very good agreement between trendlines for different
$\Pi$
and the skin friction measured. We also attempted to fit the wake function from Coles (Reference Coles1956) wake function which does not seem to affect the fitted value of
$\Pi$
. However, in this case, the coefficients in (C2) have to be altered. Regardless, this agreement shows that the method proposed in the current work is consistent with the results of Castro (Reference Castro2007), provided a suitable value of
$\Pi$
as determined from
$\Delta \beta$
is used.
Open access
