22.01.2014 Views

user's manual for corhyd: an internal diffuser hydraulics model - IfH

user's manual for corhyd: an internal diffuser hydraulics model - IfH

user's manual for corhyd: an internal diffuser hydraulics model - IfH

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.

Gradual contraction (Idelchik 1986)<br />

4<br />

3<br />

2<br />

( − .0125⋅<br />

n + 0.0224⋅<br />

n − 0.00723⋅<br />

n + 0.0044⋅<br />

n − 0.00745)<br />

3 2<br />

ζ<br />

c<br />

= 0<br />

0<br />

0<br />

0<br />

0<br />

⋅ ( β − 2πβ<br />

−10β)<br />

with A0<br />

<strong>an</strong>d β in rad<br />

n = 1.0 A<br />

0<br />

≤<br />

1<br />

For β > 50°, the <strong>for</strong>mulation <strong>for</strong> gradual exp<strong>an</strong>sion leads to a greater loss coefficient th<strong>an</strong> the one<br />

<strong>for</strong> a sudden exp<strong>an</strong>sion. There<strong>for</strong>e Idelchiks <strong>for</strong>mulas was adopted so that <strong>for</strong> β > 50° losses are<br />

equal the loss <strong>for</strong> β = 50°.<br />

Code (see files: CommonFeederPipe.m, feederpipes.m, DiffuserLosses.m,<br />

Losses_common_feeder.m).<br />

Bending<br />

(reference<br />

velocity =<br />

velocity after<br />

bending<br />

Bend (Kalide 1980)<br />

3.5<br />

⎡<br />

⎛ D ⎞ ⎤ δ<br />

ζ0<br />

= ⎢0.131+<br />

0.159⎜<br />

⎟ ⎥ ⋅<br />

⎢⎣<br />

⎝ R ⎠ ⎥⎦<br />

180°<br />

where D is the pipe diameter <strong>an</strong>d R the radius of the bend. Often applied as R = 3D. Delta is the<br />

<strong>an</strong>gle of the bend (e.g. 90° <strong>for</strong> rect<strong>an</strong>gular bends).<br />

Code (see files: CommonFeederPipe.m, feederpipes.m, DiffuserLosses.m,<br />

Losses_common_feeder.m)<br />

Division<br />

flow<br />

of<br />

Friction due to bend (Idelchik 1986)<br />

L<br />

ζ<br />

fr<br />

= λ with<br />

L δ R<br />

= π<br />

D D 180°<br />

D<br />

(Idelchik 1986)<br />

∆p<br />

ζ<br />

s<br />

c,<br />

s<br />

ζ<br />

s<br />

= =<br />

2<br />

ρV<br />

/ 2 ( / ) 2<br />

s<br />

Vs<br />

Vc<br />

∆p<br />

ζ<br />

st<br />

c,st<br />

ζ<br />

st<br />

= =<br />

2<br />

ρV<br />

( ) 2<br />

st<br />

/ 2 Vst<br />

/ Vc<br />

ζ<br />

c,s<br />

from Diagram 7.15, ζ<br />

c, st<br />

from Diagram 7.17 (Idelchik, 1986 or Annex chapter 10). Curves<br />

fitted by the following code:<br />

Determination of zeta' (in the following zeta double underline) c,s<br />

vRatio = (q(i)/Ar(i)) / ((sum_q(i-1)+q(i))/Ad(i));<br />

if Dr(i)/Dd(i)

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

Saved successfully!

Ooh no, something went wrong!