Robust optimal design method for photovoltaic cells

Abstract

This invention relates to a robust optimal design method for photovoltaic cells. Firstly, the deterministic optimal model is established, which is solved by Monte Carlo method to obtain the maximum output power value of optimization objective and its corresponding design variable value, and then the design variable value obtained from deterministic optimization is deemed as the initial point of the mean value of the robust optimal design variable. Later, the robust optimal model is solved by Monte Carlo method in order to obtain the mean value of design variable, and then appropriate materials and manufacturing techniques are selected for corresponding photovoltaic components according to the design variable obtained, so as to achieve the robust optimal design of photovoltaic cells. In fact, this invention improves the output stability and reliability of photovoltaic cells.

Claims

1. A robust optimal design method for photovoltaic cells, which is characterized that it comprises the following steps: Step 1: Take material parameters of photovoltaic components as design variable x and temperature and radiation intensity in an actual working environment as design parameter z, and a specific distribution is as shown in Table 1: TABLE-US-00006 TABLE 1 Distribution of Design Variables and Parameters of Photovoltaic Cells Parameter Identi- Distribution Standard Name fication Unit Type Mean Value Deviation .sub. R.sub.x x.sub.1 Ω Normal [0.286, 0.656] 9.596 × 10.sup.−3 Distribution .sub.  R.sub.sh x.sub.2 Ω Normal [802.24, 1602.2] 22.2298 Distribution C x.sub.3 mA/C Normal [0.0059, 0.0061]   1.2 × 10.sup.−3 Distribution n x.sub.4 / Normal [1.18, 1.6]  0.003 Distribution S z.sub.1 W/m.sup.2 Normal 600 8 Distribution T z.sub.2 K Normal 303.15 6.063 Distribution wherein R.sub.s is equivalent series resistance of photovoltaic cells, R.sub.SH is equivalent shunt resistance of photovoltaic cells, C is temperature coefficient of short circuit current, n is diode ideality factor, S is radiation intensity, and T is surface temperature of photovoltaic cells, Take a maximum output power of photovoltaic cells as optimization objective and theoretical efficiency of a studied cell as constrained performance function to establish a deterministic optimal model as follows: Find x1, x2, x3, x4, and Max P(x,z), wherein Max P(x,z) is maximum output power of photovoltaic cells, given: s.t.g(x,z)=η(x,z)−0.159≤0, wherein s.t.g is objective function corresponding to inequality constraints about conversion efficiency and η is conversion efficiency of photovoltaic cells; 0.286≤x1≤0.656, 802.24≤x2≤1602.24; 0.0059≤x3≤0.0061, 1.18≤x4≤1.6; and Solve the deterministic optimal model by Monte Carlo method to obtain the maximum output power value Max P(x,z) of the optimization objective and its corresponding design variable value x={x1, x2, x3, x4}; Step 2: Take the design variable value obtained from deterministic optimization as an initial point of a mean value of robust optimal design variable (μx′={x.sub.1, x.sub.2, x.sub.3, x.sub.4}), and then obtain a robust optimal model based on the mean value and standard deviation after conversion as listed in Table 1, which is expressed as follows: Find μ.sub.x.sub.1, μ.sub.x.sub.2, μ.sub.x.sub.3, μ.sub.x.sub.4, $\mathrm{Min}\frac{{\sigma }_p\left(x,z\right)}{{\mu }_p\left(x,z\right)},$ wherein σ.sub.p is mean value of output power and μ.sub.p is standard deviation of output power, given: s.t.G(x,z)=μ.sub.g(x,z)≤0, wherein s.t.G is objective function corresponding to inequality constraints on mean value of the conversion efficiency and g is objective function corresponding to inequality constraints about conversion efficiency, 0.286≤μ.sub.x.sub.1≤0.656, 802.24≤μ.sub.x.sub.2≤1602.24 0.0059≤μ.sub.x.sub.3≤0.0061, 1.18≤μ.sub.x.sub.4≤1.6; and Step 3: Solve the robust optimal model by Monte Carlo method to obtain a mean value μ.sub.x{μ.sub.x.sub.1, μ.sub.x.sub.2, μ.sub.x.sub.3, μ.sub.x.sub.4}, of design variables, and then select corresponding materials and manufacturing techniques of photovoltaic components according to the obtained μx, so as to achieve robust optimal design for photovoltaic cells.

Description

BRIEF DESCRIPTION OF THE FIGURES

(1) FIG. 1 Equivalent Circuit Diagram of Photovoltaic Cells;

(2) FIG. 2 Flow Chart of Robust Optimal Design for Photovoltaic Cells.

DESCRIPTION OF THE INVENTION

(3) This invention takes the reliability and stability of the output of photovoltaic cells as objective, takes into account the randomness in the actual working environment of photovoltaic cells and the internal uncertainty in the material properties of photovoltaic components, deems the minimum ratio of standard deviation to mean value of the output power of photovoltaic cells under the co-influence of randomness and interval uncertainty as optimization objective, regards the given theoretical conversion efficiency as constrained performance function, establishes the robust optimal model, and solves the model by Monte Carlo method. The flow chart is shown in FIG. 2, which includes the following steps:

(4) Step 1: Establish the model of photovoltaic cells

(5) The equivalent circuit of photovoltaic cells is composed of the constant flow source that generates photo-generated current and a series of resistances (internal shunt resistance Rsh and series resistance Rs), as shown in FIG. 1.

(6) It can be obtained as follows according to the equivalent circuit of photovoltaic cells:
I=I.sub.ph−I.sub.D−I.sub.sh  (1)

(7) It can be obtained as follows according to the principle of diode:

(8) $\begin{array}{cc}{I}_D=I_o\left(\exp \frac{q\left(V+{\mathrm{IR}}_s\right)}{nKT}-1\right)& \left(2\right)\end{array}$

(9) It can be obtained according to the Kirchhoff principle:

(10) $\begin{array}{cc}{I}_{sh}=\frac{V+{\mathrm{IR}}_s}{R_{sh}}& \left(3\right)\end{array}$

(11) It can be obtained as follows according to Simultaneous Equations (1), (2) and (3):

(12) $\begin{array}{cc}I=I_{ph}-I_o\left(\exp \frac{q\left(V+{\mathrm{IR}}_s\right)}{nKT}-1\right)-\frac{V+{\mathrm{IR}}_s}{R_{sh}}& \left(4\right)\end{array}$

(13) In Equations (1), (2), (3) and (4): I refers to the output current of photovoltaic cells; Iph refers to photo-generated current; ID refers to the current of the diode; Ish refers to the current of equivalent shunt resistance that flows through photovoltaic cells; Io refers to reverse saturation current of the diode; V refers to the output voltage of photovoltaic cells; q refers to the electron charge; n refers to the impact factor of the diode; K refers to Boltzmann constant; T refers to the surface temperature of photovoltaic cells; Rs refers to the equivalent series resistance of photovoltaic cells; Rsh refers to the equivalent shunt resistance of photovoltaic cells.

(14) In order to facilitate the follow-up calculation, it is necessary to obtain the explicit equipment of current, so that the current item will not be contained on the right side of the equation. The explicit equation of current is obtained according to Lambert W function:

(15) $\begin{array}{cc}I=\frac{R_{\mathrm{sh}}\left(I_{ph}+I_o\right)-V}{R_s+R_{\mathrm{sh}}}-\frac{{\mathrm{nV}}_{\mathrm{th}}}{R_s}W\left(X\right)& \left(5\right)\end{array}$

(16) Where,

(17) $\begin{array}{cc}{V}_{\mathrm{th}}=\frac{nKT}{q},X=\frac{I_0{R}_{sh}{R}_s}{n{V}_{\mathrm{th}}\left(R_s+R_{\mathrm{sh}}\right)}\exp \left(\frac{R_{sh}\left(R_s{I}_{ph}+R_s{I}_o+V\right)}{n{V}_{th}\left(R_s+R_{sh}\right)}\right)& \left(6\right)\end{array}$

(18) In Equation (4), Iph and Io are model parameters, which are the functions of temperature and light intensity. The calculation methods are as follows:

(19) $\begin{array}{cc}{I}_{\mathrm{ph}}=\frac{S}{S_{\mathrm{ref}}}\left[I_{\mathrm{ph},\mathrm{ref}}+C\left(T-T_{\mathrm{ref}}\right)\right]& \left(7\right)\end{array}$$\begin{array}{cc}\frac{I_o}{I_{o,\mathrm{ref}}}={\left(\frac{T}{T_{ref}}\right)}^3\exp \left(\frac{q}{nK}\left(\frac{E_{g,\mathrm{ref}}}{T_{\mathrm{ref}}}-\frac{E_g}{T}\right)\right)& \left(8\right)\end{array}$

(20) In Equations (7) and (8), S refers to radiation intensity; T refers to cell temperature; C refers to temperature coefficient of short circuit current; Eg refers to material energy gap, whose temperature characteristic is E.sub.g/E.sub.g,ref=1−1.0002677(T−T.sub.ref); the subscript ref means under standard test conditions.

(21) The output power P can be obtained as shown in the following equation:

(22) $$\begin{gathered} P=\mathrm{IV}= \\ \frac{\left\{R_{\mathrm{sh}}\left\{\left\{S\left[I_{\mathrm{ph},\mathrm{ref}}+{\alpha }_{I,\mathrm{sc}}\left(T-T_{\mathrm{ref}}\right)\right]\right\}/S_{\mathrm{ref}}+\left\{I_{o,\mathrm{ref}}{T}^3\exp \left[-\left(E_8/T-E_{g,\mathrm{ref}}/T_{\mathrm{ref}}\right)/K\right]\right\}/T_{\mathrm{ref}}^3\right\}-V\right\}}{\left(R_s+R_{sh}\right)} \end{gathered}$$$\frac{\{{\mathrm{TK}}_n\star \mathrm{lambertw}\left(\begin{array}{c}0,\left\{I_{o,\mathrm{ref}}{R}_s{R}_{\mathrm{sh}}{T}^{2}q\exp \left\{\begin{array}{c}\frac{R_{\mathrm{sh}}q\left\{\begin{array}{c}\left\{R_{s}S\left[I_{\mathrm{ph},\mathrm{ref}}+{\alpha }_{I,\mathrm{sc}}\left(T-T_{\mathrm{ref}}\right)\right]\right\}/S_{\mathrm{ref}}+\\ \left\{I_{o,\mathrm{ref}}{R}_s{T}^3\exp \left[-\left(E_g/T-E_{g,\mathrm{ref}}/T_{\mathrm{ref}}\right)/K\right]\right\}/T_{\mathrm{ref}}^3\end{array}\right\}\}}{\mathrm{Tkn}\left(R_s+R_{\mathrm{sh}}\right)}+\\ \left[\frac{-\left(E_g/T-E_{g,\mathrm{ref}}/T_{\mathrm{ref}}\right)}{K}\right]\end{array}\right\}\right\}\\ /\left[T_{\mathrm{ref}}^3\mathrm{kn}\left(R_s+R_{\mathrm{sh}}\right)\right]\end{array}\right)\}}{R_{s}q}$

(23) The conversion efficiency η of photovoltaic cells is as follows:

(24) $\begin{array}{cc}\eta =\frac{I_m{V}_m}{AS}& \left(10\right)\end{array}$

(25) Where, Vm refers to the voltage of the maximum power point; Im refers to the current of the maximum power point; A refers to cell area.

(26) The output performance of photovoltaic cells is affected by the actual working environment. At the same time, photovoltaic components such as diode, resistance and battery panel in photovoltaic cells may face uncertainties in component parameters due to inevitable factors such as manufacturing and installation errors and dispersion of materials. Under the combined action of random working environment and parameter uncertainty, the output of photovoltaic cells may change a lot. Take the material parameters of photovoltaic components as design variable x and temperature and radiation intensity in the actual working environment as design parameter z, and the specific distribution is shown in Table 1 below.

(27) TABLE-US-00002 TABLE 1 Distribution of Design Variables and Parameters of Photovoltaic Cells Param- Identi- Distribution Standard eter fication Unit type Mean value deviation .sub. R.sub.s x.sub.1 Ω Normal [0.286, 0.656]  9.596 × 10.sup.−3 .sub.  R.sub.sh x.sub.2 Ω Normal [802.24, 1602.24] 22.2298 C x.sub.3 mA/C Normal [0.0059, 0.0061]    1.2 × 10.sup.−3 n x.sub.4 / Normal [1.18, 1.6]   0.03  s z.sub.1 W/m.sup.2 Normal 600    8     T z.sub.2 K Normal 303.15 6.063

(28) Step 2: Establish robust optimal design model

(29) Take the maximum output power of photovoltaic cells as optimization objective and theoretical efficiency of the studied cell as constrained performance function to establish the deterministic optimal model as follows:
Find x1x2x3x4
Max P(x,z)
s.t.g(x,z)=η(x,z)−0.159≤0
0.286≤x1≤0.656,802.24≤x2≤1602.24
0.0059≤x3≤0.0061,1.18≤x4≤1.6

(30) Given the fact that the mean value of design variables can better reflect the real results, take the mean value of original design variables as the design variable of the new model, so that the robust optimal model based on mean value and standard deviation can be obtained after conversion according to the above-mentioned deterministic optimal design model, which is expressed as follows:
μ.sub.x.sub.1,μ.sub.x.sub.2,μ.sub.x.sub.3,μ.sub.x.sub.4,

(31) 0$\mathrm{Min}\frac{{\sigma }_P\left(x,z\right)}{{\mu }_{\rho }\left(x,z\right)}$
s.t.G(x,z)=μg(x,z)≤0
0.286≤μ.sub.x.sub.1≤0.656,802.24≤μ.sub.x.sub.2≤1602.24
0.0059≤μ.sub.x.sub.3≤0.0061,1.18≤μ.sub.x.sub.4≤1.6

(32) Step 3: Solve the robust optimal model

(33) Monte Carlo method is an effective appropriate simulation method which solves uncertain numerical values based on random sampling. The method is a statistical method that mainly uses the design variables or parameters in uncertainty model problems for random sampling and then solves the probability of these problems with random samples. As a very direct and simple estimation method, the expressions for solving the mean value and variance are as follows:

(34) $\begin{array}{cc}{\mu }_Y\approx \frac{1}{N}\underset{i=1}{\overset{n}{.Math.}}y\left(x_i,z_i\right)& \left(11\right)\end{array}$$\begin{array}{cc}{\sigma }_Y^2\approx \frac{1}{N-1}\underset{i=1}{\overset{n}{.Math.}}{\left[y\left(x_i,z_i\right)-{\mu }_Y\right]}^2& \left(12\right)\end{array}$

(35) Where, N refers to sample size, xi and zi refer to the sample point of design variable and design parameters, respectively.

(36) The integral expression of solving the failure probability of constraints by Monte Carlo method is shown in Equation (12) below.

(37) $\begin{array}{cc}{P}_f={\int }_F^{}f\left(X\right)\mathrm{dX}=\int I_{F}f\left(X\right)\mathrm{dX}& \left(13\right)\end{array}$

(38) Where, f(X) refers to the joint probability density function with variable X={x1, x2, x3, x4, z1, z2}, F s, refers to the failure domain that does not meet constraint and Pf refers to failure probability.

(39) The specific solving process is as follows:

(40) Firstly, use Monte Carlo method to solve deterministic optimal model to obtain the maximum output power value of optimization objective and its corresponding design variable value x={x1, x2, x3, x4}, as shown in Line 2, Table 2. Take the design variable value obtained from deterministic optimization as the initial point μx′={x1, x2, x3, x4} of the mean value of robust optimal design variable, then solve the robust optimal model by Monte Carlo method to obtain the mean value μx={μx1, μx2, μx3, μx4} of design variables, as shown in Line 3, Table 2. Later, select appropriate materials and manufacturing techniques for corresponding photovoltaic components according to the obtained μx, so as to achieve robust optimal design for photovoltaic cells.

Embodiment 1

(41) Step 1: Establish the model of photovoltaic cells, determine design variables x1, x2, x3 and x4 and design parameters z1 and z2 of photovoltaic cells according to engineering experience, and set the value range of the design variable x, and the specific distribution is shown in Table 1.

(42) Step 2: Take the maximum output power of photovoltaic cells as optimization objective and theoretical efficiency as constrained performance function to establish the traditional deterministic optimal model; then take the mean value of original design variables as the design variable of the new model, minimum ratio of standard deviation to mean value of the output power of photovoltaic cells as optimization objective and the mean value of theoretical efficiencies as constrained performance function, and convert the deterministic optimal design model into the robust optimal model based on mean value and standard deviation.

(43) Step 3: Complete initial sampling of design variables and parameters through simple random sampling, take design variable x as independent control parameter, use Fmincon function in Matlab to obtain the response values of deterministic optimal objective function and constrained performance function, as well as the value of the corresponding design variable {0.286, 1602.24, 0.0061, 1.6}; take the design variable x={0.286, 1602.24, 0.0061, 1.6} obtained from deterministic optimization as the initial point of the mean value of robust optimal design variable, use Fmincon function in matlab to obtain the response value of robust optimal objective function and constrained performance function, as well as the value rangeμx={0.4778, 1601.2452, 0.005996, 1.3874} of the mean value of corresponding design variable; select appropriate materials and manufacturing techniques of corresponding photovoltaic components according to the obtained μx to achieve robust optimal design for photovoltaic cells.

(44) Both x and μx obey normal distribution. According to Equations (10) and (11), the mean value and variance of deterministic optimal and robust optimal constraints can be solved separately, as shown in Table 3. According to Equation (12), the failure probability can be separately solved when the deterministic optimal and robust optimal constraints fail, as shown in Table 4. (In order to facilitate calculation, N is set to 100000). Although the output power of photovoltaic cells under the deterministic optimal model is higher than that under the robust optimal model, the standard deviation and failure probability of constraints after robust optimization are smaller than those after deterministic optimization, and the output fluctuations become smaller, which improves the stability and reliability of the output of photovoltaic cells.

(45) TABLE-US-00003 TABLE 2 Calculation Results of Optimal Model Optimization method μ.sub.x.sub.1 μ.sub.x.sub.2 μ.sub.x.sub.3 μ.sub.x.sub.4 Deterministic 0.286  1602.24     0.61 × 10.sup.−2 1.6   optimization Robust 0.4778 1601.2452 0.5996 × 10.sup.−2 1.387 optimization

(46) TABLE-US-00004 TABLE 3 Results of Robust Optimization and Deterministic Optimization Optimization method μ.sub.g σ.sub.g Deterministic optimization 0.1592 6.4944 × 10.sup.−4 Robust optimization 0.159  6.3686 × 10.sup.−4

(47) TABLE-US-00005 TABLE 4 Failure Probability of Constraints under Different Optimization Results Failure Optimization method probability P.sub.f Deterministic optimization 0.6354 Robust optimization 0.4608