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.

161<br />

Observing with attention fig. ??, it is noticeable that there are some zones of the spectrum in which<br />

an important gap is present between hybrid and direct computations. It is probable, consi<strong>de</strong>ring the direct<br />

computation, that the fluctuations of pressure coming from LES are composed by both acoustic and hydrodynamic<br />

contributions. On the other hand, presure fluctuations coming from the hybrid computation are<br />

totally due to acoustics. A suitable comparison is then not carried out, and it becomes important to be able<br />

to extract acoustics from LES computations in or<strong>de</strong>r to evaluate in a proper way the results from the hybrid<br />

computation.<br />

V. Filtering a LES pressure field to find the corresponding acoustic field<br />

Velocity fluctuations obtained by LES are composed by both hydrodynamics and acoustics<br />

Appling the operator ∂/∂t to Eq. 5 leads to<br />

u ′ i,LES = u ′ i,hyd + u ′ i,ac (5)<br />

∂u ′ i,LES<br />

= ∂u′ i,hyd<br />

+ ∂u′ i,ac<br />

∂t ∂t ∂t<br />

From linear acoustics, the momentum equation is given by<br />

(6)<br />

ρ 0<br />

∂u ′ i,ac<br />

∂t<br />

= − ∂p′ ac<br />

∂x i<br />

(7)<br />

where [] 0 and [] ′ represent respectively the mean and fluctuating flow. Combining this term of eq. 7 into<br />

Eq. 6<br />

− 1 ρ 0<br />

∂p ′ ac<br />

∂x i<br />

+ ∂u′ i,hyd<br />

∂t<br />

Finally the divergence operator to this equation is applied<br />

= ∂u′ i,LES<br />

∂t<br />

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

ac<br />

+ ∂ ( ∂u<br />

′ )<br />

i,hyd<br />

= ∂2 u ′ i,LES<br />

∂x i ρ 0 ∂x i ∂x i ∂t ∂x i ∂t<br />

(8)<br />

(9)<br />

A. Finding ∂u i,hyd<br />

∂x i<br />

Neglecting viscosity, species diffusion and heat conduction the Navier-Stokes equations for reacting flows<br />

read<br />

∂ρ<br />

∂t + ρ∂u j<br />

∂x j<br />

+ u j<br />

∂ρ<br />

∂x j<br />

= 0 (10)<br />

ρ ∂u i<br />

∂t + ρu ∂u i<br />

j = − ∂p<br />

(11)<br />

∂x j ∂x i<br />

ρc p<br />

∂T<br />

∂t + ρc pu j<br />

∂T<br />

∂x j<br />

= ˙q (12)<br />

In the low-Mach number approximation, the thermodynamic pressure P 0 only <strong>de</strong>pends on temperature.<br />

The equation of state is simply<br />

Replacing Eq. 13 in the left hand si<strong>de</strong> of the Eq. 12 leads to<br />

K 0 = ρT (13)<br />

ρc p<br />

∂K 0 /ρ<br />

∂t<br />

+ ρc p u j<br />

∂K 0 /ρ<br />

∂x j<br />

= ˙q (14)<br />

6 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!