THESE de DOCTORAT - cerfacs
THESE de DOCTORAT - cerfacs
THESE de DOCTORAT - cerfacs
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
4.1 Fundamental validation cases 61<br />
150<br />
210<br />
120<br />
240<br />
90<br />
270<br />
0.008<br />
60<br />
0.006<br />
0.004<br />
0.002<br />
180 0<br />
300<br />
30<br />
330<br />
r = 0.05<br />
r = 0.15<br />
Pressure Fluctuation (Pa)<br />
0.03<br />
0.02<br />
0.01<br />
0<br />
−0.01<br />
−0.02<br />
Analytics<br />
AVSP−f<br />
−0.03<br />
−0.2 −0.1 0 0.1 0.2<br />
X axis (m)<br />
a) | ˆp| Directivity of a monopole 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.4: Analytical solution vs Numerical solution.<br />
4.1.2 A Dipole in free space<br />
Once pure radiation is tested in section (4.1.1), the influence on the acoustic field of another<br />
monopole is tested. If this monopole radiates sound out of phase with respect to the first<br />
monopole, a substraction of acoustic fluctuations will occur leading to a diminution of the<br />
acoustic amplitu<strong>de</strong> in some regions. If this phase shift is equal to π a perfect cancelation should<br />
occur in the plane of antisymmetry between the two monopoles. Two monopoles with such a<br />
shift phase conform a dipole. The source is <strong>de</strong>fined as<br />
Ŝ(y, ω) = Ŝ 0 δ(y − y 1 )e −iφ 1<br />
+ Ŝ 0 δ(y − y 2 )e −iφ 2<br />
(4.5)<br />
where φ 1 = φ 2 + π. From Eq. (4.3)<br />
ˆp(x, ω) =<br />
i ∫<br />
4c 2 0<br />
V 0 (y)<br />
(Ŝ0 δ(y − y 1 )e −iφ 1<br />
+ Ŝ 0 δ(y − y 2 )e 2)<br />
−iφ H0(kr)dV 2 (4.6)<br />
Finally<br />
ˆp(x, ω) =<br />
i (Ŝ0 δ(y − y 1 )e 1)<br />
−iφ H0(kr 2 1 )dV 1 +<br />
4c 2 0<br />
i<br />
4c 2 0<br />
(Ŝ0 δ(y − y 2 )e 2)<br />
−iφ H0(kr 2 2 )dV 2 (4.7)