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
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)