 Open Access
 Total Downloads : 12
 Authors : A.Kavitha, S.Sutha
 Paper ID : IJERTCONV1IS06137
 Volume & Issue : ICSEM – 2013 (Volume 1 – Issue 06)
 Published (First Online): 30072018
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Reconfigurable controller design for continuous stirred tank reactor

AVITHA S.SUTHA
Department of Instrumentation Engineering, Department of Instrumentation Engineering,
MIT Campus, Anna University, MIT Campus, Anna University,
Chennai, Tamil Nadu, India. Chennai, Tamil Nadu, India. kavi25690@gmail.com
AbstractIn this paper, an optimized Eigenstructure Assignment based reconfigurable Linear Quadratic Regulator (LQR) controller is designed for linear discrete Continuous Stirred Tank Reactor (CSTR). The servo operation is achieved by augmenting the linear discrete state and ProportionalIntegral (PI) control action is incorporated in the design. LQR parameters are optimized using Particle Swarm Optimization (PSO) and Parallel implementation of Genetic Algorithm (PI GA). The closedloop performance of the normal (Faultfree) system obtained using PSO and PIGA are compared. Then, the best optimized technique is used to reconfigure the controller parameters against partial actuator failure in the process.
Key WordsEigenstructure Assignment, Reconfigurable controller, LQR, PSO, PIGA, PI, CSTR.

INTRODUCTION
With the growing complexity of modern engineering and ever increasing demand for safety and reliability there has been great interest in the development of Fault Tolerant Control Systems (FTCS) [1, 2, 6, 11, 13]. Thus FTCS are closedloop control systems which can tolerate system component faults/failures, and are capable of maintaining stability and a certain degree of performance of the system. It can be divided into passive FTCS and active FTCS. Passive methods are essentially robust control techniques which are suitable for dealing with certain types of structure failures that can be modeled as uncertainty regions around a nominal model. In contrast with passive FTCS, active FTCS react to the system component failures actively either by selecting a pre computed control law (projectionbased approaches) or by synthesizing a new control strategy online (online automatic controller redesign approaches [1].
Typically, an active FTCS consists of three parts: a reconfigurable controller, a Fault Detection and Diagnosis scheme, and a control law reconfiguration mechanism. Key issues are how to design: 1) a robust reconfigurable controller,
2) an FDD scheme with high sensitivity to faults and robustness to model uncertainties and external disturbances, and 3) a reconfiguration mechanism which can organize the reconfigured controller such that the eigenstructure of the reconfigured faulty closedloop system should recover its eigenstructure as close as possible to normal closedloop system. Several reconfigurable control techniques have been
proposed in the literatures. However, a very few techniques have dealt with multiinput and multioutput (MIMO) systems. Linear Quadratic Regulator (LQR) and Eigenstructure Assignment (EA) are the popular techniques to design controller for MIMO systems. In general, the LQR based control design guarantees the closedloop stability and certain degrees of robustness, but may not easily achieve specific system performances due to the difficulty in the selection of the weighting matrices. EAbased reconfigurable controller (RC) design can be found in the literatures [4,6].The advantage of EA is that if the specifications are given in terms of system eigenstructure, the eigenstructure can be achieved exactly for the desired stability and dynamic performance. Here EA is used for RC system design while LQR is used for normal system.
This paper focuses on design of an integrated design of reconfigurable control for MultiInput and Multioutput (MIMO) systems. The scheme uses Proportionalintegral (PI) control structure in the reconfigurable controller so as to recover both the dynamic and steadystate performance of the prefault systems. An optimized eigenstructure assignment based reconfigurable state feedback controller is developed to achieve automatic redesign of the controller that compensates partial actuator failures in the Continuous Stirred Tank Reactor (CSTR).Here reconfiguration is carried out against the availability of a perfect Fault Detection and Diagnosis (FDD) and detects the fault after a delay time of 1 sampling instant [1,4,6].
The paper is organized as follows. In Section 2, Continuous Stirred Tank Reactor is described, actuator fault is presented. Section 3 presents the problem formulation. Section 4 and 5 discusses the reconfigurable controller design and simulation respectively.

SYSTEM MODELING
A. Continuous Stirred Tank Reactor
The continuous stirred tank reactor (CSTR), also known as vat or back mix reactor is a common ideal reactor type in chemical engineering. The CSTR process assumed here is the linear time invariant multi input multi output process shown Fig 1. It is nonisothermal first order process with an
irreversible exothermic reaction (A > B) that occurs in constant volume. The heat of the reaction is removed by a
are the state, input, and output vectors respectively where
coolant medium that flows around the reactor. It has two
controlled variable namely outlet concentration (CA) and outlet temperature (T) and manipulated variable are coolant flow rate (qc) and reactant flow rate (q).
[q qC]T [100
103]T l/min and[CA
TABLE 1
T]T [0.0989
438.7763]
Figure 1. CSTR schematic diagram.
The schematic diagram of the CSTR is depicted in Figure 1 as defined by Sentil et al. (2007) [8].The dynamic behavior of CSTR can be described by mass balance and energy balance. Using mass balance and energy balance, the process is modeled by:
Physical Parameters of the CSTR
parameters
value
Feed Concentration (CAO)
1 mol/l
Feed temperature (To)
350 K
Inlet coolant temperature (Tco)
350 K
CSTR volume (V)
100 l
Heat transfer term (hA)
7 x 105 cal/(min K)
Reaction rate constant (ko )
7.2 x 1010 min1
Activation energy term (E/R)
1 x 104 k
Heat of reaction (H)
2 x 105 cal/mol
Liquid density (,c)
1 x 103 g/l
Specific heats (CP ,CPC )
1 cal/( g K)
C. Modelling of faulty system
Modelling of Actuator Faults as Loss in Control Effectiveness Factor:
In this paper, abrupt type bias fault has been
dC A t q C
t C
t k C
t exp E
considered and it is represented by control effectiveness factor
dt v A0 A
0 A RT t
() reduction. During normal operation actuator could deliver
dT t q
Hk
E
control signal without any loss. When fault occurs in the
T 0 t T t
O C A( t ) exp
actuator due to partial blockage or aging, the actuator cannot
dT v
q
c
RT ( t )
deliver control signal without loss. Such a failure condition is represented by reduced control effectiveness factor [11].The

c CC C T CO t T t 1 exp

hA
faulty model given by.
/p>
C v
qcc
xf (k+1) = Af xf (k) + Bf uf (k) (2)
In order to obtain a linear discrete model, the nonlinear system is linearized around the steadystate operating points. The linearization procedure is based on the expansion of the nonlinear function into Taylor series around the operating point and the retention of only the linear term. The physical parameters of the continuous stirred tank reactor are shown in Table 1.
B. Modeling of system
The discrete state space representation for the faultfree model (nominal) is given by
yf(k) = C xf (k)


PROBLEM FORMATION.

Augmented Dynamic System
To solve the command tracking in the presence of noise, disturbance and uncertainty, PI control is incorporated in the proposed scheme [11]. To incorporate the integral action into the controller design, one popular approach is to expand an integrator vector as an augmented state. In the discrete domain, the statespace realization of the integrator vector can be described in:
eI(k 1) eI(k) TSe(k)
(3)
where
x (k+1) = A x (k) + B u (k)
y (k) = C x (k) (1)
Where e(k) is the tracking error, Ts is the sampling time.
The augmented system is derived from Equations (1) and (3) and is represented by
A 0.2248
0.0034
B 0.0007
0.0001 ~ ~
133.3 1.5010 0.0423 0.0934
x (k 1) A x (k) B u(k) Er(k)
y(k) C ~x(k)
x(k) [CA
T]T, u(k) [q qC]T, y(k)[CA T]T
Where
~x(k1)x(k1)
s(k1)
~x(k) x(k)
s(k)
and
A A 0
TC I
performances are provided. The weighting matrices Q an R may be manipulated by natural algorithm in order to force the system to behave in the required manner by the designer, in
B
B
0
E
0
T I
CC 0
0 I
s
this case to force an eigenstructure upon the closed loop system. Provided that the controller is a solution of the Algebraic Ricatti Equation (Eq. (7)), then the performance and
s
The design of state feedback controllers based on the augmented system can allow the computation of proportional and integral controller gains for the system. PI control incorporated Reconfigurable control illustrated in Figure 2.
robustness of the LQR will be maintained. The main aim is to derive a cost function value from the difference between desired and achieved eigenstructure determined by their eigenvalues {di, ai} and their corresponding eigenvectors
{vdi , vai}, i =1,2,..n.
For this, we define the following criteria:

C1 Max i { di ai }, i 1,2,.nthe maximum
norm of the difference between the desired and achieved eigenvalues.
(ii) C2 Max i{ vdi vai 2}, i 1,2,…… .n, the
Euclidean distance between (vai) and (vdi)
Then, objective is to minimize the global criterion Cg defined
by
Figure 2 PI Integrated Reconfigurable control


Mathematical formulation of the ea and the lqr problems
In this work, the natural algorithm (PIGA and PSO) can be used to search the weighting matrices Q and R with Q is a [n x n] matrix and R is a [m x m] matrix. The design problem is described as follows:
The LQR is produced by the minimization of the following cost function:
J = (xT Q x + uT R u) dt (4)
in which Q 0 and R > 0 that defining respectively, the state and the input weighting matrices of the LQR optimization
Cg max{C1, C2}
Thus, a zero cost function value represents the ideal solution. The minimization of the cost function is achieved through manipulation of these Q and R matrices by the natural algorithm. Based on the above criteria, an algorithm for the optimized closed loop eigenstructure assignment can be given.


RECONFIGURABLE CONTROL TECHNIQUE:
The Optimized Eigen structure Assignment based LQR controller discussed in previous section can be extended to recover the performance of the closedloop system, when the system is affected by actuator fault. Faulty system can be recovered by reconfiguring the controller gain of the nominal system [11]. State feedback controller gain (Kf) is determined
by replacing (Af, Bf )with augmented faulty system (A, B ) in the Eq.(1) where
problem. According to the LQR technique, the system (Eq. (1)) is controlled by the state feedback
Af
Af – T C
0 and B
f
I
B
f
0
K
u = Q u (5)
R
s
The closedloop system representation is given by:

Particle Swarm Optimisation Via Lqr Method
X (t) A
C x(t) (A B K Q ) x(t)
(6)
Considering the social behaviour of swarm of fish, bees and other animals, the concept of PSO is developed. PSO is a
R
Via the following Algebraic Ricatti Equation (ARE)
ATP + PA PBR1BTP+Q= 0 (7)
K
We can determine the solution P and the gain Q
R
K
Q= R1BTP (8)
R
A suitable choice of the Q and R leads to the computation of
robust stochastic evolutionary computation method based on the movement of swarms looking for the most fertile feeding location [5].In general, PSO implementation is easier than GA. Indeed, PSO only has one operator; velocity calculation, so the computation time is decreased significantly and has the ability to control convergence through manipulating of the inertial weight.
the gain
Q . Therefore, the eigenvalues and eigenvectors of
K
R
All solutions in PSO can be represented as particles in a swarm. Each particle has a position and velocity vector and
the closedloop system can be chosen and the system
each position coordinate represents a parameter value. Similar to the most optimization techniques, PSO requires a fitness evaluation function relevant to the particles position. XPB and XGB are the personal best (Pbest) position and global best (Gbest) position of the ith particle. Each particle is initialized
The desired eigenvalues and vectors are designed in terms of the time domain specifications .The desired eigenvalues ( di ) and vectors (vdi) are computed as:
0.8590 0.222i
0.8590 0.222i
with a random position and velocity. The velocity of each particle is accelerated toward the global best and its own personal best based on the following equation [8].
di
0.8800
0.900
V i (new ) = w xV
i (old ) + c1 x rand () x (X PB – X i )
0.0049 – 0.0057i
0.0049 + 0.0057i
0.0032
– 0.0017
c 2 x Rand
() x (X GB – X i )
0.0011 – 0.0004i 0.0011 + 0.0004i
0.0037 –
1.0
– 0.1347 –
0.0730i – 0.1347 + 0.0730i
– 0.7568 – 0.
0046
v di
– 0.9882
0.9882
– 0.6536
0.009
Here rand() and Rand() are two random numbers in the range [0,1]; c1 and c2 are the acceleration constants and w is the inertia weight factor. The parameter w helps the particles
converge to Gbest, rather than oscillating around it.In general,
w is set according to the following equation :
w 0.5(1 rand (0,1))
The positions are updated based on their movement which is given as follows:
X i X i V i (6)
Then the fitness at each position is reevaluated. If any fitness is
greater than Gbest, then the new position becomes Gbest, and the particles are accelerated toward the new point. If the particles
p>fitness value is greater than Pbest, then Pbest is replaced by the current position. Here PSO is used to optimize the LQR tuning matrices (Q and R) in order to achieve desired eigenstructure assignment. PSO algorithm parameters are set based on trial and error as follows: Population size = 50 particles, Generation number =100 and acceleration constants: C1=1.5, C2=1.5.

Parallel implementation of genetic algorithm via lqr method
The Genetic Algorithm (G.A) is an optimisation scheme with an unconventional approach and is robust over a large range of nonconvex problems. In normal genetic algorithm, the fitness and constraint functions are evaluated separately at each member of a population. But in Parallelimplementation of genetic algorithm, the fitness and constraint evaluations are done parallel at each member of a population, which occurs once per iteration. The motivation of parallelisation of the algorithm is speedup of processing time [7]. The parameters used in PIGA algorithm are: Population size=50, Generation
number = 50, Crossover rate =0.8, Mutation rate = 0.005,
case i:Normal condition:
Optimized Eigenstructure Assignment based LQR controller using Particle Swarm Optimization and Parallel implementationGenetic Algorithm is used to control the normal system. The achieved eigenvalues (ai) and eigenvectors (vai) obtained through Particle Swarm Optimization and Parallel implementationGenetic Algorithm are given Table 1.
TABLE 2. Performance comparison between PIGA and PSO
LQR – PIGA
LQR – PSO
ISE
1.161*1005
7.7202*1006
CONTROL
EFFORT
9.976*e+007
9.978*e+007
Case II Abnormal condition:
Here control effectiveness factor (30%) is taken as =0.7.From Table 2, the closedloop performance of the normal system obtained using PSO and PIGA: PSO gives better performance than PIGA. So the PSO is used for recovery of the abnormal system. The state feedback gain matrix of the abnormal system is same as normal system.
Case III Recovered condition:
Optimized Eigenstructure Assignment based LQR controller using Particle Swarm Optimization is used to solve the reconfiguration problem. The achieved eigenvalues (ai) and eigenvectors (vai) and the recovered gain (Kf) obtained through Particle Swarm Optimization for reconfigurable control is given by:
0.8568 0.2243i
0.8568 0.2243i
Elitism number = 2.


SIMULATION RESULTS AND PERFORMANCE EVALUATION.

ai/ PSO
0.9075
0.989
– 0.0046 0.0016i
– 0.0046 – 0.0016i
– 0.0036
0.0014
vai/PSO
0.9539
0.9539
0.7183
0.0011
0.112
Abnormal
Normal
PRODUCTCONCENTRATION(mol/l)
0.11
– 0.0012 – 0.0009i
0.1618 0.2527i
– 0.0012 0.0009i 0.1618 0.2527i
– 0.0035
0.6957
0.7951
0.6065
Set point
Recovered
0.108
K f ( LQR PSO )
0.112
PRODUCT CONCENTRATION (mol/l)
0.11
0.108
0.106
1.0e + 003 * 0.0277
1.4506
0.0013
– 0.0126
– 0.0028
0.0002
0.0001
0.0059
LQR+PIGA
SET POINT
0.106
0.104
0.102
0.1
0.098
0 10 20 30 40 50 60 70 80 90 100
TIME INSTANTS
LQR+PSO
0.104
0.102
0.1
0.098
0 20 40 60 80 100
TIME INSTANTS
Figure 3. Reactor concentration (mol/l) of the normal.
Figure 5. Reactor Concentration for 30% actuator failure introduced at 65th sampling instants.
439
PRODUCT TEMPERATURE (K)
438.5
438
437.5
439
Set point LQR+PSO
LQR+PIGA
Normal Abnormal
Set point
Recovered
PRODUCT TEMPERATURE(K)
438.5
438
437.5
437
436.5
437
436.5
436
0 10 20 30 40 50 60 70 80 90 100
TIME INSTANTS
436
0 10 20 30 40 50 60 70 80 90 100
TIME INSTANTS
Figure 6. Reactor Temperature for 30% actuator failure introduced at 65th sampling instants.
Figure 4. Reactor Temperature (K) of the normal.
ACHIEVED EIGEN VALUE AND EIGEN VECTOR 
GAIN 

LQR+PSO 
0.8580 0.2219i 0.8580 – 0.2219i ai / PSO 0.9159 0.989 0.0046 0.0016i 0.0046 – 0.0046 0.0017i – 0.9557 – 0.9557 vai/PSO – 0.0012 0.0009i – 0.0012 0.0009i 0.1631 0.2449i 0.1631 0.2449i 
– 0.0038 0.7443 – 0.0034 0.6678 
0.0008 0.0006 0.7792 0.6268 
K 77.7744 0.3809 17.6834 116.7231 – 0.7949 109.6617 
– 0.3070 0.0383 
LQR+PIGA 
0.8517 0.2377i 0.8517 – 0.2377i ai / PI GA 0.9755 0.9999 – 0.0048 0.0020i – 0.0048 0.0020i 0.9659 0.9659 vai/PIGA 0.0012 0.0006i 0.0012 0.0006i 0.1659 0.1989i 0.1659 0.1989i 
– 0.0044 0.8321 – 0.0029 0.5546 
0.004 0 1.00 0.6357 
K 52.43 0.2873 15.6774 – 0.3194 – 75.453 – 0.535 – 97.45 0.52373 
A.Kavitha, S.Sutha
23
International JoTurAnBaLlEOf3.EnAgCiHneIEerVinEgD REIeGseEaNrcShTRanUdCTTUecRhEnBoYloPgSyO(IJAENRDTP)I,GICASFEOMRN2O01R3MCALonSfYeSreTnEcMe Proceedings
In Figure 5 and 6, the response shows that under normal condition the process variables are able to track the set point changes for given step input. To demonstrate the performance of Reconfigurable state feedback controller optimized by the Particle Swarm Optimization technique, the fault is introduced at 65th sampling instant in terms of loss in the control effectiveness factor (30%) i.e. =0.7. The faulty system response (represented in Green line) shows the presence of undershoot and oscillation which is reduced in the Recovered system response (represented in Red line).
CONCLUSION
In this thesis the design of an optimized Eigenstructure Assignment based reconfigurable LQR control is designed for linear discrete Continuous Stirred Tank Reactor (CSTR) against partial actuator failure and is investigated. To perform the servo operation linear discrete state is augmented and ProportionalIntegral (PI) control action is incorporated in the design. LQR parameters (Qstate weighting matrix and Rinput weighting matrix) are optimized using Particle Swarm Optimization (PSO) and Parallel implementation of Genetic Algorithm (PIGA). The closed loop performance of the normal (Faultfree) system obtained using PSO and PIGA are compared. The best closedloop performance of the nominal system obtained using PSO technique is used to reconfigure the controller parameters against partial actuator failure in the process.
REFERENCES

Youmin Zhang, Jin Jiang: Integrated Active Fault Tolerant Control Using IMM Approach. In IEEE
Transactions on Aerospace and Electronic Systems, 2001.

Y.M .Zhang, J. Jiang: Active FaultTolerant control for partial actuator failures. In IEEE ProcessControl Theory Appl., Vol. 149, No. 1, 2002.

DOUIK Ali, LIOUANE Hend, MESSAOUD Hassani: Optimised Eigenstructure Assignment By Ant System And Lqr Approaches. In International Journal of Computer Sience and Applications, Vol. 5, No. 4, pp 45 – 56, 2008.

Meng Zhang, Weiguo Zhang, Yong Sun and Ning Ning : Reconfigurable Flight Control System Design using Eigenstructure Assignment and Particle Swarm Optimisation . In IEEE Transactions on Automatic Control, 2010.

Del Valle, Y.; Venayagamoorthy. G.K.; Mohagheghi, S.; Hernandez, J.C.; Harley, R.G.; "Particle Swarm Optimization: Basic Concepts, Variants and Applications in Power Systems," Evolutionary Computation, IEEE Transactions on, vol.12, no.2, pp.171195, April 2008.

GuoSheng Wang, Qiang Lv, Bing Liang and Guang Ren Duan. Design of Reconfiguring Control Systems via State Feedback Eigenstructure Assignment. In International Journal of Information Technology, Vol. 11, No. 7, 2005.

David.R ,Clark T: Parallel implementation of genetic algorithm. In Control Engineering practices ,Vol 3,No 1,2001.

Prakash J, Senthil R. Design of observer based nonlinear model predictive controller for a continuous stirred tank reactor. Journal of Process Control 2008;18:50414 .

B.W. Bequette, Nonlinear control of chemical processes: a review, Industrial Engineering and Chemistry Research 30 (1991) 13911413.

M. Siahi, M.A. Sadrina, A. Darabi: A new method for observer design using eigenstructure assignment and its application on fault detection and isolation.In World Applied Sciences Journal 6(1):2009.

S. Sutha and T. Thyagarajan: Eigenstructure assignment based multiobjective dynamic state feedback controller for MIMO system using NSGA
II. In Proceedings of Modelling, Identification and Control, IEEE Computer Society, July, 870875, 2010.

R. Eberhart, Y. Shi, and J. Kennedy, Swarm Intelligence. San Mateo, CA: Morgan Kaufmann, 2001.

C. Y. Huang and R. F. Stengel: Restructurable Control Using ProportionalIntegral Implicit Model Following. In the J. Guidance, Control, and Dynamics, 13(2):303309, 1990.
A.Kavitha, S.Sutha
24
International Journal Of Engineering Research and Technology(IJERT), ICSEM2013 Conference Proceedings