THESE de DOCTORAT - cerfacs
THESE de DOCTORAT - cerfacs
THESE de DOCTORAT - cerfacs
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.