1. Introduction
Improving the aerodynamics of aircraft has been a continuous goal, and flow control efficiency, whose essence is to trigger local or global flow-field changes by applying physical quantities, such as force, mass, and energy, to the local flow and using the hydrodynamic interactions between the fluids, makes it stand out from other methods [
1]. According to whether external energy is required during the control, it is mainly divided into two categories: active or passive flow control [
2]. Active flow control techniques are used to improve experimental model performance by directly injecting suitable perturbation patterns into the flow environment. To date, there have been many studies conducted using active flow control techniques to improve aircraft performance [
3,
4,
5]. However, the system of active flow control is very complex, especially in terms of “robustness” and extra energy consumption. In actual flight, plasma, loop control, and other active flow control technologies that require complex additional equipment have not been applied on a large scale in aircraft, with most remaining in the wind tunnel test stage [
6]. Passive flow control does not require auxiliary energy input, and it mainly achieves the purpose of flow control and increasing lift and reducing drag by changing the flow environment, including the boundary conditions, pressure gradient, and other factors. Common passive flow control techniques include leading-edge flaps [
7], vortex generators [
8], winglets [
9], etc. Compared with active flow control, passive flow control is easier to implement and maintain; therefore, it is more suitable for engineering applications with more complex operating environments. Exploring new flow control technology to enhance the maximum exploitation of the aerodynamic performance potential of aircraft is of considerable academic research value and engineering significance.
Before the Langley Research Center found that a small rib surface could effectively reduce the downstream friction resistance of walls in 1978, it was generally believed that the smoother the surface, the lower the resistance [
10]. Researchers have conducted much of research on non-smooth surfaces and have developed many passive drag-reduction control technologies, such as rough wall surfaces [
11], superhydrophobic wall surfaces [
12], soft skins [
13], etc. Inspired by the unique feather structure of “silent flight” owls [
14], porous media have been extensively used in the field of drag reduction in recent years. Rasheed et al. carried out a boundary-layer control study on a cone with a half-cone angle of 5.06 degrees under hypersonic conditions, and their experimental results showed that porous media can effectively reduce frictional resistance compared with the experimental results without porous media; furthermore, the conditions with porous media on the wall can even delay the occurrence of transition [
15]. Venkataraman and Bottaro performed numerical simulations and found that porous media affect flow topology near the rear of the wing by spontaneously adapting to the separation flow [
16]. In 2017, Klausmann and Ruck laid porous media on the downstream side of a cylinder and measured the drag, and their results showed that porous material reduced drag by 13.2% [
17]. Joshi and Gujarathi determined that the main cause of porous media drag reduction was the Darcy flow, which transformed the no-slip boundary condition to a quasi-slip Fourier-type and changed the shear force [
18]. Li et al. studied an open-channel flow covered by a porous wall with reduced spanwise permeability through direct numerical simulations. They found that the maximum drag-reduction effect of 15.3% occurred when the depth of the porous layer was nine viscous units [
19]. Liu et al. summarized that the porous media not only has the effect of noise reduction but also can change the resistance, and by increasing the porosity, the resistance is subsequently reduced. When the porosity is close to 0.97, the highest resistance-reduction effect can be achieved [
20].
Mößner et al. installed various porous media in a wind tunnel model and examined them using pressure and PIV measurements. The results show a decrease in lift with increasing permeability. The PIV data show that the airfoil’s suction-side boundary-layer thickness increases as a consequence of the flow through the porous trailing edge. Additionally, turbulence increases together with permeability [
21]. Aldheeb et al. presented a study on porosity’s effect (honeycomb aluminum) on the aerodynamic performance of a symmetric thin airfoil and a straight half-wing based on a thin symmetric airfoil section. They found that drag and moment coefficients change as porosity changes. The drag increases due to increases in viscous stresses caused by flow through the porous region. The lift and lift slope are reduced due to a reduction in the pressure difference between the airfoil’s upper and lower surfaces [
22]. Thus, different porosity types produce different aerodynamic behaviors. At very low values of porosity, the drag was reduced at low angles of attack. In 2019, Tamaro found that the power spectral density of pressure close to the wall was reduced in the whole frequency spectrum compared with a solid case in correspondence with a porous trailing edge, while the wake structure was highly affected by the presence of porous media [
23]. In 2021, Tamaro et al. studied the flow field on solid and porous airfoils (a permeable exoskeleton) subjected to turbulence shed by an upstream cylindrical rod through PIV. Additionally, their flow-field investigation showed that porosity’s main effect is to mitigate the turbulent kinetic energy in the stagnation region, attenuating the distortion of turbulence interacting with the airfoil surface [
24]. In 2022, Du et al. found that porous media primarily reduced wall resistance because of its “micro-jet” effect, which inhibits the vortex separation frequency and reduces the wall turbulence kinetic energy [
25,
26]. Although developments in drag-reduction control methods have made considerable progress in porous media research, most studies focus on simple models, such as a flat plate or cylinder, while little research has been conducted on the wing surface, which is very sensitive to aerodynamic characteristics. The published studies only describe qualitatively that porous media has a certain control effect on wing performance, but the introduction of the influence law on porosity is relatively vague; in particular, its influence mechanisms are still unclear.
Based on this research background, our paper summarizes the effect laws of porous media parameters on the aerodynamic characteristics of airfoils. We also analyze the mechanism of influence on the surface flow-field structure through balance force and pressure measurements and using a TR-PIV experiment that includes laying porous media with different PPI on the airfoil’s suction side using a wind tunnel’s low-turbulence experimental platform. Our paper’s structural framework is as follows:
Section 1 introduces the research background of this paper;
Section 2 introduces the experimental model and conditions;
Section 3 discusses aerodynamic laws;
Section 4 analyzes the flow mechanism from the flow-field test results; and
Section 5 summarizes the research content and plans future research.
2. Experimental Scheme
2.1. Wind Tunnel Test System and Experimental Model
The experiment was carried out in the low-turbulence recirculation wind tunnel of the Nanjing University of Aeronautics and Astronautics. The wind tunnel test section measured 3.0 m along the streamwise direction and had a rectangular cross-sectional area of 1.5 m × 1.0 m. The wind tunnel’s turbulence intensity, the airflow’s deviation angle, and dynamic pressure stability coefficient were less than 0.5%, 0.2°, and 0.05, respectively, and the design’s stable wind speed ranged from 5 to 100 m/s. The test section window was equipped with optical glass to provide high-level imaging. The wind tunnel’s experimental system, which mainly involves an aerodynamic-balance force measurement system, a pressure testing system, and a TR-PIV system, is shown in
Figure 1.
We fixed the airfoil model on the balance using the support mechanism, and then we bolted the balance to the base’s attack angle mechanism, while the piezometric tube, which we pre-buried inside the model, penetrated through the round hole on the lower wall of the wind tunnel. We emitted the laser plane parallel to the incoming flow direction from downstream to upstream during the PIV experiment; in this way, we significantly reduced the reflection interference caused by the laser incident on the model’s surface. The center of the laser plane coincided with the model’s midline, while the camera was fixed at the external upper wall of the wind tunnel, with its lens perpendicular to the laser plane. Moreover, the velocity of the incoming flow was 10 m/s for this experiment, the Reynolds number was, the inlet total pressure was about 101,386 Pa, and the temperature was 288 K.
We adopted the Cartesian coordinate system, and the flow, normal, and spreading direction served as the x, y and z axes, respectively, in which
was denoted as the bottom end of the model. We chose the SD8020 airfoil for our model, which was made of Q235 steel, to ensure strength during the flow. Considering the size of the test section, we designed the chord length and wingspan to be 200 mm and
, respectively. In total, there were 25 static pressure holes of 1 mm in outside diameter on the upper and lower surfaces, which we led via a purple copper tube (wall thickness of 0.15 mm, inner diameter of 0.7 mm) from the interior and connected to the pressure-scanning valve by a plastic hose in the terminal. Furthermore, we set 14 and 10 static pressure taps on the upper and lower wing surfaces, respectively, and arranged 1 hole on the leading edge facing the incoming flow, while setting a 5° angle between the holes’ direction and the chord to avoid mutual interference. We placed slots at a position 0.2 c from the leading edge, where the maximum depth and size were 10 mm and 40 mm × 980 mm, respectively. We laid same-size porous media in the slot to ensure a consistent shape with the airfoil’s outer profile. The model design is shown in
Figure 2, and
Table 1 shows the pressure taps’ pressure distributions.
We used copper foam as the porous media, which is a multifunctional material with a large number of connected or unconnected holes uniformly distributed in the copper matrix and is commonly used for thermal conductivity, pressure buffering, and noise reduction. This material can stabilize vortices, change wall structures, and has excellent drag-reduction performance due to its unique inner structure. In addition, the copper foam material is soft and easy to cut, ensuring a perfect fit with the airfoil’s surface. We investigated three copper foam porous pads with 5 PPI, 20 PPI, and 60 PPI, whose permeabilities were 6.27 × 10
−7 mm
2, 1.31 × 10
−7 mm
2, and 1.24 × 10
−8 mm
2, respectively, with 98% porosity in all cases. Images of the porous media are shown in
Figure 3. The three selected porous media had large pore density spans and considerable permeability changes; therefore, we could obtain the general control rule through our experiment.
2.2. Experimental Equipment and Test Methods
The force measurement experiment used a six-component cassette strain balance, which separately measured six components (X, Y, Z, Mx, My, and Mz) through mechanical decomposition with a bridge circuit. The balance response frequency was greater than 80 Hz, the load range of the balance in the z direction was 50 kg, the load range in x direction was 10 kg of the resistance direction, the load range in y direction was 15 kg of the lift direction, and the calibration accuracy and precision were 0.28% and 0.041%, respectively, fully meeting the requirements of our experiment. In the formal trial, the sampling frequency of the force measurement was 2000 Hz with a sampling time of 5 s, which means we obtained the final force data from 10,000 points of time-averaged results. The airfoil lift and drag coefficient can be calculated by Equations (1) and (2), respectively.
where
is the incoming flow velocity,
is the projected area of the airfoil,
is the spread length of the airfoil, and
are the total lift and drag forces gauged by the balance, respectively.
In the experiment, we used two DSA3217-PTP/16Px with 16-channel pressure-scanning valves to gauge the static pressure on the airfoil’s surface; moreover, the pressure-scan valve’s sampling frequency was 1000 Hz, with a sampling time of 5 s. The pressure-scanning valve we used had a water-column-pressure measurement range of 10 inches, which is approximately 2540 Pa, with an accuracy within 0.05% of the full scale. The surface pressure coefficient (
) distribution can be calculated by Equation (3).
where
is the static pressure measured by the pressure-scanning valve.
We used the TR-PIV to test the flow field on the airfoil’s suction side, which consists of a double-pulse laser (Vlite-Hi-527-30, divergence angle ≤ 5 mrad), a PCO high-speed camera (2000 pixel × 2000 pixel, 12 bit), a synchronization controller, and a smoke generator. Before the experiment, we sprayed a layer of black matte paint with a thickness of approximately 20 μm on the wing model’s surface to absorb the laser and further eliminate reflection issues. In the PIV experiment, we set the laser’s luminescence frequency of to 1000 Hz, the pulse width of both A and B lasers to 50 μs, and the interval time to 60 μs. The laser plane’s intensity was uniform in the flow-field test area, with a thickness of approximately 1 mm. We equipped the camera lens with a 750 nm high-pass filter, set the camera’s shooting frame rate to 1000 fps, activated a double exposure mode, and assigned an 80 μs exposure time, while controlling the laser’s light source and the camera’s time sequence using a synchronization controller to filter out the excitation and ambient light, in addition to other sources of stray light interference. The wind tunnel had two experimental sections. Because the wind tunnel is backflowing, we placed the smoke generator upstream of the model (at another test section) to spray tracer particles through the wind tunnel observation window to observe the tracer particle situation, after which we uniformly distributed the concentration to start recording the flow-field picture. The trace particles, which were generated from smoke-oil vaporing at sizes of 1μm, displayed high followability and had negligible impact on the flow field. Before post-processing, we filtered 1000 pairs of transient raw images using the MATLAB filtering program, and then we used the PIV-view software to perform correlation analysis on the processed images. We used the standard FFT correlation algorithm with an interrogation window of 32 pixel × 32 pixel and an overlap factor of 50% to determine the correlation and thus obtain the particle displacement.
2.3. Ω Vortex Identification Method
The vorticity in the flow field is defined as the spin of the velocity, which is of great value in flows dominated by vortex structures. However, the magnitude of vorticity is not able to precisely represent the existence of an actual vortex. For example, there are vortices in the laminar boundary layer, but no rotational motion actually occurs. Liu et al. [
27] proposed a method for vortex identification: the vortex volume is further decomposed into a spin and spinless part with the introduction of the parameter
Ω, which is to define and identify the vortex, as shown in Equation (4).
where
Ω ∈ [0, 1] represents the ratio of the vortex vorticity of the total value, and the flow can be classified into translational, deformation, rotational, and non-rotational motions according to this way of decomposition.
Ω = 0 represents pure deformation motion, and
Ω = 1 means pure rotational motion. Based on the theory of Liu et al. [
28], the threshold value of
Ω is generally taken as 0.52, and the vortex core boundary selection of the threshold
Ω 0.95.
4. Flow-Field Results and Discussion
Through force and pressure analysis, we found that the porous media have a considerable effect on the airfoil’s aerodynamic performance in the entire AOA test range. Therefore, in order to improve efficiency during the following experiment, we decided to select two typical AOA for the PIV test (a small and a large AOA). The main reason for choosing the 6° AOA was that the lift–drag ratio reached its peak in this situation, while the porous media’s control effect continued to perform at a high level. Furthermore, we chose the 12° AOA mainly because of our aim to explore the control mechanisms of porous media on airfoil aerodynamics after stall, while the flow structure was obviously revealed in the separation zone under this AOA, thus making it convenient for comparison.
In
Section 4.1, we first discuss the effect of pore densities on time-averaged velocity and shear stress fields. Subsequently, we undertake an analysis of the unsteady flow field in
Section 4.2. For the sake of simplification, we examine the development law of a vortex under the highest control condition compared with the smooth airfoil, which we laid with porous media of 20 PPI and while positioning a 12° angle of attack. Finally, in
Section 4.3, we analyze the vortex’s spatial–temporal evolution processes of the transient flow field using the dynamic mode-decomposition method (DMD) to analyze the mechanisms inherent in the changes in the aerodynamic properties of porous media.
4.1. Time-Averaged Flow Field
4.1.1. Time-Averaged Velocity Field
The time-averaged velocity field describes the speed of airflow movement in a period of time, which reflects the overall flow state and facilitates the visual determination of flow characteristics. It can be obtained by averaging the transient velocities in a region over a period of time.
where
is the number of sampling points in the measurement time
;
and
are acquired from the PIV experiment data.
The time-averaged velocity field contours of the suction side are given in
Figure 7 and
Figure 8, which show models laid with different PPI values of porous media while under the 6° and 12° angles of attack. From
Figure 7, we see that the flow velocity of the airfoil laid with porous media shows a decreasing trend at a small angle (AOA = 6°) compared with the smooth airfoil. As the PPI value decreases, the low-speed area at the tail increases considerably. For example, when using 5 PPI, the trailing velocity drops considerably, indicating that the porous media somewhat hinder the airflow movement, thus reducing the circumfluence circulation and leading to reduction in the lift coefficient. From
Figure 8, we observed that the airfoil’s suction side separation region is expanded at a large angle of attack (AOA = 12°), likely due to the comparatively rough surface that characterizes porous media. As the PPI value reduces, the separation region gradually increases. Meanwhile, the vortex size in the separation region when laying 60 PPI porous media increases slightly compared with that of the smooth airfoil, while the separation vortex size when laying 20 PPI porous media greatly decreases. This indicates that although 20 PPI porous media has a certain loss of lift, it can weaken the vortex energy and inhibit large-scale vortices to a certain extent, thus improving the airfoil’s aerodynamic performance.
We captured and identifies the vortex characteristics within the separation region on the suction by means of the
Ω vortex identification method mentioned above. It can be seen from
Figure 9 that the boundary of the vortex core changes after laying the porous media. Compared with the smooth airfoil, the boundary of the vortex core in the airfoil separation area with 60 PPI porous media increases, which may lead to the increase of friction resistance. But the boundary of the vortex core in the airfoil separation area with 5 PPI and 20 PPI porous media is significantly reduced, and the effect caused by 20 PPI is the most significant, greatly improving the aerodynamics. This conclusion is consistent with
Figure 8.
4.1.2. Time-Averaged Shear Stress Field
Figure 10 presents the time-averaged shear stress contours at the 6° angle of attack, which consists of the results of 5 PPI, 20 PPI, and 60 PPI porous media laid on both the airfoil and smooth airfoil. Possibly due to the existence of experimental interference factors, the flow field on the test model’s suction side was not completely attached at the 6° angle of attack and still had a small range of separation. As seen from the figure, the surface’s shear layer is thickened due to the existence of porous media, which suppresses the airflow’s velocity on the upper airfoil. However, although the shear layer of the airfoil laid with 20 PPI porous media is thicker than that of the smooth airfoil, the wall shear force is considerably lower, thus reducing the wall friction and improving aerodynamic performance. The reason may be that the size of the 20 PPI porous media pore is similar to the wall’s spanwise vortex, which leads the spanwise vortex to be “locked” in the porous media pore. Because the rolling friction is much less than the sliding friction, other vortices flow smoothly across the wall, reducing the momentum exchange of low/high speed fluid inside/outside the shear layer, thus reducing the shear stress between the fluid and airfoil.
4.2. Unsteady Flow Field
We selected a 12° AOA for unsteady flow-field analysis because the airfoil’s flow-field structure is not obvious and is difficult to capture at a small angle of attack, while at a 12° AOA, large flow separation occurs, and the vortex structure in the separation area is obvious, which is convenient for comparison. The second reason is that we also want to study the control mechanisms of porous media on airfoil aerodynamic performance and the spatio-temporal evolution law of vortices at a deep-stall angle of attack.
Figure 11 gives three consecutive instantaneous vortex contours based on the
Ω-vortex identification criterion at 12° angles of attack for the airfoil laid with 20 PPI porous media and the clean airfoil. From the figure, we observe that the number of vortices increased on the suction side laid with 20 PPI porous media, the vortex size tended to be finer, and the
Ω value was reduced compared with the clean airfoil, indicating that the proportion of rotating vortices in this region was relatively small. Meanwhile, porous media can break the large-scale high-energy vortex into more small-scale vortices of low energy, accelerating vortex dissipation and achieving improvements in aerodynamic characteristics.
4.3. DMD Mode Analysis
The DMD method is a decomposition technique that can reduce the dimensionality of the dynamic system. The accurate identification of the vortex by modal step-down can capture the vortex frequency, wavelength, and propagation velocity of each mode, and tracking the development and evolution of the vortex structure can result in more information about the dynamics of the original system. Each mode obtained with the DMD method has a unique growth rate and frequency magnitude, which facilitate the analysis of complex flows [
29,
30,
31].
4.3.1. Mode Distribution
Modes can be classified according to their types, such as quasi-static, drifting, and conjugate modes, where the quasi-static mode, also called the first order mode, is similar to the average flow field.
Table 2 gives the first 15 order modes of the airfoil’s suction side vortex field for the airfoil laid with porous media of 20 PPI and the smooth airfoil at an angle of attack of 12°, classified according to the type of mode, noting that the modes under each type are ordered by energy magnitude.
4.3.2. Vortex–Mode Energy
Figure 12 shows the relative energy distribution of the first 15 order modes of the vorticity field on the suction side of the airfoil laid with porous media of 20 PPI and the smooth airfoil at an angle of attack of 12°. From the figure, we observed that the first six order modes’ energy accounts for approximately 70% of the total energy, which indicates that the DMD method used in this test is consistent with the physical characteristics, and its reduced-order reconstruction results are better. In comparison, we found that the airfoil’s first-order mode energy when laid with porous media was less than that of the smooth airfoil, and the first-order modal energy was similar to the average flow field, indicating that porous media can reduce the vortex energy of the average vortex field on the upper airfoil. At the same time, we found that the slope of the first nine order mode curves of the airfoil laid with porous media is considerably larger than that of the smooth airfoil, indicating that the porous media accelerates the main modal vortex dissipation rate.
Figure 13 provides the DMD amplitude and frequency distribution of the first 15 order modes of the vortex field on the smooth airfoil and the airfoil’s suction side with a 12° angle of attack for the airfoil laid with 20 PPI porous media. Because the quasi-static mode’s flow field does not change with time, it does not grow or decay, and its frequency is zero, while the drifting mode eigenvalue is zero in the imaginary part; therefore, its frequency is also zero. From the figure, we observed that the DMD amplitude of the first-order mode of the airfoil laid with porous media was considerably lower than that of the smooth airfoil, which corresponds to the results in
Figure 12. We also discovered that the conjugate mode appeared earlier after laying the porous media compared with the smooth airfoil, indicating that the porous media sped up the mode evolution process.
Comparing the mode frequencies, we found that the frequency increased with the increase in the mode order, indicating that the vortex energy gradually dissipated, and that the low-order-mode vortex had a large scale, high energy, and low frequency, while the high-order-mode vortex had a small size, low energy, and high frequency. After laying the porous media, the mode frequency considerably increased compared with the smooth airfoil, which indicated that the porous media can break the large-scale vortex structure into more small-scale vortex structures, effectively weakening the energy of different modes and at the same time making the high-order modes stabilize more quickly by absorbing high-frequency disturbances, accelerating the vortex evolution process, and thus improving the airfoil’s aerodynamic performance.
4.3.3. Conjugate Mode of the Vorticity Field
Figure 14 shows the first three pairs of conjugate modes of the airfoil’s suction side vortex volume field modal amplitude ordering for the airfoil laid with 20 PPI porous media and the smooth airfoil at the 12° angle of attack. Only one pair of conjugate modes is shown in the figure because the structure of the pair of conjugate modal flow fields is the same.
We found that the vorticity gradually decreased with the increase in the mode order, indicating that the vortex gradually lost energy during the evolution process. Comparing the results in
Figure 14, we discerned that the vorticity values of the corresponding modes were lower, with the vortex size relatively smaller after laying with the porous media, indicating that the vortex strength was weaker. At the eighth and ninth modes, the smooth airfoil’s vorticity was obviously distinguishable, while the vorticity of the airfoil laid with porous media was almost dissipated, further indicating that the porous media promoted the vortex evolution process.