13.02.2013 Views

Mechanics of Fluids

Mechanics of Fluids

Mechanics of Fluids

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

Compressible flow in pipes <strong>of</strong> constant cross-section 531<br />

To account for the changes <strong>of</strong> fluid properties we must consider initially an<br />

infinitesimal length <strong>of</strong> the pipe. For adiabatic flow the relation involving the<br />

friction force may be integrated over a finite length <strong>of</strong> the pipe by introducing<br />

the continuity equation and the energy equation. For isothermal conditions<br />

integration is made possible by inserting the condition T = constant.<br />

Applying the steady-flow momentum equation to an element with crosssectional<br />

area A, perimeter P and length δl (Fig. 11.26), we have<br />

pA − (p + δp)A − τ 0Pδl − ρgAδl sin α = ρAuδu<br />

Substituting τ 0 = 1<br />

2 f ρu2 (eqn 7.3), δl sin α = δz and dividing by ρgA we<br />

obtain<br />

−δz = δp uδu f δl<br />

+ +<br />

ρg g A/P<br />

u 2<br />

2g<br />

(11.62)<br />

The final term corresponds to Darcy’s formula for head lost to friction.<br />

Equation 11.62 applies to any steady one-dimensional flow where A is constant,<br />

but, for integration, the relation between density and pressure must<br />

be known and this depends on the degree <strong>of</strong> heat transfer. We consider first<br />

the case <strong>of</strong> zero heat transfer.<br />

11.10.1 Adiabatic flow in a pipe: Fanno flow<br />

Before embarking on the integration <strong>of</strong> eqn 11.62 for adiabatic conditions it<br />

will be instructive to look at the problem in general thermodynamic terms.<br />

When no external work is done and changes <strong>of</strong> elevation may be neglected,<br />

the steady-flow energy equation (11.12) and the continuity equation m =<br />

ρAu together give<br />

Constant h0 = h + 1<br />

2 u2 = h + 1<br />

� �2 m<br />

(11.63)<br />

2 ρA<br />

For given values <strong>of</strong> m, A and the stagnation enthalpy h0, curves <strong>of</strong> h against<br />

ρ could be plotted from eqn 11.63. A more significant relation, however, is<br />

that between h and specific entropy s. For any pure substance, s, like h and<br />

ρ, is a function <strong>of</strong> state and so may be determined from values <strong>of</strong> h and ρ.<br />

In particular, for a perfect gas s − s1 = cv ln(h/h1)(ρ1/ρ) γ −1 . Starting from<br />

a specified state (point 1 on Fig. 11.27) the curve <strong>of</strong> h against s traces the<br />

Fig. 11.26

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

Saved successfully!

Ooh no, something went wrong!