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.
6.3 The 1D linearized Euler equations for compact systems 99<br />
∂ρ<br />
∂t + ∂ (ρu) = 0<br />
∂x<br />
(6.9)<br />
∂ρu<br />
+ ∂<br />
∂t ∂x (ρu2 ) = − ∂p<br />
∂x + Ϝ (6.10)<br />
∂ρh t<br />
+ ∂<br />
∂t ∂x (ρuh t) = ˙Q + W k (6.11)<br />
Integrating over the 1D domain between positions x 1 and x 2 results in<br />
∫ x2<br />
x 1<br />
∫<br />
∂ρ<br />
∂t dx + x2<br />
∫ x2<br />
x 1<br />
x 1<br />
∂<br />
∂x (ρu)dx = 0<br />
∂ρ<br />
∂t dx + (ρu) ∣ ∣∣∣<br />
x 2<br />
x 1<br />
= 0<br />
(6.12)<br />
∫ x2<br />
x 1<br />
∂ρu<br />
∂t dx + ∫ x2<br />
x 1<br />
∫ x2<br />
∫<br />
∂<br />
∂p<br />
∂x (ρu2 )dx +<br />
x 1<br />
∂x dx − x2<br />
Ϝdx = 0<br />
x 1<br />
∫ x2<br />
x<br />
∂ρu<br />
x 1<br />
∂t dx + 2<br />
x 2<br />
(ρu2 )<br />
∣ + p<br />
∣ − F = 0<br />
x 1 x 1<br />
(6.13)<br />
∫ x2<br />
x 1<br />
∂ρh t<br />
∂t dx + ∫ x2<br />
x 1<br />
∂<br />
∂x (ρuh t)dx −<br />
∫ x2<br />
x 1<br />
∫ x2<br />
x 1<br />
∂ρh t<br />
∂t dx + (ρuh t)<br />
∫ x2<br />
˙Qdx − W k dx = 0<br />
x 1<br />
x 2<br />
∣ − ˙q − W = 0<br />
x 1<br />
(6.14)<br />
where the following notations have been introduced: ∫ ˙Qdx = ˙q, ∫ Ϝdx = F and ∫ W k dx =<br />
W. After linearizing and recalling that h t = c p T t , Eqs. (6.12), (6.13) and (6.14) become<br />
( ¯ρu ′ + ρ ′ ū ) ∣ ∣ ∣∣<br />
x 2<br />
(<br />
p ′ + ρ ′ ū 2 + 2 ¯ρūu ′) ∣ ∣ ∣∣<br />
x 2<br />
[<br />
cp ¯T t ( ¯ρu ′ + ρ ′ ū) + ¯ρū(c p T ′ + ūu ′ ) ] ∣ ∣ ∣∣<br />
x 2<br />
x 1<br />
= −<br />
x 1<br />
− F ′ = −<br />
x 1<br />
− ˙q ′ − W ′ = −<br />
∫ x2<br />
x 1<br />
∂ρ ′<br />
∂t<br />
∫ x2<br />
x 1<br />
∫ x2<br />
x 1<br />
dx (6.15)<br />
∂ ( ¯ρu ′ + ρ ′ ū ) dx (6.16)<br />
∂t<br />
∂ [<br />
ρ ′ c p ¯T t + ¯ρ ( c p T ′ + ūu ′)] dx<br />
∂t<br />
(6.17)