Issue 
Int. J. Simul. Multisci. Des. Optim.
Volume 6, 2015



Article Number  A3  
Number of page(s)  6  
DOI  https://doi.org/10.1051/smdo/2015003  
Published online  01 July 2015 
Research Article
Nonlinear continuoustime generalized predictive control of solar power plant
Process Control Laboratory, Control Engineering Department, École Nationale Polytechnique, 10 avenue Pasteur, HassanBadi, 16131
Algiers, Algeria
^{*} email: khoukhi.billal@gmail.com
Received:
3
April
2015
Accepted:
30
May
2015
This paper presents an application of nonlinear continuoustime generalized predictive control (GPC) to the distributed collector field of a solar power plant. The major characteristic of a solar power plant is that the primary energy source, solar radiation, cannot be manipulated. Solar radiation varies throughout the day, causing changes in plant dynamics and strong perturbations in the process. A brief description of the solar power plant and its simulator is given. After that, basic concepts of predictive control and continuoustime generalized predictive control are introduced. A new control strategy, named nonlinear continuoustime generalized predictive control (NCGPC), is then derived to control the process. The simulation results show that the NCGPC gives a greater flexibility to achieve performance goals and better perturbation rejection than classical control.
Key words: Predictive control / Solar power plant / Nonlinear system description
© B. Khoukhi et al., Published by EDP Sciences, 2015
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
In recent years, the use of renewable energies for reducing environmental impact of fossil energies has grown considerably. There is a pressing need to find inexhaustible energy sources. Solar energy is by far the most abundant energy source of renewable energy. One of the main characteristics of a solar power plant is that the primary energy source (solar radiation) cannot be manipulated. Besides, the solar radiation intensity depends on daily and seasonal cycle variations, like clouds, atmospheric humidity, and air transparency. This justifies the relevance of solar power plant control. The distributed collector field consists of a series of parabolic mirrors that reflect solar radiation on a pipe where oil gets heated while circulating. The objective of the control system in a distributed collector field is to maintain the outlet oil temperature at a desired level in spite of disturbances, such as changes in the solar irradiance level (caused by clouds), mirrors' reflectance, or inlet oil temperature. This is accomplished by varying the flow of the fluid through the field. The field exhibits a variable delay time that depends on the control variable (flow). The transfer function of the process varies with factors such as irradiance level, mirrors’ reflectance, and oil inlet temperature. Those changes affect the dynamics of the process. In order to reduce the impact of disturbances, a feedforward controller is needed [1].
Control of solar plants has been addressed in many ways. In Igreja et al. [2], a nonlinear control using exact feedback linearization and Lyapunov’s approach is presented. In Camacho et al. [3], a survey on advanced control strategies tested in the Acurex experimental plant at the PSA (Solar Plant of Almería in Spain) is described. In Limon et al. [4], a robust tracking model predictive control (MPC) is presented. In GarciaGabin et al. [5], a sliding model predictive control is developed. In Gallego & Camacho [6], an adaptive statespace model predictive control is presented.
This paper applies nonlinear continuoustime model predictive control into solar power plants. The main benefit of the proposed technique is to achieve optimal performance as defined by a performance index without the need of online optimization that is normally required for implementing a model predictive control strategy. The online nonlinear optimization imposes a heavy computational burden, which requires extensive computing power and a long sampling time in engineering implementation, this is why it is natural to develop a generalized predictive control (GPC) in a continuoustime setting by pursuing an analytic approach, where it is intended to develop a closed form nonlinear model predictive control (NMPC) where online optimization is not required.
The paper is organized as follows: Section 2 describes the mathematical models of a collector loop. In Section 3, the new NCGPC algorithm and the obtaining of the control law are developed. In Section 4, simulations and results are described. Finally, the paper draws to a close with some concluding remarks.
2 Solar power plant
2.1 Process description
The main objective of the solar plant is to use solar radiation to heat oil passing through the field. The field consists of collectors that are arranged in rows that form parallel loops. The inlet oil is extracted from the bottom of a storage tank and is passed through the field using a pump located in the field inlet. This fluid is heated and then introduced back into the storage tank to be used for electricity generation.
The considered solar power plant is located in Almería, Spain, which is a tested and validated model of a trough collector loop [1]. As shown in Figure 1, the Acurex field consists of 480 collectors arranged in 20 rows that form 10 parallel loops and lie along an eastwest axis [3].
Figure 1. Schematic diagram of the Acurex distributed collector field. 
2.2 Mathematical model
In the case of control applications, the most important consideration is to describe the outlet oil temperature as a function of the flow rate of the oil. Generally, such a model could be formed by applying the energy conservation law for a length control volume of dl in the collector over a time interval dt. The actual equation can be set up as:(1)where f the subscript refers to the oil fluid. Also t is the time (s); x is the position (m); ρ is the oil density (kg/m^{3}); c is the specific heat of oil (J/kgK); A is the crosssection of the pipeline (m^{2}); T is the oil temperature (°C); I is the solar radiation; no is the geometric efficiency; H _{ l } is the global coefficient of thermal losses (W/m^{2} °C); G is the collector aperture (m); T _{ a } is the environment temperature C; q is the oil flow rate (m^{3}/s); K _{opt} is the optical efficiency.
Parameters ρ _{ f } and c _{ f } are function of temperature and so function of time and position.
Notice that the partial differential equation represented in equation (1) when discretized in space the following expression is obtained:(2)where n is the nth segment in which the pipe is divided. The system has four inputs (the ambient temperature T _{ a }, the effective solar radiation InoK _{opt}, the inlet temperature T(0, t), and the oil flow q) and one output (the outlet temperature T(N, t)) where N is the last segment that the pipe is divided. Figure 2 shows the schematic of discretization of the parabolic trough crosssection.
Figure 2. Parabolic trough crosssection. 
Let us put T(n, t) = x _{ n }, so that the output of the plant is then Y = x _{ N } = T(N, t) and the inputs are: U _{ q }= q, U _{ Ta } = T _{ a }, U _{irr } = IGnoK _{opt}, and U _{Tin} = x _{0 } = T(0, t). We can then write:
The control objective is to maintain the outlet temperature Y = x _{ N } = T(N, t) at constant value using the oil flow U _{ q } = q. As explained previously, the plant has many control difficulties: (1) the main source of energy (solar irradiation) cannot be manipulated. (2) The solar irradiation is a fast perturbation that acts directly to the output. (3) The system is highly nonlinear mainly at high temperatures of the oil. This is why a classical control law cannot deal with these difficulties, and an advanced control is necessary to encounter those difficulties.
3 Nonlinear generalized predictive control
Model predictive control (MPC) is a control strategy which uses a mathematical model of the plant to predict its evolution and, taking into account this evolution, find the control signal which minimizes a cost function [8].
Various nonlinear predictive control methods for discretetime systems have been developed. The main shortcoming of those methods is that online dynamic optimization is required, which, in general, is nonconvex, and its computational burden grows exponentially with the decision variables. Hence, computational issue is an active subject in predictive control. There are several model predictive controllers implemented in the Acurex field. For example, in Silva et al. [9] a variable sample time statespace adaptive controller is proposed. In Igreja et al. [10], a statespace receding horizon controller, which explicitly explores the distributed parameter character of the solar field, is developed.
Solar power plant is described by continuoustime differential equations and it is natural to develop GPC in a continuoustime setting. The interest in developing a continuoustime predictive control in nonlinear environment is that the online optimization involved in NMPC is nonconvex, which implies that the optimization procedure might be terminated at a local minimum. Consequently, poor performance even instability might result after the control sequence yielded by the online optimizer is implemented [7]. The heavy computing power needed in the online optimization requires to develop GPC in a continuoustime setting by pursuing an analytic approach.
The NCGPC algorithm solves those constraints and has been applied several times for other different systems as in Chen [12]. Alike, this paper uses the advantages of this control law to achieve optimal performance of solar power plants. However, to apply the NCGPC strategy, the following assumptions should be respected [11, 12]:

the dynamic of zeros exists and is asymptotically stable,

all stats are measurable,

the relative degree ρ is defined,

the output y(t) and the reference signal are continuously differentiable in time [13].
Once the assumptions are verified we could start the development of the output. In the NCGPC algorithm the prediction of the output is given by its development of Taylor series:(3)where t represents the present time, t + τ represents the prediction instant, ρ represents the relative degree of the plant, and R(τ ^{ ρ }) is the high order element when neglected (3) becomes:(4)
In the same way as (4) the prediction of the reference is obtained:(5)
The objective is to find a control law that makes the tracking error tend to zero when t tends to ∞. The prediction error is defined as:(6)
The cost function is then given as:(7)
Theorem: Consider the solar plant model (2) and the performance criterion J (7). Then, the control law U _{ q } given in equation (8) minimizes J. (8)
Proof:
The output in a solar power plant is the outlet temperature that represents the last segment that the pipe is divided; its derivative is obtained by subtracting n by N in equation (2):(9)
The relative degree of the plant of the plant is ρ = 1. By subtracting ρ by its values in equations (4) and (5):
The tracking error is defined as the difference between the output and the reference signal:
Let us define vectors Λ(τ), Y(t), Ω(t) so that:
Let us define another vector E(t) so that:
When replacing in equation (7):
The formula to express the cost function can be obtained.
When replacing in the criterion by its expression given in equation (10):
The control law can be obtained by minimizing the cost function (7):
The output of the NCGPC controller can be found by resolving the equation (11):(12)
When reference signal x(t) is canceled, the control law is the inputoutput linearization (FL) without poles placement terms. For SISO nonlinear systems, when the relative degree is almost equal to 4, the closed loop stability is ensured [12]. For ω(t) equal to zero, the control law (12) is exactly that of the FL [16]. The only difference is that the coefficients K that characterize the stability of the closed loop system are different from those of the FL. For NCGPC they result directly from the minimization of the cost function, whereas in FL, they are arbitrarily chosen to ensure the closed loop stability. The NCGPC is then, an optimal FL in a new coordinate system [17].
4 Simulation results
4.1 NCGPC response
Due to the degree of nonlinearities, linear methods cannot perform well in this plant under all operation conditions, as can be seen in Figure 3, the fixed parameters PID with feedforward control developed in Camacho et al. [14] cannot deal with the high nonlinearities at high temperature.
Figure 3. NCGPC response compared with PID + feed forward. 
The NCGPC developed deals well with these nonlinearities and possesses a lower time response, this is due to the predictive and the nonlinear characteristic of the control law.
4.2 Influence of prediction horizon
To study the influence of the prediction horizon let us take the general form of firstorder system similar to our system:
The resulting system is given by:
By replacing u by the control law given in equation (12) we have:
We can then obtain the characteristic polynom of the closed loop system:
Let us consider D _{1}(p) the denominator of the firstorder transfer function H _{1}(p): G _{1} and θ represent the static gain and the time constant, respectively. Then D _{1}(p) = 1 + θp.
The transfer function H _{1}(p) has a unique pole p = −1/θ. On the other hand the polynomial P _{1T } (λ) has unique root λ = −K _{10}. By equalizing those two equations we obtain:
But K _{10} = 3/2T = 1/θ. And so:
We can conclude that the response time is function of the prediction horizon:
When T increases the root λ = −3/2T tends to zero but remains in the negative side of the complex plan, and the response time increases linearly. The stability is then guaranteed and the time response is set by manipulating the prediction horizon. To illustrate this result we have simulated the influence of the prediction horizon in Figure 4.
Figure 4. Influence of the prediction horizon in NCGPC. 
As demonstrated, by increasing the prediction horizon, the response time increases. We can notice also, that the control signal is more energetic when the control horizon is small; we should then determine a compromise between the response time and constraints.
4.3 Robustness tests
To ensure that the control law gives proper results, robustness tests are conducted. The two main different disturbances affecting the robustness of the solar power plant system are the change in system parameters and the solar irradiation.
The change in system parameters is due to the synthetic oil used in the system. The characteristic of the synthetic oil used in the solar power plant can change during the entire heating process. Also, parameters ρ _{ f } and C _{ f } can change with time and so influence the response of the entire plant. It is necessary that the control law developed deals with those variations. Figure 5 shows the influence of a variation of 5% of those parameters.
Figure 5. Response of NCGPC with parameter variations. 
As can be seen, the variation introduced affects the response of the control law by acting directly on the outlet temperature, this can be seen by a deviation from the reference signal, this deviation is too small to be taken into account, and the response can be considered as satisfying.
As explained, the second test deals with outside perturbation. The solar irradiation is the main perturbation; it can directly affect the outlet temperature and acts with high frequency in an unpredictable way. In order to simulate a fast variation of this perturbation, we have used a model [15] of solar irradiation to simulate the day of 04/05/2013 at Ouargla (Algeria) city. The results of this simulation are given in Figure 6.
Figure 6. Response of NCGPC to the irradiation of the 04/05/2013 at Ouergla (Algeria). 
As can be seen, the very strong variation of the solar irradiation of that day directly affects the outlet temperature, the NCGPC reacts correctly to this variation by adjusting the oil flow, and this makes the tracking process work well with small deviation in the presence of the strong variation of the irradiation. This response makes the NCGPC a robust controller that deals correctly with perturbations and parameter variations.
5 Conclusion
In this paper a nonlinear continuoustime generalized predictive controller has been developed. The simulations, which have been carried out, demonstrate that the performances of the controller are very good in most cases. The NCGPC controller surpasses the performance of finetuned PID + feed forward. The control law shows its optimality due to its cost function minimization and its robustness due to the use of the real output for predicting the future output. The optimality is, however, approximate due to the approximating system dynamic with its Taylor expansion. The continuous characteristic of the control law developed allowed resolving the computational issue and the nonconvex optimization. In general, the outcome can be considered satisfactory.
References
 Camacho EF, Berenguel M, Rubio FR. 1997. Advanced control of solar plant. SpringerVerlag: London, UK. [CrossRef] [Google Scholar]
 Igreja JM, Lemos JM, Barao M, Silver RN. 2003. Adaptive nonlinear control of a distributed collector solar field, in Proc. of European Control Conference 2003, ECC03. [Google Scholar]
 Camacho EF, Rubio F, Berenguel M, Valenzuela L. 2007. A survey on control schemes for distributed solar collector fields. Part I: modeling and basic control approaches. Solar Energy, 81, 1240–1251. [CrossRef] [Google Scholar]
 Limon D, Alvarado I, Alamo T, Ruiz M, Camacho EF. 2008. Robust control of the distributed solar collector field ACUREX using MPC for tracking, in Proceedings of 17th IFAC World Congress, Seoul, Korea, Vol. 17, p. 958–963. [Google Scholar]
 GarciaGabin W, Zambrano D, Camacho EF. 2009. Sliding mode predictive control of a solar air conditioning plant. Control Engineering Practice, 17, 652–663. [CrossRef] [Google Scholar]
 Gallego AJ, Camacho EF. 2012. Adaptative statespace model predictive control of a parabolictrough field. Control Engineering Practice, 20, 904–911. [CrossRef] [Google Scholar]
 Chen W, Ballance DJ, Gawthrop PJ. 1999. Analytic approach to generalized predictive control of nonlinear systems. IEE: Savoy Place, London WCPR OBI, UK. [Google Scholar]
 Camacho EF, Bordons C. 2004. Model predictive control, 2nd edition. Springer Verlag: London. DOI: 10.1007/9780857293985. [Google Scholar]
 Silva RN, Lemos JM, Rato LM. 2003. Variable sampling adaptive control of a distributed collector solar field. IEEE Transactions on Control Systems Technology, 5(11), 765–771. [CrossRef] [Google Scholar]
 Igreja J, Lemos J, Silva J. 2005. Adaptative receding horizon control of a distributed solar collector field, in Proceedings of the 44th IEEE conference on decision and control, and the European control conference. p. 1282–1287. DOI: 10.1109/CDC.2005.1582335. [Google Scholar]
 Chen WH. 2001. Analytic predictive controllers for nonlinear systems with illdefined relative degree. IEE Proceedings – Control Theory and Applications, 148(1), 9–16, ISSN: 13502379. [CrossRef] [Google Scholar]
 Chen WH. 2003. Optimal control of nonlinear systems: a predictive control approach. Automatica, 39(1), 633–641. [CrossRef] [MathSciNet] [Google Scholar]
 Nijmeijer H, van der Schaft AJ. 1990. Nonlinear dynamical control systems. Springer Verlag: New York. [CrossRef] [Google Scholar]
 Camacho EF, Berenguel M, Rubio FR, Martínez D. 2012. Control of solar energy systems, in Advances in industrial control. SpringerVerlag London Limited. DOI: 10.1007/9780857299161_2. [CrossRef] [Google Scholar]
 ASHRAE. 1999. ASHRAE handbook: HVAC applications. Atlanta (GA). [Google Scholar]
 Isidori A. 1995. Nonlinear control systems, 3rd edition. SpringerVerlag, Englewood Cliffs: New York. [CrossRef] [Google Scholar]
 Dabo M, Chafouk H. 2010. Commande prédictive non linéaire a temps continu de systèmes complexes. Editions universitaires européennes. ISBN13: 9786131543395. [Google Scholar]
Cite this article as: Khoukhi B, Tadjine M & Boucherit MS: Nonlinear continuoustime generalized predictive control of solar power plant. Int. J. Simul. Multisci. Des. Optim., 2015, 6, A3.
All Figures
Figure 1. Schematic diagram of the Acurex distributed collector field. 

In the text 
Figure 2. Parabolic trough crosssection. 

In the text 
Figure 3. NCGPC response compared with PID + feed forward. 

In the text 
Figure 4. Influence of the prediction horizon in NCGPC. 

In the text 
Figure 5. Response of NCGPC with parameter variations. 

In the text 
Figure 6. Response of NCGPC to the irradiation of the 04/05/2013 at Ouergla (Algeria). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.