Pin Unit 2
Pin Unit 2
VUT
Vaal University of Technology
2/10
EIPINI Chapter 4: Level Measurement Page 4-1
4. LEVEL MEASUREMENT
This chapter aims to introduce students to some of the level measurement technologies
used to measure the level of liquids and granulars in a container.
4.1 INTRODUCTION
On the 28th of March, 1979, thousands of people fled from Three Mile Island when
the cooling system of a nuclear reactor failed. This dangerous situation developed
because the level controls turned off the coolant flow to the reactor when they detected
the presence of cooling water near the top of the tank. Unfortunately, the water
reached the top of the reactor vessel not because there was too much water in the tank,
but because there was so little that it boiled and swelled to the top. From this example,
we can see that level measurement is more complex than simply the determination of
the presence or absence of a fluid at a particular elevation.
Most level measurement techniques can be categorized into one of two groups –
direct and indirect (inferred) methods. Direct methods involve measuring the height of
fluid directly as for example with a dipstick, overflow pipe, float or sight glass. With
indirect methods, another variable is measured that correlate to the liquid level.
Measuring the weight of a substance in a container or the pressure exerted on the
bottom of a tank or transmitting an ultrasonic beam to the level surface and measuring
the time of flight of the transmitted and received signal, is indirectly related to the
level height.
Level measurement may also be continuous or “point-level”. A continuous method
refers to a technique whereby the device measures level on a constant basis, displaying
or transmitting the actual level of the liquid as it changes. Point-level devices measure
liquid at specific points within the tank. As the liquid level rises and falls, it passes
through definite points during it's transit. Continuous methods may however be
programmed to output alarms at specific points also, in addition to transmitting a
continuous level measurement.
Level measuring devices are also described as either a contact or non-contact type of
instrument. A contact type device, such as a float, makes physical contact with the
liquid in the container, in order to determine the level. A non-contact device, such as
ultrasonic or radar, does not require contact with the material in the container to
measure the level.
4.2 THE SIGHT GLASS
The sight glass consists of a strengthened glass tube, attached to the container as
shown in Figure 4-1, through which the fluid level in the container can easily be
observed by the operator. Monitoring the level from a distance, is facilitated when a
magnetic float inside the tube, is allowed to rise and fall with the liquid level, causing
metallic flags (or louvers) to flip and expose a different colour, indicating the level.
These devices are available with continuous monitoring equipment, allowing both a
local visible indication as well as an external signal for remote monitoring.
EIPINI Chapter 4: Level Measurement Page 4-2
Valve
Sight glass
Magnetic
float
Metallic
flaps
Valve
Container with sight glass Sight glass with magnetic indicator
Figure 4-1
4.3.3 Magnetic Float (Magnetic Coupled Float and Follower) Level Meter
Torque tube
Archimedes’s Principle
Torque tube
flange
Level indicator
Torque rod
Torque arm
Chain
Displacer
Figure 4-6
Example 4-1
During calibration of a displacer type level meter, it is found that the torque registered
by the meter when the tank is empty, is 10 N-m. If the torque arm is 0.1 m, calculate
the torque that will be generated when the displacer, with volume 0.002 m3, is fully
immersed in the process fluid with density 1000 kg/m3.
Remembering that torque is given by T = Fr, the weight of the displacer in air is
10/0.1 = 100 N. Using Archimedes’ law, the apparent weight of the displacer when it
is completely immersed in the liquid, is 100 – 0.002×1000×9.81 = 80.38 N. The torque
generated then, will be 80.38×0.1 = 8.038 N-m.
EIPINI Chapter 4: Level Measurement Page 4-5
Perhaps the most frequently used technique of measuring level, is the method of
measuring differential pressure. The primary benefit of this method is that the
equipment measuring the differential pressure, can be externally installed or
retrofitted to an existing vessel. It can also be isolated safely from the process
using block valves for maintenance and testing.
There are, however, some disadvantages to this method. One vessel penetrations
near the bottom of the vessel is needed, where leak paths could be the cause of
many problems. Measurement errors also occur due to changes in liquid density.
Density variations are caused by temperature changes or change of product. These
variations must always be compensated for, to maintain accurate measurements.
H Patm H
z HP LP z HP LP
Zero line h
DP
transmitter
U tube
manometer
Example 4-2
A mercury (density 13600 Patm Patm
kg/m3) u-tube manometer is
used to measure the level of a
liquid (density 1000 kg/m3) in 5m
a 5 meter high container, with 3m
the zero level of the
manometer, exactly in line h Zero line
with the bottom of the tank.
Calculate the manometer X Y
reading h, if the tank level is 3
meter.
Equating pressures on the XY line: U-tube
P X = PY
∴Patm + 1000×(3 + ½h)×9.81 = Patm + 13600×h×9.81
∴3 + ½h = 13.6h
∴13.1h = 3
∴h = 0.229 m
= 229 mm
We will now use method a), b) and c) to solve this problem, but we will prefer the
approximate but easier method c) for the remaining well type manometer problems
and also for closed tanks involving well type manometers.
a) Comparing pressures on the XY line:
PATM + 1000×3×9.81 + 1000×d×9.81 = PATM + 13600×h×9.81 + 13600×d×9.81
∴1000×3×9.81 + 1000×d×9.81 = 13600×h×9.81 + 13600×d×9.81
∴29430 + 9810d = 133416h + 133416d
∴29430 = 133416h + 123606d………………………………………………..(a)
But from Equation (2), Chapter 2, Page 2-7: d = (A2/A1)h
∴d = 0.01h……………………………………………………………………..(b)
(b) in (a):
29430 = 133416h + 123606×(0.01h)
∴29430 = 133416h + 1236.06h
∴134652.06h = 29430
∴h = 0.21856 m
b) P1 = PATM + 1000×3×9.81 + 1000×d×9.81
= PATM + 29430 + 9810×d = PATM + 29430 + 9810×(0.01h)
∴P1 = PATM + 29430 + 98.1h …………………………………….…………… (a)
And P2 = PATM ………………………………………………………………… (b)
From (a) and (b):
P1 – P2 = 29430 + 98.1h ………………………………………………………. (c)
But from Equation 2-5 (Chapter 2):
P1 – P2 = ρhg(1 + A2/A1) …………………………………...…………………. (d)
(c) in (d):
29430 + 98.1h = 13600×h×9.81(1 + 0.01)
∴29430 + 98.1h = 133416×h×1.01
∴29430 + 98.1h = 134750.16h
∴134652.06h = 29430
∴h = 0.21856 m
c) P1 = PATM + 1000×3×9.81 (neglecting d)
∴P1 = PATM + 29430 ………………………………………….…………….… (a)
And P2 = PATM …………………………………………………….……...…… (b)
From (a) and (b):
P1 – P2 = 29430 ………………………………………………….……………. (c)
And from Equation 2-5 (Chapter 2):
P1 – P2 = ρhg(1 + A2/A1) …………………………………...…………………. (d)
(c) in (d):
29430 = 13600×h×9.81(1 + 0.01)
∴29430 = 134750.16h
∴h = 0.2184 m
EIPINI Chapter 4: Level Measurement Page 4-8
Example 4-4
A DP transmitter must be calibrated to measure the level of a liquid in an open tank.
The density of the liquid is 1000 kg/m3. The DP transmitter will be mounted one meter
below the bottom of the tank. The tank is full when the height of the liquid in the tank
is 5 meter and it is empty when there is only liquid in the high pressure line connected
to the DP transmitter. Determine the necessary calibration specifications for an output
signal of 4 to 20 mA.
Example 4-5
If, in example 4-4, the DP cell is replaced by a Empty Patm
Patm
u-tube manometer using mercury (density 13600
kg/m3), with its zero level 1 meter below the
bottom of the container, calculate the full and
empty readings, hempty and hfull. 1m
ZL hempty
Empty: X Y
Equating pressures on the XY line:
Patm+1000×(1+½h)×9.81 = Patm+13600×h×9.81 U tube
∴1 + ½h = 13.6h
∴hempty = 0.07634 m = 76.34 mm. Full Patm
Patm
Full: 5m
Equating pressures on the XY line:
Patm+1000×(6+½h)×9.81 = Patm+13600×h×9.81
1m
∴6 + ½h = 13.6h ZL hfull
∴hfull = 0.458 m. = 458 mm.
X Y
U tube
EIPINI Chapter 4: Level Measurement Page 4-9
Example 4-7
A mercury (density 13600 kg/m3) u tube manometer Patm
is used to measure the level of a liquid (density
Patm
1000 kg/m3) in a 5 meter high container, with the
zero level of the manometer, 1 meter below the 5m
H
bottom of the tank. Calculate the level in the tank if
the manometer reading is 300 mm (0.3 m).
Comparing pressure on the XY line: 1m
Patm+1000×(H+1+0.15)×9.81=Patm+13600×0.3×9.81 ZL 0.3 m
∴H + 1.15 = 13.6×0.3 X Y
∴H = 4.08 – 1.15
= 2.93 meter U tube
Example 4-8
A mercury (density 13600 kg/m3) well type
Patm
manometer with A2/A1 = 0.01, is used to measure
the level of a liquid (density 1000 kg/m3) in a
5 meter high container, with the zero level of the 5m
manometer, 1 meter below the bottom of the tank. P2=Patm
H
Calculate the level in the tank if the manometer
reading is 200 mm (0.2 m).
The input pressure, P1, on the well side of the 1m
P1 0.2 m
manometer, with respect to the zero line, is:
P1 = Patm + 1000×(H + 1)×9.81 ZL
and P2 = Patm
∴P1 – P2 = Patm + 1000×(H + 1)×9.81 – Patm Well type
d manometer
= 1000×(H+1)×9.81 = 9810(H+1)
Furthermore, from Equation 2-5 (Chapter 2):
P1 – P2 = ρhg(1 + A2/A1) = 13600×0.2×9.81×(1 + 0.01) = 2747×9.81 = 26948
Therefore 9810(H+1) = 26948 ⇒ H+1 = 2.747
∴H = 1.747 m.
Note: If we did wish to include d, then d = (A2/A1)h = 0.01×0.2 = 0.002 m
P1=Patm+1000×(H+1+d)×9.81=Patm+1000×(H+1+0.002)×9.81=9810(H+1.002)
And P1 – P2 = 9810(H+1.002) = 26948 ⇒ H + 1.002 = 2.747
∴H = 1.745 m
EIPINI Chapter 4: Level Measurement Page 4-10
L L
H H
z LP HP
LP HP
z
DP transmitter
ZL h
U tube
manometer
Figure 4-8 (a) Figure 4-8 (b)
EIPINI Chapter 4: Level Measurement Page 4-11
Example 4-9
A DP transmitter is used to measure the level of a liquid, with density 1000 kg/m3 in a
closed tank, that can store liquid to a maximum level of 5 meter. The DP transmitter is
installed 1 meter below the bottom of the tank. Calculate the input differential pressure
to the DP transmitter when the tank is empty and when it is full.
Pt Pt
Pt Pt
5m 5m
1m 1m
P2 P1 P2 P1
5m 5m
1m 1m
hempty ZL hfull ZL
u-tube u-tube
X Y X Y
Empty: Full:
Example 4-11
The level of a liquid (density 1000
Pt Pt
kg/m3) in a closed tank, 5 meter
high, is measured with a mercury
(density 13600 kg/m3) well type 5m
H Please note that
manometer (A2/A1=0.01). The the well (high
manometer is installed with its pressure input)
zero line 1 meter below the bottom is connected to
the outside tube
of the tank. If the reading on the 1m P2
manometer is 0.2 m, calculate the
0.2 m P1
level H of the fluid in the tank. d=.002
Zero line
P1 = Pt + 1000×(5+1)×9.81 = Pt + 58860
X Y
P2 = Pt + 1000×[H+(1-0.2)]×9.81 = Pt + 9810×(H+0.8)
∴P1 – P2 = 58860 – 9810(H+0.8) = 58860 – 9810H – 7848 = 51012 – 9810H…(1)
And P1 – P2 = ρhg(1 + A2/A1) = 13600×0.2×9.81×(1 + 0.01) = 26950 ………… (2)
From (1) and (2): 51012 – 9810H = 26950 ⇒ 9810H = 51012 – 26950 = 24062
∴H = 2.453 m.
(The exact value of H may be obtained if we equate pressures on the XY line:
Pt+1000×[H+(1–0.2)]×9.81+13600×(0.2+0.002)×9.81=Pt+1000×(5+1+0.002)×9.81
∴[H+0.8] + 13.6×0.202 = 6.002 ⇒ H + 0.8 + 2.747 = 6.002 ⇒ H + 3.547 = 6.002
∴H = 2.455 m.)
Example 4-12 Pt
The level of a liquid, with density ρℓ, stored in a Pt
closed tank, L meters high, is measured with a u
tube manometer, using a manometer liquid, with
L
density ρm (ρm > ρℓ). The manometer is installed H
with its zero line a distance z meter below the
bottom of the tank. When the fluid level in the ρℓ z
container is H meter, the reading obtained from
the manometer is h meter. Derive an equation h ZL
that can be used to calculate the liquid level H, in X Y
terms of the manometer reading h and in terms of
the constants ρℓ, ρm, L, and z. ρm
Equating pressures at points X and Y:
Pt + ρℓ×(H + z - ½h)×g + ρm×h×g = Pt + ρℓ×(L + z + ½h)×g
∴ρℓ×(H + z - ½h)×g + ρm×h×g = ρℓ×(L + z + ½h)×g
∴ρℓ×(H + z - ½h) + ρm×h = ρℓ×(L + z + ½h)
∴ρℓH + ρℓz - ½ρℓh + ρmh = ρℓL + ρℓz + ½ρℓh
∴ρℓH = ρℓL + ρℓz + ½ρℓh - ρℓz + ½ρℓh - ρmh = ρℓL + ρℓh - ρmh
ρ ⎛ ρ ⎞
⎜ ⎟
∴H = L + h - m h = L + ⎜1 − m ⎟h (interestingly, z doesn’t play a role)
ρ ⎜ ρ ⎟
l ⎝ l ⎠
EIPINI Chapter 4: Level Measurement Page 4-13
Air/gas
Pressure Filter supply
regulator
Bubbler
sight glass
Dip tube
Air/gas
Pressure Filter supply
regulator
Bubbler
sight glass
Dip tube
5. TEMPERATURE MEASUREMENT
The purpose of this chapter is to introduce students to fundamental concepts related to
temperature and to discuss important industrial methods used to measure temperature.
5.1 INTRODUCTION
Galileo Galilei (picture to the left), was born in the year
1564 in Pisa, Italy. This brilliant scientist and astronomer is
most famous for his early development of the telescope,
being the first to see the moons of Jupiter and other celestial
objects. Through his work in astronomy, Galileo supported
the theory originated by Copernicus, that the earth moved
around the Sun. Eventually, and under immense pressure
from the Church, he publicly retracted his support for the
Copernicus theory. Even so, Galileo was placed under house
arrest and spent the last 10 years of his life in almost
complete seclusion, having dared to offend the Church.
Galileo was indeed an incredible man. He was the first to realise how the swinging
of a pendulum, could be used to measure time and he is also credited to have
developed the first device to indicate changes in temperature.
This instrument did not measure temperature as such, and is
therefore called a thermoscope. (Vincenzo Viviani states in
his Vita di Galileo, that the thermoscope was conceived by
Galileo in 1597.) Galileo’s thermoscope was based on the
fascinating idea that air will expand when heated and contract
when cooled. The instrument consisted of a glass bottle about
the size of an egg, with a long glass neck. Air in the bottle
was heated with the hands and the open end side of the tube
immersed partially in a vessel containing water or wine. As
the air in the tube cools, it contracts and the liquid is drawn Cold Hot
into the tube. Once an equilibrium point is established, a rise
in temperature increases the volume of the air, forcing fluid back down; a fall in
temperature reduces the volume of the air, drawing more fluid into the tube.
Example 5-1
Convert the absolute zero of temperature (–273.15 °C), to degrees Fahrenheit.
F = (9/5)C + 32 = (9/5)×(-273.15) + 32 = -459.67 °F ≈ -460 °F
Example 5-2
Convert 86 °F into degrees Celsius, Rankine and Kelvin.
C = (5/9)(F – 32) = (5/9)×(86 – 32) = 30 °C
R = F + 460 = 86 + 460 = 546 R
K = C + 273 = 30 + 273 = 303 K
Example 5-3
Determine the temperature in degrees Fahrenheit that would translate into the same
value in degrees Celsius.
In Equation 5-2 (or Equation 5-1), set C = F:
∴F = (5/9)(F – 32) ⇒ 9F = 5F – 160 ⇒ 4F = -160 ⇒ F = -40 °F
Example 5-4
Convert 0 °F to degrees Celsius.
C = (5/9)(F – 32) = (5/9)(0 – 32) = -17.78 °C.
• The oxygen point: The boiling point of liquid oxygen: –182.97 °C.
• The ice point: The melting point of pure ice: 0 °C.
• The boiling point: The boiling point of pure water: 100 °C.
• The sulphur point: The boiling point of pure sulphur: 444.6 °C.
• The silver point: The melting point of silver: 961.78 °C.
• The gold point: The melting point of gold: 1064.18 °C.
The background of the glass tube is covered with white enamel and the front of
the glass tube forms a magnifying glass that enlarges the liquid column and
facilitates with reading the temperature.
In Figure 5-1 (a), an all glass thermometer is depicted, with its scale etched
into the stem. Liquid in glass thermometers are fragile and for industrial use, the
thermometer is mounted in a protective housing and the scale is engraved on a
separate plate that is part of the protective case. An industrial thermometer is
shown in Figure 5-1 (b).
Protective case
Scale Scale
(etched) (on plate)
Bore Bore
Socket
Bulb
Figure 5-1 (a) Figure 5-1 (b) Heat conducting
medium
Liquids used in glass thermometers
Liquid Temperature range (Celsius)
Mercury -35 to +510
Alcohol -80 to +70
Pentane -200 to +30
Toluene -80 to +100
Creosote -5 to +200
EIPINI Chapter 5: Temperature Measurement Page 5-6
Shaft
Figure 5-5 (a) Figure 5-5 (b) Guide
A bi-metal thermometer uses a bi-metal strip, shaped in
Stem
a helix or spiral form, as shown in Figure 5-6. The one
end is fixed and the other end is free to rotate as the helix
curls in or out with changing temperature. A shaft and Helical
pointer is linked to the rotating helix, to indicate the bi-metal
temperature. The stem is filled with silicone fluid, to element
provide damping and thermal conductivity between the Figure 5-6
stem and bi-metal strip.
EIPINI Chapter 5: Temperature Measurement Page 5-8
R (ohm)
Partial temperature/
800 resistance table for
Nickel
platinum (Ro=100Ω)
600
Temperature Resistance
400 Copper (°C) (Ω )
-200 17.26
200 Platinum -100 59.64
t (°C) -80 67.83
-60 75.96
-200 0 200 400 600 800 -40 84.03
-20 92.04
Figure 5-8 0 100.00
20 107.92
40 115.78
An excerpt from the resistance table for a Minco 60 123.60
platinum thermometer (code PB), is given in Table 5-1. 80 131.38
100 139.11
Example 5-6 200 177.04
In the temperature range 0 °C ≤ t < 850 °C, the 300 213.81
resistance / temperature relation of a platinum 400 249.41
500 283.84
thermometer with ice point resistance of 100 Ω, is
600 317.09
described by the second order expression: 700 349.18
Rt = 100×[1 + 0.0039692×t – (5.8495×10-7)×t2]. Table 5-1
a) Calculate the resistance of the thermometer at 50 °C.
b) Calculate the temperature when a resistance of 110 Ω is measured.
a) R50 = 100×(1 + 0.0039692×50 – 5.8495×10-7×502)
= 100×(1 + 0.19846 - 0.001462375) = 119.7 Ω.
b) Rt = 100[1 + (3.9692×10-3)×t – (5.8495×10-7)t2].
∴110 = 100[1 + (3.9692×10-3)×t – (5.8495×10-7)t2]
∴110 = 100 + 0.39692t – 5.8495×10-5t2
∴(5.8495×10-5)t2 + (– 0.39692)t + 10 = 0 [ax2+bx+c=0]
−5
0.39692 ± ( 0.39692 ) 2
− 4 × ( 5.8495 × 10 ) × 10 b 2 − 4ac
∴t =
−5
[x= − b ± ]
2 × ( 5.8495 × 10 ) 2a
0.39692 ± 0.39396
= = 25.301 °C (ignoring 6760 °C.)
−6
116.99 × 10
EIPINI Chapter 5: Temperature Measurement Page 5-10
Example 5-7
The Wheatstone bridge in Figure 5-10, is supplied from a 10 V battery, and used to
measure temperature with a platinum resistance RT. The value of each of the fixed
resistors R1, R2 and R3, is 100 Ω (in a practical circuit, R1 and R2 would be chosen
much higher to improve bridge linearity and minimize errors caused by I2R losses in
the thermometer). The meter M measures the thermometer output voltage V, which
is given by: V = VAB = VAC - VBC. At 0 °C, the resistance RT, of the platinum RTD,
is 100 Ω, while at 100 °C, the resistance RT increases to 139.1 Ω. Calculate the
output voltage V, of the thermometer, when the temperature is 0 °C and when the
temperature is 100 °C. Assume that the measuring device draws negligible current
from the circuit.
0 °C: I2 I1
R2 R1
I1 = 10/(100 + 100) = 50.00 mA.
∴VAC = RT×I1 = 100×(50×10-3) 10V 100Ω V 100Ω
= 5 volt B A
I2 = 10/(100 + 100) = 50.00 mA. R3
VBC VAB RT VAC
∴VBC = R3×I2 = 100×(50×10 ) -3
100Ω 100Ω
= 5 volt
∴V = VAB = VAC – VBC (using Kirchoff’s C 0 °C C
= 5 – 5 = 0 volt voltage law: VAB-VAC+VBC=0)
We of course expected this result, as the bridge is indeed balanced at 0 °C.
100 °C:
I1=10/(100+139.1)=41.82 mA. I2 I1
R2 R1
∴VAC = RT×I1 = 139.1×(41.82×10-3)
= 5.817 volt 10V 100Ω V 100Ω
I2 = 10/(100 + 100) = 50.00 mA. B A
∴VBC = R3×I2 = 100×(50×10-3) VBC
R3 RT
VAC
= 5 volt 100Ω 139.1Ω
∴V = VAB = VAC – VBC
= 5.817 – 5.000 = 0.817 volt C 100 °C C
Example 5-8
The Wheatstone bridge in Figure 5-10, is supplied from a 10 V battery, and used to
measure temperature with a platinum resistance RT. The value of each of the fixed
resistors R1, R2 and R3, is 100 Ω. The meter M measures the thermometer output
voltage V, which is given by: V = VAB = VAC - VBC.
a) Derive an expression for the resistance RT of the platinum element, in terms of
the thermometer output voltage V.
b) Derive an expression for the temperature t in terms of RT, if in the temperature
range 0 °C ≤ t < 850 °C, the relationship between the resistance RT of the
thermometer and the temperature t, is given by:
RT = 100[1 + (3.9692×10-3)×t – (5.8495×10-7)t2].
c) Calculate the measured temperature when the bridge output voltage is 0, 0.2, 0.4,
0.6, 0.8, and 1.0 volt.
EIPINI Chapter 5: Temperature Measurement Page 5-12
10
a) I1 = ampere
RT + 100 I2 I1
10 × RT R 2 R1
∴VAC = I1×RT = volt
RT + 100 10V 100Ω V 100Ω
10 B A
and I2 = = 50 milliamp.
100 + 100 R3
∴VBC = I2×R3 = (50×10-3)×100 VBC
100Ω
VAC
RT
= 5 volt
∴V = VAB = VAC - VBC C C
10 × RT 10 × RT
∴V = -5⇒ = V + 5 ⇒ 10×RT = (RT + 100)×(V + 5)
RT + 100 RT + 100
∴10×RT = V×RT + 100×V + 5×RT + 100×5
∴10×RT – 5×RT - V×RT = 100×V + 100×5
∴(10 – 5 – V)×RT = 100×(V + 5) ⇒ (5 – V)×RT = 100×(V + 5)
100(V + 5)
∴RT = …………………………….………………. Equation (a)
5-V
b) RT = 100[1 + (3.9692×10-3)×t – (5.8495×10-7)t2]
∴RT = 100 + 0.39692×t – (5.8495×10-5)×t2
∴(5.8495×10-5)×t2 - 0.39692×t + RT – 100 = 0
5.8495 ×10 - 5 2 0.39692 1
∴ ×t - ×t + ×(RT-100) = 0
5.8495 ×10 - 5 5.8495 ×10- 5 5.8495 ×10- 5
∴ t2 - 6786t + [17095×(RT-100)] = 0 [x2+bx+c=0]
6786 − (6786) 2 − 4 × 17095 × (RT − 100) b 2 − 4c
∴t = [x= − b ± ]
2 2
6786 − 46.05 × 106 − 68380(RT − 100)
∴t = …….………….. Equation (b)
2
(Note: From Example 5-6 we already know that we can ignore the second
solution associated with the positive square root.)
c) V = 0: From Equation (a): RT = 100(0 + 5)/(5 - 0) = 100 Ω
From Equation (b): t = {6786 - √[46.05×106 - 68380(100 - 100)]}/2
= [6786 - 6786]/2 = 0 °C
V = 0.2: From Equation (a): RT = 100(0.2 + 5)/(5 - 0.2) = 108.33 Ω
From Equation (b): t = {6786 - √[46.05×106-68380(108.33-100)]}/2
= [6786 - √45.48×106]/2 = 21.05 °C
V = 0.4: RT = 100(0.4 + 5)/(5 - 0.4) = 117.39 Ω
t = {6786 - √[46.05×106 - 68380(117.39 - 100)]}/2 = 44.09 °C
V = 0.6: RT = 100(0.6 + 5)/(5 – 0.6) = 127.27 Ω
t = {6786 - √[46.05×106 - 68380(127.27 - 100)]}/2 = 69.40 °C
V = 0.8: RT = 100(0.8 + 5)/(5 - 0.8) = 138.1 Ω
t = {6786 - √[46.05×106 - 68380(138.1 - 100)]}/2 = 97.37 °C
EIPINI Chapter 5: Temperature Measurement Page 5-13
R1
R2 c Rlead
V
E
A B b Rlead
a Rlead RT
R3
C Figure 5-11
The four wire method, depicted in Figure 5-12, alleviates many problems
associated with the Wheatstone bridge. Instead of using a Wheatstone bridge
configuration, a current source is employed to supply a constant current I, to the
thermometer, through wires a and d. A high impedance voltmeter M, measures the
voltage developed across the thermometer, via wires b and c. The measured
voltage is directly proportional to the resistance of the thermometer, so only the
conversion from resistance to temperature is necessary. Wires b and c only act as
voltage sense leads, and with M a high impedance meter, virtually no current flows
in wires b and c, and therefore no voltage drop in these leads and thus no lead
resistance error in the measurement. The disadvantage of the four wire method is
that we need one more wire than with the three wire method, but it is a small price
to pay if we are at all concerned with the accuracy of the temperature
measurement. The four wire method is however gaining in popularity in industry.
a Rlead
b Rlead
I M
RT
c Rlead Figure 5-12
d Rlead
Example 5-9
A Wheatstone bridge, with fixed resistors R1, R2 and R3, of 100 Ω and supplied from a
10 V battery, is used to measure a temperature of 100 °C with a platinum resistance
RT, that has a resistance of 139.1 Ω at 100 °C . The thermometer is connected to the
bridge with leadwires that have a resistance Rlead of 10 Ω at 20 °C and the thermometer
is correctly calibrated for an ambient temperature of 20 °C. Assume that the lead
resistance changes to 11 Ω when the ambient temperature rises to 30 °C.
EIPINI Chapter 5: Temperature Measurement Page 5-15
10V V Process
10V V Process
10Ω 11Ω
temberature temberature
100°C 100°C
b 139.1Ω b 139.1Ω
100Ω 100Ω
a RT a RT
10Ω 11Ω
I1 = 10/(100+10+139.1+10) = 0.0386 A I1 = 10/(100+11+139.1+11) = 0.0383 A
I2 = 10/(100+100) = 0.05 A I2 = 10/(100+100) = 0.05 A
V = (10+139.1+10)×0.0386 - 100×0.05 V = (11+139.1+11)×0.0383 - 100×0.05
= 6.141 – 5 = 1.141 V (calibrated 100°C) = 6.17 – 5 = 1.17 V (105°C)
Therefore, if the bridge was correctly calibrated for an ambient temperature of 20 °C,
an error of 0.029 V will occur when the ambient temperature rises to 30 °C. So instead
of reading 100 °C, it would measure approximately 105 °C, an error of 5 °C or 5%.
b) Measuring with three leads
Calibrated bridge voltage at Bridge output voltage at
20 °C ambient temperature 30 °C ambient temperature
I2 I1 100Ω I2 I1 100Ω
100Ω Process
100Ω Process
10Ω temberature
100°C 11Ω temberature
100°C
c c
10V V 10Ω 10V V 11Ω
b b
139.1Ω 139.1Ω
100Ω a RT 100Ω a RT
10Ω 11Ω
I1 = 10/(100+10+139.1+10) = 0.0386 A I1 = 10/(100+11+139.1+11) = 0.0383 A
I2 = 10/(100+100) = 0.05 A I2 = 10/(100+100) = 0.05 A
V = (139.1 + 10)×0.0386 - 100×0.05 V = (139.1 + 11)×0.0383 - 100×0.05
= 5.755 – 5 = 0.755 V (calibrated 100°C) = 5.749 – 5 = 0.749 V (99.13°C)
Therefore, measuring with the three wire method, an error of only 0.006 V will
occur. If the thermometer was correctly calibrated to read 100 °C at an ambient
temperature of 20 °C, then if the ambient temperature rises to 30 °C, the reading of
0.749 V would correspond to a value of 138.76 Ω for RT with no rise in ambient
temperature, or in terms of temperature, 99.13 °C, an error of 0.87 °C or 0.87%. This
is clearly a remarkable improvement over the two wire method.
EIPINI Chapter 5: Temperature Measurement Page 5-16
5.10 THERMOCOUPLES
5.10.1 The Seebeck effect
If two dissimilar metals are joined together to form a closed loop, and if one
junction is kept at a different temperature from the other, an electromotive force is
generated and electric current will flow in the closed loop.
Metal A
Metal B Metal B
Seebeck emf
This very important discovery in the field of thermometry, was made by T.J. Seebeck
in the year 1821. Experiments by Seebeck and others, have shown that the generated
emf (called the Seebeck voltage in his honor), is relative, in a predictable manner, to
the difference in temperature between the two junctions. So, if the temperature of one
junction is kept at a known value, the temperature of the other junction can be
determined by the amount of voltage produced. This discovery resulted in the
temperature sensor that we know today as the thermocouple.
The Seebeck voltage is made up of two components: the Peltier voltage generated at
the junctions, plus the Thomson voltage generated in the wires by the temperature
gradient. The Peltier voltage is proportional to the temperature of each junction while
the Thomson voltage is proportional to the square of the temperature difference
between the two junctions. It is the Thomson voltage that accounts for most of the
observed voltage and non-linearity in thermocouple response.
T3
T1 T2
T4 Figure 5-13
In Figure 5-13, a thermocouple is shown with junction temperatures at T1 and T2.
Along the thermocouple wires, the temperature is T3 and T4. The thermocouple emf is,
however, still a function of only the temperature gradient T2 – T1.
EIPINI Chapter 5: Temperature Measurement Page 5-17
Metal A
T1 T2
T1 T2 T2 T3 T1 T3
Cold junction
Ice bath
Figure 5-19
EIPINI Chapter 5: Temperature Measurement Page 5-20
Temperature (°C) versus emf (μV) for type K thermocouple with 0 °C reference.
deg C 0 1 2 3 4 5 6 7 8 9
0 0 39 79 119 158 198 238 277 317 357
10 397 437 477 517 557 597 637 677 718 758
20 798 838 879 919 960 1000 1041 1081 1122 1163
30 1203 1244 1285 1326 1366 1407 1448 1489 1530 1571
40 1612 1653 1694 1735 1776 1817 1858 1899 1941 1982
50 2023 2064 2106 2147 2188 2230 2271 2312 2354 2395
60 2436 2478 2519 2561 2602 2644 2685 2727 2768 2810
70 2851 2893 2934 2976 3017 3059 3100 3142 3184 3225
80 3267 3308 3350 3391 3433 3474 3516 3557 3599 3640
90 3682 3723 3765 3806 3848 3889 3931 3972 4013 4055
100 4096 4138 4179 4220 4262 4303 4344 4385 4427 4468
Table 5-3
We could store these look-up table values in a computer and use the table to convert
between emf and temperature. A more viable approach used by manufacturers
however, is to approximate the table values using a power series polynomial (see
Example 5-14) and allow the instrument’s microprocessor or process computer, to
calculate temperature from emf or the emf from temperature (inverse polynomial).
The ice bath cold junction is not considered practical anymore. Instead the terminals
connecting the thermocouple to the measuring device, are now assumed to play the
role of the reference junction or 'cold junction', as it is still called today. The reference
junction temperature may now be kept, for example, at room temperature, where the
junction temperature is measured with an auxiliary temperature sensor, such as a RTD.
According to the law of intermediate temperatures, the thermocouple voltage that
corresponds to the cold junction temperature, may be added to the measured
thermocouple voltage. The true temperature of the hot junction, with respect to 0 °C,
can then be determined from this augmented voltage.
Example 5-10
Calculate the average sensitivity (μV/°C) of a type K thermocouple in the temperature
range 0 °C to 100 °C.
From Table 5-3 the change in emf developed by a type K thermocouple from 0 °C to
100 °C, is 4096 μV. The average sensitivity is therefore 4096/100 = 40.96 μV/°C.
EIPINI Chapter 5: Temperature Measurement Page 5-21
Example 5-11
The cold junction of a type K thermocouple is kept at 0 °C. Use Table 5-3 to
determine the temperature if the measured voltage is a) 798 μV and b) 2602 μV.
a) 20 °C b) 64 °C
Example 5-12
The relationship between emf and temperature for a certain (imaginary) thermocouple,
is described by the relation: v = t2, where v is the generated thermocouple emf in
microvolt (μV), and t the temperature difference in °C, between the hot junction and
0 °C. If the thermocouple emf reading is 3000 microvolt and the temperature of the
cold junction is 25 °C, calculate the temperature of the hot junction.
Emf corresponding to (25 – 0) °C = 252 = 625 μV 0°C 25°C 25°C T
Total emf (T - 0) = 3000 + 625 = 3625 μV 625μV 3000μV
According to the law of intermediate temperatures:
Hot junction temperature T = √v = √3625 = 60.21 °C 0°C T
(T is NOT = √3000 + 25 = 54.77 + 25 = 79.77 °C) 3625μV
Example 5-13
An unknown temperature is measured with a type K thermocouple. A thermocouple
voltage of 2602 μV is measured. If the cold junction temperature is 20 °C, calculate
the process temperature, measured by the hot junction side of the thermocouple.
From Table 5-3, the cold junction (20°C) emf is
0°C 20°C 20°C T
798 μV. According to the law of intermediate
798μV 2602μV
temperatures, the correction voltage of 798 μV
should be added to the measured voltage of 0°C T
2602 μV, to obtain 3400 μV. The corrected voltage 3400μV
represents the thermocouple emf that would be
obtained, if the reference junction was kept at 0 °C.
Again from Table 5-3, the temperature that 84°C
corresponds to 3400 μV, is somewhere between
83 °C and 84 °C. To find the correct temperature,
1°C
we must use linear extrapolation between these
?
two values. The difference between 3433 μV
(84 °C) and 3391 μV (83 °C) is 42 μV, while X
3400 μV is 9 μV more than 3391 μV (83 °C). 83°C
9μV
Therefore the temperature we are looking for is
42μV
83 °C plus (9/42) °C, which is 83.2143 °C.
The way NOT to calculate the hot junction
3391μV 3400μV 3433μV
temperature, is to look up the measured voltage
(2602 μV) as 64 °C and then to add the cold X 1 9
junction temperature of 20 °C, to obtain 84 °C. = ∴X = ×1
9 42 42
This is NOT CORRECT.
EIPINI Chapter 5: Temperature Measurement Page 5-22
Example 5-14
The following equations are provided for a type K thermocouple, to calculate
temperature (t in °C) from emf (v in μV) in the temperature range 0 °C to 500 °C:
t = (2.508355×10-2)v + (7.860106×10-8)v2 - (2.503131×10-10)v3
+ (8.315270×10-14)v4 – (1.228034×10-17)v5 + (9.804036×10-22)v6
- (4.413030×10-26)v7 + (1.057734×10-30)v8 – (1.052755×10-35)v9 ,
Equation (a)
and to calculate emf (v in μV) from temperature (t in °C) in the temperature range
0 °C to 1372 °C:
v = -17.600413686 + 38.921204975t + (1.8558770032×10-2)t2
- (9.9457592874×10-5)t3 + (3.1840945719×10-7)t4
- (5.6072844889×10-10)t5 + (5.6075059059×10-13)t6
- (3.2020720003×10-16)t7 + (9.7151147152×10-20)t8
−4
– (1.2104721275×10-23)t9 + 118.5976 e − 1.183432 ×10 (t − 126.9686) 2
Equation (b)
a) Use Equation (a) to calculate the temperature when the emf is 2602 μV.
b) Use Equation (b) to calculate the emf when the temperature is 20 °C.
a) t = (2.508355×10-2)×2602 + (7.860106×10-8)×26022 - (2.503131×10-10)×26023
+ (8.315270×10-14)×26024 – (1.228034×10-17)×26025 + (9.804036×10-22)×26026
- (4.413030×10-26)×26027 + (1.057734×10-30)×26028 – (1.052755×10-35)×26029
= 65.26740 + 0.5321609 – 4.409664 + 3.811584 – 1.464694 + 0.3042627
- 0.03563592 + 0.002222464 – 0.05755631
= 63.95 °C
b) v = -17.600413686 + 38.921204975×20 + (1.8558770032×10-2)×202
- (9.9457592874×10-5)×203 + (3.1840945719×10-7)×204
- (5.6072844889×10-10)×205 + (5.6075059059×10-13)×206
- (3.2020720003×10-16)×207 + (9.7151147152×10-20)×208
−4 2
– (1.2104721275×10-23)×209 + 118.5976 e − 1.183432 ×10 (20 − 126.9686)
= -17.600414 + 778.4241 + 7.423508 – 0.7956607 + 0.05094551
- 0.001794331 + 35.88804×10-6 – 409.8652×10-9 + 2.487069×10-9
- 6.197617×10-12 + 30.61899
= 798.1 microvolt
(Compare these answers to those of Example 5-11)
5.10.6 Compensating leads
Thermocouple thermometers are normally installed some distance away from the
voltmeter or computer that measures the emf generated by the thermocouple. For
this purpose, cheaper and lower grade thermocouple wires, called extension wire or
compensating leads, are used to connect the thermocouple to the measuring device
at the reference junction. Compensating leads have the same thermoelectric
EIPINI Chapter 5: Temperature Measurement Page 5-23
properties as the thermocouple and do not introduce a significant error into the
temperature measurement. Compensating leads must be matched to the
thermocouple and for each type of thermocouple, corresponding extension leads
are available.
5.11 THERMISTORS
A thermistor is similar to a resistance thermometer, but a semiconductor material is
used instead of a metal. A distinct advantage of thermistors over resistance
thermometers, is that their temperature coefficient of resistance is approximately
ten times higher than that of a resistance thermometer, causing thermistors to be
much more sensitive temperature detectors. The resistance change with
temperature is however very nonlinear, and unlike RTD’s, most thermistors have a
negative temperature coefficient of resistance, that is, the resistance of a thermistor
decreases with increasing temperature.
In Figure 5-20, the resistance curve of a typical thermistor is shown. The following
equation, called the Steinhart and Hart equation, is used to describe the resistance /
temperature relation for thermistors:
1
= A + Bln(R) + ClnR 3 ,
T
where T is absolute temperature in Kelvin, R is the thermistor resistance at
temperature T and A, B and C are coefficients that describe a given thermistor.
Resistance R
Figure 5-20
Temperature T
The resistivity of thermistors is also much larger than that of RTD’s, and
thermistors can therefore be made very small which means they will respond
quickly to temperature changes.
Thermistors cannot be used to measure high temperatures, compared to RTDs. In
fact, the maximum temperature of operation is sometimes only 100 or 200 oC. The
high resistivity as well as high temperature coefficient of resistance of the
thermistor, makes it feasible to use the simpler two wire technique when measuring
thermistor resistance.
EIPINI Chapter 5: Temperature Measurement Page 5-24
For example, a common thermistor value is 5000 ohms at 25°C. With a typical
temperature coefficient of resistance of 0.04/°C, a measurement lead resistance of
10 Ω, produces only 0.05 °C error. This error is a factor of 500 times less than the
equivalent RTD error. Thermistors have not gained nearly the popularity of RTDs or
even thermocouples in industry due to their limited span as well as other
disadvantages. Since thermistors are semiconductor devices, they are quite susceptible
to permanent decalibration when exposed to high temperatures. In addition,
thermistors are quite fragile and great care must be taken to mount them so that they
are not exposed to shock or vibration.
Thermistors
A thermistor is a resistive device composed of metal oxides formed into a bead and
encapsulated in epoxy or glass. A typical thermistor shows a large negative
temperature coefficient. Resistance drops dramatically and non-linearly with
temperature. Sensitivity is many times that of RTDs but useful temperature range is
limited. There are wide variations of performance and price between thermistors from
different sources. Typical benefits are:
• Low sensor cost: Basic thermistors are quite inexpensive. However, models with
tighter interchangeability or extended temperature ranges often cost more than
RTDs.
• High sensitivity: A thermistor may change resistance by tens of ohms per degree
temperature change, versus a fraction of an ohm for RTDs.
• Point sensing: A thermistor bead can be made the size of a pin head for small area
sensing.
6. PROCESS CONTROL
The purpose of this chapter is to introduce students to industrial process control
methods used in industry to ensure that manufactured products meet predetermined
quality requirements. This is accomplished by continuously monitoring the production
process and automatically correcting or minimizing any deviations from the required
specifications, that may be detected during the manufacturing process.
6.1 INTRODUCTION
The history of automatic control goes back many centuries. Water level control may
already be identified in water clocks used in the middle ages while mechanical clocks,
making their appearance in the 1200’s, used an escapement mechanism that may be
described in terms of feedback control. With the advent of the steam engine, it was
evident that some sort of speed control was needed. In what is generally considered a
major event in the history of feedback control, James Watt, completed
the design of the centrifugal flyball governor for regulating the
speed of the rotary steam engine, in 1788. This device
employed two pivoted rotating flyballs which were flung
outward by centrifugal force. As the speed of rotation increased,
the flyweights swung further out, operating a steam flow throttling
valve which slowed the engine down. Thus, a constant speed was
achieved automatically. From this point onwards, many new developments followed
to improve control systems as well as mathematical research to understand physical
processes. A major invention, however, was the introduction of proportional, integral
and derivative (PID) control, which was formulated in 1922 by Nicholas Minorsky
(1885-1970). Observing the way in which a helmsman steered a ship and compensate
for the disturbances from the ocean, motivated this threefold control strategy.
Although advanced tools, such as model predictive controllers, are available today,
the Proportional, Integral, Derivative (PID) control strategy, is still the most widely
used in modern industry, controlling more than 95% of closed loop industrial
processes.
If the water level reaches the desired value (also called the set point), the operator
must try to keep the inflow equal to the outflow. This manual control example
illustrates essentially how we would expect an automatic control system to behave.
In Figure 6-2, an
arrangement for Inflow
automatic control
Inflow Level Level Outflow
of the water level
valve controller detector valve
in the container,
is shown. A
level measuring
instrument senses Figure 6-2
the water level in Water inlet Water level Outflow
the container, and control
feeds this value back to the controller. The basic function of the controller is to
calculate the difference between the desired water level and the measured water level,
and to determine a sensible control action, that will minimize this difference. Once the
controller has calculated an appropriate control action, this value is transmitted as a
pneumatic signal to the inflow valve, which is adjusted accordingly. A well-designed
controller should respond quickly to changes in the water level and always steer the
level in the water container back to the required level, in a disciplined way.
Example 6-1
A control system for a heat exchanger is shown in Figure 6-4. Cold water enters the
container at point 1 and is heated by steam, entering the steam line 3, at point 2. The
heated water leaves the system at point 7. The temperature of the hot water is
measured by a thermometer at point 6 and presented to the temperature controller,
block 5. The controller operates a steam valve 4, such that when the hot water is below
the required temperature, the steam flow rate through the steam line is increased,
allowing more heat transfer to the cold water, and when the hot water is above the
required temperature, the steam flow is reduced. Identify: a) the controlled variable, b)
the measured value, c) the manipulated variable and d) some disturbance variables.
Feedforward control systems, on the other hand, will measure disturbance variables
that directly influence the controlled variable. It is conceivable, for example, to
measure the outflow in the water level control system in Figure 6-2, and to base the
inflow control on the value of this disturbance variable. Because the source of the
level disturbance is monitored and acted upon, a change in water level is anticipated
and corrected before it even occurs.
Whereas feedback systems must wait for a deviation to occur before corrective
action is taken, the principle advantage of feedforward control systems is that they act
before deviations occur. It is clear that feedforward control alone would not be
adequate for the water level control system, as the water level would eventually drift
away from set point because of measuring and modelling errors. In general,
feedforward control is used to complement feedback control and to enhance system
performance.
Feedback control: Measure the controlled variable to determine the control strategy.
Feedforward control: Measure disturbance variables to determine the control
strategy.
Figure 6-6 shows the typical response to a step input of a process such as a control
valve receiving a command to open from 50 % to 60 %, resulting in the increase in
inflow from 50 % to 60 %.
Input
60 %
50 %
Time
Output
60 %
Dead time and first order lag
50 %
Time
td tf
Figure 6-6
The response curve in Figure 6-6 illustrates two kinds of time delays that may be
identified. Firstly a dead time td, when the system does not respond at all, and
secondly a first order delay tf, when the system output changes to its new final value
in an exponential fashion.
The dead time is associated with a delay in the transportation of material or
information. For example, it may take 20 minutes to send a signal to the Rover on
Mars. During this time, the vehicle can not respond to this control command, because
the control information is still being transported to the vehicle. The dead time is also
known as a distance-velocity lag, and the reason for this is illuminated in
Example 6-2. Another term used for the dead time, is transportation lag.
The first order lag is associated with a delay in the transfer of energy. For example,
when you step on the accelerator of a car, the car will take time to convert the
chemical energy in the fuel, into mechanical kinetic energy, and the speed increase
will be exponential in nature. The first order lag is also known as a transfer lag,
exponential lag or resistance-capacitance lag.
Dead time: Delay due to the time it takes information or material to be transported
from one point to another.
First order lag: Delay due to the time it takes energy to be transferred from one
point or form to another.
Example 6-2
In the system depicted Figure 6-7, a conveyor belt is used to transport material that
was deposited by a feeder, to the delivery station. A request for more material is given
to the feeder. Calculate the time it will take the sensor at the delivery site to detect a
reaction to the increased demand for material. The velocity of the conveyor belt is v
meter per second and the distance between the feeder and delivery point is d meter.
EIPINI Chapter 6: Process Control Page 6-7
Feeder
Sensor
d
Any increase in material delivered, must travel a distance d with a velocity v, before it
reaches the delivery point, where the sensor will detect the increase. The time it takes
the material to be transported to the sensor is t = distance/velocity = d/v seconds. The
dead time lag in this case is called a distance-velocity lag and it is a problem not easily
handled by a PID controller.
Therefore, if the controller output C is 0 %, the water inflow QIN will immediately
cease (QIN = 0 m3/sec.) and when C is 100 %, the water inflow will immediately rise
to its maximum value of QIN = 0.01 m3/sec.
The outflow demand or load QOUT will be represented by the percentage variable L.
L will be 0 % when the outflow QOUT = 0 m3/sec and L will equal 100 % when the
outflow is the maximum value of QOUT = 0.01 m3/sec.
The complete water level control system is shown in Figure 6-9.
Water inflow C
100%
Controller
Output C
Inflow C (0% – 100%) 100%
C Water level M
Inflow Level M Level
(0% – 100%)
valve controller detector
Figure 6-9
Water supply Water level Outflow L (0 – 100%)
control system
Before we can use our water tank, we must find out how it works. This means that
we need to understand how the inflow and outflow influence the water level in the
container. Clearly when C = L, the water level M will stay constant. When C > L, the
water level will rise and when C < L, the water level will drop. It is therefore easy to
ΔM
argue that the rate at which the water level changes, , is proportional to C – L. In
Δt
the example given at the end of this chapter in Appendix 6-1, the exact relationship
between M, C and L, was found to be:
dM
= 0.002×(C – L) percent per second.
dt
Finally, the error value will play a fundamental part in the operation of on-off and
PID control. The percentage error value E is defined as:
the inflow is switched on fully (100%) when the water level drops below a certain
level, say 40%, or it is closed completely (0%), when the water level rises above a
certain level, say 60%. With the aid of Equation 6-1, we could express the control law
for the water level on-off control as:
0% if E ≤ -10% (or M ≥ 60%)
C= unchanged if -10% < E < 10% (or 60% > M > 40%)
100% if E ≥ 10% (or M ≤ 40%)
A possible response of the water level to on-off control as well as the controller
action, is shown in Figure 6-10. On average, the water level remains 50%.
Water level (M)
60%
40%
Time
Water inflow (C)
100%
0% Time
Figure 6-10 Typical on/off control performance
On-off control will also work well for a room temperature control – the heater is
fully on, until the room temperature increases to a prescribed value, after which the
heater is switched off, and so on. Not all systems are suited to on-off control. Imagine
a speed cruise control in a car, set to 120 km/h, using on-off control. 118 km/h: full
acceleration, 122 km/h: no throttle. Surely not very comfortable driving.
A disadvantage of on-off control is the wear on the final control element, such as the
control valve in the present example that continually moves from one extreme position
to the other.
Example 6-3
A water container is 5 meter high and has a base area of 1 m2. The maximum water
inflow rate, as well as the maximum water outflow rate, is 0.01 m3/sec. On-off control
is used to regulate the water level in the tank, between 40% and 60% of its maximum
level (5 meter). Assume that the rate of water outflow, is regulated and constant at
60% of the maximum rate of outflow. Calculate:
a) The time period that the inflow valve is closed
b) The time period that the inflow valve is open.
Solution:
a) Water level when the inflow valve must close: 60% of 5 m = 0.6×5 = 3 meter.
Water level when the inflow valve must open: 40% of 5 m = 0.4×5 = 2 meter.
Constant demand from tank: 60% of 0.01 = 0.6×0.01 = 0.006 meter3/second.
EIPINI Chapter 6: Process Control Page 6-10
The inflow valve closes when the Max inflow (100%) = 0.01 m3/s
water level reaches the 3 m mark, and
stays closed until the water level Max level=60%=3 m
drops to the 2 m mark. Total volume 5m Min level=40%=2 m
of water that must leave the container
before the 2 m mark is reached:
Controlled and
V = Ah = 1×1 = 1 m3. 2
1m constant outflow
Rate of outflow Q, is 0.006 m3/s.
But volume = flow rate×time Outflow=60% of max (0.01) = 0.006 m3/s
∴V = Qt [(m3/sec)×(sec)]
∴tclosed = V/Q = 1/0.006 M
= 166.7 sec = 2.778 min 60%
b) The inflow valve opens when the
water level reaches the 2 m mark and 40%
stays open until the water level rises
to the 3 m mark. Total volume of t
water that must enter the container, is C
1 m3. The resultant rate at which
water flows into the tank is the rate of 100%
inflow, minus the outflow rate.
Q = 0.01–0.006 = 0.004 m3/sec. 0%
t
∴topen = V/Q = 1/0.004 2.778 4.167
min min
= 250 sec = 4.167 min.
On-off control: A control strategy in which the controller output switches the final
control element fully on or off to keep the controlled variable near set point.
A simple and logical proportional control law for the water level system can now be
formulated. The inflow equals 100% when the tank is empty, the inflow equals 50%
(bias) when the tank is half full (equal to the desired value) and the inflow equals 0%
when the tank is full (overflowing). This control strategy is depicted in Figure 6-11.
Water inflow C (%)
100 %
Figure 6-11
Basic
proportional
50 % control law
R (Bias)
Error (%)
-50 % 0% 50 % (E=50-M)
Tank full, Tank half full, Tank empty,
M = 100% M = 50% M = 0%
This is obviously a reverse acting control strategy, because the controller output
increases when the measured value decreases. A negative slope of the C-E line would
indicate a direct acting control strategy. It is now a simple matter to write down an
expression for this proportional control law, as the graph of C versus E, is a straight
line of the form: y = mx + c. From Figure 6-11: C = E + 50. The slope of this line is
clearly 1, because the controller output C, changes by 100% when the error E, changes
by 100%. The slope of the line is called the proportional gain KP, of the controller.
For the controller characteristic in Figure 6-11, the gain is therefore KP = 1.
If the controller gain, KP, is changed, the slope of the straight line, describing the
relationship between controller output C, and error value E, will also change. Graphs
of C versus E, for three different controller gains, KP = 2, KP = 1 and KP = 0.5, are
shown in Figure 6-12. (Note: Although a gain of 0.5, prescribes a theoretical error of
100% before the output reaches 100 %, such an error could of course not occur in our
water level control system.)
C (%) C (%) C (%)
100% 100% 100%
75%
R=50% R=50% R=50%
E(%) E(%) 25% E(%)
-25% 0% 25% -50% 0% 50% -100% -50% 0% 50% 100%
KP = 2 KP = 1 KP = 0.5
KP = 2 KP = 1 -100% KP = 0.5
pivoted beam, represents changes in the error E, while the right hand side, sweeps out
corresponding values of the controller output C. For KP = 2, an error range of 50%
(from -25% to 25%), sweeps out the complete output range of 100%. For KP = 1, an
error range of 100% results in an output change of 100%. For KP = 0.5, the error value
must change 200% to cover the whole output range. The error range that results in the
total change in output range, is called the proportional band, of the controller, and is
denoted by the percentage variable PB.
For a high gain (narrowband) system (KP =2 in Figure 6-13 for example), a small
error would cause a large reaction from the controller in his effort to correct the error.
A low gain (wideband) system (KP = 0.5 in Figure 6-13 for example), would act more
gently, because large errors, will cause a milder controller reaction. The relationship
between proportional band PB (in percent) and proportional gain KP, is given by:
100
PB = percent. …………...……………………….. Equation 6-2
K
P
We can now obtain a general expression for the output C, of a proportional
controller, in terms of the proportional band PB, the set point S, the measured value M
and the bias R. In general, for any gain KP and bias R, the relationship between the
controller output C, and the process error E, is a straight line of the form, y = mx + c,
with the slope of the line equal to the controller gain KP and the y intersect, equal to
the controller bias R.
C (%)
∴C = KPE + R ………………..………...….… (1)
From Equation 6-1: E = S – M …...............….. (2)
R (%) Slope=KP
100
And from Equation 6-2: KP = ………..…. (3)
PB
(2) and (3) in (1): E (%)
100
C = KPE + R = (S – M) + R ……………………..…..…. Equation 6-3
PB
100
In a typical application, S = 50 and R = 50, therefore C = (50 – M) + 50.
PB
EIPINI Chapter 6: Process Control Page 6-13
Example 6-4
A process is controlled by a proportional controller. The controller is programmed for
a positive gain (reverse acting controller) and proportional band of 80 %, a set point of
50 % and a bias of 50 %. The system error is calculated from the difference between
set point and the measured value i.e. E = 50 – M.
a) Calculate the proportional gain of the controller.
b) If the measured value of the process is indicated as 36 %, calculate the output of
the controller at that instant.
c) If the measured value of the process is indicated as 65 %, calculate the output of
the controller at that instant.
d) For the given settings, give the proportional control law (C as a function of E).
e) Draw a graph of the controller control law (C versus E).
a) From Equation 6-2: d) C = KPE + Bias
KP = 100/PB = 100/80 = 1.25. ∴ C = 1.25E + 50
b) From Equation 6-3: e)
C
C=
100
×(S – M) + R 100% ΔC
Slope =
PB ΔE
100
= ×(50 – 36) + 50 = 67.5%. =
100
80 80
c) From Equation 6-3: 50%
= 1.25
100
C= ×(S – M) + R
PB E
100
= ×(50 – 65) + 50 = 31.25%. -40% 0% 40%
80
Example 6-5 (students must study this important example very carefully)
The water level in a container is controlled by a proportional controller with KP = 1
and bias R = 50%. At a given moment, the outflow demand L, is 50% and the inflow
C, is 50%. The water level in the tank is also stable at its set point of 50% (the error is
zero, therefore the inflow is 50% which is equal to the outflow). The outflow demand
suddenly increases to L = 60%. Determine the new stable water level in the tank.
Solution: The moment C=50% C=60%
the outflow increases, the
water level will start to
drop, and it will continue S=50% S=50% Offset
M=50% (E=10%)
to drop until the inflow C,
M=40%
also increases to 60%.
The inflow is driven by
the error value according L=50% L=60%
to C = KPE + R = E + 50.
Therefore the controller output will reach 60% when E = 10%. This means, from
E = 50–M, that the water level M, must drop to 40%. Although trying to keep the
water level close to 50%, this example clearly illustrates that a proportional controller
is not able to force the measured value back to set point, when the demand is
EIPINI Chapter 6: Process Control Page 6-14
different from the bias value. Decreasing the proportional band (increasing KP), will
however keep the level closer to set point. For example, if KP = 2, the controller will
need an error of 5 % to increase its output to 60 %. The water level will thus drop to
only 45%.
The difference between the set point and the measured value that may occur in
proportional control, is called the offset, and it is a serious disadvantage of
proportional only control. It is also important to note that we could have restored the
level to 50 %, by changing the bias value of the controller to 60 %, if the outflow
remained constant at 60 %.
Example 6-6
The water level in a container is controlled by a proportional controller and control
valve. The controller gain is 1 and the bias is 50 %. Assume that the water inflow
delivered by the valve is equal (no time delay) to the control signal C, and that the
level detector instantly reflects the level M, of the water in the container. The
relationship between the water level M, and the inflow C and outflow L, is given by:
dM
= 0.002×(C – L).
dt
a) Sketch the installation.
b) Draw a block diagram of the system.
c) The system is initially in the state: C = 50%, L = 50%, M = 50% and S = 50%.
When time equals zero, the outflow demand L is suddenly increased to 60%.
Determine a time expression for the water level M. (This will demonstrate the
system’s response to a disturbance and is also called regulatory control action).
a)
S=50 Inflow = C
Inflow C Level M Level M
valve controller detector
KP=1 & R=50
Water supply
Outflow demand = L
b)
Water container M
S=50 E=50-M Control unit C Control C
dM
C=E+50 valve =0.002(C–L)
dt
L
M Controller
Sensor
dM
∴ +0.002M = 0.08 and with M(0) = 50 gives the solution M = 40 + 10e-t/500
dt
A graph of the response of the water level M to the disturbance, is shown below.
Demand L
60 %
50 %
Time (sec.)
Water level M
M = 40 + 10e-t/500
50 %
43.7% Offset
40 %
Time (sec.)
t=0 t=500 sec.
The behaviour of the water level may be compared with our findings in Example 6-5.
The graphs from point B onwards, suggest a possible procedure with which we
could eliminate the offset (perhaps someone gave us a hosepipe with running water
and the task to restore the water level to 50 %). Firstly, we must keep the total inflow
more than 60 % so that the level starts returning to 50 %. While we are doing this, we
must remember that the error is getting smaller and we will get less and less assistance
from the proportional controller. Secondly, after the level is restored to 50 %, the
inflow from the proportional controller will only be 50 % and we need to supply the
extra 10 % to equal the demand of 60 % so that the level remains at 50 %.
Integral action gives an output proportional to the time integral of the error. This
action can accomplish the two tasks mentioned above, namely, restoring the measured
value to set point and providing the extra bias to match the demand. Integral control is
rarely used alone but proportional control together with integral control, or PI control,
is widely used in industry because it provides the important practical advantage of
eliminating the offset. The controller output C, of a PI controller, can be expressed as:
∫
C = KPE + KI Edt + R, ……………….……………………… Equation 6-4
where KP is the proportional gain, KI is defined as the integral gain and R is the bias.
(Note: For PI control mode the bias term is optional and may be disabled.)
∫
I=0.0005 Edt
10%
0% Time
Bias R Bias R
50% 50%
Time Time
Controller output Controller output
C=E+R C=E+I+R
60% 60%
50% 50%
Time Time
Figure 6-15 P and PI control comparison
Repeat time and reset time
Many commercial PI controllers do not directly allow for the adjustment of the
integral gain KI, but rather the adjustment of the repeat time or reset time. Imagine an
error signal that suddenly changes from 0 to 10 % and remains constant after that, as
shown in Figure 6-16.
Error E
Area = 10×TR
10%
0% Time
∫ ∫
PI controller output
Integral effect = KI Edt = KI 10dt = 10KITR
∫
C=KPE+KI Edt
Proportional effect = KP×E = KP×10
50%
Time
TR Figure 6-16
(Repeat or reset time) Reset time
If this step in error signal is presented to the input of a PI controller, that initially
EIPINI Chapter 6: Process Control Page 6-18
50
65
49
48 60
47
55
46
45 50
0 500 1000 1500 2000 2500 3000 3500 4000 t (sec) 0 500 1000 1500 2000 2500 3000 3500 4000 t (sec)
a) Water level Figure 6-17 b) Controller output
System response to demand change from 50% to 60% with PI control (KP=1.25, KI =0.002083)
contribution of the proportional term to equal the derivative contribution is called the
rate time, and it gives an indication of how much faster the derivative action is
compared to the proportional action.
It is immediately clear from Figure 6-18 that 10×KITR = 10×KP = 10×KD/TD or
KPTD = KD. The rate time is therefore given by:
K
TD = D or KD = KPTD ……..……….…..………………….. Equation 6-7
KP
∫
dE
C = KPE + KI Edt + KD + R. …………….……………… Equation 6-8
dt
Controller tuning
Adding integral and derivative action to a controller increases the complexity of the
system. Choosing suitable P, I and D gains for a PID controller (tuning the controller),
becomes complicated. Even for the very simple water level control system in Figure
6-9, analysis of the complete PID system is difficult. We will however conclude with
one example of PID control of the water level control system. The controller gains
chosen were KP = 5, KI = 0.08 (reset time = 1.04 min.) and KD = 500 (rate time = 1.7
min.). The plots in Figure 6-19, were obtained with Matlab. The response curves in
Figure 6-19, show the water level (although oscillating), not even dropping below
49% and returning back to setpoint in less than seven minutes.
M(%) 50.6 C(%) 65
50.4
50.2
50
60
49.8
49.6
49.4
49.2 55
0 500 1000 1500 2000 2500 3000 3500 4000 t (sec) 0 500 1000 1500 2000 2500 3000 3500 4000 t (sec)
a) Water level b) Controller output
System response to demand change from 50 % to 60 %
Figure 6-19 with PID control (KP = 5, KI = 0.08 and KD = 500)
EIPINI Chapter 6: Process Control Page 6-21
Controller
Pilot
output C
Set point S relay
(20-100 kPa)
50% = 60 kPa Restriction Air
supply
Proportional
S bellows C
Beam
C Controller
Pilot
output C
Set point S relay
(20-100 kPa)
50% = 60 kPa Restriction Air
supply
Proportional
S bellows C
Beam
The restriction will prevent the pressure R in the reset bellows to immediately
follow the output pressure C and causes a delay in the pressure build up (or decay)
in the reset bellows which will depend on the needle valve setting. This delay in
reset pressure will cause integral action to be added to the controller and it can be
∫
a ra
KP = and K = , this equation then becomes C = K E + K Edt .
I P I
b b
The controller thus provides a proportional plus integral output. The dynamic
integral term effectively replaces the fixed bias with an ‘automatic bias’. Commercial
controllers do however allow a fixed bias to be enabled for PI control, if so required.
Again it can be shown that for the configuration in Figure 6-23, the controller output
a a d(S - M) a a
is given by C = (S – M) + + R. Using E = S - M, KP = and KD = ,
b bd dt b bd
dE
this equation then becomes C = KPE + KD + R.
dt
It is now clear that inserting a delay in the proportional pressure line, causes
derivative action to be added to the proportional action.
Adding delay actions simultaneously in both the reset and the feedback line, will
result in a complete PID controller and such a structure is shown in Appendix 6-2.
EIPINI Chapter 6: Process Control Page 6-24
Proportional
RPF action
E RPI VP 100k
Op-amp
R
VP = - PF E
Gnd R
PI
= -KPE
Error
100k Integral
100k
M action
CI Adder
E E
Op-amp 100k
RI VI 100k
S 100k E=S-M Op-amp
Op-amp
C
100k 1
Gnd Gnd
VI= -
R C ∫Edt
I I
∫
= -KI Edt Gnd
C = - VP – VI – VD
∫
dE
=KPE + KI Edt + KD
dt
Derivative
RD action P I D
E
CD VD 100k
Op-amp
dE
VD= -RDCD
Gnd dt Figure 6-24
dE
= -KD Electronic PID
dt
controller
EIPINI Chapter 6: Process Control Page 6-25
For instance, referring to Figure 6-26, at time instant 4, just after M4 was sampled,
the error E4 may be obtained from 50–M4, given that the set point is 50%. From E4 the
proportional action can be calculated as P4 = KP×E4. If we assume that the controller
started its operation at time equals zero, then a very simple method to calculate the
integral component I4 is to approximate the area A, under the error curve with
4×T
rectangles to obtain A= ∫ S Edt ≈ T ×E + T ×E + T ×E + T ×E . The integral
s 0 s 1 s 2 s 3
0
action I4 may then be expressed as I4 = KI×(Ts×E0 + Ts×E1 + Ts×E2 + Ts×E3) or
I4=KI×Ts×(E0+ E1+E2+E3). The derivative component D4, depends on the slope of the
E −E
error curve at time instant 4, and may be approximated with D4 = KD 4 3
T
s
M E E
S=50%
E3
E0 E2
E1 E1 E3×Ts
E2 E4
M0 E3 E4 E2×Ts E3
E0 E1×Ts
M1
M2 M3 E0×Ts
M4
0 1 2 3 4 Time 0 1 2 3 4 Time 0 1 2 3 4 Time
Ts Ts Ts Ts Ts Ts Ts Ts Ts Ts Ts Ts
E −E
C4 = KPE4 + KI×TS ×(E0 + E1 + E2 + E3) + KD× 4 3
Ts
Figure 6-26 P I D
EIPINI Chapter 6: Process Control Page 6-26
6.4.1 Actuators
Actuators may be categorized
as ‘direct acting’ or ‘reverse
acting’ and some configurations
are indicated in Figure 6-27. In
a reverse acting actuator an
increase in the pneumatic
pressure applied to the
diaphragm lifts the valve stem
(in a normally seated valve this
will open the valve and is
called ‘air to open’). In a direct- Reverse acting
acting actuator, an increase in
the pneumatic pressure applied
to the diaphragm extends the
valve stem (for a normally
seated valve this will close the
valve and is called ‘air to
close’). The choice of valve
action is dictated by safety
considerations. In one case it
may be desirable to have the
valve fail fully open when the
pneumatic supply fails. In
another application it may be
Direct acting
considered better if the valve
fails fully shut.
Figure 6-27
EIPINI Chapter 6: Process Control Page 6-27
Spring
Diaphragm
plate
Actuator
(Motor)
Diaphragm
Air Stem connector
supply (stroke adjust
20 – 100 kPa)
Open
Output ½
Travel
indicator
Close
Stem
Valve
positioner Connector Yoke
arm
Gland and
packing
Instrument
signal Bonnet (yoke
(20-100 kPa) hold down) nut
Bonnet
Gasket
Plug
Valve
body
Seat
Ball valve
Gate valve
1. Globe valves
Globe valves are the most widely used in industry for flow control in both on/off and
throttling service. They typically have rounded bodies from
which their name is derived. They are linear motion valves
with a tapered plug or disk attached to the stem that closes onto
a seating surface to act as flow control element, as shown in
Figure 6-31. The liquid must make two 90° turns when passing
through the valve and because of that, the pressure drop in the
globe valve is significant, even when fully open. The high Figure 6-31
pressure drop is the main disadvantage of the globe valve.
Globe valve
2. Gate valves
Gate valves (also known as knife valves or slide valves) are the most common valve
used for on/off service but could be used for throttling. They are linear motion valves
and a flat disk or wedge slides into the flow stream
Parallel
to act as flow control element, as shown in Figure disk
6-32. The direction of fluid flow is not changed by
the valve. When fully open, the wedge completely
clears the flow path creating minimum pressure drop.
Gate valves are advantageous in applications
involving slurries, as their ‘gates’ can cut right Wedge
through the slurry. They are also used in applications Figure 6-32
that involve viscous liquids such as heavy oils, light
Knife Gate valve
grease, varnish, molasses, honey and cream.
3. Needle valves
A needle valve, shown in Figure 6-33, is used to make relatively
fine adjustments to the fluid flow and can be used for on/off and
throttling service. They are linear motion valves with a tapered
‘needle like’ cone shaped plug that acts as the flow control
element. The fluid enters from below into the seat which forms
part of the flow path. The needle plug permits the flow opening
in this channel to be increased or decreased very gradually.
Needle valves are widely used for steam, air, gas, water or other
Figure 6-33
non-viscous liquids. Disadvantages of the needle valve are a
large pressure loss and possible clogging of the flow orifice. Needle valve
4. Pinch valves
The relatively inexpensive pinch valve, illustrated
In Figure 6-34, is the simplest of all the valve
designs and may be used for on/off and throttling
service. Pinch valves are linear motion valves
that use a flexible tube or sleeve to effect valve
closure. Pinch valves are ideally suited for the
handling of slurries, liquids with large amounts Figure 6-34
of suspended solids and corrosive chemicals. Pinch valve
EIPINI Chapter 6: Process Control Page 6-30
5. Diaphragm valves
A diaphragm valve, illustrated in Figure 6-35, is related to pinch
valves and one of the oldest types of valve known. Leather
diaphragm valves were used by Greeks and Romans to control
the temperature of their hot baths. Diaphragm valves use linear
stem movement for on/off and throttling and are excellent for
controlling fluid flow containing suspended solids. A stud moves
a flexible and resilient diaphragm into the flow thereby acting as Figure 6-35
flow control element. The diaphragm valve is used primarily for Diaphragm valve
handling viscous and corrosive fluids as well as slurries.
6. Plug valves
Plug valves, illustrated
in Figure 6-36, are also
called plug cock or stop
cock valves and they
also date back to ancient Close Throttle Open
times, when used by Romans in plumbing systems. Today they Figure 6-36
remain one of the most widely used valves for both on/off and Plug valve
throttling services. A plug valve is a rotary moving valve that uses a cylindrical or
tapered plug as the flow control element. The opening through the plug may be
rectangular or round. Flow is regulated from closed to open during a 90° turn. If the
opening is the same size or larger than the pipe’s inside diameter, it is referred to as a
full port otherwise as standard round port. Full port plug valves offer little resistance
to flow when fully open, resulting in small pressure loss. A disadvantage off plug
valves is its poor throttling characteristics.
7. Ball valves
Ball valves, illustrated in Figure 6-37, are related to plug valves
and are used in situations where tight shut-off is required. Ball
valves are rotational motion valves and are used for on/off and
throttling service. The flow control element is a sphere with a
round opening rotating in a spherical seat. Flow is regulated
from closed to open during a 90° turn. Full port ball valves Figure 6-37
create a minimum pressure loss when fully open. A disadvantage
of the ball valve is its poor throttling characteristics. Ball valve
8. Butterfly valves
A butterfly valve, illustrated in Figure 6-38, is a rotary
movement valve used for on/off flow control and especially
in throttling applications. The flow control element in
butterfly valves is a circular disk with its pivot axis at right
angles to the direction of flow. A 90° turn of the axis moves
the valve from closed to open. Butterfly valves are well
suited to handle large flows of liquids or gasses as well as Figure 6-38
slurries and liquids with large amounts of suspended solids. Butterfly valve
EIPINI Chapter 6: Process Control Page 6-31
Δp Δp Δp Δp Δp
80
60
40
20
Δp Δp Δp Δp Δp
0 psi 0 psi 0 psi 0 psi 0 psi
4 4 4 4 4
Figure 6-40
EIPINI Chapter 6: Process Control Page 6-32
The special value of C that corresponds to the valve fully open, is called the
characteristic flow coefficient of the valve and denoted by CV. The characteristic flow
coefficient CV, is an extremely important parameter of a specific valve, and is used
extensively when choosing the correct valve for a certain application. For the valve
characteristics in Figure 6-40, for example, CV = 40.
The way q varies when the pressure is kept constant (normally the rated pressure
drop across the valve for maximum flow) and the way C changes as the valve opening
x changes, may be available for a particular valve, provided by the manufacturer in
tabulated format. For example, the graph of q (with Δp = 4 psi) as a function of x and
the graph of C as a function of x, are shown in Figure 6-41 (a) and (b) respectively, for
the flow characteristics of our hypothetical valve, depicted in Figure 6-40.
Figure 6-41(a) Figure 6-41(b)
q (gpm) C
80 40
60 30
40 20
Valve Valve
20 10
opening opening
x x
0 0 1
¼ ½ ¾ 1 ¼ ½ ¾
(0%) (25%) (50%) (75%) (100%) (0%) (25%) (50%) (75%) (100%)
The different values of C obtained as the valve travels through its full range (or
stroke as it is also called) from closed to open, is generally expressed in the following
format, by valve manufacturers:
C = CV×f(x) Equation 6-10
where f(x) is called the inherent valve characteristic of a particular valve. For
example, from the graph of C, in Figure 6-41(b), we could express C as:
C = 40×x,
and we conclude therefore from Equation 6-10, that CV = 40 and f(x) = x.
The function f(x) varies from 0 (valve closed with x equal to 0) to 1 (valve fully
open with x equal to 1). The inherent valve characteristic f(x), is intimately linked to
the flow rate q through a valve, as we can see if we replace C with CVf(x) (as per
Equation 6-10), to rewrite the valve Equation 6-9 in the form:
ΔP
q = CVf(x) Equation 6-11
G
It is indeed clear from Equation 6-11 that if we keep Δp constant (for a given liquid, G
is constant as well), q and f(x) will have exactly the same shape.
EIPINI Chapter 6: Process Control Page 6-33
A valve for which f(x) = x, is called a linear valve because the flow rate will change
linearly with valve opening x. Depending however on the valve design, f(x) may also
reflect quick opening or equal percentage (slow opening) characteristics, as shown in
Figure 6-42.
q f(x) Quick
opening
100% 1
Quick
Graphs also opening Figure 6-42
represent Inherent valve Linear
flow rate q characteristics
Linear
(% of max.
flow) if Δp
is constant. (Ratio of
Equal valve travel Equal
percentage x to maximum percentage
0% 0 valve travel)
0 1
(0% open) (100% open)
The flow behaviour of a valve is dictated by the manner in which the flow areas A1
and A2 changes with valve position and therefore by the style and design of the valve
trimming and in particular the design of the valve seat and closure member (plug). The
quick opening flow characteristic provides for maximum change in flow rate at low
valve travels with a nearly linear relationship. Additional increases in valve travel
gives sharply reduced changes in flow rate. The linear flow characteristic curve allows
the flow rate to be directly proportional to the valve travel (Δq/Δx equals a constant) or
in terms of the inherent valve characteristic, f(x) = x. An equal percentage valve starts
initially with a slow increase in flow rate with valve position which dramatically
increases as the valve opens more.
The term equal percentage for a slow opening characteristic curve may at first be
confused with the description of a linear characteristic curve. However, for an equal
percentage valve, Δq/Δx at any stage, is proportional to the flow rate q at that moment.
This is in contrast with a linear characteristic for which Δq/Δx is constant. That Δq/Δx
is proportional to q, may be rephrased as Δq/q is proportional to Δx. This means the
percentage change Δq with respect to the current flow rate q (that is (Δq/q)×100), is
equal at every valve travel position x for the same change in valve travel Δx, hence the
term ‘equal percentage’. The inherent valve characteristic for an equal percentage
valve is exponential in nature and is normally given by valve manufacturers in the
x-1
form f(x) = R , where R is a constant for the valve. The exponential behaviour of an
equal percentage valve is explored further in Example 6-7.
Example 6-7
An equal percentage valve delivers 8 gpm of water when the valve is 50% open
(x=0.5). When the valve is 60% open (x=0.6), the flow rate increases to 16 gpm.
EIPINI Chapter 6: Process Control Page 6-34
Estimate the flow rate through the valve when it is 70% open (x=0.7). Assume that the
pressure drop across the valve remains constant.
Valve sizing
The characteristic flow coefficient CV of a valve and the inherent valve characteristic
f(x), completely describe the flow-pressure profile of the valve and play an important
role during the valve selection or valve sizing process. Although the characteristic
flow coefficient (CV) is the most fundamental and important parameter influencing the
user when selecting a valve, it is not the only consideration. Conditions such as piping
particulars, fluid type, laminar or turbulent flow, may be included as additional factors
in Equation 6-9. Nevertheless, to select an appropriate valve for a certain
implementation, the user must specify the maximum flow rate required and the
pressure drop expected with the valve fully open. With this information, the necessary
characteristic valve flow coefficient CV, can be calculated from Equation 6-9 by
setting C = CV, q = QRated (rated maximum flow rate with valve fully opened) and
Δp = ΔPFull (pressure across fully opened valve at maximum rated flow rate):
ΔP
CV = Q Full ………………...…………………. Equation 6-12
Rated G
Example 6-8
A system is pumping water from one tank to another through a piping system with
total pressure drop of 150 psi. The maximum design flow rate is 150 gpm. Calculate
the characteristic flow coefficient required for a control valve to be used in this system.
Solution: The usual rule of thumb is that a valve should be designed to use 10-15% of
the total pressure drop or 10 psi, whichever is greater, when fully open. For this
system 10% of the total pressure drop is 15 psi, which is what we must use. From
Equation 6-12 (remember that specific gravity G is the same as relative density or in
other words, G = δ):
CV = QRated / ΔPFull /G = 150/√(15/1) = 38.72 ≈ 39
(Note: CV is normally given as a dimensionless number but it really is a dimensional
quantity with units [gallon-inch/minute-poundforce½] from Equation 6-12.)
Example 6-9
The user in Example 6-8 has a choice between a valve with CV = 35 which would be
too small and the next one with CV = 45 which could be too large. So he decides to
check what the valve travel positions will be for the larger valve if he wants to control
the water flow between a maximum flow rate of 150 gpm and a minimum flow rate of
40 gpm. He expects a pressure drop of 15 psi across the valve at full flow and a rise in
pressure drop to 25 psi, when the flow rate drops to 40 gpm. The valve he intends to
purchase is a linear valve with inherent valve characteristic, f(x) = x.
Solution: From Equation 6-11 (q = CVf(x) ΔP/G ),
with qmax = 150 gpm. and Δpqmax = 15 psi: 150 = 45×x × 15/1
∴x = 150/(45×√15) = 0.8607 (86%)
with qmin = 40 gpm. and Δpqmin = 25 psi: 40 = 45×x × 25/1
∴x = 40/(45×√25) = 0.1778 (18%)
This seems to be good enough, as the aim should be to keep valve movement between
20% and 80% of maximum travel.
EIPINI Chapter 6: Process Control Page 6-36
Valve actuator
Elastic Cam
Force-
balance Valve stem
beam
Figure 6-43 Valve positioner
Pivot Instrument
bellows Controller signal
(Instrument signal) Actuator
output
20-100 kPa
Pilot
Flapper and
relay
nozzle
Air supply
Restriction
Figure 6-43 shows a typical pneumatic valve positioner, assisting a reverse acting
valve. The operation revolves around a force balance flapper and nozzle
arrangement. The valve stem position is communicated to the force balance beam
by means of a lever and cam. If the controller signal increases, with the intention to
open the valve more, the instrument bellows will push the flapper (a flexible plate)
towards the nozzle. This will increase the valve actuator pressure and the valve
stem will start to move upwards. This will turn the cam which will subsequently
push the flapper away from the nozzle, until balance is reached so that the valve
remains in the new more opened position. If the controller signal decreases, the
flapper will move away from the nozzle and the actuator pressure will be reduced.
The valve will begin to close and the valve stem will begin to move downwards,
relieving cam pressure on the beam that kept the flapper away from the nozzle. The
flapper will tend to move towards the nozzle thereby restoring balance with the
valve more closed than before.
EIPINI Chapter 6: Process Control Page 6-37
APPENDIX 6-1
A water container is 5 meter high and has a base area of 1 meter2. The rate, measured
in cubic meter per second, at which water flows into the container at any instant, is
QIN. The rate, measured in cubic meter per second, at which water flows out of the
container at any instant, is QOUT. The level, measured in meter, of the water in the
container at any instant, is H. The maximum value of QIN and QOUT, is 0.01 cubic
meter per second, and the maximum value of H is 5 meter.
a) Express H in terms of a percentage M, of the maximum level.
b) Express QIN in terms of a percentage C, of the maximum inflow.
c) Express QOUT in terms of a percentage L, of the maximum outflow.
d) Calculate the volume V, of the water in the container, in terms of M.
e) Calculate the resultant rate Q, at which water accumulates into the tank or drains
out of the container. Express your answer in terms of C and L.
ΔV
f) Calculate the rate , at which the volume changes, in terms of C and L.
Δt
ΔM
g) Calculate the rate , at which the water level changes, in terms of C and L.
Δt
C (inflow in %)
Solution: 0.01m3/s
Qin (inflow in m3/s)
M Max
a) H = ×5 = 0.05×M meter …...……… (1)
100
C 5m max
3
b) Qin = ×0.01 = 0.0001×C m /s ….….. (2) M (%)
100 H (m)
L 3
c) Qout = ×0.01 = 0.0001×L m /s ….…. (3) L (outflow in %) 0.01m3/s
100 3
Qout (outflow in m /s) max
d) Volume V, of water in tank:
V = Height×Area = H×1
Q=0.0001×(C-L)
From (1): V = 0.05×M×1 C
∴V = 0.05×M cubic meter ……..…..….. (4)
e) Resultant rate Q = Qin - Qout
M V=0.05M
From (2) and (3):
∴Q = 0.0001×C–0.0001×L
∴Q = 0.0001(C – L) m3/s ….....……....... (5) Area = 1 m2 L
f) The rate at which the water volume in the container increases or decreases, is the
same as the rate at which water is accumulated in or drained from the tank.
ΔV ΔV
∴ = Q. Therefore from (5): = 0.0001(C - L) cubic meter per second.
Δt Δt
∴ΔV = 0.0001(C - L)×Δt meter3 ….….... (6)
g) From (4): V = 0.05×M
∴ΔV = 0.05×ΔM meter3 ……….…......... (7)
Substituting ΔV from (7) into (6): 0.05×ΔM = 0.0001(C – L)×Δt
ΔM 0.0001 dM ΔM
∴ = ×(C – L), or, using the proper notation for ,
Δt 0.05 dt Δt
dM
= 0.002×(C – L) percent per second.
dt
EIPINI Chapter 6: Process Control Page 6-38
APPENDIX 6-2
Pd C
a(r + d) a d(S - M)
∫
ar
C= (S-M) + (S - M)dt (S-M) +
bd b bd dt
a(r + d) ar a
Using E = S-M, KP = , KI = and KD = :
bd b bd
dE
∫
C = KPE + KI Edt + KD
dt