Dec 2017 PDC Solutions
Dec 2017 PDC Solutions
Q.1.a) A liquid surge tank has the following transfer function of [05]
H ( s) 10
Qi ( s ) 50s 1
The system is operating at the steady state with qis=0.4m3/s, and hs=4m when the inlet flow rate
fluctuates as a sine wave with an amplitude of 0.1m3/s and a period of 500 sec. what is the maximum
and minimum value of the level after 10 min?
2t
Qi (t ) 0.1sin Qi (t ) 0.1sin0.01256t
500
1.256 * 10 3
Qi s
0.1 * 0.01256
s 2 0.012562 s 2 1.577536* 10 4
10 10 1.256 *103
H ( s) * Q ( s) *
50s 1 i 50s 1 s 2 1.577536*104
We know the Laplace inverse of first order process for sine input as t=600 secs the first part of inverse
can be neglected.
10 * 0.1
H (t ) sin(0.01256t )
1.577536*104 * 502 1
10 * 0.1
A 0.846854
4
1.577536*10 * 50 1 2
The maximum and minimum value can be found as The given height is 4m at steady state then
Maximum value=Hs+A=4+0.846854=4.846854 m
Minimum Value=Hs-A=4-0.846854=3.153146 m
Q.1.b) A second order system is found to have a peak amplitude ratio of 1.1547 at a frequency of
0.7071 rad/min. What are the values of time constant and the damping coefficient of the system? [05]
Overshoot=0.1547
overshoot exp
1 2
2
0.35897
1 2
0.85974
2 2
1 2
0.722sec
Q.1.c) The following reaction takes place in a CSTR at a constant temperature
2 A
k1
B rA k1C A2
where CA is the concentration of A in the reactor. Derive the transfer function relating the outlet
concentration CA to the inlet concentration CAi. Assume volume is constant. [05]
V
dC A
dt
F C Ai C A k1C A2V
Has we have nonlinear term we have to linearize the reaction term in the above balance equation
By Taylor series
V
dC A
dt
F C Ai C A k1C A2,sV 2k1VC A,s C A
V
dC A, s
dt
F C Ai , s C A, s k1C A2 , sV 2k1VC A, s C A, s
Divide by V
dC A F
dt
V
C Ai C A 2k1C A, s C A
dC A F F
C A,i (2k1C A, s )C A
dt V V
dC A F F
C A,i (2k1C A, s )C A
dt V V
F
CA V
C A ,i
s ( 2 k1 C A , s F )
V
Q.1.d) If the input ym to a PI controller changes stepwise (ym(s) = 2/s) and the controller output changes
initially as in figure below, what are the values of the controller gain and integral time? [05]
Intercept =6
p (t ) 6 1.2t
6 1.2
p ( s )
s s2
2 0 .6
p ( s ) 3
s s
Given that 2/s is ym(s)
p ( s) 0.6
3
y m ( s) s
p ( s) 0.2
31
y m ( s) s
p ( s) 1
31
y m ( s) 5s
Q.2.a) The variation of liquid level in a spherical tank with the inlet flowrate qi and outlet discharging
through the valve can be described as
dh 1
qi C v h
dt D h h
Derive the transfer function relating the changes in the liquid level h to the changes in the inlet flow
rate qi. The diameter of the tank is D and Cv is the constant of the valve in the outlet line. [10]
dh 1
qi C v h
dt D h h
D 2hs
h hs
1 1 1
D h h D hs hs D hs hs 2
Cv
Cv h Cv hs h hs
2 hs
qi q D 2hs C
i h hs Cv hs v h hs
dh
dt D hs hs 2
D hs hs 2 hs
dhs qi s q D 2hs C
i s hs hs Cv hs v hs hs
dt D hs hs 2
D hs hs 2 hs
1
K1
D hs hs
qi s D 2hs C
K2 v
D hs hs 2 hs
2
dH H ( s) K1
K1Qi K 2 H
dt Q1 ( s) s K 2
2.5K c
GOL
s s2
2
for what value of Kc will the point -1 be encircled? Will it be clockwise direction? Will the closed
loop system be stable?
2.5K c
G|OL
s 2s s 2
2
1.25K c
G|OL
( s 1)(0.5s 1)
1.25K c
The amplitude ratio is given by A.R
2 1 0.25 2 1
A.R 1.25
phase angle tan 1 ( ) tan 1 (0.5 )
Kc 1 0.25 1
2 2
We get ωco=0
Substituting in A.R/Kc equation we known A.R=1 at cross over frequency thus Kc=0.8
0.2
0.15
0.1
0.05
Imaginary Axis
-0.05
-0.1
-0.15
-0.2
-0.25
-1 -0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 0
Real Axis
Q.3.a) A stirred-tank reactor has an internal cooling coil to remove heat liberated in the reaction. A
proportional controller is used to regulate coolant flow rate so as to keep the reactor temperature
reasonably constant. The controller has been designed so that the controlled reactor exhibits typical
underdamped second-order temperature response characteristics when it is disturbed, either by feed
flow rate or by coolant temperature changes.
(a) The feed flow rate to the reactor changes suddenly from 0.4 to 0.5 kg/s, and the temperature of the
reactor contents, initially at 100 °C, changes eventually to 1020C. What is the gain of the transfer
function (under feedback control) that relates changes in reactor temperature to changes in feed flow
rate? (Be sure to specify the units.).
(b) The operator notes that the resulting response is slightly oscillatory with maxima estimated to be
102.50C and 102 0C occurring at times 1000 and 3060 s after the change is initiated. What is the
complete process transfer function?
(c) The operator failed to note the rise time. Predict t, based on the results in (a) and (b). [10]
Solution
(a) The gain is obtained by dividing the steady-state change in temperature by the feed flow rate
(disturbance) change:
102 100 0
C
K 20
0 .5 0 .4 kg s
(b) The oscillatory characteristics of the response can be used to find the dynamic elements in the
transfer function relating temperature to feed flow rate. Assuming the step response is due to an
underdamped secondorder process, Figs. 5.8 and 5.11 can be used to obtain estimates of ' and T.
Alternatively, analytical expressions can be used, which is the approach taken here.
Overshoot can be calculated as
ln OS 2
0.4
2 ln OS 2
1 2
P 300s
2
tr ( cos1 )
1 2
tr=649 s
T ( s ) 20
where the process gain has unit 0C/kg/s
W ( s ) 300 s 2(0.4)(300) s 1
2 2
Q.3.b) Two streams w1 and w2 each at a constant density of 900 kg/m3, and carrying solute of mass
fraction x1 and x2 respectively, enter a continuous stirred tank of 2m3 capacity. At steady-state,
w1s=500 kg/min, w2s=200 kg/min, x1s=0.4, and x2s=0.75. Suddenly the inlet flow rate w2 decreases to
100 kg/min and remains there. Determine an expression for the mass fraction of the solute x(t).
Assume that liquid hold up is constant
w1 w2 w
dV 1
dt
dVx
w1 x1 w2 x 2 wx
w1 x1 w2 x2
At steady state condition x
dt w
Substituting value of w
dx w1
x1 x w2 x 2 x
dt V V
dx w1 w2 wx w x
x 1 1 2 2
dt V V V
V dx w1 x1 w2 x2
x
w dt w w
dx wx w x
3 x 1 1 2 2
dt w w
dx
3 x 0.45833
dt
t t
x(t ) 0.5e 3
0.45833(1 e )
3
80
Q.4.a) A unit feedback control system has: G ( s ) Generate Bode plot and
s( s 2)(s 20)
comment on the stability.
Bode plot
0
10
-1
10
AR/K
-2
10
-3
10
-2 -1 0 1 2
10 10 10 10 10
0
phaseangle(degrees)
-50
-100
-150
-200
-250
-2 -1 0 1 2
10 10 10 10 10
frequency (rad/time)
A.R/K PHASE ANGLE (ϕ)0C Frequency (ω)
2.5000 -90.3151 0.0100
0.8984 -90.8768 0.0278
0.3226 -92.4388 0.0774
0.1154 -96.7655 0.2154
0.039 -108.4026 0.5995
0.0115 -134.5975 1.6681
0.0021 -169.7553 4.6416
0.0002 -204.050 12.9155
0.0000 -237.7182 35.9381
0.0000 -257.5443 100.0000
5e 0.2 s
G p (s)
2s 2 s 1
Design a PI controller for the negative feedback loop of the process, based on the Zeigler and Nicholos
tuning rules.
5
A.R
1 2
2 2
2
tan 1 2
0.2
1 2
We have to determine the cross frequency and ultimate gain by substituting ϕ=-180 in the above
equation
5
The A.R 1.0445
20.169 2.7455
Q.5.a) Using Routh Criteria determine the limits of Kc for the stability of the system with following
Kc
open loop transfer function. G ( s ) [10]
s 1(0.5s 1)( s 3 1)
Solution
Kc
G( s)
s 1(0.5s 1)( s 3 1)
The closed loop transfer function is given by
s 1 K c
0.5 3 5.5
1 G (s) s s2
3 3
Solution
Number of zeroes(m)=0
Center of Gravity
n m
si s j
i 1 j 1
2
nm
s 3 6 s 2 8s K 0
1 8 0
6 K 0
K 48
0 0
6
6s2+K=0
s=0
s 8j
Q.6.a) A continuous stirred-tank heating system shown in the figure below. The liquid inlet stream
consists of a single component with a mass flow rate wi and a inlet temperature Ti. The tank contents
are agitated and heated using a electrical heater that provides a heating rate Q, density and specific heat
capacity of the mixture is constant. Carry out the degrees of freedom analysis, also state the disturbance
and manipulated variables, for following two cases (i) holdup is constant,
(ii) holdup is varying [10]
Solution
Carrying out the general mass balance and energy balance for the given system
dV
wi w
dt
dVT
Cp wi C p Ti wC p T Q
dt
dV wi w
dt
dVT wi Ti wT Q
dt C p
Case(1)
No of parameters (3)=V,ρ,Cp
No of Variables Nv (4)=T,Ti,w,Q
Degrees of freedom=4-1=3
1 manipulated variable Q.
Case(2)
Volume is varying
dV
wi w
dt
dT w Q
i (Ti T )
dt V C pV
No of parameters (02), ρ, Cp
No of equations NE (02)
Degrees of freedom=NV-NE=04
1 manipulated variable Q.
Q.6.b) Consider a liquid phase, irreversible, first order reaction taking place in a CSTR, where the
reactant A gets converted to product B. The reaction is cooled by coolant passing through a coil at a
temperature Tc. Develop a state space model, assuming volume is constant, if both, the concentration
of A and the reaction temperature T, are required to be monitored. Following is the mass & energy
balance equations. [10]
qC Ai C A VkC A
dC A
V
dt
Solution
x Ax Bu Ed
y Cx
dC A
x dt
dT
dt
C
x A
T
dC A q q
C Ai k C A
dt V V
dT qC qC UA ( H R ) UA
Ti T T VkC A Tc
dt VC V C V C V C V C
dT q q UA ( H R ) UA
Ti T T kC A Tc
dt V V VC C VC
dT q q UA (H R ) UA
Ti T kC A Tc
dt V V VC C VC
q
V k 0
A
H R k q UA
C V VC
q
V 0 0
q UA
B 0
V VC
C Ai
u Ti
Tc
E=[0], as d=[0]
C
y A
T
1 0
C
0 1