22.10.2014 Views

THESE de DOCTORAT - cerfacs

THESE de DOCTORAT - cerfacs

THESE de DOCTORAT - cerfacs

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

6.6 Transmitted and Reflected Waves through an i<strong>de</strong>al Compressor: the enthalpy jump case 113<br />

π c = p t,2 /p t,1 where the indices [1] and [2] represent the flow upstream and downstream of the<br />

compressor respectively. It is then useful to rewrite Eq. (6.13) as function of the parameter π c .<br />

F = (ρu 2 )<br />

∣<br />

2<br />

1<br />

+ p<br />

∣<br />

2<br />

1<br />

∫ x2<br />

+<br />

x 1<br />

∫ x2<br />

F = ∆p t + ∆e k +<br />

x 1<br />

∂<br />

(ρu) dx (6.49)<br />

∂t<br />

∂<br />

(ρu) dx (6.50)<br />

∂t<br />

where p t = p + ρu 2 /2 and e k = ρu 2 /2. From Eq. (6.50) it is observed that the compressor F<br />

acts as an element that exerts a change in the total pressure p t of the flow as well as in its kinetic<br />

energy e k . The parameter π c is now inserted into the expression leading to<br />

∫ x2<br />

F = p t1 (π c − 1) + ∆e k +<br />

x 1<br />

∂<br />

(ρu) dx (6.51)<br />

∂t<br />

Equation (6.51) is now linearized and expressed in the frequency domain. It yields<br />

∫ x2<br />

̂F = ¯p t ˆπ c + ˆp t,1 ( ¯π c − 1) + ∆ê k − iω ( ¯ρû + û ˆρ) dx (6.52)<br />

x 1<br />

where ê k = ˆρū 2 /2 + ¯ρūû and ˆp t = ˆp + ˆρū 2 /2 + ¯ρūû. For a compact system the integral can be<br />

neglected and Eq. (6.52) results in<br />

ˆ F = ¯p t ˆπ c + ˆp t,1 ( ¯π c − 1) + ∆ê k (6.53)<br />

In Appendix A, it is shown that π ′ c is a cumbersome term that <strong>de</strong>pends mainly on the mean and<br />

fluctuation value of the enthalpy both upstream and downstream of the compressor. Neglecting<br />

this term might be an strong assumption. The conditions of a compressor such that π ′ c ≈ 0<br />

are discussed in Appendix A. Nevertheless, it is practical as a first approximation and it is useful<br />

to see the influence of the mean value of the total pressure ratio ¯π c . After this assumption,<br />

the momentum equation (Eq. 6.22) becomes<br />

ˆp t,1 ¯π c = ˆp t,2 (6.54)<br />

Finally, the system of equations to resolve is expressed as function of the Mach Number<br />

The continuity equation (Eq. 6.31) and the momentum equation (Eq. 6.54) read<br />

¯M.

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

Saved successfully!

Ooh no, something went wrong!