22.10.2014 Views

THESE de DOCTORAT - cerfacs

THESE de DOCTORAT - cerfacs

THESE de DOCTORAT - cerfacs

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

162 Chapter B: Publications<br />

The temporal <strong>de</strong>rivative of the <strong>de</strong>nsity in Eq. 10 is combined with Eq. 14. After some algebra, the<br />

divergence of the velocity can therefore be expressed as:<br />

∂u j<br />

= 1 ˙q (15)<br />

∂x j c p K 0<br />

This velocity field is supposed to be composed only by hydrodynamics, due to the fact that in the low<br />

Mach number mo<strong>de</strong>l the acoustic wave length is infinitely long. One can state that the divergence of the<br />

fluctuating velocity is<br />

B. Finding the Acoustic Pressure<br />

Injecting Eq. 16 into Eq. 9 leads to<br />

− ∂ ( 1 ∂p ′ )<br />

ac<br />

= ∂ ( ∂u<br />

′ )<br />

i,LES<br />

− ˙q′<br />

∂x i ρ 0 ∂x i ∂t ∂x i c p K 0<br />

or in the frequency domain<br />

Finally, multiplying everywhere by γP 0<br />

∂u ′ j,hyd<br />

∂x j<br />

= 1<br />

c p K 0<br />

˙q ′ (16)<br />

(17)<br />

( )<br />

∂ 1 ∂ ˆp ac<br />

= iω ∂û i,LES<br />

− iω ˆ˙q<br />

(18)<br />

∂x i ρ 0 ∂x i ∂x i c p K 0<br />

(<br />

∂<br />

c 2 ∂ ˆp )<br />

ac<br />

∂x i ∂x i<br />

= iωγP 0<br />

∂û i,LES<br />

∂x i<br />

} {{ }<br />

T 1<br />

− iω γP 0ˆ˙q<br />

c p K 0<br />

} {{ }<br />

T 2<br />

It has been found by the author that the contribution of T 2 can be neglected. As a consequence, one can<br />

state that for a reactive flow the divergence of the hydrodynamic velocity field ( ∂u hyd<br />

∂x i<br />

) can be consi<strong>de</strong>red<br />

zero as typically done for non-reactive flows. Equation 19 simplifies<br />

F(ˆp ac ) =<br />

∂ (<br />

c 2 ∂ ˆp )<br />

ac<br />

∂û i,LES<br />

= iωγP 0 (20)<br />

∂x i ∂x i ∂x i<br />

(19)<br />

VI.<br />

LES Vs Hybrid Results<br />

Figures 7 and 8 shows the Sound Pressure Level (SPL) given by the solution of Eq. 20 (LES ‘Filtered’)<br />

and the hybrid computation for microphones 5, 6 and 7.<br />

SPL (dB) − micro 5<br />

180<br />

160<br />

140<br />

120<br />

100<br />

80<br />

Hybrid Computation<br />

Direct Computation: LES (Filtered)<br />

SPL (dB) − micro 6<br />

180<br />

160<br />

140<br />

120<br />

100<br />

80<br />

Hybrid Computation<br />

Direct Computation: LES (Filtered)<br />

60<br />

0 500 1000 1500 2000 2500<br />

Frequency (Hz)<br />

60<br />

0 500 1000 1500 2000 2500<br />

Frequency (Hz)<br />

Figure 7.<br />

Sound Pressure Levels from the direct and hybrid approaches<br />

After extracting the acoustic field from the complete pressure fluctuation field, it is seen that results match<br />

pretty well with the values predicted by the hybrid method, especially for microphones 6 and 7. There is a<br />

high contain of hydrodynamic fluctuations at low frequencies (before 400 Hz) that has been removed. Also<br />

7 of 9<br />

American Institute of Aeronautics and Astronautics

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

Saved successfully!

Ooh no, something went wrong!