1. Introduction
A portion of the parasitic drag on an aircraft is due to irregularities and deviations from a smooth (clean) external surface. This contribution is commonly termed excrescence drag. The percentage of the total drag due to excrescence drag varies significantly based on the plane design, size, and operation. The study of surface imperfections dates back to the 1930s (e.g., [
1]) and a review of early studies is provided in Young and Paterson [
2]. More recently, studies have been performed to better understand excrescence effects on transition and the use of flow-control techniques to mitigate their impact [
3,
4,
5,
6]. A common source of excrescence drag is the creation of a step due to a mismatch in height between two surfaces. Steps are frequently formed at skin joints, around windows, and control surfaces. In addition, a backward-facing step (BFS) is frequently created when protective films are added to the leading edge of airfoils for various applications including wind turbine and helicopter blades. While these protective films do protect the integrity of the leading edge, the backward-facing step created by the film is detrimental to the performance [
7,
8,
9]. The issue of a backward-facing step on airfoils is not new and has been studied extensively [
10,
11,
12,
13]. The current study is a preliminary examination of a serrated (i.e., saw-tooth) pattern being added to the downstream edge of a BFS as a means of mitigating its adverse effects on airfoil aerodynamic performance.
Before the previous study on the sBFS is presented, low-profile flow-control methods will be discussed. A common passive flow-control device is a vortex generator (VG), which is commonly used on aircraft to control the wing aerodynamic boundary layer. It does so by re-energizing the boundary layer to delay flow separation and aerodynamic stalling [
14]. Some VGs are designed for specific flow regimes. For example, Yao et al. show that low-profile vortex generators outperform traditional VGs at high angles of attack [
15]. The Gurney flap can significantly increase lift at low angles of attack without any penalty in terms of drag [
16]. While the advantages of VGs are plenty, passive flow-control methods such as these require tuning of the geometry with the flow conditions, typically resulting in a relatively narrow range where they are beneficial, which is why they are generally designed for specific flow regimes. An extensive review of low-profile VGs is presented by Lin [
17]. Additionally, the micro vortex generators that are used in supersonic regimes improve shockwave-induced boundary layer separation [
18], or to eliminate the flap separation of take-off and landing configurations of aircrafts and increase the lift-to-drag ratio at various angles of attack [
19]. Since these micro VGs are mostly used in supersonic regimes, the literature on subsonic regimes is scarce.
There is evidence that the addition of various patterns to the trailing edge of a BFS is an effective passive flow-control technique for mitigating the adverse effects (e.g., [
20]) and even reduction in noise [
21]. Serration (i.e., making a saw-tooth pattern) in the BFS is a pattern that has shown potential for mitigating the BFS impact. For the current study, a serrated BFS (sBFS) was created by applying a thin film around the leading edge of an airfoil and then cutting triangular shapes into the downstream edge that would form the BFS. This mimics the use of protective leading edge tapes or coatings used for helicopter and wind turbine blades. Protective tapes without the serration have previously been studied and shown to protect the rotor blades while significantly increasing the profile drag [
22]. Qualitative investigations adding the serration suggested that it could mitigate the added drag. The National Renewable Energy Laboratory (NREL) performed a comparative study of the performance of wind turbines with a serrated protective film relative to a non-serrated protective film [
23]. The turbine was a two-bladed Westinghouse WWG-0600 with a nominal rotor diameter of 44 m. Preliminary data from that study show the KWH power production increased with the sBFS relative to the BFS by 6% and 7.4% at wind speeds of 7 m/s and 11 m/s, respectively. However, this study only looked at the bulk performance (i.e., turbine efficiency) and no information was provided about the impact of the aerodynamic performance due to the sBFS. Consequently, there is little insight into the relationship between the sBFS shape and flow parameters, which is required to develop an optimized design for a given application.
Most reported aerodynamic work on sBFS has been related to their application on airplane wings. Commercial flight testing with sBFS (termed a conformal vortex generator or CVG) on a Boeing 737 was performed by Safair, which reportedly showed a 1% reduction in fuel consumption [
24]. The first assessment of the aerodynamic impact of the sBFS on the Boeing 737 suggested that shock stabilization on the wing could be a cause for the reduced fuel consumption [
25]. However, flight tests on a subsonic aircraft (Piper Cherokee) that provided wall shear distribution on the wings indicated that sBFS redistributes the near-wall momentum. This was apparent as a low-shear diamond pattern formed downstream of the sBFS valleys (i.e., the furthest upstream step location) [
26]. A computational model was developed to further investigate the induced flow pattern observed on the subsonic aircraft [
24]. The model confirmed the low-shear diamond pattern from the subsonic flight tests. Further evaluation computationally [
10] and experimentally [
26] showed that these diamonds were created by the sBFS peak (i.e., furthest downstream step location), inducing the laminar-to-turbulent transition while delaying the transition downstream of the valley. In addition, these studies showed that low-profile (i.e., step height was approximately 5% of the boundary layer height) sBFS could produce strong coherent structures that persisted downstream with a height that was comparable to the boundary layer thickness.
The first part of the current study experimentally investigates the sBFS applied to an airfoil (LA203A) mounted in a wind tunnel and operated at subsonic conditions. The aim was to quantify the aerodynamic impact of the sBFS relative to a BFS as well as a clean (i.e., no BFS) wing. The total drag was quantified from wake surveys in the far-wake region. We present evidence that the sBFS-induced coherent structures that persisted into the far-wake region led to further examination of the induced vorticity and its dependence on sBFS size (length and width) for a subset of conditions. All results are compared relative to the BFS and clean conditions. Then, a computational model was created and validated against the available experimental results to further examine the flow behavior on the airfoil that produced the coherent structures. The remainder of the article describes the experimental and computational methods used in
Section 2;
Section 3 presents drag coefficient and vorticity results from both the experiment and computations;
Section 4 provides a qualitative and quantitative discussion of the sBFS’s performance; and
Section 5 provides a brief summary and concluding remarks.
2. Methods
2.1. Experimental Methods
2.1.1. Test Facility and Model
The experiments were conducted in the custom built Flexible-use Wind Tunnel (Diehl Aero-Nautical) located at Oklahoma State University [
27,
28,
29]. It was an open loop, draw down tunnel powered with a 125 hp centrifugal fan, which can achieve speeds of up to 30.5 m/s (100 ft/s). The tunnel test section can be varied, and, for the current study, it was 2.44 m (8 ft) long with a 0.91 m (3 ft) square cross section. The freestream turbulent intensity within the test section was measured to be below 1.4% [
28]. A schematic of the wind tunnel test section is provided in
Figure 1 with the test model located 0.762 m downstream from the test section inlet and vertically centered. The test model was an airfoil with an LA203A profile and a chord length (
c) of 0.197 m (7.75 in). The model airfoil angle of attack (defined in
Figure 2a) was then varied for a given test speed.
2.1.2. Backward-Facing Steps
A film (i.e., thin tape) was applied around the leading edge of the airfoil model and extended downstream to the desired chord location. If a BFS was being tested, the end of the film sheet was cut straight, while for sBFS, a serrated pattern, as illustrated in
Figure 2b, would be cut into the film. As previously mentioned, the upstream point and downstream point of the sBFS are referred to as the valley and peak locations, respectively. The height of the BFS and sBFS was fixed at 1.02 mm. The foil thickness varied with a maximum thickness of 15.6% of the chord length (30.7 mm). The step height was selected to correspond to that of a Boeing 737 aircraft following the scaling proposed in Lucido et al. [
26], though the current study was subsonic. Four different sBFS configurations were tested, termed sBFS-V1, -V2, -V3, and -V4. The size and position of each sBFS and BFS, as well as their position on the airfoil, are listed in
Table 1.
The sBFS-V1 length (
L = 70 mm), width (
W = 50 mm), and height (
H = 1.02 mm) were three times those tested on the Boeing 737 aircraft, which was selected to scale to flight conditions at cruise altitude [
13]. The sBFS-V1 valley was located at 10% chord length as the flight tested on the Boeing 737. To separate geometric and kinematic dependence, the length and width were halved for sBFS-V2, which produced the same
L/
W ratio without changing boundary thickness or step height. The sBFS-V2 valley location was matched with sBFS-V1 at 10% chord length, but the peak location was changed due to the decreased length. Consequently, sBFS-V3 had the same dimensions as sBFS-V2 but shifted downstream such that its peak location matched the downstream position of sBFS-V1. For the final configuration (sBFS-V4), the length matched that of sBFS-V1, while the width was halved (i.e., matched sBFS-V2 and -V3) to assess the impact of the
L/
W ratio on the sBFS performance. Both sBFS-V3 and -V4 configurations were tested at the same time with half of the airfoil model span having one of the configurations, which was carried out since initial configurations showed minimal spanwise variation for scales larger than the sBFS width. In addition to the sBFS, the BFS was tested with the step height fixed at 1.02 mm and its location at 10%, 27%, or 44% chord lengths, which corresponds to peak and/or valley locations for the sBFS configurations. Pressure measurements at these locations were not recorded for the experiments; however, they could be tabulated using the simulations. The corresponding pressure gradients were as follows: −8 kPa/m at 10% chord length, 2.5 kPa/m at 27% chord length, and 5 kPa/m at 44% chord length.
2.1.3. Instrumentation and Data Analysis
The primary measurements for the current study were wake surveys via a Pitot-static tube mounted five chord lengths downstream of the airfoil model (see
Figure 1). The Pitot-static probe was traversed through the test section cross-section on an automated traverse (Dantec) with each wake survey acquiring data at nominally 100 locations. For all conditions, the streamwise velocity component (
u) was measured, and for a subset of conditions, a 5-hole probe was installed. The 5-hole probe measures the velocity magnitude (
V), angle of attack (α), and the slip angle (β), which can be converted to the streamwise
, spanwise
, and vertical
velocity components.
In addition, another Pitot-static probe was located in the test section inlet to measure the free-stream velocity. The turbulent intensity at the test section inlet was also measured with a hot wire. Finally, the barometric pressure and temperature of the air at the test section inlet were recorded with a data-logger (SD700, Extech Instruments, Nashua, NH, USA), which was used to determine the fluid properties (e.g., density, kinematic viscosity). The signals were sampled at 1000 Hz via an analog data acquisition module (NI-9220, National Instruments, Austin, TX, USA) and chassis (cDAQ-9188, National Instruments, Austin, TX, USA). A commercial software package (LabView, National Instruments, Austin, TX, USA) created a Virtual Instrument (VI) that was used to handle the data with the mean voltage from every 100 samples recorded, which produced an effective sample rate of 10 Hz. The pressure and temperature data-logger did not have an external output, so that data was manually inputted into the VI.
The vertical spacing of the velocity profiles within the wake was 3.18 mm, which was half the diameter of the probe head. Integration of the streamwise velocity is a well-established method for obtaining the total drag on 2D objects that can be extended to 3D models [
30]. Three-dimensional considerations were required due to some of the sBFS producing significant coherent structures that persisted into the far-wake region where the wake survey was acquired. A brief overview of the drag analysis is provided here for completeness as well as its relationship to the measurement uncertainty. Assuming a uniform inlet velocity and nearly horizontal streamlines sufficiently far from the model, integration of the streamwise component of velocity in the wake gives the drag coefficient
as
Here,
is the total drag,
is the planform area,
is the freestream speed,
is the fluid density,
is the chord length,
is the vertical distance, and
is the mean streamwise velocity in the wake. The vertical wake velocity profiles were integrated using the trapezoidal rule to determine
. No correction for blockage was performed since the maximum solid body blockage was 7%, and for the majority of conditions tested it was below 5%. Corrections for measurements within the wake of a streamlined body are minimal prior to separation [
31]. The uncertainty of
calculation was determined using standard propagation of uncertainty procedures. While the uncertainty for instruments and geometric measurements were readily available, the final uncertainty was heavily dependent on the covariance between the mean velocity measurements and the free-stream speed. With a covariance of 0.395, the resulting uncertainty in
was 3.6% at the maximum test speed. Given the uncertainty of the covariance and weak variation between operation conditions, error bars on
plots have been fixed at 5% for reference.
Due to the observation of coherent structures in the far-wake region for a subset of conditions, the induced vorticity within the wake was computed as a measure of the strength of the coherent structures produced by the sBFS. Vorticity is defined as the curl of the flow velocity vector,
Since the traverse system could not move in the streamwise direction, none of the streamwise gradients could be computed experimentally. Preliminary results from the computational study indicate that the streamwise component of vorticity
is dominant [
19]. This is also consistent with the expectation that a coherent structure persisting to the far-field would have significantly stronger gradients in the cross-stream directions relative to the streamwise. Note that since the vorticity is dependent on resolving the velocity vector, only the subset of conditions that used the 5-hole probe could analyze the vorticity. The streamwise vorticity was calculated at each measurement location and then averaged to create vorticity profiles of the wake deficit for each test configuration.
2.1.4. Test Conditions
The airfoil configurations were tested at free-stream speeds of 20 m/s and 26 m/s, which have a corresponding chord length-based Reynolds number
of nominally
and
, respectively. The current study focuses on the
results since no significant variation was observed between test speeds. This was expected due to the relatively narrow Reynolds number range but tested to confirm the insensitivity to Reynolds number over this range. For each test configuration, the angle of attack (α) was typically varied between 0 and 10 degrees in 2-degree increments. The test model configurations were clean wing (i.e., no film), three BFS, and four sBFS (see
Table 1 for geometries and locations). For all configurations, the vertical increment (∆
z) for the wake survey was 3.18 mm, which was half of the probe head diameter. For the clean wing and 7 BFS configurations, the spanwise increment (∆
y) was increased to 25 mm since there was negligible spanwise variation. Conversely, the sBFS configurations produced coherent structures that created spanwise variations; the spanwise increment was reduced to 12.5 mm. In addition, the measurement locations with the sBFS were selected such that measurements were acquired directly downstream of peaks and valleys. The subset of conditions tested with the 5-hole probe were clean wing, BFS at 10% chord, BFS at 44% chord, and sBFS-V1.
2.2. Computational Methods
STAR-CCM+ v2022.1.1, a commercially available CFD software, was used for the computational study. The complex 3D models were created in SolidWorks and imported to STAR-CCM+, which has its own built in 3D-CAD, where the models were further edited to facilitate the simulation. STAR-CCM+ provides built-in automated mesh capabilities that include tetrahedral, trimmer, and polyhedral unstructured meshes. These mesh algorithms also have the capability of utilizing prism layers, which is essential in resolving the near-wall region of a boundary layer. STAR-CCM+ also has the capability to run a large range of physics- and flow-specific models that include RANS, DES, LES, and qDNS, with sub models for each. For the current study, RANS, or Reynolds-Averaged Navier–Stokes equations, which are often used to study the mean characteristics of the flow, were chosen.
The experimental setup was replicated in the computational model. A 3D rectangular domain was created for all simulations. The side walls were given periodic boundary conditions, while the top and bottom walls were subjected to slip conditions. The inlet was a velocity inlet, while the outlet was treated as a pressure outlet. The side walls had periodic boundary conditions. The domain total length was 3.75 m, width 15.24 (width of 3 sBFS), and the height was 0.91 m to match the height of the wind tunnel in the experimental study. In this study, the span of the airfoil was reduced to 3 sBFS widths (152.4 mm). The sBFS-V1 configuration was chosen since it provided the most significant wake structures. The length, width, and thickness were 69.9 mm, 50.8 mm, and 1.1 mm, respectively, with the valley and peak locations at 10% and 45%, respectively. The turbulence model used for this study was k-omega turbulence and the turbulent intensity at the inlet was set to 1%. Surface, polyhedral and prism layer meshers were chosen for this study. The mesh was compared by studying 3 independent meshes comprising 2.5, 14, and 25 million cells. All 3 meshes produced similar results and the one with the highest number of cells was used for all of the analysis.
4. Discussion
The lowest drag results for the BFS (i.e., 45% chord) and sBFS are compared along with the clean wing results in
Figure 12. As expected, the sBFS did not reduce the drag relative to the clean wing, though at the highest angle of attack tested, the two configurations are within the measurement uncertainty. Conversely, sBFS-V3 had lower drag than BFS at 45% chord for nearly all angles of attack tested, and at zero angle of attack, the increase was comparable to the measurement uncertainty. This was in spite of the fact that sBFS-V3 had the valley location (27% chord) closer to the leading edge than the BFS at 45% chord. Note that the peak location of sBFS-V3 was at 45%, which matches the BFS location compared in
Figure 12. As previously discussed, airfoils regularly have BFSs such as that created by the slat step on the wing on the Boeing 737, which is nominally located at the 10% chord length. These results indicate that the sBFS can mitigate the negative drag effect created by BFS. These findings are consistent with recent work studying microsteps, which have shown that they could have a significant impact on boundary layer flows at scales much larger than their step height [
37,
38]. It should also be noted that minimal effort has been dedicated to date on the sizing of the sBFS relative to the flow parameters, which is critical, as most passive flow-control methods require a tuning of the geometry given the flow parameters.
The computational study showed evidence that these coherent structures are formed by the sBFS redirecting the near-wall flow at the valley towards the peaks. Then, downstream of the valley, the flow remains attached even when the flow separates downstream of the sBFS peak. The pattern then persists downstream, forming the spanwise variation observed in the far-wake region. In addition, the streamwise vorticity in the computational study showed that the coherent structures in the wake are counter-rotating pairs.
5. Conclusions
Wake surveys in the far-wake region of a LA203A airfoil model were used to measure the drag associated with the application of a protective film around the model’s leading edge. A BFS was created on the model with a thin film. The use of a serrated (i.e., saw-tooth) pattern on the trailing edge of the BFS was explored to investigate its potential for mitigating negative impacts. For reference, the clean wing (i.e., no film) and three standard BFSs were tested. These results showed that the BFS increased the drag relative to the clean wing, and that the drag increased the closer the BFS was to the leading edge. Four sBFS configurations, each with the same step height, were tested. In general, over most angles of attack, the sBFS produced lower drag relative to the BFS located at the sBFS valley location. In addition, vorticity and contour plots showed that the sBFS produced strong coherent structures that could persist into the far-wake region.
Computational simulations showed that these structures were produced by the sBFS, with the flow near the surface of the airfoil being entrained in the direction of the sBFS peak. This induced motion modified the downstream flow structure, which, in the sBFS condition computationally examined, impacted separation at the trailing edge. This pattern at the trailing edge likely produced such a strong pattern in the far-wake region, though they were induced by the sBFS as evidenced by the structure alignment with the sBFS peak and valley geometry. The strength of the structures does appear to be dependent, in part, on the length-to-width ratio of the sBFS; however, the possible influence of the width itself was not studied and should be explored further in connection with the impact on the trailing edge separation.
Overall, there is evidence that sBFS could mitigate the adverse effects on drag associated with BFS on airfoils, which are frequently created due to the addition of protective leading edge films or mismatch in height between adjoining surfaces. Mitigating these effects could lead to improved performance across certain flight regimes, and improve the efficiency of wind turbines. The next steps are to develop an LES-based simulation to enable a study of the coherent structures on the airfoil and how they are modified by the sBFS; these insights will enable the tuning of the step geometry based on the operation condition.