Effect of Upstream Side Flow of Wind Turbine on Aerodynamic Noise: Simulation Using Open-Loop Vibration in the Rod in Rod-Airfoil Conﬁguration

: Adaptive and ﬂexible control techniques have recently been examined as methods of controlling ﬂow and reducing the potential noise in vertical axis wind turbines. Two-Dimensional (2D) ﬂuid ﬂow simulation around rod-airfoil is addressed in this study as a simple component of the wind turbine by using Unsteady Reynolds Averaged Navier–Stokes (URANS) equations for prediction of noise using Ffowcs Williams-Hawkings (FW-H) analogy. To control the ﬂow and reduce noise, the active controlling vibration rod method is utilized with a maximum displacement ranging from 0.01 C to 1 C (C: airfoil chord). Acoustic assessment indicates that the leading edge of the blade produces noise, that by applying vibration in cylinder, blade noise in 0.1 C and 1 C decreases by 22 dB and 35 dB, respectively. Applying vibration is aerodynamically helpful since it reduces the ﬂuctuations in the airfoil lift force by approximately 48% and those in the rod by about 46%. Strouhal assessment (frequency) shows that application of control is accompanied by 20% increase. Applying vibration in the rod reduces the ﬂow ﬂuctuations around the blade, thus reduces the wind turbine blade noise. This idea, as a simple example, can be used to study the incoming ﬂow to turbines and their blades that are affected by the upstream ﬂow.


Introduction
Loss of fossil energy sources and global warming has led to the little economic acknowledgment of carbon as an international strategy for sustainable development [1]. Therefore, there has been an increase in the demand for green and renewable energies. Among different types of renewable energies, wind has had incredible growth all over the world and plays a significant role in the future economy [2]. As a sustainable energy source, wind energy can be an important factor in increasing energy production in different countries and international politics. With the worldwide spread of wind plants and wind farms, it would be necessary to investigate different aspects of wind energy, especially the environmental aspect, along with finding appropriate ways to increase its efficiency [3]. The aerodynamic noise generated by wind turbines, as one of the main sources of noise pollution, has always seriously challenged engineers, because the generated noise has reduced the quality of life near residential areas.
The effect of noise produced by rotating blades is perceivable also by citizens living in the area [4,5]. So, identifying aspects such as the factors involved in generation or reduction of noise and the appearance are required along with examining the flow structure to meet public interests. This has resulted in a researcher's attempt to study and understand aerodynamic noise generated by vertical and horizontal axis turbines. Oerlemans et al. [6] conducted an experimental study to predict noise in a three-bladed horizontal axis turbine. The results indicated that beside a minor source at the rotor hub, all noise is practically produced during the downward movement of the blades. Their experiment also showed that the main source of the noise was broadband trailing edge noise. In a numerical investigation around vertical axis turbines, Botha et al. [7] found that increasing the velocity entering the turbine and blades increases the noise generated by the blades. Maizi et al. [8] used Unsteady Reynolds Averaged Navier-Stokes (URANS) and direct numerical simulation (DNS) models to conduct a numerical and Three-Dimensional (3D) study of horizontal axis turbine blades. Among their results, 5 dB aerodynamic noise reduction after varying the shape of blades be stated. Ghasemian and Nejat [9] employed improved delayed detached eddy simulation (IDDES) turbulence (hybrid RANS-large eddy simulation (LES) (k-ω-Shear-Stress Transport (SST)_LES) model to conduct a 3D study of flow inlet velocity into the horizontal axis turbine. The results indicated that increasing the velocity from 7 m/s to 10 m/s increases the aerodynamic noise by 8 dB. Tadmasa and Zangeneh [10] studied the effect of wind speed and blade rotation on the horizontal axis blade noise by using k-ω-SST turbulence model and Ffowcs Williams-Hawkings (FW-H) analogy. The results showed that increase in the rotational speed increases the aerodynamic noise. Cho et al. [11] examined NREL model turbine acoustics by employing a receiver array. They identified the blade tip as the cause of noise. They also found that the maximum recorded noise occurs in stall. Mo and Lee [12] investigated the aerodynamic noise of the blades in a rotating turbine by using LES turbulence model. They modeled the aerodynamic noise of the field by using FW-H analogy for frequencies lower than 500 Hz. They identified the vortices created around the blade as a result of the flow passing the trailing edge of the adjacent blades as the cause of noise generation. Mohammad [13] conducted a Two-Dimensional (2D) study of the effect of the blade cross-section and the rotational speed on noise by using URANS equations. His results indicated that highvelocity rotation entails more noise than normal conditions. He also examined the three double-airfoil acoustic turbine afterwards and achieved noise reduction by increasing the distance between the two airfoils [14].
Noise generated by the turbulent wake and the blade's trailing edge (airfoil) is an important source of noise in turbomachinery. In case the airfoil is under the influence of a turbulent flow, the noise on the leading edge may be dominant [15]. The type of inlet flow to the turbine blade affects its performance and aerodynamic characteristics, so the turbines used downstream of other bluff bodies are exposed to turbulent flows, and the behavior of the flow around them changes and may also have destructive effects [16]. To adapt such a phenomenon, a combination of rod-airfoil has been employed for investigation of noise and flow turbulence structure, namely, the flow passing the rod at the upstream side penetrates towards the airfoil at the downstream side of the rod. The distance between the rod and the airfoil entails alterations in the flow structure and noise and intensifies the fluctuations. Thus, investigation of rod-airfoil noise and controlling it has recently attracted the researchers' attention. Munekata et al. [17] examined the effect of altering the distance between the rod and the airfoil on aerodynamic aspects experimentally. They found that von Kármán vortex street does not develop by reducing the distance. Furthermore, the noise resulting from the interaction of the flow with the rod and the airfoil is reduced. Moreover, in another study, they found that increasing the angle of attack reduce the noise generated due to the structure of the diffused wake [18]. Jacob et al. [19] obtained the noise spectrum around the airfoil and experimentally found that the interaction of the object and the flow largely accounts for noise radiation. Li et al. [20] accomplished noise reduction in an experimental study by controlling the rod-airfoil aerodynamic noise by simultaneous application of blowing behind the rod and covering the airfoil leading edge with a thick cover. The Rod-Airfoil combination has also been numerically investigated. For instance, Magagnato et al. [21], Boudet et al. [22], Greschner et al. [23], Agrawal et al. [24], Giret et al. [25], and Jacob et al. [19] measured the rod-airfoil noise with LES model numerical method. Their results had a relatively precise correspondence with experimental data. Casalino et al. [26] and Jacob et al. [19] examined noise using 2D k-ω-SST turbulence model and Caraeni et al. [27], Gerolymos et al. [28], also investigated the rod-airfoil noise by employing Detached eddy simulation (DES) model. Jiang et al. [29] conducted a study on noise and flow fluctuations by altering the distance between the rod and the airfoil. They found that if this distance is twice the rod diameter, the pressure fluctuations significantly decrease. Chen et al. [30] achieved noise reduction in a 3D numerical study by wavy leading edges the airfoil leading edge. Rousoulis et al. [31] conducted a numerical 3D and k-ω-SST turbulence model study by rotating the rod with different frequencies to examine noise. They found that with rod rotation with a frequency higher than the natural shedding frequency, the noise is reduced. They also examined the results of 2D model, which showed that 2D results were also reliable [32]. Abbasi and Suori [33] achieved noise reduction by applying blowing and suction in rods and airfoils with a 2D URANS model. Li et al. [34] used natural rod-base blowing to reduce aerodynamic noise. By changing the diameter of the internal slot, they reduced the noise to 12 dB at best.
Researchers have always tried to suggest wind turbines more to the environmentallyfriendly than ever before. It is possible to increase the efficiency of wind turbines by controlling the inlet flow. Furthermore, controlling the inlet flow helps to control the aerodynamic noise generated by the flow. There have been great efforts in the aerodynamic noise of vertical-centered wind turbines, but largely it is considered that the turbine is not located downstream of a bluff body; so, the noise from the blade trailing edge being known as the dominant noise source. However, if the blade is located downstream of a bluff body, the aeroacoustic condition is affected, so in this study, the rod-airfoil configuration is used to apply this phenomenon. A review of the literature indicates that employing the vibrating rod active control method has not been investigated so far for controlling the flow and noise. In fact, other researchers have not used this configuration as a simple component of a wind turbine to study the flow around the blade. Hence, the present study investigates the effect of the control method on Rod-Airfoil noise by applying FW-H equations. In addition, the details of the 2D flow structure have been addressed with the k-ω-SST turbulence model to achieve a deeper understanding of the process of improving the flow and aerodynamic coefficients.

Material and Methods
Simulation of the results was achieved using a two-equation URANS turbulence model. The equations governing the turbulence model must be explained and solved. Moreover, FW-H equations have been utilized to predict the noise generated by the incompressible flow, and the equations of this analogy need to be explained and solved like the turbulence model. Reliable results are predictable with an appropriate solution domain and boundary conditions and the grid must not have radical alterations. It should be noted that in the next equations, Einstein's summation rules are employed. In this notation, the indices repeating twice in a term are summed over the range of the index. In this model, continuity [36][37][38][39] and momentum equations are stated as: S M denotes the sum of body forces, µ e f f is the effective viscosity, and P " is defined as: where µ e f f = µ + µ t . It is assumed in this model that µ t depends on the flow kinetic energy and energy loss according to the following relation.
This model does not have complicated nonlinear equilibrium functions in the k-ε model. k and ω are directly obtained from the differential equations of turbulent kinetic energy transfer and turbulent frequency transfer (Equations (5) and (7)) [40,41]: where Here the general description for each of the terms in Equations (5) are the usual terms for accumulation, convection, diffusion, production, and dissipation of ω. Last term (VI) is called a "cross-diffusion" term, an additional source term which has a role in the transition of the modeling from ε to ω [40]. γ, σ k and β are constants.
Here parameters can be defined as ( In 1950 s, Lighthill [42] extracted an accurate form of the heterogeneous wave equation based on Navier-Stokes and continuity equations. However, the extracted equation is limited since it estimates the sound waves only at unbounded domains and is incapable of considering reflection effects, absorption wave failure and extended wave dispersion over solid boundaries. Later, Curle [43] developed Lighthill's model to include the presence of a solid object. Williams and Hawkings [44] improved the previous theories (Lighthill [42] and Curle [43]).
To analyze the noise generated by the flow, FW-H equation (Equation (8)) is used in addition to solving Navier-Stokes equations. This analogy includes all sound sources (monopole, dipole, and quadrupole) which is shown in Equation (8). This equation is basically a heterogeneous wave equation and is derived from the continuity equation and Navier-Stokes equations. Equation (9) is Lighthill stress tensor.
where c 0 denotes the speed of sound in a far field. H( f ) is Heaviside function and δ( f ) denotes the Dirac Delta function. u n and v n are the fluid velocity in the direction normal to the integration surface and the surface velocity in the direction normal to the integration surface, respectively. Furthermore, τ ij , δ ij , ρ and P are respectively the viscous stress, Kronecker delta, the density and the sound pressure in the far field in P = p − p 0 . In FW-H equation, f = 0 indicates the source surface (cylinder and splitter plate), f > 0 is unbounded space, n i represents the vertical unit vector towards outside the area, and P ij is the compressive stress tensor, which is defined as Equation (10) for a Stokes flow.
The first term of the right of Equation (8) refers to quadrupole sources that is caused by unsteady shear stresses, and these are outside the source surface and have the least contribution in each period of sound generation, so it can be neglected. The quadrupole term is the distribution of the volume of sources defined by the Heaviside function H( f ). The second term, known as the dipole or loading source, is related to the forces applied to the surface of the body. The last term of the right of Equation (8) is for monopole sound source, also known as the thickness noise, caused by fluid displacement by the body surface in the flow field. In Equation (8), δ( f ) represents the mono-and dipole sources. To solve the FW-H equation, the free-space Green's function is utilized to compute the sound pressure at the receiver location x. Acoustic pressure is defined by Equation (11), where the Q, L, and T components represent the quadrupole, loading (dipole), and thickness (monopole), respectively. These parameters are described in Equations (12)- (14), respectively [44,45]. Since quadrupole sources have not an effective role in noise generation, they are ignored.
Here the subscript (r et ) means for the retarded time, and M is local Mach number. The radiated sound, based on unsteady force on the body in the calculations, is written as [46]: where r is the observer distance from the sound source and ∂F i ∂t is the time gradient of unsteady force on the body.

Geometric Model
In the present study, the geometric model is adopted from the experimental model proposed by Jacob et al. [19]. Thus, an airfoil is of NACA 0012 type with a 0.1 m (C = 0.1 m) chord is located at the downstream side of a rod at a distance equal to the airfoil cord (L = 1 C) with a diameter of 0.1 times the airfoil chord (D = 0.1 C). The Reynolds number is 48,000 based on the rod diameter, and 480,000 based on the airfoil chord, and the Mach number is 0.2. Furthermore, between the rod and the airfoil in the vertical direction (y), there is a y = 2 mm shift between the rod axis and the airfoil chord. The x-axis is parallel to the flow, and the y-axis is perpendicular to the flow in downward and upward directions. The origin of coordinates is also located at the leading edge of the airfoil. A literature review shows that the rectangular domain with the intended dimensions is suitable for solving this problem ( Figure 1) [30,31]. To measure the aerodynamic noise of the flow, a microphone is embedded at the top and center of the airfoil at a distance of 18.5 C. Figure 1 illustrates a schematic of computational domain along with the boundary conditions.

Meshing Process
Considering the fact that this paper aims at investigation of the noise generated through the transient fluid flow, a precise estimation of flow turbulence is of high importance. To accurately estimate extreme velocity gradients, a precise meshing is deemed necessary and vital especially near the wall. Given that the rod has vibrational movements in the vertical direction to the flow, meshing of the calculation field is depicted as triangles in Figure 2a. The unstructured mesh has a lower computational cost than structured mesh. However, they have almost the same accuracy. Moreover, since the recorded noise is generated from variations of pressure on a solid surface, the meshing around the rod and the airfoil must be finer, and to avoid sudden variations in mesh density, the meshing of the field around the rod and the airfoil is quadrilaterals as in Figure Figure 2b,c. The circumference of rod and airfoil is discretized by 200 and 300 nodes, respectively. To assess the quality of generated mesh, a validation study is presented in Table 2. In the three examined density levels, y + is investigated to determine the best mesh for simulation because it plays an important role in estimating the flow behavior in the boundary layer. y + denotes the non-dimensionalized distance of the grid points from the wall and defined as y + = u τ y /v (u τ : friction velocity and y: distance from the wall).  Rousoulis et al. [29] y + < 3

Solver Settings
ANSYS-Fluent commercial software has been used for the numerical analysis of the flow. URANS have been used for solving the governing equations. Two-equation turbulence models are the simplest forms of turbulence models. Two-equation turbulence models allow the determination of both the turbulent velocity and length scales by solving two separate transfer equations [47]. The k-ω-SST turbulence model has been used to accurately estimate viscosity eddies. Stability, economy, and appropriate precision for a wide range of turbulent flows explain its popularity in the simulation of turbulent heat transfer and flow in the industry [10,31,40,45,[48][49][50][51]. To provide for a comparison of the present numerical solution results with the empirical research conducted by Jacob et al. [19], all the applied dimensions, parameters, and boundary conditions have been considered exactly according to the conditions of the empirical research. A constant flow at the speed of 72 m/s has been assumed at the inlet boundary. For solution, simple pressure-velocity coupling with Least Squares Cell-Based Gradient discretization was used. Furthermore, second-order spatial discretization approach was used for pressure and momentum. Moreover, first-order spatial was employed for turbulent kinetic energy and dissipation rate. For the unsteady simulation process, a time step of 0.0001 s is used and the total simulation time is 0.8 s with 50 inner iterations at each time interval. After reaching the aerodynamic force and convergence diagram to an iterative pattern and ensuring its stability (0.35 s) (Figure 3), the FW-H analogy is coupled with a flow solution, and the data are extracted at the end of the analysis (0.8 s). The convergence condition in every governing equation is for residuals to reach 10 −5 . Moving mesh is implemented for the rigid body (rod) and near region, whereas the outer of near region is kept stationary (deforming region) ( Figure 1). As previously mentioned, FW-H solution has been used in this study in order to extract the noise pressure level and the aerodynamic noise. Moreover, a frequency resolution, equal to the sampling frequency (i.e., 10 kHz), is employed for the FW-H results. Since the literature review has shown that 2D models are also reliable, this model is used in this study to extract noise and to simulate flow structure [26,31,33,45,52]. The numerical procedure outlined above for predicting the generation and propagation of aeolian tone noise is summarized in Figure 4 [53].

Results and Discussion
First, the results' becoming independent of the grid is addressed. Then validation of acoustics and aerodynamic parameters is made to examine the validity of the current study. Eventually, acoustic control and flow disturbances are addressed by using the rod vibration.

Examination of Independence from the Grid
To make sure of the conducted meshing and the number of calculation cells considered, dependence from the grid is made for the results of the numerical analysis. To this aim, the meshing of the measurement zone has been performed in multiple cases around the Rod-Airfoil, the sound pressure level (SPL), and the peak amplitude Strouhal number, and the mean pressure parameter around the rod and airfoil and mean velocity parameter are obtained as the most important results. Meshing has been carried out with different numbers of computational cells. In Table 3, the maximum value of SPL (i.e., the peak value of the noise spectrum (NSP) recorded by embedded microphone) is calculated according to the contribution of noise generated by the rod and airfoil in the total of NSP for different numbers of calculation cells. Strouhal number that occurs in the NSP is also recorded. Equation (16) presents the calculation of SPL and Equation (17) demonstrates the calculation of the Strouhal number. In Equation (16), P rms denotes the root mean square of pressure disturbances, and P re f is equal to 2 × 10 −5 (pa). It is evident that following the variations in the meshing of the solution field from case I to case III, the NSP and the peak amplitude Strouhal number show significant changes; however, by increasing the number of cells from case III to case IV, the results do not undergo much change. Maximum SPL recorded from the Rod-Airfoil configuration is equal to the dominant sound source. It is clear that airfoil is the predominant source.  Figure 5 depicts the results of mesh study for the average pressure coefficient around the rod and the airfoil and velocity at x/C = − 0.25 and x/C = + 0.25. Clearly, in case (I) and (II), the pressure coefficient and velocity are of relatively greater value; however, with a finer meshing in case (III), this value is decreased. By increasing the number of trihedral and tetrahedron grids in case (IV), no significant change occurs in the results in comparison with case (III). According to the results of Table 2 and Figure 5, it is clear that the third and fourth cases are slightly different in predicting the results, which can be ignored. On the other hand, fining the mesh lead to increasing the computational cost. Therefore, considering the acoustic and the aerodynamic results and reduce the computational cost, the grid in case (III) is appropriate for the present investigation.

Validation of Results
Penetration of vortex structures at the airfoil leading edge results in the generation of powerful acoustic sources such as the noise generation mechanism in turbomachinery in a study by Homicz and George [15]. Hence, the leading edge significantly has the most contribution to the generation of noise in comparison with the rod and the trailing edge. Validation of the acoustic results are achieved using the SPL and the Strouhal number ( Figure 6) in comparison with experimental results of Jacob et al. [19]. The peak amplitude Strouhal number, equal to 0.2120, is stated as Equation (17) at a frequency of 1526 Hz. The obtained value has a trivial difference with that of the empirical study, which is equal to 0.192. Furthermore, the obtained NSP in the numerical method is 99 dB, the value which is 6.5 dB higher than that in the results of Jacob et al. [19]. The discrepancy is due to the nature of the present modeling. 2D models do not take into account the features of experimental and 3D models, so they make mistakes in correctly estimating the separation position, which leads to more estimates of the amount of lift force and its fluctuations [54,55]. On the other hand, due to the significant contribution of lift in noise production, this leads to a higher estimate of SPLs than experimental and 3D studies. Another feature of 2D models is that it does not take into account the spanwise effect, which leads to a higher estimate of the frequency (St) at NSP, according to what is obtained in Figure 6. The difference in SPL between the present model and the experimental model, in other Strouhal numbers, is due to the URANS method, which does not consider smaller eddies in the near wake region [33,52,56,57]. Although the flow of Rod-Airfoil in high Reynolds numbers is inherently 3D, 2D approaches can still illustrate the key aspects of vortex shedding with significant cost reductions. Furthermore, despite the differences and shortcomings of this model, its prediction capability is high, and the simulated SPL is reasonably consistent with the experimental results. Evaluation of the average velocity is an important parameter for examining the flow pattern. Consequently, this study attempts to validate this parameter as one of the most important aerodynamic parameters. The average velocity validation has been achieved at two points in the flow: X C = −0.25 and X C = +0.25. X C = +0.25 and X C = −0.25 have been devised as perpendicular to the airfoil and at the rod wake respectively. In Figure 7a, the minimum value of u mean U 0 · X C = −0.25 in the experiential study equals 0.701, and its minimum value in the numerical study is 0.721; the relative error is 1.6%. As shown in Figure 7b, the average velocity difference at X C = +0.25 between the empirical and the numerical methods is 5.9%. It is evident that the present numerical method has a trivial difference with LES, and both patterns are the same as that of the empirical study with very little difference. After reviewing the results and the acceptable validation provided, Computational Fluid Dynamics (CFD) prediction is possible for further sufficiently accurate investigation of the flow pattern for the present study.

Applying Control and Study of Its Effect on Aerodynamic Noise
According to Homicz and George study [15], it was found that the leading edge is the most important source of noise generation. Therefore, in the present study, the effect of the open-loop vibrating rod control method is investigated. In other words, the rod will have an open-loop sinusoidal motion perpendicular to the flow whose range varies from 0.02C to 1C. The frequency of harmonic motion is 5 Hz.
The mesh motion is shown in Figure 8 by moving the rod, while A = 1C is considered as the displacement.
In the following, the role of the introduced control method on aeroacoustic characteristics is examined. Figure 9 illustrates the SPL variations by applying vibration to the rod and increasing the amount of rod displacement. As is observed, applying vibration has resulted in a drop in the aerodynamic noise. Since the dominant noise in the Rod-Airfoil combination is due to flow disturbance at the downstream side of the rod on the airfoil leading edge, the rod motion perpendicular to the flow causes the flow at the downstream side of the rod to also distance away from the airfoil leading edge and lead to flow disturbance reduction around the edge. Eventually, the pressure fluctuations at the airfoil attack surface also decrease. Hence, it is expected that noise also decreases with rod motion perpendicular to the flow. The extracted charts demonstrate the dominance of the airfoil noise. Clearly, the recorded airfoil noise equals that of Rod-Airfoil and the rod noise without applying control is 10 dB less than the Rod-Airfoil noise. Applying oscillation to the cylinder in the range of A ≤ 0.1C generally reduces the SPL, caused by the Rod-Airfoil, airfoil and rod; by about 22 dB (Figure 10). By applying oscillation to the rod, the vortex shedding phenomena are changing which leads to an increase in the peak amplitude Strouhal number (about 20%). To evaluate the role of aerodynamic forces on a bluff body in noise generation, magnitude of root mean square (rms) sound pressure ( P = P rms/ρU 2 ∞ ) are presented due to the lift and drag fluctuations for the 0 C-0.1 C variations range in Figures 11 and 12. These plots show the share of the acoustic pressure of the lift and drag from the total radiated pressure from configuration based on Equation (15), drawn separately to better understand the values. According to Figure 11, it is evident that the sound waves radiates by the lift fluctuations at the direction ±90 • from the sound source. It is also shown in Figure 12 that the direction of the sound field waves radiates by the drag fluctuations are 0 • and 180 • . From Figures 11 and 12, it can be seen that the noise emitted by the body is dipolar, which is consistent with Curl's [36] result. The peaks of the graphs follow their respective forces. The major changes in P are related to changes in aerodynamic forces. It is clear that in all cases, the magnitude of P for the sound source from the drag signal is small than from the lift signal. This is due to the small amplitude of drag fluctuation, as compared to the amplitude of lift fluctuation. It can be found that in the interval 30 • -150 • and 210 • -330 • , the highest proportion of noise produced is due to lift fluctuations and at angles outside this range, the highest contribution is related to drag fluctuations. This result may indicate that if the microphone approaches 0 • and 180 • angles, the drag contribution to noise production will increase. Since the noise received from the drag fluctuations is far below the lift, it is expected that a lower SPL will be recorded. Sun et al. [58] have shown that with the microphone approaching angles approaches 0 • and 180 • the SPL is significantly reduced. Moreover, with respect to Figure 11, where the noise of the rod is less than that of the Rod-Airfoil, the two graphs also show that the rod has less sound pressure fluctuations.   The changes in NSP are examined by applying vibration to the cylinder for intervals greater than 0.1C (0.1C < A ≤ 1C). Figure 13 shows that increasing displacement decreases the NSP. Increasing the displacement reduces the role of the leading edge in noise generation, and along with that, the trailing edge also becomes a factor for noise generation. In A = 1C, the NSP is reduced by nearly 35 dB. The noise level changes in the range of 0.1C < A ≤ 1C are not significant, so it is expected that the change pattern of noise will be repeated by increasing the displacement rate in 1C ≤ A. The reason for this behavior is that the rod is far enough away from the chord at A = 0.2C, and the transient flow of the rod has low effect on the airfoil, so it is expected that increasing the amplitude more than 2C has no significant effect on the noise level.

Effect of Flow Control on Aerodynamic Coefficients
Hereafter, the control method's effect on dimensionless coefficients of lift and drag will be discussed. Hence, lift and drag coefficient values of the Rod-Airfoil have been produced. The average values are not indicative of the generated sound, But the generated sound field is due to fluctuations of the forces applied to the bluff body. According to the results of Figures 11 and 12, the current position of the microphone, drag has the lowest domain of noise generation. Given the fact that dipole noise is a function of C L rms then variations in C L rms are examined by applying control. Mean drag, known as a structural destructive force, and mean lift have also been investigated. rms and mean values have been reviewed as absolute values according to Equations (18)- (21). The values have been recorded in the figures studied as compared to the case with no control.
C Lmean = C L mean X C L mean X=0 (20) where C L rms X=0 , C D mean X=0 , and C L mean X=0 represent the fluctuations of the lift, the mean drag, and lift in the uncontrolled case, respectively. Moreover, C L rms X , C D mean X and C L mean X respectively indicate the lift fluctuations, mean drag, and lift for the studied cases (0.02C < A ≤ 1C). Figure 14a has illustrated the airfoil and the rod lift coefficient fluctuation values. It is observed that C Lrms has generally decreased. Since variations in lift coefficient fluctuations are dominant, then it is expected that they will be similar to variations in the sound pressure level. It may be observed in this figure that the minimum recorded value of the lift coefficient fluctuations occurs in the airfoil at 0.2C ≤ A ≤ 1C, which has decreased by about 48% compared to the case with no control. The amount of reduction in the rod's lift coefficient fluctuations at A C = 1 has decreased by about 46% in comparison to the case with no control. Figure 14b illustrates the averaged drag coefficient. The results generally indicate increase in the airfoil's average drag coefficient after applying control by 7% at 0.02C ≤ A ≤ 0.1C and 2% at 0.2C ≤ A ≤ 1C, and reduction in the rod's average drag coefficient by 16%. Moreover, vibration in the rod increases the mean lift coefficient, which is aerodynamically desirable (Figure 14c).  Then, further investigation of flow pattern and the force applied to the system will be made. Since the flow characteristics in 0.2C ≤ A ≤ 1C are similar to 0.02C ≤ A ≤ 0.1C, so the flow characteristics are only examined for the second range. Figure 15 has presented the average static pressure. It is highly important to study static pressure contours because it clearly demonstrates the greatest static pressure variation domains is surrounded the airfoil. Pressure variations around the airfoil lead to variations in drag and lift forces. It is clear from the figure that the harmonic motion of the rod causes the pressure distribution pattern to change, which increases the pressure difference in both sides of the airfoil. Unsteady pressures which support lift dipole are the most important and the most common factor in the generation of noise [43]. This unsteadiness is mainly created in the vicinity of the solid boundary and its surrounding. Figure 16 has depicted pressure disturbance contour (P rms ). These values are obtained from P rms = 1 N ∑ N n=1 P − P 2 . As is evident in the figure, the highest RMS pressure level occurs at the airfoil leading edge. Nevertheless, this pressure stress significantly helps the noise, for this stress is located adjacent to the solid boundary. Hence, it can support dipole sources, and it is in fact not relevant to the quadrupole shape. The maximum recorded intensity can indicate the source of the noise. it is necessary to mention the emitted noise does not depend only on the amplitude of the pressure fluctuations, i.e., the coherence of the fluctuations and the shape of the close by solid surface may play a role, but pressure fluctuations are predominant It is evident from Figure 16 that the airfoil leading edge is the cause of sound generation and this zone of tension is located adjacent the airfoil due to the motion of the rod. It should be noted that after the application of control, the area of P rms decreases by increasing the amount of rod displacement so that it would decrease by nearly 50% at A/C = 0.1. Considering the previously discussed issues, if the airfoil is not located at the downstream side of a bluff body, its trailing edge noise might be dominant. It is known that at x/C = 0.1, the intensity of P rms is also generated at the edge of the airfoil trailing edge. This indicates that the noise recorded in Figure 10 for A/C = 0.1, in addition to the leading edge, is also due to the trailing edge. Vorticity contours can perfectly illustrate the boundary layers, the shear layer growth, and the instability of vortices that eventually turn into von Kármán street. The occurrence of this periodic structure is the main cause of the pressure difference, the fluctuations of the forces applied to the system, and the noise generation. Figure 17 show the vorticity contours after applying the control method (0C ≤ A ≤ 0.1C). Vorticity contours are extracted 0.45 s after the start of the analysis. According to Figure 17, in the case of A = 0C, it is clear that the turbulent flow has been applied to the leading edge, which has caused the leading edge to be dominant in noise generation. Vibration in the rod causes the vortex shedding phenomena to change due to vortex roll-up distance. According to Figure 17 which shows vibration apply in the rod, the created vortexes occur smaller and closer to the back surface of the rod, resulting in a dramatic increase in the St [59] due to what happened in Figure 9. It is clear that the vibration in the rod has caused the downstream vortex of the rod not to penetrate on the leading edge as in an uncontrolled case. So, the fluctuations of the flow near the leading edge are reduced. The flow characteristics are relatively similar in different cases, so the noise reduction is relatively same for different cases. An important point in A = 0.1C is increasing the vortex size created at the trilling edge which increases its role in noise generation. This result can be seen in Figure 16, with increasing pressure fluctuations in the mentioned region. Accordingly, similar results are expected to occur at A = 0.1C, in other words, increasing the role of the trailing edge in increasing noise along with reducing the role of the leading edge.

Conclusions
This study demonstrates that the vibrating rod method is capable of reducing the generated noise. This method is potentially used to reduce the noise in vertical-axis wind turbines when they are located downstream of a bluff body. As a general and final conclusion, it may be stated that any factor contributing to control of the airfoil leading edge vortices will result in noise reduction. For example, corrugating the airfoil leading edge and rotating the rod (according to the literature), controlled or deviated vortices. Noise reduction capability is advantageous when these practical programs are truly taken into consideration. Currently, the experimental data on the effect of the vibrating rod on the Rod-Airfoil aeroacoustics are not available; therefore, reliable evaluation of the proposed approach is optimistic for real cases. Moreover, 3D simulation for considering the effects of spanwise with the LES turbulence model can be a useful idea to ensure the results.
In this study, the performance of 2D models was evaluated and the results were sufficiently reliable. However, to better understand and improve these models, our team plans to evaluate Gaussian function and 3D models for more complex configurations in future studies. Funding: This research received no external funding.

Conflicts of Interest:
The authors declare no conflict of interest.