Controlling Liquid Height within Two Tanks in Series Using a PID Controller
Controlling Liquid Height within Two Tanks in Series Using a PID Controller
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 is the height of tank in meters, (set to 1 here) is the area of each tank in , is the inlet flow rate (the manipulated variable) into the first tank, and and are the outlet flow rates (expressed in /s) from the first and second tank, respectively. The initial height of the liquid in each tank is assumed to be 0.25 meters.
h
i
i
A
2
m
F
0
F
1
F
2
3
m
The flow equation is given by = where or and is the valve constant expressed in /(s). It is assumed that =2.0 and =1.5.
F
i
K
i
h
i
i=1
2
K
i
3
m
0.5
m
K
1
K
2
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: =1+e(t)+e(t)dt+, where is the error, is the proportional gain, and and are the integral and differential time constants, respectively.
F
0
K
p
1
τ
i
τ
d
de(t)
dt
e(t)=(3-h(t))
K
p
τ
i
τ
d
For very large and =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).
τ
i
τ
d
PI control is achieved when 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.
τ
d
The most general case is when ≠0and is not too large; one gets PID control of the second tank's height.
τ
d
τ
i
The Demonstration plots the height of both tanks as a function of time; the blue curve displays (t) and the magenta curve represents (t).
h
1
h
2