Print version ISSN 0100-6916
Eng. Agríc. vol.32 no.4 Jaboticabal July/Aug. 2012
Mathematical analysis of maximum power generated by photovoltaic systems and fitting curves for standard test conditions1
Modelo matemático da potência máxima gerada por sistemas fotovoltaicos e curvas de ajuste às condições de teste em campo
Luís R. A. Gabriel FilhoI; Daniel dos S. Viais NetoII; Camila P. CremascoIII; Odivaldo J. SeraphimIV; Fernando de L. CaneppeleV
IMestre em Matemática, Doutor e Pós-Doutor em Agronomia, Prof. Doutor, UNESP - Univ Estadual Paulista, Tupã-SP, e FCA, UNESP, Botucatu - SP, firstname.lastname@example.org
IIMestre em Matemática, Prof. Associado, FATEC - Faculdade de Tecnologia de Presidente Prudente, email@example.com
IIIMestre em Matemática, Doutora em Agronomia, Pós-Doutora em Biometria, Profa. Doutora, FATEC - Faculdade de Tecnologia de Presidente Prudente, e FCA/UNESP, Botucatu - SP, firstname.lastname@example.org
IVProfessor Adjunto III, Faculdade de Ciências Agronômicas, Botucatu - SP, email@example.com
VEngo Eletricista, Prof. Doutor, Faculdade de Zootecnia e Engenharia de Alimentos, USP, Pirassununga - SP, firstname.lastname@example.org
The rural electrification is characterized by geographical dispersion of the population, low consumption, high investment by consumers and high cost. Moreover, solar radiation constitutes an inexhaustible source of energy and in its conversion into electricity photovoltaic panels are used. In this study, equations were adjusted to field conditions presented by the manufacturer for current and power of small photovoltaic systems. The mathematical analysis was performed on the photovoltaic rural system I-100 from ISOFOTON, with power 300 Wp, located at the Experimental Farm Lageado of FCA/UNESP. For the development of such equations, the circuitry of photovoltaic cells has been studied to apply iterative numerical methods for the determination of electrical parameters and possible errors in the appropriate equations in the literature to reality. Therefore, a simulation of a photovoltaic panel was proposed through mathematical equations that were adjusted according to the data of local radiation. The results have presented equations that provide real answers to the user and may assist in the design of these systems, once calculated that the maximum power limit ensures a supply of energy generated. This real sizing helps establishing the possible applications of solar energy to the rural producer and informing the real possibilities of generating electricity from the sun.
Keywords: rural electrification, photovoltaic cells, numerical methods.
A eletrificação rural é caracterizada pela dispersão geográfica da população, baixo consumo, alto investimento por consumidor e elevado custo operacional. Por outro lado, a radiação solar constitui-se numa inesgotável fonte energética, e para sua conversão em energia elétrica são utilizados painéis fotovoltaicos. Neste trabalho, foram determinadas equações ajustadas às condições de campo apresentadas pelo fabricante para corrente e potência de sistemas fotovoltaicos de pequeno porte. A análise matemática foi feita sobre o sistema fotovoltaico rural I-100 da ISOFOTON, com potência de 300 Wp, localizado na Fazenda Experimental Lageado da FCA/UNESP. Para o desenvolvimento de tais equações, procurou-se estudar o circuito elétrico de células fotovoltaicas, aplicar os métodos numéricos iterativos para a determinação de parâmetros elétricos e analisar possíveis erros na adequação das equações existentes da literatura à realidade. Sendo assim, foi proposta a simulação de um painel fotovoltaico por meio de equações matemáticas que foram ajustadas segundo os dados de radiação locais. Os resultados apresentam equações que fornecem respostas reais para o usuário e podem auxiliar no dimensionamento destes sistemas, uma vez que a potência máxima calculada garante um limite de fornecimento da energia gerada. Este real dimensionamento auxilia a determinação das possíveis aplicações da energia solar para o produtor rural, informando as possibilidades reais de geração de energia elétrica pelo Sol.
Palavras-chave: eletrificação rural, células fotovoltaicas, métodos numéricos.
The amount of energy the sun provides the Earth is considered very large, and in just one hour, energy, which is consumed by humanity during a year, reaches our planet. As the sun is a free renewable energy source and does not pollute the environment, the use of its energy for electricity supply represents a viable alternative, particularly where access to conventional electricity is difficult or uneconomic (MACIEL et al., 2008).
In this sense, global environmental concerns and growing demand for energy, along with the steady progress in renewable energy technologies are opening new opportunities for the use of renewable energy resources. The solar energy is the most abundant, inexhaustible and clean of all renewable energy resources until this date. The power of the sun intercepted by the Earth is about 1.8 × 1011 MW, which is much bigger than the current rate of all energy consumption (PARIDA et al., 2011).
According to MARINI & ROSSI (2005), the energy supply to isolated rural communities and remote areas has been a constant challenge, especially in underdeveloped or developing countries, due to the enormous economic and social disparities, and often the geographical dimensions of these countries. In agricultural engineering, mathematical models with real application in the area are being developed by using correction equations of numerical temporal variables (GABRIEL FILHO et al., 2011), whether by developing systems based on artificial intelligence (GABRIEL FILHO et al., 2011 be PEREIRA et al., 2008), and analytical equations for the optimization of photovoltaic systems (GABRIEL FILHO et al., 2011).
Small farms in Brazil, particularly those located in rural homes, according to their characteristic location often in remote areas of the distribution networks of conventional electricity, plus the fact that it has a low demand, do not justify, in most cases, the investment network extension: in this context, they present a suitable space for the application of the photovoltaic power generation, which provide a local solution as a renewable and reliable energy source and it can promote the self-sustainable social and economic development of rural populations, improving the quality of life of its inhabitants (FRAGA, 2009).
Excluding the manufacturing process of the modules, according ALVES & CAGNON (2010), another important factor in expanding the use of photovoltaic technology is the environmental appeal, because these systems are characterized not to change the thermal balance of the earth, operating without the need of greenhouse gases or acid gases for the formation of lakes reservoirs and the deforestation only occurs in large installations, which are often built in semi-arid or desert areas.
According to CASTRO (2008), photovoltaic cells (Figure 1) comprise a semiconductor material - silicon - which is added substances, named doping to create a suitable means for establishing the photovoltaic effect, i.e. direct conversion of power linked to the solar radiation into DC electrical power. Furthermore, CASTRO (2008) states that the cell is the smallest element of the photovoltaic system, producing electrical power typically of the order of 1.5 Wp (corresponding to a voltage of 0.5 V and a current of 3 A).
According to FRAGA (2009), photovoltaic module (Figure 2) is the basic unit of the entire system. The module is composed of solar cells connected in arrangements producing sufficient voltage and current for energy use. The number of cells connected in a module and their arrangement, which can be set and/or parallel, depends on the operating voltage and electric current desired. Once having the desired configuration, the assembly is encapsulated with special material that protects it from possible external damage.
Photovoltaic modules are also possible in the cluster so that the available power is adequate for the intended application. These panels associated with other devices that transform or store electrical energy so that it can be used, is called the photovoltaic system.
Applications of electricity generated by PV systems are small, as the system used in this study are generally very useful for farmers. The energy produced by these systems may be used among other ways to: pump water for human consumption, animals and small irrigation systems, operate electric motors of lower power, lighting, telephony, etc. These services are generally of great importance for the farmer, therefore equations of such systems for an adequate sizing of use is required.
The equations governing the operation of PV systems generally refer to some difference in the results of field measurements, possibly due to the fact of the theoretical approaches already has in its conception or may not contain the details of the module manufacturing equipment, or due to local environmental factors.
There are studies (CELIK & ACIKGOZ, 2007) searching for means of electrical parameters to compare the results of modeling and experimental verification, indicating which option best fits the model to the actual operating conditions.
However, an appropriate model to the actual operating conditions without the need for testing parameters combinations would be of fundamental importance for the correct use of an analytical model / mathematical models for PV.
The objective of this study is to determine the analytical forms of the curves of current and photovoltaic power, and calculate the points of power, voltage and maximum current, as well as create a mathematical model to field conditions related to the PV system located in the Experimental Farm Lageado of FCA/UNESP.
MATERIAL AND METHODS
The photovoltaic system used in this work has been installed at the Laboratory of Rural Power of the Engineering Department of UNESP, Faculdade de Ciências Agronômicas, Lageado, located in Botucatu, State of São Paulo.
For the PV system formation , concerning the generation of electricity from solar energy incident on the site, it was used three single-crystalline silicon photovoltaic modules of the I-100 Isofoton with nominal power of 100 Wp each (a total of 300 Wp installed ), the technical task of cells 36 in series, two cells in parallel, 0.72 m2 of cells, and TONC (800 Wm-2, 20 °C, AM 1.5, 1 ms‑1) 47 °C (Figure 3 and Table 1). The data supplied by the manufacturer of the panels for simulating the analytical models / preparation of curves developed and also to compare the results of the equations are given in Table 1.
The methodology of this study was developed in accordance to the proposed objectives, concerned with the development of deterministic mathematical models and obtaining practical results via numerical methods.
To develop governing equations of the strength of the photovoltaic systems, the circuitry of photovoltaic cells was studied, according to CASTRO (2008), it can be described according to Figure 4, where it can be seen the series resistance Rs[Ω], representing the voltage drop in circuit to the external contacts, and the parallel resistance Rp[Ω], representing the leakage current.
The current lighting (or photogerated) IL resulting from the loads generated by the illumination of the material represents the electric current generated by the beam of light radiation consisting of photons reaching the active surface of the cell (photovoltaic effect). This unidirectional current is constant for a given incident radiation. The p-n junction acts as a diode which is traversed by an inward-way ID which depends on the voltage V at the terminals of the cell. According to LORENZO (1994), the relationship between current and voltage V in the photovoltaic cell, passing through the load, is given by
I = IL - ID - IP , (1)
- current that goes through diode, A;
Io - dark current, A;
- thermal voltage, V;
m - diode ideality factor, dimensionless;
k = 1.38 10-23 - Boltzmann constant, J K-1;
Tc - temperature of the PV cell, K;
q = 1.6 10-19 - electron charge, C, and
- represents the leakage currents, A.
Equation (2) can appropriately represent the operation of a solar cell, but it has some inconveniences in practical use. The term has an implicit character; the current appears on both sides of the equation, forcing its resolution by iterative methods. Furthermore, the parameters, IL, IO, RS and RP are difficult to determine.
For the application of numerical methods is required operating parameters of the modules provided by the manufacturer measured in conditions called Standard Conditions (STC - Standard Test Conditions) or referred to herein by Field Test Conditions, defined by the irradiance of 1000 W m-2 (solar radiation received at the Earth's surface in daylight, at noon), 25 °C in the cell (cell efficiency is reduced with increasing temperature), the wind speed of 1 m s-1 and spectral distribution (BF) 1.5.
The numerical method (iterative) used was Newton's method (ARENALES & DAREZZO, 2008), which consists in the estimation of the roots of a function. For its application, it becomes initially, as shown in Figure 5, a point x0 in the field of function possibly close to the solution r1, and calculates the intercept of the tangent line to the function at that point with the x-axis whose abscissa is denoted by x1. If f(x1) is near zero, the point x1 is elected as the solution (root) of the equation. Otherwise, the process x1 is repeated until f(xi) is sufficient to almost zero, thus establishing the convergence of the process. If there is no convergence, or elects another value for the initial condition x0, or should look for another iterative method. In mathematical notation, this process can be written as follows:
Where in n indicates the nth algorithm interaction and f(xn) is the derivate function f in xn.
Thus, the application of numerical methods to the conditions STC, it is possible to obtain the points of maximum power (Pmax), and also the current (Imax) and voltage (Vmax) that produce such power, in which Pmax = Imax . Vmax. Because these numerical values having certain error when compared to the Field Test conditions, curves were obtained from the equations of the PV system adjusted to these conditions using analytical methods.
RESULTS AND DISCUSSION
According to LORENZO (1994), the effect of parallel resistance is negligible, then we can consider Ip = 0. The generated current and short circuit are equal, i.e. IL = Isc. Thus, in any work conditions, and . Therefore, the current (I) may be shown by means of:
In which VOC is the open circuit tension. As P = VI, so the Power (P) in the PV cell is given by:
The theorems below show ways of obtaining the series resistance and the maximum power of a solar cell that these results can be extended to modules by laws in series and parallel association.
Theorem 1: The series resistance of solar photovoltaic cell is given by:
By eq.(4), is given:
concluding the demonstration.
In the next theorem, we present an equation involving the maximum power for photovoltaic cells. This equation can be extended to modules and photovoltaic generators according to GABRIEL FILHO (2007).
Theorem 2: The maximum power (Pmax) of a photovoltaic cell, and also the voltage (Vmax) and current (Imax) in the cell that produce such power is given by:
As , then:
The maximum power point Pmax may be determined by doing . Therefore,
Therefore, the maximum power point occurs when .
This implies that . The demonstration of (7) and (9) follows (4) and (5).
It is necessary to determine some characteristic values of the photovoltaic system in relation to the photovoltaic cell, the photovoltaic module and solar generator used. The generator consists of 3 modules in parallel (Npm = 3), absence of modules connected in series (Nsm = 1 ). The association of the cells of the PV system was made with module of 72 cells, grouped in the form of 36 series cell (Nsc = 36 ) and 2 parallel cell (Npc = 2 ).
The data for the photovoltaic cell shown in Table 2, can be calculated through the laws of associations in series and parallel resistors with voltage and current associated (Table 3), except parameters and Voltage Thermal Resistance Series. These parameters are then calculated for the photovoltaic cell, and later obtained the values equivalent to the modules through associations with such laws.
In field testing conditions, the cell temperature is 25 °C and thus Tc = (273.16 + 25) K = 298.16 K. Observing that the Boltzmann constant is given by k = 1.3806505 10-23 J K-1 the electron charge e = 1.60205.10-19 C and it is given the diode ideality factor as (the optimum condition), it was calculated, respectively, the thermal voltage of the solar cell, and the generator module as follows:
By theorem 1 is given:
To determine the error between the data from Field Test Conditions (Table 3) and those obtained by the mathematical models, the maximum power points of photovoltaic solar cell will be determined through of Theorem 2. Using the values Rs and Vt presented in Table 3, and also some data provided by the manufacturer, it was possible to apply Newton's method to obtain the value Vmax of the following equation:
in which it is considered f(x) equal to the first member (12), x = Vmax and x0 = 0.4 obtaining the maximum voltage value to the cell Vmax c = 0.461705, according to Table 4.
By Theorem 1 and in accordance to the laws in series and parallel associations of resistances with associated voltage and current, it is defined Imax c = 3.0976 and pmax c = 1.4302.
By eqs.(4) and (5) and to the values listed in Table 3 for the standard conditions of operation of each photovoltaic cell module under consideration, it is possible to obtain the characteristic curve and the power curve for the cell (Figure 6). The characteristic curve and the power curve for the photovoltaic cell that currently exist in the model are given by:
In Table 5, we can see different values in quantities specified by the manufacturer and the values calculated using the mathematical model, with respective errors between them in STC conditions. Such differences are not due to numerical model applied as in Table 4, the error was null and not assigned to the observations performed under standard conditions by the manufacturer, compared to tests performed in laboratories.
The errors are due to analytical procedures used in the model described because the equations in the literature that have been used for the development of theorems hypotheses have already demonstrated that they do not consider certain electrical phenomena. This can be featured in eq.(4), which is a simplification of eq.(2) even with the setting of one of its variables, and also in eq.(10), which is considered as the optimum condition for the ideality factor LED (m = 1), which does not always occur under field conditions.
Thus, it can be assumed that all possible values being obtained by the proposed calculation will always have an error. Therefore, it is necessary when creating graphics for current and power to consider this approximation error for the adequacy of existing equations from literature to reality.
Shortly, it will be seen that the highest values for any possible errors of P, V and I, are observed errors (Table 5) to the maximum power point (Pmax), and voltage (Vmax) and current (Imax) that produce such power. For adaptation to Field Test Conditions (STC) is then established by Table 5 that the errors in the data provided by the manufacturer of voltage, current and power are -4.47%, 7.93% and 2.97% , respectively. Thus, the variables V*, I* and P* are created such magnitudes in field test conditions and therefore have V = 0.9553.V*, I = 1.0793.I* and P = P*/1.0297. Substituting these values in equations (13), have new equations for current and Power:
Figure 7 shows a comparison between the current and power curves of the photovoltaic cell voltage as a function of standard operating conditions using eq.(13) developed in this study by eq.(14) adjusted on the Field Test Conditions supplied by the manufacturer.
This study presented a methodology for the determination of curves and points of occurrence of maximum power. The scientific contribution of the methods examined is the adequacy of mathematical models for real results concerning the Field Test Condition.
Thus, the curves created from this study provide real answers for the user. The methods and theorems of the photovoltaic systems showed may be helpful for the further sizing applications of these systems, since the maximum power limit calculated ensures a supply of power output.
Moreover, this methodology is a tool to enable analysis of the scaling in real situations established by the manufacturer, helping to determine the potential applications of solar energy for rural producers, informing the real possibilities of generating electricity from the sun.
The authors thank the Graduate Program in Agronomy FCA/UNESP/ Botucatu supported by scientific and structural availability, Prope - Dean of Research at UNESP, for the financial support which enabled this publication and also to the anonymous collaborators, whose comments and suggestions contributed to important improvements for this study.
ALVES, A.F.; CAGNON, J.A. Desenvolvimento e avaliação de desempenho de um sistema de posicionamento automático para painéis fotovoltaicos. Energia na Agricultura, Botucatu, v.25, n.2, 2010, p.1-19. Disponível em: <http://126.96.36.199/index.php/energia/article/viewFile/52/34>. Acesso em: 1 jun. 2012. [ Links ]
ARENALES, S.; DAREZZO, A. Cálculo numérico: aprendizagem com apoio de software. São Paulo: Editora Tomson Learning, 2008. [ Links ]
CASTRO, R.M.G. Introdução à energia fotovoltaica. ed. 2.2. Lisboa: Universidade Técnica de Lisboa, Instituto Superior Técnico - DEEC/Secção de Energia, Energias Renováveis e Produção Descentralizada, 2008. [ Links ]
CBSOLAR. Centro Brasileiro para Desenvolvimento da Energia Solar Fotovoltaica. Porto Alegre: Faculdade de Física PUC-RS. Disponível em: <http://www.pucrs.br/cbsolar/ntsolar/index.htm>. Acesso em: 1 jun. 2012. [ Links ]
CELIK, A.N.; ACIKGOZ, N. Modelling and experimental verification of the operating current of mono-crystalline photovoltaic modules using four- and five-parameter models, Applied Energy, London, v.84, n.1, p.1-15, jan. 2007. Disponível em: <http://www.sciencedirect.com/science/article/pii/S0306261906000511>. Acesso em: 1 jun. 2012. [ Links ]
CRESESB/CEPEL. Centro de Referência para Energia Solar e Eólica Sérgio de Salvo Brito/Centro de Pesquisas de Energia Elétrica. Energia Solar - princípios e aplicações, 2009. Disponível em: <http://www.cresesb.cepel.br/tutorial/tutorial_solar.pdf>. Acesso em: 1 jun. 2012. [ Links ]
FRAGA, J.R.C.P. Análise do comportamento da bateria utilizada em sistemas fotovoltaicos de pequeno porte. 2009. 137 f. Tese (Doutorado em Energia na Agricultura) - Faculdade de Ciências Agronômicas, Botucatu, 2009. [ Links ]
GABRIEL FILHO, L.R.A. Análise e modelagem geométrica da potência gerada por um sistema híbrido solar fotovoltaico eólico. 2007. 136 f. Tese (Doutorado em Agronomia) - Faculdade de Ciências Agronômicas, Universidade Estadual Paulista, Botucatu, 2007. [ Links ]
GABRIEL FILHO, L.R.A.; CREMASCO, C.P.; SERAPHIM, O.J. Análise diferencial da potência máxima gerada por um sistema solar fotovoltaico. Energia na Agricultura, Botucatu, v.25, n.2, 2010. Disponível em: <http://188.8.131.52/index.php/energia/article/view/60>. Acesso em: 1 jun. 2012. [ Links ]
GABRIEL FILHO, L.R.A.; CREMASCO, C.P.; SERAPHIM, O.J.; CANEPPELE, F.L. Caracterização analítica e geométrica da metodologia geral de determinação de distribuições de Weibull para o regime eólico e suas aplicações. Engenharia Agrícola, Jaboticabal, v.31, n.1, 2011. Disponível em: <http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0100-69162011000100006&lng=en&nrm=iso>. Acesso em: 1 jun. 2012. [ Links ]
GABRIEL FILHO, L.R.A.; PEREIRA, D.F.; BARATO, F.F.; MAGALHÃES, M.M. Method of numerical correction of errors occasioned by delay of records during the monitoring of environmental variables of interest for animal production. Engenharia Agrícola, Jaboticabal, v.31, n.5, 2011. Disponível em: <http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0100-69162011000500001&lng=en&nrm=iso>. Acesso em: 1 jun. 2012. [ Links ]
GABRIEL FILHO, L.R.A.; CREMASCO, C.P.; PUTTI, F.F.; CHACUR, M.G.M. Application of fuzzy logic for the evaluation of livestock slaughtering. Engenharia Agrícola, Jaboticabal, v.31, n.4, 2011. Disponível em: <http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0100-69162011000400019&lng=pt&nrm=iso>. Acesso em: 1 jun. 2012. [ Links ]
LORENZO, E. Eletricidade solar: ingenieria de los sistemas fotovoltaicos. Madrid: Artes Gráficas Galas, 1994. 338 p. [ Links ]
MACIEL, N F.; LOPES, J.D.S.; LIMA, F.Z. Energia solar para o meio rural: fornecimento de eletricidade. Viçosa: CPT, 2008. 254 p. [ Links ]
MARINI, J.A.; ROSSI, L.A. Sistematização do dimensionamento técnico e econômico de sistemas fotovoltaicos isolados por meio de programa computacional. Engenharia Agrícola, Jaboticabal, v.25, n.1, 2005. Disponível em: <http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0100-69162005000100008&lng=en&nrm=iso>. Acesso em: 1 jun. 2012. [ Links ]
PARIDA, B.; INIYAN, S.; GOIC, R. A review of solar photovoltaic technologies. Renewable and Sustainable Energy Reviews, v.15, n.3, 2011. Disponível em: <http://www.sciencedirect.com/science/article/pii/S1364032110004016>. Acesso em: 1 jun. 2012. [ Links ]
PEREIRA, D.F.; BIGHI, C.A.; GABRIEL FILHO, L.R.A.; CREMASCO, C.P.C. Sistema fuzzy para estimativa do bem-estar de matrizes pesadas. Engenharia Agrícola, Jaboticabal, v.28, n.4, 2008. Disponível em: <http://www.scielo.br/scielo.php?script=sci_arttext&pid=S0100-69162008000400002&lng=pt&nrm=iso>. Acesso em: 1 jun. 2012. [ Links ]
Recebido pelo Conselho Editorial em: 26-5-2011
Aprovado pelo Conselho Editorial em: 1º-3-2012
1 Elaborado a partir da tese de doutorado do primeiro autor