OJFD  Vol.7 No.4 , December 2017
A Numerical Approach to Possible Identification of the Noisiest Zones of a Wall Surface with a Flow Interaction
This paper examines the use of proper orthogonal decomposition (POD) and singular value decomposition (SVD) to identify zones on the surface of the source that contribute the most to the sound power the source radiates. First, computational fluid dynamics (CFD) is used to obtain the pressure field at the surface of the blade in a subsonic regime. Then the fluctuation of this pressure field is used as the input for the loading noise in the Ffowcs Williams and Hawkings (FW&H) acoustic analogy. The FW&H analogy is used to calculate the sound power that is radiated by the blade. Secondly, the most important acoustic modes of POD and SVD are used to reconstruct the radiated sound power. The results obtained through POD and SVD are similar to the acoustic power directly obtained with the FW&H analogy. It was observed that the importance of the modes to the radiated sound power is not necessarily in ascending order (for the studied case, the seventh mode was the main contributor). Finally, maps of the most contributing POD and SVD modes have been produced. These maps show the zones on the surface of the blade, where the dipolar aeroacoustic sources contribute the most to the radiated sound power. These identifications are expected to be used as a guide to design and shape the blade surface in order to reduce its radiated noise.

1. Introduction

Reducing the noise produced by the interaction between a turbulent subsonic flow and a solid’s surface can be difficult, especially in industrial configurations where the noise level is a quality and selection criterion. When considering subsonic turbomachines, this noise can be reduced if the geometry of the profile is appropriate. However, there is no general solution, and each fan configuration needs a specific design. Therefore, most of the time acoustic engineers use a case by case approach, either experimentally or numerically [1] [2] . In this context, identifying and characterizing the region responsible for the noise interaction becomes a challenge. Since the acoustic source is connected to the flow some studies for example separated the turbulent field into acoustically radiating and non-radiating components using Navier-Stokes equations [3] [4] . Other approaches are based this time on the decomposition of the sound radiated itself in order to identify the most radiating zones as the case in our investigation.

This paper proposes a method of identifying and analyzing these zones based on the combination of Ffowcs Williams and Hawkings’ acoustic analogy (FW&H) [5] and either proper orthogonal decomposition (POD) or singular value decomposition (SVD) [6] . The SVD approach used in this paper is the application of the POD to a matrix formed by the product of the correlation matrix obtained from the elementary acoustic pressure to a receiver and it transposed. This study will develop the POD approach in more detail.

The pressure fluctuations obtained from a CFD approach are used as input data for the FW&H acoustic analogy [5] . For that, the large eddy simulation (LES) method [7] [8] proposed by OpenFOAM [9] is used. OpenFOAM’s popularity has been growing rapidly for several applications, with an increasing number of users over the last few decades. Numerous investigations have been published in the CFD field. For example, de Villiers [8] investigated the influence of subgrid scale models in the method implemented in OpenFOAM. This investigation was carried out on several flow configurations around a rigid solid. The noise generated by the interaction between an airflow and the rearview mirror of a car was also studied. The results are in agreement with the experimental results. Furthermore, Dimtry et al. [10] modeled the flow around a cylinder with Smagorinsky’s subgrid models and the one-equation turbulent kinetic energy model implemented by de Villiers. The results are in agreement with experimental data [11] [12] . Several other studies in the literature have demonstrated OpenFOAM’s effectiveness in modeling complex configurations [13] [14] . In acoustics, in addition to de Villiers’ work, Olivier’s work [15] addressed the noise of the trailing edge of NACA 0012. It used data from the OpenFOAM CFD to supply the acoustic source of the Curle analogy. These acoustic results are in agreement with Herr et al.’s results [16] . The subgrid model of one of de Villiers’s turbulent kinetic energy equations will be used for this study for the CFD.

POD was first introduced by Lumley in 1967 [17] in fluid mechanics as an objective method to identify and extract coherent turbulent flow structures. Building on this idea, Bonnet et al. [18] used POD in conjunction with a stochastic estimation method to identify the key characteristics of a turbulent field on the basis of a reduced number of measurements. POD was also employed either for optimal flow control [19] [20] or reduction model development [21] [22] . In the latter application, the Navier-Stokes equations are projected on the eigenmodes using the Galerkin projection. The POD method was used in many aeroacoustics investigations, mostly to describe the noise radiated by turbulent flows (e.g. jet noise, wake noise) [23] [24] . The optimal decomposition basis velocity field is then used to estimate acoustic noise. Druault et al. [25] [26] employed this method in the far field region to separate the acoustic contribution of the most energetic structures from the residues caused by the clean flow noise (i.e. the noise from velocity fluctuations). They noticed that 99.8% of the acoustic energy comes from the dominant modes and 0.2% from acoustic residues. Hekmati and Ricot [27] reached the same conclusion by applying the Druault et al.’s method [25] [26] for wake noise generated by the blades of an axial fan. Furthermore, Gleggs and Devenport [28] used POD to describe the input to the turbulence using a set of statistically independent modes of velocity. The wall pressure is linked to the input turbulence for the evaluation of the acoustic interaction sound pressure.

SVD is generally used to search for the propagation operators such as the Green function [6] [29] , or to solve inverse problems in acoustics [30] [31] [32] [33] . To the authors’ knowledge, SVD has not yet been used to locate radiating areas on a moving surface due to its interaction with flow.

In the light of previous works, this study applies the POD and SVD methods to the problem of aeroacoustic noise generated by the interaction between a stationary blade and a turbulent flow in a channel. The objective is to understand the link between the decomposition modes of POD or SVD and the noisiest zones of the blade’s surface. The methodology relies on a three-step methodology: 1) The internal flows of the centrifugal fan are modeled using LES method [7] [8] . OpenFOAM extend 3.2 [9] software is used. The objective is to estimate the wall pressure fluctuation on the blade. 2) The previous wall pressure fluctuation of the blade is used to estimate the loading noise from the FW&H analogy [5] [34] . 3) Finally, POD and SVD are used to extract the most important acoustic modes and visualize them in order to identify the zones that radiate the most on the blade’s surface. This paper is organized as follows. First, the estimation of the acoustic field based on the FW&H analogy is introduced. Then both the POD and SVD methods are developed. Finally, an application is demonstrated using a stationary blade.

2. Theory

2.1. Estimation of the Acoustic Field (FW&H)

The acoustic approach is based on the FW&H analogy [5] and is used as the reference approach. The subsonic regime is considered in this paper. The finite thickness of the blades is neglected; only loading noise ( p L ) generated by the fluctuating wall pressure on the blade is considered in the following example.

To alleviate the computation time problem, Formulation 1A proposed by Farassat [34] can be used. In this approach, the receiver time derivative in Formulation 1 [34] is transformed into a retarded time derivative. This has also the great advantage of permuting the time derivative and the integral. Moreover, no derivation of the integral is needed. For the case in which the source and the receiver are stationary and the propagation medium is at rest, the Farassat Formulation 1A [34] in the far field becomes:

p L ( x , t ) = 1 4 π c S y [ l ˙ r r ] τ d S y (1)

where S y is the source surface, r = r = y x is the distance between the source position y on the blade’s surface and the receiver position x , c is the speed of sound of the acoustic medium at rest, l r = p n r / r with n is the unit normal vector to the source’s surface, and p is the wall pressure fluctuation of the blade’s surface obtained by CFD calculation. M is the Mach vector number and M r = M r / r and [ ] τ indicates that all the integrands should be evaluated at the retarded time τ = t r / c , with receiver time t.

As mentioned previously, Formulation 1, proposed by Farassat, has the main advantage of avoiding the spatial derivatives; however, the receiver time derivative is maintained on l r and M r . Implementing this operation is complex, and the computation time increases. To evaluate Equation (1), two computational approaches are available in the literature [35] [36] using either retarded time or advanced time. The latter is usually chosen when the aerodynamic data comes from CFD computations, as in this study. The received acoustic pressure must then be determined while an irregular receiving time discretization appears despite the regular emission time. An interpolation is thus necessary to obtain the received sound pressure at regular time.

The advanced time approach and the Lagrange interpolation are used in this paper. The calculation of the sound pressure at receiver x in the far field and in a free medium allows for the calculation of root mean square sound pressure. When considering the loading noise (Equation (1)), the mean square acoustic pressure reads:

p L 2 ( x ) = p L ( x , t ) p L ( x , t ) T 0 (2)

where T 0 is the temporal average over the time period T 0 . Thus, the radiated sound power estimated from far field microphones (i.e., receivers) located on a spherical surface encompassing the source can be written:

P = S x p L 2 ( x ) ρ c d S x 1 ρ c x p L 2 ( x ) Δ S x (3)

where S x is the receiving surface, Δ S x is the elementary surface associated with receiver x , and ρ is the density of the surrounding fluid medium.

2.2. Proper Orthogonal Decomposition (POD) Approach

Generally, POD used in aeroacoustics is not based on acoustic analogies. The approach developed in this paper is a combination of the dipole term (loading noise) of the FW&H analogy and the POD theory. If one considers a dipole located at y i , the sound pressure received at point x at time t according to Equation (1) can be written:

p L ( x , y i , t ) = 1 4 π c [ l ˙ r r ] τ Δ S y i (4)

where Δ S y i is the i-th elementary surface of the source located at y i . Here, this elementary surface comes from the discretization of the source’s surface for the LES calculation. Since the sound pressure received at one point corresponds to the contribution of all the dipoles located on the wall surface, two matrices A and W o b s are defined by

W o b s = 1 N A A T (5)


A = ( p L ( x , y 0 , t e 0 ) p L ( x , y 0 , t e 1 ) p L ( x , y 0 , t e N 1 ) p L ( x , y 1 , t e 0 ) p L ( x , y 1 , t e 1 ) p L ( x , y 1 , t e N 1 ) p L ( x , y m 1 , t e 0 ) p L ( x , y m 1 , t e 1 ) p L ( x , y m 1 , t e N 1 ) )

where t e j represents the emission time for the dipole source located at y i ( i = 0 , 1 , , m 1 ) at time step j ( j = 0 , 1 , , N 1 ) . Each column vector of the matrix A is the sound contribution of all dipoles at a given reception time, while each row vector represents the sound pressure of one single dipole source along the receiving time. Matrix W o b s in Equation (5) is the correlation matrix of the sources for receiver x . It is symmetric, real, positive definite, and spatial. Its eigenvalues (modes) are thus real, positive, and space-dependent. When developing a modal basis for matrix W o b s , one defines λ = diag ( λ 0 , λ 1 , , λ m 1 ) and ϕ = [ ϕ 0 , ϕ 1 , , ϕ m 1 ] , the diagonal matrix of eigenvalues and the matrix of eigenvectors at receiver x respectively. Each column ϕ i is the eigenvector associated with the eigenvalue λ i at receiver x . Then for every receiver x , the problem to be solved is the following eigenvalue problem:

W o b s ϕ = λ ϕ (6)

If ones multiplies each member of Equation (6) to the right by the transpose of the matrix of eigenvectors and considers normalized modes so that ϕ ϕ T = I , the expression for the correlation matrix is given by:

W o b s = k = 0 m 1 ( λ k Φ k ) Φ k T (7)

The correlation matrix W o b s is then written as a sum of independent matrices defined as spatial autocorrelation patterns with proper modes as components. Since the eigenvectors form an orthonormal basis of the source space, p L ( x , y i , t ) can be written as:

p L ( x , y i , t ) = k = 0 m 1 α k ( t ) Φ k , i (8)

where α k ( t ) = α k ( x , t ) = i = 0 m 1 p L ( x , y i , t ) Φ k , i are the temporal modal amplitudes or the projection coefficients on the modal basis. Φ k , i = ϕ k ( x , y i ) is the i-th component of the k-th eigenvector ( ϕ k ). Thus, coefficients are the root mean square of the acoustic pressure projected on the Φ k ( y ) axis in the source space (i.e., the blade). According to Merces’ theory [37] , the projection coefficients form an orthogonal basis of the temporal space, and its root mean square α k ( t ) α k ( t ) T 0 corresponds to the eigenvalue λ k . The eigenvalues represent the square of the sound pressure projected on the Φ k ( y ) axis in the source space. The total sound pressure from the loading noise radiated to the receiver x at time t is evaluated from Equations 1 and 8. It is given by:

p L ( x , t ) = k = 0 m 1 ( α k ( t ) i = 0 m 1 Φ k , i ) (9)

Considering the orthogonal eigenvectors, the quadratic pressure at receiver x is:

p L 2 ( x ) = k = 0 m 1 ( i = 0 m 1 Φ k , i ) 2 λ k (10)

In this proper orthogonal decomposition, each mode does not contribute equally to the total quadratic pressure. The latter could then only be evaluated by taking into account the modes that contribute the most. This is reported by the accumulated acoustic energy of the first q modes, E A c q , divided by the total acoustic energy of all modes:

E A c q = k = 0 q 1 λ k k = 0 m 1 λ k (11)

Once the modes that contribute the most are identified, summations in Equations (3) and (10) are limited up to q 1 instead of m 1 . Thus, Equations (9) and (10) become respectively:

p L ( x , t ) k = 0 q 1 ( α k ( t ) i = 0 m 1 Φ k , i ) (12)

p L 2 ( x ) k = 0 q 1 ( i = 0 m 1 Φ k , i ) 2 λ k (13)

In the previous equations, the spatial eigenvectors Φ k , i give information on the acoustic radiation of all dipole sources distributed over the surface S. Thus, considering Equation (13), the acoustic power defined by Equation (3) becomes:

P 1 ρ c l = 0 N o b s 1 ( Δ S l k = 0 q 1 λ l k ( i = 0 m 1 Φ l k , i ) 2 ) (14)

where Δ S l = Δ S x is the elementary surface of the ℓ-th receiver. In this study, Δ S l is a constant. N o b s is the number of receivers, λ l k is the k-th eigenvector of the ℓ-th receiver, and Φ l k , i is the i-th component of the eigenvector associated with the k-th eigenvalue λ l k .

2.3. Singular Value Decomposition (SVD) Approach

SVD investigates sound generation (i.e., the eigenvalues and the eigenvectors resulting from the POD) independent of the receiver x . A global matrix W S V D gathering all the correlation matrices W k ( 0 k < N o b s ) of the N o b s receivers is built. The global matrix of dimension ( N o b s × m ) × m is defined by:

W S V D = [ W 0 W 1 W N o b s 1 ] (15)

Here, the receivers are distributed over a sphere of radius R around the source. The radius is large enough that the far field assumption is verified. The position of the receivers over the sphere is done according to ISO 3745 [38] . SVD consists of decomposing the W S V D matrix in the following form:

W S V D = U σ V T (16)

where σ is the diagonal matrix of singular eigenvalues, and U and V are the matrices of left and right eigenvectors respectively. U accounts for the difference of acoustic radiation of the sources in the free field of all receivers, and V provides the average information on acoustic radiation of the sources in the free field of all receivers. One can show that U and V can be obtained by applying the POD to matrices W S V D T W S V D and W S V D W S V D T respectively. By retaining the first q energetic modes only, the expression of sound power (Equation (3)) based on the SVD can finally be expressed as:

P 1 ρ c k = 0 q 1 ( σ k i = 0 m 1 V k , i l = 0 N o b s 1 ( Δ S l j = m l m ( l + 1 ) 1 U k , j ) ) (17)

where σ k is the k-th eigenvalue of the matrix W S V D , V k , i and U k , i respectively represent the right and left i-th component eigenvectors, respectively associated with the eigenvalue σ k .

3. Application

In Figure 1, a centrifugal fan blade based on the modified NACA12 profile is used to illustrate our approach. The chord and wingspan equals l = 0.056 m and l 1 = 12 / 25 l respectively. The blade is placed in a l 1 height periodic channel whose periodic faces are separated by the width l 2 = 17 / 50 l . The cross section area of the channel equals l 1 × l 2 , where l 1 and l 2 are the length and the width respectively. This area remains constant from the inlet until the

Figure 1. Geometry and boundary conditions.

outlet (Figure 1). The numerical domain comprises three volumes separated by interfaces: the inlet, the blade, and the outlet volumes, whose lengths equal l 3 = 9 / 5 l , l 4 = 7 / 5 l and l 5 = 17 / 5 l respectively. The fluid is pure air and its physical properties are estimated at 25˚C ( ρ = 1.2 kg / m 3 and μ = 1.831 × 105 Pa・s). An airflow of U = 70 m / s is imposed on the inlet. The periodicity boundary conditions are applied to the periodic face and the no slip boundary conditions are imposed on the walls (blade, top and bottom). The purpose here is to investigate the possibility of both the POD and SVD approaches for identifying zones of the blade that are responsible for noise generation. Only frequencies in the 50 Hz to 10 kHz range are considered.

The calculation is performed on a hybrid spatial discretization realized with the free software Salome [39] . The discretization is composed of a structured and inhomogeneous mesh in the input and output volumes. In the blade volume and far from the blade, the mesh is unstructured .The mesh near the blade is structured with eight layers of geometrical progression 1.1 and total thickness δ = 0.0001 m . Around the blade, the dimensionless variables x + = 54 (flow direction), z + = 26 (spanwise direction), and y + = 1 (direction normal to the walls) are used. These values are sufficient according to Sagaut criteria [40] to correctly predict the boundary layer behavior when using LES method. This is also in accordance with de Villiers investigation [8] that employed a mesh satisfying these criteria and obtained numerical data close to the experiments. Recently, in one of our accepted articles for publication [41] , the influence of discretization on POD and SVD approaches has been studied. The coarse, medium and fine meshes describe a same noisiest area of the surface of the source. Thus, the mesh is composed of 1,100,000 cells, including m = 30400 cells on the blade.

The large eddy simulation method is used to simulate the internal flow channel with the free software OpenFOAM [9] . The one-equation subgrid model for turbulent kinetic energy is used [8] . The linear system is solved with the iterative preconditioned conjugate gradient method (PCG) with preconditioning DIC (diagonal incomplete Cholesky) for the pressure field, and the stable biconjugate gradient method (BiCGStab) with preconditioning DILU (diagonal incomplete LU) for the velocity field and turbulent kinetic energy. The backward second order scheme is used for the temporal resolution with a time step of Δ t C F D = 2.95 × 10 7 s for a Courant-Friedrichs-Lewy (CFL) number lower than 0.2. The gradient terms are calculated with a Gauss linear scheme, and the divergence terms are calculated with that of a Gauss vanLeerV 0.5 scheme for velocity and a Gauss upwind scheme for turbulent kinetic energy. Three weeks were necessary to reach convergence using 24 processors on Compute Canada-Sherbrooke [42] .

The second invariant of the velocity gradient tensor, named Q-criterion, was introduced by Hunt et al. [43] to better visualize the coherent structures of the flow characterized by a positive Q-criterion value. One observes in our configuration that there are no coherent structures in the region upstream the blade (Figure 2), despite the 5% intensity level imposed at the inlet. Coherent structures appear from a region located at 30% of the blade chord and develop gradually until the trailing edge. In this location, a strong intensification of the production of coherent structures that are then transported by the flow is particularly noticeable.

Once the calculation is converged, the wall pressure fluctuations on the blade are saved at time interval Δ t a c o u s t = 10 Δ t C F D . N = 2500 samples of the pressure fluctuations are then the input data of both POD and SVD methods. The sound pressure, correlation matrix and modes are calculated for receivers placed on a sphere of radius R around the source, according to ISO 3745. The distance R is such that the far field hypothesis is verified (i.e. k R 1 k l where l is a characteristic length of the source and k = 2 π f c / c is the wave number based on the smallest cut-off frequency that is f c = 50 Hz ). For this investigation, N o b s = 20 receivers were placed on the sphere of radius R = 6 m encompassing the source (Figure 3(a)). To better see the directivity of the radiation of the blade, 50 receivers were placed along a circle having the same radius as the


Figure 2. Isocontour Q-criterion from 0 to 1000: (a) All volume, and (b) Zoom on the middle volume and contour lines at the trailing edge of the blade.

(a) (b)

Figure 3. Receivers positions: (a) Receivers on the sphere and (b) Receivers position on the circle.

sphere. The circle is in the plane perpendicular to the chord of the blade (i.e. in the plane ( y , z ) . The 50 receivers are identified by their polar coordinates ( R , φ ) . The angular origin is from the z axis. The positions of 4 of the 50 receivers are illustrated in Figure 3(b).

3.1. POD Analysis and Interpretation

The POD principle is to search an orthonormal basis of the m elementary sources of the discretized radiating blade’s surface. From this orthogonal basis, the q eigenvectors that contribute the most to the acoustic radiation at a given receiver are identified. The mapping of these q eigenvectors on the blade’s surface allows identifying its most radiant zones. Consequently, the application of the POD method on our stationary blade in the channel aims to minimize the number of modes necessary to understand the most radiant zones of the blade due to its interaction with the turbulent flow.

Indeed, Figure 4 presents the accumulated acoustic energy of the q first modes (Equation (11)) at receivers φ = 0 , φ = 28 , φ = 56 , and φ = 91 . One observes that the first 10 eigenvectors capture 99.60% of the total acoustic energy, whatever the receivers. As a result, these 10 eigenvectors could then be employed for the reconstruction of the total loading noise. This leads to a reduction in the number of modes (from m = 30400 to q = 10 ) to be analyzed in order to identify the noisiest zones.

For example, Figure 5 shows the temporal evolution of the sound pressure of the loading noise reconstructed with q = m (Equation (1)) and (Equation (9)), and q = 10 (Equation (12)) at receivers φ = 0 , φ = 28 , φ = 56 , and φ = 91 . Acoustic pressures constructed with the first 10 modes (Equation (12)) of each receiver are in accordance with the reference acoustic pressure (Equation (1)). The maximum relative error equals 1.7% at receiver φ = 91 . The relative

Figure 4. POD percentage accumulated acoustic energy (eigenvalues, Equation (11)) of the first q modes divided by the total energy of all modes for receivers φ = 0 , φ = 28 , φ = 56 , and φ = 91 .

Figure 5. Comparison between total loading noise ( q = m ) and loading noise reconstruction using the first 10 eigenvalues at receivers: (a) φ = 0 , (b) φ = 28 , (c) φ = 56 , and (d) φ = 91 .

error remains less than 0.1% for the other receivers.

One also observes that the temporal evolution of the sound pressure depends on the receiver location (Figure 5), which characterizes the directivity of radiation. It results in the POD being dependent on the receiver. For example, the maximum amplitude of sound pressure is almost equal to 0 for the φ = 0 receiver, and is equal to 4 × 10 3 Pa for the receiver φ = 91 . The amplitude of the sound pressure signal is very small for receiver located at 0˚ in comparison with the others. In order to better see this signal a subfigure has been added on Figure 5(a).

Consider the receivers in Figure 3(b). We calculated the mean quadratic sound pressure with all modes (Equations (2) or (10)) and with the first 10 dominant modes (Equation (13)). First, One observes that the acoustic radiation is dipolar (Figure 6). The lowest sound pressure is noticed for receivers located at φ = 0 and φ = 180 , while the highest radiations are for receivers at φ = 90 and φ = 270 (Figure 6). This observation is in agreement with the directivity defined by Farassat [34] . According to this formulation (Equation (1)) the scalar product between the normal surface of the blade and the direction with the receiver is considered in the loading noise evaluation. Thus this product becomes null for the two former positions (i.e. φ egals 0˚ and 270˚) while equals unity for the two latter (i.e. φ egals 0˚ and 270˚).

In addition, the two curves of the mean quadratic sound pressure are similar, since the relative error remains less than 3% (Figure 6). This confirms the analysis of the reconstruction of sound pressure with the first 10 most important modes (Figure 5). When considering each term of the mean quadratic sound pressure (Equation (13)). We noted that the contribution of the modes does not necessarily appear in the ascending order (Figures 4-7. For example, considering receiver φ = 91 , the contribution of mode 1 is greater than the contribution of mode 0, even if the latter accumulates 71% of the energy, while mode 1

Figure 6. Directivity estimation in dB between mean quadratic sound pressure ( q = m ) and mean quadratic sound pressure using the first ten eigenvalues at receivers of the Figure 3(b).

Figure 7. Mean quadratic sound pressure (product of eigenvalues and sum of the components of the eigenvector, Equation (13) for each of the first 15 modes of receivers φ = 0 , φ = 28 , φ = 56 , and φ = 91 .

only accumulates 23% (Figure 4). This contradiction is due to the phase shift of the components of an even eigenvector in calculating the mean quadratic sound pressure (Equation (13)). As a consequence the modes must be classified in ascending order when calculating either the mean quadratic sound pressure or the acoustic power. Thus, the reclassification of the modes was carried out for each receiver in Table 1.

One considers the reconstructed sound power (Equation (14)) with q = 10 and q = 5 . As shown in Table 2, the relative error of the acoustic power of the first five modes (Equation (14)) with respect to all modes (Equation (3)) is low and even lower with the first 10 modes (Equation (14)). This confirms that all useful acoustic information needed here is contained in the first 10 modes.

It is possible to project the eigenvectors that contribute the most to the acoustic power onto the blade surface, since the quadratic mean pressure is proportional to the sound power (Equations (13) and (14)). For example, when considering receivers φ = 0 , φ = 28 , φ = 56 , and φ = 91 , one observes that the POD method highlights the trailing edge as the region characterized by a great amplitude of the component Φ k , whatever the side of the blade (Figure 8 and Figure 9). The other regions of the blade do not seem to be affected by such amplitudes. Hence, mode 2 for receivers φ = 0 , 56˚, and 91˚, and mode 3 for receiver φ = 28 highly contribute to the acoustic power. The same observation is made for modes 3, 2, 4 and 9 for receivers φ = 0 , φ = 28 , φ = 56 , and φ = 91 respectively (last column of Figure 8 and Figure 9).

These previous observations are in agreement with the Q-criterion distribution (Figure 2), since according to this parameter, the turbulence is high near the trailing edge. The high velocity and pressure fluctuations in this region

Table 1. Classification of modes according to their contribution to the calculation of mean square acoustic pressure.

Table 2. Results of the sound power reconstructed with POD and relative error.

Figure 8. Mapping POD eigenvector on the overpressure side of the blade according to the classification given in Table 1: line 1, receiver φ = 0 ; line 2, receiver φ = 28 ; line 3, receiver φ = 56 ; line 4, receiver φ = 91 , and columns rank from most to least important modes from left to right.

would therefore make this region responsible for the loading noise.

It was shown here that both the radiated sound pressure and the sound power can be reconstructed with only a few POD modes. Instead of using all the modes

Figure 9. Mapping POD eigenvector on the underpressure side of the blade according to the classification given in Table 1: line 1, receiver φ = 0 ; line 2, receiver φ = 28 ; line 3, receiver φ = 56 ; line 4, receiver φ = 91 , and columns rank from most to least important modes from left to right.

proposed by the POD method (i.e. m = 30400 ), the first 10 dominant POD modes only were sufficient. In addition, the dipole character of the blade of radiation was noticed. The mapping of the major eigenvectors for several receivers demonstrated that the trailing edge is the region that radiates the most due to the high level of turbulence.

3.2. SVD Analysis and Interpretation

Like the POD, SVD methods based on CFD calculations and acoustic analogies were applied to identify the zones on a stationary blade in a channel that contribute the most to the sound power radiated in a subsonic regime. Using the SVD method, radiated sound power and pressure can be recovered using only the first few modes.

Figure 10 represents the acoustic energy that the q first modes accumulate divided by the total energy (Equation (11)). It emerges that 99.90% of the total acoustic energy is contained within the first 10 modes, which causes a reduction in the number of modes to be analyzed in order to understand the noisiest zones of the blade. The reconstructed sound power (Equation (17)) using the first 10 modes only is given in Table 3 and is compared with the sound power calculated

Figure 10. SVD percentage accumulated acoustic energy of the first q modes E A c q divided by the total energy of all modes.

Table 3. Results of the sound power reconstructed with SVD and relative error.

using all the modes. The relative error equals 0.43%.

As shown previously for the POD method, the most energetic eigenvalues do not necessarily contribute the most to the reconstruction of the acoustic power. Modes 6, 10, 4, and 9 in order of their importance contribute the most (Figure 11).

In addition, the acoustic power estimated with only mode 6 equals 1.15 × 10 7 W (i.e., 50.6 dB) representing 93% of the total acoustic power. When adding the contribution of mode 10 to the contribution of mode 6, 96% of the total acoustic power is recovered. Thus, the number of modes to be analyzed to comprehend the most radiant zones of the blade regardless of receiver changes from m = 30400 modes to two modes (modes 6 and 10).

As for the POD method, it is possible to use singular value decomposition to map the right eigenvectors V k of the four main modes (i.e. modes 6, 10, 4, and 9) on the blade surface as shown in Figure 12. One notices that the amplitude for mode 6 is quite homogenous along the surface. As mentioned previously, this must be due to the turbulence that develops quite homogenously in the blade region; the entire surface therefore contributes equally to the loading noise. When considering the distribution for mode 10, one observes that the trailing edge is characterized by greater amplitudes, which is in accordance with the strong intensification of the turbulence in this region as shown by the

Figure 11. SVD acoustic power radiated by each mode.

Figure 12. Mapping SVD eigenvectors on the overpressure and underpressure side of the blade: line 1, overpressure side; line 2, underpressure side, and columns rank from most to least important modes.

Q-criterion (Figure 2). As a result, this region contributes largely to the radiated acoustic power through mode 10.

As a result, based on this observation, one could expect that the radiated sound power could be reduced by a proper modification of the blade geometry, or surface treatment, in these regions. This conclusion is yet to be validated in future work.

4. Conclusions

A CFD calculation was performed using OpenFOAM in order to estimate the wall pressure fluctuations on the surface of a blade located in a channel. The loading noise was then evaluated using the FW&H acoustic analogy and decomposed using both the POD and SVD methods. It was observed that the sound power and the radiated pressure were recovered using these two methods, even if only a few modes were considered. In the configuration studied here, 10 POD modes or 1 SVD mode can be sufficient to predict the radiated sound pressure. Whatever the approach, the trailing edge of the blade is distinguished. An acoustic treatment or a geometrical modification of this zone can therefore influence the blade’s acoustic radiation.

It was shown that the POD approach provides information about the directivity of the source, which is not the case when using the SVD approach. However, the latter approach would be employed when a global noise reduction is sought, no matter the direction of the sound. On the contrary, if a decrease is preferred in a particular direction, the POD method should be used.


This work was supported by Natural Sciences and Engineering Research Council of Canada (NSERC). The authors also wish to thank Compute Canada in Sherbrooke, Quebec, for its assistance.

Cite this paper
Kone, T. , Marchesse, Y. and Panneton, R. (2017) A Numerical Approach to Possible Identification of the Noisiest Zones of a Wall Surface with a Flow Interaction. Open Journal of Fluid Dynamics, 7, 525-545. doi: 10.4236/ojfd.2017.74036.
[1]   Guedel, A. (2002) Bruit des ventilateurs—Parte 2 [fan noise—part 2]. In: T.I., Ed., Techniques de l’ingénieur, Ref. Num. bm4178, 1-25.

[2]   Embleton, T.F.W. (1963) Experimental Study of Noise Reduction in Centrifugal Blowers. The Journal of the Acoustical Society of America, 35, 700-705.

[3]   Goldstein, M.E. (2005) On Identifying the True Sources of Aerodynamic Sound. Journal of Fluid Mechanics, 526, 337-347.

[4]   Snayoko, S., Agarwal, A. and Hu, Z. (2011) Flow Decomposition and Aerodynamic Sound Generation. Journal of Fluid Mechanics, 668, 335-350.

[5]   Ffowcs Williams, J.E. and Hawkings, D.L. (1969) Sound Generation by Turbulence and Surfaces in Arbitrary Motion. Philosophical Transactions of the Royal Society of London. Series A, Mathematical and Physical Sciences, 264, 321-342.

[6]   Borgiotti, G.V. (1990) The Power Radiated by a Vibrating Body in an Acoustic Field and Its Determination from Boundary Measurements. Journal of the Acoustical Society of America, 88, 1884-1893.

[7]   Da Silva, C.B. and Pereira, J.C.F. (2007) Analysis of the Gradient Diffusion Hypothesis in Large Eddy Simulations Based on Transport Equations. Physics of Fluids, 19, Article ID: 035106.

[8]   De Villiers, E. (2006) The Potential of Large Eddy Simulation for the Modeling of Wall Bounded Flows. PhD Thesis, Imperial College of Science, Technology and Medicine, London.

[9]   OpenFoam (2016) The Open Source CFD Toolbox.

[10]   Lysenko, D.A., Ertesvag, I.S. and Rian, K.E. (2012) Large Eddy Simulation of the Flow over a Circular Cylinder at Reynolds Number 3900 using the OpenFoam Toolbo. Flow Turbulence Combust, 89, 191-518.

[11]   Parnaudeau, P., Carlier, J., Heitz, D. and Lamballais, E. (2008) Experimental and Numerical Studies of the Flow over a Circular Cylinder at Reynolds Number 3900. Physics of Fluids, 20, Article ID: 085101.

[12]   Son, J.S. and Hanratty, T.J. (1969) Velocity Gradients at the Wall for Flow around a Cylinder at Reynolds Numbers from 5 × 103 to 104. Journal of Fluid Mechanics, 35, 353-368.

[13]   Lysenko, D.A., Ertesvag, I.S. and Rian, K.E. (2013) Modeling of Turbulent Separated Flows using OpenFOAM. Computers & Fluids, 80, 408-422.

[14]   Robertson, E., Choudhury, V., Bhushana, S. and Walters, D.K. (2015) Validation of OpenFOAM Numerical Methods and Turbulence Models for Incompressible Bluff Body Flows. Computers & Fluids, 123, 122-145.

[15]   Verhoeven, O. (2011) Trailing Edge Noise Simulations: Using IDDES in OpenFOAM. Master’s Thesis, Delft University of Technology.

[16]   Herr, M., Appel, C., Dierke, J. and Ewert, R. (2010) Trailing-Edge Noise Data Quality Assessment for CAA Validation. Proceeding of 16th AIAA/CEAS Aeroacoustics Conferences, Stockholm.

[17]   Lumley, J.L. (1967) The Structure of Inhomogeneous Turbulent Flows. In: Yaglom, A.M. and Tartarsky, V.I., Eds., Atmospheric Turbulence and Radio Wave Propagation, 166-177.

[18]   Bonnet, J.P., Cole, D.R., Delville, J., Glauser, M.N. and Ukeiley, L.S. (1994) Stochastic Estimation and Proper Orthogonal Decomposition: Complementary Techniques for Identifying Structure. Experiments in Fluids, 17, 307-314.

[19]   Camarri, S. and Iollo, A. (2010) Feedack Control of the Vortex-Shedding Instability Based on Sensitivity Analysis. Physics of Fluids, 9, 94-102.

[20]   Ravindran, S.S. (2000) Reducced-Order Adaptive Controllers for Fluid Flows using POD. Journal of Scientific Computing, 15, 457-477.

[21]   Gloerfelt, X. (2008) Compressible Proper Orthogonal Decomposition/Galerkin Reduced-Order Model of Self-Sustained Oscillation in a Cavity. Journal of Physics of Fluids, 20, 105-122.

[22]   Iollo, A., Lanteri, S. and Dsidri, J.A. (2000) Stability Properties of POD-Galerkin Approximations for Compressible Navier-Stockes Equation. Theoretical and Computational Fluid Dynamics, 13, 377-396.

[23]   Hekmati, A. (2011) Analyse des évènements aeroacoustiques à l'origine des éemissions sonores à partir de simulations numériques. [Analysis of Aerodynamic Events at the Origin of Noise Emissions from Numerical Simulations.] PhD Thesis, Universit\'e Pierre et Marie Curie, Paris.

[24]   Arndt, R.E.A. and George, W.K. (1974) Investigation of the Large Scale Coherent Structure in a Jet and Its Relevance to Jet Noise. Tech. Rep. 74N27505 NASA-CR-138908.

[25]   Druault, P., Yu, M. and Sagaut, P. (2010) Quadratic Stochastic Estimation of Far-Field Acoustique Pressure with Coherent Structure Event in a 2D Compressible Plane Mixing Layer. International Journal for Numerical Method in Fluids, 62, 906-926.

[26]   Druault, P., Hekmati, A. and Ricot, D. (2013) Discrimination of Acoustic and Turbulent Components from Aeroacoustic Wall Pressure Field. Journal of Sound and Vibration, 332, 7257-7278.

[27]   Hekmati, A. and Ricot, D. (2009) Aeroacoustic Analysis of the Automotive Ventilation Outlets using Extended Proper Orthogonal Decomposition. Proceedings of the 15th AIAA/CEAS Aeroacoustics Conference, Miami, 11-13 May 2009.

[28]   Glegg, S.A.L. and Devenport, W.J. (2001) Proper Orthogonal Decomposition of Turbulent Flows for Aeroacoustic and Hydroacoustic Applications. Journal of Sound and Vibration, 239, 767-784.

[29]   Photiadis, D.M. (1990) The Relationship of Singular Value Decomposition to Wave-Vector Filtering in Sound Radiation Problems. Journal of the Acoustical Society of America, 88, 1152-1159. https://doi.org/10.1121/1.399811

[30]   Grace, S.P., Atassi, H.M. and Blake, W.K. (1996) Inverse Aeroacoustic Problem for a Streamlined Body, Part 1. Basic Formulation. American Institute of Aeronautics and Astronautics Journal, 34, 2233-2240.

[31]   Grace, S.P., Atassi, H.M. and Blake, W.K. (1996) Inverse Aeroacoustic Problem for a Streamlined Body, Part 2. Accuracy of Solutions. American Institute of Aeronautics and Astronautics Journal, 34, 2241-2246.

[32]   Nelson, P.A. and Yoon, S.H. (2000) Estimation of Acoustic Source Strength by Inverse Methods: Part I. Conditioning of the Inverse Problem. Journal of Sound and Vibration, 233, 643-668. https://doi.org/10.1006/jsvi.1999.2837

[33]   Nelson, P.A. and Yoon, S.H. (2000) Estimation of Acoustic Source Strength by Inverse Methods: Part II. Experimental Investigation of Methods for Choosing Regularisation Parameters. Journal of Sound and Vibration, 233, 669-705.

[34]   Farassat, F. (2007) Derivation of Formulations 1 and 1A of Farassat. Tech. Rep. TM-2007-214853 NASA/TM-2007-214853.

[35]   Fedala, D., Kouidri, S., Bakir, F. and Rey, R. (2007) Modelling of Broadband Noise Radiated by an Airfoil-Application to an Axial Fan. International Journal of Vehicle Noise and Vibration, 3, 106-117.

[36]   Casalino, D. (2003) An Advanced Time Approach for Acoustic Analogy Predictions. Journal of Sound and Vibration, 261, 583-612.

[37]   Mercer, J. (1909) Functions of Positive and Negative Type and Their Connection with the Theory of Integral Equations. Philosophical Transactions of the Royal Society A, 209, 415-446.

[38]   ISO 3745 (2003) Acoustics-Determination of Sound Power Levels and Sound Energy Levels of Noise Sources using Sound Pressure-Precision Methods for Anechoic Rooms and Hemi-Anechoic Rooms. Tech. Rep. European Standard.

[39]   Salome (2015) The Open Source Integration Platform for Numerical Simulation.

[40]   Wagner, C., Huttl, T. and Sagaut, P. (2007) Large-Eddy Simulation for Acoustics. Cambridge Aerospace Series. Cambridge University Press, Cambridge.

[41]   Kone, T.C., Marchesse, Y. and Panneton, R. (2016) Numerical Approach for Possible Identification of the Noisiest Zones on the Surface of a Centrifugal Fan Blade. In: Nobrega, J. and Jasak, H., Eds., OpenFoam: Selected Papers of the 11th Workshop, Springer, Guimaraes.

[42]   Compute-Canada (2016) Compute Canada in Sherbrooke in Sherbrooke, Quebec.

[43]   Hunt, J., Wray, A. and Moin, P. (1988) Eddies, Stream, and Convergence Zones in Turbulent Flows. Proceeding of the Summer Program in Center for Turbulence Research, 193-208.