Zener reverse breakdown is not modelled only avalance breakdown

[mikofski]

Current method is to use avalanche breakdown

Irbd = a_rbd*(Vd/Rsh)*(1-Vd/Vrbd)^-nRBD

The proposed method is to add a quadratic term to soften the shoulder of the breakdown like a Zener diode. It's a hybrid approach that should be easy to implement in PVMismatch.

Zener breakdown\left(1-\frac{V{diode}}{V{rbd}}\right)^{-n_{rbd}}")

Irbd = Isc*(a_rbd*(Vd/Rsh/Isc) + b_rbd*(Vd/Rsh/Isc)^2)*(1-Vd/Vrbd)^-nRBD

New 2^nd order term: b_rbd*(Vd/Rsh/Isc)^2

Another even better option would be to allow any expression to be used as the breakdown current, and implement this as an object that the user can pass to the PVcell class instead of specifying the parameters aRBD, bRBD, nRBD and VRBD.

The plot below shows the measured and simulated breakdown response.

Equation rendering by Online LaTeX Equation Editor API

[mikofski]

closed by @nzweibaum in 333b1bb6c443baf8f286ce85190f3daa1e2f5a8e