Simulate the behaviour of a grid connected PV system under different
conditions of irradiance and temperature. This function is used by the
prodGCPV function.
fProd(inclin, module, generator, inverter, effSys)If inclin is data.table or Gef object, the result
is a data.table object with these components:
cell temperature, \(^{\circ}{\rm C}\).
open circuit voltage, short circuit current, MPP voltage and current, respectively, in the conditions of irradiance and temperature provided by Inclin
voltage and current at the input of the inverter. If no voltage limitation occurs (according to the values of inverter$Vmax and inverter$Vmin), their values are identical to Vmpp and Impp. If the limit values are reached a warning is produced
power at the input of the inverter, W
power at the output of the inverter, W
efficiency of the inverter
A Gef object, a data.table object. In case of being data.table it must include a component named Gef (effective irradiance, W/m²) and another named Ta (ambient temperature, \(^{\circ}{\rm C}\)).
list of numeric values with information about the PV module,
Vocnopen-circuit voltage of the module at Standard Test Conditions (default value 51.91 volts.)
Iscnshort circuit current of the module at Standard Test Conditions (default value 14.07 amperes.)
Vmnmaximum power point voltage of the module at Standard Test Conditions (default value 43.76 volts.)
ImnMaximum power current of the module at Standard Test Conditions (default value 13.03 amperes.)
Ncsnumber of cells in series inside the module (default value 24)
Ncpnumber of cells in parallel inside the module (default value 6)
CoefVTcoefficient of decrement of voltage of each cell with the temperature (default value 0.0049 volts per celsius degree)
TONCnominal operational cell temperature, celsius degree (default value 45).
list of numeric values with information about the generator,
Nmsnumber of modules in series (default value 22)
Nmpnumber of modules in parallel (default value 130)
list of numeric values with information about the DC/AC inverter,
Kivector of three values, coefficients of the efficiency curve of the inverter (default c(0.002, 0.005, 0.008)), or a matrix of nine values (3x3) if there is dependence with the voltage (see references).
Pinvnominal inverter power (W) (default value 1.5e6 watts.)
Vmin, Vmaxminimum and maximum voltages of the MPP range of the inverter (default values 822 and 1300 volts)
Gumbminimum irradiance for the inverter to start (W/m²) (default value 20 W/m²)
list of numeric values with information about the system losses,
ModQualaverage tolerance of the set of modules (%), default value is 3
ModDispmodule parameter disperssion losses (%), default value is 2
OhmDCJoule losses due to the DC wiring (%), default value is 1.5
OhmACJoule losses due to the AC wiring (%), default value is 1.5
MPPaverage error of the MPP algorithm of the inverter (%), default value is 1
TrafoMTlosses due to the MT transformer (%), default value is 1
Displosses due to stops of the system (%), default value is 0.5
Oscar Perpiñán Lamigueiro, Francisco Delgado López.
Jantsch, M., Schmidt, H. y Schmid, J.: Results on the concerted action on power conditioning and control. 11th European photovoltaic Solar Energy Conference, 1992.
Baumgartner, F. P., Schmidt, H., Burger, B., Bründlinger, R., Haeberlin, H. and Zehner, M.: Status and Relevance of the DC Voltage Dependency of the Inverter Efficiency. 22nd European Photovoltaic Solar Energy Conference, 2007.
Alonso Garcia, M. C.: Caracterización y modelado de asociaciones de dispositivos fotovoltaicos. PhD Thesis, CIEMAT, 2005.
Perpiñán, O, Energía Solar Fotovoltaica, 2015. (https://oscarperpinan.github.io/esf/)
Perpiñán, O. (2012), "solaR: Solar Radiation and Photovoltaic Systems with R", Journal of Statistical Software, 50(9), 1-32, tools:::Rd_expr_doi("10.18637/jss.v050.i09")
fInclin,
prodGCPV,
fTemp.
library("data.table")
setDTthreads(2)
inclin <- data.table(Gef = c(200,400,600,800,1000),Ta = 25)
#using default values
fProd(inclin)
#Using a matrix for Ki (voltage dependence)
inv1 <- list(Ki = rbind(c(-0.00019917, 7.513e-06, -5.4183e-09),
c(0.00806, -4.161e-06, 2.859e-08),
c(0.02118, 3.4002e-05, -4.8967e-08)))
fProd(inclin, inverter = inv1)
#Voltage limits of the inverter
inclin <- data.table(Gef = 800,Ta = 30)
gen1 <- list(Nms = 10, Nmp = 11)
prod <- fProd(inclin,generator = gen1)
print(prod)
with(prod, Vdc * Idc / (Vmpp * Impp))
Run the code above in your browser using DataLab