Assignment 2
Assignment 2
CONTINUOUS ASSESSMENT
Assignment: Part 2
Fi (L/min), Ti (oC)
Thermocouple
_ T
Q
TSP +
F (L/min), T (oC)
e
Condensate
Integral control
Fst (kg/s)
Steam
Ti
Ti,S
T0 Time (s)
1
Table 1 STH operating conditions and parameters
operating conditions
Inlet temperature 20 oC
Process parameters
During operation, some process varibles (disturbances) may change and affect operation of the
STH. One such variable is the feed temperature. For a particular case where feed temperature
suddenly increases from steady state 𝑇𝑖𝑠 = 20 oC to 50 oC, it is desired to control the STH
temperature. This can be achieved by performing the following.
a) Predicting the temperature response when the STH is controlled using a proportional-
integral control.
b) For a given integral-time constant, investigate the effect of the process time constant on
behaviour of STH.
c) For a given time constant, investigate the effect of the controller tuning parameter
‘controller gain and integral time’ on the STH’s response.
Assignment guidelines
The format should follow that of a technical report (Introduction, problem statement,
hypothesis, objectives, methodology, results & discussion, conclusion).
The objectives should be in line with what is given above in the assignment.
Additional MATLAB code is provided for for I-only controlled STH. The code can be run
in MATBLAB to get a general idea of how MATLAB works.
2
MATLAB software can be downloaded from MathWorks MATLAB R2021a Free
Download (getintopc.com) and a video containing instruction on how to install can be
found from the same website or ‘ReadMe Win’ in the downloaded MATLAB folder.
Integral control
𝑑(𝑇 − 𝑇𝑠 )
𝑉𝜌𝐶𝑝 = 𝐹𝜌𝐶𝑝 [(𝑇𝑖 − 𝑇𝑖,𝑠 ) − (𝑇 − 𝑇𝑠 )] + 𝑄 − 𝑄𝑠 (1)
𝑑𝑡
𝐾𝑐
𝑄 = 𝑄𝑠 + ∫ 𝑒(𝑡)𝑑𝑡 (2)
𝜏𝐼
Where 𝑒(𝑡) = 𝑇𝑠 − 𝑇
𝑑𝑥
𝑉𝜌𝐶𝑝 = 𝐹𝜌𝐶𝑝 [(𝑇𝑖 − 𝑇𝑖,𝑠 ) − 𝑥] + 𝑄 − 𝑄𝑠 (3)
𝑑𝑡
𝑑𝑥 1 𝑄 − 𝑄𝑠
= [(𝑇𝑖 − 𝑇𝑖,𝑠 ) − 𝑥] + (4)
𝑑𝑡 𝜏 𝜌𝑉𝐶𝑝
Or
𝑑𝑥 1 1 𝑄 − 𝑄𝑠
= (𝑇𝑖 − 𝑇𝑖,𝑠 ) − 𝑥 + (5)
𝑑𝑡 𝜏 𝜏 𝜌𝑉𝐶𝑝
Where 𝜏 = 𝑉⁄𝐹
Since 𝑒(𝑡) = 𝑇𝑠 − 𝑇 = −(𝑇 − 𝑇𝑠 ) = −𝑥, then the integral control law can also be written as
𝐾𝑐
𝑄 = 𝑄𝑠 − ∫ 𝑥𝑑𝑡 (6)
𝜏𝐼
And
𝐾𝑐
𝑄 − 𝑄𝑠 = − ∫ 𝑥𝑑𝑡 (7)
𝜏𝐼
3
Let 𝑦 = 𝑄 − 𝑄𝑠 = − ∫ 𝑥𝑑𝑡 then equations (5) and (7) become
𝑑𝑥 1 1 𝑦
= (𝑇𝑖 − 𝑇𝑖,𝑠 ) − 𝑥 + (8)
𝑑𝑡 𝜏 𝜏 𝜌𝑉𝐶𝑝
𝐾𝑐
𝑦=− ∫ 𝑥𝑑𝑡 (9)
𝜏𝐼
Differentiating (9) with respect to t
𝑑𝑦 𝐾𝑐
=− 𝑥 (10)
𝑑𝑡 𝜏𝐼
Equation (8) and (10) are the two differential equation in 𝑥 and 𝑦 and can be solved simultaneously
using MATLAB.
%differential equations
4
dy1dt = -Kc*T1/I1;
Script file
tspan = [0 10000];
x0 = [0 0 0 0 0 0]; % initial values of T and Heat supply
[t,x] = ode45(@I_ctrl_STH,tspan,x0);
p1 = x(:,1);
p2 = x(:,3);
p3 = x(:,5);
plot(t,p1,'r','linewidth',1.5);
hold on
plot(t,p2,'g','linewidth',1.5);
hold on
plot(t,p3,'b','linewidth',1.5);
xlabel('Time (s)','Fontname','times new
roman','fontsize',16,'fontweight','bold');
ylabel('Temperature (oC)','Fontname','times new
roman','fontsize',16,'fontweight','bold');
legend('I = 0.5','I = 2','I = 10','fontname','times new
roman','fontsize',12,'fontweight','bold');
title('load disturbance response for first-oder PI-controlled
STH','fontname','times new roman','fontsize',16);
grid on