Next Article in Journal
Crack Identification in Cantilever Beam under Moving Load Using Change in Curvature Shapes
Previous Article in Journal
Color Variability Constrains Detection of Geometrically Perfect Mirror Symmetry
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Model of Series-Parallel Photovoltaic Arrays Using Double-Diode Model and Parallel Computing

by
Juan David Bastidas-Rodríguez
1,†,
Carlos Andrés Ramos-Paja
2,† and
Sergio Ignacio Serna-Garcés
3,*,†
1
Facultad de Ingeniería y Arquitectura, Universidad Nacional de Colombia, Manizales 170003, Colombia
2
Facultad de Minas, Universidad Nacional de Colombia, Medellín 050034, Colombia
3
Departamento de Electrónica y Telecomunicaciones, Instituto Tecnológico Metropolitano, Calle 54A No. 30-01, Medellín 050013, Colombia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 13 April 2022 / Revised: 14 May 2022 / Accepted: 16 May 2022 / Published: 18 June 2022

Abstract

:
Several applications require to estimate the power production of photovoltaic (PV) systems under partial shading conditions. For example, dynamic reconfiguration of the array connections is needed to maximize the power production under partial shading conditions, which requires estimating the power generated by the PV array in several possible configurations. Therefore, a fast and accurate modeling technique is needed to perform those calculations in practical times and with low estimation errors. To address those kinds of problems, this paper proposes a modeling approach based on the double-diode model to provide high accuracy at low voltage and low irradiance conditions, which are important for partial-shading analysis. Moreover, the proposed modeling technique is designed to be implemented in parallel processing devices; thus, the calculation time is much shorter in comparison with classical serial solutions. The proposed model is tested in terms of accuracy and speed, obtaining satisfactory results. Finally, the applicability of the parallel model in reconfiguration applications is demonstrated using an application example.

1. Introduction

Photovoltaic (PV) power generation is a growing option for providing electricity, where large PV fields are more common every year due to multiple reasons such as the need of reducing costs [1,2], fossil fuel shortage, greenhouse effect, the integration of generation systems to buildings, and the deployment of stand-alone systems for remote/rural areas. The energy harvest from PV fields depends on several factors [3], some intrinsic (panel degradation) and other extrinsic (solar irradiance, partial shading), where the fast estimation of the power production is crucial for control decisions like dynamic reconfiguration in the presence of partial shading [4].
Determine the power produced by a PV array requires the calculation of the voltage (V) and current (I) of the modules, which is needed to reconstruct the current vs. voltage (I–V) curves of the PV system; those curves can be later used to interpolate the power production at a given operation point [4]. Alternatively, a mathematical model of each PV section can be obtained using an Equivalent Electrical Circuit (EEC), where the single-diode and double-diode models are the most widely adopted [5]. The main difference between those models concerns the current produced by the recombination mechanism, which is only considered by the double-diode model. Comparing those models using criteria like root mean square error (RMSE), normalized RMSE, absolute error, relative error (RE), and mean RE, demonstrate the higher accuracy provided by the double-diode model. On the other hand, the number of parameters to be estimated for the single-diode model is lower, hence it is a simpler model [6]. However, the recombination current, which is only considered by the double-diode model, occurs in the depletion region (a low voltage region); this characteristic allows the double-diode model to provide better reproduction of the PV module behavior at low voltages [7], which is particularly important under partial shading conditions.
Several works highlight the differences between the single-diode and double-diode models. For example, in [8,9,10,11,12,13,14,15,16,17,18,19,20] different optimization methods are used to calculate the parameters of the models, such as FBIA (Forensic-Based Investigation Algorithm), GBO (Gradient-based optimizer), COA (Coyote Optimization Algorithm), CSO (Cat Swarm Optimization), BHCS (Biogeography-based Heterogeneous Cuckoo Search), PSO (Particle Swarm Optimization), GA (Genetic Algorithm), MLSHADE (Multi-strategy Success-History based Adaptive Differential Evolution), and IJAYA (Improved Jaya Algorithm); those works put into evidence the higher accuracy provided by the double-diode model over the single-diode one. Moreover, the double-diode model is also more accurate for different types of PV modules and under different operating conditions, such as low irradiance conditions and temperature variations [17,19,21,22].
The calculation of the power generated by a PV array requires information about the model and the interconnection topology of the PV modules. Regardless of the adopted model (single-diode or double-diode), the mathematical expression describing the I-V curve of the module is a transcendent equation, which is commonly solved using numerical methods. On the other hand, the modules-interconnection topologies are diverse, but the series-parallel structure is the most widely adopted for both small (urban) and large-scale installations [23,24,25].
There are several numerical solutions for the transcendental equations, where one simple method concerns the use of the Lambert-W function to transform those equations into explicit expressions [26]; however, the Lambert-W function is not available in several programming languages, and its accurate implementation requires complex numerical methods. Another strategy adopted to solve transcendental equations is the use of recursive numerical methods, where the most common options are the Newton-Raphson method, the bisection method, and the false position method [3,27,28]. However, the previous methods only seek one solution of the transcendental equation in a given range; hence, those methods are not recommended to solve systems with more than one solution. Other methods that have been used, due to their implementation in specialized software like Matlab, are the trust-region-dogleg [23,29], trust-region and Levenberg-Marquardt. However, only the trust-region-dogleg algorithm is designed to solve nonlinear equations, while the other ones are mainly aimed at minimizing the squares of a function.
The implicit equation of the single-diode model has been used to evaluate the power production of series-parallel PV arrays under mismatching conditions [29]. That method avoids evaluating the Lambert-W function, which decreases the calculation time, but keeps the relatively low accuracy of the single-diode model. A similar approach is reported in [23,29], both adopting the implicit expression of the PV current, but in this case, based on the double-diode model. In those works, the equation system representing the series-parallel PV array is solved using the trust-region-dogleg algorithm, and it is tested under non-homogeneous conditions, but only considering a series (classical) model execution. Similarly, the work reported in [3] proposes a methodology to solve the equation system of a series-parallel PV array based on the double-diode model. As in the previous case, this method is tested under non-homogeneous conditions, but the Newton-Raphson method used for the solution is not parallelized.
The Lambert-W function is a mathematical definition that has been used in the last years to solve the transcendental equation of PV models. For example, the work reported in [30] uses the Lambert-W function to estimate the power of a PV panel in a power system with a DC/DC converter, battery, load, and MPPT controller. In this case, the PV model is based on a single equation, and the authors do not provide the calculation analysis. Another approach is presented in [31], which uses the bisection method and a modified version of the Lambert-W function to speed up the calculation time by, approximately, a factor of 2 concerning other methods [32,33,34]; however, this solution is only tested with the fifth-order Schröder’s equation. A different approach is reported in [35], which uses an explicit expression for the double-diode model based on the Lambert-W function, and it is applied to solve a series-parallel PV array, including the action of the bypass diode. This solution provides an accurate solution, but the calculation algorithm is not concurrent.
A different solution was presented in [36], which proposes a discrete representation based on the Lambert-W function to calculate the power vs. voltage (P-V) and I-V curves of a PV array in different configurations and under homogeneous and non-homogeneous conditions. This solution, based on the single-diode model, is evaluated point-to-point to reach the MPPT, analyzing both regions of the diode conduction (negative and positive voltage) using the Lambert-W function [37,38]. The accuracy of the model is validated under partial shading, but there is not an analysis of performance and speed.
The difficulty of dealing with transcendent equations has been faced by transforming those expressions; for example, the authors of [39] propose a power-law model to calculate the I-V curve under real operating conditions for different types of PV modules, which avoids the use of the single-diode model transcendental equation. Although some parameters of this model lack physical meaning, the resulting expression has an analytical solution, high accuracy, and low computational burden. Another approach, reported in [40], proposes an explicit model based on the sine function, which only has a shape parameter. This model is aimed at constructing the I-V and P-V curves for the PV cells under different temperature and solar irradiance conditions, but the model is not tested under mismatched conditions. Similarly, the authors of [41] present a transformation of the transcendental equation of the PV cell model into a non-transcendental expression that reduces the computation time without impacting the accuracy, but this is applied only to the single diode model.
Other modern approaches to solve transcendental equations include artificial intelligence and advanced numerical methods. For example, in [42] is proposed a multilayer neural network to solve a transcendental equation; this method is tested in some problems modeled by transcendental expressions such as the diode equation. However, the authors do not provide a performance evaluation or a comparative analysis against other solutions. Similarly, Ref. [43] proposes some approximate analytical solutions for a transcendental equation similar to the PV model, where the goal is to obtain a simpler equation with an acceptable numerical evaluation. Despite some cases are evaluated, it is not confirmed that such an approximation improves the accuracy or speed in comparison with other solutions adopted to calculate the PV model.
To take advantage of the concurrent execution provided by FPGAs, the work in [44] uses an SoC formed by an ARM processor and an FPGA based on the Artix-7 logic, all in a single chip, for the simulation of a mismatched PV string. This solution is based on the system of equations reported in [45], which is formed by the transcendental equations derived from the single-diode model of the PV modules. Such a technique requires constructing a Jacobian matrix representing the nonlinear equations system, which can be symbolically inverted [45,46,47]; hence, the numerical method to solve the system is reduced to algebraic matrix operations. The authors split the simulation algorithm into three parts: the model update, the guess solution, and the Newton-Raphson numerical method. To perform a comparative analysis, the authors implemented three software versions. First, the complete solution runs on the ARM processor, which is the baseline; second, using the SIMD feature (Single-Instruction Multiple-Data) of the NEON coprocessor integrated into the ARM processor, the instructions are parallelized to calculate the currents and the differential resistances necessary to obtain the inverse of the Jacobian matrix; third, the numerical solution of the Lambert-W function [32] was described in VHDL to be implemented on the FPGA. The first two versions were software-only solutions, while the third one is a hardware/software hybrid solution. The second software version improved the execution time by 68% concerning the baseline (first solution), and the third solution improves the execution time by an additional 3%. In this work, the parallelization of the algorithm is the main element used to improve the execution time, since the performance of the hybrid solution depends mainly on the communication protocol between the hardware and the software stages. Finally, it is worth noting that the work uses the single-diode model and an approximation of the Lambert-W function, which prevents a complete parallelization of the system solution.
This paper proposes a solution aimed at providing high accuracy and speed in the calculation of the power production of a PV array since many PV installations are based on multiple strings. First, this new solution is based on the more accurate double-diode model of the PV modules, thus improving the accuracy over solutions based on the single-diode model for mismatched (partial-shading) conditions. Second, the proposed solution is focused on performing a parallel calculation of the strings, thus providing a fast calculation over solutions based on classical (series) calculation algorithms. Finally, the use of this parallel and accurate solution is tested under a real-world condition, which considers the reconfiguration of a partially-shaded PV array. The paper is organized in the following way: the PV array model considering the double-diode model for the PV modules is presented in Section 2. Section 3 reports the results and discusses the different test conditions. Finally, Section 4 presents the conclusions of the work.

2. Proposed Model

The general structure of the most widely used PV array configuration, named series-parallel (SP), is illustrated in Figure 1. The array is formed by M columns or strings, and each string is formed by a set of modules and a blocking diode connected in series. In turn, each module is formed by one or more series-connected submodules, where each submodule is formed by N s cells connected in series and a protection diode connected in anti-parallel, named bypass diode; such a submodule structure is depicted in Figure 2, where the double diode model represents the N s series-connected cells. Hence, each string of a SP array can be analyzed as a set of N submodules and one blocking diode connected in series. It is worth noting that the array voltage ( V a r r in Figure 1) is usually fixed by the power converter where the array is connected to; therefore, it is assumed as a known value.
This section introduces the mathematical descriptions of the main elements that form the PV array (submodules and blocking diodes) and the system of non-linear and implicit equations that describe each string, as proposed in [23]. Moreover, this section also presents the procedure to calculate the array current by using parallel computing.

2.1. Submodule and Blocking Diodes Model

The submodule is the main element of the PV array and it is represented by using the equivalent electric circuit illustrated in Figure 2, where the current source ( I p h ) represents the current generated by the photovoltaic effect and the diodes D 1 and D 2 represent the cells diffusion and recombination currents, respectively. Furthermore, the resistors R s and R h represent the ohmic losses and leakage currents, respectively, and the bypass diode ( B D ) limits the power dissipation in the submodule cells when they are under partial shading conditions [48].
The equation that relates the submodule’s current ( I s m ) and voltage ( V s m ) is shown in (1), which is a nonlinear and implicit equation. In (1) the parameters I s 1 and η 1 are the inverse saturation current and ideality factor of diode D 1 , respectively; while I s 2 , η 2 , I s b d and η b d correspond the same parameters for diodes D 2 and B D . Moreover, V t d is the thermal voltage of the PV cells and it is defined as V t d = k · T / q , where k is the Boltzmann constant, q is the electron charge, T is the cells temperature. Additionally, I b d is the current through the diode B D , as shown in (), and its thermal voltage is defined as V t b d = k · T b d / q , where T b d is the bypass diode temperature. It is worth noting that the submodule’s parameters depend on the irradiance and temperature, which is discussed in detail in [49].
f ( V s m , I s m ) = I s m + I p h I s 1 · e V s m + ( I s m I b d ) · R s N s · η 1 · V t d 1 I s 2 · e V s m + ( I s m I b d ) · R s N s · η 2 · V t d 1
V s m + ( I s m I b d ) · R s R h + I b d = 0
I b d = I s b d · e V s m n b d · V t b d 1
The other key element in the SP arrays is the blocking diode, which prevent opposite current flows through the string coming from the other strings of the array. The relation between the blocking diode current ( I b l k ) and voltage ( V b l k ) is introduced in (3), where I s b l k is the inverse saturation current, η b l k is the ideality factor, and V t b l k is the thermal voltage. Furthermore, V t b l k is defined as V t b l k = k · T b l k / q , where T b l k is the blocking diode temperature.
I b l k + I s b l k · e V b l k η b l k · V t b l k 1 = 0

2.2. Calculation of the String Current and Voltages

In a SP array formed by M strings connected in parallel, where each string has N submodules, i.e., an M × N array, each string is formed by N + 1 elements (N submodules and one blocking diode). Therefore, there are N + 2 unknown electrical variables that correspond to: the voltages of the N submodules ( V s m 1 V s m N ), the blocking diode voltage, and the string current ( I s t r ), which is shared by all the elements in the string.
To calculate the N + 2 unknowns, it is necessary to define the system of N + 2 nonlinear and implicit equations given in (4) [23], where the first N equations are obtained by applying (1) for each submodule, the N + 1 equation is obtained by using (3), and the last equation is obtained by applying KVL to the string including the array voltage. In expression (4), V is a vector formed by N + 1 voltages, where the first N elements correspond to the submodules voltages from top to bottom (i.e., V ( i ) = V s m , i i [ 1 N ] ) and the last element is the blocking diode voltage with the opposite polarity (i.e., V ( N + 1 ) = V b l k ). Moreover, I s t r is the string current and V a r r is the array voltage, which is shared by all the strings and it is assumed known. The system of equations F ( V , I s t r ) can be solved by using a numerical method or an optimization algorithm to obtain all the electrical unknowns in the string [50,51]. In this paper the Trust-Region Dogleg optimization algorithm [51] is used to solve F ( V , I s t r ) for each string, but any other suitable method can be adopted.
F ( V , I s t r ) = f 1 ( V ( 1 ) , I s t r ) f 2 ( V ( 2 ) , I s t r ) f N ( V ( N ) , I s t r ) I s t r + I s b l k · e V ( N + 1 ) η b l k · V t b l k 1 V ( 1 ) + + V ( N + 1 ) V a r r = 0

2.3. Calculation of the Array Current Using Parallel Computing

To calculate the array current it is necessary to solve (4) for each string of the array. At that point, all the strings currents are known and the array current can be calculated by using (5).
I a r r = k = 1 M I s t , k
It is worth noting that the system of nonlinear equations of each string does not depend on the electrical variables of the other strings; hence, all the strings in an M × N array can be solved independently. This particular characteristic is used in this paper to reduce the model solution time by using parallel computing, since each string can be solved by a different core or processing worker. This type of parallel computation can be performed using synchronous or asynchronous algorithms [52], where the workers in synchronous solutions wait for some synchronization condition, while the workers in asynchronous solutions can be used for other task when they finish the initial work. As it was previously described, the array current is calculated as the sum of all the strings currents; hence, that condition is an inherent synchronization point for the parallel calculation of the array current. Therefore, the solution proposed in this paper is a synchronous parallel algorithm.
The flow chart of the proposed model to solve SP arrays using parallel computing is shown in Figure 3, where the input data are: the array voltage ( V a r r ), N, M, and the model parameters of all the submodules and blocking diodes. Those parameters are organized in eleven N × M matrices ( M I p h , M I s 1 , M I s 2 , M η 1 , M η 2 , M V t d , M R s , M R h , M I s b d , M η b d , and M V t b d ), for the DDM parameters, and three 1 × M matrices for the blocking diode parameters ( M I s b l k , M η b l k , M V t b l k ). Moreover, N w is the number of available workers and S f and S l are auxiliary variables that define the first and last string to be evaluated, respectively.
The first step is to send the data of N w strings to the available workers. Then, each worker constructs the vector V and solves the system of nonlinear and implicit equations of the corresponding string. Once the string is solved, the current is saved to calculate I a r r later. When the N w workers have solved the N w strings, it is necessary to verify if all the strings has been solved; otherwise, the process is repeated until the M strings are solved. Finally, the array current is calculated by using (5).
At this point all the array is solved for a given value of V a r r , which can be used for different model based applications. For example, to generate the array I-V and P-V curves it is necessary to perform a voltage sweep of the array voltage from 0 V to the array open-circuit voltage and save the values of I a r r . Moreover, the proposed model can also be used to perform dynamic simulations of SP arrays in conjunction with power converters and maximum power point tracking techniques.

3. Result and Discussion

The proposed model is validated using experiments and detailed simulations, and those results are contrasted with the model proposed in [23], which is denominated Serial model. The details of the experimental results are introduced in the corresponding Section 3.4; hence, the following paragraphs provide information that apply for all the simulation scenarios.
This section includes, first, three evaluation scenarios to show the performance of the proposed model in terms of the errors in the I-V curve reproduction, the improvement in the calculation time for SP arrays with different number of strings, and an application example (model-based reconfiguration analysis). Moreover, all the simulations of the proposed model, from here on denominated Parallel model, are performed by using only two processing workers (i.e., N w = 2 ) to consider the worst case scenario for a parallel computing application. However, higher number of workers can be used depending on the computational platform.
The Trina Solar TSM-PD05 270 W module is used as a reference to construct the arrays considered in this section, whose main electrical characteristics in standard test conditions (STC) are: short-circuit current 9.22 A, open-circuit voltage 37.9 V, Maximum Power Point (MPP) current 8.73 A, MPP voltage 30.9 V, shot-circuit current temperature coefficient 0.05 %/K, open-circuit voltage temperature coefficient −0.32 %/K, number of cells 60, and number of submodules 3; therefore, each submodule has 20 cells (i.e., N s = 20 ). From those characteristics, and considering an irradiance of 1 kW/m 2 ( G = 1 kW/m 2 ) and a cell temperature of 55 C ( T = 328.15 K), it is possible to calculate the following DDM parameters, of each submodule, by using the procedure proposed in [53]: I p h = 9.3583 A, I s 1 = 18.846 η A, I s 2 = 16.316 μ A, η 1 = 1 , η 2 = 2 , R s = 0.1002 Ω , R h = 307.487 Ω . Moreover, the parameters of both bypass and blocking diodes models are calculated from the GF3045T bypass diode datasheet information, obtaining I s b d = I s b l k = 1 μ A, and η b d = η b l k = 0.2694 .
For sake of simplicity, it is assumed that all the submodules and blocking diodes have the same parameters, and also it is considered that the cells, bypass diodes, and blocking diodes operate at almost the same temperature (i.e., T T b d T b l k ) [47,54]. In any case, under experimental conditions, each parameter can be calculated independently using the procedure reported in [55,56]. Finally, the partial shading conditions are defined by using the matrix M I p h , since I p h is directly proportional to the irradiance reaching each submodule [57].

3.1. Error Evaluation under Different Partial Shading Conditions

This analysis considers a 15 × 2 (i.e., N = 15 and M = 2 ) array with two partial shading conditions ( C 1 and C 2 ) to evaluate the errors in the current prediction of both the Series and Parallel models regarding the array equivalent circuit implemented in Simulink. Condition C 1 represents a shading profile with two irradiance levels ( G = 1 kW/m 2 and G = 0.8 kW/m 2 ), while C 2 represents a random shading condition with submodule irradiance varying from G = 0.0286 kW/m 2 to G = 0.7765 kW/m 2 .
The I-V and P-V curves of the SP array, for operating conditions C 1 and C 2 , are introduced in Figure 4, where the Simulink results are plotted in red lines, while Serial and Parallel models are plotted in magenta and blue discontinuous lines, respectively. As expected, Serial and Parallel models provide the same results for both operating conditions, since the blue dashed line is over the magenta dashed line.
The errors in the current calculation of the Series and Parallel models are shown in Figure 5. As expected, the errors obtained with the two models are exactly the same, since both models solve the same systems of equations with the same algorithm; hence, both models provide the same Normalized Sum of Squared Errors (NSSE) for both operating conditions: 0.0139% for C 1 and 0.0032% for C 2 . In the left plot of Figure 5, it can be observed that the largest errors for condition C 1 are close to the open-circuit voltage (169 V) and the inflection point (71 V), where the errors are + 0.91 A (+93%) and −29.2 mA (−0.17%), respectively. Around those points the derivative of the current with respect to the voltage increases; therefore, the small differences between the equivalent array circuit and the models are translated into increments in the current errors. Additionally, it is important to mention that the high percentage in the current error close to the open-circuit voltage is expected since in this region of the I-V curve the currents are close to 0 A. Nevertheless, those current errors do not considerably affect the I-V curve reproduction, as observed in Figure 4. The increment of the current errors around the inflection points of the I-V curve is also evidenced in the right plot of Figure 5 for condition C 2 . In such a plot each increment in the magnitude of the current error corresponds to an inflection point in the I-V curve; however, those errors do not affect significantly the I-V reproduction. Moreover, for condition C 2 the errors move from −0.103 A (−0.75%) close to 16 V to +0.253 A (6.53%) close to 127 V; while the largest percentage error is again close to the open-circuit voltage (163 V) with 0.160 A (74.6%).

3.2. Calculation Time for Different Number of Strings

The main advantage of the Parallel model is the reduction of the calculation time without sacrificing the solution accuracy. Hence, the analysis performed in this section evaluates the calculation time of both Series and Parallel models to generate the I-V curve for a SP array with fifteen rows ( N = 15 ), and the number of columns varying from one to twenty (i.e., M = [ 1 20 ] ). All the I-V curves are generated with voltage steps of 2 V in V a r r , and each I-V curve is calculated six times to evaluate the average computation time and standard deviation for each array. In addition, the shading profiles for all the strings are randomly generated to represent a complex partial shading condition with a high number of inflection points in the I-V curves.
Figure 6 shows the calculation time and standard deviation of the I-V curves for both Serial and Parallel models considering SP arrays from 15 × 1 to 15 × 20 . It is worth noting that those results are obtained with two processing workers (i.e., N w = 2 ) for the Parallel model, which is the worst case scenario regarding the number of workers. Moreover, the exact mean values and standard deviations are presented in Table 1 as well as the average time reduction obtained with the Parallel model and its standard deviation.
For arrays with from 1 to 3 the calculation times are close but there is a reduction below between 2.5% and 28.0% with the Parallel model. Then, the calculation time of both models increases linearly regarding the number of strings in the array. However, the calculation time of the Serial model increases with a higher slope (1.73 s/string) than the one of the Parallel model (1.058 s/string). Hence, in general, the proposed model reduces the calculation time between 20.6% and 39.5% with respect to the Serial model with an average standard deviation of 2.13%. Furthermore, the standard deviations of the calculation times of both models confirm their repeatability.

3.3. Application Example: Exhaustive Model-Based Reconfiguration

The reduction in the calculation time provided by the Parallel model is useful in model-based reconfiguration systems. In those applications, the I-V curves must be calculated as fast as possible to reduce the reconfiguration time (period), since the shorter the reconfiguration period, the higher the power generation under partial shading conditions [58,59]. This is due to the continuous change in the shading profiles generated by the movements of the sun in the sky along the day and year.
The reconfigurable array considered in this application example is illustrated in Figure 7, which is a 15 × 4 array where the first two rows are reconfigurable (i.e., 2 × 4 ) and the rest of the array is fixed (i.e., 13 × 4 ). The eight reconfigurable submodules can be connected to any string through the reconfiguration matrix, which is implemented with switches (e.g., relays or MOSFETs) and controlled by the reconfiguration algorithm. This type of systems with a fixed and a reconfigurable parts have been previously reported in the literature because, sometimes, only part of the array is affected by partial shading conditions [60].
In this example, it is assumed that the reconfiguration algorithm evaluates all the possible configurations, i.e., it performs an exhaustive search. For the array explained before, it is possible to construct 105 different configurations with the eight reconfigurable submodules and the reconfiguration matrix. For each configuration the model generates the P-V curve, with voltage steps of 2 V, to identify the global maximum power point (GMPP). Then, the GMPP of each configuration is saved to determine the one that provides the highest GMPP. Once the best configuration is identified, the reconfiguration algorithm sets the reconfiguration matrix to implement such a best configuration into the PV array.
The analysis is performed for three different shading profiles defined through the M I p h matrix. The values of M I p h for the reconfigurable part are reported in Table 2, while the values of M I p h for the fixed part of the array are defined as follows: 5.6150 A for rows 3 to 9 (i.e., M I p h ( a , b ) = 5.6150 A a [ 3 9 ] b [ 1 4 ] ) and 2.8075 A for rows 10 to 15 (i.e., M I p h ( c , b ) = 2.8075 A c [ 10 15 ] b [ 1 4 ] ).
The 105 possible configuration were evaluated six times for each shading profile, with both the Serial and Parallel models, in order to calculate the average processing time values and the standard deviations shown in Table 3. The average time required by the Serial and Parallel models to evaluate the 105 configurations are 719.96 s (or 11.99 min) and 518.58 s (or 8.64 min), respectively, which implies a time reduction between 26.6% and 29.3% (with standard deviations lower than 0.23%) in the reconfiguration algorithm period when the Parallel model is used. Such a reduction indicates that the reconfiguration system using the Parallel model is able to enhance the array power production under mismatching conditions. Once again, it is important to highlight that these results were obtained with two processing workers ( N w = 2 ); hence, the calculation time can be further reduced by using as many workers as strings in the array (i.e., N w = M ).
In addition, in this application example the improvement in the maximum power extracted, comparing the best and worst configurations, varies from 7.39% for shading profile 1 to 4.36% for shading profile 3 (see Table 3). This indicates that the increment in the power production obtained by a reconfiguration system depends on the shading profile; however, such increment is translated into a higher energy production, a reduction in the project’s payback time, and a mitigation of the power generation variability of PV generators.

3.4. Experimental Results

The proposed model is also validated with experimental results for a 3 × 3 array formed by three modules Trina Solar TSM-PD05 270W connected in parallel, where each module has three series-connected submodules. The I-V curves of each string and the array are measured with a BK Precision electronic load, which is controlled from a PC through serial communication. The parameters of the DDM parameters used for the experimental validation are introduced in Table 4 and Table 5, where the first one shows the I p h values for each submodule, while the second table presents the rest of the DDM parameters. It is worth highlighting that the parameters of the three submodules in one string are assumed equal; that is why Table 5 shows one DDM parameter for each string.
Serial and Parallel models are used to reproduce the experimental results. The I-V curves of each string and the array are introduced in the left plot of Figure 8; while the error in the estimation of the array current is presented in the right plot of Figure 8. As in the previous examples, the results obtained with both models are the same, which entails the same errors. Moreover, the error varies from −0.62 A (−36.64%) to + 0.41 A (+4.79%), where the highest values of the errors are obtained close to the open-circuit voltage, as in the simulation results. However, both models correctly reproduce the electrical behavior of the PV array according to the results shown in Figure 8 and the NSSE obtained for both models ( N S S E = 0.012%).
Additionally, Serial and Parallel (with N w = 2 ) models are executed six times to validate the repeatability of the results obtaining that the average time required by the Serial and Parallel models are 2.1464 s and 2.0908 s. Hence, in this case, the time improvement obtained with the Parallel model is relatively low (2.590%), which is expected considering the low number of strings and rows in the array.
Finally, the simulation and experimental results introduced in this section show the applicability of the proposed model and the consistent time reduction obtained in different evaluating scenarios.

4. Conclusions

This paper has presented a parallel modeling strategy for SP arrays, which can reduce the processing time required to calculate the power production of the array without impacting the solution accuracy. This parallel solution is based on the double-diode model, it providing a higher accuracy in comparison with the classical single-diode model at low voltage and low irradiance conditions, which are important for partial shading analysis. Therefore, this solution is very useful for reconfiguration systems, since the parallel model will enable a fast modification of the array connections to improve the power production, and such a reconfiguration will be based on accurate predictions.
The previous claims were tested in several scenarios. The first test demonstrated the high accuracy of the parallel model since the prediction errors were very low in comparison with a commercial solution based on Simulink. The second test demonstrated that the proposed parallel model significantly improves the calculation speed over the classical serial model, which in this case reaches 40%; however, the calculation speed can be further improved by increasing the number of workers used to parallelize the calculation process. The third test demonstrated the applicability of the parallel model to reconfiguration systems, where the calculation speed was significantly improved, thus enabling the reduction of the reconfiguration period. Therefore, the PV array can be reconfigured more times during the day, which enables the optimization of the array connections to more partial shading conditions; hence, more power can be generated from the same reconfigurable array. The final test, based on experiments made on commercial PV panels, confirmed the model accuracy for real-world applications.
Another possible application for this parallel solution concerns the estimation of the potential power production of PV plants, which can be subjected to partial shading during some hours of the day (caused by seasons or unavoidable obstacles). This is critical for the economic viability analysis of PV plants, since partial shading could significantly decrease the power production; hence, the lifetime production of the PV plant is a decisive factor to consider. However, estimating the power production of a partially-shaded PV array requires a lot of calculations (and time), thus the proposed parallel solution could significantly reduce the processing time required. Such a fast calculation could enable to test much more possible installation sites and architectures since each test will take less time.
The main disadvantage of the parallel model is the requirement of multiple cores or workers, which is not needed to implement the classical serial model. However, processors available inside commercial computers have multiple cores (4, 6, 8, and more), thus the parallel model can take advantage of that hardware, which is neglected by the serial model. Another scenario concerns the embedded implementation adopted for real-time applications, which usually provide a single core; in any case, newer SoCs provide multiple cores [44], hence the parallel model can be implemented in embedded applications using those newer devices. Of course, multi-core embedded systems are costly in comparison with single-core devices, but the increased power production provided by the parallel model in reconfiguration systems, under partial shading conditions, could justify the increased cost of the embedded system. Otherwise, in the cases when the power increment is not enough, the serial solution will be more profitable; thus, an economic viability analysis is mandatory.

Author Contributions

Conceptualization, J.D.B.-R.; methodology, J.D.B.-R., C.A.R.-P. and S.I.S.-G.; software, J.D.B.-R. and S.I.S.-G.; validation, J.D.B.-R., C.A.R.-P. and S.I.S.-G.; formal analysis, J.D.B.-R., C.A.R.-P. and S.I.S.-G.; resources, J.D.B.-R., C.A.R.-P. and S.I.S.-G.; data curation, J.D.B.-R.; writing, J.D.B.-R., C.A.R.-P. and S.I.S.-G.; visualization, J.D.B.-R., C.A.R.-P. and S.I.S.-G.; funding acquisition, J.D.B.-R., C.A.R.-P. and S.I.S.-G. All authors have read and agreed to the published version of the manuscript.

Funding

The APC were covered by the Instituto Tecnológico Metropolitano under the project P21101 “Fortalecimiento y Consolidación del Grupo Automática, Electrónica y Ciencias Computacionales para responder a las necesidades de las industrias 4.0”.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the data generated in this study are reported in both figures and tables of the paper.

Acknowledgments

The authors thank the support of the Instituto Tecnológico Metropolitano, the Universidad Nacional de Colombia, and Minciencias (Fondo nacional de financiamiento para la ciencia, la tecnología y la innovación Francisco José de Caldas), under the projects P21101 “Fortalecimiento y Consolidación del Grupo Automática, Electrónica y Ciencias Computacionales para responder a las necesidades de las industrias 4.0” and “Estrategia de transformación del sector energético Colombiano en el horizonte de 2030—Energetica 2030”—“Generación distribuida de energía eléctrica en Colombia a partir de energía solar y eólica” (Code: 58838, Hermes: 38945).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. REN21. Renewables 2020 Global Status Report; Technical Report; Renewable Energy Policy Network for the 21st Century: Paris, France, 2020. [Google Scholar]
  2. Shukla, A.K.; Sudhakar, K.; Baredar, P. Recent advancement in BIPV product technologies: A review. Energy Build. 2017, 140, 188–195. [Google Scholar] [CrossRef]
  3. Kermadi, M.; Chin, V.J.; Mekhilef, S.; Salam, Z. A fast and accurate generalized analytical approach for PV arrays modeling under partial shading conditions. Sol. Energy 2020, 208, 753–765. [Google Scholar] [CrossRef]
  4. Serna-Garcés, S.; Bastidas-Rodríguez, J.; Ramos-Paja, C. Reconfiguration of Urban Photovoltaic Arrays Using Commercial Devices. Energies 2015, 9, 2. [Google Scholar] [CrossRef] [Green Version]
  5. Petrone, G.; Ramos-Paja, C.A.; Spagnuolo, G. Photovoltaic Sources Modeling; John Wiley & Sons, Ltd.: Chichester, UK, 2017. [Google Scholar]
  6. Abbassi, R.; Abbassi, A.; Jemli, M.; Chebbi, S. Identification of unknown parameters of solar cell models: A comprehensive overview of available approaches. Renew. Sustain. Energy Rev. 2018, 90, 453–474. [Google Scholar] [CrossRef]
  7. Hassan, H.M.; Sulaiman, I.A.; Singh, M.; Gupta, P.; Sahni, M.; Rajoriya, M. Exploratory modelling and simulation of InGaN as a basis of comparison between single and double diode models of PV devices. Mater. Today Proc. 2020, 49, 3112–3116. [Google Scholar] [CrossRef]
  8. Shaheen, A.M.; Ginidi, A.R.; El-Sehiemy, R.A.; Ghoneim, S.S.M. A Forensic-Based Investigation Algorithm for Parameter Extraction of Solar Cell Models. IEEE Access 2021, 9, 1–20. [Google Scholar] [CrossRef]
  9. Ismaeel, A.A.K.; Houssein, E.H.; Oliva, D.; Said, M. Gradient-Based Optimizer for Parameter Extraction in Photovoltaic Models. IEEE Access 2021, 9, 13403–13416. [Google Scholar] [CrossRef]
  10. Yang, X.; Gong, W. Opposition-based JAYA with population reduction for parameter estimation of photovoltaic solar cells and modules. Appl. Soft Comput. 2021, 104, 107218. [Google Scholar] [CrossRef]
  11. Yang, B.; Wang, J.; Zhang, X.; Yu, T.; Yao, W.; Shu, H.; Zeng, F.; Sun, L. Comprehensive overview of meta-heuristic algorithm applications on PV cell parameter identification. Energy Convers. Manag. 2020, 208, 112595. [Google Scholar] [CrossRef]
  12. Diab, A.A.Z.; Sultan, H.M.; Do, T.D.; Kamel, O.M.; Mossa, M.A. Coyote Optimization Algorithm for Parameters Estimation of Various Models of Solar Cells and PV Modules. IEEE Access 2020, 8, 111102–111140. [Google Scholar] [CrossRef]
  13. Hao, Q.; Zhou, Z.; Wei, Z.; Chen, G. Parameters Identification of Photovoltaic Models Using a Multi-Strategy Success-History-Based Adaptive Differential Evolution. IEEE Access 2020, 8, 35979–35994. [Google Scholar] [CrossRef]
  14. Abdel-Basset, M.; Mohamed, R.; Mirjalili, S.; Chakrabortty, R.K.; Ryan, M.J. Solar photovoltaic parameter estimation using an improved equilibrium optimizer. Sol. Energy 2020, 209, 694–708. [Google Scholar] [CrossRef]
  15. Chen, X.; Yu, K. Hybridizing cuckoo search algorithm with biogeography-based optimization for estimating photovoltaic model parameters. Sol. Energy 2019, 180, 192–206. [Google Scholar] [CrossRef]
  16. Pourmousa, N.; Ebrahimi, S.M.; Malekzadeh, M.; Alizadeh, M. Parameter estimation of photovoltaic cells using improved Lozi map based chaotic optimization Algorithm. Sol. Energy 2019, 180, 180–191. [Google Scholar] [CrossRef]
  17. Et-Torabi, K.; Nassar-Eddine, I.; Rmaily, R.; Obbadi, A.; Errami, Y.; Sahnoun, S.; El Fajri, A.; Agunaou, M. Photovoltaic parameter extraction methods for single and double diode models under real changing environmental conditions. In Proceedings of the 2018 6th International Renewable and Sustainable Energy Conference (IRSEC 2018), Rabat, Morocco, 5–8 December 2018. [Google Scholar] [CrossRef]
  18. Guo, L.; Meng, Z.; Sun, Y.; Wang, L. Parameter identification and sensitivity analysis of solar cell models with cat swarm optimization algorithm. Energy Convers. Manag. 2016, 108, 520–528. [Google Scholar] [CrossRef]
  19. Tamrakar, V.; Gupta, S.; Sawle, Y. Study of characteristics of single and double diode electrical equivalent circuit models of solar PV module. In Proceedings of the 2015 International Conference on Energy Systems and Applications, Pune, India, 30 October–1 November 2015; pp. 312–317. [Google Scholar] [CrossRef]
  20. Alrahim Shannan, N.M.A.; Yahaya, N.Z.; Singh, B. Single-diode model and two-diode model of PV modules: A comparison. In Proceedings of the 2013 IEEE International Conference on Control System, Computing and Engineering, Penang, Malaysia, 29 November–1 December 2013; pp. 210–214. [Google Scholar] [CrossRef]
  21. Elgohary, R.; Elela, A.A.A.; Elkholy, A. Electrical Characteristics Modeling for Photovoltaic Modules Based on Single and Two Diode Models. In Proceedings of the 2018 Twentieth International Middle East Power Systems Conference (MEPCON), Cairo, Egypt, 18–20 December 2018; pp. 685–688. [Google Scholar] [CrossRef]
  22. Suthar, M.; Singh, G.K.; Saini, R.P. Comparison of mathematical models of photo-voltaic (PV) module and effect of various parameters on its performance. In Proceedings of the 2013 International Conference on Energy Efficient Technologies for Sustainability, Nagercoil, India, 10–12 April 2013; pp. 1354–1359. [Google Scholar] [CrossRef]
  23. Pérez Archila, L.M.; Bastidas Rodríguez, J.D.; Correa, R. Implicit modelling of series-parallel photovoltaic arrays using double-diode model and its solution. Sol. Energy 2021, 214, 131–137. [Google Scholar] [CrossRef]
  24. Shankar, N.; SaravanaKumar, N. Reduced Partial shading effect in Multiple PV Array configuration model using MPPT based Enhanced Particle Swarm Optimization Technique. Microprocess. Microsyst. 2020, 2020, 103287. [Google Scholar] [CrossRef]
  25. Dhimish, M.; Holmes, V.; Mehrdadi, B.; Dales, M.; Chong, B.; Zhang, L. Seven indicators variations for multiple PV array configurations under partial shading and faulty PV conditions. Renew. Energy 2017, 113, 438–460. [Google Scholar] [CrossRef] [Green Version]
  26. Hayes, B. Why W? Am. Sci. Mag. 2005, 90, 104–108. [Google Scholar] [CrossRef]
  27. Md Golam, M.; Chakraborty, A.; Bhuiyan, M.A.A. A Robust Method for Solving Transcendental Equations. IJCSI Int. J. Comput. Sci. Issues 2012, 9, 413–419. [Google Scholar]
  28. Li, Z. Numerical evaluations of transcendental equations for transient experiments. Int. J. Hydrogen Energy 2012, 37, 8118–8122. [Google Scholar] [CrossRef]
  29. Bastidas-Rodriguez, J.D.; Cruz-Duarte, J.M.; Correa, R. Mismatched Series–Parallel Photovoltaic Generator Modeling: An Implicit Current–Voltage Approach. IEEE J. Photovolt. 2019, 9, 768–774. [Google Scholar] [CrossRef]
  30. Tripathy, M.; Kumar, M.; Sadhu, P. Photovoltaic system using Lambert W function-based technique. Sol. Energy 2017, 158, 432–439. [Google Scholar] [CrossRef]
  31. Fukushima, T. Precise and fast computation of Lambert W-functions without transcendental function evaluations. J. Comput. Appl. Math. 2013, 244, 77–89. [Google Scholar] [CrossRef]
  32. Veberič, D. Lambert W function for applications in physics. Comput. Phys. Commun. 2012, 183, 2622–2628. [Google Scholar] [CrossRef] [Green Version]
  33. Chapeau-Blondeau, F.; Monir, A. Numerical evaluation of the Lambert W function and application to generation of generalized Gaussian noise with exponent 1/2. IEEE Trans. Signal Process. 2002, 50, 2160–2165. [Google Scholar] [CrossRef] [Green Version]
  34. Fritsch, F.N.; Shafer, R.E.; Crowley, W.P. Solution of the transcendental equation we w = x. Commun. ACM 1973, 16, 123–124. [Google Scholar] [CrossRef]
  35. Lun, S.; Wang, S.; Yang, G.; Guo, T. A new explicit double-diode modeling method based on Lambert W-function for photovoltaic arrays. Sol. Energy 2015, 116, 69–82. [Google Scholar] [CrossRef]
  36. Fathabadi, H. Lambert W function-based technique for tracking the maximum power point of PV modules connected in various configurations. Renew. Energy 2015, 74, 214–226. [Google Scholar] [CrossRef]
  37. Psarros, G.N.; Batzelis, E.I.; Papathanassiou, S.A. Partial Shading Analysis of Multistring PV Arrays and Derivation of Simplified MPP Expressions. IEEE Trans. Sustain. Energy 2015, 6, 499–508. [Google Scholar] [CrossRef]
  38. Batzelis, E.I.; Routsolias, I.A.; Papathanassiou, S.A. An Explicit PV String Model Based on the Lambert $W$ Function and Simplified MPP Expressions for Operation Under Partial Shading. IEEE Trans. Sustain. Energy 2014, 5, 301–312. [Google Scholar] [CrossRef]
  39. Zhang, C.; Zhang, Y.; Su, J.; Gu, T.; Yang, M. Modeling and Prediction of PV Module Performance Under Different Operating Conditions Based on Power-Law I–V Model. IEEE J. Photovolt. 2020, 10, 1816–1827. [Google Scholar] [CrossRef]
  40. Oulcaid, M.; El Fadil, H.; Ammeh, L.; Yahya, A.; Giri, F. One shape parameter-based explicit model for photovoltaic cell and panel. Sustain. Energy Grids Netw. 2020, 21, 100312. [Google Scholar] [CrossRef]
  41. Mahmoud, Y.; El-Saadany, E.F. A Photovoltaic Model with Reduced Computational Time. IEEE Trans. Ind. Electron. 2014, 62, 3534–3544. [Google Scholar] [CrossRef]
  42. Jeswal, S.K.; Chakraverty, S. Solving transcendental equation using artificial neural network. Appl. Soft Comput. J. 2018, 73, 562–571. [Google Scholar] [CrossRef]
  43. Barsan, V. Siewert solutions of transcendental equations, generalized Lambert functions and physical applications. Open Phys. 2018, 16, 232–242. [Google Scholar] [CrossRef]
  44. Guarino, A.; Monmasson, É.; Spagnuolo, G. SoC-based embedded real-time simulation of mismatched photovoltaic strings. Math. Comput. Simul. 2021, 184, 267–281. [Google Scholar] [CrossRef]
  45. Petrone, G.; Ramos-Paja, C. Modeling of photovoltaic fields in mismatched conditions for energy yield evaluations. Electr. Power Syst. Res. 2011, 81, 1003–1013. [Google Scholar] [CrossRef]
  46. Orozco-Gutierrez, M.; Ramirez-Scarpetta, J.; Spagnuolo, G.; Ramos-Paja, C. A technique for mismatched PV array simulation. Renew. Energy 2013, 55, 417–427. [Google Scholar] [CrossRef]
  47. Petrone, G.; Spagnuolo, G.; Vitelli, M. Analytical model of mismatched photovoltaic fields by means of Lambert W-function. Sol. Energy Mater. Sol. Cells 2007, 91, 1652–1657. [Google Scholar] [CrossRef]
  48. Satpathy, R.; Pamuru, V. Manufacturing of crystalline silicon solar PV modules. In Solar PV Power; Elsevier: Amsterdam, The Netherlands, 2021; pp. 135–241. [Google Scholar] [CrossRef]
  49. Gong, Y.; Wang, Z.; Lai, Z.; Jiang, M. TVACPSO-assisted analysis of the effects of temperature and irradiance on the PV module performances. Energy 2021, 227, 120390. [Google Scholar] [CrossRef]
  50. Lindfield, G.; Penny, J. Solution of Non-Linear Equations. In Numerical Methods; Elsevier: Amsterdam, The Netherlands, 2019; pp. 157–190. [Google Scholar] [CrossRef]
  51. Nocedal, J.; Wright, S.J. Trust-Region Methods. In Numerical Optimization, 2nd ed.; Springer: New York, NY, USA, 2006; Chapter 4; pp. 66–100. [Google Scholar] [CrossRef]
  52. Li, S.E.; Wang, Z.; Zheng, Y.; Sun, Q.; Gao, J.; Ma, F.; Li, K. Synchronous and asynchronous parallel computation for large-scale optimal control of connected vehicles. Transp. Res. Part C Emerg. Technol. 2020, 121, 102842. [Google Scholar] [CrossRef]
  53. Hejri, M.; Mokhtari, H.; Azizian, M.R.; Ghandhari, M.; Söder, L. On the parameter extraction of a five-parameter double-diode model of photovoltaic cells and modules. IEEE J. Photovolt. 2014, 4, 915–923. [Google Scholar] [CrossRef]
  54. Lee, C.G.; Shin, W.G.; Lim, J.R.; Kang, G.H.; Ju, Y.C.; Hwang, H.M.; Chang, H.S.; Ko, S.W. Analysis of electrical and thermal characteristics of PV array under mismatching conditions caused by partial shading and short circuit failure of bypass diodes. Energy 2021, 218, 119480. [Google Scholar] [CrossRef]
  55. Accarino, J.; Petrone, G.; Ramos-Paja, C.A.; Spagnuolo, G. Symbolic algebra for the calculation of the series and parallel resistances in PV module model. In Proceedings of the 2013 International Conference on Clean Electrical Power (ICCEP), Alghero, Italy, 11–13 June 2013; pp. 62–66. [Google Scholar] [CrossRef]
  56. Kumari, P.A.; Geethanjali, P. Parameter estimation for photovoltaic system under normal and partial shading conditions: A survey. Renew. Sustain. Energy Rev. 2018, 84, 1–11. [Google Scholar] [CrossRef]
  57. Bhol, R.; Dash, R.; Pradhan, A.; Ali, S.M. Environmental effect assessment on performance of solar PV panel. In Proceedings of the 2015 International Conference on Circuits, Power and Computing Technologies [ICCPCT-2015], Kumaracoil, India, 19–20 March 2015; pp. 1–5. [Google Scholar] [CrossRef]
  58. Yousri, D.; Babu, T.S.; Beshr, E.; Eteiba, M.B.; Allam, D. A Robust Strategy Based on Marine Predators Algorithm for Large Scale Photovoltaic Array Reconfiguration to Mitigate the Partial Shading Effect on the Performance of PV System. IEEE Access 2020, 8, 112407–112426. [Google Scholar] [CrossRef]
  59. Kumar, N.; Hussain, I.; Singh, B.; Panigrahi, B.K. Rapid MPPT for Uniformly and Partial Shaded PV System by Using JayaDE Algorithm in Highly Fluctuating Atmospheric Conditions. IEEE Trans. Ind. Inform. 2017, 13, 2406–2416. [Google Scholar] [CrossRef]
  60. Velasco-Quesada, G.; Guinjoan-Gispert, F.; Pique-Lopez, R.; Roman-Lumbreras, M.; Conesa-Roca, A. Electrical PV Array Reconfiguration Strategy for Energy Extraction Improvement in Grid-Connected PV Systems. IEEE Trans. Ind. Electron. 2009, 56, 4319–4331. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Structure of a Series-Parallel (SP) photovoltaic array.
Figure 1. Structure of a Series-Parallel (SP) photovoltaic array.
Computation 10 00100 g001
Figure 2. Equivalent circuit of the double-diode model used to represent each submodule.
Figure 2. Equivalent circuit of the double-diode model used to represent each submodule.
Computation 10 00100 g002
Figure 3. Flow chart of the solution of an M × N Series-Parallel array using parallel computing.
Figure 3. Flow chart of the solution of an M × N Series-Parallel array using parallel computing.
Computation 10 00100 g003
Figure 4. I-V (left) and P-V (right) curves of a 15 × 2 array under two different partial shading conditions.
Figure 4. I-V (left) and P-V (right) curves of a 15 × 2 array under two different partial shading conditions.
Computation 10 00100 g004
Figure 5. Current errors in A and % for a 15 × 2 array under two different partial shading conditions. (Left): current errors for condition C 1 . (Right): current errors for condition C 2 .
Figure 5. Current errors in A and % for a 15 × 2 array under two different partial shading conditions. (Left): current errors for condition C 1 . (Right): current errors for condition C 2 .
Computation 10 00100 g005
Figure 6. Calculation time scalability analysis for an SP array with 15 rows ( N = 15 ) and columns varying from 1 to 20 ( M = [ 1 20 ] ).
Figure 6. Calculation time scalability analysis for an SP array with 15 rows ( N = 15 ) and columns varying from 1 to 20 ( M = [ 1 20 ] ).
Computation 10 00100 g006
Figure 7. Reconfiguration system used in the application example.
Figure 7. Reconfiguration system used in the application example.
Computation 10 00100 g007
Figure 8. I-V curves (left) and current errors in A and % (right) of a 3 × 3 array under partial shading conditions.
Figure 8. I-V curves (left) and current errors in A and % (right) of a 3 × 3 array under partial shading conditions.
Computation 10 00100 g008
Table 1. Average calculation time and standard deviation (std) for an SP array with 15 rows ( N = 15 ) and columns varying from 1 to 20 ( M = [ 1 20 ] ).
Table 1. Average calculation time and standard deviation (std) for an SP array with 15 rows ( N = 15 ) and columns varying from 1 to 20 ( M = [ 1 20 ] ).
Columns (M)Serial ModelParallel ModelAverage Time Reduction (%)Std (%)
Average Time (s)Std (s)Average Time (s)Std (s)
14.9501.4853.8700.91420.5606.943
26.5130.6324.6550.65028.03012.026
38.0150.0746.2570.43521.9235.417
49.9860.3246.2490.06737.3821.621
512.5560.9278.5000.24132.1083.425
614.2510.5298.6040.01339.5552.224
716.3890.60310.5300.03435.6812.250
817.8000.44110.7910.04539.3461.395
919.2040.80512.7180.04633.6802.753
1020.1690.06213.2610.03034.2490.191
1122.0100.25215.1080.02731.3530.804
1223.7080.06215.2280.04635.7660.300
1325.6860.06717.3230.04332.5590.252
1427.1640.13216.9420.11237.6290.554
1529.2890.85418.7950.08835.7941.530
1630.6670.04519.1490.02937.5580.145
1732.7170.06521.1930.03035.2230.163
1834.6190.07521.4630.06438.0020.182
1936.4270.13123.2970.03836.0440.260
2038.3010.10023.4440.04238.7890.246
Table 2. I p h currents of the reconfigurable submodules in the array (the first two rows).
Table 2. I p h currents of the reconfigurable submodules in the array (the first two rows).
Shading ProfileString
M = 1M = 2M = 3M = 4
11.17997.16625.99163.1579
7.26653.63401.06246.8555
24.57450.71313.55091.5638
5.08605.12870.54753.0708
33.58420.58993.58312.9958
3.61193.63331.81700.5312
Table 3. Average reconfiguration times and standard deviations (std) to evaluate 105 configuration for an SP array with 15 rows ( N = 15 ) and 4 columns ( M = 4 ).
Table 3. Average reconfiguration times and standard deviations (std) to evaluate 105 configuration for an SP array with 15 rows ( N = 15 ) and 4 columns ( M = 4 ).
ProfileSerial ModelParallel ModelAverage Time Reduction (%)Std (%)Power Increment (%)
Average Time (s)Std (s)Average Time (s)Std (s)
1733.0451.236517.957517.95729.3410.2297.390
2715.4550.654515.705515.70527.9190.0954.793
3711.3990.307522.080522.08026.6120.0914.363
Table 4. I p h currents of the experimental 3 × 3 array.
Table 4. I p h currents of the experimental 3 × 3 array.
RowString
M = 1M = 2M = 3
N = 13.72498.35646.2051
N = 23.72498.35646.1083
N = 33.72498.35646.1083
Table 5. Double diode model parameters of the 3 × 3 array.
Table 5. Double diode model parameters of the 3 × 3 array.
ParametersString
M = 1 and N = [1, 2, 3]M = 2 and N = [1, 2, 3]M = 3 and N = [1, 2, 3]
I s 1  ( η A)0.3633.8290.270
I s 2  (mA)0.1040.3580.089
N s · η 1 · V t d  (V)0.5480.5770.544
N s · η 2 · V t d  (V)1.0961.1541.089
R s  ( Ω )0.0110.0450.039
R h  (k Ω )6.00613.36310.015
I s b d  (mA)8.5158.5158.515
η b d · V t d  (mV)26.89284.91065.028
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bastidas-Rodríguez, J.D.; Ramos-Paja, C.A.; Serna-Garcés, S.I. Model of Series-Parallel Photovoltaic Arrays Using Double-Diode Model and Parallel Computing. Computation 2022, 10, 100. https://0-doi-org.brum.beds.ac.uk/10.3390/computation10060100

AMA Style

Bastidas-Rodríguez JD, Ramos-Paja CA, Serna-Garcés SI. Model of Series-Parallel Photovoltaic Arrays Using Double-Diode Model and Parallel Computing. Computation. 2022; 10(6):100. https://0-doi-org.brum.beds.ac.uk/10.3390/computation10060100

Chicago/Turabian Style

Bastidas-Rodríguez, Juan David, Carlos Andrés Ramos-Paja, and Sergio Ignacio Serna-Garcés. 2022. "Model of Series-Parallel Photovoltaic Arrays Using Double-Diode Model and Parallel Computing" Computation 10, no. 6: 100. https://0-doi-org.brum.beds.ac.uk/10.3390/computation10060100

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop