Handbook For Control Valve Sizing: Bulletin 1-I
Handbook For Control Valve Sizing: Bulletin 1-I
Handbook For Control Valve Sizing: Bulletin 1-I
C v1 , F L1
Fpv
(p F )
Cv1= Cv 2
Fpv
(p F )
L 1
L2
pm ax (1) = F
p1
pmax (2) =F
C v 2 ,F L 2
FL1 > FL 2
pv
pvc
p2I
p2II
p2III
p2IV
HANDBOOK
FOR CONTROL
VALVE SIZING
* PARCOL
HANDBOOK
FOR CONTROL
VALVE SIZING
CONTENTS
NOMENCLATURE
VALVE SIZING AND SELECTION
1
PROCESS DATA
VALVE SPECIFICATION
FLOW COEFFICIENT
3.1 KV coefficient
3.2 Cv coefficient
3.3 Standard test conditions
SIZING EQUATIONS
4.1 Sizing equations for incompressible fluids (turbulent flow)
4.2 Sizing equations for compressible fluids
(turbulent flow)
4.3 Sizing equations for two-phase fluids
4.4 Sizing equations for non turbulent flow
-1-
* PARCOL
Symbols
Description
Units (note)
Cd
various
Cv
Flow coefficient
U.S. gallons/min
mm
mm
Fd
dimensionless
FF
dimensionless
FL
dimensionless
FP
dimensionless
FR
dimensionless
dimensionless
dimensionless
Kc
dimensionless
Kv
Flow coefficient
m 3/h
K1 and K2
dimensionless
kg/kmole
pc
bar
pv
bar
pvc
bar
p1
bar
p2
bar
bar
p max
bar
qm
kg/h
qv
m 3/h
qm(max)
kg/h
qv(max)
m 3/h
Rev
dimensionless
T1
m/s
FLP
dimensionless
-2-
* PARCOL
Symbols
Description
Units
Specific volume
m 3/kg
dimensionless
xcr
dimensionless
xFZ
dimensionless
xT
dimensionless
xTP
Expansion factor
dimensionless
dimensionless
dimensionless
kg/m 3
1
r
Kinematic viscosity ( = /)
Centistoke = 10 -6 m 2/s
Dynamic viscosity
Centipoise = 10 -3 Pa s
-3-
dimensionless
* PARCOL
SIZING AND SELECTION OF CONTROL
VALVES
- VALVE SPECIFICATION
- valve rating
- body and valve type
- body size, after having calculated the maximum flow coefficient Cv with the appropriate
sizing equations
- type of trim
- materials trim of different trim parts
- leakage class
- inherent flow characteristic
- packing type
- type and size of actuator
- accessories
- PROCESS DATA
The following data should at least be known:
a - Type of fluid and its chemical-physical and
thermodynamic characteristics, such as pressure p, temperature T, vapour pressure
pv, thermodynamic critical pressure pc,
specific mass , kinematic viscosity or
dynamic viscosity , specific heat at constant pressure Cp, specific heat at constant
volume Cv , specific heat ratio , molecular mass M, compressibility factor Z, ratio
of vapour to its liquid, presence of solid particles, inflammability, toxicity.
- FLOW COEFFICIENT
K v = qv
p( Kv )
(1)
where:
p(Kv) is the static pressure drop of 105 Pa
p is the static pressure drop from upstream to
downstream in Pa
is the specific mass of fluid in kg/m3
o is the specific mass of water in kg/m3
-4-
* PARCOL
4
C v = qv
p( Cv )
(2)
4.1 - SIZING EQUATIONS FOR INCOMPRESSIBLE
FLUIDS (TURBULENT FLOW)
where:
p(Cv) is the static pressure drop of 1 psi (see
above)
p is the static pressure drop from upstream to
downstream expressed in psi.
is the specific mass of the fluid expressed in
Ib/ft3
o is the specific mass of the water expressed in
Ib/ft3
6D
p1
20D (*)
- SIZING EQUATIONS
p2
(*)
10D (*)
-5-
* PARCOL
approximation is justified by simplicity purposes
and by the fact that it is not practically important
to predict the exact flow rate in the hatched area;
on the other hand such an area should be
avoided, when possible, as it always involves
vibration and noise problems as well as mechanical problems due to cavitation.
Basic equation
Valid for standard test conditions only.
qv = K v
p
/
with qv in m3/s
p in bar (105 Pa)
p
/
qv = Cv
with qv in gpm
p in psi
CV =
CV =
p <
F 2
pmax = LP
Fp
(p F pv )
1
qm
865 FRp p r
1.16 q v
FpR
p
r
qm
p max = FL p1 FF p v
IEC limit flow
approximation of IEC
equations
2%
normal flow
semi-critical flow
p =
K c (p1 pv ) flow
f rate affected by cavitation
beginning of cavitation
p = x FZ (p1 p v )
p
Fig.2 -Flow rate diagram of an incompressible fluid flowing through a valve plotted versus downstream
pressure under constant upstream conditions.
-6-
* PARCOL
Limit flow
It is individuated by the relationship:
Cv =
Cv =
q m (max )
865 FLP
2
2
FLPP
(p1 FF pv )
p pmax =
F
p
(p1 FF p v )r
1.16 q v (max )
FLP
(p1 FFp v )
r
The Fig. 3 shows the flow rate diagram of a compressible fluid flowing through a valve when
changing the downstream pressure under constant upstream conditions. The flow rate is no
longer proportional to the square root of the pressure differential p as in the case of incompressible fluids. This deviation from linearity is
due to the variation of fluid density (expansion)
from the valve inlet up to the vena contracta.
Normal flow
It is individuated by the relationship
x < F xT
limit flow
-7-
2/3 < Y 1
Cv =
qm
27.3 Fp Y x p1 1
Cv =
qv
M T1 Z
2120 Fp p1 Y
x
limit flow
density variation effect
or
* PARCOL
Limit flow
It is individuated by the relationship
x F xTP
Cv =
Cv =
and/or
A second physical model overcomes this limitation assuming that the two phases cross the vena
contracta at the same velocity.
Y = 2/3 = 0.667
The mass flow rate of a gas (see above) is proportional to:
q m( max )
18.2 Fp F x TP p1 1
Y x 1 = Y
x
= x/ V eg
V g1
q v (max )
M T1 Z
1414 Fp p1
F x TP
V g1 /Y 2
In other terms this means to assume that the
mass flow of a gas with specific volume Vg1 is
equivalent to the mass flow of a liquid with specific volume Veg under the same operating conditions.
Assuming :
Ve = f g
Vgl
Y
+ f liq Vliq1
q m = 27. 3 Fp C v
x p1
Ve
Therefore:
C v = C v g + C v liq
-8-
* PARCOL
The most reliable explanation of such results is
that the two phases flow at quite different velocities, though mutually exchanging mass and
energy.
C v = C v liq + C v vap
Cv =
- experimental data are available on liquid-vapour mixtures flowing through orifices at flow
rates 1012 times higher than the ones resulting from calculation when considering the fluid
as compressible with a specific mass equal to
the one at the valve inlet.
qm
27 .3 Fp
x p1
Ve
Temperature
1
Fig. 4 - Thermodynamic transformations of a water / vapour
mixture inside a valve.
1
2
Vc
2
Vc
Enthropy
-9-
* PARCOL
The effect of fittings attached to the valve is
probably negligible in laminar flow condition and
it is presently unknown.
In equations applicable to compressible fluid the
correcting factor p 1+p 2/2 was introduced to
account for the fluid density change.
Re = u d
CV =
qm
865 FR p r
incompressible fluid
CV =
1.16 q v
FR
CV =
p
r
qm
T1
67 FR
p ( p 1 + p 2 ) M
FL =
compressible fluid
CV =
qv
M T1
1500 FR
p ( p 1 + p 2 )
p1 p2
p1 pvc
- 10 -
* PARCOL
longer and this flow rate is assumed as q v(max).
- Accuracy in determination of FL
It is relatively easier determining the critical flow
rate qv(max) for high recovery valves (low FL) than
for low recovery valves (high FL ). The accuracy
in the determination of FL for values higher than
0.9 is not so impor tant for the calculation of the
flow capacity as to enable to correctly predict
the cavitation phenomenon for services with high
differential pressure.
- Determination of FL
c
FL =
1.16q v (max )
C v p1 0 .96 p v
C v1 , F L 1
Fpv
(p F )
Cv1= Cv 2
Fpv
(p F )
L 1
pm ax (1) = F
p1
L2
pmax (2) = F
C v2 ,F L 2
FL1 >FL 2
Cv1FL1
Cv2FL2
pv
pvc
- 11 -
p2 I
p2II Fig. 5 - Comparison between two valves with
equal flow coefficient but
with different recovery facIII
p2 tor, under the same inlet
fluid condition, when varyIV
p2 ing the downstream pressure. At the same values
of Cv, p1 and p2 valves
with higher F L can accept
higher flow rates of fluid.
* PARCOL
1
1
0.95
4
2, 4
Coefficiente di recupero FL
Pressure recovery factor FL
0.9
7
0.85
5
3
0.8
0.75
0.7
7
0.65
3
1
2
3
4
5
6
7
0.6
10
50
1-8110
1-6911
1-6911
1-6933
1-6600
1-2471
1-8110
% del Cv max
% of rated Cv
100
Fig. 6 - Typical FL values versus % value Cv and flow direction for different PARCOL valve types.
5.2 - COEFFICIENT OF INCIPIENT CAVITATION
XFZ AND COEFFICIENT OF CONSTANT CAVITATION Kc
Index of resistance to cavitation
stellite gr. 6
chrome plating
17-4-PH H900
AISI 316/304
monel 400
gray cast iron
chrome-molybdenum alloyed steels (5% chrome)
carbon steels (WCB)
bronze (B16)
nickel plating
pure aluminium
20
(5)
2
1
(0.8)
0.75
0.67
0.38
0.08
(0.07)
0.006
- 12 -
* PARCOL
the formula Kc = 0.8 FL2. Such a simplification
is however only acceptable when the diagram of
the actual flow rate versus p , under constant
upstream conditions, shows a sharp break point
between the linear/proportional zone and the
horizontal one. If on the contrary the break point
radius is larger (i.e. if the p at which the deviation from the linearity takes place is different from
the p at which the limit flow rate is reached)
then the coefficient of proportionality between
Kc and FL2 can come down to 0.65. Since the
coefficient of constant cavitation changes with
the valve opening, it is usually referred to a 75%
opening.
Critical conditions are obviously reached gradually. Moreover the velocity profile in the vena
contracta is not completely uniform, hence may
be that a part only of the flow reaches the vaporization pressure. The FL recovery factor is
determined in proximity of fully critical conditions,
so it is not suitable to predict an absolute absence of vaporization. In order to detect the beginning of the constant bubble formation, i.e. the
constant cavitation, the coefficient Kc was defined. This coefficient is defined as the ratio p/
(p1 - pv ) at which cavitation begins to appear in
a water flow through the valve with such an intensity that, under constant upstream conditions,
the flow rate deviation from the linearity versus
p exceeds 2%. Usually the beginning of cavitation is identified by the coefficient of incipient
cavitation xFZ. The xFZ coefficient can be determined by test using sound level meters or accelerometers connected to the pipe and relating
noise and vibration increase with the beginning
of bubble formation. Some informations on this
regard are given by standard IEC 534-8-2 Laboratory measurement of the noise generated by a
liquid flow through a control valve, which the Fig.
8 was drawn from. A simple calculation rule uses
p/(p1-pv)
Fp =
XFZ
ptr
x FZ =
p1 pv
being:
1
K C v
1+
0.00214 d 2
K = K 1+ K 2 + K B1 K B2
- 13 -
* PARCOL
In case of reducers:
d 2
K1 = 0.5 1
D
it becomes:
2
FLp
(p1 FF pv )
Fp
In case of expanders:
d 2
=
K 2 1.01
D
(see Fig. 9)
It is determined by test, like for the recovery factor FL (see point 5.1).
FLP =
d 2
K1 + K 2 = 1.51
D
FLP =
5.4 -
(p1 F F p v ) , but
FL
FL2
(K)1 C2v
1+
0.00214
d
p max = FL p1 FF pv
q max FL Cv
q max FLP Cv
p1 FF pv
pC
v
FLP
Fp
q
F
q
C
p max =
p
Fig. 9 - Effect of reducers on the diagram of q versus
- 14 -
* PARCOL
dition and the vapour pressure of the liquid at
inlet temperature:
FF = pvc /pv
When the flow is at limit conditions (saturation)
the flow rate equation must no longer be expressed as a function of p = p1-p2, but of pvc
= p1 -pvc(differential pressure in vena contracta).
Starting from the basic equation (at point 4.1):
qv = Cv
q mY x
A mathematic procedure allows to calculate the
value of Y which makes maximum the above
function (that means finding the point where the
rate dqm/ dx becomes zero.
p1 p 2
r
and from:
FL =
3
q m (1 ax ) x = x a x
p1 p2
p1 pvc
By setting
dq m
1
3a x
=
=0
dx
2
2 x
q v = FL C v
p1 p vc
r
1
= 3a x
x
q v = FL C v
i.e.:
pv
pc
1
3a
2
1
a =
3a
3
Y =1
pv
pc
x
3x T
thus taking into account also the third hypothesis. As a matter of fact xT is an experimental
value to be determined for each valve type. Finally the second hypothesis will be taken into
account with an appropriate correction factor:
x=
FF = 0.96 0.28
Y = 1
hence:
EXPANSION FACTOR Y
F = /1.4, which is the ratio between the exponent of the adiabatic transformation for the actual gas and the one for air.
Y = 1
- 15 -
x
3Fy x T
* PARCOL
Cv/d2
(d in mm)
15 x 10-3
20 x 10-3
25 x 10-3
30 x 10-3
35 x 10-3
40 x 10-3
FL
.5 .6 .7 .8 .9
.5 .6 .7 .8 .9
.5 .6 .7 .8 .9
.5 .6 .7 .8 .9
.5 .6 .7 .8 .9
.5 .6 .7 .8 .9
d/D
FLP
FLP
FLP
FLP
FLP
FLP
.25
.33
.40
.50
.66
.75
Fig. 10 - Values of FLP for valves with short type reducer at the inlet with abrupt section variation
1
0.96
FF
0.9
0.8
0.7
FF = 0.96 0.28
0.6
0
0.1
0.2
0.3
0.4
0.5
0.6
0.7
0.8
0.9
pv
pc
Pv
Pc
1
0,96
0,95
0,90
0,85
0,80
221.2
0,75
0,70
0,68
0,65
0,60
pv
FF = 0.96 0.28
221.2
50
100
150
pv = Vapour pressure (bar abs.)
200
pc
250
1
Y
0.9
0.8
XT
ing
eas
incr
0.7
0.667
x=
0.6
0.2
0.4
0.6
0.8
- 16 -
1.0
p
p1
* PARCOL
Therefore the maximum flow rate is reached
when x = F . xT (or F xTP if the valve is supplied with reducers) ; correspondently the expansion factor reaches the minimum value of 0.667.
5.7 -
p
x cr =
p1 cr
xTP is the same coefficients xT however determined on valves supplied with reducers or installed not in according to the standard set up.
y 1
x cr = F 1
+ 1
2
L
Cd
xT
d/D
.80
.75
.67
.60
.50
.40
.33
.25
10
15
xTP
.40 .49 .59 .69 .78
.40 .50 .59 .69 .78
.40 .50 .60 .69 .78
.41 .51 .60 .70 .79
.41 .52 .61 .70 .80
.42 .52 .62 .71 .80
.43 .53 .62 .72 .81
.44 .53 .63 .73 .83
x TP =
20
Fp
.99
.98
.98
.97
.96
.95
.94
.93
xTP
.40 .49 .58 .67 .75
.40 .49 .58 .67 .75
.41 .50 .59 .68 .76
.42 .52 .61 .69 .78
.44 .53 .63 .71 .79
.44 .55 .65 .74 .82
.46 .56 .66 .75 .83
.48 .58 .67 .76 .85
xTP
.39 .48 .56 .64
.40 .49 .57 .65
.42 .51 .59 .67
.43 .53 .61 .69
.46 .55 .64 .72
.49 .58 .67 .75
.50 .60 .69 .78
.52 .62 .71 .79
1+
x T (K1 + K B1 ) C v
2
0.0024
d
25
Fp
.98
.97
.95
.93
.91
.89
.88
.87
xT
(Fp )2
30
Fp
.96
.94
.91
.89
.85
.82
.81
.79
xTP
.21 .30 .39 .47
.22 .31 .40 .48
.24 .33 .43 .51
.25 .36 .45 .54
.28 .39 .49 .58
.30 .42 .53 .62
.31 .44 .55 .64
.33 .46 .57 .67
.15
Fp
.94
.91
.87
.84
.79
.76
.74
.72
.20
.25
xTP
.17
.18
.19
.21
.24
.26
.27
.27
.21
.23
.25
.27
.30
.33
.34
.37
.26
.27
.30
.32
.36
.40
.40
.44
Fp
.91
.88
.83
.79
.73
.70
.69
.65
Fig. 14 -Calculated values of xTP and Fp for valves installed between two commercial concentric reducers (with abrupt section variation)
Cd = Cv / d 2 (d expressed in inches).
Example: For a 2" valve is: Cv = 80 and xT = 0.65
The valve is installed in a 3" pipe between two short type reducers.
Cd = Cv / d2 = 20 d / D = 2/3 = 0.67
A linear interpolation between xT = 0.6 and xT = 0.7 results in xTP = 0.63
- 17 -
* PARCOL
FR
Cd=10
Cd=15
5.9 -
Cd=20
laminar flow
transitional flow
turbulent
Rev
Flow direc
ectti on
Max. opening
90
60
Cage valve
1-6933, 1-4433,
1-6971, 1-4471
Number of holes
50
100
200
Double seat
1-8110
Parabolic
V-port
1. 00
Flow-to-open
Flow-to-close
0.10
0.20
0.46
1.00
Whatever
0.20
0.20
0.7
0.5
Whatever
0.45
0.32
0.22
0.14
0.10
0.07
Between
seats
0.10
0.10
0.32
0.28
Fig. 16 -Typical Fd values for PARCOL control valves. More accurate values on request
- 18 -
* PARCOL
This data sheet was derived from IEC 60534-7 with some improvements not affecting the numbering of the original items.
- 19 -
Telephone: +39 0331 413 111 - Fax: +39 0331 404 215
e-mail: [email protected] - http://www.parcol.com
1120047 Studio Trevisan - Gallarate -