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

    Effects of geometry variations on tandem airfoil interaction noise

    2020-08-08 02:46:28ZHANGChiSONGWenbin
    空氣動力學學報 2020年3期
    關鍵詞:接收點遠場聲壓級

    ZHANG Chi, SONG Wenbin

    (School of Aeronautics and Astronautics, Shanghai Jiao Tong University, Shanghai 200240, China)

    Abstract:Interaction noise of a tandem airfoil configuration is studied using high-order large eddy simulations (LES) methods on a multi-block structured grid. The sound could be generated from aerodynamic interactions between the two airfoils positioned one after the other. The flow field is first calculated, and spectral properties of pressure fluctuations at a few given observation points are analyzed. The far-field noise obtained using FW-H integration is presented for the base configuration. The interaction from the wake of the airfoil in the front could interact with the leading edge of the trailing airfoil for specific configuration settings. The geometric features of both airfoils are also varied simultaneously to study the combined effects on the noise. The impact of geometric characteristics of trailing-edge and leading-edge of the front and rear airfoil, respectively, is the focus of the research. A typical surrogate-based approach is deployed in this study to reduce computational costs. The results show that, by careful tuning of the geometry and configuration, aeroacoustic characteristics of the tandem airfoil configuration could be modified to achieve lower acoustic levels.

    Key words:aeroacoustics; high order scheme; tandem airfoil; interaction noise; FW-H method; response surface

    0 Introduction

    The noise standards for commercial aircraft operating in civil airspace surrounding airports are becoming increasingly stringent. The proposal from the Advisory Council for Aviation Research in Europe (ACARE) aims to achieve a 10 dB reduction by 2020[1]. Compared with the present Chapter 4 standard, the noise standard introduced by the International Civil Aviation Organization (ICAO) in February 2010 requires civil aircraft airworthiness noise emission to be reduced by 7 EPNdB[1]. The goal of N+3 generation NASA research aircraft wished to gain a 52 dB reduction in combined noise levels by 2025 compared to that of aircraft entering service in 2000[1]. Huge challenge is faced by engineers considering the gap between the state of the art and ambitious objectives set by the relevant organizations.

    Aircraft noise comes primarily from two sources, including airframe noise and engine noise. Both need to be reduced to a comparable level to achieve the overall objective of noise reduction. Aerodynamically generated noise is a major mechanism for aircraft noise. Two primary sources for aircraft noise are landing gear and high-lift devices. A number of different noise generation mechanisms have been identified in the literature. The interaction noise between two airfoils placed in tandem configuration has gained interest because it represents one of many typical configurations for aerodynamic noise generation. For example, high bypass ratio turbofan engine and contra-rotating open rotor (CROR) engine may lead to further reductions in fuel consumption, but will likely lead to higher noise levels compared to the current turbofan engines. The levels of noise generated from the interactions between rotor/stator in the turbofan engine and front/aft rotor in CROR can be higher than that of the self-noise of the rotor. In the take-off and landing configurations, the wake and vortex of the main wing would also interact with the horizontal tails. In these cases, interaction noise becomes one of the dominant sources. Such a configuration, though slightly different type, can also be found in typical landing gear configurations, in which, the wake shed from components in the upstream interacts with components in the downstream, generating a major source of noise for the take-off and landing configurations. The interactions of the wake from the front airfoil with the second airfoil are the primary research focus of the current paper.

    Changes in the geometry of leading-edge and trailing edge ofairfoils, rotors or stators and their relative position will alter the surrounding flow field, hence leading to changes in the noise generation as well as aerodynamic characteristics of the configuration. To accurately predict the noise level and study possible means for noise reduction, simplified configurations of two cylinders or rod and airfoil placed in tandem were studied widely in the literature. However, for tandem airfoil configurations, most related studies were focused on the aerodynamic performance. Studies on acoustic characteristics with numerical methods are relatively insufficient to the authors’ knowledge. There are several mechanisms of the sound generation of an airfoil in the flow field. Such noise can be divided into two categories: tonal noise and broadband noise. Broadband noise generated from the turbulence of flow is highly disordered. Therefore, the noise covers a broader spectrum of frequencies. The tonal noise is typically related to the flow feature with a certain stable periodic pattern. The noise sources of isolated airfoil include the following four items according to the paper by Migliore and Oerlemans[2]:

    a) Turbulent boundary layer-trailing edge noise.

    b) Laminar boundary layer-vortex shedding noise.

    c) Separation stall noise.

    d) Trailing-edge bluntness noise.

    The presence of the rear airfoil in the tandem configuration leads to interaction noise, and also self-noise of rear airfoil into the domain. The wake and the vortex shedding from the front airfoil will possibly impinge on the leading edge of the rear airfoil, hence generating interaction noise, which is affected by the geometry of both airfoils.

    This paper presents a computational study into the effects of changes in the geometry of the leading edge of the rear airfoil and trailing edge of the front airfoil could have on far-field noise. LES (Large Eddy Simulation) method from an in-house code is adopted in this paper to calculate the flow field. The far-field noise is obtained by FW-H method. The methodology used in the current paper is validated against benchmark cases to demonstrate its capability for noise prediction. NACA 0012-63 is chosen as the baseline airfoil in this paper. Mach number is kept constant as the takeoff and approach speed of aircraft is nearly 0.2 with a Reynolds number of 105. The combined effects of geometry variation from two airfoils are studied. The conclusions are drawn on the noise implications. This could be coupled with aerodynamic analysis in further studies.

    This paper is organized as follows: Section 1 presents the related work about the models and numerical methods. Section 2 introduces the methodology used for modeling the flow and calculating the sound. Validation is given by comparing results with benchmark cases. Section 3 presents the calculation result of the case using NACA 0012-63. Furthermore, a surrogate-based approach is used to study the geometry variations of the trailing edge of the front airfoil and the leading edge of the rear airfoil. Section 4 concludes the results obtained from the last section.

    1 Related work

    1.1 Typical models and configurations

    Tandem cylinder and rod-airfoil models have been investigated by many researchers using numerical and experimental methods. The widely mentioned tandem cylinder experiments were conducted in NASA Langley Research center at the Basic Aerodynamic Research Tunnel (BART) and Quiet Flow Facility (QFF). Experimental results revealed some flow features in this configuration with the two spaced at a diameter ratio of 1.435 and 3.7[3], respectively. Fifteen tandem cylinder solutions on the First Workshop on Benchmark problems for Airframe Noise Computations were summarized in the paper by Lockard[4]. There were also other studies related to the configuration, many of which are focused on flow control technology.

    Rod-airfoil interaction noise generated at low Mach number was investigated numerically by Casalino[5]with 2D URANS method. The spanwise effects were introduced into the time-domain formulation of the FW-H analogy by a statistical model. The results with spanwise correction compared very well to the accompanying experiment. Jacob conducted a low Mach number test on rod-airfoil in the large anechoic wind tunnel of the Ecole Centrale de Lyon, and a 2D URANS combined with a spanwise stochastic model was compared with the experiment[6]. Two CAA (Computational AeroAcoustics) methods were compared and validated with Jacob’s experiment by Schell[7]. A direct noise calculation with compressible DES (Detached Eddy Simulation) and a hybrid method with DES and Lighthill analogy were applied with StarCCM+. The calculations were three dimensional with a span of 0.3 m, which were the same as in the experiment. An incompressible LES coupled with Lighthill’s equation was used to predict the noise of rod-airfoil by Eltaweel[8]and again compared with Jacob’s[6]experiment. Most of these works identified that the leading edge is the dominant sound source.

    The effects of the geometry of isolate airfoil on the turbulent noise were studied numerically by Gill[9]. The dominant noise reduction mechanism was found to be related to the thickness and leading-edge radius of NACA series airfoil. The result was compared with an analytic flat plate that over-predict the noise by 3 dB at high frequencies. Paruchuri[10]extended Gill’s work to experiment with airfoil printed by 3D printer. Their work focused on the interactions of an isolated airfoil with the gust of a specific wavelength.

    As to the tandem airfoil configuration, most of the related studies were based on the background application in turbomachinery. McGlumphy[11]investigated the feasibility of using a tandem rotor in the rear stage of a core compressor. Laurent[12]studied the aerodynamic performance of tandem airfoil under a high angle of attacks. Shirsath[13]explored the effect of the angle of attack of the rear airfoil on unsteady aerodynamic quantities. It is to the author’s best knowledge that the aeroacoustics characteristics of two airfoils in a tandem configuration have not been sufficiently studied in the reported literature.

    It can be seen from the above that most of the researches focused on the flow condition and the geometry effects of an isolated airfoil. And the turbulence interaction is modeled with isotropic turbulence or generated by the wake of a rod. The wake of an airfoil and the geometry effect on interaction noise have not been studied numerically yet. In this paper, the background application the configuration corresponds to the take-off and approach flight conditions of aircraft. The flow conditions in the calculation stay unchanged while geometries of the two airfoils were varied.

    1.2 Numerical study

    In a typical CAA approach, three distinct phases are often used, including the calculation, extraction and characterization of the noise source, propagation calculation and analysis of noise levels at the receiving points. Accurate identification of noise sources using either numerical and experimental methods is a challenge. Numerical methods have a long term potential value as these can be integrated into the design process to study the effects of relevant factors such as geometry, configuration, and the effectiveness of noise reduction measures. It has become common practice to use RANS in aerodynamic design in industry settings on complex configurations; this has not been possible for noise calculations on realistic configurations due to the extremely high computational cost. Significant efforts are still required from the three areas of noise generation, propagation and characterization of observed noise. The current study focused on the problem of broadband noise generation for tandem airfoil configurations in order to further the state-of-the-art in the fields, development of hybrid approach appropriate for the simulation, improved understanding, characterization of noise generation mechanisms and more efficient numerical algorithms on massively parallel high-performance computing facilities.

    Jacob[5]also compared 2D URANS with 3D URANS and concluded that though 2D URANS over predicted the shedding frequency, the 3D URANS did not improve the shedding frequency significantly because of the intrinsic deficiency of the URANS formulation. Yi Jiang studied the rod-airfoil interaction noise with high order LES method[14]. The spacing between the rod and airfoil was studied as well. It was found that three modes of interaction exist for the variation of spacing.

    Liu[15]studied the trailing edge serration technology on the interaction noise of tandem airfoil without looking at other geometry variations. The chosen airfoil is the NACA-65 series, which was widely used in the propeller design. The aft airfoil is placed in the mirrored vertical direction. Chen[16]investigated the geometry optimization of a tandem airfoil at a low Mach number on aerodynamic performance. The geometry and relative positions are varied with 23 design variables in total.

    2 Numerical methods and validation

    2.1 Sound source and acoustic analogy model

    Far-field noise is obtained by the acoustic analogy method based on a two-step approach. First, the unsteady flow field around the sound source is simulated with CFD method; Second, the scattering of the sound wave is calculated by solving the wave equation. In this paper, the flow field for acoustic calculation is obtained from the result of LES method. Far-field noise is calculated with the sound analogy model FW-H.

    In the current work, LES is used to investigate the noise mechanism of tandemairfoil. A fifth-order weighted essentially non-oscillatory (WENO) scheme is used in the study. Roe method is used for the flux splitting. The temporal integration used is LU-SGS approach. The subgrid-scale (SGS) turbulence model used in this paper is the Smagorinsky-Lilly model. Far-field boundary condition is obtained with the method based on Riemann invariant from the 1-D inviscid flow. Side boundary is set as a symmetric boundary condition.

    2.2 Validation using turbulent rod configuration

    The experimental data fromBoudet for turbulent rod configuration is used to validate the numerical methods[17]. This configuration is a cylinder flow with an inflow speedU0=72 m/s, temperatureT=293 K andρ=1.2 kg/m3. Reynolds number based on rod diameter is 46000. A circular computational domain with a radius of 80D is used. The computation mesh used is a 181×181×31 sized structured grid. The first layer thickness is chosen as 0.05D(Dis the diameter of the rod andD=0.01m). A growth rate of 1.0359 is used.

    The first layer thickness in this paper is varied to investigate its effect on far-field noise. Three values were chosen for the first layer thickness, 0.05D, 0.013Dand 0.005D, which correspond toy+(yplus value) of 3.8, 1 and 0.38, respectively. Circumferential points at Case 2 is increased to study the effect ofx+(relative to the definition ofy+). The non-dimensional time step is set as Δt=0.002. The total flow time is chosen as 100,000 time steps. Acoustic data is sampled every time step, and the sampling frequency is 360 kHz. Table 1 shows the aerodynamic results of different first layer thickness. It can be seen that lowery+value gives a better shedding frequency. The drag coefficients and lift coefficients from these cases did not show much difference.

    Figure 1 shows the power spectra density (PSD) of pressure fluctuation in 70D. Seo[18]obtained the result from a two-dimensional LPCE (Linearized Perturbed Compressible Equations) computation with input from LES result. A sixth-order finite difference scheme was used. Strohal number is 0.19 in the experiment. It can be seen that controlling they+value at 1 can obtain better frequency result in comparison with the results obtained by others. For case 2, though the results slightly over predicted the peak frequency, the value agrees reasonably well with the aerodynamic result. It should be noted that only numerical results of pressure at location (0, 70D) was found in the published paper. Therefore, the overprediction at high frequency may result from that smaller scales flow structure was resolved by the finer grid.

    表1 不同網格的圓柱氣動結果Table 1 Aerodynamic results of rod using different grids

    圖1 點(0, 70D)處的壓力脈動自相關頻譜圖Fig.1 Pressure fluctuation auto-correlation spectra at a receiving point of (0,70D)

    Figure 3 depicts the span-wise coherence function of surface pressure at the given shedding frequencyf0. It can be seen that the harmonic frequency of shedding frequency decay more slowly than other frequencies. This means the span length is enough for capturing the flow scale apart from those corresponding shedding frequency. Therefore, corrections should be made on the final result. The correction method used here is introduced by Seo[18].

    圖2 點(0, 70D)處的壓力脈動自相關頻譜圖(放大)Fig.2 Enlarged auto-correlation spectra of pressure at receiving point of (0, 70D)

    (a) Spanwise coherence coefficient of surface pressure at Shedding frequency over the span (90°)

    (b) Two-point coherence function of surface pressure over the rod at 90° (z=0 and z=1.5D)

    Far-field noise calculated from the FW-H surface overlapped with the rod surface is shown in Figure 4. It can be seen that both the scale and frequency of the tonal noise are well predicted. The harmonics of the main tonal noise are also captured. The high fluctuation of the spectra is the common phenomenon that appears in the calculation of noise, which results from the relatively short length of the total sampling time.

    圖4 觀測點(0,1.85 m)處的聲壓級頻譜圖Fig.4 Sound pressure level at receiver point (0,1.85 m)

    2.3 Validation using rod-airfoil configuration

    The same rod airfoil configuration, as used by Jacob[6], is considered in this validation, as depicted in Figure 5. Freestream flow conditions are set asU0=72 m/s,T=293 K andρ=1.2 kg/m3. The reference chord length is chosen asc=100 mm, which results in a Reynolds number of 4.8×105.

    Three grids were compared to validate the meshused for the following studies.yplus values of the first layer in all three grids were kept at 1. The differences were the growth ratio in the normal direction and the number of points along the circumference direction. Table 2 shows the nodes number allocated on the surface of the rod and airfoil.

    Total computational domain is set as (-8c,8c) in vertical direction and (-10c,20c) in the horizontal direction. Computational methods are the same as the rod configuration but the non-dimensional time step is reduced toΔt=0.001. The sampling time for acoustic data was 20,000 time-steps and sampled every time step. Sampling frequency was then defined as 720 kHz. Results from RANS of 20,000 iterations were used as the initial values for LES calculation.

    The calculation results are shown in Figure 6. Results of three cases are compared with reference LES and experiment results. Turbulent intensity is obtained based on streamwise velocity fluctuations. Considering the rod was shifted 2 mm vertically in the experiment,the velocity profile at position A (x/c=-0.255) shows the greatest velocity deficit at the locationy/c>0. Apart from the result of grid No.2, all simulated cases exhibit a larger averaged velocity deficit compared to experimental results. And turbulent intensities are also larger along the middle line. Despite these differences, it is believed that the results obtained here are acceptable. These phenomena appear in nearly all the numerical simulations of rod airfoil configuration.

    圖5 圓柱機翼構型示意圖[6]Fig.5 Rod-airfoil configuration for validation[6]

    By comparing the calculated results of the two cases above, it can be seen that the method used in this paper is feasible in the current study on the geometry variation of tandem airfoil configurations.

    表2 三個圓柱機翼網格的細節(jié)參數(shù)Table 2 Mesh details of three cases

    (a) Averaged meanstreamwise velocity profile at x/c=-0.255

    (a) Averaged meanstreamwise velocity profile at x/c=0.25

    圖8 機翼前緣上方(0,1.85 m)處的遠場SPL頻譜圖Fig.8 Far-field SPL spectra at the location (0,1.85 m) above the airfoil leading edge

    3 Aerodynamic and aeroacoustics calculation of tandem airfoil

    3.1 Geometry and computational model

    The tandem airfoil model used here is a model used in the 1st International Workshop on High-Order CFD Methods[19], as shown in Figure 9. Thedoffis the distance of vertical displacement relative to the chordc.doffis set as 0.24 in this paper and the angle of attack of the front airfoil is set as 10°. Parameterdsepis the horizontal distance between the two airfoils. The airfoil used is NACA0012 and the methods used for generating the airfoil are described in the following.

    圖9 計算所選用的串聯(lián)機翼模型示意圖Fig.9 Tandem airfoil configuration chosen for computation

    Geometry variations are introduced for the leading edge radius and the maximum thickness of the aft airfoil. As for the front airfoil, maximum thickness and boat-tail angle are varied. Class/Shape Transformation method[20]is used to generate varied geometry. Leading-edge radius is related to the maximum thickness through the following expression:

    (1)

    wheretis the maximum thickness of airfoil, andIstands for non-dimensional parameter that controls the shape of leading edge. The standard NACA 4-series airfoil is defined withI=6.

    The model is shown as Figure 9. In this paper, the relative position of two airfoils did not changed, only geometries were varied. The chosen parameters are listed below in Table 3.

    表3 串聯(lián)機翼構型的外形參數(shù)Table 3 List of parameters for tandem airfoil configuration

    Two NACA0012 airfoils with the same chord ofc=100 mm were placed in tandem configuration with twice the chord length apart between the leading edge of the two airfoils. The front airfoil was rotated around point quarter-chord (0.25c, 0) to achieve an angle of attack of 10°.

    O-type blocks are used to discretize the computational domain.Each airfoil surface is meshed with 538 nodes along the circumferential direction. The normal direction growth ratio is set as less than 1.05 near the wall. Total nodes number is 128 in the normal direction near the wall.

    It should be notedthat the current noise result is calculated from the finite span 0.2c, as the length of span would directly influence the total noise results. Spanwise length correction has been made for the tandem airfoil cases.

    Q-criterion Iso-surfaces of tandem airfoil calculated using LES are shown in Figure 10. Periodic pressure fluctuations are generated from the leading edge shedding vortex on the upper face of the front airfoil. Thick low-speed zone occurs in the wake domain of the trailing edge. Then the vortices are pushed downward by the presence of the aft airfoil, impinging on the upper face of the aft airfoil. Finally, the vortices are mixed with the wake of aft airfoil downstream of the aft airfoil.

    圖10 Q準則等值面圖(速度云圖)Fig.10 Q-criterion iso-surfaces for velocity

    The turbulence flow is mainly developed at the upper face and trailing edge of the front airfoil, then the upper face of the rear airfoil. These areas are likely to be the primary regions for noise sources. The SPL (sound pressure levels) spectra are shown in Figure 11. Three FW-H integration surfaces comprised of the front airfoil and aft airfoil, and tandem airfoil were used to obtained the far-field results of the three lines, respectively. It can be seen that aft airfoil is the dominant noise source in the overall model. The maximum SPL frequency of aft airfoil in the tandem configuration is totally different from the isolated condition, which has a max SPL frequency of about 2700 Hz[21]. This is similar to the rod airfoil configuration. Therefore, the following work in this paper focused on the trailing edge of the front airfoil and the leading edge of aft airfoil.

    圖11 接收點(0,1.85 m)處的串聯(lián)機翼遠場聲壓級頻譜圖Fig.11 Sound pressure level of the tandem airfoil at the receiver (0, 1.85 m)

    3.2 Surrogate-based study of shape parameters

    This section adopted aresponse surface methodology to study the combined effects of variations in these geometric features on noise levels[22]. A Latin hypercube sampling method is used to generate a set of sample points. A series of Kriging based surrogate models were generated in order to find the best combinations of these geometry parameters regarding noise level reduction at the observation points. The response surface will demonstrate the combined effects of these geometry features and provide an insight into the design of future lower noise configurations.

    Four parameters were chosen in the study of the combined effect, including maximum thickness and boat-tail angle of the front airfoil, leading-edge radius and the maximum thickness of the aft airfoil. The ranges of these four parameters are listed in Table 4.

    表4 幾何外形參數(shù)變化范圍Table 4 Range of variations for geometry parameters

    A total of 30 sample points were obtained for the generation of the response surface. TheR2value is defined as the following to describe the accuracy of the model.

    (2)

    As Figure 12(a) shows, OASPL is increased with the maximum thickness when the boat-tail angle is low. Moreover, the two curves become closer when the maximum thickness is greater than 15% of the chord. That means that a low boat-tail angle has very few effects on OASPL when the maximum thickness is large. Judging from the different trend at low and high boat-tail angle, it is reasonable to presume that there is a certain range of ratios between maximum thickness and boat-tail angle that give the lowest sound emission. This ratio makes sure that there is a smooth transition from the maximum thickness point to the trailing edge. Figure 12(b) shows a similar trend as in Figure 12(a) apart from the fact that two curves are overlapped.

    (a) The effect of maximum thickness at different boat-tail angles

    (b) The effect of boat-tail angle

    (c) The effect of leading edge radius

    (d) The effect of maximum-thickness for different airfoils

    As to Figure 12(c), OASPL increases withI(leading edge radius) at small maximum thickness. Sharper leading edge can reduce the direct impingement of vortex on the surface of the aft airfoil, driving the vortex downstream more smoothly. While at larger maximum thickness, OASPL is smaller with largerI. In this case, the front of aft airfoil in inevitably large, leading to higher blockage at the leading edge. Sharper leading edge would make the pressure at the static zone higher. Whereas rounded leading-edge could relieve this phenomenon. Figure 12(d) shows a similar trend as Figure 12(c).

    The response surface is plotted with HAT (Hierarchical Axes Technology) graph. The effects of four parameters on OASPL are plotted. Variables in the outer axis are front airfoil maximum thickness in thexdirection and boat-tail angle in theydirection. In the inner coordinate, two variables are aft airfoil leading-edge radius in thexdirection and maximum thickness in theydirection.

    Far-field OASPL of the front airfoil is shown in Figure 13, wheretmaxmeans front airfoil thickness, andβmeans front airfoil boat-tail angle. It can be seen that front airfoil OASPL is less sensitive to the effect of aft airfoil geometry. The highest OASPL lies in the subplot (1,3) and (3,1). The counterpart front airfoil geometry is the combination of smalltmaxwith smallβ, and largetmaxwith largeβ. When an airfoil is featured with slim rear shape, namely low boat-tail angle, noise is decreased with the increase of maximum thickness. This conclusion is similar to that of James[9]. But the trend is reversed when the boat-tail angle is high. In this case, geometry has a blunt trailing edge. Because the scale of vortex shedding from the trailing edge is increased, making the airfoil thickness larger may worsen this phenomenon.

    圖13 接收點(0,1.85 m)接收到的前翼OASPL圖

    Figure 14 depicts the overall noise of the front airfoil and aft airfoil. It is obvious that tandem airfoil OASPL is more sensitive to the aft airfoil geometry. Two trendscan be observed from the HAT graph. The first one is that the trend is totally opposite in subplot (3,1) and (1,3). When front airfoiltmaxandβare small, tandem airfoil noise is larger if aft airfoiltmaxorIbecome larger. The second is that in the rest of the subplot, OASPL is lower when aft airfoil maximum thickness and leading-edge radii are either smaller or larger.

    The primary areas with low OASPL lie on the line between Figure 14 subplot (1,3) and (3,1). But the counterpart of this area in Figure 13 is the lowest noise areas of the front airfoil. It should be noted that aft airfoil takes the dominant position in tandem airfoil configuration. Figure 13 and Figure 14 show the opposite trend in subplot (3,1) and (1,3), since the response surfaces were generated from single front airfoil and tandem airfoil, respectively. The effect of the front airfoil is overwhelmed by aft airfoil.

    Figures 15 and 16 show the frequency obtained at max SPL. It is easy to find that the frequency has a very similar trend as max SPL. That means louder noise is accompanied at higher frequencies.

    圖15 接收點(0,1.85 m)接收到的前翼最大聲壓級對應的頻率圖Fig.15 Max SPL frequency of front airfoil at the receiver (0,1.85 m)

    圖16 接收點(0,1.85 m)接收到的串聯(lián)機翼最大聲壓級對應的頻率圖Fig.16 Max SPL frequency of tandem airfoil at the receiver (0,1.85 m)

    4 Conclusions

    Far-field noise characteristics of typical tandem airfoil configuration have been studied numerically with 3D LES combined with FW-H analogy. Geometric features of both airfoils are varied simultaneously to study the combined effects on noise. The results of four geometry parameters: front airfoil maximum thickness, front airfoil boat-tail angle, aft airfoil leading-edge radius and aft airfoil maximum thickness, have been investigated. The key findings are listed as follows:

    1) Geometries of both airfoils in tandem airfoil configuration should be considered for noise reduction. Since the reduction of noise from the front airfoil alone could worsen the overall noise level of the tandem airfoil configuration.

    2) The ratio between maximum thickness and boat-tail angle of front airfoil should bekept in a certain medium range. Too large or too small would worsen the aeroacoustics performance. This can avoid the sharp change in geometry that may induce vortex generation. Keeping the transition of flow from the maximum thickness point toward trailing edge could also avoid large trailing edge vortex.

    3) Sharp leading edge should only be used in a thin airfoil. Larger leading-edge radius should be used under the circumstance that maximum thickness islarge. Namely,Iandtmaxshould be both large or small at the same time in most cases.

    4) Noise frequency of max SPL has a similar trend as the max SPL. The frequency for the Max SPL should be carefully considered in order to avoid undesirable effects such as acoustic fatigue.

    The complex interactions between two airfoils in a tandem configuration make it necessary to consider the case using an integrated approach. The current study should be further extended to include spanwise effects for the same tandem configuration in the next steps of the research, along with comparative studies using an experimental approach. In reality, aerodynamic aspects should also be considered, which further complicates the problem.

    From the results, it is also suggested that configuration parameters such as the vertical and horizontal distance, and spanwise variations, should be included in the future, three-dimensional studies.

    猜你喜歡
    接收點遠場聲壓級
    機器噪聲平均聲壓級計算方法差異性實證研究
    電動工具(2024年1期)2024-02-29 01:40:24
    一種計算消聲室聲壓級的新方法
    計量學報(2020年2期)2020-04-11 04:33:22
    全新DXR mkll有源揚聲器
    演藝科技(2019年4期)2019-03-30 03:21:46
    基于仿真與實測的列車遠場氣動噪聲分析
    測控技術(2018年9期)2018-11-25 07:44:12
    更正
    動態(tài)網絡最短路徑射線追蹤算法中向后追蹤方法的改進*1
    某種陣列雷達發(fā)射通道遠場校準簡易方法
    淺海波導界面對點源振速方向的影響?
    應用聲學(2015年3期)2015-10-27 02:52:49
    Diodes1.9W D類音頻放大器提供高聲壓級水平并延長電池壽命
    戰(zhàn)斗部遠場水下爆炸對艦船沖擊損傷評估
    午夜精品久久久久久毛片777| 国产精品 欧美亚洲| 精品无人区乱码1区二区| 中亚洲国语对白在线视频| 亚洲黑人精品在线| 十八禁人妻一区二区| 精品免费久久久久久久清纯| 亚洲五月婷婷丁香| 操出白浆在线播放| 国产精品 欧美亚洲| 欧美日韩瑟瑟在线播放| 日本欧美视频一区| 亚洲精品美女久久av网站| 精品久久久久久久毛片微露脸| 后天国语完整版免费观看| 国产97色在线日韩免费| 一卡2卡三卡四卡精品乱码亚洲| 狂野欧美激情性xxxx| 亚洲一区高清亚洲精品| 国产精品秋霞免费鲁丝片| 欧美成人免费av一区二区三区| 成人三级做爰电影| 老司机福利观看| 老司机午夜十八禁免费视频| 亚洲第一欧美日韩一区二区三区| 国产亚洲av嫩草精品影院| 午夜激情av网站| 激情在线观看视频在线高清| 咕卡用的链子| 成年女人毛片免费观看观看9| 少妇的丰满在线观看| 日本在线视频免费播放| 中文字幕高清在线视频| 亚洲精品一区av在线观看| 日韩 欧美 亚洲 中文字幕| 国产在线观看jvid| 国产高清激情床上av| 91在线观看av| 91老司机精品| 九色国产91popny在线| 91精品国产国语对白视频| 丝袜在线中文字幕| 亚洲精品粉嫩美女一区| 国内精品久久久久久久电影| 国产免费男女视频| 久久国产乱子伦精品免费另类| 亚洲人成77777在线视频| 女警被强在线播放| 久久精品国产综合久久久| 18禁观看日本| 露出奶头的视频| 99久久国产精品久久久| 欧美乱色亚洲激情| 在线观看免费视频网站a站| 高清黄色对白视频在线免费看| 欧美中文日本在线观看视频| 国产乱人伦免费视频| 久久亚洲精品不卡| 久久久久久久久久久久大奶| 精品第一国产精品| ponron亚洲| 美女高潮喷水抽搐中文字幕| 大型av网站在线播放| 亚洲五月天丁香| 午夜精品在线福利| 一级作爱视频免费观看| 91精品国产国语对白视频| 男人舔女人的私密视频| 亚洲狠狠婷婷综合久久图片| 91麻豆av在线| 午夜激情av网站| 日韩欧美一区二区三区在线观看| 久久天躁狠狠躁夜夜2o2o| 欧美国产日韩亚洲一区| 欧美一级a爱片免费观看看 | 午夜成年电影在线免费观看| 麻豆成人av在线观看| 婷婷六月久久综合丁香| 久久九九热精品免费| 久久久久久久久免费视频了| 伊人久久大香线蕉亚洲五| 91成年电影在线观看| 亚洲人成伊人成综合网2020| 韩国av一区二区三区四区| 涩涩av久久男人的天堂| 免费人成视频x8x8入口观看| 两个人看的免费小视频| 久久久久久人人人人人| 性色av乱码一区二区三区2| 亚洲成国产人片在线观看| 久久精品亚洲熟妇少妇任你| 男女下面进入的视频免费午夜 | 一区在线观看完整版| 亚洲色图综合在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 18禁裸乳无遮挡免费网站照片 | 韩国av一区二区三区四区| 日韩成人在线观看一区二区三区| 高清黄色对白视频在线免费看| 国产成人啪精品午夜网站| 变态另类成人亚洲欧美熟女 | 午夜日韩欧美国产| 亚洲 国产 在线| 男女床上黄色一级片免费看| 大陆偷拍与自拍| 免费观看人在逋| www.www免费av| 操美女的视频在线观看| 久久久久久久久免费视频了| 美国免费a级毛片| 国产精品国产高清国产av| 亚洲欧美日韩无卡精品| 电影成人av| 色av中文字幕| 啦啦啦观看免费观看视频高清 | 黄色视频,在线免费观看| 十分钟在线观看高清视频www| 超碰成人久久| 日韩欧美一区二区三区在线观看| 午夜精品久久久久久毛片777| 在线国产一区二区在线| 亚洲精品在线观看二区| 男女做爰动态图高潮gif福利片 | 午夜福利在线观看吧| 久久九九热精品免费| 国产高清videossex| 精品国产超薄肉色丝袜足j| 久久香蕉国产精品| 亚洲欧美精品综合久久99| 97人妻天天添夜夜摸| 99热只有精品国产| 国产亚洲欧美98| 国产精品国产高清国产av| 久久天堂一区二区三区四区| 亚洲人成电影免费在线| 激情在线观看视频在线高清| 亚洲精品国产精品久久久不卡| 性少妇av在线| 手机成人av网站| 国产精品免费一区二区三区在线| 日本vs欧美在线观看视频| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲欧美激情综合另类| 亚洲 国产 在线| 亚洲av第一区精品v没综合| 欧美成人免费av一区二区三区| 黑人巨大精品欧美一区二区mp4| 亚洲男人天堂网一区| 国产麻豆成人av免费视频| 欧美绝顶高潮抽搐喷水| 亚洲男人天堂网一区| 非洲黑人性xxxx精品又粗又长| 亚洲一区二区三区不卡视频| 国产成年人精品一区二区| 国产国语露脸激情在线看| 一本综合久久免费| 国产成年人精品一区二区| av网站免费在线观看视频| 国产精品秋霞免费鲁丝片| 757午夜福利合集在线观看| 他把我摸到了高潮在线观看| 久久亚洲真实| 欧美一级毛片孕妇| 亚洲欧美日韩高清在线视频| 老司机深夜福利视频在线观看| 一本综合久久免费| 亚洲中文字幕一区二区三区有码在线看 | 亚洲国产欧美日韩在线播放| 黄色视频,在线免费观看| 免费久久久久久久精品成人欧美视频| 国产一级毛片七仙女欲春2 | 午夜成年电影在线免费观看| 99国产极品粉嫩在线观看| 亚洲aⅴ乱码一区二区在线播放 | 丰满人妻熟妇乱又伦精品不卡| 国产成人欧美| 欧美在线一区亚洲| 波多野结衣巨乳人妻| 亚洲第一av免费看| 亚洲av日韩精品久久久久久密| 亚洲色图 男人天堂 中文字幕| 久久精品国产99精品国产亚洲性色 | 女警被强在线播放| 国产高清videossex| 国产精品一区二区精品视频观看| 国产一区二区三区综合在线观看| 国产激情久久老熟女| 悠悠久久av| 欧美日韩乱码在线| 国产高清videossex| 国产熟女xx| 老司机靠b影院| 美女国产高潮福利片在线看| 亚洲男人天堂网一区| 桃色一区二区三区在线观看| 免费看十八禁软件| 国产免费男女视频| 琪琪午夜伦伦电影理论片6080| 亚洲一区高清亚洲精品| 国产伦一二天堂av在线观看| 神马国产精品三级电影在线观看 | 欧美一级毛片孕妇| 看片在线看免费视频| av天堂在线播放| 国产精品爽爽va在线观看网站 | 色综合站精品国产| 国产精品久久视频播放| 成人国语在线视频| 亚洲国产欧美网| 精品电影一区二区在线| 国产欧美日韩精品亚洲av| 久9热在线精品视频| 欧美 亚洲 国产 日韩一| 在线免费观看的www视频| 成人国语在线视频| 午夜福利高清视频| 青草久久国产| 亚洲国产精品合色在线| 欧美激情极品国产一区二区三区| 麻豆国产av国片精品| 日本在线视频免费播放| 伊人久久大香线蕉亚洲五| 中亚洲国语对白在线视频| 亚洲国产欧美网| 18禁国产床啪视频网站| 午夜福利成人在线免费观看| 乱人伦中国视频| 久久精品91无色码中文字幕| 超碰成人久久| 国产黄a三级三级三级人| 欧美日韩亚洲国产一区二区在线观看| 悠悠久久av| 久久久久久久午夜电影| 18禁黄网站禁片午夜丰满| 国产av一区在线观看免费| 精品少妇一区二区三区视频日本电影| 欧美激情 高清一区二区三区| 嫁个100分男人电影在线观看| 国内久久婷婷六月综合欲色啪| 国产亚洲av高清不卡| 欧美激情 高清一区二区三区| 好男人在线观看高清免费视频 | 国产精品1区2区在线观看.| or卡值多少钱| 脱女人内裤的视频| 他把我摸到了高潮在线观看| 精品欧美一区二区三区在线| 亚洲伊人色综图| 午夜福利视频1000在线观看 | 午夜福利视频1000在线观看 | 99在线视频只有这里精品首页| 亚洲午夜精品一区,二区,三区| 成人免费观看视频高清| 我的亚洲天堂| av欧美777| 中国美女看黄片| 欧美一级毛片孕妇| www日本在线高清视频| 国产精品香港三级国产av潘金莲| 一级,二级,三级黄色视频| 精品国产国语对白av| 亚洲男人天堂网一区| 欧美日韩亚洲国产一区二区在线观看| 欧美av亚洲av综合av国产av| 18禁黄网站禁片午夜丰满| 国产亚洲精品久久久久久毛片| 日韩欧美免费精品| 亚洲色图综合在线观看| 十八禁网站免费在线| 中出人妻视频一区二区| av中文乱码字幕在线| 欧美成狂野欧美在线观看| 国产精品国产高清国产av| 国产成人精品在线电影| 91精品国产国语对白视频| 亚洲,欧美精品.| 欧美一级毛片孕妇| 亚洲国产毛片av蜜桃av| 性少妇av在线| av天堂在线播放| www.www免费av| 精品久久久久久久久久免费视频| 成人特级黄色片久久久久久久| 亚洲成av片中文字幕在线观看| 熟女少妇亚洲综合色aaa.| 妹子高潮喷水视频| 露出奶头的视频| 国产高清videossex| 丝袜在线中文字幕| 国产av精品麻豆| 男人操女人黄网站| 亚洲av成人一区二区三| 婷婷精品国产亚洲av在线| 日本免费a在线| 久久午夜综合久久蜜桃| 两个人免费观看高清视频| 十分钟在线观看高清视频www| 999久久久国产精品视频| 给我免费播放毛片高清在线观看| 亚洲色图av天堂| 制服丝袜大香蕉在线| 真人做人爱边吃奶动态| 久久久久久国产a免费观看| 欧美久久黑人一区二区| 久久精品国产清高在天天线| 正在播放国产对白刺激| 亚洲视频免费观看视频| 亚洲国产中文字幕在线视频| 一a级毛片在线观看| 在线av久久热| 欧美av亚洲av综合av国产av| 少妇的丰满在线观看| 99国产极品粉嫩在线观看| 在线观看日韩欧美| 一区二区三区高清视频在线| 99香蕉大伊视频| 亚洲精品一卡2卡三卡4卡5卡| 欧美激情 高清一区二区三区| 首页视频小说图片口味搜索| 国产免费男女视频| 伊人久久大香线蕉亚洲五| 日韩三级视频一区二区三区| 国内精品久久久久久久电影| 曰老女人黄片| 久久婷婷人人爽人人干人人爱 | 欧美精品啪啪一区二区三区| 香蕉丝袜av| 亚洲精品av麻豆狂野| 亚洲熟妇中文字幕五十中出| 黄色a级毛片大全视频| 少妇被粗大的猛进出69影院| 精品久久久久久成人av| 啦啦啦免费观看视频1| 欧美丝袜亚洲另类 | 日本撒尿小便嘘嘘汇集6| 久久人妻av系列| 天天躁夜夜躁狠狠躁躁| 亚洲av电影在线进入| 国产aⅴ精品一区二区三区波| 国产精品香港三级国产av潘金莲| 色av中文字幕| 久热这里只有精品99| 搡老熟女国产l中国老女人| 精品国产美女av久久久久小说| 国产亚洲欧美精品永久| 日韩欧美一区视频在线观看| 亚洲成国产人片在线观看| 亚洲精品在线观看二区| 最新在线观看一区二区三区| 村上凉子中文字幕在线| 久久草成人影院| 国产三级在线视频| 中出人妻视频一区二区| 亚洲国产精品999在线| 夜夜夜夜夜久久久久| av中文乱码字幕在线| aaaaa片日本免费| 国产精品久久久久久亚洲av鲁大| 国产成人av教育| 欧美+亚洲+日韩+国产| 免费在线观看黄色视频的| 每晚都被弄得嗷嗷叫到高潮| 国产麻豆成人av免费视频| 亚洲专区中文字幕在线| 美女 人体艺术 gogo| 午夜精品国产一区二区电影| 男人操女人黄网站| 亚洲全国av大片| 欧美日韩一级在线毛片| 大型黄色视频在线免费观看| 精品国产一区二区久久| 好男人电影高清在线观看| 无人区码免费观看不卡| 欧美日本亚洲视频在线播放| 少妇的丰满在线观看| 一级a爱片免费观看的视频| 久久久久国产一级毛片高清牌| 久久国产精品男人的天堂亚洲| 十分钟在线观看高清视频www| 久久人人精品亚洲av| 免费不卡黄色视频| 久久国产精品男人的天堂亚洲| 精品午夜福利视频在线观看一区| 日韩精品免费视频一区二区三区| 高清在线国产一区| 国产精品一区二区三区四区久久 | 熟妇人妻久久中文字幕3abv| 久久精品国产亚洲av高清一级| 天天一区二区日本电影三级 | 悠悠久久av| 亚洲人成电影免费在线| 国产av在哪里看| 亚洲一区高清亚洲精品| av在线播放免费不卡| 国产激情久久老熟女| 国产av一区二区精品久久| 正在播放国产对白刺激| 免费高清视频大片| 91精品国产国语对白视频| 淫秽高清视频在线观看| 国产精品亚洲美女久久久| 日本a在线网址| 欧美激情久久久久久爽电影 | 天天躁夜夜躁狠狠躁躁| 欧美黄色淫秽网站| 午夜福利欧美成人| 伦理电影免费视频| АⅤ资源中文在线天堂| 看黄色毛片网站| 久久久精品欧美日韩精品| 国产成+人综合+亚洲专区| 国产亚洲精品久久久久5区| 在线十欧美十亚洲十日本专区| 狠狠狠狠99中文字幕| 国语自产精品视频在线第100页| 18美女黄网站色大片免费观看| 国产亚洲精品久久久久久毛片| 欧美乱妇无乱码| 亚洲国产欧美日韩在线播放| 日本撒尿小便嘘嘘汇集6| 美女大奶头视频| 国产亚洲欧美精品永久| 亚洲,欧美精品.| 国产99久久九九免费精品| 久久久久久大精品| 99国产精品一区二区三区| 免费人成视频x8x8入口观看| 夜夜看夜夜爽夜夜摸| 在线观看日韩欧美| 精品久久久久久,| 国产精品久久电影中文字幕| 人人妻人人爽人人添夜夜欢视频| 热re99久久国产66热| 国产精品野战在线观看| 国产黄a三级三级三级人| 午夜日韩欧美国产| 大型av网站在线播放| 国产视频一区二区在线看| 欧美色视频一区免费| 亚洲av成人不卡在线观看播放网| 成年女人毛片免费观看观看9| 国产高清激情床上av| 国产精品一区二区免费欧美| 国产精品九九99| 咕卡用的链子| 最近最新免费中文字幕在线| 12—13女人毛片做爰片一| aaaaa片日本免费| 国产不卡一卡二| 亚洲精品国产色婷婷电影| 一二三四社区在线视频社区8| 91精品国产国语对白视频| 少妇裸体淫交视频免费看高清 | 黄色丝袜av网址大全| 99国产精品免费福利视频| 国产1区2区3区精品| 国产av精品麻豆| 亚洲av五月六月丁香网| 男人操女人黄网站| 亚洲成人精品中文字幕电影| 黄色丝袜av网址大全| 天堂√8在线中文| 真人做人爱边吃奶动态| 国产成人精品久久二区二区91| 亚洲欧美激情在线| 国产精品秋霞免费鲁丝片| 国产一区在线观看成人免费| 久久影院123| 欧美日韩中文字幕国产精品一区二区三区 | 久久人人97超碰香蕉20202| 国产单亲对白刺激| 精品乱码久久久久久99久播| 国产亚洲精品综合一区在线观看 | 国产精品av久久久久免费| 久久久久国产精品人妻aⅴ院| 一本久久中文字幕| 午夜免费成人在线视频| 国产精品 欧美亚洲| 欧美日本中文国产一区发布| 免费在线观看视频国产中文字幕亚洲| 18禁裸乳无遮挡免费网站照片 | 一区福利在线观看| 看免费av毛片| 狠狠狠狠99中文字幕| 黄网站色视频无遮挡免费观看| 国产成人av教育| 国产又色又爽无遮挡免费看| 黄色a级毛片大全视频| 俄罗斯特黄特色一大片| 国产精品亚洲一级av第二区| x7x7x7水蜜桃| 日本 av在线| 精品国内亚洲2022精品成人| 一夜夜www| 欧美成人午夜精品| 老司机午夜十八禁免费视频| av视频免费观看在线观看| 免费在线观看日本一区| 精品免费久久久久久久清纯| 久久久久久免费高清国产稀缺| 操美女的视频在线观看| 91麻豆精品激情在线观看国产| 一区二区三区国产精品乱码| 99国产精品一区二区蜜桃av| 亚洲国产欧美日韩在线播放| 91麻豆av在线| 天堂√8在线中文| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品久久久久久精品电影 | 国产av在哪里看| 日本在线视频免费播放| 成人亚洲精品一区在线观看| 国产精品98久久久久久宅男小说| 国产麻豆成人av免费视频| 在线观看一区二区三区| 日日爽夜夜爽网站| av视频在线观看入口| 999精品在线视频| 国产精品久久久久久亚洲av鲁大| 欧美中文综合在线视频| av有码第一页| 777久久人妻少妇嫩草av网站| 中文字幕人妻熟女乱码| 男人舔女人下体高潮全视频| 窝窝影院91人妻| 69av精品久久久久久| 成人av一区二区三区在线看| 亚洲精品国产一区二区精华液| 变态另类丝袜制服| 色婷婷久久久亚洲欧美| 在线观看日韩欧美| 亚洲片人在线观看| 中文字幕高清在线视频| 久久精品成人免费网站| 精品人妻1区二区| 久久午夜综合久久蜜桃| 亚洲九九香蕉| 日韩中文字幕欧美一区二区| 桃色一区二区三区在线观看| 午夜福利,免费看| 日韩国内少妇激情av| 久久性视频一级片| 一级a爱片免费观看的视频| 日韩大尺度精品在线看网址 | 麻豆久久精品国产亚洲av| 大码成人一级视频| 亚洲精品粉嫩美女一区| 国产精品99久久99久久久不卡| 国产麻豆成人av免费视频| 久久婷婷人人爽人人干人人爱 | 国产精品 国内视频| 好男人电影高清在线观看| 国产精品永久免费网站| 日韩视频一区二区在线观看| 国语自产精品视频在线第100页| 无人区码免费观看不卡| 亚洲视频免费观看视频| 成人永久免费在线观看视频| 波多野结衣一区麻豆| av在线播放免费不卡| 波多野结衣巨乳人妻| 丝袜人妻中文字幕| 精品久久久久久,| 日韩大尺度精品在线看网址 | 免费在线观看影片大全网站| 精品久久久久久成人av| 亚洲av电影在线进入| 中文亚洲av片在线观看爽| 在线av久久热| 久久久水蜜桃国产精品网| 777久久人妻少妇嫩草av网站| 国产麻豆69| 91麻豆精品激情在线观看国产| 成人永久免费在线观看视频| 麻豆国产av国片精品| 日本 av在线| 91在线观看av| 如日韩欧美国产精品一区二区三区| 丁香欧美五月| 91精品国产国语对白视频| 日本免费a在线| 国产精品一区二区精品视频观看| 亚洲五月婷婷丁香| 此物有八面人人有两片| 桃色一区二区三区在线观看| 亚洲成人精品中文字幕电影| 色综合欧美亚洲国产小说| 在线永久观看黄色视频| 高潮久久久久久久久久久不卡| 国产亚洲欧美在线一区二区| 两性午夜刺激爽爽歪歪视频在线观看 | 国产伦一二天堂av在线观看| 日韩av在线大香蕉| 久热爱精品视频在线9| 69精品国产乱码久久久| 亚洲精品在线观看二区| 成人亚洲精品av一区二区| 久久精品aⅴ一区二区三区四区| 中文字幕最新亚洲高清| 免费在线观看完整版高清| 日日摸夜夜添夜夜添小说| 久久久久久大精品| 超碰成人久久| a在线观看视频网站| 12—13女人毛片做爰片一| 很黄的视频免费| 亚洲精品久久成人aⅴ小说| 午夜免费鲁丝| 叶爱在线成人免费视频播放| 国产精品99久久99久久久不卡| 国产欧美日韩一区二区精品| 午夜福利18| 丰满人妻熟妇乱又伦精品不卡| 日韩三级视频一区二区三区| 天天一区二区日本电影三级 |