In [5]:
!pip install pvmismatch pvlib
Requirement already satisfied: pvmismatch in /usr/local/lib/python3.10/dist-packages (4.1)
Collecting pvlib
Downloading pvlib-0.10.4-py3-none-any.whl (29.5 MB)
━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━━ 29.5/29.5 MB 42.5 MB/s eta 0:00:00
Requirement already satisfied: numpy>=1.13.3 in /usr/local/lib/python3.10/dist-packages (from pvmismatch) (1.25.2)
Requirement already satisfied: matplotlib>=2.1.0 in /usr/local/lib/python3.10/dist-packages (from pvmismatch) (3.7.1)
Requirement already satisfied: scipy>=1.0.0 in /usr/local/lib/python3.10/dist-packages (from pvmismatch) (1.11.4)
Requirement already satisfied: future>=0.16.0 in /usr/local/lib/python3.10/dist-packages (from pvmismatch) (0.18.3)
Requirement already satisfied: six>=1.11.0 in /usr/local/lib/python3.10/dist-packages (from pvmismatch) (1.16.0)
Requirement already satisfied: pandas>=1.3.0 in /usr/local/lib/python3.10/dist-packages (from pvlib) (2.0.3)
Requirement already satisfied: pytz in /usr/local/lib/python3.10/dist-packages (from pvlib) (2023.4)
Requirement already satisfied: requests in /usr/local/lib/python3.10/dist-packages (from pvlib) (2.31.0)
Requirement already satisfied: h5py in /usr/local/lib/python3.10/dist-packages (from pvlib) (3.9.0)
Requirement already satisfied: contourpy>=1.0.1 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=2.1.0->pvmismatch) (1.2.1)
Requirement already satisfied: cycler>=0.10 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=2.1.0->pvmismatch) (0.12.1)
Requirement already satisfied: fonttools>=4.22.0 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=2.1.0->pvmismatch) (4.51.0)
Requirement already satisfied: kiwisolver>=1.0.1 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=2.1.0->pvmismatch) (1.4.5)
Requirement already satisfied: packaging>=20.0 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=2.1.0->pvmismatch) (24.0)
Requirement already satisfied: pillow>=6.2.0 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=2.1.0->pvmismatch) (9.4.0)
Requirement already satisfied: pyparsing>=2.3.1 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=2.1.0->pvmismatch) (3.1.2)
Requirement already satisfied: python-dateutil>=2.7 in /usr/local/lib/python3.10/dist-packages (from matplotlib>=2.1.0->pvmismatch) (2.8.2)
Requirement already satisfied: tzdata>=2022.1 in /usr/local/lib/python3.10/dist-packages (from pandas>=1.3.0->pvlib) (2024.1)
Requirement already satisfied: charset-normalizer<4,>=2 in /usr/local/lib/python3.10/dist-packages (from requests->pvlib) (3.3.2)
Requirement already satisfied: idna<4,>=2.5 in /usr/local/lib/python3.10/dist-packages (from requests->pvlib) (3.6)
Requirement already satisfied: urllib3<3,>=1.21.1 in /usr/local/lib/python3.10/dist-packages (from requests->pvlib) (2.0.7)
Requirement already satisfied: certifi>=2017.4.17 in /usr/local/lib/python3.10/dist-packages (from requests->pvlib) (2024.2.2)
Installing collected packages: pvlib
Successfully installed pvlib-0.10.4
In [6]:
import numpy as np
from matplotlib import pyplot as plt
from pvmismatch.contrib.gen_coeffs import gen_two_diode
from pvmismatch import *
def last_guess(sol):
isat1 = np.exp(sol.x[0])
isat2 = np.exp(sol.x[1])
rs = sol.x[2] ** 2.0
rsh = sol.x[3] ** 2.0
return isat1, isat2, rs, rsh
example generation of 2-diode coefficients¶
In [7]:
electrical_data = {
550: {"vmp": 37.2, "imp": 14.78, "voc": 44.8, "isc": 15.76}}
beta_voc = -0.0025 # per degC
alpha_isc = 0.0004 # per degC
nseries = 66
nparallel = 2
tc = 25.0 # degC
In [8]:
isc, voc = electrical_data[550]['isc'], electrical_data[550]['voc']
imp, vmp = electrical_data[550]['imp'], electrical_data[550]['vmp']
x, sol = gen_two_diode(isc, voc, imp, vmp, nseries, nparallel, tc)
sol
Out[8]:
message: The solution converged.
success: True
status: 1
fun: [ 5.430e-10 -6.277e-11 -4.501e-10 1.139e-14]
x: [-2.446e+01 -1.350e+01 6.605e-02 3.028e+00]
nfev: 8
njev: 1
fjac: [[-7.774e-01 -4.070e-02 -6.277e-01 1.692e-06]
[ 6.266e-01 -1.365e-01 -7.673e-01 -5.690e-06]
[ 5.448e-02 9.898e-01 -1.316e-01 5.477e-05]
[-1.897e-06 5.492e-05 -3.907e-06 -1.000e+00]]
r: [ 9.096e+00 1.509e+00 1.440e+02 -5.465e-02 6.668e-01
1.676e+02 -2.647e-03 1.486e+01 3.546e-02 3.714e-06]
qtf: [ 1.547e-09 -1.096e-08 -6.425e-10 1.440e-13]
In [9]:
x
Out[9]:
(2.3722544230539216e-11, 1.3665379526271328e-06, 0.004362363709123517, 9.168352828127766)
In [10]:
last_guess(sol)
Out[10]:
(2.3722544230539216e-11, 1.3665379526271328e-06, 0.004362363709123517, 9.168352828127766)
In [11]:
pvc = pvcell.PVcell(
Isat1_T0=x[0],
Isat2_T0=x[1],
Rs=x[2],
Rsh=x[3],
Isc0_T0=isc/nparallel,
alpha_Isc=alpha_isc)
In [12]:
pvc.Voc*nseries
Out[12]:
44.816799849648454
In [13]:
mpp = np.argmax(pvc.Pcell)
In [14]:
pvc.Vcell[mpp][0]*nseries
Out[14]:
37.18592977422984
In [15]:
pvc.Icell[mpp][0]*nparallel
Out[15]:
14.78556291473152
In [16]:
halfcell132 = pvmodule.crosstied_cellpos_pat([22, 22, 22], 2, partial=True)
In [17]:
pvm = pvmodule.PVmodule(cell_pos=halfcell132, pvcells=[pvc]*132)
In [18]:
# make some comparison plots
pvm.plotMod() # plot the pvmm module
plt.tight_layout()
f = plt.gcf()
ax = f.axes
ax[0].set_ylim([0, 20])
ax[0].plot(0, isc, 'or')
ax[0].plot(voc, 0, 'or')
ax[0].plot(vmp, imp, 'or')
ax[0].plot(0, pvm.Isc.mean()*nparallel, 'xk')
ax[0].plot(pvm.Voc.sum()/nparallel, 0, 'xk')
mpp = np.argmax(pvm.Pmod)
ax[0].plot(pvm.Vmod[mpp], pvm.Imod[mpp], 'xk')
ax[1].plot(vmp, imp*vmp, 'or')
ax[1].plot(pvm.Vmod[mpp], pvm.Pmod[mpp], 'xk')
Out[18]:
[<matplotlib.lines.Line2D at 0x7ebb5c8f4bb0>]
In [19]:
mpp = np.argmax(pvm.Pmod)
pvm.Vmod[mpp], pvm.Imod[mpp]
Out[19]:
(37.22876225039524, 14.765733698218524)
In [19]:
In [20]:
pvm.Vcell[:, mpp].shape
Out[20]:
(132,)
In [20]: