Printer Friendly

Computation of Aerodynamic Noise Radiated from Ducted Tail Rotor Using Boundary Element Method.

1. Introduction

With wide applications of helicopter, its noise level has become one of substantial indexes of helicopter performance and acceptability [1-5]. Helicopter tail rotor is a major noise source; therefore, investigation of tail rotor noise characteristic becomes distinguishably valuable and important to control the helicopter noise level. Tail rotor is an important aerodynamic assembly of conventional configuration helicopter, and its aerodynamic characteristics could have a significant impact on performance and flying qualities of the helicopter. For this reason the analysis of the tail rotor has been one of the research priorities in the field of helicopter technology [6-8].

Due to the effect of the duct, the slipstream pattern and aerodynamic characteristic of ducted tail rotor are different from the conventional open-type rotor. It can reduce the risk of component damage and enhance the operational safety. In addition, both the duct and rotor of the ducted tail rotor can produce thrust. Thrust which is produced by the duct is due to negative pressure on the duct inlet, which contributes maximum 50% to the total. Furthermore, the noise of ducted tail rotor is lower than the conventional open-type rotor [9,10]. Consequently, the analysis of the aerodynamic noise of the ducted tail rotor is the main topic for this paper.

In order to predict the aerodynamic noise, aerodynamic characteristics of the ducted tail rotor should be obtained firstly. Along with the development of the ducted tail rotor, many researchers have studied the aerodynamic characteristics of various types of ducted tail rotors [11-13]. In recent years, the main method which investigates the aerodynamic performance of ducted tail rotor is CFD (Computational Fluid Dynamic) technique. Among those CFD techniques, momentum source method is one of the commonly used methods. The momentum source with functional relationship to the local flow conditions represents the spanning rotor in flow field and is added to the N-S equations by UDF procedure. Cao and Yu [14] used it to solve the numerical simulation of turbulent flow around ducted tail rotor. Song et al. [15] applied it to analyze the aerodynamic characteristics of ducted tail rotor. Even though momentum source method is quite successful for the prediction of the aerodynamic performance of ducted tail rotor, detailed three-dimensional flow features cannot be accurately predicted. Because of the approximation of the blade shape parameters, the accuracy of the momentum source method is still limited. In addition, it is difficult to obtain the accurate noise source information by using the momentum source method. In this work, a CFD technique which is considered the blade shape parameters has been proposed to investigate the aerodynamic characteristics of ducted tail rotor, and the method is based on moving reference frame (MRF) model instead of momentum source method.

Calculation of the acoustic field radiated by rotating sources is a meaningful problem in the prediction of noise of aircraft rotors. The main techniques for noise prediction are based on acoustic analogy Curle equation, FW-H (Ffowcs Williams-Hawkings) equation, and the generalized treatment of Goldstein [16-18]. The inhomogeneous acoustic wave equation divides the aeroacoustic source into three types: monopole source, dipole source, and quadrupole source.

Although the aerodynamic characteristics of ducted tail rotor have been discussed by many researchers [19, 20], the investigations of the aerodynamic noise are relatively few. The prediction of rotating blade aerodynamic noise can be acquired easily by solving acoustic analogy Curle equation, but the scattering effect of the duct wall on the propagation of the sound wave is difficult to discuss and the scattering mechanism is rather complex. The acoustic boundary element method (BEM) is usually applied to predict the sound radiating and scattering field in the exterior and interior closed domain [21, 22]. Hu et al. [23] proposed a thin-body BEM, in which an imaginary interface is constructed to divide the domain into interior and exterior subdomains, and the imaginary surface is not discretized in the numerical implementation. In the present study, the thin-body boundary element method is also used to consider the scattering effect of the solid wall and to predict the far-field aerodynamic noise of the ducted tail rotor.

2. Governing Equations

Because the rotary speed is constant, in order to study the noise of the ducted tail rotor, we focus on the flow field in stable state. So the unsteady Navier-Stokes equation is applied here to calculate the flow field. In rotating coordinate system, the unsteady Navier-Stokes equation can be expressed as follows [24]:

[mathematical expression not reproducible], (1)

where V is control body, [partial derivative]V is the boundary area, and n is the normal vector. W is conservation vector, and W = [[[rho], [rho][u.sub.0], [rho][v.sub.0], [rho][w.sub.0], [rho]E].sup.T], p is the density of the fluid, [w.sub.0], [v.sub.0], [w.sub.0] are projection of absolute velocity in relative coordinate, and E is the internal energy of the unit fluid. [H.sub.w] and [H.sub.v] are the flux tensors. [S.sub.v] is source term.

Equation (1) is discretized by using Jameson cell-centered finite volume method. An explicit time integration algorithm based on the multistep Runge-Kutta is used to advance the solution in time. The linear system obtained by using the above method is solved by applying Gauss-Seidel method. The standard k - [epsilon] model is adopted here because of its reasonable accuracy and fewer resources. The moving reference frame is also used to solve the problem of rotating. The far-field boundary is treated as the nonreflecting boundary condition. In order to accelerate the convergence, a local time step and implicit residue smoothing have been used. The detailed procedure for solving the governing equations can be found in [25].

3. Free Field of Rotor Noise

Lighthill aeroacoustic analogy method is used to investigate the generation of the aerodynamic noise. The general equation is given by [16]

[mathematical expression not reproducible], (2)

where [c.sub.0] is the speed of sound and p' is the acoustic pressure. The terms on the right-hand side of (2) are interpreted as source terms. q denotes the volume source per unit volume per unit time, and [partial derivative]q/[partial derivative]t is the monopole source. [f.sub.i] denotes the force per unit volume, -[partial derivative][f.sub.i]/[partial derivative][x.sub.i] means the dipole source. [T.sub.ij] is the Lighthill stress tensor and [[partial derivative].sup.2][T.sub.ij]/[partial derivative][x.sub.i] [partial derivative][x.sub.j] denotes the quadrupole source.

After the CFD calculation, the aerodynamic noise at observation point is calculated using (2) by referring to the time history of the blade loading of the CFD results. Since a low rotational speed is considered for the tail rotor, the monopole source and dipole source are dominant, and quadrupole source could be negligible. A brief description will be given later to explain that the loading noise generated by dipole source is much bigger than the thickness noise generated by the monopole source for the tail rotor. Then (2) can be reduced to

[mathematical expression not reproducible]. (3)

The free space solution of (3) is [26]

[mathematical expression not reproducible], (4)

where G(x, y, t-[tau]) is Green function given by

[mathematical expression not reproducible], (5)

where r = [absolute value of x - y] denotes the distance between the observer location x and source location y. The relation between observe time t and source time [tau] is [tau] = t - r/[c.sub.0]. The square bracket means evaluation at the retarded time r.

Using the equality

[mathematical expression not reproducible], (6)

(4) can be written as

[mathematical expression not reproducible]. (7)

By using the Fourier transform and its inverse transform, force [f.sub.i], sound pressure p'(x, t), and G(x, y, t-[tau]) can be given in the frequency domain as

[mathematical expression not reproducible], (8)

[mathematical expression not reproducible], (9)

[mathematical expression not reproducible], (10)

where k = [omega]/[c.sub.0] is called the wave number and j = [square root of - 1] the complex unit. Substituting (10) into (7), we obtain

[mathematical expression not reproducible]. (11)

Comparing (11) with (9), we have

[mathematical expression not reproducible]. (12)

Differentiating (12) with regard to the direction of the normal unit vector n(x), we get

[mathematical expression not reproducible]. (13)

where P(y, [omega]) is the static pressure on the source surface in frequency domain. The function G(x, y, [omega]) and its derivations [partial derivative]G(x, y, [omega])/[partial derivative]n(y) and [[partial derivative].sup.2]G(x, y, [omega])/[partial derivative]n(x)[partial derivative]n(y) in above equations are expressed as follows:

[mathematical expression not reproducible]. (14)

4. Thin-Body Acoustic Boundary Element Method (BEM)

In this section, a thin-body boundary integral formulation is applied to calculate the far-field sound pressure. Due to the scattering effect of the solid wall in the duct, the total sound pressure is acquired as the sum of the incident and scattered pressure

P' (x, [omega]) = [P'.sub.I] (x, [omega]) + [P'.sub.S] (x, [omega]), (15)

where [P'.sub.I] and [P'.sub.S] are incident and scattered sound pressure in the frequency domain, respectively. The incident sound pressure can be obtained by using (12), and the scattered sound pressure will be got by using BEM. The calculated domain is shown in Figure 1; the surface of the duct is denoted by S and an imaginary surface s is added to close the duct and divide the computational domain into an interior subdomain [D.sup.+] and an exterior subdomain [D.sup.-]. The sound pressure on the outside of the surface S + s is denoted by [P'.sup.-] and that on the inside is denoted by [P'.sup.+]. The integral equation can be used to each subdomain [23] [mathematical expression not reproducible], (16)

where [n.sub.1] and [n.sub.2] are normal unit vectors at the two sides of the wall and [C.sup.+](x) and [C.sup.-](x) are the two constants that depend on the position of x

0,[mathematical expression not reproducible]. (17)

Adding (16) gives the thin-body boundary integral equation

[mathematical expression not reproducible], (18)

where [partial derivative]/[partial derivative][n.sub.1] = -[partial derivative]/[partial derivative][n.sub.2] = [partial derivative]/[partial derivative]n, and the continuous boundary conditions of the pressure and its partial derivation on the imaginary surface s are used

[mathematical expression not reproducible], (19)

x [member of] s.

The assumption of acoustic rigid boundary conditions are applied over the entire surface S

[mathematical expression not reproducible]. (20)

Then (18) reduce to

[mathematical expression not reproducible], (21)

[mathematical expression not reproducible], (22)

x [member of] [D.sup.+] [union] [D.sup.-] [union] s.

Equation (22) is not sufficient to obtain the two unknowns [P'.sub.+](x, [omega]) and [P'.sup.-](x, [omega]). Differentiating (21) with regard to the direction of normal vector n(x), it can be transformed into

[mathematical expression not reproducible]. (23)

The problem of scattering by the duct can be dealt with by initially solving (23) to calculate the sound pressure jump [P'.sup.+](y, [omega]) - [P'.sup.-](y, [omega]) on the surface S, and afterwards evaluating the acoustic pressure at any field point. Then the acoustic pressure on both sides of the duct could be easily got.

The integral Equations (21)-(23) cannot be directly calculated; a discretized scheme based on BEM should be used to calculate the unknown value [P'.sup.+](y, [omega]) - [P'.sup.-](y, [omega]). The simplest constant boundary element is applied in this paper and the thin-body surface Sis discretized into N elements. Each element has one node which is located in the center of the element. Then, (23) can be transformed into a system of algebraic equations.

[mathematical expression not reproducible], (24)

where [mathematical expression not reproducible] are the center of element i, j, respectively. d[S.sub.j] is the area of the element j. Therefore, the form of matrix of (24) is written as

[[[phi]].sup.T] = [B][[[phi]].sup.T]. (25)

Equation (25) may be solved easily by using the software Mathematica. When the unknown [phi] is calculated, the acoustic pressure at any field point can be obtained through (21) and (22).

5. Numerical Results and Discussions

5.1. Geometric Configuration and Mesh Generation. To demonstrate the effectiveness of the present method, we consider a model of ducted tail rotor at TsAGI, which has been installed on the Ka-60 helicopter. Figure 2 shows CATIA geometric model of TsAGI ducted tail rotor. Its geometric dimensions are tabulated in Table 1 as reported in [27].

A structured mesh is generated around the duct of the TsAGI model, and the surface quadrilateral for the all configuration is presented in Figure 3. The tail rotor of the TsAGI model is meshed with an unstructured mesh which is shown in Figure 4, and an adaptive encryption technique is adopted to mesh the leading edge and the trailing edge of the tail rotor.

5.2. Aerodynamic Performance of the Ducted Tail Rotor. An initial condition is made for the ducted rail rotor at a tip Mach number of 0.22 and a collective pitch angle of 40 degree at the blade root. At this moment, the predicted thrust from the present calculation is 9.547 kg. The rotor thrust in [28] and [29] are 9.54 kg and 9.69 kg, respectively. The coincidence indicates that our method can investigate aerodynamic performance of ducted tail rotor effectively. The predicted sectional thrust distribution along the rotor in Figure 5 is also compared with the momentum source theory and the vortex theory [27], which shows a good agreement.

In Figures 6 and 7, spanwise distribution of the induced downwash and the circumferential induced velocity are compared with the experiment and the results obtained by vortex theory [27]. It is also shown that the induced downwash and the circumferential induced velocity are generally in good coincidence with the experiment results. The results obtained by our method also compare well with those of the vortex theory. In addition, the present results indicate better comparison than those of the vortex theory when approaching toward the tip.

Figure 8 shows the distribution of the pressure on the surface of the ducted tail rotor in hover. It is obvious that the pressure changes violently on the inlet of the duct while few on the other zones remain the same.

A grid sensitivity analysis is given for a steady-state CFD simulation. The induced downwash, the circumferential induced velocity, and the induced attack angle are analyzed for the hybrid grids which are obtained by steadily increasing the meshing parameters. Their mean error ([e.sub.mean]) and maximum error ([e.sub.max]) with respect to the reference grid are established in Table 2. D denotes the max grid cell size.

From Table 2, we can see that the errors for all three flow variables are smaller and smaller when the grid cells increase. The small errors also display a converging trend.

To verify the effectiveness and the accuracy of the CFD method, the analytical solution of a dipole source has been investigated. A dipole has an axis, which is the line connecting the centers of the two monopoles (pulsating sphere with radius of a) in Figure 9. The two monopoles are positioned in the z direction at a distance of d = 2.5a from the center (d is the distance between the two monopole points). The pressure fluctuation induced by the dipole source for an observer positioned at (x, y, z) is

[mathematical expression not reproducible], (26)

where [mathematical expression not reproducible].

To perform the flow field characteristic of the dipole source, we put the dipole source into the center of a sphere (the radius [r.sub.s] of the sphere to be 3.25a). The speed of sound c0 is 340 m/s. The density for medium is 1.2 kg/[m.sup.3]. The angular velocity of the source is 1020 rad/s. The other parameters are A = A[phi][a.sup.2]U, a = 0.01 m, and U = 8m/s. The value of the spherical pressure at any point can be acquired by using the CFD method and (26), respectively.

The calculation is performed for this case, and the spherical pressure at any point is calculated by using CFD method and analytic method. The CFD results are compared with those of the analytical solution in Figure 10. From Figure 10, we can see that the CFD results are in perfect coincidence with the analytical solution.

5.3. Aerodynamic Noise of the Ducted Tail Rotor. The sound pressure is expressed as dB (decibels) and the predicted SPLs (sound pressure levels) are given by the following:

SPLs = 20 lg [P.sub.e]/[P.sub.r], (27)

where [P.sub.e] denotes predicted pressure and [P.sub.r] denotes the reference pressure and equals 2 x [10.sup.-5] Pa.

Before investigating the aerodynamic noise of the TsAGI ducted tail rotor, we firstly consider the radiation sound field generated by its isolate tail rotor. Figures 11 and 12 are the directivity of SPLs of the thickness noise and the loading noise, respectively. They can be obtained by using Farassat 1A formulation [26]. From the two figures, we can find that the thickness noise generated by the monopole source is much smaller than the loading noise generated by the dipole source for the isolate tail rotor. For the ducted tail rotor, the loading noise is dominant in rotation noise. The main reason is that the number of tail rotors is too high, and the size of tail rotors is too small. Therefore, the loading noise generated by the dipole source is considered only for the following discussions.

Through the analysis of aerodynamic performance of TsAGI ducted tail rotor, we can acquire the unsteady flow properties and the enough pressure data of the surface of the rotor. Then they are transformed into data in frequency domain by applying fast Fourier transform method. The farfield sound pressure is calculated using Curle's analogy in the frequency domain. Before using the BEM method to study the radiation and propagation of the sound source, we should use a simplified configuration which is shown in Figure 13 to consider the scattering effect of the duct. From Figure 13, we know that the surface is constructed by quadrilateral mesh. The maximum size of the mesh is 0.023 m and the number of elements is 2584. In this section, we use a point force method which regards the rotors as a series of points to investigate the scattering effect of the duct. The method can be applied in the prediction of the noise when the force fluctuations of the rotors are obtained.

In the present article, we just focus on rotation noise of the ducted tail rotor. The dominant noise in rotation noise is blade passing frequency (BPF) noise and its harmonic noise [30]. And the BpF of the ducted tail rotor is 2400 x 11/60 = 440 Hz.

Acoustic pressure in frequency domain is calculated by combining Curle equation and the thin-body BEM. For different observation distance, the directivity of SPLs at the BPF and 2BPF are shown in Figures 14 and 15, respectively. Free field is also rotor alone. The more interesting result is the comparison of the total field and the free field. This provides information about the scattering effect of the duct. For the angle (180 to 360 degrees), upstream of the rotor (180 to 270 degrees) is slightly louder due to scattering. And at the angle (270 to 360 degrees), it is quieter. From Figure 14(a), we can see clearly that the sound pressure at inlet is bigger than that of the outlet. From Figures 14(b), 14(c), and 15, we can also find the directivity of the incident sound pressure is symmetrical as the properties of the dipole sound source. These characteristics are consistent with the fact.

For 30 m/s forward flight velocity, the directivity of SPLs at the BPF is shown in Figure 16. The polar plot is made at a distance of ten times the rotor radius (10k) away from the center of the duct. Compared with the SPLs in Figure 14, the SPLs for forward flight are bigger than the SPLs for hover, and the noise directivity is also changed.

The directivity of SPLs at the BPF for rotational speed 2100 rpm and 2400 rpm are displayed in Figure 17. The SPLs of rotational speed 2400 rpm are bigger than that of 2100 rpm; this can be noticed by taking a closer look at Figure 17. The directivity of SPLs at the BPF for various observation distances is presented in Figure 18. From Figure 18, we can find that the SPLs become more and more small when the observation distance increases.

In order to illustrate the validation of the scheme, the results obtained in this paper are compared with the results obtained by commercial software Virtual.Lab Acoustics (VLA) [30]. The comparisons are shown in Figure 19.

In Figure 19, at the majority of azimuth of the ducted tail rotor, the SPLs obtained by our method are in good agreement with those obtained by using VLA.

6. Conclusions

The aim of this paper is to develop the Curle/thin-body BEM method for predicting the aerodynamic noise of the ducted tail rotor numerically. A CFD method based on MRF model is used to study the aerodynamic performance of ducted tail rotor. Calculations are made for TsAGI ducted tail rotor, and the results are compared with the known results and the analytical solution of a dipole source for validation.

The unsteady blade force is also obtained and transferred to the Curle equation to calculate the sound source without the duct. In order to consider the scattering effect of the duct, a thin-body BEM is proposed to predict the far-field aerodynamic noise. The calculated SPLs at BPF accord well with the results obtained by VLA and show that our method can effectively predict the far-field SPLs of the ducted tail rotor. The numerical results also show that the scattering effect of the duct should be included in the prediction of the noise. The control of the inlet noise is advisable for the ducted tail rotor noise reduction.
Nomenclature

V:                                   Control body

[partial derivative]V:               Boundary area

n:                                   Normal vector

W:                                   Conservation vector

[rho]:                               Fluid density

[u.sub.0], [v.sub.0], [w.sub.0]:     Projection of absolute velocity
                                     in relative coordinate

E:                                   Internal energy of the unit fluid

[H.sub.w], [H.sub.v]:                Flux tensors

[S.sub.v]:                           Source term

[omega]:                             Rotational angular velocity

[c.sub.0]:                           Sound speed

G(x, y, t - [tau]):                  Green function in time domain

G(x, y, [omega]):                    Green function in frequency
                                     domain

[f.sub.i]:                           Unsteady force

[F.sub.i]:                           Unsteady force in frequency

t:                                   Observe time

[tau]:                               Source time

P':                                  Sound pressure in time domain

P':                                  Sound pressure in frequency
                                     domain

[P'.sub.I]:                          Incident sound pressure in
                                     Frequency domain

[P'.sub.S]:                          Sound pressure in frequency
                                     domain

x:                                   Observer position

y:                                   Source position

r:                                   Distance between observer
                                     and source

k:                                   Wave number

S:                                   Duct surface

s:                                   Imaginary surface

n(x):                                Normal unit vector on the
                                     observer surface

n(y):                                Normal unit vector on the
                                     source surface

[P.sub.e]:                           Predicted pressure

[P.sup.r]:                           Reference pressure.


http://dx.doi.org/10.1155/2017/5726253

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work was supported by the National Natural Science Foundation of China under Grant no. 51307004.

References

[1] P. Cyril, Z. Joelle, R. Olivier, and K. Martin, "Helicopter rotor noise prediction using ONERA and DLR Euler/Kirchhoff methods," Journal of the American Helicopter Society, vol. 44, no. 2, pp. 121-131, 1999.

[2] E. Greenwood and F. H. Schmitz, "Effects of ambient conditions on helicopter rotor source noise modeling," Journal of Aircraft, vol. 51, no. 1, pp. 90-103, 2014.

[3] M. J. Kingan, "Advanced open rotor noise prediction," Aeronautical Journal, vol. 118, no. 1208, pp. 1125-1135, 2014.

[4] W. Johnson, Helicopter Theory, Princeton University Press, 1980.

[5] Y. Mao, Q. Zhang, C. Xu, and D. Qi, "Two types of frequency-domain acoustic-velocity formulations for rotating thickness and loading sources," AIAA Journal, vol. 53, no. 3, pp. 713-722, 2015.

[6] R. Mouille, "Fenestron, shrouded tail rotor of the SA341 gazelle," Journal of the American Helicopter Society, vol. 15, no. 4, pp. 3137, 1970.

[7] R. Mouille and F. D. Ambra, "The fenestron a shrouded tail rotor concept for helicopters," in Proceedings of the 38th Annual Forum ofAHS, June 1986.

[8] M. Lowson, "Acoustic kinematics of helicopter BVI noise," International Journal of Aeroacoustics, vol. 14, no. 1-2, pp. 229256, 2015.

[9] C. Keys, M. Sheffler, S. Weiner et al., "LH wind tunnel testing: key to advanced aerodynamic design," in Proceedings of the 47th annual forum of American Helicopter Society, pp. 77-87, Phoenix, Ariz, USA, 1991.

[10] R. A. Desjardins, W. L. Noehren, and A. N. Bertolazzi, "Design and flight test evaluation of the Fantail antitorque system," in Proceedings of the 47th Annual Forum of American Helicopter Society, pp. 857-867, Phoenix, Ariz, USA, May 1991.

[11] R. G. Rajagopalan and C. N. Keys, "Detailed aerodynamic analysis of the RAH-66 FANTAIL using CFD," Journal of the American Helicopter Society, vol. 42, no. 4, pp. 310-320, 1997

[12] S. Bandoh, M. Fudamoto, and T. Akiyama, "The ducted tail rotor system of the New Observation (XOH-1)," in Proceedings of the AHS International Meeting on Advanced Rotorcraft Technology and Disaster Relief, Gifu, Japan, 1998.

[13] A. Vuillet and F. Morelli, "New aerodynamic design of the Fenestron for improved performance," in Proceedings of the 12th European Rotorcraft Forum, pp. 57-69, 1986.

[14] Y. Cao and Z. Yu, "Numerical simulation of turbulent flow around helicopter ducted tail rotor," Aerospace Science and Technology, vol. 9, no. 4, pp. 300-306, 2005.

[15] C. H. Song, Y. F. Lin, W. X. Chen, and Y. T. Liu, "CFD analysis for the ducted tail rotor based on momentum source method," Helicopter Technique, vol. 1, pp. 6-11, 2009.

[16] J. E. F. Williams and D. L. Hawkings, "Sound generation by turbulence and surfaces in arbitrary motion," Philosophical Transactions of the Royal Society, vol. 264, no. 1151, pp. 321-342, 1969.

[17] F. Farassat and K. S. Brentner, "Uses and abuses of the acoustic analogy in helicopter rotor noise prediction," Journal of the American Helicopter Society, vol. 33, no. 1, pp. 29-36, 1988.

[18] F. Farassat and K. S. Brentner, "Supersonic quadrupole noise theory for high-speed helicopter rotors," Journal of Sound and Vibration, vol. 218, no. 3, pp. 481-500, 1998.

[19] Y. Cao, K. Chen, and J. Wang, "Analyses on ducted tail rotor and airfoil aerodynamic characteristics with CFD," Aircraft Engineering and Aerospace Technology, vol. 77, no. 1, pp. 62-67, 2005.

[20] Z.-W. Yu and Y.-H. Cao, "CFD simulation and validation of ducted tail rotor," Journal of Aerospace Power, vol. 21, no. 1, pp. 19-24, 2006.

[21] R. D. Ciskowski and C. A. Brebbia, Boundary Element Method in Acoustics, Elsevier Science Publishing, Norw York, NY, USA, 1991.

[22] T. W. Wu, Boundary Element Acoustics: Fundamentals and Computer Codes, WIT Press, Southampton, UK, 2000.

[23] B.-B. Hu, H. Ouyang, Y.-D. Wu, G.-Y. Jin, X.-Q. Qiang, and Z.H. Du, "Numerical prediction of the interaction noise radiated from an axial fan," Applied Acoustics, vol. 74, no. 4, pp. 544-552, 2013.

[24] X.-H. Li and J.-J. Xiong, "Navier-Stokes equation-based aerodynamic interaction between rotor blades in hover phase," Acta Aerodynamica Sinica, vol. 25, no. 2, pp. 271-276, 2007

[25] Y. Cao, J. Wang, and Y. Su, "Mixed Jameson/total-variation-diminishing scheme applied to simulating rotor airfoil flow field," Journal of Aircraft, vol. 40, no. 1, pp. 213-216, 2003.

[26] F. Farassat, "Derivation of formulations 1 and 1A of Farassat," Tech. Rep. NASA/TM-2007-214853, L-19318, NASA, 2007

[27] B. N. Bourtsev and S. V. Selemenev, "Fan-in-Fin performance at hover computational method," in Proceedings of the 26th European Rotorcraft Forum, The Hague, The Netherlands, September 2000.

[28] W. R. Tan, S. L. Huang, and G. H. Xu, "Analysis on aerodynamics characteristics of helicopter ducted tail rotor, the 3th science and technology society of Chinese aviation," pp. 204-210, 2007.

[29] H. D. Lee, O. J. Kown, and J. Joo, "Aerodynamic performance analysis of a helicopter shrouded tail rotor using an unstructured mesh flow solver," in Proceedings of the 5th Asian Computational Fluid Dynamics Forum, pp. 73-85, 2003.

[30] K. Kucukcoskun, "Prediction of free and scattered acoustic fields of low-speed fans," 2012.

Yunpeng Ma, Mingxu Yi, Lifeng Wang, and Jun Huang

School of Aeronautic Science and Engineering, Beihang University, Beijing 100191, China Correspondence should be addressed to Mingxu Yi; yimingxu88@163.com

Received 23 November 2015; Revised 26 April 2016; Accepted 25 September 2016; Published 3 January 2017 Academic Editor: Yuri Vladimirovich Mikhlin

Caption: Figure 1: A diagram of acoustic scattering by a thin-body.

Caption: Figure 2: Geometric model of TsAGI ducted tail rotor.

Caption: Figure 3: Surface quadrilateral for the duct.

Caption: Figure 4: Surface quadrilateral for the tail rotor.

Caption: Figure 5: Comparison of the spanwise distribution of the blade loading.

Caption: Figure 6: Comparison of the spanwise distribution of the induced downwash.

Caption: Figure 7: Comparison of the spanwise distribution of the circumferential induced velocity.

Caption: Figure 8: Pressure distribution for the surface of the ducted tail rotor in hover.

Caption: Figure 9: Dipole source and spherical are shown.

Caption: Figure 10: The spherical pressure at any point for different method.

Caption: Figure 11: Directivity of SPLs of the thickness noise at different distance.

Caption: Caption: Figure 12: Directivity of SPLs of the loading noise at different distance.

Caption: Figure 13: Mesh of duct used in BEM.

Caption: Figure 14: Noise directivity pattern at the BPF: (a) R =1 m; (b) R = 5 m; (c) R = 3 m.

Caption: Figure 15: Noise directivity pattern at the 2BPF: R = 3 m.

Caption: Figure 16: Noise directivity pattern at the BPF: 30m/s forward flight.

Caption: Figure 17: Directivity of SPLs of the ducted tail rotor at different rotational speed.

Caption: Figure 18: Directivity of SPLs of the ducted tail rotor at different distance.

Caption: Figure 19: The comparison of the SPLs at different distance: (a) R =1 m; (b) R =3 m.
Table 1: Geometric dimensions of the ducted tail rotor.

             Radius                    297 mm
Rotor   Number of blades                 11
            Solidity                   0.4951
           Twist angle                 -12 deg
            Tip speed                 74.6 m/s
             Airfoil         NACA23012 (r/R = 0.35-1.0)
Duct    Tip clearance gap               0.01R
          Inlet radius                  0.2R
         Diffuser length                0.7R
         Diffuser angle               4 degrees

Table 2: Mean and maximum values of the induced downwash, velocity,
and attack angle for different grid densities.

Meshing parameters            Induced downwash

Grid cells     D(mm)          [e.sub.mean]   [e.sub.max]

48 k           1.335          1.028          3.374
96 k           0.861          0.914          3.021
258 k          0.742          0.857          2.853
480 k          0.604          0.777          2.634
712 k          0.359          0.513          2.028
837 k          0.147          0.434          1.634
968 k          0.082          0.125          1.008

Circumferential induced       Induced attack angle
velocity

[e.sub.mean]   [e.sub.max]    [e.sub.mean]   [e.sub.max]

0.799          1.376          0.542          0.835
0.723          1.224          0.407          0.801
0.583          1.011          0.380          0.788
0.371          0.969          0.236          0.728
0.126          0.808          0.105          0.519
0.101          0.761          0.095          0.353
0.098          0.486          0.088          0.104
COPYRIGHT 2017 Hindawi Limited
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2017 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Title Annotation:Research Article
Author:Ma, Yunpeng; Yi, Mingxu; Wang, Lifeng; Huang, Jun
Publication:Mathematical Problems in Engineering
Article Type:Report
Date:Jan 1, 2017
Words:5261
Previous Article:Design, Analysis, and Experiment for Rescue Robot with Wheel-Legged Structure.
Next Article:Novel Stability Analysis for Uncertain Neutral-Type Lur'e Systems with Time-Varying Delays Using New Inequality.
Topics:

Terms of use | Privacy policy | Copyright © 2020 Farlex, Inc. | Feedback | For webmasters