WOLFRAM|DEMONSTRATIONS PROJECT

Controlling Liquid Height within Two Tanks in Series Using a PID Controller

​
proportional gain
2.4
integral time constant
1.5
differential time constant
0.1
Consider two tanks in series; the dynamic behavior of their heights is governed by the ODEs:
​​
A
d
h
1
dt
=
F
0
-
F
1
,
A
d
h
2
dt
=
F
1
-
F
2
,
where
h
i
is the height of tank
i
in meters,
A
(set to 1 here) is the area of each tank in
2
m
,
F
0
is the inlet flow rate (the manipulated variable) into the first tank, and
F
1
and
F
2
are the outlet flow rates (expressed in
3
m
/s
) from the first and second tank, respectively. The initial height of the liquid in each tank is assumed to be 0.25 meters.
The flow equation is given by
F
i
=
K
i
h
i
where
i=1
or
2
and
K
i
is the valve constant expressed in
3
m
/(s
0.5
m
). It is assumed that
K
1
=2.0
and
K
2
=1.5
.
The setpoint for the height of the liquid (the process variable) in the second tank is chosen to be 3 meters.
The inlet flow rate is varied in order to achieve the desired setpoint value using P, PI, or PID (proportional–integral–derivative) control:
F
0
=1+
K
p
e(t)+
1
τ
i
e(t)dt+
τ
d
de(t)
dt
, where
e(t)=(3-h(t))
is the error,
K
p
is the proportional gain, and
τ
i
and
τ
d
are the integral and differential time constants, respectively.
For very large
τ
i
and
τ
d
=0
, one recovers the usual proportional control, which is usually characterized by a small offset value (i.e., the final steady state height is not exactly equal to the setpoint value).
PI control is achieved when
τ
d
is taken equal to zero. PI control can show an overshoot and dumped oscillations around the setpoint. No offset is observed and the final steady state tank height is equal exactly to the setpoint value.
The most general case is when
τ
d
≠0
and
τ
i
is not too large; one gets PID control of the second tank's height.
The Demonstration plots the height of both tanks as a function of time; the blue curve displays
h
1
(t)
and the magenta curve represents
h
2
(t)
.