Modelling and Simulation of A Photovoltaic System
Modelling and Simulation of A Photovoltaic System
Modelling and Simulation of A Photovoltaic System
CERTIFICATE
This is to certify that the thesis entitled “MODELLING AND SIMULATION OF A
PHOTOVOLTAIC SYSTEM WITH MPPT CONTROLLER” submitted by
RAJEEV RANJAN (110EE0411) in partial fulfilment of the requirements for the award
of Bachelor of Technology degree in Electrical Engineering at National Institute of
Technology, Rourkela is an authentic work carried out by them under my supervision
and guidance. To the best of my knowledge the content of this thesis has not been
submitted to any other university/Institute for the award of any Degree.
NIT Rourkela
i
ACKNOWLEDGEMENTS
We would also like to express our deepest gratitude to Prof. J K SATPATHI, Head
of the Electrical Engineering Department, National Institute of Technology,
Rourkela for his invaluable advice and for providing us the platform to successfully
complete our project.
We would like to thank all others for their technical as well as moral support and
their significant contribution in successful completion of this project.
Lastly, we would like to appreciate our family and friends, for their advice and
support during the course of this project.
ii
DEDICATED TO
Our beloved family, friends & the
almighty
iii
ABSTRACT
Solar energy has become one of the most popular renewable energies due to the drastic increase
of global energy demand and rapid depletion of conventional fossil fuel resources. Advanced
power electronic systems are required to utilize and develop renewable energy sources. In solar
system, utilizing of maximum power from the source is one of the most important functions of the
power electronic converters. The output power characteristics of a PV array are influenced by the
environmental factors so that the conversion efficiency is low. Therefore a maximum power
tracking (MPPT) technique is necessary to track the peak power to maximize the efficiency.
This project focuses on mathematical piece wise linear modelling and simulation of three diode
model of a PV solar cell. The shunt and series resistance is calculated by using Newton Raphson
iterative method by varying solar irradiance and temperature. The MPPT algorithm is implemented
by using buck converter in MATLAB/SIMULINK.
iv
CONTENTS
ACKNOWLGEMENT ii
ABSTRACT iv
TABLE OF CONTENETS v
LIST OF FIGURES vi
ABBREVIATIONS AND ACRONYMS viii
CHAPTER – 1
INTRODUCTION
1.1 INTRODUCTION 1
1.2 MOTIVATION 1
1.3 OBJECTIVE OF PROJECT WORK 2
1.4 LITERATURE SURVEY 2
1.5 ORGANIZATION OF THESIS 3
CHAPTER – 2
PV MODELLING
v
CHAPTER – 3
MAXIMUM POWER POINT TRACKING (MPPT)
CHAPTER –5
5.1. CONCLUSION 24
FUTURE WORK
REFERENCE
vi
LIST OF FIGURES
LIST OF TABLES
vii
ABBREVIATIONS AND ACRONYMS
PV Photovoltaic
DC Direct Current
viii
CHAPTER – 1
1.1. INTRODUCTION
Over the past decade, solar energy source has shown the potential by robust and continuous
growth. In universal, P-V or I-V curve have distinct point called Maximum power point, Where
With maximum efficiency the total P-V system is been operated and gives its maximum power
output. For this, Maximum Power Point Tracker (MPPT) algorithm is used to extract the maximum
available power at every operating point. Various MPPT algorithms are available in literature such
as Perturb and observe (P & O), Hill climb search (HCS), Advanced P & O, Fuzzy logic, neural
network algorithms etc. MPPT techniques are implemented by using different types of power
electronics converter like buck converters, boost converters, SEPIC converter, and buck converter
etc. Due to the PV system being quite expensive, optimization of the system performance are
widely used by predictive performance tools. Under the Optimized Test Conditions (OTCs), PV
manufacturers are providing limited tabular data which points to a cell of 25°C at 1000 W/m2 and
at 1.5 air mass spectral distributions. But PV cells work under standards very far from the OTCs.
Research needs failed due to this reason, when the data in the datasheet shown. Usually, A physical
image which express PV current with the environmental parameters and most basic technical
features depending on this there are analytical equations, such as the operating voltage, the
surrounding temperature, and the irradiance. By mathematical formulae Analytical ways are
common methods in estimating the parameters. Furthermore, analytical methods use the curve
characteristics or semiconductor variables which are absent in the datasheet. Recently, PV
parameter prediction and analysis is considered as a multidimensional parametric problem. Several
solving ways, such as Newton Raphson iterative technique, Genetic Algorithms (GA), Chaos
Particle Swarm Optimization (CPSO), Firefly, and Pattern Search (PS) were issued. These
algorithms bring out the appropriate parameters by reducing the Root Mean Square Error (RMSE)
as the fundamental relation in the optimization system.
1.2. MOTIVATION
It is clear that the world energy demands cannot be make my conventional source .Taking into
account the growing population it is very important to find as sustainable source of power. In this
regard renewable energy presents unique potential. Solar power is an important of the renewable
energy mix in the world today but several factor have limited its widespread use and adoption.
1
One of the major issue is the relatively costlier installation of PV array .That has rendered it
uncompetitive is the energy market. So need of the ours is more research into making solar cheaper
modelling and simulation of PV arrays is done to estimate its characteristics batter and extract
maximum power possible .The I-V and P-V characteristics of a PV array for varying atmospheric
condition enable as to make better design interfacing of a synchronous buck converter with PV
array acts as impedance matching device and make it possible to get optimum power from it
.
1.3 OBJECTIVE
Development of a piece wise linear model for PV module consisting three diodes connected
in parallel with voltage controlled resistance for high performance.
Estimation of the series resistance and shunt resistance by Newton Raphson iteration
technique for matching of power and I-V characteristics for the variation of temperature and
irradiation.
Development of MPPT controller by using buck converter.
1.4.LITERATURE SURVEY
The researches on “Solar Energy Systems” is going on since a long time all around the world by
many researchers and engineers. To improve the efficiency and control the output of solar cell
within economic cost, advanced control mechanism are required as it comes with low efficiency.
This modern advanced control mechanism based on power electronics is known as “Maximum
Power Point Tracking (MPPT)”.
A PV cell operation is analysed at various solar radiation intensity and temperature over a day [1].
The current vs voltage and power vs voltage characteristics is obtained by mathematical modelling
and MATLAB simulation of PV system [2].
The cause of variation in the temperature and solar radiation in PV cell’s characteristics curve and
effect of partial shading over PV cell are assessed by an accurate analytical modelling of a two
diode PV cell [3].
Various PV cells are simulated and compared with the manufacturing datasheet for error analysis
[4]. By using MATLAB/Simulink, generalised modelling of PV cell is obtained and partial shading
effect over different MPPT techniques is used to obtain optimum power of system and algorithm
2
is studied for optimised system [5]. For partial shading variation of peak points in power vs voltage
characteristics is shown [6].
For different temperature from STP is discussed for mathematical and theoretical analysis of a PV
cell [7]. At various solar radiation intensity, a relative analysis of different maximum power point
tracking (MPPT) is done to make the output power nearly constant in each cases
3
CHAPTER – 2
Modelling of PV cell
4
Rs I
Id Ish
IL
V
Rsh
AkT
qv
I d I 0 e 1 2.1.2
Where
5
2.3. R sh Model
Similar method to find load current and diode current as in classical model but determined
by the Newton Raphson method.
And RS 0 (2.3.3)
Where
𝑉𝑚𝑎𝑥 - Voltage at maximum power
6
The advantages of R sh and R s model over classical model is matching of R sh and R s for maximum
Power output by iterative method. Hence for high accurate modelling, using datasheet information
those parameters can be calculated which is open circuit voltage, maximum power and short circuit
current. Newton-Raphson method is used for computing these parameters and computed, by
equation (2.3.1). Newton-Raphson method has the good and a quick convergence with an initial
values near the root, as given in equation (2.3.2). So, within a few iteration a better solution of R sh
function of R s by equation (2.3.1) is obtained until the optimum experimental power and voltage
values (𝑃𝑚𝑎𝑥 ,𝑉𝑚𝑎𝑥 ) matching with datasheet values of different module.
Rs
I
D1 D2 D3
V
IL R1 R2 R3 Rsh
V1 V2 V3
To control the power (constant power) voltage controlled resistor in series with diodes is used near
the maximum power point voltage (𝑉𝑚𝑎𝑥 ) although there is little change in voltage near the 𝑉𝑚𝑎𝑥
(between.9Vmax to 1.1 Vmax with 0.1 Vmax interval) by suitable switching ON and OFF of the
diodes connected in parallel rather than switching of all the diodes at the same time.
7
The detailed analysis of switching ON and OFF of the diodes connected in parallel and
determination of the value of the voltage controlled register are given below. The switching
voltages are taken at 𝑉 1 = 0.9Vmax, V2 = Vmax, V3 =1.1Vmax.
(V1 = 0.9 Vmax). When the generated voltage from the PV panel is less than 𝑉 1 , then all the diodes are
in OFF state, and there will be no flow of current through those diodes, So all the light generated
current from the PV panel will flow through the load and a very small portion of the IL will pass through
the shunt resistance Rsh, that’s why the current is almost remains constant in this interval
In this case diode D1 will be ON and the value of PV current fails from IL to IL-Id1
VMAX 0.9VMAX
R1 (2.4.2)
I D1
1.1VMAX 0.9VMAX
I D1 (2.4.3)
R1
The current passing through the diode D2 becomes:
I D 2 I D I D1 I (2.4.4)
𝑅2 Computed as follows
8
In this case all the three diodes will be turned ON and the current passing through the
load will be zero at Voc. The current passing through the diode D1 is given in the
equation (2.4.6)
VOC 1.1VMAX
I D1 (2.4.6)
R1
Similarly the current passing through the diode D2 is given as:
VOC 1.1VMAX
I D2 (2.4.7)
R2
I D3 I L I D1 I D 2 I (2.4.8)
In this proposed model Rsh and Rs are calculated by Newton – Raphson iteration method same to
the procedure proposed given in equation [2.3-1, 2, 3].
The Classical model, Rsh model, proposed piece wise linear model are designed in
MATLAB/Simulink for different model and their comparative study is given below.
Current vs Voltage (I-V) and Power vs Voltage (P-V) characteristics are shown in Fig 2.4 and
Fig 2.5 for various models at temperature 250c and 500c and the variation is observed.
9
Fig. 2.5 Power-Voltage characteristics of PV module at different temperature
From the on top of I-V and P-V characteristics of various models (classical, Rsh and planned
model) at completely different temperature, the planned model has its characteristics a great
deal near the plot obtained from the producing datasheet
I-V and P-V characteristics area unit shown in Fig two.6 and Fig two.7 for numerous models at
irradiation 800W/m2 and 1000W/m2 and also the variation is discovered. The legends utilized
in higher than figure is additionally
Applicable here
35
0
Fig. 2.6 Current vs Voltage Characteristics of PV module at different irradiation
10
VOLTAGE (V)
Fig. 2.7 Power vs Voltage Characteristics of PV module at different irradiation
From the higher than i-v and p-v characteristics of various models (classical, Rsh and planned
model) at totally different irradiation, the planned model has its characteristics a great deal near
the plot obtained from the producing datasheet.
At low temperature the classical model gives desired values but significant error is
introduced as temperature and insolation increases.
The shunt resistance model gives higher results compared to classical model as
approximation is finished by Newton Raphson methodology but it's conjointly not
satisfactory for higher temperature and insolation amount.
The three diode model gives economical results with significant accuracy overcoming the
restrictions of 1𝑠𝑡 and 2𝑛𝑑 models.
With increase in insolation the current, voltage and corresponding power output will
increase.
As temperature will increase current will increase to a small degree however voltage
decreases considerably leading to decrease in power output.
11
2.5.3 Relative error in power w.r.t. temperature
6
CLASSICAL MODEL R sh
5
MODEL PROPOSED PIECEWISE LINEAR
RELATIVE ERROE INPmp (%)
4 MODEL
0
25 50 75
TEMPRATURE (*C)
6
CLASSICAL MODEL
RELATIVE ERROR IN Voc
5 Rsh - MODEL
1
(%)
0
25 50
75
TEMPRATURE (*C)
12
CHAPTER 3
3.1 Introduction to Maximum Power Point Tracking
The maximum power is extracted for the optimal utilisation of panel as installation cost of PV
panel is high. Over the period of a day, temperature and solar isolation changes. So, main moto
is to obtain voltage that corresponds to maximum power output for the different radiation
intensity and temperature. Solar panel without optimised technique leads to wastages of solar
energy. The variation in solar radiation limits the optimum utilisation of solar panel which can
be overcome by Maximum Power Point Tracking (MPPT) technique to obtain optimum power
from the panel at every instant of time. The aim of MPPT technique is to get operating voltage
of panel corresponding to maximum power at the load. MPPT is an electronic circuit with an
ability to track voltage (VMPP) or current (IMPP) periodically to get optimal power from PV
panel at all environmental conditions.
Some required characteristics of MPPT techniques are:
• Cheaper
• Easy to implement and to install.
• Rapid tracking response for the dynamic analysis and tracking.
• Ability to track the MPP for varying solar radiation and temperature over wide range.
v. Current sweep
13
3.3 Perturb and Observe Algorithm
14
2
Table 3.1 Result of Different Models at 25°C and 1000 W/m
15
CHAPTER 4
DESIGN OF BUCK CONVERTER
The basic operation of the buck converter with an inductance is controlled by two switches
(Usually a semiconductor device and a diode). With controlled switching of converter, output
power can be controlled for desired value. It is assumed that input and output voltage remains
constant over the cycle.
When switch is turned off, inductor act as voltage source which supplies power to the load and
diode turns on. . The voltage across the output will be the summation of input voltage and the
voltage across the inductor. Thus stepped up voltage is obtained across the load. Capacitor
maintains the constant voltage across the load. This operation is called OFF-Mode operation.
In this way output power control is achieved.
16
.4.1. MPPT controller
To obtain the optimistic power by interlinked PV or solar powered system from transmission
line which is indeed optimised system, MPPT methods are used in grids linked with inverters,
battery chargers powered by solar power system and other same type of device. At no-load,
relationship between PV voltage and other parameter like radiation, resistance and temperature
is too much complex. Hence dependent efficiency on P-V and I-V characteristics is non- linear.
To resolve this problem an appropriate load is connected at output side so that desirable power
output can be obtained at any condition with changing solar radiation and temperature. The
MPPT device is incorporated with the converter to convert electric power and supplies
converted voltage or current, regulates and filters for operation of loads such as motors,
batteries or power grids.
Out of many MPPT methods and techniques used in PV system one method is discussed
here:
The main concept of the algorithm is that PV voltage at Maximum Power Point is directly
proportional to the open circuit voltage Voc. The constant of proportionality depends on the
manufacturing of the cell like solar cell technology, fill factor and climate conditions.
𝐾1 = 𝑉𝑚𝑝𝑝 *𝑉𝑜𝑐
The constant, K1 varies in between 0.71 and 0.78. The constant k depends on the sort and
configuration of the electrical phenomenon panel. The circuit voltage should be measured and
therefore the MPP determined in how for different close conditions. Usually, the system
disconnects the load sporadically to live FOC and calculate the in operation voltage. This
methodology has some clear disadvantages, temporary power loss is obvious. An alternate
Methodology would be to use one or a lot of watching cells, however they additionally ought
To be chosen and placed terribly rigorously to mirror truth open circuit voltage of the system.
17
Although this methodology is kind of easy and sturdy and does not need a microcontroller, the
constant solely allows a crude approximation of the MPP. Other algorithms can considerably
increase the highest power drawn from a similar PV installation
18
4.3 design parameters:
Vin 32.9
Vout 26.3
I load 7.61A
Fsw 200kHZ
Vout
K 0.8
Vin
Iripple 2.283
(a)Inductance Calculation:
In DC-DC converter Inductor and Capacitor act as filter to filter out higher order harmonics.
Hence Inductance removes the ripples from the output current.
For an inductor,
L Vin Vout *
D (4.3.1)
Fsw * I
ripple
L=11.56 µH
(b)Output Capacitor Calculation:
The output capacitor voltage is resultant sum of effective series resistance ripple voltage, load
current sag voltage supplied by capacitor as the inductor is discharged and capacitor’s effect
series inductance ripple voltage. As ESL is not rated on capacitor, so let’s assume ESL = 0.
But with increase in switching frequency, ESL value is taken under consideration.
I * T
Cout (4.3.2)
V I * ESR
C out 321.5F
This equation shows that it depends on ESR, C and ESL. With suitable approximation, we can
remove the other parameter from the equation and depends only on ESR which can be
controlled. ESR is taken from rating of amps rated capacitors.
19
(c)Input Capacitor
For critical condition of duty cycle, D = 80% and under this condition buck converter’s ripple
input current is half of the output load current. Input capacitor is selected based on ESR to
meet desired ripple voltage. Output ripple voltage is required more stringent compared to
input. Here, the maximum input voltage ripple was defined as 200 millivolts. For input
capacitor, ripple current rating is the important factor. Many times input ripple is more than
output.
T
Cin (4.3.3)
V
ripple ESR
I
ripple
C in 511F
20
SIMULATION MODEL
21
SIMULATION RESULT
30
25
20
voltage
15
10
0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35
time
5
current
0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35
time
22
300
250
200
power
150
100
50
0
0 0.05 0.1 0.15 0.2 0.25 0.3 0.35
time
23
CHAPTER 5
5.1. CONCLUSION
The project process simplified piecewise linear model for a PV array. In this work we have
chosen the KC-200GT solar panel as the power source and simulated its I-V and P-V
characteristics. The final simulation of PV-array along with MPPT and an impedance matching
device. A synchronous buck converter was used to implement the MPPT algorithm because of
its better efficiency in the fractional open-circuit voltage MPPT algorithm. The value of K was
set to 0.75, as usually K varies between 0.7 & 0.8. A PI controller is used to precisely control
the transient and study state response of the model under varying load as well as source
variation. The MPPT algorithm used was fractional open circuit voltage algorithm,
Which is very Simple and can be implemented cheaply. In open loop simulation of
synchronous buck converter there was significant error in the output value. This was due to the
losses in the circuit elements and the inability of the system to compensate for it. Switching
losses depend upon switching frequency chosen. In closed loop simulation the transient
response, steady state error etc. depend on the proportional and integral constant values.
FUTURE WORK
Study and simulation of advanced MPPT algorithms
Hardware implementation of the set-up
Simulation of effects of PI controller on transient response for different K p and K i
values
24
5.2 Reference
[1] Phang, J. C. H., D. S. H. Chan, and J. R. Phillips. "Accurate analytical method for the
extraction of solar cell model parameters," Electronics Letters, vol.20, no.10, 1984, pp.
406-408.
[2] Villalva, Marcelo Gradella, and Jonas Rafael Gazoli. "Comprehensive approach to
modeling and simulation of photovoltaic arrays." IEEE Trans. Power Electron. 24.5
(2009): 1198-1208.
[3] Altas, I. H., and A. M. Sharaf. "A photovoltaic array simulation model for matlab-simulink
gui environment." International Conference on Clean Electrical Power, ICCEP, IEEE,
2007.
[4] Tsai, Huan-Liang, Ci-Siang Tu, and Yi-Jie Su. "Development of generalized photovoltaic
model using matlab/simulink." Proceedings of the World Congress on Engineering and
Computer Science, 2008.
[5] Ishaque, Kashif, Zainal Salam, and Hamed Taheri. "Modeling and simulation of
photovoltaic (PV) system during partial shading based on a two-diode model. “Simulation
Modelling Practice and Theory, vol.19, no.7, 2011, pp. 1613-1626.
[6] Kawamura, Hajime, et al. "Simulation of I & V characteristics of a pv module with shaded
pv cells." Solar Energy Materials and Solar Cells, vol.75, no.3, 2003, pp. 613-621.
[7] G. E. Ahmad, H. M. S. Hussein, and H. H. El-Ghetany, “Theoretical analysis and
experimental verification of pv modules,” Renewable Energy, vol. 28, no. 8, pp. 1159–
1168, 2003.
[8] https://en.wikipedia.org/wiki/Buck_converter
[9] H. Altas and A. M. Sharaf, “A photovoltaic array simulation model for matlab-simulink
gui environment,” Proc. Of International Conference on Clean Electrical Power, ICCEP
’07, May 21-23, 2007.
[10] S. Rahman, M. Khallat, and B. Chowdhury, “A discussion on the diversity in the
applications of photo-voltaic system,” IEEE Trans., Energy Conversion,, vol. 3, pp. 738–
746, Dec. 1988.
25