13.02.2013 Views

Mechanics of Fluids

Mechanics of Fluids

Mechanics of Fluids

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.

12.3.4 The method <strong>of</strong> characteristics<br />

Since the velocity and pressure in a pipe subject to pressure transients are<br />

continuous functions <strong>of</strong> position and time, they are described essentially<br />

by partial differential equations. The method <strong>of</strong> characteristics is a general<br />

technique in which such partial differential equations are converted into particular<br />

simultaneous total differential equations which, after being expressed<br />

in finite-difference form, can be solved by computer. The effects <strong>of</strong> friction<br />

and difference <strong>of</strong> elevation can be retained in the equations and complex<br />

pipe system can be dealt with. This method is not the only one possible for<br />

solving the equations but it is increasingly used for the study <strong>of</strong> problems<br />

involving pressure transients, so we shall here derive the finite-difference<br />

equations and outline the method to the point where a computer program<br />

can be written.<br />

Pressure transients 575<br />

Figure 12.9 shows a short length δx <strong>of</strong> the pipe through which a pressure The equation <strong>of</strong> motion<br />

wave travels upstream. The velocity <strong>of</strong> the fluid is assumed uniform over<br />

each cross-section and so the fluid in the length δx may be regarded as a<br />

single particle to which Newton’s Second Law can be directly applied. Then<br />

the force in the direction <strong>of</strong> flow<br />

�<br />

pA − p + ∂p<br />

∂x δx<br />

��<br />

A + ∂A<br />

∂x δx<br />

�<br />

−τ0Pδx − mg sin α<br />

= Mass × Acceleration = m<br />

�<br />

+ p + 1<br />

2<br />

�<br />

u ∂u<br />

�<br />

∂u<br />

+<br />

∂x ∂t<br />

∂p<br />

∂x δx<br />

�<br />

∂A<br />

∂x δx<br />

(12.14)<br />

Here the third term on the left-hand side consists <strong>of</strong> the average pressure at<br />

the sides <strong>of</strong> the element multiplied by the projected area <strong>of</strong> the sides perpendicular<br />

to the flow direction. The term therefore represents the component in<br />

the flow direction <strong>of</strong> the force exerted by the pipe walls on the fluid particle.<br />

Fig. 12.9

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

Saved successfully!

Ooh no, something went wrong!