• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Experimental investigation of typical surface treatment effect on velocity fluctuations in turbulent flow around an airfoil

    2023-11-10 02:15:48BkELHAMIMojtDEHGHANMANSHADIAliAkrDEHGHAN
    CHINESE JOURNAL OF AERONAUTICS 2023年10期

    Bk ELHAMI,Mojt DEHGHAN MANSHADI,*,Ali Akr DEHGHAN

    a Faculty of Mechanical Engineering, Malek Ashtar University of Technology, Isfahan 177415875, Iran

    b Mechanical Engineering Faculty, Yazd University, Yazd 8915818411, Iran

    KEYWORDS

    Abstract Noise reduction of different airfoils is important because these sections are used in wind turbines,propellers,and aircraft wings.Several methods are used for passive noise reduction of sections.One of these methods is the use of surface treatment.In this research,the effect of the typical surface treatment element (finlets) on the vortex structure at different frequencies in the turbulent flow created on the NACA2412 section is investigated.For this purpose,one-dimensional hot wire probe is used.The used surface treatment in this research has a special geometry and the distance between two consecutive finlets is 6 mm (S=6 mm).This study shows that this surface treatment element can be used for noise reduction in high frequencies.Another result of this research is the suggestion of the most suitable position to install this special surface treatment element to reduce vortex energy in all frequency ranges.This installation location is determined based on a dimensionless parameter (Xaft/h).

    1.Introduction

    In the past decades,a lot of studies on trailing edge noise have been done in analytical,numerical,and experimental research.This is due to the importance of trailing edge noise in a wide range of applications, such as aircraft, submarines, wind turbines, fans, etc.Oerlemans et al.conducted noise reduction research on a standard blade, a blade with trailing edge serrations,and a blade with an optimized airfoil shape used in wind turbines.1Both modified blades show a significant trailing edge noise reduction at low frequencies but the serrated blade showed more noise reduction.On the average, overall noise reduction of 3.2 dB is observed for the serrated blade.

    The reason for the wide frequency range of the trailing edge noise of the turbulent boundary layer is the presence of vortices with different sizes in this area.When the vortices create in the turbulent boundary layer near the trailing edge,pressure fluctuations in different frequency spectra occur due to the different sizes of the vortices, which causes the production of acoustic waves in the wide frequency spectrum.The physical process of creating trailing edge noise in a turbulent boundary layer was explained by Roger and Moreau.2They derived a correction and further extensions of Amit’s trailing edge noise model.

    As mentioned, the transient pressure field in the turbulent boundary layer near the trailing edge is one of the main causes of airfoil noise.Studies in 2008 showed that a part of the aircraft noise during the landing process was due to the intensity of pressure fluctuations in the turbulent boundary layer near the trailing edge of the wings and flaps.These studies were done by Brouwer and Rienstra.3Knowing the relationship between the velocity field of the turbulent boundary layer and the transient pressure field,and especially the relationship between the strength of the vortices and the pressure gradient on a wall is very important for controlling the flow.Andreopoulos and Agui did this research in 1996.4They used four high-frequency response pressure transducers in their research.Tsuji et al.measured pressure fluctuations in the turbulent boundary layer with a new technique in 2007.5In this technique,mean and pressure fluctuation and stream-wise velocity near the wall have been measured simultaneously in the turbulent boundary layer.In 2012, Ghaemi et al.conducted a transient three-dimensional velocity field measurement by employing a laser graphic imaging method in order to obtain the pressure field in the turbulent boundary layer.6In this experiment, the pressure field near the wall was also measured by capacitive microphones.The comparison showed agreement between the two pressure signals obtained from the tomo-PIV and the microphones with a cross correlation coefficient of 0.6 while their PSDs overlap up to 3 kHz.In 2013,Ghaemi and Scarano repeated the experiment on a flat plate with a new laser imaging technique.7The accuracy of experiments based on velocity field measurement within the turbulent boundary layer is more than pressure field measurement.Its reason is the lack of robust pressure measurement techniques in this area.They investigated the positive and negative high amplitude pressure peaks in a turbulent boundary layer.This research showed the positive high amplitude pressure peaks to the shear layer structures formed by an upstream sweep of high speed flow opposing a downstream ejection event.The hairpin category of vortices causes the ejection of near wall formation.It initiates the sweep event in the outer layer.The negative high amplitude pressure peaks are linked to both the span-wise and quasi-stream-wise vortices.

    With the growth of the aviation industry in the early 1990 s and the development of wind turbines, research on the use of passive noise reduction methods increased.These methods include airfoil trailing edge serrations,8–11trailing edge brushes,12,13porous trailing edge,14–16airfoil shape optimization,17equipping the surface with a 3D surface treatment element18–26and installing a vortex generator.

    Gstrein et al.investigated the effect of surface treatment with typical finlet spacing (4 mm) on a flat plate in 2020.20In this research, the Reynolds number was about 9.9 × 105near the trailing edge and the angle of attack was 0°.They suggested the best installing position of this surface treatment for reducing overall pressure PSD near the trailing edge.In 2021,Gstrein et al.investigated the effect of biologically inspired finlets on the reduction of trailing edge noise of a NACA0012 airfoil at Reynolds number of 400000.21They observed noise reduction of about 6 dB in the frequency range from 1000 Hz to 4000 Hz.In 2022, Fiscalleti et al.investigated the effect of stream-wise fences with different spacing on far-field noise experimentally.22When increasing the spacing between the fences from 2 mm to 4 mm, they observed a lower noise reduction at high frequencies and a higher noise reduction at low frequencies.Afshari et al.investigated the effect of the novel three-dimensional surface treatments on trailing edge noise reduction in 2019.26They showed that their aeroacoustic performance is better than the standard two-dimensional surface treatments.

    In the present study,the effect of the surface treatment with S=6 mm(finlet spacing)on velocity fluctuations in turbulent boundary layer is investigated.The effect of noise manipulation with the surface treatment with S = 6 mm has received less attention in previous researches.This element is installed slightly upstream of the trailing edge of the airfoil NACA2412 section so that the velocity field downstream of the finlets can be analyzed.Velocity fluctuations in turbulent boundary layer are measured by a hot wire sensor.In previous researches, it has been shown that if the surface treatment element is installed at the end of the trailing edge, the noise at high frequencies can be reduced; however, it will be shown that it is possible to find a mounting location on the airfoil that reduces the energy of the vortices in all frequency ranges near the trailing edge.Obviously, this mounting location depends on the characteristics of the surface treatment element, especially S(finlet spacing) and free stream conditions.

    2.Experimental equipment

    The tests were performed in a closed-circuit wind tunnel with an opened test section (Fig.1).The test section is a cylinder whose diameter is about 500 mm.Data acquisition in the velocity field was done by a one-dimensional hot wire sensor.In this sensor, a tungsten wire with a diameter of 5 μm and a length of 1.25 mm is used.

    The used model in the experiments is the NACA2412 section(Fig.2).Its chord is 430 mm and its span is 150 mm.Both sides of the model are confined by flat surfaces made of Plexiglas for preventing vortex shading (Fig.3).

    3.Test conditions

    Fig.1 Photo of wind tunnel.

    Fig.2 Model of NACA2412 section installed in test section of wind tunnel.

    Fig.3 Top view of model and installed trip and their dimensions.

    Fig.4 Views and dimensions of used surface treatment element.

    In all the tests, the free stream velocity in the test section is 30 m/s.Based on the free stream velocity and the model’s chord, the Reynolds number near the trailing edge is about 850000.The experiments were performed at three angles of attack(0°,4°,and 8°).At the highest angle of attack,the tunnel blockage ratio is less than 2%.At all angles of attack,the section was tested both without the surface treatment element and with its presence, so that the results could be compared.The hot wire anemometry was calibrated before each test.The frequency of data acquisition by the hot wire anemometry was 20 kHz and the data acquisition duration was 10 s.Therefore,200,000 data were saved at each point.The cutting frequency in most experiments was about 17 kHz.The maximum free stream Mach number was about 0.09 based on the lowest temperature of the test conditions.Therefore, the flow in the test section was incompressible.The surface treatment element must be used in the incompressible flow regime.A trip was installed near the leading edge of the airfoil (0.05C–0.1C) to ensure turbulent flow creation.The height of the used trip is 1 mm and its length is 6 mm.The end of the trip is located at 0.1C.

    4.Surface treatment geometry and its installation location

    The views of the used surface treatment element in this experiment are shown in Fig.4.In Fig.4, the dimensions of the designed element are specified.The length of this element in the x direction is 106.07 mm, and the height of each finlet is 7.02 mm.The total width of the surface treatment element is 150 mm.The distance between two consecutive finlets (finlet spacing) is 6 mm (S = 6 mm).The thickness of each finlet is 1.5 mm.The bottom of the constructed element has a slight curvature, which is determined by the installation location of the surface treatment on the NACA2412 section.

    The installation location of the surface treatment on the NACA2412 section is shown as follows.The beginning of the surface treatment element is located at x = 249 mm(x/C = 0.579) and the end of it is located at x = 355.07 mm (x/C = 0.826).The reason for installing the surface treatment element in this location is that the downstream area of the finlets can be investigated by a hot wire sensor.

    5.Data acquisition points

    Fig.5 shows the top view of the data acquisition point’s location which is near the wall in the z-direction(about 0.05 mm).The distance between two consecutive points in width is 1 mm.Points 0, 7, 14, 15, 16, 17, 28, and 34 are along the middle of two consecutive finlets, and Points 3,4,10,11,20,21,37,and 38 are located along behind the finlet approximately.Meanwhile, from 3 mm after the surface treatment element, exactly from the middle of two consecutive finlets,data were acquired in each 0.5 mm in the z-direction (perpendicular to the flow).

    6.Velocity field investigation

    At first, it should be noted that all experiments were repeated three times at different ambient temperatures and the obtained results were identical approximately.

    Fig.5 Top view of data acquisition points downstream of surface treatment element.

    Meanwhile,at Point 0(at x=37 cm or x/C=0.86),without the presence of the surface treatment element, the boundary layer thickness was measured at α = 0° and its value was about 18 mm.Considering that the height of the surface treatment element is 7.02 mm,after installing on its location on the NACA2412 section, it stayed on turbulent flow in all the angles of attack completely (as mentioned in the previous sections, a trip is used about 0.05C-0.1C for turbulent flow creation).Meanwhile, the Reynolds number is about 730300 at this Point 0.

    The velocity profiles were compared at three different points,with and without the presence of the surface treatment element.Meanwhile,these were compared to the results of the experiment that was done by Afshari et al.27They investigated the effect of two surface treatment elements in the turbulent flow installed on the flat plate at α=0°.In one of the surface treatment elements used by them,the distance between the two consecutive finlets is 2 mm(S=2 mm),and in the other 8 mm(S = 8 mm).Velocity profile diagrams have been dimensionless to allow comparison between the results of the presented study and Afshari et al.In these diagrams, obtained mean velocity profiles are validated by the 1/7 power law turbulent boundary layer profile in the absence of the surface treatment element.The small difference between the mean velocity and 1/7 power law profiles shows that the turbulent boundary layer has occurred at Xaft/h=2.5,Xaft/h=6.8,and Xaft/h=8.6.In these diagrams,Xaftis the distance from the end of the surface treatment element to the data acquisition points in the xdirection (stream-wise) and h is the height of the finlets.According to Fig.6, the smallest S creates the lowest momentum after the surface treatment element near the wall.In this case, there is a possibility of flow separation.It was reported by Afshari et al.27The momentum reduction by S = 6 mm is between S = 2 mm and S = 8 mm.It shows that the aerodynamic behavior of the surface treatment elements depends on their characteristics, especially S (finlet spacing).However,Afshari et al.showed that by changing the finlet spacing, the aeroacoustic behavior of the surface treatment element also changes.

    In this research, the energy of vortices in different frequencies at data acquisition points has also been investigated.For this purpose, the velocity Power Spectral Density (PSD)parameter has been measured at different points.To obtain the velocity PSD,first,velocity fluctuations relative to the time that are acquired by the hot wire sensor must transmit to the frequency domain.It is done by the Fourier transform.

    where U(ω,T) is the Fourier transform of u′(t).Then, velocity PSD can be calculated by28

    In the above relation,U*(ω,T)is a complex conjugation of U(ω,T) and E is the expected value that is a discrete random variable.To compare the energy of vortices in different frequencies, the following equation has been used in both cases,with and without the surface treatment element.In other words, ΔL?(f-C)is the logarithmic velocity PSD difference between the cases with and without the surface treatment element at one point.

    To amplify the above value, it is multiplied by 10.The above relation is a logarithmic velocity PSD difference.In this research, this equation is used to compare the energy of vortices at different frequencies in turbulent flow when the surface treatment element presents to that when it absents.

    Fig.7 shows the energy of vortices at different frequencies with the presence of surface treatment relative to that without it, at Points a, b, c, d, e, f, and g which are located at Xaft/h=0.43 and three angles of attack(0°,4°,and 8°).This figure indicates that the energy of high-frequency small vortices reduces by the surface treatment.But at α=0°,4°,the surface treatment increases the energy of vortices at low and intermediate frequencies,exactly after it(between two consecutive finlets).The reason for this phenomenon is presented in the next section.At these angles of attack,when points are close to the wall of the finlet, the energy of low-frequency large vortices increases strongly because the turbulent boundary layer of the finlet interacts downstream of its flow.At α = 8°, lowfrequency large vortex energy reduces.The energy of vortices in all frequency ranges is reduced at points located behind the finlet at all angles of attack because this area is a lowenergy region.A phenomenon happens similar to cavitation behind the finlet.

    Fig.8 shows the energy of vortices at different frequencies with the presence of surface treatment relative to that without it, at Points 0, 1, 2, 3, 4, 5, and 6 which are located at x/C = 0.86 or Xaft/h = 2.5.This figure is drawn at three angles of attack (0°, 4°, and 8°).This figure indicates that the surface treatment with S = 6 mm has reduced the energy of highfrequency small vortices.This reduces the intensity of highfrequency pressure fluctuations near the wall and also decreases the Sound Pressure Level (SPL) in far-field at high frequencies.Of course,at these points,the energy of large vortices with low frequency is almost under control and has not changed much.Meanwhile, the surface treatment with S = 6 mm reduces the vortex energy at intermediate frequencies approximately.The cause of these phenomena is presented in the next section.Generally, at Xaft/h = 2.5, the energy of vortices is decreased at high and intermediate frequencies by the surface treatment element with S = 6 mm and is almost controlled at low frequencies.

    Fig.6 Velocity profiles at α = 0° and different Xaft/h.

    Fig.7 10 × ΔL?(f-C) at Points a, b, c, d, e, f and g at Xaft/h = 0.43.

    Fig.8 10 × ΔL?(f-C) at Points 0, 1, 2, 3, 4, 5 and 6 at Xaft/h = 2.5.

    Fig.9 (Xaft/h = 6.8), Fig.10 (Xaft/h = 8.6), and Fig.11(Xaft/h=10.11)show that the energy of the vortices is increased at low and intermediate frequencies by the surface treatment element relative to points that are located at Xaft/h=2.5.The cause of this phenomenon is presented in the next section.Of course,in all of these diagrams, the reduction of high-frequency vortex energy can be seen at all angles of attack.Of course, there are exceptions at some points.At Xaft/h = 6.8 and α = 0°, Points 11 and 12 (these points are along behind the finlet approximately),at Xaft/h=6.8 and α=0°,Points 20 and 21(these points are along behind the finlet approximately),at Xaft/h=10.11 and α=0°,Points 33 and 34(these points are along the middle of two consecutive finlets),the energy of vortices does not change and remains constant at low and intermediate frequencies near the wall compared to the state that there is no surface treatment element approximately.

    Fig.9 10 × ΔL?(f-C) at Points 7, 8, 9, 10, 11, 12 and 13 at Xaft/h = 6.8.

    Fig.12 shows the energy of vortices at different frequencies near the wall, with the presence of the surface treatment element, relative to its absence, at points that are along the middle of two consecutive finlets.These graphs are also drawn at three angles of attack (0°, 4°, and 8°).Figs.13 and 14 show the energy of vortices at different frequencies near the wall, with the presence of the surface treatment element, relative to its absence, at points that are along behind the finlet approximately.In Figs.12–14, it is observed that the lowest energy of vortices in all frequency ranges occurs at Xaft/h = 2.5 by the surface treatment element with S = 6 mm.

    Fig.10 10×ΔL?(f-C)at Points 17,18,19,20,21,22,23,24,25,26 and 27 at Xaft/h = 8.6.

    Fig.11 10×ΔL?(f-C)at Points 29,30,31,32,33,34,35,36,37,38 and 39 at Xaft/h = 10.11.

    The following equation has been used to better determine which points have the least vortex energy in all frequency ranges with the presence of the surface treatment element with S=6 mm.ΔL?fin(P2-P1)compares the energy of the vortices in all frequency spectra between two points, while the surface treatment element is the installation.To amplify the bellow value, it is multiplied by 10 in figures.

    Fig.12 10 × ΔL?(f-C) at points along the middle of two consecutive finlets.

    Fig.15 compares the energy of vortices, at different frequencies, of points located at Xaft/h = 2.5 to that of points at Xaft/h = 0.43.This figure shows that the energy of vortices in all frequency spectra is reduced by the surface treatment with S = 6 mm at Xaft/h = 2.5 relative to Xaft/h = 0.43 approximately.Of course, some exceptions are seen.Since the cavitation occurs behind the finlet approximately,this area is a low-energy region.This causes the vortex energy of some points at Xaft/h = 2.5 to be more than that of the points located behind the finlet.This is shown in Fig.15, with black dashed and green plots.Fig.16 also compares the energy of vortices, at different frequencies, of points located at Xaft/h = 6.8 to that of points at Xaft/h = 2.5.These comparisons are done at three angles of attack(0°,4°,and 8°).In this figure,it is seen clearly that the surface treatment with S = 6 mm increases the energy of vortices in all frequency ranges at Xaft/h = 6.8 relative to points located at Xaft/h = 2.5.Therefore, the least vortex energy in all frequency spectra happens close to points located at Xaft/h = 2.5 clearly.

    Fig.13 10 × ΔL?(f-C) at some points along behind finlets at Points 3, 10, 20, 37.

    Based on the above, it can be concluded that in order to reduce velocity and pressure PSD in all frequency spectra near the wall, the surface treatment with S = 6 mm should be moved slightly towards the leading edge relative to the end of the trailing edge.This displacement is about a dimensionless distance of Xaft/h = 2.5.Gstrein et al.also reported the accuracy of this information for the surface treatment with S = 4 mm.20Their research was done on a flat plate at α=0°.In this research,the Reynolds number near the trailing edge was about 9.9 × 105.Some results of this research are illustrated in Fig.17.In Fig.17, L, LF, and PFare flat plate length, finlet length, and finlet position, respectively.Gstrein et al.showed that the least overall pressure PSD (in all frequency ranges) occurs at PF= 0.835L by the surface treatment with S = 4 mm.Overall pressure PSD was measured by a microphone installed at x/L = 0.996.

    Fig.14 10 × ΔL?(f-C) at some points along behind finlets at Points 4, 11, 21, 38.

    Fig.18 shows the velocity profiles at Point a,with the presence of the surface treatment element.As shown in Fig.18,flow separation does not occur after the surface treatment element with S = 6 mm at all angles of attack (α = 0°, 4°, and 8°).However,according to Afshari’s report27,flow separation occurs after the surface treatment element with S = 2 mm at α = 0°, on the flat plate.

    Fig.15 10×ΔL?fin(P2-P1)between Points 0,1,2,3,4,5 and 6 at Xaft/h = 2.5 and Points a, b, c, d, e, f and g at Xaft/h = 0.43.

    Data acquisition was taken at Point a in the z-direction,with and without the surface treatment element.The results are shown in Fig.19.According to the presented information in Fig.19,the energy of the vortices is clearly reduced in all frequency spectra by the surface treatment element with S = 6 mm at 0.57 ≤z/h ≤0.855 and at α = 0°, 4°.This decrease happens at low frequencies near the wall at α=8°.In the same way, it is clear that the surface treatment element causes the increase in the vortex energy in all frequency spectra at z/h ≥1 strongly.This increase is greater at intermediate frequencies.Its reason is flow injection from between finlets into the turbulent flow that passes above them.The growth of the boundary layer is caused to reduce in the mass flow rate between two consecutive finlets and flow injection into the stream that passes through the top of the finlets.

    Fig.16 10×ΔL?fin(P2-P1) between Points 7,8, 9,10,11,12 and 13 at Xaft/h=6.8 and Points 0,1,2,3,4,5 and 6 at Xaft/h=2.5.

    Finally, it should be noted that the behavior of finlets changes with geometric properties, especially the S parameter.Afshari et al.27and Gstrein et al.20showed this in their researches.Afshari et al.investigated the effect of the surface treatment with various S (finlet spacing) on turbulence intensity profiles at different Xaft/h.According to Afshari’s research,turbulence intensity near the wall is reduced by the surface treatment that has S = 8 mm with Xaft/h increasing.But in the case of the surface treatment with S=2 mm,the opposite is true.However, in this article, it is proved that the surface treatment with S = 6 mm reduces turbulence intensity near the wall until Xaft/h = 2.5 and then increases it after this location.

    Fig.17 Surface pressure fluctuation PSD for conventional finlets with h=12 mm,S=4 mm and LF=65 mm applied at various positions at ReL = 9.9 × 105.20

    Fig.18 Velocity profiles at 3 mm after surface treatment element with S = 6 mm, exactly between two consecutive finlets.

    7.Analysis of some phenomena

    In this section, the reasons for some of the phenomena are explained.

    Fig.19 10 × ΔL?(f-C) at points in z-direction.

    On the walls of two consecutive finlets, the growth of the boundary layer occurs.This causes the increase in the kinetic energy of the flow between the finlets at the end of them.It increases the energy of vortices at low and intermediate frequencies at Xaft/h = 0.43 near the wall.Fig.7 confirms the occurrence of this phenomenon.It has also been shown by Afshari et al.27.Their research has been done on a flat plate at α = 0°.Meanwhile, Bodling et al.have also mentioned this phenomenon.29

    Fig.20 Effect of surface treatment with various S (finlet spacing) on turbulence intensity profile at different Xaft/h.27

    After the finlets (after Xaft/h = 0.43), a decrease in kinetic energy occurs in the turbulent flow structure.This causes a decrease in the energy of vortices to be seen in all frequency ranges at Xaft/h = 2.5 near the wall.This is shown in Fig.15.Meanwhile, this phenomenon has been reported by Afshari et al.27and Bodling et al.29

    Due to the growth of the boundary layer between the finlets, the mass flow rate decreases in this region.This causes some flow to be injected into the turbulent flow passing through the top of the finlets.By injecting the flow,the energy of vortices increases in all frequency spectra strongly.This phenomenon can be well seen in Fig.19.

    Part of the turbulent flow above the finlets hits the wall of the airfoil after Xaft/h=2.5 probably.Vortex shedding on the wall increases the energy of vortices in all frequency ranges after Xaft/h=2.5 near the wall.Figs.9-14 and Fig.16 show it.

    Right behind the finlets,low-energy regions are created similar to the cavitation phenomenon.This is well seen in Figs.7 and 15 (black dashed and green plots).

    8.Conclusions

    Noise reduction of different airfoils is important.Several passive noise reduction methods are used on airfoils.One of these methods is the use of surface treatment.In this research, the effect of the surface treatment element with S = 6 mm on the vortex structure at different frequencies in the turbulent flow that was created on the NACA2412 section has been investigated.Based on the performed experiments, the following results can be presented:

    (1) Flow separation does not occur exactly after the surface treatment element with S = 6 mm up to α = 8°.However, Afshari et al.reported that flow separation happens after the surface treatment with S = 2 mm on a flat plate at α = 0°.27It is shown in Fig.18.

    (2) A noticeable momentum decrease occurs after the surface treatment element with S = 6 mm, near the wall at Xaft/h=2.5,Xaft/h=6.8,and Xaft/h=8.6.According to Afshari’s research,27the momentum decreasing with the presence of a surface treatment element with S = 6 mm is less than S = 2 mm and more than S=8 mm.According to items 1 and 2 as well as Figs.6 and 20, the performance of the surface treatment element on the flow behind it depends on its characteristics,especially S (finlet spacing).

    (3) The energy of the vortices at the intermediate and high frequencies, with the presence of the surface treatment element that has S = 6 mm, decreases more at Xaft/h = 2.5 compared to Xaft/h = 0.43 and remains almost constant at low frequencies.The reason may be the occurrence of a fully developed flow between the finlets, which causes the flow velocity to increase at Xaft/h = 0.43 near the wall.

    (4) The energy of vortices in all frequency spectra is increased after Xaft/h=2.5 by the surface treatment element with S=6 mm.It is shown in Fig.16.The reason for that is the injection of flow and vortices on the top of the finlets and then a part of them shedding after Xaft/h = 2.5, at all angles of attack.The increase in the energy of the vortices in all frequency ranges above the finlets is well seen in Fig.19.This phenomenon is explained in Section 7.

    (5) The main effective factor of SPL in the far-field is the intensity of pressure fluctuations that occur near the wall at different frequencies.In order to reduce velocity and pressure PSD in all frequency spectra near the wall, the surface treatment with S = 6 mm should be moved slightly towards the leading edge of the NACA2412 airfoil relative to the end of the trailing edge.This displacement is about a dimensionless distance of Xaft/h = 2.5.

    (6) Just after the surface treatment element with S=6 mm and at 0.57 ≤z/h ≤0.855, the energy of vortices in all frequency spectra is clearly reduced at α = 0°, 4°.It happens in the middle of two consecutive finlets.It is seen in Fig.19.

    (7) Right after the surface treatment element with S = 6 mm and near the wall, the energy of vortices is reduced at low frequencies and α = 8°.It is shown in Fig.19.

    (8) At z/h ≥1,the energy of vortices in all frequency ranges is strongly increased, at all angles of attack.This increase is further at intermediate frequencies.The reason may be the injection of flow from between the finlets to their top.

    Declaration of Competing Interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    成人国产麻豆网| 99久久久亚洲精品蜜臀av| 看片在线看免费视频| 91狼人影院| 国产精品乱码一区二三区的特点| 美女被艹到高潮喷水动态| 国产高清有码在线观看视频| 激情 狠狠 欧美| 欧美性感艳星| 欧美日韩乱码在线| 日本爱情动作片www.在线观看 | 中文字幕精品亚洲无线码一区| 亚洲成a人片在线一区二区| 老司机午夜福利在线观看视频| 日本在线视频免费播放| 国产在线精品亚洲第一网站| 成人精品一区二区免费| 日韩成人伦理影院| 久久久久国产网址| 全区人妻精品视频| 亚洲精品日韩在线中文字幕 | 最新中文字幕久久久久| 亚洲自偷自拍三级| 我要搜黄色片| 给我免费播放毛片高清在线观看| 国产人妻一区二区三区在| 免费看av在线观看网站| 真人做人爱边吃奶动态| 欧美+亚洲+日韩+国产| АⅤ资源中文在线天堂| 三级毛片av免费| 国产探花在线观看一区二区| 成年女人看的毛片在线观看| 长腿黑丝高跟| 秋霞在线观看毛片| 一个人看视频在线观看www免费| 又粗又爽又猛毛片免费看| 久久人人爽人人爽人人片va| 午夜激情福利司机影院| 中文资源天堂在线| 亚洲成人中文字幕在线播放| 久久韩国三级中文字幕| 国产大屁股一区二区在线视频| 人妻制服诱惑在线中文字幕| 国产高清激情床上av| 亚洲欧美日韩卡通动漫| 久久久久精品国产欧美久久久| 国产白丝娇喘喷水9色精品| 国产一区二区在线观看日韩| 毛片一级片免费看久久久久| 真实男女啪啪啪动态图| 欧美三级亚洲精品| 国产一区亚洲一区在线观看| 亚洲最大成人手机在线| 国产精品电影一区二区三区| av天堂在线播放| 激情 狠狠 欧美| 久久精品国产亚洲网站| 少妇熟女aⅴ在线视频| 少妇人妻一区二区三区视频| 久久6这里有精品| 女的被弄到高潮叫床怎么办| 在线观看66精品国产| 久久精品国产亚洲av香蕉五月| 亚洲五月天丁香| 高清毛片免费看| 波野结衣二区三区在线| 日韩欧美三级三区| 国产黄色视频一区二区在线观看 | 又爽又黄a免费视频| 变态另类成人亚洲欧美熟女| 九色成人免费人妻av| 别揉我奶头~嗯~啊~动态视频| 啦啦啦韩国在线观看视频| 悠悠久久av| 男女边吃奶边做爰视频| 一级av片app| 国内精品一区二区在线观看| 免费电影在线观看免费观看| 国产黄色视频一区二区在线观看 | 国产精品一区二区免费欧美| 91午夜精品亚洲一区二区三区| 久久久久久国产a免费观看| 乱系列少妇在线播放| 真人做人爱边吃奶动态| 3wmmmm亚洲av在线观看| 男插女下体视频免费在线播放| 午夜久久久久精精品| 日本黄色视频三级网站网址| a级毛片免费高清观看在线播放| 一级黄片播放器| 亚洲性久久影院| 成人三级黄色视频| 最后的刺客免费高清国语| 丝袜喷水一区| 婷婷精品国产亚洲av在线| 国产精品女同一区二区软件| 久久精品国产亚洲av天美| 在线a可以看的网站| 丰满的人妻完整版| 久久九九热精品免费| 老师上课跳d突然被开到最大视频| 欧美三级亚洲精品| 欧美性猛交╳xxx乱大交人| 男人狂女人下面高潮的视频| 九九在线视频观看精品| 久久精品综合一区二区三区| 午夜精品国产一区二区电影 | 1000部很黄的大片| 免费大片18禁| 蜜桃久久精品国产亚洲av| av天堂中文字幕网| 久久久精品94久久精品| 女的被弄到高潮叫床怎么办| 不卡一级毛片| 免费大片18禁| 午夜久久久久精精品| 麻豆久久精品国产亚洲av| 丰满乱子伦码专区| 老司机午夜福利在线观看视频| 三级经典国产精品| 日本精品一区二区三区蜜桃| 免费在线观看成人毛片| 精品欧美国产一区二区三| 99热精品在线国产| 欧美bdsm另类| 久久国内精品自在自线图片| 麻豆一二三区av精品| 国产精品福利在线免费观看| 18禁在线无遮挡免费观看视频 | a级毛片a级免费在线| 欧美3d第一页| 欧美区成人在线视频| 十八禁网站免费在线| 久99久视频精品免费| 尾随美女入室| 国内揄拍国产精品人妻在线| 简卡轻食公司| 精品福利观看| 九九热线精品视视频播放| 干丝袜人妻中文字幕| 色综合色国产| 女生性感内裤真人,穿戴方法视频| 在线免费十八禁| 国产在视频线在精品| 免费看光身美女| 国产三级中文精品| 亚洲综合色惰| 欧美一区二区国产精品久久精品| 国产黄色小视频在线观看| 亚洲国产色片| 老司机福利观看| 啦啦啦观看免费观看视频高清| 欧美激情国产日韩精品一区| 亚洲精品亚洲一区二区| 哪里可以看免费的av片| 久久午夜福利片| 久久精品国产99精品国产亚洲性色| 在线国产一区二区在线| 亚洲成人久久性| 亚洲av.av天堂| 久久6这里有精品| 国产单亲对白刺激| 国产精品人妻久久久久久| 成人av一区二区三区在线看| 日韩欧美一区二区三区在线观看| 狠狠狠狠99中文字幕| 最新在线观看一区二区三区| 成人美女网站在线观看视频| 国产高清不卡午夜福利| 久久精品国产自在天天线| 亚洲av一区综合| 美女高潮的动态| 内射极品少妇av片p| av中文乱码字幕在线| 日韩欧美免费精品| 成年女人毛片免费观看观看9| 国产精品99久久久久久久久| 大型黄色视频在线免费观看| 久久人人爽人人爽人人片va| 日本一二三区视频观看| 精品久久久久久久久久久久久| 黄色日韩在线| 亚洲熟妇中文字幕五十中出| 最近在线观看免费完整版| 国产aⅴ精品一区二区三区波| 尾随美女入室| 国内少妇人妻偷人精品xxx网站| 精品乱码久久久久久99久播| 一本久久中文字幕| 成人综合一区亚洲| 欧美+亚洲+日韩+国产| 日本爱情动作片www.在线观看 | 久久99热这里只有精品18| 午夜老司机福利剧场| 国产欧美日韩精品一区二区| 亚洲乱码一区二区免费版| 搡老妇女老女人老熟妇| 午夜精品国产一区二区电影 | 真人做人爱边吃奶动态| 亚洲精品日韩av片在线观看| 国产精品久久视频播放| 中文在线观看免费www的网站| av.在线天堂| 国产精品一区二区三区四区久久| 1024手机看黄色片| 日日摸夜夜添夜夜添小说| 亚洲精品日韩在线中文字幕 | 麻豆精品久久久久久蜜桃| 黄色日韩在线| 嫩草影院新地址| 男女那种视频在线观看| eeuss影院久久| 国产淫片久久久久久久久| 国产探花极品一区二区| 免费不卡的大黄色大毛片视频在线观看 | 天堂网av新在线| 欧美日韩在线观看h| 春色校园在线视频观看| 少妇高潮的动态图| 99九九线精品视频在线观看视频| 色av中文字幕| 人妻制服诱惑在线中文字幕| 久久精品夜夜夜夜夜久久蜜豆| 男女视频在线观看网站免费| 亚洲不卡免费看| 亚洲国产高清在线一区二区三| 国产国拍精品亚洲av在线观看| 99热这里只有是精品在线观看| 欧美绝顶高潮抽搐喷水| 成人国产麻豆网| 18禁裸乳无遮挡免费网站照片| 中文字幕久久专区| 久久久久久伊人网av| 99热这里只有是精品在线观看| 国产熟女欧美一区二区| 男女视频在线观看网站免费| 伊人久久精品亚洲午夜| 国产精品爽爽va在线观看网站| 直男gayav资源| 波多野结衣高清无吗| videossex国产| 中文资源天堂在线| a级毛片免费高清观看在线播放| 尤物成人国产欧美一区二区三区| 欧美日韩精品成人综合77777| 熟女人妻精品中文字幕| 身体一侧抽搐| 男插女下体视频免费在线播放| 两性午夜刺激爽爽歪歪视频在线观看| 欧美一区二区国产精品久久精品| 黄色欧美视频在线观看| 高清日韩中文字幕在线| 亚洲不卡免费看| av中文乱码字幕在线| 99久国产av精品国产电影| 国产精品国产高清国产av| 三级经典国产精品| 亚洲无线在线观看| 啦啦啦观看免费观看视频高清| 最新中文字幕久久久久| 欧美成人精品欧美一级黄| 日本色播在线视频| 久久久久久九九精品二区国产| 色哟哟·www| 亚洲精品乱码久久久v下载方式| 亚洲真实伦在线观看| 亚洲四区av| 色综合色国产| 岛国在线免费视频观看| 欧美绝顶高潮抽搐喷水| 淫秽高清视频在线观看| 看片在线看免费视频| 97超视频在线观看视频| 级片在线观看| 久久久成人免费电影| 久久精品国产99精品国产亚洲性色| ponron亚洲| 18禁裸乳无遮挡免费网站照片| 久久鲁丝午夜福利片| 久久久久久国产a免费观看| 最近最新中文字幕大全电影3| 久久欧美精品欧美久久欧美| 亚洲成人久久爱视频| 高清午夜精品一区二区三区 | 麻豆国产av国片精品| 亚洲经典国产精华液单| 丝袜美腿在线中文| 欧美激情久久久久久爽电影| 一个人免费在线观看电影| 天堂av国产一区二区熟女人妻| 在线观看一区二区三区| 免费看av在线观看网站| 国产av在哪里看| 欧美色视频一区免费| 国产高清不卡午夜福利| 久久天躁狠狠躁夜夜2o2o| 成人永久免费在线观看视频| 我要看日韩黄色一级片| 亚洲成av人片在线播放无| 九色成人免费人妻av| 久久久a久久爽久久v久久| 成人国产麻豆网| 老司机午夜福利在线观看视频| 欧美又色又爽又黄视频| 麻豆乱淫一区二区| 国产三级中文精品| 中文字幕精品亚洲无线码一区| 精品一区二区三区人妻视频| 亚洲欧美日韩高清专用| 插逼视频在线观看| 97热精品久久久久久| 人妻制服诱惑在线中文字幕| 国产黄色小视频在线观看| 久久久久久久久中文| 九色成人免费人妻av| 亚洲在线自拍视频| 亚洲欧美成人综合另类久久久 | 在线免费观看不下载黄p国产| 欧美三级亚洲精品| 一本一本综合久久| 两个人的视频大全免费| 精品一区二区三区视频在线| 国产高潮美女av| 俺也久久电影网| av福利片在线观看| 一个人观看的视频www高清免费观看| 久久精品国产亚洲av涩爱 | 桃色一区二区三区在线观看| 婷婷精品国产亚洲av在线| 国产精品美女特级片免费视频播放器| 22中文网久久字幕| 亚洲成人久久爱视频| 12—13女人毛片做爰片一| 国内精品宾馆在线| 免费在线观看影片大全网站| 国产真实乱freesex| 久久精品91蜜桃| 国产成人影院久久av| 三级男女做爰猛烈吃奶摸视频| 99国产精品一区二区蜜桃av| 国语自产精品视频在线第100页| 91在线精品国自产拍蜜月| 欧美+日韩+精品| 欧美性感艳星| 成人永久免费在线观看视频| 亚洲七黄色美女视频| 欧美人与善性xxx| 深爱激情五月婷婷| 一个人看视频在线观看www免费| 偷拍熟女少妇极品色| 尾随美女入室| 一个人免费在线观看电影| 成人永久免费在线观看视频| а√天堂www在线а√下载| 久久久久久久久久成人| 亚洲七黄色美女视频| 熟妇人妻久久中文字幕3abv| 狠狠狠狠99中文字幕| 亚洲最大成人中文| 少妇猛男粗大的猛烈进出视频 | 国产三级在线视频| 免费看日本二区| 日韩欧美精品v在线| 啦啦啦啦在线视频资源| 日韩中字成人| 亚洲欧美精品综合久久99| 五月伊人婷婷丁香| 亚洲av不卡在线观看| 日本与韩国留学比较| a级毛色黄片| 婷婷精品国产亚洲av在线| 亚洲人与动物交配视频| av在线观看视频网站免费| 日本撒尿小便嘘嘘汇集6| 99国产极品粉嫩在线观看| 麻豆精品久久久久久蜜桃| 人人妻人人澡欧美一区二区| 欧美成人精品欧美一级黄| 黄片wwwwww| 一级av片app| 亚洲成人中文字幕在线播放| 亚洲美女搞黄在线观看 | 99热精品在线国产| a级毛片a级免费在线| 18+在线观看网站| 激情 狠狠 欧美| 欧美日本视频| 日本爱情动作片www.在线观看 | 欧美色视频一区免费| 亚洲美女搞黄在线观看 | 亚洲天堂国产精品一区在线| 国产高清视频在线观看网站| 精品久久久久久久人妻蜜臀av| 丰满乱子伦码专区| 欧美激情国产日韩精品一区| 91久久精品国产一区二区成人| 一个人看的www免费观看视频| 嫩草影院入口| 日本黄色视频三级网站网址| 国产人妻一区二区三区在| 又粗又爽又猛毛片免费看| 久久久久性生活片| 久久国产乱子免费精品| 久久精品国产亚洲av天美| 久久久久久久久大av| 日本一本二区三区精品| 久久精品综合一区二区三区| 欧美人与善性xxx| 岛国在线免费视频观看| 精品一区二区免费观看| 精品人妻一区二区三区麻豆 | 欧洲精品卡2卡3卡4卡5卡区| 亚洲av.av天堂| 精品久久久噜噜| 九色成人免费人妻av| 中文字幕av在线有码专区| 噜噜噜噜噜久久久久久91| av在线亚洲专区| а√天堂www在线а√下载| 看十八女毛片水多多多| 午夜福利视频1000在线观看| 日韩av不卡免费在线播放| 日本撒尿小便嘘嘘汇集6| 亚洲成人中文字幕在线播放| 级片在线观看| 色哟哟·www| 亚洲国产日韩欧美精品在线观看| 女人十人毛片免费观看3o分钟| 黄色日韩在线| 99在线视频只有这里精品首页| 午夜福利高清视频| 一级av片app| 草草在线视频免费看| 男女那种视频在线观看| 亚洲无线在线观看| 久久久久国产精品人妻aⅴ院| 美女内射精品一级片tv| 精品一区二区三区人妻视频| a级毛色黄片| 一a级毛片在线观看| 午夜精品国产一区二区电影 | 淫秽高清视频在线观看| 人妻制服诱惑在线中文字幕| 最近手机中文字幕大全| 亚洲激情五月婷婷啪啪| 亚洲av不卡在线观看| 高清日韩中文字幕在线| 精品久久久久久久人妻蜜臀av| 免费看美女性在线毛片视频| 亚洲av二区三区四区| 国产精华一区二区三区| 国产黄色小视频在线观看| 99久久九九国产精品国产免费| 性色avwww在线观看| 身体一侧抽搐| 国产亚洲欧美98| 五月伊人婷婷丁香| 又黄又爽又刺激的免费视频.| 亚洲精品粉嫩美女一区| 熟妇人妻久久中文字幕3abv| 日韩欧美精品免费久久| 一进一出抽搐gif免费好疼| 麻豆av噜噜一区二区三区| 久久精品久久久久久噜噜老黄 | 内地一区二区视频在线| 亚洲自偷自拍三级| 亚洲无线观看免费| 亚洲高清免费不卡视频| 欧美一区二区精品小视频在线| 色视频www国产| 久久精品国产亚洲av香蕉五月| 午夜精品一区二区三区免费看| 激情 狠狠 欧美| 伦理电影大哥的女人| 禁无遮挡网站| 自拍偷自拍亚洲精品老妇| 高清毛片免费看| 长腿黑丝高跟| 两性午夜刺激爽爽歪歪视频在线观看| 一本久久中文字幕| 人人妻,人人澡人人爽秒播| 欧美极品一区二区三区四区| 国产精品一区二区三区四区久久| 色5月婷婷丁香| 亚洲婷婷狠狠爱综合网| 天天躁夜夜躁狠狠久久av| 国产v大片淫在线免费观看| 免费电影在线观看免费观看| 欧美高清成人免费视频www| 亚洲七黄色美女视频| 日本黄色视频三级网站网址| 国产视频内射| 国产成人91sexporn| 亚洲av中文字字幕乱码综合| 欧美日韩乱码在线| 久久精品91蜜桃| 99热这里只有是精品在线观看| 男女做爰动态图高潮gif福利片| 最近最新中文字幕大全电影3| 综合色av麻豆| 三级男女做爰猛烈吃奶摸视频| 可以在线观看的亚洲视频| 搡老岳熟女国产| 亚洲av成人精品一区久久| 18禁在线播放成人免费| 伦理电影大哥的女人| 久久久国产成人免费| 午夜影院日韩av| 亚洲天堂国产精品一区在线| 中国国产av一级| 成人亚洲欧美一区二区av| 少妇人妻精品综合一区二区 | 国产一区二区在线av高清观看| 免费看光身美女| 午夜精品在线福利| 亚洲精品日韩av片在线观看| 精品久久久噜噜| 天美传媒精品一区二区| 亚洲中文字幕一区二区三区有码在线看| aaaaa片日本免费| 丰满的人妻完整版| 97超碰精品成人国产| 一区二区三区免费毛片| 欧美+亚洲+日韩+国产| 欧洲精品卡2卡3卡4卡5卡区| 欧美一区二区国产精品久久精品| 国语自产精品视频在线第100页| 成人av在线播放网站| 97人妻精品一区二区三区麻豆| 亚洲第一电影网av| 亚洲精品国产成人久久av| 波野结衣二区三区在线| 秋霞在线观看毛片| 国产成人91sexporn| 免费av观看视频| 日韩高清综合在线| 2021天堂中文幕一二区在线观| av黄色大香蕉| 成人漫画全彩无遮挡| 熟女电影av网| 又爽又黄无遮挡网站| 亚洲18禁久久av| 日韩,欧美,国产一区二区三区 | 丰满乱子伦码专区| 在线播放无遮挡| 亚洲av一区综合| 国内精品一区二区在线观看| 一个人观看的视频www高清免费观看| 亚洲欧美日韩东京热| 精品免费久久久久久久清纯| 久久久久久久久中文| 色综合站精品国产| 国产美女午夜福利| 中文资源天堂在线| 在线观看美女被高潮喷水网站| 国产精品99久久久久久久久| 亚洲av第一区精品v没综合| 色综合色国产| videossex国产| 免费观看在线日韩| 国产精品一区www在线观看| 国产av在哪里看| 精品午夜福利视频在线观看一区| 最新中文字幕久久久久| 一本精品99久久精品77| 亚洲av成人av| 欧美高清成人免费视频www| 亚洲av第一区精品v没综合| 波野结衣二区三区在线| 亚洲四区av| 亚洲激情五月婷婷啪啪| 性欧美人与动物交配| 男人舔女人下体高潮全视频| 亚洲国产精品国产精品| 国产精品,欧美在线| a级一级毛片免费在线观看| 中国美女看黄片| 欧美日韩在线观看h| 大又大粗又爽又黄少妇毛片口| 久久韩国三级中文字幕| 国产av在哪里看| 淫妇啪啪啪对白视频| а√天堂www在线а√下载| 在线观看av片永久免费下载| 国产在视频线在精品| 国产黄色小视频在线观看| 老熟妇仑乱视频hdxx| 在线播放国产精品三级| 我要搜黄色片| 老熟妇仑乱视频hdxx| 国产乱人偷精品视频| 免费观看人在逋| 国内少妇人妻偷人精品xxx网站| 尤物成人国产欧美一区二区三区| 在线观看美女被高潮喷水网站| 成人av一区二区三区在线看| 直男gayav资源| 免费观看人在逋| ponron亚洲| 一本精品99久久精品77| 狠狠狠狠99中文字幕| 国产三级中文精品| 国产在线精品亚洲第一网站| 亚洲国产精品成人综合色| av天堂在线播放| 人妻丰满熟妇av一区二区三区| 91久久精品国产一区二区成人| 日产精品乱码卡一卡2卡三| 黄色欧美视频在线观看| 精品久久久久久久末码| 亚洲乱码一区二区免费版| 成人欧美大片| 亚洲精品国产av成人精品 | 真实男女啪啪啪动态图|