The Kronecker SummationMethod for Robust Stabilization Applied to a Chemical Reactor

The paper focuses on robust stabilization where the suitable parameters of a simple continuous-time PI controller are determined through a combination of the Kronecker summation method, sixteen plant theorem, and an algebraic approach to control design in the ring of proper and stable rational functions. The initial theoretical background is followed by an illustrative experiment which includes computation of the controller and verification of control results for a continuous stirred tank reactor with exothermic reaction modelled as a fourth-order interval system.


Introduction
A popular and practically preferred approach to appropriate control synthesis for plants with complex properties consists of construction of an uncertain model, which should cover all possible operating points, parameter variations and nonlinear behavior, and consequent robust control design [1,2].Although the resulting controllers can have a simple structure and fixed parameters, which is even acknowledged from the practical application viewpoint, both process modelling and control design are generally nontrivial tasks.
Chemical reactors, which belong among the most interesting and critical processes in all chemical engineering, represent the class of systems suitable for robust control applications.Their control is usually affected by very complex behavior and, moreover, bounded with potential safety problems.A common type of reactor is known as a continuous stirred tank reactor (CSTR) [3].The mathematical model of CSTR, robustly stabilized in this paper, has been constructed in [4].Moreover, the same work has presented stabilization of the CSTR using technique from [5,6] embellished with a polynomial control.Besides, robust static output feedback control has been utilized to this CSTR in [7].The idea of robust stabilization applied in this paper is similar to [4], but the contribution is mainly in the use of alternative methods, that is, combination of the Kronecker summation method [8] and an algebraic approach to control design under the ring of proper and stable rational functions (R PS ) allowing the elegant tuning [9][10][11][12][13][14] which has been already investigated in [15,16].
This paper deals with design of robustly stabilizing continuous-time PI controllers for a continuous stirred tank reactor (CSTR) in which exothermic reaction occurs.The controlled plant is assumed as a fourth-order interval plant, and easily tunable PI controller is designed in order to robustly stabilize the closed control loop.As a synthesis method, the combination of Kronecker sum method, sixteen plant theorem and an algebraic approach, is utilized.This compound and its application to the chemical reactor model represent the key contribution of the work.The robust stabilization is verified on a simulative example.
The paper is organized as follows.In Section 2, the Kronecker summation method is described.Section 3 then extends the idea for an interval plant.Subsequently, Section 4 briefly outlines the algebraic approach to controller design itself.Next, a chemical reactor description and specific control experiment are provided in the extensive Section 5.And finally, Section 6 offers some conclusion remarks.

C(s) G(s)
Figure 1: The classical closed control loop.

Kronecker Summation Method
An interesting technique for computation of stabilizing PI controllers based on the Kronecker summation has been presented in [8].The main purpose of the method is to find possible variations of PI controller parameters which ensure stability of the classical one-degree-of-freedom closed control loop according to Figure 1, where is transfer function of a controlled system with fixed coefficients, and is a PI controller.Remember that Kronecker summation of two general square matrices Q (with size q-by-q) and R (r-by-r) is defined as [17] where I q , I r are identity matrices of size q-by-q and r-by-r, respectively, and where ⊗ stands for the Kronecker product.For example, The significant feature of the final square matrix (3) (qrby-qr) is that it has qr eigenvalues which are pair-wise combinatoric summations of the q eigenvalues of Q and r eigenvalues of R. In other words, the operation of Kronecker summation induces the "eigenvalue addition" feature to the matrices.One can exploit this property to obtain the equation for which all pairs (k P , k I ) leading to purely imaginary roots comply.The characteristic equation of the closed control loop (Figure 1) has the form (5) Now define and transform (5) into matrix differential equation where M is matrix of size n-by-n and x n ] T .Equations ( 5) and ( 7) are connected by relation Evidently, the same complex variable s is both the root of ( 5) and the eigenvalue of M. Thanks to the fact that M is a constant matrix, the complex conjugates of s must satisfy also (9), that is, det(s * I − M) = 0. (10) Consequently, as it was presented in [8], if s = jω is the root of (5) it must be the eigenvalue of M. Furthermore, s * = − jω is also the root of ( 5) and the eigenvalue of M. As the summation of two eigenvalues s = jω and s * = − jω is equal to zero, the Kronecker summation of two matrices must be singular when such correspondence of k P , k I , and ω occurs.Thus, determines the stability boundary in (k P , k I ) plane, because every couple (k P , k I ) satisfying (11) means that the same couple inserted into (5) will lead to the pair of conjugate purely imaginary roots or zero roots.These positions together with the line k I = 0 are the only ones where the system stability can shift.Generally, the stability boundary splits the (k P , k I ) plane into the stable and unstable areas.The selection of the stabilizing region(s) can be performed through a test point and corresponding representative polynomial within each area.

Robust Stabilization of Interval Systems Using PI Controller
So far, we could apply the Kronecker sum method to calculate a region of stabilizing PI controller parameters for plant with fixed coefficients.However, the papers [5,6,8] extended this (or an alternative) stabilization technique also for interval systems.The simple idea consists of the combination with the so-called sixteen plant theorem [1,18,19].This proposition says that a first-order controller (such as PI controller) robustly stabilizes an interval plant if and only if it stabilizes its sixteen Kharitonov plants.
The values b − i , b + i , a − i , and a + i represents, respectively, lower and upper bounds for parameters in numerator and denominator.
Remember that the Kharitonov plants are defined as where i, j ∈ {1, 2, 3, 4}; B 1 (s) to B 4 (s) and A 1 (s) to A 4 (s) are the Kharitonov polynomials for the numerator and denominator of the interval system (12), that is [20], and analogically Thus, robust stabilization of interval system directly follows from the simultaneous stabilization of all sixteen fixed Kharitonov plants.The final stability region for the original interval system is given by the intersection of sixteen partial regions obtained individually via the Kronecker summation method.

Algebraic Tuning of PI Controller
Recapitulate that for now we are able to compute all possible robustly stabilizing variations of proportional and integral parts in PI compensator.However, the final choice of the controller from the obtained stability region represents another task.A very good survey on PI(D) control issues is provided, for example, in [21].For the purpose of this paper, a simple but effective solution offers an algebraic approach to control synthesis [9][10][11].This fractional design is grounded in general solutions of Diophantine equations in R PS , Youla-Kučera parameterization of controllers, and conditions of divisibility in R PS .One of main advantages of this method can be seen in the existence of the only tuning parameter m > 0 which serves for influencing the control behavior.This paper is not going to explain the details of this approach.Interested readers can find them, for example, in [11][12][13][14].This work only takes advantage of one simple tuning rule.The parameters of PI controller (2) can be computed in compliance with where a 0 and b 0 come from the first-order nominal controlled plant and where the tuning parameter m > 0 can be selected, for example, using the recommendation from [12] Suitable coefficient k depends on the size of first overshoot.For example, the choice k = 2.14, which is applied in the following simulation experiments, leads to 3% overshoot.

CSTR Description.
The controlled process adopted from [4,7,22] represents hydrolysis of propylene oxide to propylene glycol in a CSTR.More specifically, the chemical reaction of the process is Besides propylene oxide and water, methanol is also added to the CSTR in order to improve the solubility of propylene oxide in water.The excess of water ensures higher selectivity to propylene glycol and eliminates consecutive reactions of propylene oxide with nascent propylene glycol.Dependence of the rate constant of chemical reaction on the temperature can be described by the well-known Arrhenius equation where k means reaction rate constant, k ∞ is the preexponential factor, E represents the activation energy, R signifies the universal gas constant, and T r is the temperature of the reaction mixture.
Under assumption of ideal mixing in the CSTR, constant reacting volume, and the identical volumetric flow rates of the inlet and outlet streams, the mass balance of the system can be given by where V r stands for the reacting volume, c i means the molar concentration of the ith component, c i0 is the feed molar concentration of the ith component, q r represents the volumetric flow rate of the reaction mixture, ν i determines the stoichiometric coefficient of the ith component, and r = kc C3H6O is the molar rate of the chemical reaction.Further, independency of the specific heat capacities, densities, and volumetric flow rates on temperature or mixture composition has been supposed.Moreover, the mixing volume and the heat of mixing have been neglected.So, the simplified enthalpy balance of the reaction mixture and the simplified enthalpy balance of the cooling medium introduced in the monograph [23] and subsequently in the papers [4,7] can be formulated as, respectively, where T means the temperature, ρ represents the density, c p is the specific heat capacity, Δ r H 0 stands for the reaction enthalpy, U is the overall heat transfer coefficient, and A is the heat exchange area.Furthermore, meaning of subscripts is as follows: 0 is for the feed, c for the cooling medium, and r for the reaction mixture.Equation ( 22) represents a standard in CSTR design.Interested readers can find the specific values of all constant parameters and steady-state inputs of the CSTR in tables in [4,7].
In fact, the physical parameters (such as reaction enthalpy, pre-exponential factor, and overall heat transfer coefficient) of this CSTR are not known exactly but they are supposed to vary within some intervals (see again [4,7]).However, this paper particularly takes advantage of the final mathematical model of the CSTR introduced in [4], where it is obtained in the linearized form of an interval transfer function: with parameters which can vary within the following bounds: (24)

Control Experiments.
The interval transfer function (23) with parameters (24) describing the CSTR is considered to be the controlled plant.The first of its sixteen Kharitonov plants ( 13) is G 1,1 (s) = −0.0245s 2 − 0.0199s − 0.000574 s 4 + 0.903s 3 + 0.2299s 2 + 0.0062s + 0.0001094 . ( Corresponding closed-loop characteristic (5) can be computed as which means that the matrix (8) takes the following form: where The stability boundary is determined by (11).The positions of such pairs (k P , k I ) which fulfill (11) are plotted in Figure 2. The decision on area of stability and instability can be simply done with the assistance of an arbitrary testing  point from the appropriate set.Moreover, the figure is supplemented with the half line k I = 0 bordering the "upper" part of the stabilizing area.Further, we must repeat the analogical procedure for all sixteen Kharitonov plants.The stability areas for the Kharitonov plants are depicted in Figure 3, and its zoomed version with highlighted intersection of all particular areas of stability is shown in Figure 4.The highlighted area represents the final region of robustly stabilizing PI controller parameters for the original interval transfer function of CSTR (23), (24).
The following question is how to find the practically convenient PI controller from the obtained robust stability region.This paper utilizes the algebraic-based approach outlined in the Section 5. Nevertheless, this method requires a first-order fixed nominal model (17) of controlled plant in order to compute the final controller of appropriate (first) order, that is, with PI structure.So, the paper employs very simple but effective method of obtaining such model.In the first step, the fourth-order model with fixed parameters has been computed using the average values of interval parameters in ( 23), (24): G A (s) = −0.0268s 2 − 0.0163s − 0.00046445 0.74155s 3 + 0.16505s 2 + 0.0102s + 0.0001753 . (29) Next, the simplest possible approximation takes advantage of neglecting the higher than zero-order powers of s in numerator and higher than first-order powers of s in denominator, that is, the corresponding coefficients are supposed to be zero.This approximation is easy but efficient enough as can be seen in the following steps.Obviously, it results in a first-order nominal model suitable for the applied synthesis: Further, (16) has been used for the calculation of PI controller parameters while the tuning parameter m has been chosen according to (18) in order to obtain 3% first overshoot for the nominal case.Consequently, the resulting controller is given by This controller is located inside the stability region as can be seen in Figure 5.It means that the regulator robustly stabilizes the CSTR ( 23), (24).Besides, the controller lies on the curve hypothetically connecting the other potential controllers tuned by various parameters m > 0.
Finally, the robust stability and control performance is confirmed and demonstrated in Figure 6.It shows the output signals of the control loop with designed PI controller (32) and 128 "representative" systems from the interval family describing the CSTR ( 23), (24).The minimum and maximum values of each interval parameter have been used.It results in 2 7 = 128 systems for simulation.On top of that, the red curve represents the output for the nominal system (30).Furthermore, Figure 7 depicts corresponding 128 + 1 control (actuating) signals.As far as control quality is considered, the results from Figure 6 might not be impressive at first sight, but there is a tradeoff between the simplicity of applied control algorithm and the performance here.As can be seen, the only one off-line tuned feedback PI controller with fixed coefficients has been utilized for controlling the chemical reactor with all the possible variations of uncertain parameters.We can observe that the CSTR is robustly stabilized successfully.

Conclusion
The main aim of the paper has been to present a possible approach to robust stabilization of a CSTR with exothermic reaction modelled as a fourth-order interval system.The developed easy but effective method of PI controller design has combined the Kronecker summation method with sixteen plant theorem and the algebraic tools.The proposed technique is applicable to a wide range of real processes, provided that they can be expressed by means of interval system and subsequently temporarily approximated by a first order model.The demerit of the method can be seen in a missing guarantee of coincident nominal performance and robust stability before the design process itself.They have to be verified during or after the design.Nevertheless, the applicability has been clearly demonstrated on the example where the CSTR has been successfully robustly stabilized.