Modelling of Photovoltaic Water Pumping System Based on Artificial Intelligence

Modelling of Photovoltaic Water Pumping System Based on Artificial Intelligence

Mensia NawelTalbi Mourad Bouaїcha Mongi 

Photovoltaic Laboratory, Center of Researches and Technologies of Energy of Borj Cedria, Borj Cedria B.P N 95 2050 - Hammam Lif, Tunisia

Laboratory of Semiconductors, Nanostructures and Advanced Technology, Center of Researches and Technologies of Energy of Borj Cedria, Borj Cedria B.P N 95 2050 - Hammam Lif, Tunisia

Corresponding Author Email:
28 October 2018
27 January 2019
31 March 2019
| Citation



In this work, we develop and simulate an accurate model of a gird connected photovoltaic (PV) system composed by a PV generator, a buck converter and a motor-driven pump. The purpose of this study is to represent closely the system behavior under different climatic conditions (temperature and insolation), without using any order reduction technique and without identifying the operating point. The proposed modeling method is based on the multi-model technique for representing the PV system by a set of simple models. Those latters are derived from a polytopic transformation of premises that represent non-constant variables in PV system state equations. Every premise is bounded by two extreme values depending on temperature and insolation. Each model is associated to an activate function which estimates the model contribution in the global functioning of the system. The combination of all these models builds a flexible model which can accurately describe the behavior of the PV system. The simulation results prove that the obtained model can reproduce with high precision the conduct of the system state variables under different conditions (temperature and insolation). This model permits to develop an intelligent law of control to effectively attain Maximum Power Point (MPP) in a PV system.


convex polytopic transformation, vertex matrixes, extreme values, temperature, insolation, multi-model approach, photovoltaic system

1. Introduction

Traditional energy sources are not able to meet the growing demand for energy worldwide. Therefore, alternative energy sources like biomass, sunlight and wind, come into mentality. In this context, photovoltaic energy presents a source of attractive energy; it is non-polluting, sustainable and renewable. Today, photovoltaic energy is increasingly used in various terrestrial applications such as lightning, water pumping and telecommunication [1-2]. Though, its exploitation is less than conventional energy sources due to the high cost and the low efficiency of photovoltaic systems. The rational employment of solar energy imposes the MPP tracking of Photovoltaic Generators (PVGs) and estimation of the produced energy.

The exact identification of MPP is closely related to the efficiency of controller that provides it. Today, intelligent controllers are among the most powerful and robust controllers. Their synthesis is based on a model of system to command. Consequentially, it is necessary to have flexible and reliable model of the photovoltaic system, in order to perform an intelligent command for getting its performance as maximum as possible.

In literature, many techniques have been proposed for modeling and simulating (PVG) based on diode model. Although, a great number of this techniques are based on the adjustment of the curve I(V) at the points of the experimental curve, which requires the graphic extraction of curve slope at determined points or require an experimental analysis of photovoltaic device. These techniques have drawbacks, for the experimental curves I(V) are not permanently accessible in manufacturers’ data sheets, and devices for measurement of experimental data are not usually available as well. Others modeling techniques based on the modification of the electronic characteristic of a diode model of photovoltaic cell [3]. The model resulting from these techniques does not precisely reflect the real system behavior. Others research works proposed techniques of modeling tools based on the adjustment of the curve P(V). Although, this latter is based only on the characteristic I(V) and the peak power of the resulting model is not always peak power of the real device. Also, there are few modeling studies based on artificial intelligence such as multi-model approach [4]. These studies are confronted with the problem of determination of PV system operating points for constructing the model base.

This problem is overcoming in this work by extracting a basis of eight linear models from a convex polytopic transformation of scalar functions origin of nonlinearity in the system state equations. Linear basis models are calculated according to the extreme values of temperature and insolation, and the nonlinearity of the global system is rejected in the weighting functions. A flexible model which can accurately describe the behavior of the PV system is provided by combination of models assigned by the weighting functions. Consequently the PV system model is obtained without use of an order reduction technique which causes an information loss, and without requirement of the functioning point identification.

The rest of this paper is organized as follow: in section 2, is presented the photovoltaic pumping system to be study. Its model state, is given in section 3. The elaboration of the multi-model representation of the system is detailed in section 4. In section 5, are given the simulation results and interpretations. Finally, the conclusion and perspectives are given in section 6.

2. Study of the Photovoltaic Pumping System

A photovoltaic pumping system is generally composed by a photovoltaic generator, a motor-driven pump and a water tank to ensure the continuity of the supply when solar energy is not sufficient. Two types of systems can be used. In the first type, the Photovoltaic Generator (PVG) is directly coupled to the motor-driving pump. In the second type the PVG can be coupled either to an AC motor-driven pump through a DC-DC converter connected to a DC-AC converter [5], or to a DC motor-driven pump through a DC-DC converter. The latter type is the subject of our interest; indeed we will be modeling the photovoltaic system consists of a PV generator coupled to a DC permanent magnet motor through a buck convertor which equivalent circuit diagram is the following:


Figure 1. Electrical diagram of water pumping photovoltaic system

It’s consisted of the following components:

·    Photovoltaic Generator (PVG):

The PVG consist of Np arrays connected in parallel. Each array is composed of np panels associated in series, and each panel is formed of nc cells related in series. Therefore, the array is composed by Ns=np×nc cells. So, if the cell current is Ipvc and its voltage is Vpvc then the PVG current and voltage are successively: Ip=NpIpvc and Vp = NsVpvc.

·    Buck converter:

The chopper is controlled by the switched control variable u and its input is connected to the capacitor C, which allows transforming the current source Ip to a voltage source Vp.

·    The Diode:

The diode ensures continuity of the current in the motor at the moment of disconnecting power to avoid peak overvoltage.

·    The motor-driven pump:

we consider a DC permanent magnet motor with Constant flux, where the armature reaction and the switching phenomenon are assumed to be negligible.

3. State Model of the Photovoltaic Pumping System

The electrical behavior of a PV cell can be described with good accuracy by the single diode equivalent circuit [6]. The ideal equivalent circuit is a current source in parallel with a single-diode which governed by the following equation:

${{I}_{pvc}}\,=\,{{I}_{phc}}-{{I}_{sd}}\left( \exp \left( \frac{q{{V}_{pvc}}}{Kn{{T}_{C}}} \right)-1 \right)$      (1)


$I_{p v c}$: PV cell current,

$I_{p h c}$: Light current,

$I_{sd}$: Diode reverse saturation current,

$V_{pvc}$: PV cell voltage,

q=1.6. 10-19C, is telectron charge,

K=1.38∙10-23J/K, is the Boltzmann’s constant,

n is the ideal factor,

$T_c$: the actual cell temperature.

From the last equation, we can deduce the expression of the PVG output current $I_P$ according to its voltage $V_P$:

${{I}_{p}}\,=\,{{N}_{P}}{{I}_{phc}}-{{N}_{p}}{{I}_{sd}}\left( \exp \left( \frac{q{{V}_{p}}}{{{N}_{S}}Kn{{T}_{C}}} \right)-1 \right)$      (2)

If we consider: $I_{p h}=N_{p} I_{p h c}$ and $\cdot I_{S}=N_{p} I_{s d}$ then we can write the following:

${{I}_{p}}\,=\,{{N}_{P}}{{I}_{phc}}-{{N}_{p}}{{I}_{sd}}\left( \exp \left( \frac{q{{V}_{p}}}{{{N}_{S}}Kn{{T}_{C}}} \right)-1 \right)$    (3)

The motor-pump functioning is described by the following equations:

${{V}_{m}}\,=\,{{k}_{e}}\Omega +L\frac{d{{I}_{m}}}{dt}+R\cdot {{I}_{m}}$  (4)

$J\frac{d\Omega }{dt}=\,{{k}_{b}}{{I}_{m}}-f\cdot \Omega $  (5)


$k_e$ is the counter-electromotive force constant,

Ω is the angular velocity,

L is the electric inductance,

R is the electric resistance,

f is the motor viscous friction constant,

$k_b$ is the torque equation constant,

J is the total inertia moment.

The current $I_c$ and the voltage $V_m$ are expressed according the switched control variable u as following:

${{I}_{C}}\,=\,{{I}_{p}}-u\cdot {{I}_{m}}$ (6)

${{V}_{m}}\,=\,\,\,u\cdot {{V}_{p}}$  (7)

Combining the above equations gives:

$\left\{ \begin{matrix}   \frac{d{{V}_{P}}}{dt}=\frac{1}{C}{{I}_{P}}-\frac{1}{C}{{I}_{m}}\cdot u  \\   \frac{d{{I}_{m}}}{dt}=\frac{1}{L}V{}_{P}\,\cdot u-\frac{R}{L}{{I}_{m}}-\frac{{{k}_{e}}}{L}\Omega   \\   \frac{d\Omega }{dt}=\frac{k{}_{b}}{J}{{I}_{m}}-\frac{f}{J}\Omega   \\\end{matrix} \right.$  (8)

If we consider the state vector $\mathrm{x}=\left[\begin{array}{l}{\mathrm{x}_{1}} \\ {\mathrm{x}_{2}} \\ {\mathrm{x}_{3}}\end{array}\right]$ where: $\mathrm{x}_{1}=\mathrm{V}_{\mathrm{p}} ; \mathrm{x}_{2}=\mathrm{I}_{\mathrm{m}} \cdot \mathrm{and} \cdot \mathrm{x}_{3}=\Omega$. Then the model state of the water pumping photovoltaic system presented in Figure (1) can be written as following:

$\left\{ \begin{matrix}   {{{\dot{x}}}_{1}}=\frac{1}{C}{{I}_{P}}-\frac{1}{C}{{x}_{2}}\cdot u  \\   {{{\dot{x}}}_{2}}=\frac{1}{L}{{x}_{1}}\,\cdot u-\frac{R}{L}{{x}_{2}}-\frac{{{k}_{e}}}{L}{{x}_{3}}  \\   {{{\dot{x}}}_{3}}=\frac{k{}_{b}}{J}{{x}_{2}}-\frac{f}{J}{{x}_{3}}  \\\end{matrix} \right.$   (9)

When we suppose that $\alpha\left(x_{1}\right)=\frac{I_{p}}{C}$, then system of equations (9) can be written as follow:

$\left\{ \begin{matrix}   {{{\dot{x}}}_{1}}=\alpha \left( {{x}_{1}} \right)-\frac{1}{C}{{x}_{2}}\cdot u  \\   {{{\dot{x}}}_{2}}=\frac{1}{L}{{x}_{1}}\,\cdot u-\frac{R}{L}{{x}_{2}}-\frac{{{k}_{e}}}{L}{{x}_{3}}  \\   {{{\dot{x}}}_{3}}=\frac{k{}_{b}}{J}{{x}_{2}}-\frac{f}{J}{{x}_{3}}  \\\end{matrix} \right.$  (10)

where $α(x_{1})$ varies depending on the temperature and insolation.

The system non-linearity comes from the exponential function which describes the internal diode behavior in the expression of Ip. This function prevents to isolate Vp by a reliable linearization. Consequently, the dependency between Ip and Vp constitutes a major obstacle for the command synthesis. To solve this problem, we will use in this paper, a transformation procedure of system in order to writing the model state in multi-model form. That makes it possible to calculate easily a robust control for the system, under different climatic conditions of temperature and insolation.

4. Multi-Model Representation of the Photovoltaic Pumping System

4.1 Multi-model principle

A multi-model is a system representation composed by a set of simple models easy to manipulate, each of which is valid in a well-defined feature space corresponding to a part of global system’s behavior [7]. One degree of validity is associated to each of these models forming a base of model. An activation function which tends to 1 in the feature space and tends to zero outside, specifies the validity degree which estimates the model contribution in the global functioning of the system. The global system’s behavior can then be described by the combination of all models of the base. In literature, we can mention three approaches which are largely used for obtaining a multi-model representation: by identification, by linearization and by non-linear sectors transformation [8]. In this paper, we will be interested in the last approach. It is a convex polytopic transformation of scalar functions origin of nonlinearity.

4.2 Multi-model elaboration 

In order to represent a pumping system in a multi-model form, we will first transform the non-linear model state (10) into a quasi-Linear model with Variable Parameters (quasi-LPV) having the following form:

$\dot{x}=A(\rho(x, u)) \cdot x+B(\rho(x, u)) \cdot u$     (11)

It deserves mentioning that the parameters vector ρ(x,u) will be different in the matrixes A and B. However, for simplicity, a single notation will be used in those two matrixes.

Since the state variable $x_1$ can vary depending the temperature and insolation between a non-zero minimum and a maximum value, therefore we can reformulate the system of equation (10) as follow:

$\left\{ \begin{matrix}   {{{\dot{x}}}_{1}}=\alpha \left( {{x}_{1}} \right)\cdot \left( \frac{{{x}_{1}}}{{{x}_{2}}} \right)-\frac{1}{C}{{x}_{2}}\cdot u  \\   {{{\dot{x}}}_{2}}=\frac{1}{L}{{x}_{1}}\,\cdot u-\frac{R}{L}{{x}_{2}}-\frac{{{k}_{e}}}{L}{{x}_{3}}  \\   {{{\dot{x}}}_{3}}=\frac{k{}_{b}}{J}{{x}_{2}}-\frac{f}{J}{{x}_{3}}  \\\end{matrix} \right.$   (12)

This makes it possible to establish the following quasi-LPV form of the photovoltaic pumping system:

$\left[ \begin{matrix}   {{{\dot{x}}}_{1}}  \\   {{{\dot{x}}}_{2}}  \\   {{{\dot{x}}}_{3}}  \\\end{matrix} \right]=\left[ \begin{matrix}   {\alpha \left( {{x}_{1}} \right)}/{{{x}_{1}}}\; & 0 & 0  \\   0 & -{R}/{L}\; & -{{{k}_{e}}}/{L}\;  \\   0 & {{{k}_{b}}}/{J}\; & -{f}/{J}\;  \\\end{matrix} \right]\left[ \begin{matrix}   {{x}_{1}}  \\   {{x}_{2}}  \\   {{x}_{3}}  \\\end{matrix} \right]+\left[ \begin{matrix}   -{{{x}_{2}}}/{C}\;  \\   {{{x}_{1}}}/{L}\;  \\   0  \\\end{matrix} \right]\cdot u$    (13)

Thereafter, we define the set of non-constant variables in matrixes A and B as the set of premises $V_z$ [9]:

$V_{z}=\left\{z_{1}(x), z_{2}(x), z_{3}(x)\right\},$ wherez$_{1}(x)=\alpha\left(x_{1}\right) / x_{1}$

$z_{2}(x)=-x_{2} / c, z_{3}(x)=x_{1} / L$

Each of these premises can be bounded by two extreme values depending on temperature and insolation:

$\underline{z_{1}}<z_{1}(x)<\overline{z_{1}}, z_{2}<z_{2}(x)<\overline{z_{2}}$


Then, we apply to these bounded premises a polytopic transformation based on the following lemma [10].

Lemma: Any nonlinear function: g(x):IR→IR satisfying $\underline{g}<g(x)<\overline{g} ; \forall x$ can be written as:

$g\left( x \right)={{G}_{1}}\left( x \right)\,\underline{g}+{{G}_{2}}\left( x \right)\,\bar{g}$


${{G}_{1}}\left( x \right)\,=\frac{\overline{g}-g\left( x \right)}{\overline{g}-\underline{g}},\,\ {{G}_{2}}\left( x \right)=\frac{g\left( x \right)-\underline{g}}{\overline{g}-\underline{g}}$

and the functions $\mathrm{G}_{\mathrm{i}}(\mathrm{x}), \mathrm{i} \in\{1,2\}$ satisfy the convex sum property i.e. $\mathrm{G}_{1}(\mathrm{x})+\mathrm{G}_{2}(\mathrm{x})=1,$ and $\ 0<\mathrm{G}_{\mathrm{i}}(\mathrm{x})<1, \forall  \mathrm{x}$.

After premises transformation, we obtain:      

${{z}_{i}}\left( x \right)\,={{F}_{i1}}\left( {{z}_{i}}\left( x \right) \right)\,\,\overline{{{z}_{i}}}+{{F}_{i2}}\left( {{z}_{i}}\left( x \right) \right)\underline{\,\,{{z}_{i}}},\,\ i\in \left\{ 1,\,2,\,3 \right\}$       (14)


${{F}_{i1}}\left( {{z}_{i}}(x) \right)\,=\frac{{{z}_{i}}\left( x \right)-\underline{{{z}_{i}}}}{\overline{{{z}_{i}}}-\underline{{{z}_{i}}}}\,,\,\,{{F}_{i2}}\left( {{z}_{i}}\left( x \right) \right)=\frac{\overline{{{z}_{i}}}-{{z}_{i}}\left( x \right)}{\overline{{{z}_{i}}}-\underline{{{z}_{i}}}},\,\ i\in \left\{ 1,\,2,\,3 \right\}$    (15)

${{F}_{i1}}\left( {{z}_{i}}\left( x \right) \right)\,\,+{{F}_{i2}}\left( {{z}_{i}}\left( x \right) \right)=1$   (16)

We will denote by $F_{ij}$ instead of $\mathrm{F}_{\mathrm{ij}}\left(\mathrm{z}_{\mathrm{i}}(\mathrm{x})\right), \mathrm{j} \in\{1,2\}$.

Our aim consists in evaluating the matrixes A and B from the polytope vertex matrixes which defined from the partitions of the premises that configure in these matrixes.


$\mathrm{A}\left(\mathrm{z}_{1}(\mathrm{x})\right)=\left[\begin{array}{ccc}{\mathrm{z}_{1}(\mathrm{x})} & {0} & {0} \\ {0} & {0} & {0} \\ {0} & {0} & {0}\end{array}\right], \mathrm{A}_{0}=\left[\begin{array}{ccc}{0} & {0} & {0} \\ {0} & {-\mathrm{R} / \mathrm{L}} & {-\mathrm{k}_{\mathrm{e}} / \mathrm{L}} \\ {0} & {\mathrm{k}_{\mathrm{b}} / \mathrm{J}} & {-\mathrm{f} / \mathrm{J}}\end{array}\right]$

Then the matrix A(ρ(x,u)) can be expressed as follow:

$A\left( \rho \left( x,u \right) \right)=A\left( {{z}_{1}}\left( x \right) \right)+{{A}_{0}}$   (17)

When $z_1 (x)$ is replaced by its expression of equations (14) and (15), we can obtain the next equation: 

$A\left( {{z}_{1}}\left( x \right) \right)={{F}_{11}}\cdot \left[ \begin{matrix}   \overline{{{z}_{1}}} & 0 & 0  \\   0 & 0 & 0  \\   0 & 0 & 0  \\\end{matrix} \right]+{{F}_{12}}\cdot \left[ \begin{matrix}   \underline{{{z}_{1}}} & 0 & 0  \\   0 & 0 & 0  \\   0 & 0 & 0  \\\end{matrix} \right]$     (18)

The next step is to multiply the members of equation (18) by the unit term, $(F_{21}+F_{22})(F_{31}+F_{32})$, which gives the equation (19).

$\begin{align}  & A\left( {{z}_{1}}\left( x \right) \right)={{\mu }_{1}}{{A}_{1}}+{{\mu }_{2}}{{A}_{1}}+{{\mu }_{3}}{{A}_{1}}+{{\mu }_{4}}{{A}_{1}}+{{\mu }_{5}}{{A}_{2}} \\  & +{{\mu }_{6}}{{A}_{2}}+{{\mu }_{7}}{{A}_{2}}+{{\mu }_{8}}{{A}_{2}} \\ \end{align}$   (19)


$\mu_{1}=F_{21} F_{31} F_{11}, \mu_{2}=F_{21} F_{32} F_{11}, \mu_{3}=F_{22} F_{31} F_{11}, \mu_{4}=F_{22} F_{32} F_{11}$

$\mu_{5}=F_{21} F_{31} F_{12}, \mu_{6}=F_{21} F_{32} F_{12}, \mu_{7}=F_{22} F_{31} F_{12}, \mu_{8}=F_{22} F_{32} F_{12}$

$A_{1}=\left[\begin{array}{ccc}{\overline{z_{1}}} & {0} & {0} \\ {0} & {0} & {0} \\ {0} & {0} & {0}\end{array}\right]$ and $A_{2}=\left[\begin{array}{ccc}{\frac{Z_{1}}{}} & {0} & {0} \\ {0} & {0} & {0} \\ {0} & {0} & {0}\end{array}\right]$

Consequently, the matrix $A(z_{1} (x))$ is written as a function of two vertex matrixes $A_1$ and $A_2$.

We proceed in the same way for B in order to express it according to the premises $z_2$ and $z_3$. We have:

$B\left( \rho \left( x,u \right) \right)=B\left( {{z}_{2}}\left( x \right),{{z}_{3}}\left( x \right) \right)=\left[ \begin{matrix}   {{z}_{2}}\left( x \right)  \\   {{z}_{3}}\left( x \right)  \\   0  \\\end{matrix} \right]$  (20)

First, B(ρ(x,u))  is decomposed with respect to $z_{2} (x)$, then with respect to $z_{3} (x)$:

 $\begin{align}  & B\left( {{z}_{2}}\left( x \right),{{z}_{3}}\left( x \right) \right)={{F}_{21}}\cdot {{F}_{31}}\left[ \begin{matrix}   \overline{{{z}_{2}}}  \\   \overline{{{z}_{3}}}  \\   0  \\\end{matrix} \right]+{{F}_{22}}{{F}_{31}}\left[ \begin{matrix}   \underline{{{z}_{2}}}  \\   \overline{{{z}_{3}}}  \\   0  \\\end{matrix} \right] \\  & +{{F}_{21}}{{F}_{32}}\left[ \begin{matrix}   \overline{{{z}_{2}}}  \\   \underline{{{z}_{3}}}  \\   0  \\\end{matrix} \right]+{{F}_{22}}{{F}_{32}}\left[ \begin{matrix}   \underline{{{z}_{2}}}  \\   \underline{{{z}_{3}}}  \\   0  \\\end{matrix} \right] \\ \end{align}$     (21)

After that we multiply the members of the equation (21) by the unit quantity (F_{11}+F_{12}). So we find the subsequent expression of B(ρ(x,u)):

$\begin{align}  & B\left( \rho \left( x,u \right) \right)={{\mu }_{1}}{{B}_{1}}+{{\mu }_{2}}{{B}_{3}}+{{\mu }_{3}}{{B}_{2}}+{{\mu }_{4}}{{B}_{4}}+{{\mu }_{5}}{{B}_{1}} \\  & +{{\mu }_{6}}{{B}_{3}}+{{\mu }_{7}}{{B}_{2}}+{{\mu }_{8}}{{B}_{4}} \\ \end{align}$   (22)


$B_{1}=\left[\begin{array}{c}{\overline{z}_{2}} \\ {\overline{z}_{3}} \\ {0}\end{array}\right]$, ${{B}_{2}}=\left[ \begin{align}  & {{z}_{2}} \\  & {{\overline{\overline{z}}}_{3}} \\  & 0 \\ \end{align} \right]$, ${{B}_{3}}=\left[ \begin{align}  & {{\overline{z}}_{2}} \\  & {{z}_{3}} \\  & \overline{0} \\ \end{align} \right]$, $B_{4}=\left[\begin{array}{l}{z_{2}} \\ {\overline{z}_{3}} \\ {\overline{0}}\end{array}\right]$

Therefore, matrix B(ρ(x,u)) is then written as a function of four vertex matrixes $B_1$, $B_2$, $B_3$ and $B_4$.

Finally the state model of the photovoltaic pumping system can be reformulating according the vertex matrixes as following:

$\dot{x}=\left(A_{0}+v_{1} A_{1}+v_{2} A_{2}\right) x+\left(v_{3} B_{1}+v_{4} B_{2}+v_{5} B_{3}+v_{6} B_{4}\right) u$     (23)


$\dot{x}=\left(A_{0}+v_{1} A_{1}+v_{2} A_{2}\right) x+\left(v_{3} B_{1}+v_{4} B_{2}+v_{5} B_{3}+v_{6} B_{4}\right) u$

$v_{1}=\mu_{1}+\mu_{2}+\mu_{3}+\mu_{4}, v_{2}=\mu_{5}+\mu_{6}+\mu_{7}+\mu_{8}, v_{3}=\mu_{1}+\mu_{5}$

$v_{4}=\mu_{3}+\mu_{7}, v_{5}=\mu_{2}+\mu_{6}, v_{6}=\mu_{4}+\mu_{8}$

It deserves mentioning that each vertex of the polytope defines a linear subsystem, and that the nonlinearity of the global system is rejected in the weighting functions $ν_i$.

5. Simulation Results

In order to show the reliability of the developed model, we will consider a photovoltaic pumping system which has the following characteristics:

- PVG: PV module: SPM (P) 240-260W

$I_{SC}=8.3 A,  V_{mp}=31.58 V,  I_{mp}=7.6 A; V_{OC}=36.6 V, α_{0}=0.05%/°C , N_{S}=60$.

$β_{0}=-0.34%/°C$, NOCT=45∓2°C.

Motor-pump DC: It is characterized by a nominal functioning point:

u_{n}=24V and I_{n}=12A; ω_{n}=2000tr/mn.

Its parameters are: $\mathrm{R}=1.072 ; \mathrm{L}=0.05 \mathrm{H} ; \mathrm{f}=8814.10^{-7} \mathrm{N} . \mathrm{m} . \frac{\mathrm{s}}{\mathrm{rd}} ; \mathrm{k}_{\mathrm{e}}=0.5 ; \mathrm{k}_{\mathrm{b}}=0.1 ; \mathrm{J}=476.10^{-6} \mathrm{Kg} . \mathrm{m}^{2}$ and $\mathrm{C}=4.7 \mathrm{mF}$.

Concerning the experimental conditions, we choose the maximum and minimum values of temperature and irradiance as follow:

$\mathrm{G}_{\max }=1200 \mathrm{W} / \mathrm{m}^{2} ; \mathrm{G}_{\min }=400 \mathrm{W} / \mathrm{m}^{2} ; \mathrm{T}_{\max }=60^{\circ} \mathrm{C} ; \mathrm{T}_{\min }=0^{\circ} \mathrm{C}$

The extreme values of the PVG voltage $V_p$ and the current $I_p$ are calculated according to the temperature and irradiance by the following formulas [11-12]:

$V_{p}=V_{m p}\left[1+0.0539 \ Ln \left(\frac{G}{G_{0}}\right)\right]+\beta_{0}\left(T_{C}-T_{C 0}\right)$         (24)

$\mathrm{I}_{\mathrm{p}}=\mathrm{I}_{\mathrm{ph}}-\mathrm{I}_{\mathrm{S}}\left(\mathrm{e}^{\frac{\mathrm{qV}_{\mathrm{P}}}{\mathrm{N}_{\mathrm{S}} \mathrm{KT}_{\mathrm{c}} \mathrm{n}}}-1\right)$       (25)


$\mathrm{T}_{\mathrm{C}}=\mathrm{T}_{\mathrm{a}}+\left[\frac{\mathrm{NOCT}-20}{800}\right] . \mathrm{G}$


$T_C$: Temperature of the PV module;

$T_a$: is the ambient temperature;

$T_{C0}$: Temperature of the PV module at STC;

NOCT: Nominal Operating cell Temperature;

$V_{mp}$: The maximum PVG voltage under standard conditions;

$α_0$: Coefficient of current according the temperature;

$β_0$: Coefficient of voltage according the temperature;

$\mathrm{T}_{0}=298^{\circ} \mathrm{k} ; \mathrm{G}_{0}=1000 \mathrm{Watt} / \mathrm{m}^{2}$.

The ideality factor n depends on PV cell technology and can be chosen in [13].

$I_{ph}$ is expressed according temperature and insolation by the following equation [14]:

$\mathrm{I}_{\mathrm{ph}}=\frac{\mathrm{G}}{\mathrm{G}_{0}}\left(\mathrm{I}_{\mathrm{SC}}+\alpha_{0}\left(\mathrm{T}_{\mathrm{C}}-\mathrm{T}_{\mathrm{C} 0}\right)\right)$         (26)

and $I_S$ is given by the equation (27).

$\mathrm{I}_{\mathrm{S}}=\mathrm{I}_{\mathrm{S} 0}\left(^{\mathrm{T}_{\mathrm{C}}} / \mathrm{T}_{\mathrm{C} 0}\right)^{3} \mathrm{e}^{\frac{\mathrm{qEg}}{\mathrm{N}_{\mathrm{S}} \mathrm{nK}}\left(\frac{1}{\mathrm{T}_{\mathrm{C} 0}}-\frac{1}{\mathrm{T}_{\mathrm{C}}}\right)}$       (27)

where Eg is the bang-gap energy of the semiconductor used in the cell, and $I_{S0}$ is given by:

$\mathrm{I}_{\mathrm{S0}}=\frac{\mathrm{I}_{\mathrm{SC}}}{\mathrm{e}^{\frac{\mathrm{qV} \mathrm{OC}}{\mathrm{N}_{\mathrm{S}} \mathrm{n} \mathrm{K} \mathrm{T}_{\mathrm{CO}}}}-1}$      (28)

After calculating the weighting functions $γ_i$ and the vertex matrixes relating to the model developed according to the preceding parameters, this latter is simulated in Matlab environment. The result is given in the following figure, which illustrate the evolution of system’s PVG voltage and that of model PVG voltage for different temperature and irradiance pairs chosen arbitrarily, and for two values of command randomly selected.


Figure 2. Evolution of the system’s PVG voltage and that of the model’s PVG voltage, for two different values of command and different couples of temperature and irradiance


Figure 3. Evolution of the angular velocity Ω of the system and that of the model. for two different values of command and different couples of temperature and irradiance

We note that the error between voltage generate by the system’s PVG and that generate by the model’s PVG under diverse climatic conditions of temperature and illumination, and for different command decreases up to 10-9. This result leads to a model able to imitate with a great precision the PVG voltage of the system.

To test the model aptitude to replicate the behavior of all system state variables, we have simulated the angular velocity Ω of the system and that of the model. As well as the motor current $I_m$ of the system and that of the model. The simulation result is given respectively in Figures 3 and 4.


Figure 4. Evolution of motor current $I_m$ of the system and that of the model, for two different values of command and different couples of temperature and irradiance

We notice according to the simulation results given in Figures 3 and 4, that the modeling error reaches 10-12. These results demonstrate the model success to reproduce with high precision the conduct of the system state variables. That attests again the high quality of the elaborate model.

6. Conclusion

The performance of a Photovoltaic pumping system highly depends on the weather conditions, such as temperature and insolation, so does the performance of the PV generator degrade with increasing temperature and decreasing insolation. To provide water continuously throughout the year, a PV pumping system must be led by intelligent control which depends heavily on the quality of the system model. The model must guarantee a correct description of the expected variation in the system. This objective is reached by the model structure proposed in this paper. Indeed, we have developed an accurate model that represents the behavior of the system under different climatic conditions of temperature and insolation. The multi-model approach is used to represent the system by a set of simple models, which calculates according the extreme values of temperature and insolation. The modeling begins with the establishment of the quasi-LPV form, then the identification of the premise variables, next the premise decomposition and finally the development of matrixes and weighting functions. The simulation results of the method proved the good quality of the established model. Therefore we recommend for future researches to control the system with a multi-model command based on the elaborated model. Indeed, each model of the base can allow the generation of a partial controller which is used to derive the global control to be applied to the system. The simplicity of the eight models can involves a powerful and easy control law which can be implemented in real time.

In perspective, this work will be continued with a practical implementation of the model by an Arduino board for increasing the efficiency of the system to extract the maximum energy.


We would like to thank all the people who contributed in some way to this work which was supported by the Center of Researches and Technologies of Energy of Borj Cedria, Tunis and Ministry of Higher Education and Scientific Research


[1] Pa, J., Gaurag, S.H. (2013). Modeling and Simulation of Solar Photovoltaic Module using Matlab/Simulink. International Journal of Research in Engineering and Technology, 02(03).

[2] Abdulkadir, M., Samosir, A.S.A., Yatim, H.M. (2013). Modeling and simulation of a solar photovoltaic system, its dynamics and transient characteristics in LABVIEW. International Journal of Power Electronics and Drive System (IJPEDS), 3(2): 185-192.

[3] Hongmei, T., Mancilla-Davida, F., Kevin, E., Eduard, M., Peter, J. (2012). A cell- to – module –to-a cell detailed model for photovoltaic panels. Solar Energy, 86(9): 2695-2706.

[4] Ravinder, K.K., Shimib, S.L., Chatterjib, S., Fahim, A. (2014). Modelling of solar PV module and maximum power point tracking using Anfis. Renew. Sustain. Energy Rev., 33: 602-612.

[5] Ait-Cheikh, S.M. (2007). Etude, Investigation et conception d’algorithmes de commande appliqués aux systèmes photovoltaïque. PhD thesis. Ecole Nationale Polytechnique.

[6] Petibon, S. (2009). Nouvelles architectures distribuées de gestion et de conversion de l’énergie pour les applications photovoltaïques. Ph.D. thesis, Toulouse university, France.

[7] Smith, R.M.T.A. (1997). Johansen Multiple Multimodel Approachs to Modelling and Control. Edition, Taylor and Francis.

[8] Nesrine, E. (2010). Approche neuronale de la representation et de la commande multimodèle de processus complexe. Ph.D. Thesis. Ecole Doctorale Science pour l’Ingénieur Université Lille, France. 

[9] Anca, M.N. (2010). Analyse et synthèse de multimodèle pour le diagnostic. Application à une station d'épuration. Ph.D. Thesis. National Polytechnic Institute of Lorraine, France.

[10] Kazuo, T., Hua, O.W. (2001). Fuzzy Control Systems Design and Analysis: A Linear Matrix Inequality Approach. A Wiley-Interscience Publication.

[11] Belhaj, M., Benouaz, T., Cheknane, A., Bekkouche, S.M.A. (2010). Estimation de la puissance maximale produite par un générateur photovoltaique. Revue des Energies Renouvelables, 13(2): 256-257.

[12] Habbati, B., Ramdani, Y., Moulay, F. (2014). A detailed modeling of photovoltaic module using MATLAB. NRIAG Journal of Astronomy and Geophysics, 3: 53-61.

[13] Bogdan, S.B., Salameh, Z.M. (1996). Methodology for optimally sizing the combination of a battery bank and PV array in a wind/PV hybrid system. IEEE Transactions on Energy Conversion, 11(2): 367-375.

[14] Tsai, H.L., Tu, C.S., Su, Y.J. (2008). Development of generalized photovoltaic model using MATLAB/Simulink. Proceedings of the World Congress on Engineering and Computer Science. San Francisco, USA.