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.
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