 Open Access
 Total Downloads : 456
 Authors : Arathi George, Riya Mary Francis
 Paper ID : IJERTV4IS060683
 Volume & Issue : Volume 04, Issue 06 (June 2015)
 DOI : http://dx.doi.org/10.17577/IJERTV4IS060683
 Published (First Online): 20062015
 ISSN (Online) : 22780181
 Publisher Name : IJERT
 License: This work is licensed under a Creative Commons Attribution 4.0 International License
Model Reference Adaptive Control of Binary Distillation Column Composition using MIT Adaptive Mechanism
Arathi George
Department of Electrical and Electronics engineering, TKM College of Engineering,Kollam, India
Prof. Riya Mary Francis Department of Electrical and Electronics engineering, TKM College of Engineering,Kollam, India
AbstractThis paper presents a model reference adaptive control for binary distillation column based on the adaptation mechanism calledMIT rule.The goals are: to present a theoretical calculation procedure of a binary distillation column for simulation. Second, the controller design: a reducedorder linear model is derived suchthat it best reflects the dynamics of the distillation process and used as the reference model.Third, selecting an adaptation mechanism suitable for the system.Simulation results using MATLAB shows that the responses are satisfactory and stability is achieved.
Index TermsModel Reference Adaptive Control, Adaptation mechanism,MIT Rule, binary distillation column.

INTRODUCTION
Distillation is the most common method of separation of two or more components from a mixture using the relative volatility of the components.It is widely used in chemical industry to separate components with desired purity.Distillation columnsare made up of several components, each of which is usedeither to transfer heat energy or enhance mass transfer.
The distillation column model is computed isbased on a real petroleum project gas processingplant.The feed to the column is a pseudomixture of ligas and naphthas and the distillate product is gasoline and bottoms product is LPG.The input variables are feed flow rate F and feed compositionxF. The manipulated variables are reflux flow rate L and vapor boilup rate V. The quality of the outputproducts is the purityof the distillate, xD, higher/equal than98% and the impurity of the bottoms, xB, less/equal than 2%.The basic feed stock data and its actual compositions are basedon the reference [1].
The distillation column has 14 trays and each tray is modeled using mass balance equations of LV control configuration. The control of composition of distillate and bottoms product is not easy as the feed flow cannot be constant. Thus proper controlling techniques should be used.
Fig.1: Distillation equipment
The mathematical modelling of distillation column is obtained from Vu Trieu Minh and John Pumwa [2], which is based on physical laws such as mass balance equation. The nonlinear set of equations are linearized and order is reduced for create the reference model for MRAC.The controller design and control theories for adaptive controller proposed by Astrom was found very helpful [4].

MATHEMATICAL MODELLING OF BINARY DISTILLATION COLUMN
The purity of distillate or bottom is affected by twofundamental variables: feed split and fractionation. The feed split variable refers to the fraction of the feed that is takenoverhead or out the bottom. The fractionation variable refersto the energy that is put into the column to accomplish theseparation.Feed split variables are used as those have much stronger effect on product compositions than fractionation variables.
The LV control structure, which is called energy balance structure, used as standard thecontrol structure for dual composition control of distillation [5] has been adopted in this paper. In this control structure,he reflux flow rate L and the boilup flow rate V are used to control the primary outputsassociated with the product specifications.
The nonlinear system is modelled with the mass balance and component balance equations.The model is simplified under assumptions in [7].Constant relative volatility throughout the column and thevapourliquid equilibrium relation can be expressed as:
= (1)
1+( 1)
wherexn : liquid concentrations on nth stage; yn : vapour concentration on nth stage; : relative volatility.
The dynamic model can be expressed by the following equations [8]:
Condenser (n=N+2):
= + 1 (2) Tray n (n=f+2 to N+1):
= + (1 ) + (+1 ) (3) Tray above the feed flow (n=f+1):
= (1 ) + (+1 ) ( )
(4)
Tray below the feed flow (n=f):
= (1 ) + (+1 ) ( )
(5)
Tray n (n=2 to f1):
= (1 ) + (+)(+1 ) (6) Reboiler (n=1):
1 = + 2 1 1 (7)
The dynamic model is represented bynonlinear differential equations where x1=xB: liquidconcentration in bottom; x2: liquid concentration in the 1st tray, x3: liquid concentration in the 2nd tray; ;x15: liquidconcentration in the 14th tray; and x16=xD: liquidconcentration in the distillate. The distillation processsimulation is shown in Figure 3.
If there is no disturbance in the operating conditions, thesystem reaches steady state such that the purity of thedistillate product xDis 0.9654 and the impurity of the bottomsproduct xBis 0.0375 (Fig 2).
Table. I indicates the steady state values of concentration ofxnand ynon each tray.
Fig.2: Steady values of concentration of xn on each tray
Table I: Steady values of concentration of xn on each tray
The changes in feed flow rates and feedcompositionsof the feed stream can be considered disturbances.Simulations with disturbances indicate that the qualityof the output products deteriorate with more disturbances (Table II).
Table II: Product quality depending on the change in feed flow rate
The designed system does not achieve the operationalobjective of the product quality (xD 0.98and xB 0.02) andthe product quality will get worse dealing with disturbances.Hence we will use an adaptive controller – MRAC to take thesystem from these steady state outputs to the desired outputtargets (xD 98%andxB 2% ).

LINEARISED DISTILLATION COLUMN MODEL
The nonlinear set of differential equations is linearized using Taylors series. In order to obtain a linear mathematical model for a nonlinear system, it is assumed that thevariables deviate only slightly from some operating condition [6].The higherorder terms in theseries can be neglected when the variation is small. The linearization leads to 16 order linear model in state space.
= + ()
= , where;
The model reference adaptive control is composed of two loops. The inner loop consists of an ordinary feedback
1 1
2 2
control loop. The outer loop consists of the adaptation mechanism and also looks like a feedback loop. The model
z(t) =
16
. ,
.
16
output plays the role of the setpoint while the process output is the actual measurement [9].
The mechanism for adjusting the controller parameters in a model reference adaptive system can be obtained in two
=
=
u(t) = ,
1 1 =
ways: by using gradient method or by applying stability theory [4].The desired behaviour of the system is specified
by the reference model,and the parameters of the controller
16
y(t) =
16
= (8)
are adjusted on the error,which is the difference between the output of the closed loop system and the reference
The full order model can be reduced to second order as
follows:
= 1 (0) ; (9)
model
We found that the reducedorder linear model in(11) can also maintain the similar steady state outputs as thebasic
1+
nonlinear model. Now we use this model as reference
where: G(0) is the steady state gain: G(0) =
1, isthe time constant:
model MRAC totake the process plant from these steady state outputs to thedesired targets amid the disturbances
= + (1 ) + (1 )
(10)
and the plantmodelmismatches as the influence of the feed
ln
stock disturbances.
The MIT Rule is adopted as the adaptation mechanism is
where(kmole): the total holdup of liquid inside the column; (kmole): liquid holdup in the condenser;
(kmole): liquidholdup in the reboiler; : the impurity sum; and S : theseparation factor.
As the result of calculation, the reducedorder linearmodel of the plant is a first order system with a time constant
of =1.9588 is: [3]
this paper which is also called as the gradient method. According to MIT Rule there will be one controller parameter that can be adjusted. The parameter can be adjusted in such a way that the loss function
= 1 2 (12)
2
is minimized where e is the error between the process
output and the model output ym. To minimize J,the parameters have to be changed in the direction of negative
= 1
0.0042 0.0060
(11)
gradient of J, that is,
= = (13)
1+1.9588
0.0050 0.0072

MODEL REFERNCE ADAPTIVE CONTROL
SYSTEM
The two main reasons for applying adaptive control widely in petroleum industriesare:First, most of chemical processes arenonlinear and the linearized models are used to design thecontrollers, so that the controller must change and adapt to themodelplant mismatch; Second, most of the chemical processes arenonstationary or their characteristics are changed with time,this leads again to adapt the changing control parameters.
Fig 3: MRAC block diagram
The feedforward control law using MIT Rule is:
= (14) and the parameter update rule is:
= (15)
The feedback control law using MIT Rule is:
= 1 2(16)
and the parameter update rule is:
1 = (17)
2 = (18)

SIMULATION RESULTS AND DISCUSSIONS
The reducedorder linear model in (11) can also maintain thesimilar steady state outputs as the basic nonlinear model. This model is used as anMRAC totake the process plant from these steady state outputs (xD =0.9654 andxB= 0.0375) to the desired targets (0.98 xD 1and 0 xB 0.02) amid the disturbancesand the plantmodel mismatches as the influence of the feed stock disturbances.
Fig 4: Distillate purity of MRAC FF System
Fig 5: Bottoms impurity of MRAC FF System
Simulation in Fig 4, 5 shows that the controlled outputs of MRAC system with feedforward controller using MIT Rule xD and xB are stable and tracking to the model outputs and the reference setpoints amid thedisturbances and the plant model mismatches.
While the Fig 6, 7 shows that the outputs xD and xB are stable and tracking to the model outputs and the reference setpoints are more accurately controlled by using feedback controller. The same can be inferenced from table III.
Table III: Comparision of MRAC Syatem, Non Linear System, Non Linear System with disturbance
Purity of distillate Product(mole fraction)
Impurity of bottoms Product(mole fraction)
NonLinear Model
.9654
.0375
NonLinear Model with disturbance
.9225
.1087
Feedforward MIT
.9708
.03612
Feedback MIT
.9719
.03703

CONCLUSION
In this paper, an MRAC controller is designed using the MIT Rule for both feedforward gain and feedback gain even under disturbances. Then both feedback and feedforward controller performance are compared and found feedback controller is more efficient. The empirical evidence shows that the system is stable for small values of adaptation gain in the case of MIT rule but high gains lead to instability.

FUTURE SCOPE
The main drawback of MIT feedforward and feedback controller is both the distillate product and bottoms product composition are not controlled to desired value simultaneously. By designing MRAC PID the performance can be improved so that the desired composition is achieved. Also decoupling can incorporated so that distillate and bottoms product purity can be controlled separately and more accurately.
Fig 6: Distillate purity of MRAC FB System
Fig 7: Bottoms impurity of MRAC FB System
REFERENCES

PetroVietnam Gas Company,Condensate Processing Plant Project Process Description Document No. 8203602BM 01.PetroVietnam,pp.154, 1999

Vu Trieu Minh and John Pumwa ,Modeling and Adaptive Control Simulation for a Distillation Column,14th International Conference on Modelling and Simulation, 2012

Skogestad, S., and Morari, M. The Dominant Time Constant for Distillation Columns. Computers & Chemical Engineering., pp. 607617, 1987.

K,J.Astrom and Bjorn Wittenmark, Adaptive ControlPearson Education, NJ, USA, 1995.

Skogestad, S., and Morari, M. Control Configuration Selection for Distillation Column. AIChE., Vol 33, No. 10, pp. 16201653, 1987.

Katsuhiko Ogata, Model Control Engineering. PrenticeHall International , Englewood Cliffs, NJ, USA,1982.

Papadouratis, A., Doherty, M. F., and Douglas, J. M. Approximate Dynamic Models for Chemical Process Systems. Ind. Eng. Chem. Res.,pp. 546522 ,1989.

George, S. (1986). Chemical Process Control: An Introduction to Theory and Practice, PrenticeHall, ISBN: 9780131286290, New Jersey.

Stephanopoulos, G.Chemical Process Control. Prentice Hall International, 1984.