Assessment of aerodynamic roughness parameters of turbulent boundary layers over barnacle-covered surfaces

Full-scale drag penalty predictions of flows over rough walls require surface roughness characterisation from laboratory experiments or numerical simulations. In either approach, it is necessary to determine the so-called equivalent sand-grain roughness height (ks\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k_s$$\end{document}). There are several steps involved in determining this aerodynamic roughness lengthscale, but its procedure typically includes a combination of measurement of wall-shear stress (τw\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tau _w$$\end{document}) using direct or indirect methods as well as analysis of velocity profiles. Indirect methods usually rely on assumptions made about flow and its scaling including the validity of universal outer-layer similarity. However, the implications of the underlying assumptions involved in full-scale drag prediction are unclear. In this work, we carry out wind tunnel measurements over a realistic rough surface (from a fouled ship-hull) to evaluate the impact of different methods with an emphasis on using the outer-layer similarity hypothesis for full-scale drag predictions. Wall-shear stress is measured using an in-house floating-element drag balance (DB), and velocity profiles are obtained using particle image velocimetry (PIV), allowing the evaluation of ks\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k_s$$\end{document}, and the associated wake parameters through several methods. The aerodynamic roughness parameters hence obtained are used for full-scale drag penalty calculations. It is observed that the predicted drag penalty can vary by over 15%\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\%$$\end{document} among the different methods highlighting the care that should be taken when employing such methods.


Introduction
The skin friction drag is well-established as one the most valuable quantities characterising smooth-and rough-wall bounded flows, as it embodies the result of the conversion of the free stream momentum into shear force by means of the boundary layer.Considering the ever-increasing interest in reducing the carbon footprint and the environmental impact of fuel-based systems, accurate measurement and prediction of this quantity remain vital to a whole range of engineering applications (e.g.design of aerodynamically efficient vehicles).Skin friction is equally important for researchers as it influences strongly the fluid motion near the wall but is also crucial in understanding wall-turbulence dynamics due to its scaling capabilities for both the near-and outer-wall regions (Chung et al. 2021).
Regardless of whether the frictional drag is directly or indirectly measured, the use of similarity laws is required

Background
Above the roughness sublayer, the mean flow over a rough wall is sufficiently captured through the log-wake law, similar to a smooth wall.Therefore, the mean velocity profile in inner and outer forms can be, respectively, written as, where y is the wall-normal location whose origin is taken as the lowest roughness height (h min ) , U ∞ is the freestream velocity, U is the mean streamwise velocity, U = √ w ∕ is the skin-friction velocity ( w is the wall-shear-stress and is the density of the fluid) and the boundary layer thickness is defined as the wall-normal height where U ≈ 0.99U ∞ , is the fluid kinematic viscosity.In the above equations, the inner-normalised velocity deficit in the log region ΔU∕U ≡ ΔU + is known as the roughness function (Hama  1954) and depends on the nature of the surface roughness.d is a virtual origin (also referred to as the zero-plane displacement) commonly defined as the height at which the mean momentum sink occurs (Jackson 1981), and similarly depends on the roughness.Both ΔU + and d are equal to zero for a smooth wall.The logarithmic region is described by its slope , paired with its smooth-wall intercept B. However, it has been shown that these "constants" vary depending on the flow geometry with ranging from 0.37 to 0.41 between turbulent channel, boundary layer, and pipe flows (Nagib and Chauhan 2008).
Beyond the log layer, an overshoot of the velocity from the logarithmic behaviour occurs and is believed to emanate from a dual effect of the velocity jumps in the uniform momentum (1) zones as well as the jump across the turbulent/non-turbulent interface (Krug et al. 2017).This is described by the term W(y∕ ) known as the Coles wake function (Coles 1956), which is characterised by its wake strength parameter Π .There have been numerous studies in the literature that proposed analytical expressions to capture this region, and depending on the expression, the same flow yields different estimates of Π (see for example Jones et al. (2001); Guo et al. (2005); Nagib et al. (2007) for reference).In fact, for a given and B, Π in a zero- pressure gradient turbulent boundary layer is reported to be around 0.55, 0.65, and 0.45 for the expressions used by Hinze (1975), Jones et al. (2001) and Nagib et al. (2007), respectively.Therefore, comparing the wake strength parameter is only meaningful when the same expression is used.Moreover, for a given wake function, Π is also sensitive to pressure gradients, with Π increasing and decreasing with adverse and favourable pressure gradients, respectively.By adding equation 1 and 2, an expression for the skinfriction coefficient can be obtained, and is given by, with the superscript "+" denoting quantities non-dimensionalised by U and .Assuming and B universal constants, expression 3 clearly shows that C f = f + , ΔU + , Π , with + ≡ Re = U ∕ being the friction Reynolds number.Therefore, the accuracy in determining the constitutive quantities of C f can have serious implications on the pre- diction of skin friction and drag penalty at high Reynolds numbers.In this study, ΔU + and Π will be evaluated through different methods, allowing us to quantitatively assess their resulting predictions.

Outer-layer similarity
Over rough walls, outer-layer similarity (OLS) emerges as an extension of Townsend's wall-similarity hypothesis, which states that the influence of viscosity is limited to the near-wall region (Townsend 1976).In this scenario, roughness effects are hypothesised to be confined to the near-wall region, provided the roughness length scale (k -representative roughness height) is a small portion of (i.e.k∕ is small).In other words, a large separation between scales is required ( 1 ≪ k + ≪  + ).This means that in the wake region, the mean and the turbulence motions can be assumed to have a universal form between smooth-and rough walls.This is generally expressed in a functional form as This means that the role of roughness is limited to the adjustment of the friction velocity ( U ) and the boundary layer thickness ( ), while the mean and turbulence structure functions (f) and (g ij ) , respectively, remain unaffected.From the early work of Perry and Abell (1977) on roughened pipe flows, numerous studies have shown the presence of outer-layer similarity in rough-wall flows over various surface conditions (Andreopoulos and Bradshaw 1981;Acharya et al. 1986;Perry and Li 1990;Raupach 1992;Schultz and Flack 2005;Castro 2007; Wu and Christensen 2007;Squire et al. 2016).However, other researchers have also reported cases where rough-walls do not behave similarly to smooth-wall flows under certain circumstances (Krogstad et al. 1992;Krogstad and Antonia 1994;Volino et al. 2009;Chan et al. 2018;Chung et al. 2018;Medjnoun et al. 2018Medjnoun et al. , 2020;;Womack et al. 2022).
Intuitively, the relative height of the roughness ( k∕ with k being a representative roughness height) is suggested to be particularly important for a flow to follow outer-layer similarity, with several works proposing thresholds with a minimum of ∕k > 10 (Jimenez 2004; Schultz and Flack 2005).However, the results outlined by Flack et al. (2007) and later by Amir and Castro (2011) clearly indicated that outer-layer similarity can still be valid even at ∕k < 10 suggesting that a single critical roughness height beyond which outer-layer similarity fails does not necessarily exist.They argued that the outer flow is only gradually altered with increasing roughness height.Nonetheless, there is growing evidence from both recent experimental and numerical studies that suggest similarity remains strongly dependent not only on roughness amplitude but also on other parameters such as its phase and directionality (Chung et al. 2021), suggesting that universal similarity might not be possible for some types of surface roughness.Therefore, it is important to examine the implications of this hypothesis when drag penalty predictions are desired and provide corrections if necessary.

Roughness function
The skin-friction law in equation 3 can be written for a rough-and a smooth-wall ( ΔU + = 0) and subtracted from each other to determine the roughness function expression: (5) with the superscripts S and R identifying quantities measured over smooth and rough wall flow conditions, respectively.Under the current hypotheses, equation 5 shows that the roughness function is characterised by three main contributions.The term I) represents the roughness function contribution stemming from the change in frictional drag between surfaces, the term II) captures the contribution arising due to the difference in the frictional Reynolds numbers, whereas the term III) shows the roughness function contribution caused by the alteration in the wake strength parameter.
In equation 5, the function governing the wake W in the outer flow is assumed to be universal (i.e.universal form), however, the wake strength (i.e. the amplitude of the function) can be different between surfaces, and it is captured by its wake strength parameter.In fact, changes in the outer flow are usually reflected in different values of Π for differ- ent surface conditions.Therefore, an accurate assessment of the roughness function requires accurate knowledge of the different constitutive terms that characterise the roughness function.
If Π R = Π S , then there is universal outer-layer similarity.In this instance, ΔU + can still be determined by substituting the two different values of + in equation 5.This is identical to the procedure used by Granville (1987), and with some further manipulations, the integrated skin-friction coefficient C F over a plate of length L can be obtained.If also + R and + S are matched in the presence of universal outerlayer similarity, then ΔU + reduces to the difference in the frictional drag, and can be obtained from direct wall-shear stress measurements, and is expressed as: with the subscript, DB representing the quantity measured directly using the drag balance technique, employed in this study.The skin friction coefficient C S f for the smooth wall can be separately measured or obtained from the smoothwall expression (i.e.equation 3).
In the absence of direct skin-friction measurements, velocity profile data are required.It is a common practice for studies on rough-wall flows to infer the aerodynamic quantities from a fitting procedure of the logarithmic velocity data, or the total shear stress profile (Perry and Li 1990;Krogstad et al. 1992;Flack et al. 2005;Schultz and Flack 2005;Flack et al. 2007;Wu and Christensen 2007).Alternatively, the use of Townsend's similarity hypothesis enables estimating U from the outer flow, however, with the assumption that similarity does exist.It is only in some isolated facilities such as in towing tanks and a few wind tunnels where direct drag information has been used to determine ΔU + (Cheng and  Castro 2002; Schultz and Myers 2003; Krogstad and Efros  2010; Baars et al. 2016; Ferreira et al. 2018).Their limited (6) use is essentially due to their manufacturing cost as well as the delicate attention needed when operating such devices.However, while it is possible to use equation 6 to determine ΔU + under certain conditions (matched + and pres- ence of OLS), in practice, both conditions cannot be a-priori assumed without examining the flow.Hence, in the absence of a guarantee that outer-layer similarity can be valid for a given rough surface, velocity profile measurements are essential to measuring all relevant quantities and ensuring appropriate predictions are possible.When velocity information is available, the roughness function can be determined from the shift in the log region, relative to the smooth-wall profile, directly from the inner-scaled velocity distribution.In fact, the value of ΔU + hence determined represents the momentum loss or gain (in case of drag reduction) and is the correct estimate of the roughness function, provided there is an independent (direct) estimate of U .
Once ΔU + is determined from drag, velocity measure- ments, or both, the rough surface is exposed to different freestream speeds to obtain the variation of ΔU + as a func- tion of the aerodynamic roughness lengthscale.By evaluating ΔU + as a function of k + s , the equivalent sand-grain roughness height k s can be determined using the expression: where, C is an additive constant, determined empirically to be equal to 3.5 (Nikuradse 1933).Combining equation 7 together with 3, an expression for fully-rough flows can be written as: Equation 8 shows that k s depends on the boundary layer thickness (with being also a function of x for boundarylayers), the rough wake strength parameter Π R (with Π R = Π S if OLS is satisfied), the rough-wall intercept B R (with B R = B + C ), the logarithmic slope κ and the skin-friction coefficient C f .The quantity k s does not represent a geometri- cal lengthscale of the roughness but rather an aerodynamic property of the rough-wall flow investigated.Hence, it can only be determined with the assessment of the flow, and as such, equations 7 or 8 require the flow to be exposed to a range of Reynolds numbers to detect whether the data has reached an asymptotic behaviour.If the roughness function follows the trend in either equations 7 or 8, then the flow over the surface is considered to be "fully rough", and k s can subsequently be used to predict the drag penalty at higher Reynolds numbers.For the purposes of this study, we assume the flow to be "fully rough" if the skin-friction coefficient C f is invariant with Reynolds number (i.e.data (7) can fit equation 7).It is also worth recalling that boundary layer flows are always developing, and their C f and (per- haps even Π ) remain a function of the streamwise distance x.Therefore, the assessment of k s in this study is associated with a fixed streamwise location.

Contribution and outline
It is evident that in the absence of direct drag measurements, determination of the skin friction hence ΔU + (k s ) can only be achieved using velocity profile data at the price of either assuming the validity of outer-layer similarity or other alternative ways such as the modified Clauser-chart (which requires fittings U , d and ΔU + while fixing and B in equation 1) or total shear stress methods.The choice of either will most certainly affect the drag prediction at higher Reynolds numbers.However, the implications of these steps and assumptions on the final prediction remain unclear.
In this work, direct drag measurements are carried out using an in-house floating element drag balance previously validated by Ferreira et al. (2018), while velocity field measurements are performed using particle image velocimetry.The experiments are performed over a realistic rough surface with the objective of examining the limitations of indirect methods in determining the equivalent sand-grain roughness height.Once the value of k s is determined, it is possible to make predictions of the total frictional drag ( C F ) at higher Reynolds numbers using the integral boundary layer evolution equations as outlined by Monty et al. (2016), which in turn will allow the estimation of the total ship resistance.We also compare the results from indirect methods in determining the aerodynamic roughness parameters and examine the implications of using them for drag penalty predictions at high Reynolds numbers over moderately fouled ship hull surfaces.The different methods used in this study are summarised in Table 1.
The remainder of the manuscript is presented in three sections.The experimental methodology is described in Sect. 2 depicting the roughness topography, drag, and flow field measurement techniques.The results and discussion

Facility and biofouling roughness
Measurements are conducted in an open-circuit suctiontype wind tunnel at the University of Southampton.The test section is 0.6 m × 0.9 m in the (y, z) cross-plane and spans 4.5 m in the x-wise direction, following a 7:1 contraction.The test section was designed with a weak diverging cross section to allow a constant free stream along the streamwise direction and the growth of a turbulent boundary layer with a nominally zero pressure gradient.The dx was shown to range between 1 and 5 ×10 −8 from various recent studies per- formed in this facility (Ferreira and Ganapathisubramani 2021;Medjnoun et al. 2021).The turbulent boundary layer develops over a flat surface composed of five equally-sized wooden boards onto which the biofouled surface roughness is mounted.In order to assess the skin friction, a 0.2 m × 0.2 m slot was made in one of the boundary layer plates at the measurement location to accommodate the floating element drag balance.The boundary layer plates are preceded by a ramp of 0.2 m long inclined by four degrees to the horizontal, ensuring a smooth transition of the flow from the bottom floor of the test section.The free stream speed can reach up to 30 ms −1 , with a turbulence level less than 0.5% , and is monitored and acquired using a micromanometer FC510.To account for air density variations, the temperature is also acquired, and its standard deviation for an average run was less than ±0.5 o Celsius.
The roughness was made of a tiled surface manufactured from a moderately fouled coupon.The coupon is taken from a dry-dock silicone mould obtained directly from a ship hull.The coupon consists of some tubeworm fouling as well as one barnacle.The size of the coupon was 25 mm × 50 mm and this was repeated in two directions to generate a repeating-unit tile of 100 mm × 100 mm.The surface properties of the coupon are also listed in the left panel of Table 2 while the surface topography is shown on the right panel.To manufacture these surfaces, four master tiles are initially 3D printed using a Formlab Resin printer with a height resolution of 10 m equivalent to one viscous scale at the highest Reynolds number of the facility.These master tiles were used to create several silicone moulds, which were subsequently used to cast over 400 replicas of polyester resin needed to cover the entire floor of the wind tunnel test section.The boards are mounted in two configurations in the wind tunnel generating two different roughness patterns from the same tiles, which 169 Page 6 of 20 can be depicted by the directionality of effective slopes ES (Napoli et al. 2008) as shown in Table 2. Therefore, the spanwise characteristic wavelength of Config 1 is equal to the streamwise wavelength of Config 2, and vice-versa ( S z,1 = S x,2 = 50 mm and S z,2 = S x,1 = 25 mm).

Floating-element drag balance
Wall shear stress is directly measured using the floating-element drag balance developed and validated by Ferreira et al. (2018).Its design is based on the parallel-shift linkage, featuring dedicated pairs of bending-beam transducers to monitor not only the streamwise load but also the induced pitching moment, as a means to decouple extraneous loads.The balance itself is a 200-mm-side square that is small enough in relative terms to ensure local measurements of wall shear stress.The balance is flush mounted with the wind tunnel floor, roughly 3.2 m downstream of the contraction.Special care was taken to mounting the surrounding tiles of the balance, given the tight tolerance of the fix-to-floating surface joint that is only 0.5 mm wide.It is emphasised that the dimensions of the balance are a multiple of a tile size, which is 100-mm-side square.Therefore, on the floating element, four complete tiles were placed.
In order to determine C R f as a function of the Reynolds number and determine the aerodynamic roughness lengthscale of the rough surface, the floating element is subjected to a series of nine free stream velocities ranging from 9 up to 27 ms −1 .Each acquisition lasted 30 s sampled at 150 Hz, equivalent to 2500 boundary-layer eddy turnover times ( eddy = ∕U ∞ ) at the lowest operating speed, with a total of five repetitions per velocity.Pre-and post-calibrations were conducted for each configuration with a change in the calibration coefficient less than 0.5% .The acquisition procedure and uncertainty analysis are detailed in Ferreira et al. (2018).

Particle image velocimetry
In order to examine the flow, velocity fields are obtained from both planar and stereoscopic particle image velocimetry (PIV) measurements at a similar location to that of the drag measurements.All planar and stereoscopic PIV measurements are performed in the (x, y)-and (y, z)-planes, respectively.For the first surface configuration, seven PIV measurements are taken at various free stream speeds ranging from 7 up to 23 ms −1 .This included six planar-PIV cases with five above the valley (between two successive barnacle centres) and one additional plane at the centre of barnacles, while one stereo-PIV plane was also taken above the barnacles.For the second surface configuration, another seven PIV measurements are also acquired at similar free stream speeds to the first configuration.This consisted of five stereo-PIV planes at the top of the barnacles with two additional planar-PIV cases; one at the top of the barnacles and another between them.A summary of the different cases is provided in Table 3.
To trace the flow, vaporised glycerol-water particles produced by a Magnum 1200 fog machine are used, then illuminated with a laser-light sheet sourced by a two-pulse Litron Nd:YAG laser operating at 200 mJ.A LaVision optical system for the beam focus/expansion of the light sheet is employed, which is comprised of convex and concave lenses in order to focus the beam and a cylindrical lens in order to expand the sheet with relatively constant thickness in the measurement plane (roughly 1 mm thickness).The particle images are recorded by high-resolution LaVision Imager LM 16 MP CCD cameras fitted with 200 and 300 mm AF Micro Nikkor lenses for the planar-and stereo-PIV planes, respectively.One camera was used for the planar-PIV setup and is positioned at roughly 1 m away from the object plane, whereas two other cameras needed for stereo-PIV planes are mounted on Scheimpflug adapters to account for the oblique view angle ( ±42 o ), and placed at approximately 1.3 m from the object plane at either side of the test section.A single and a double-sided dual-plane calibration target aligned with the laser light sheet were used to determine the mapping function for each setup, using a third-order polynomial fit.This resulted in a field of view of nearly 0.9 × 1.5 in the (x, y)-plane and 1.3 × in the (y, z)-plane for the planar and stereoscopic PIV, respectively.
For both configurations, 500 to 3000 independent realisations of image pairs are acquired for the Reynolds number sweep at 0.6 Hz, with a time delay ranging between 25 and 70 s resulting in an average particle displacement of 20 pixels/s.The velocity vector fields are then obtained by interrogating particle images using a decreasing multipass scheme starting from 48 pixels × 48 pixels down to a final pass of 16 pixels × 16 pixels for the planar PIV and 24 pix- els × 24 pixels for the stereo PIV.Using a 50% interroga- tion window overlap, the resulting effective vector spacing ranged between 0.2−0.4mm, equivalent to an inner-normalised volume resolution of Δx + × Δy + × Δz + ≈ 7.4 × 7.4 × 32 for the lowest Reynolds number up to 35.7 × 35.7 × 87.4 for the highest Reynolds number (also depending on the different planar and stereoscopic PIV measurements).

Results
This section focuses on results reported from the analysis of the aerodynamic roughness parameters evaluated using different methods.Section 3.1 discusses the results obtained from the floating element and predicts the roughness lengthscale obtained from direct drag measurements together with equation 6.The second Sect.3.2 discusses the aerodynamic properties obtained using both direct drag and velocity measurements.Section 3.3 examines the effect of using indirect methods in the assessment of the aerodynamic properties of the roughness.Finally, Sect.3.4 attempts to reconcile these results by showing how discrepancies from different methods affect skin friction predictions at higher Reynolds numbers.

Method 1: DB + Eq. 6
Figure 1 shows the variation in skin-friction coefficient for the two different configurations, alongside smooth wall data obtained using oil-film interferometry (OFI) in the same facility.These measurements were obtained at various free stream velocities and therefore the abscissa in Fig. 1a represents the Reynolds number Re x at a fixed streamwise fetch (3.2 m downstream the leading edge of the surface roughness), while Fig. 1b shows the same data where PIV measurements are performed and plotted against the friction Reynolds number Re .It is clear that both configurations of rough walls yield a much higher skin-friction coefficient than a smooth wall with an increase of 320% and 280% at Re x ≈ 5 × 10 6 for Config 1 and Config 2, respectively.Both figures show that C f has reached an asymptotic limit for each case as it appears to not change as the Reynolds number increases.This suggests that both surfaces have reached a fully-rough regime whereby the skin friction is invariant with the Reynolds number (for a given streamwise fetch).
Furthermore, the skin friction produced by surface Config 1 is 9 .5% higher than that generated by Config 2. The difference can be attributed to the directional alignment present in the biofouled rough surface, as indicated by the effective slope (ES) in both streamwise and spanwise directions, despite both surfaces sharing similar statistical characteristics.Following the surveys of Napoli et al. (2008) and Schultz and Flack (2009) on turbulent boundary layers over irregular and regular rough-walls, respectively, both current configurations can be considered to be within the roughness regime (roughness regime if ES > 0.35 and waviness regime otherwise).It is also seen that Config 1 has an effective slope 15% larger than that of Config 2. This result suggests that this drag increase could be due to contributions from an increase in the flow separation around steeper roughness features.It is also intuitive to assume that the drag can easily be influenced by the change in the spanwise/streamwise characteristic wavelengths S x and S z .However, the latter parameters have both changed from Config 1 to Config 2, and it would be impossible to ascertain which parameter had more impact in increasing the drag.
The skin-friction coefficient can be used to estimate the roughness function and hence the equivalent sand-grain roughness height.Following equation 6 which assumes outer-layer similarity, the roughness function ΔU + DB can be determined from solely the friction velocity and boundary layer thickness information.The latter can easily be obtained with a quick survey of the mean flow using Pitot, hot-wire anemometry, or as in this study, PIV measurements.Since outer-layer similarity implies the wake strength parameters between smooth and rough walls to be matched ( Π R = Π S ), the roughness function ΔU + DB can be obtained by either matching + of both the smooth and rough-wall measurements or by simply adding the 1 ln + R + S term to equation 6 if the + 's are unmatched.In fact, our results showed both methods yield identical ΔU + DB estimates.
Results of the roughness function ΔU + DB under the assumptions of equation 6 (matched + and Π between smooth and rough-wall flows) are depicted in Fig. 2(a) showing its variation with respect to k + rms .The outcome of choosing k rms or any other representative roughness height remains the same in the determination of the equivalent sand grain roughness height k s,DB .The results of the rough- ness function for both cases are shown to vary logarithmically with this length scale as a consequence of the fullyrough regime observed in Fig. 1.
Using the appropriate value of k s,DB , the same data can be collapsed to the fully-rough line represented by equation 7 as shown in Fig. 2b.The results show that even at the lowest Reynolds number, all measurement points exhibit a k + s,DB > 300, well within the fully-rough regime.The large values in the reported k + s,DB are a consequence of the corresponding large values in ΔU + DB (over 11 at the lowest Reynolds number) which can be considered as a strong roughness in both configurations.These surfaces are shown to lead to considerable values of k s,DB when compared with their geo- metrical attributes with k s,DB ≈ 2.6k p and k s,DB ≈ 1.8k p or approximately 15% and 10% when expressed as a fraction of the boundary layer thickness for Config 1 and Config 2, respectively.The subsequent analysis will investigate how the roughness function and other aerodynamic parameters obtained from velocity measurements compare with the values presented above.Although PIV fields are not shown here for brevity, results for both configurations indicate that the mean and turbulent flow heterogeneity remain confined within the roughness sublayer ( y < 0.1 ) in both the (x, y)-and (y, z)-planes, demonstrating that these surfaces (flows) behave as homogeneous rough surfaces (flows).Hence, the subsequent results remain applicable within the homogeneous roughness framework.The following analysis presents the aerodynamic roughness parameters assessed from direct drag together with mean velocity profile measurements which are spatially averaged over from 2 to 4 repeating units (depending on the case), to obtain statistically converged data.For the sake of clarity, only one profile per configuration, from the sPIV data at matched Re x ( Re x ≈ 4 × 10 6 ) are presented for comparison.Two inner-normalised velocity profiles scaled with the friction velocity from the drag balance at the same Re x for both surface configurations are presented in Fig. 3 and compared with a smooth-wall turbulent boundary layer profile from the DNS of Sillero et al. (2013).The velocity profiles are shown to exhibit a substantial shift in the vertical axis indicative of the momentum loss caused by the higher wall shear stress experienced by the flow in comparison with the smooth wall.Config 1 is observed to have a relatively larger downward shift in comparison with Config 2 despite having a similar free stream velocity.However, this is expected since the skin-friction coefficient of Config 1 is larger than that of Config 2.
The zero-plane displacement (d) and the roughness function ( ΔU + ) which characterise the shift in the velocity pro- files observed in Fig. 3 are obtained by assuming the loglaw.The former is determined by taking the derivative of the logarithmic part of equation 1 with respect to y, whereas the latter is obtained by simply taking the difference between the inner-normalised measured velocity profile and the log-law.These are expressed as and Equation 9 is traditionally employed as a means to infer the logarithmic slope (Von-Kármán constant ) and to provide a conservative indication of the log region; thus Ξ is gener- ally referred to as the indicator function (Österlund et al. 2000).Consequently, by assuming to be universal between smooth-and rough walls, d can be determined by measuring the velocity gradient and the friction velocity (as done in this study).The zero-plane displacement is found as the value that minimises the difference between the left-hand-side ( Ξ ) and the right-hand-side ( 1∕ ) of equation 9.The results shown in Fig. 4a indicate that the best solution that minimises this difference (the value of d + that yields Ξ ≈ 1∕ ) occurs around 0.15 , leading to a good collapse of both profiles.
The results have shown that the values of d have a degree of scatter, however, this is believed to be caused by the uncertainty in locating the wall-normal origin due to the nature of these rough surfaces.The average zero-plane displacement as a fraction of the maximum roughness height d∕k p = 0.28 and 0.57 for Config 1 and Config 2, respectively.Nonetheless, the variability in d was shown to have only a marginal influence on the roughness function.For instance, there is only a 1 % change in ΔU + if the value of d is to be fixed at d = k p ∕2 , compared to the values returned from the fit described in equation 9.The roughness function ΔU + is subsequently determined by examining the Ψ profiles.Fig- ure 4b presents the profiles of both surfaces and the results show the optimum fit solution that returns the best plateau occurs within the overlap region (highlighted grey region) 2k p < Ψ fit < 0.2 .The profiles are also compared with the smooth-wall DNS data of Sillero et al. (2013) and the figure clearly indicates that the rough profiles have a weaker wake region.In fact, the wake strength parameter can be deduced through the modified-log function Ψ , which quanti- fies the maximum departure from the logarithmic behaviour; Π = ∕2 × max(Ψ) .This method is sometimes favoured as it circumvents the use of a particular function.This also means that it can easily be inferred that the weaker the maximum value of Ψ gets, the weaker the wake strength parameter becomes.
With the roughness function values as a function of Reynolds number known, the equivalent sand-grain roughness height can be deduced using the fully-rough relation (equation 7).The variation of ΔU + with k + s obtained from direct drag and velocity measurements is shown in Fig. 5a for both configurations.As expected, the data are shown to line up very well with the fully-rough asymptote already hinted at from the C f (Re x ) plot.The determined values of k s are reported to be 11.6 mm and 8.05 mm, corresponding to 1.85k p and 1.28k p for Config 1 and Config 2. As fractions of the boundary layer thickness, they represent 10.3% and 7.4% for Config 1 and Config 2, respectively.These values are seen to be lower than those reported by direct drag measurements and using equation 6, indicating that k s,DB overesti- mates the roughness lengthscale by as much as 38% .Using these values, the inner-normalised mean velocity profiles can be re-plotted against the wall-normal distance normalised by k s as depicted in Fig. 5b.
The results exhibit a good collapse of the two configurations within a range of 1 < (y − d)∕k s < 2 corresponding to the overlap region where ΔU + was estimated.They are observed to follow a logarithmic behaviour represented by the black-dashed line, with a fully rough intercept B r ≈ 8 , as opposed to the classically reported 8.5.It should be recalled that the latter value is derived from pipe flow data.In fact, the logarithmic slope for boundary layers (currently taken as 0.39) is smaller than that observed in pipe flows, generally reported as being 0.41 (Marusic et al. 2013), hence the associated intercept B r being lower to the classically used value in pipes and channel flows.It is worth noting that the selected value for C remains applicable to boundary layers, as evidenced from the current intercepts for both smoothand rough-walls, calculated as B r − B = 8 − 4.5 = 3.5 = C .The profiles are also shown to have a weaker departure from the logarithmic distribution in the outer region, but this is perhaps clearer when comparing the Ψ profiles with respect to the smooth wall as presented in Fig. 4b.
Using the values of k s obtained from the profiles and drag-balance, we can also predict the change in Π R .This prediction can subsequently be checked against the value of Π R determined directly from the velocity profile data.Using equation 7 together with 5 and 6 at matched values of + , it can be shown that: With the values of k s and k sDB , the average values of Π R are predicted to be 0.41 for both Config 1 and Config 2, respectively, as shown in Fig. 6a as blue-dashed lines.This (11) prediction confirms the profiles are expected to have a weaker wake parameter in comparison with the smooth wall.
To examine Π R , it is possible to use either formulation of the velocity profiles expressed by equations 1 or 2.Although different analytical expressions for W yield different values of Π (as stated in the introduction), the conclusion drawn in the current study remains unchanged when employing alternative expressions.Therefore, instead of using an analytical expression, equation 3 is currently preferred as it directly quantifies Π R without the need for a fit, thus eliminating additional fitting uncertainties.Rearranging the terms in equation 3 yields the expression, Π R is subsequently determined directly given all quanti- ties ( U + ∞ + ΔU + − B ) have been measured.The results of the wake strength parameter with respect to Reynolds number are depicted in 6a together with the prediction (blue line) showing the values to be substantially lower than that of a smooth-wall of Π ≈ 0.57 (determined using Sillero et al. (2013) DNS smooth-wall data).Despite the scatter in the results, both surfaces are shown to yield relatively similar mean values ranging between 0.38 and −0.44 (across the different speeds and cases), with an average value of 0.41 nearly a 30% reduction in wake strength.In hindsight, these results provide an explanation for the higher values observed in ΔU + DB (hence k s,DB ) obtained using drag information with the assumption of outer-layer similarity without prior knowledge of the flow behaviour.Examining equations 5 and 6 shows this difference can only be attributed to differences in the value of Π between smooth and rough walls.This would suggest that velocity profiles would not follow universal outer-layer similarity.
Additionally, the values of Π R determined from equa- tion 11 are shown to be in agreement with the ones measured from the profiles.The latter is shown to have a degree of variability ( ±5% ) but this is expected given the change in Reynolds number.The values of Π R are also reported in Table 4 alongside other relevant aerodynamic parameters for both configurations.The similarity in Π R for Config 1 and Config 2 despite their difference in ∕k s reveals that perhaps outer-layer similarity between Config 1 and Config 2 can still be valid.
In a similar fashion to Figure 3 presented by Squire et al. (2016), Fig. 6b presents the current data in the ( ∕k s , k + s ) -plane along with data from other rough-wall turbulent boundary layer experiments; the sand-grain roughness and woven mesh of Flack et al. (2007), the sand-grain roughness of Squire et al. (2016) and the urban roughness of Ferreira and Ganapathisubramani (2021), with the latter being performed in the same facility as in the current study.The horizontal axis reflects the ratio of the aerodynamic roughness to the viscous length scales and is generally indicative of the flow regime (i.e.fully-rough regime for k + s > 100 , aerody- namically-smooth for k + s < 6 and transitionally-rough regime in between).On the other hand, the vertical axis indicates a scale separation between the boundary layer thickness and the aerodynamic roughness height equivalent to the friction to roughness Reynolds numbers ratio ( ∕k s ≡ + ∕k + s ).In this example, the plane is divided into four quadrants delineating regions of; mild roughness in transitionallyrough regime ( Q 1 with ∕k s > 10 and k + s < 100 ), mild roughness in fully-rough regime ( Q 2 with ∕k s > 10 and k + s > 100 ), strong roughness in fully-rough regime ( Q 3 with ∕k s < 10 and k + s > 100 ) and finally strong rough- ness in transitionally-rough regime ( Q 4 with ∕k s < 10 and k + s < 100 ).According to this chart, the current roughness seems to generate a relatively strong aerodynamic roughness lengthscale compared with the P36-sand-grain roughness tested by Squire et al. (2016).It instead compares well with the medium-mesh and P24-sand-grain cases examined by Flack et al. (2007) but remains weaker than the cube roughness generally used to model urban roughness as demonstrated by Ferreira and Ganapathisubramani (2021).
As stated in the introduction, the outer-layer similarity hypothesis implies that the roughness sublayer confines itself to the inner region while a large portion of the flow remains unaffected by the surface condition.To assess the validity of this hypothesis, the mean velocity deficit along with the streamwise turbulence intensity profiles are examined in the wake region (which nearly represents the upper 80% of the boundary layer).The results from both configurations are illustrated in Fig. 7a and b, and compared with their equivalent smooth-wall profiles from Sillero et al. (2013).It should be noted that despite the unmatched + 's of the different profiles (the DNS date has a friction Reynolds number of + ≈ 2000 ), the wake region is only weakly dependent on Reynolds number, hence the comparison remains valid.The results presented in Fig. 7a clearly show the absence of similarity between the smooth-and rough-wall profiles, with current profiles presenting a smaller deficit compared with the smooth wall (solid line).The lack of similarity extends almost up to 0.7 -0.8 , which represents nearly 7 k s and 10k s for Config 1 and 2, respectively.
The figure also compares the current results with roughwall profiles from surfaces made of regularly distributed cubes with random (C10R) and uniform (C10U) heights, examined by Ferreira et al. (2018) within the same facility.These types of surfaces are commonly utilised to model urban flows and have been documented in the literature (Castro 2007;Amir and Castro 2011;Castro et al. 2013).Both cube rough surfaces yield different behaviours in the outer flow in comparison with the biofouled surfaces.
In fact, these surfaces maintain a good degree of similarity (in the mean flow) with the smooth-wall flow for y∕ > 0.2 , with a wake strength parameter of Π ≈ 0.42 and 0.51 for C10R and C10U, respectively.These are associated with a ∕k s ≈ 5 and 7 for C10R and C10U, which can be considered a relatively stronger roughness than the current ones, also seen in the ( ∕k s , k + s )-plane of Fig. 6b.Moreover, the blockage due to roughness in both experiments is also similar meaning their corresponding acceleration parameters are very similar.This suggests that the absence of outer-layer similarity cannot be attributed to a small value of ∕k s , nor to the weaker value in the wake strength parameter alone (if caused by a favourable pressure gradient).However, their influences cannot be fully ignored.
Figure 7b shows the inner-normalised streamwise turbulence intensity profiles also compared with the smooth-wall profile.These profiles have been corrected to account for the energy attenuation caused by the filtering of the subresolution scales in the PIV measurements, using the method outlined by Lee et al. (2016).It is seen that the turbulence intensities have only been very marginally affected in the outer region (less than 1 % attenuation in the outer region.Overall, the results clearly indicate that the turbulence intensities have a smaller magnitude in comparison with their equivalent smooth-wall profile.This observation reveals that the increase in frictional drag caused by the roughness is not systematically accompanied by a proportional increase in the turbulence intensity, as evidenced by the absence of outer-layer similarity.

Method 3 and 4: OLS and CSS
In the absence of drag balance measurements, U has to be indirectly determined and then used to obtain ΔU + , k s , and Π .As opposed to direct methods, indirect methods  (2013).The profiles are also compared with the fully-rough profiles (C10R and C10U) from Ferreira and Ganapathisubramani (2021) measured in the same facility.The turbulence intensity profiles have been corrected using the method proposed by Lee et al. (2016) to account for the energy attenuation due to the filtering that resulted from the PIV spatial resolution Page 13 of 20 169 essentially rely on mean and turbulence flow assumptions to infer the wall shear stress.For example, Clauser chart method (Clauser 1956) is widely used to infer U in smooth- walls, while its modified version applied to rough-walls (Perry and Li 1990), determines U , ΔU + and d altogether from the log region for a given and B. Alternatively, U can be determined from the total shear stress method, which assumes a constant shear stress layer in the overlap region to exist, assumed to be equal to the wall shear stress (Schultz and Flack 2005).Other methods rely on the flow's streamwise growth, such as the Von-Kármán momentum integral method.In this work, two methods at two ends of complexity are evaluated: 1) assumption of outer-layer similarity for streamwise mean and turbulence intensity (Townsend 1976) and 2) comprehensive shear-stress method recently proposed by Womack et al. (2019).Results from both methods are presented in the following and compared to the direct measurement results.
Assuming the validity of the outer-layer similarity hypothesis implies that both f and g ij in equa- tion 4a are universal for the smooth-and rough-wall flows.This is clearly not the case for both currently investigated surfaces, however, it is important to know the effect of this assumption as this method has been previously employed to determine full-scale frictional drag and predict the drag penalty of ships at high Reynolds numbers (Monty et al. 2016).Forcing collapse of profiles is achieved by simultaneously minimising the difference in both velocity deficit and turbulence intensity between the smooth-and rough-wall data in the outer region ( y > 0.2 ).This results in a new U value which is subsequently used to evaluate ΔU + from the new inner-scaled profiles, for which an appropriate k s is determined using equation 7.
In the second method, the aerodynamic parameters are estimated using the comprehensive shear-stress method proposed by Womack et al. (2019).This method builds on the procedure initially outlined by Volino and Schultz (2018) and follows an iterative error minimisation procedure to determine U , d, y 0 (equivalent to ΔU + or k + s ) and Π using integral momentum equations -total shear stress balance together with the law-of-the-wall.Womack et al. (2019) demonstrated by using a single streamwise velocity and a Reynolds shear stress profile, U can be retrieved to within ±1 − 4% when compared with direct drag meas- urements.To allow a consistent comparison between the direct method and this indirect method, the fitting region adopted for the log-law is similar to the previously used, 2k p ≈ 0.1 < y∕ < 0.2 .For the Reynolds shear stress pro- files, the optimisation was carried out within a range of 0.2 < y∕ < 0.3 .The interactive graphical user interface software made available by Womack et al. (2019) was used to carry out this analysis.
Results of the aerodynamic parameters from the different methods are presented in Table 4, and their differences are illustrated in Fig. 8 for both configurations.Figure 8a compares the percentage in the relative difference of the friction velocity ( E U (%) = (U − U ,indirect )∕U ) obtained from the two indirect methods (OLS and CSS) with the measures from the drag balance.Despite the noticeable scatter, the results indicate that in both methods, U is underestimated for both configurations.The CSS method is shown to give consistently relatively better results as the Reynolds number Fig. 8 Percentage in the relative difference of the friction velocity E U (%) = (U − U ,indirect )∕U ) obtained from the two indirect methods (OLS and CSS) with respect to the direct measures from the drag balance.b Percentage in the relative difference of the roughness function ( E ΔU + (%) = (ΔU + − ΔU + indirect )∕ΔU + ) obtained from the three indirect methods (OLS and CSS in addition with using equation 6) compared with the measures using U ,DB and the log-law profile 169 Page 14 of 20 increases, with the difference reaching less than 3 % , which falls within the accuracy range results provided by Womack et al. (2019).On the other hand, the OLS method is shown to underpredict the friction velocity by over 10% at low speeds then the difference reduces nearly to 5 % as the speed increases.
The relative difference in estimating the roughness function ( E ΔU + (%) = (ΔU + − ΔU + indirect )∕ΔU + ) is shown in Fig. 8(b) and compares the three different methods; OLS, CSS and equation 6 (which uses direct drag measurements but assumes Π S = Π R ), against the log-law method which is expected to be the most accurate method (since it uses U measured with the drag balance without any mean and turbulence flow assumptions).It is clear from this figure that equation 6 overestimates the roughness function when compared with the correct ΔU + , while the other indirect methods tend to underestimate the value of ΔU + .
Similarly to the degree of scatter observed in Fig. 8a, the relative difference in the roughness function shown in Fig. 8b seems to reduce as the Reynolds number increases.Nonetheless, these differences are enough to cause a considerable underestimation in k s as shown in Table 4, with k s,OLS and k s,CSS underpredicting k s by 32% and 15% for Config 1 and 31% and 23% for Config 2, respectively.Hence, this could have an important effect when predicting the drag penalty at higher Reynolds numbers.
It is shown in Table 4 that the OLS method results in a wake strength parameter value of roughly 0.48 which is relatively larger than the previous estimate, however, remains about 15% smaller than the smooth-wall Π S of 0.57.This is believed to result from the absence of a complete simultaneous collapse when imposing similarity in the mean and turbulence statistics, but also due to subtle differences in the shape of the wake region itself.This behaviour is also observed in the C10R and C10U profiles, in which, despite the relatively good collapse of the mean profiles with the smooth-wall data, the overlap region remains noticeably deviated from that of the smooth wall.For the CSS methods, the wake strength parameter showed a degree of scatter  4 Aerodynamic parameters determined from four different methods namely; method 1 using U ,DB and equation 6, method 2 using U ,DB and log-law profile, method 3 using outer-layer similarity hypothesis, and method 4 using the comprehensive shear stress method of Womack et al. (2019), respectively.The estimates of the equivalent sand-grain roughness height from the different methods for Config 1 are: k s,DB = 16.11mm, k s = 11.59 mm, k s,OLS = 7.88 mm and k s,CSS = 9.76 mm associated with the methods 1, 2, 3 and 4, respectively.The results for Config 2 are: k s,DB = 11.13 mm, k s = 8.05 mm, k s,OLS = 5.56 mm and k s,CSS = 6.13 mm associated with the methods 1, 2, 3 and 4, respectively across different speeds.However, this is inherently related to the way Π is determined, since Womack et al. (2019) defined it as the wake's maximum departure from the extrapolated log-law, to avoid assuming a universal wake function.
Figure 9 illustrates the different results in using indirect methods in estimating C f in these laboratory-scale surfaces.Estimates of the skin-friction coefficient from both OLS and CSS methods, compared with the direct drag measurements are exemplified for Config 2. The results clearly show that both methods underpredict C f , with the CSS method under- estimating C f by close to 6 % at the highest Reynolds number, as opposed to 10% when assuming the outer-layer similarity hypothesis.

Drag penalty predictions
Now that direct and indirect methods have been assessed, a prediction of the spatially-averaged frictional drag C F at high Reynolds numbers Re x (with fixed U ∞ and but increasing x) can be made through each method.To this end, the procedure proposed by Monty et al. (2016) is employed.It is an integral boundary layer evolution method based on numerical integration of the skin friction over a given length of the developing flow, assuming a velocity profile composed of a logarithmic and wake region (see Monty et al. (2016) for more details).The profile can be fully described by its freestream velocity U ∞ , kinematic viscosity , and wake strength parameter Π along the equivalent sand-grain rough- ness height k s .For the sake of consistency with the above analysis, the logarithmic region uses a = 0.39 and B = 4.5.
Figure 10 shows the variation of C F as a function Re L , highlighting full-scale predictions for each method (colour-coded curves) and depicting the influence of the wake strength parameter (dashed curves).The graph is shown for Config 2, but the same behaviour is also observed for Config 1.This figure is obtained for a unit Reynolds number U ∞ ∕ ≈ 2 × 10 7 .The results show that the relative change in frictional drag with respect to the smooth-wall ( ΔC F,s = (C F − C F,s )∕C F,s , with C F,s being the smooth- wall frictional drag) is over 230% at a Reynolds number of Re x ≈ 10 9 (typical of a Reynolds number of a 50 m long vessel) when using the direct method ( k s , Π R ).This rela- tively large value in ΔC F,s is not unexpected as according to Schultz (2007), the equivalent sand-grain roughness of the current surface should be regarded as heavy calcareous fouling ( k s ∼ O(10) mm).Similarly, the same result can be achieved when using k s,DB and Π S as shown by the relatively good collapse of both predictions (solid black and blue curves).The reason for this is trivial and is explained by equation 11 which shows that the overestimated value in k s,DB is compensated for by the assumption of Π S .Hence, the skin friction prediction using ( k s,DB , Π S ) is similar to using ( k s , Π R ).
Furthermore, it's worth noting that both indirect methods, OLS and CSS, tend to provide underestimated values for C F ; however, the CSS method shows relatively better accuracy.On the other hand, the assumption of Π = Π S clearly shows that the drag is systematically underpredicted.To highlight these differences, three different fetch Reynolds numbers have been selected for comparison as seen in Fig. 10 (vertical dashed lines).At Re L = 10 8 , the prediction using ( k s,DB , Π S ) overestimated the predicted value by less than 1.5% , whereas using ( k s , Π S ) the prediction worsens to 6-7% below the direct prediction.As observed from the graph, the substantial undervaluation scenario occurs if the outer-layer similarity hypothesis is assumed even at these high Reynolds numbers, as this results in a drag underprediction of around 14-15% .On the other hand, the CSS is seen to give relatively better predictions with underestimates around 6-7% when using Π R but worsen to over 11-12% when using a Π S .
These different estimates remain consistent with increasing Reynolds numbers, albeit with a clear tendency of reduction in the relative difference due to the asymptotic nature of the frictional drag at very high Reynolds numbers.This can be explained by the fact that at very high Reynolds numbers, most of the drag is expected to be captured by the logarithmic region which increases proportionally with Re L , while the wake region (hence changes in Π ) will have a propor- tionally lesser impact as observed.Quantitative figures are provided in Table 5 for both configurations, summarising the influence of using the different methods on drag prediction.Nevertheless, it is important for the reader to be cautious when interpreting these relative differences, as they are specific to the particular roughness under investigation.In fact, other rough surfaces will likely produce different values.This was observed after carrying out the same analysis described in Sect.3.3 on the surfaces examined by Ferreira et al. (2018).The results showed that the OLS and CSS method only underpredicted the value of C F by roughly 3 % at the highest Reynolds number.Therefore, the OLS and CSS methods remain potentially excellent alternatives to assess the aerodynamic roughness parameters and predict the drag penalty at high Reynolds numbers over a wide range of rough surfaces, when direct drag measurements are not available (Monty et al. 2016;Womack et al. 2019).It is also recalled that these predictions are only valid for flat surfaces, without pressure gradients, and care should be taken when interpreting these results for practical engineering applications (e.g.flows under favourable and/or adverse pressure gradients due to surface curvatures).
Finally, it is essential to keep in mind that the above analysis predicts drag within the classical framework of homogeneous roughness.However, in practical applications, this tool does not account for the inherent patchiness or heterogeneity found in several real-world applications.The lack of comprehensive data on this topic poses significant challenges for our current predictive framework, as noted by Chung et al. (2021).This limitation primarily stems from our limited understanding of the drag behaviour over heterogeneous rough surfaces (Medjnoun et al. 2018(Medjnoun et al. , 2020)).In their recent work, Hutchins et al. ( 2023) have proposed a power-mean approach for defining an equivalent homogeneous roughness length that takes into account patchiness or the heterogeneous distribution of roughness on ship hulls.This approach seems promising, especially when dealing with patch sizes significantly larger than the boundary layer thickness.It has the potential to address some challenges posed by heterogeneous roughness, however, its application in existing full-scale drag prediction methods highlights the need for further investigations in this field.Obtaining highfidelity data across a broad range of Reynolds numbers is particularly crucial in advancing our understanding of this complex problem.

Conclusions
A turbulent boundary-layer flow over two configurations of a realistic biofouled rough surface scanned, scaled and replicated from a ship hull is experimentally investigated.Detailed laboratory wind-tunnel experiments have been carried out, which consisted of direct drag measurements using an in-house floating element drag balance together with PIV to assess the flow field.
Data from direct drag measurements revealed that the skin-friction coefficient reached the fully-rough regime under both surface configurations, shown by the invariance of C f with respect to the Reynolds number.In these labora- tory-scale surfaces, the frictional drag is shown to exceed 300% that of a smooth wall at matched Re x .Using this information along with velocity measurements, limitations of similarity laws have been discussed.When considering only direct drag and boundary layer thickness data, it becomes evident that the aerodynamic roughness lengthscale tends to overestimate the actual value, sometimes by as much as 40% .This overestimation is shown to be caused by the reduction in the wake strength parameter, which itself is shown to be 30% weaker than that of a smooth wall.
Lack of outer-layer similarity in both mean flow and turbulence intensities is reported, despite following the classically admitted conditions (  + > 3000 , k + s > 200 and k s ∕ ≈ 0.1 ).Although flow fields were not shown for brev- ity, the mean and turbulent structures remain generally unperturbed in the outer region, with an excellent degree of spanwise homogeneity.The canopy layer harboured intense shear layer regions generated by the barnacles which seemed to propagate towards the outer layer, albeit confined below y∕ < 0.15 .This means that the lack of universal outer- layer similarity can potentially be caused by the generation of roughness-scale structures that emanate from the barnacles, leading to changes in the uniform momentum zones and the wake intermittent activity resulting in a weaker wake parameter.
By assuming the validity of the outer-layer similarity hypothesis and using the comprehensive shear stress (Womack et al. 2019) methods, indirect estimates of the aerodynamic parameters have been obtained and compared with results from the direct measurements.The comparison highlighted the limitation of using these flow assumptions, which resulted in underrating the equivalent sand-grain roughness on average for both configurations by over 30% and 20% for the OLS and CSS methods, respectively.Using the integral boundary layer evolution method outlined by Monty et al. (2016) results from the different indirect methods have been used to assess the frictional drag penalty at high Reynolds and contrasted with the predictions from direct measurements.
The current results also showed that accurate prediction of skin friction remains possible using an indirect estimate of the roughness function ΔU + DB , provided direct drag measurements at laboratory-scale data are available.In fact, for a nominally zero-pressure gradient turbulent boundary layer developing over a flat surface, k s,DB (obtained from direct drag measurements) together with a wake strength parameter Π = Π S result in a similar frictional drag prediction as when using directly measured estimates of k s and Π R .How- ever, when indirect drag estimates are used to characterise the roughness, frictional drag predictions at high Reynolds numbers can be underestimated.At a Reynolds number of Re x ≈ 10 9 (equivalent to a Reynolds number of a 50 m long vessel at full-speed), assuming outer-layer similarity underestimated the frictional drag by over 12% while the comprehensive shear stress method resulted in approximately 6% difference.However, as the Reynolds number increases the differences become smaller since the inner layer contributes substantially to the momentum deficit, and "erroneous" outer-layer information does not substantially affect the predictions.
Finally, in spite of highlighting the importance of using a direct method to assess drag as opposed to relying on indirect techniques, it is necessary to recall that the quantitative results of the indirect methods analyses and their predictions of drag penalty at high Reynolds numbers remain tied to the surfaces investigated in this study.The fact that the analysis of both surfaces showed different rates of drag penalty predictions (see Table 5) means that other surfaces will likely produce other figures and in some situations, OLS and CSS methods will remain excellent options to assess the aerodynamic roughness parameters and predict the drag at high Reynolds numbers when direct drag measurements are not available (Monty et al. 2016;Womack et al. 2019).
169 Page 18 of 20 case), and reduces with increasing speed, as reported in Table 6.
The uncertainty in the roughness function is deduced by applying a linear error propagation to equation 1, which is rewritten in its log-law form without the wake extension, as: Assuming all quantities are uncorrelated, the combined uncertainty ΔU + can be expressed as: Equation 14 factors in the most significant error contributions, including the uncertainty in the logarithmic slope ( ), the wall-normal location ( y + ), the zero-plane displacement ( d + ), the log-law intercept ( B ), the streamwise mean velocity ( U +).
Despite an ongoing debate about the universality of the Von-Kármán constants -universality between smooth and rough-wall flows, in the current study, they are assumed to be "universal".Admitting an error of 1 % ( = 0.39 ± 0.004 ) as reported by Marusic et al. (2013), the uncertainty in the loglaw intercept becomes B = 4.5 ± 0.13 meaning B (%) = 3% .The uncertainty in the wall-normal distance ( y + ) depends on the uncertainty in the friction velocity, air viscosity (itself dependent on room pressure and temperature), and PIV grid resolution.The uncertainty of the zero-plane displacement ( d + ) depends on the mean velocity gradient, the friction velocity, air viscosity, and the logarithmic slope.Besides the friction velocity, the uncertainty in the mean streamwise velocity ( U + ) also depends on both systematic and precision errors due to the limitation of the PIV reconstruction algorithm and the finite number of snapshots, respectively.Finally, the roughness function is prone to the fitting procedure, with ΔU + fit essentially being influenced by the extent of the log region and the number of fitted data points.
Once the uncertainty in the roughness function is determined, the error in the equivalent sand-grain roughness height can be evaluated, using either equation 7 for k + s or equation 8 to directly assess k s .A similar linear uncertainty propagation procedure is carried out on equation 7, which yields the following expression: Equation 15 shows the sensitivity coefficients are expected to produce more significant errors given the expression of k s is based on an exponential function.This expression requires the knowledge of the overall uncertainty in the roughness function ΔU + , but also C which stems from the difference between the smooth-rough-wall log-law intercepts.
Finally, to complete the uncertainty analysis, the error in the wake strength parameter is evaluated using equation 12 which yields the following expression:  In this instance, it is important to note that besides the uncertainties , B and ΔU + , equation 16 shows that Π R also depends on the uncertainty in + but more importantly the error in U + ∞ (i.e.C f ).Hence, Π R is inherently expected to produce larger uncertainties.
The results of the different quantities mentioned above are tabulated in 7, for both surface configurations and at all speeds examined in the wind tunnel.The results show that the roughness function is accurately measured to within ±3.5% and with the error systematically dropping down to nearly ±1.5% at the highest speeds.The overall error in the roughness Reynolds number (normalised equivalent sandgrain roughness height) is shown to be an order magnitude larger than the roughness function, with k + s being around ±15% at low speeds, and marginally decreases to around ±10% at the highest Reynolds numbers.
Ultimately, the largest error is observed with the wake strength parameter, which exhibits errors larger than ±20% at low speeds but reduces near or below ±15% .However, it is important to recall these uncertainty estimates do not conflict with our conclusions, and they also do not necessarily imply the wake strength parameters discussed in earlier sections are erroneous.In fact, in addition to using equation 16, we have also determined the uncertainty Π R ,fit by fitting Cole's wake function and the results clearly indicate smaller errors.The error Π R ,fit is shown to be within ±6% or less for Config 1 and ±3% Config 2, respectively.This observation provides an additional note of caution as evidenced by how small errors from drag can propagate and can substantially impact other aerodynamic roughness parameters.Therefore, these need to be carefully considered, given the drag penalty predictions at high Reynolds numbers are shown to be influenced by the equivalent sand-grain roughness height but also the wake strength parameter.

Table 2
(left) Key surface parameters from the scanned coupon, where h ′ represents the surface elevation of the topography about the mean height (h � (x, z) = h(x, z) − h) .(right) The surface topography of a single repeating unit was used for the experiments.Measurements are carried out in two different configurations as shown with Config 1 having a streamwise characteristic wavelength ( S x,1 ) between bar-nacles of 50 mm while Config 2 has a wavelength ( S x,2 ) of 25 mm.Note that S x,1 = 2S z,1 and S x,2 = S z,2 ∕2 .The influence of surface directionality is highlighted in the effective slope ES with ES x,1 and ES x,2 representing the streamwise effective slopes for Config 1 and Config 2, respectively.Hence, the streamwise effective slope of Config 1 is equal to the spanwise effective slope of Config 2 and vice-versa Parameters Value Units Definition

Fig. 1 a
Fig. 1 a Variation of the skin-friction coefficient C f measured using the drag balance with respect to Reynolds number Re x for the two configurations, with the red-filled markers representing the interpolated values where the s/PIV measurements are performed.The black solid line indicates Schlichting's power law for smooth-wall with a ≈ −1∕5 and b ≈ 0.058 .The blue star symbols represent oil-film

Fig. 2 a
Fig. 2 a Variation of the roughness function ΔU + DB as a function k + rms with ΔU + DB determined from drag balance measurements and boundary layer thickness, using equation (1.5).b Variation of ΔU + DB with respect to k + s,DB determined using equation 7 with the constant

Fig. 3 Fig. 4 a
Fig.3Wall-normal distribution of the inner-normalised streamwise velocity profiles for both configurations at Re x ≈ 4 × 10 6 , compared with the DNS turbulent boundary layer smooth-wall ofSillero et al. (2013).The dashed line represents the logarithmic distribution with the constants and B being 0.39 and 4.5 respectively

Fig. 5 a
Fig. 5 a Variation of the roughness function ΔU + against k + s for two configurations determined from the direct drag measurements and velocity profiles.b Inner-normalised streamwise velocity profile with wallnormal distance normalised by k s .The additive constant B r = B + C ≈ 8 for this data

Fig. 6 a
Fig.6a Variation of the wake strength parameter Π R with respect to the friction Reynolds number Re ( + ) estimated from the composite profile, compared with the smooth-wall Π S (solid black line) and the predicted Π R (dashed blue lines) using equation 11.Π S is determined from the DNS smooth-wall profile ofSillero et al. (2013).b Com-

Fig. 7 a
Fig. 7 a Wall-normal distribution of the velocity profiles in defect form for both configurations.b Variation of the inner-normalised streamwise turbulence intensity profiles.The black solid line in both figures represents the DNS smooth-wall data from Sillero et al.(2013).The profiles are also compared with the fully-rough profiles

Fig. 9
Fig. 9 Comparison of the skin-friction coefficient determined from the indirect methods with the estimates from the drag balance for Config 2

Fig. 10
Fig. 10Comparison between the different methods of the predicted full-scale average skin-friction coefficient C F as a function of the fetch Reynolds number Re L using the procedure presented byMonty et al. (2016) for Config 2. The used unit Reynolds number is U ∞ ∕ ≈ 2 × 10 7 with U ∞ ≈ 18 ms −1 and ≈ 9 × 10 −7 m 2 s −1 y + − d + + B − U + .(14)

Table 1
Methods used for the current case study

Table 3
Experimental parameters and flow conditions for both configurations.The PIV and sPIV data refer to the (x, y) and (z, y) measurement planes while the subscripts

Table 5
Assessment

Table 7
Overall uncertainty estimates of the aerodynamic parameters as measured directly using the floating-element drag balance and the velocity profiles from PIV/sPIV, for both surface configurations.The uncertainty