22.10.2014 Views

THESE de DOCTORAT - cerfacs

THESE de DOCTORAT - cerfacs

THESE de DOCTORAT - cerfacs

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

4.1 Fundamental validation cases 63<br />

150<br />

210<br />

120<br />

240<br />

90<br />

270<br />

0.03<br />

0.02<br />

0.01<br />

180 0<br />

60<br />

300<br />

30<br />

330<br />

r = 0.05<br />

r = 0.15<br />

Pressure Fluctuation (Pa)<br />

0.04<br />

0.02<br />

0<br />

−0.02<br />

Analytics<br />

AVSP−f<br />

−0.04<br />

−0.2 −0.1 0 0.1 0.2<br />

X axis (m)<br />

a) | ˆp| Directivity of a dipole b) Pressure fluctuation<br />

for two different radii.<br />

along the horizontal e 1 axis<br />

[ —— Analytics, (◦△) AVSP-f ] [ —— Analytics, (◦) AVSP-f ]<br />

Figure 4.7: Analytical solution vs Numerical solution.<br />

turbulent flame is acoustically equivalent to a set of different monopole sources, each of them<br />

radiating at different strength, phase and frequencies. Therefore, for a last test of this first<br />

validation, it is interesting to evaluate the noise radiation and directivity that results from three<br />

monopoles with arbitrary positions and phase shifts. The contribution of a third pole is just<br />

ad<strong>de</strong>d as a third term in Eq. (4.6). After <strong>de</strong>veloping, the analytical solution reads<br />

ˆp(x, ω) = i (Ŝ0<br />

4c 2 δ(y − y 1 )e 1)<br />

−iφ H0(kr 2 1 )dV 1 +<br />

i (Ŝ0<br />

0<br />

4c 2 δ(y − y 2 )e 2)<br />

−iφ H0(kr 2 2 )dV 2<br />

0<br />

+ i<br />

(4.8)<br />

(Ŝ0<br />

4c 2 δ(y − y 3 )e 3) −iφ H0(kr 2 3 )dV 3<br />

0<br />

where Ŝ 0 = 1E10 W/m 3 s; y 1 = [−0.015, 0.015] m, y 2 = [0.015, 0.015] m and y 3 = [0.0, −0.015]<br />

m; φ 1 = 0 rad, φ 2 = 1 rad and φ 3 = 2 rad; ω = 2π · 5000 rad/s. Fig. (4.8) shows the solution<br />

given by the co<strong>de</strong> AVSP-f. The reconstruction of the acoustic field at t 0 = 0 is evaluated by<br />

p ′ = | ˆp| cos(arg( ˆp)). This acoustic field is shown in Fig. (4.9). Fig. (4.10) shows that a good<br />

agreement is obtained between both analytical and numerical results. The directivity pattern<br />

produced by these three monopoles is satisfactorily well recovered.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!