Research Article | Open Access
Saad Motahhir, Abdelaziz El Ghzizal, Souad Sebti, Aziz Derouich, "Modeling of Photovoltaic System with Modified Incremental Conductance Algorithm for Fast Changes of Irradiance", International Journal of Photoenergy, vol. 2018, Article ID 3286479, 13 pages, 2018. https://doi.org/10.1155/2018/3286479
Modeling of Photovoltaic System with Modified Incremental Conductance Algorithm for Fast Changes of Irradiance
The first objective of this work is to determine some of the performance parameters characterizing the behavior of a particular photovoltaic (PV) panels that are not normally provided in the manufacturers’ specifications. These provide the basis for developing a simple model for the electrical behavior of the PV panel. Next, using this model, the effects of varying solar irradiation, temperature, series and shunt resistances, and partial shading on the output of the PV panel are presented. In addition, the PV panel model is used to configure a large photovoltaic array. Next, a boost converter for the PV panel is designed. This converter is put between the panel and the load in order to control it by means of a maximum power point tracking (MPPT) controller. The MPPT used is based on incremental conductance (INC), and it is demonstrated here that this technique does not respond accurately when solar irradiation is increased. To investigate this, a modified incremental conductance technique is presented in this paper. It is shown that this system does respond accurately and reduces the steady-state oscillations when solar irradiation is increased. Finally, simulations of the conventional and modified algorithm are compared, and the results show that the modified algorithm provides an accurate response to a sudden increase in solar irradiation.
The energy generated by the PV systems depends on various parameters, either environmental as temperature and irradiance or internal parameters of the PV panel, namely, the series and shunt resistors [1, 2]. Thus, the load imposes its own characteristic on the output power . Therefore, in order to predict and analyze the effect of these parameters on the PV power, the model of the PV panel should be previously studied and achieved, and this model should be in accordance with the real comportment of the PV panel. Therefore, different models were proposed in the literature. In , a single diode model is used; in , a two-diode model is proposed to illustrate the influence of the recombination of carriers; and in , a model of three diodes is used to present the effects which are neglected by the two-diode model. However, the single-diode model is the most adopted due to its good simplicity and accuracy . Moreover, manufacturers of PV panels offer only some characteristics. But other characteristics required to model PV panel are missed in the datasheet, as the photocurrent, the diode saturation current, the series and shunt resistors, and the ideality factor . Hence, in [7, 9, 10], researchers have proposed different methods to extract the lacking characteristics based on the datasheet values, but these methods require an implementation and this can increase the time spent in the development of a PV application. Therefore, this paper aims firstly to extract lacking parameters in the manufacturers’ datasheet by using a simple tool provided by MathWorks  and then model the PV panel. The single-diode model is used in this work because it gives a high compromise between accuracy and simplicity , and several researchers have used it in their works [13, 14]. In addition, this work shows the effect of parameters that may change the performance of the PV panel. Then, this model will be a platform to design a PV array.
On the other side, the maximization of the PV power always remains a major challenge. Researchers have proposed different MPPT algorithms to maximize PV power, namely, FSCC, FOCV, fuzzy logic, neural network, P&O, and INC [15, 16]. FSCC and FOCV are the simplest MPPT algorithms, which are based on the linearity of short-circuit current or open-circuit voltage to the maximum power point current or voltage. However, these techniques isolate the PV panel to measure the short-circuit current or open-circuit voltage. Therefore, the loss of energy is increased due to the periodic isolation of the panel . Alternatively, fuzzy logic and neural network obtain a consistent MPPT technique due to their ability to treat the nonlinearity of the PV panel. But the main drawback of fuzzy logic is that the effectiveness of this algorithm depends a lot on choosing the right error computation and an appropriate rule base . In addition, neural network presents many disadvantages like the fact that the data needed for the training process have to be specifically acquired for every PV array and location, and also, the PV characteristics change with time, so the neural network has to be periodically trained . On the other hand, P&O and INC are mostly used. These techniques use the () characteristic of the PV panel. For P&O, steady-state oscillations occur after the MPP is found due to the perturbation made by this technique to maintain the MPP, which in turn increases the loss of power . For INC, it is founded in the fact that slope of characteristic is zero at the maximum power, and theoretically, there is no perturbation after the MPP is found. Therefore, oscillations are minimized. However, during implementation, the zero value is hardly found on the slope of the characteristic due to the truncation error in digital processing. Thus, the INC technique can make an inaccurate response when the irradiation is suddenly increased . Therefore, this work aims also to propose and implement a modified INC algorithm, which can overcome the wrong response made by the conventional INC algorithm when the irradiance is suddenly increased. Therefore, this work proposes a new technique to detect the increase in solar irradiation. The variation of voltage (∆V) and current (∆I) are used to identify the increase in irradiation instead of the slope (∆P/∆V) of the characteristic. The modified algorithm detects the increase of irradiance and makes a correct decision. Moreover, a mini error is accepted to admit that the slope is near to zero and minimize the steady-state oscillations.
This paper is structured as follows. Following the introduction, Section 2 presents the modeling of PV panel and array and presents the impact of different environmental and internal parameters. Section 3 presents the design of the boost converter, the conventional algorithm, and the proposed algorithm.
2. Modeling of PV Panel and Array
2.1. Model of PV Panel
As shown in Figure 1, the single-diode model of PV panel can be presented by a photocurrent source and a diode connected with series and shunt resistances. The mathematical model of the PV panel can be presented by the following equations :
Hence, the physical behavior of the PV panel depends on the shunt and series resistances, solar irradiation, and temperature. Therefore, in this work, the impact of these parameters on the output of the PV panel is investigated.
The panel used in this work is MSX-60 panel, and as presented in Table 1, the datasheet of PV panel provides only some characteristics of PV panel. Thus, other characteristics required to model PV panel are missed in the datasheet, as the photocurrent, the diode saturation current, the series and shunt resistors, and the ideality factor. Hence, researchers have proposed different methods to extract the lacking characteristics based on the datasheet [7, 9, 10]. But these methods require an implementation and this can increase the time spent in the development of a PV system. Therefore, this work aims firstly to extract these parameters by using a simple tool provided by MathWorks which is “PV array;” the latter is available in Simulink 2015 or later . Hence, as shown in Figure 2, we only have to set the datasheet values and automatically it will generate the lacking parameters.
Figure 4 shows the and curves of experimental and PSIM model under STC. The experimental data P(V) and I(V) are taken from the manufactured datasheet . And as presented, the model data are in accordance with the experimental data both in the current and power curves.
2.2. Effect of Solar Irradiation Variation
Figure 3 contains the model of the three equations: one of these equations computes the photocurrent based on temperature and irradiance (2). The model of this equation is presented in Figure 5, and Figure 6 presents and curves for different values of solar irradiation.
As presented in Figure 6, the PV panel current depends heavily on solar irradiation. However, the voltage increases just by 1 V once the irradiance is increased from 400 W/m2 to 1000 W/m2. Therefore, the irradiance change affects heavily the PV panel current.
2.3. Effect of the Temperature Variation
Figure 3 contains also the modeling of (3), which computes the diode saturation current based on the temperature. The model of this equation is shown in Figure 7, and Figure 8 shows the and curves for different values of temperature.
Generally, as shown in Figure 8, for a fixed solar irradiation and when the temperature increases, the open-circuit voltage decreases and the short-circuit current increases with a little value. Therefore, the temperature change affects strongly the PV panel voltage.
2.4. Effect of Series Resistor Variation
The series resistor value is very small, and it may be neglected in some cases. Nevertheless, to make the appropriate model for any PV panel, it is recommended to make a variation of this resistor and show its effect on the PV panel output. As shown in Figure 9, the change of the series resistor results on the deviation of the MPP.
The simulation was made for three values of series resistance (1 mΩ, 4 mΩ, and 8 mΩ). Moreover, as shown in Figure 9, the upper values of series resistance decrease the output power. In addition, the fill factor presented by (4) decreases as series resistance increases .
2.5. Effect of Shunt Resistor Variation
As presented in Figure 10, Rsh should be quite large for a good fill factor. In fact, when Rsh is small, the current collapses more strongly, then the loss of power is high, and the fill factor is low. Therefore, Rsh of any PV panel should be large enough for a good efficiency.
2.6. Effect of Shading
Partial shading also presents a major impact on PV output power. When the insolation received by a part of the PV panel (shaded cells) is less than the insolation received by another part (illuminated cells), the current generated by the illuminated cells is greater than the current produced by the shaded cells; this mismatch makes the diode of shaded cells reverse biased; consequently, the power will be lost in the shaded cells and that may cause a hot spot problem which is the reason of permanent damage to the PV panel . Hence, in order to overcome this problem, the bypass diodes can be connected in parallel with PV cells .
To simulate the effect of shading, a bypass diode is associated with each string of the panel, and it should be mentioned that the panel used includes two strings and each string is a set of 18 cells. Thus, as shown in Figure 11, the first string is exposed by 1000 W/m2 and the second string by 700 W/m2.
Under uniform irradiation, the bypass diodes have no impact because they are reversely biased. But under shading, the current flows through the bypass diode instead of the shaded string because the bypass diode is directly biased, and as a result, no power will be lost in the shaded cells and only the illuminated cells generate power. Figure 12 shows the effect of bypass diodes on the PV panel characteristics and as presented, multiple peaks may occur on the curve; as in this case, there are two peaks, point A which is the global peak and point B which is the local peak. Therefore, conventional MPPT algorithms are unable to track the global peak which is the real MPP .
2.7. Photovoltaic Array
To get benefit from the model developed, a PV array of 18 PV panels has been built in order to supply a solar pumping station, not studied in this paper. Therefore, as shown in Figure 13, three strings of six PV panels have been linked in parallel and each group is composed of six panels in series.
The model of the PV array has been achieved on PSIM, and simulation result obtained is presented in Figure 14.
As presented in Figure 14, the panels connected in parallel increase the current and the panels connected in series increase the voltage. However, as discussed in the previous section, this connection between panels can lead to hot spot problem when the insolation received by a part of the PV array (shaded panels) is less than the insolation received by another part (illuminated panels).
3. Modified INC Algorithm with Boost Converter
PV panel provides and curves presented in Figure 15; these curves highlight one point where the power is maximum. As presented above, this point depends on solar irradiation and temperature. Moreover, as presented in Figure 15, in general, load’s characteristic is different from the MPP. Therefore, the boost converter controlled by a duty cycle (α) generated by the MPPT controller is put between the panel and the load . The interest of this addition is to remove the mismatch between the panel and the load, and then the PV panel can operate at MPP.
3.1. Boost Converter Design
The operation principle of this converter is described by the following equations :
By using (5), (6) is obtained, which is the relationship between the resistance seen by the PV panel (Req) and the load resistance (R). Hence, based on this equation, the MPPT controller can find the optimum α to remove the mismatch between the load and MPP. Therefore, The boost converter is required to get maximum power available from the panel.
3.1.1. Selection of the Inductor
The choice of the inductor can directly influence the performance of the boost converter. Moreover, the selection of the inductance is a trade-off between its cost, its size, and the inductor current ripple. A higher inductance value results in a minor inductance current ripple; however, that results in a higher cost and larger inductor’s size, which means a larger PCB surface.
By the way, the inductance value can be given as follows:
During TON state, where ∆IL can be computed as below, and r is the inductor current ripple ratio, which is optimal in the range [0.3, 0.5] :
Therefore, the optimum inductor value can be computed by using
Based on Figure 17, in order to guarantee the performance of boost converter in the continuous conduction mode, the following equation must be verified:
3.1.2. Selection of the Output Capacitor
The choice of the output capacitor is made by using output voltage ripple as follows:
Therefore, the output capacitor value can be calculated as below, where the desired ∆VO equals to 2% of output voltage :
3.1.3. Selection of the Input Capacitor
An input capacitor is used to decrease the input voltage ripple and to deliver an alternative current to the inductor. The input voltage ripple matches to the charge voltage during the charge phase of the capacitor, and during this phase, ICin is greater than zero, so this phase is illustrated by the blue area in Figure 18; therefore, this area is used to calculate the input capacitor as follows:
The design of the used boost is presented in Table 2.
3.2. Problem with the Conventional INC Algorithm
A good MPPT algorithm balances between the tracking speed and steady-state performance. In accordance with these requirements, the INC algorithm can be used even if it can fail in some cases  and in this study, it will be modified in order to improve its performance. INC algorithm is founded in the fact that slope of characteristic is zero at the MPP . Therefore, this algorithm can be modeled as follows:
The flowchart of the INC algorithm is presented in Figure 19 . This algorithm measures the current and voltage of the panel. If (21) is met, the duty cycle is increased, and vice versa if (20) is met. Then, there is nothing to do if (19) is met. Therefore, theoretically, if MPP is reached, there is no more perturbation of α; consequently, steady-state oscillations are decreased, and that is the main advantage of INC algorithm.
However, the conventional INC algorithm fails to make a good decision when the irradiance is suddenly increased . As presented in Figure 20, once the solar irradiance is at 500 W/m2 and the PV system operates at load_2, the INC technique controls the PV system in order to reach the MPP (point B). When the irradiance is increased to 1000 W/m2, load_2 will lead the system to point G in characteristic, which matches to point C in characteristic. The INC technique calculates the slope between point C and point B which is positive. Therefore, the INC algorithm will decrease the duty cycle and consequently, the PV panel voltage will be increased. But since the MPP of 1000 W/m2 is at point A, and the slope between point A and C is negative, then the PV panel voltage should have been decreased in order to reach point A, instead of increase voltage and recede from point A as made by the conventional INC algorithm. In addition, as presented in Figure 6, generally when the solar irradiance increases, the MPP moves to the right and consequently, the same problem will occur.
Conversely, this weakness does not happen if the solar irradiation is decreased. Because as shown in Figure 20, the slope is positive between point A and D, and also between point B and D.
3.3. Modified INC Algorithm
Based on the above analysis, it is noted that when the solar irradiance increases, both the voltage and the current are increased. Therefore, the sudden increase in solar irradiation can be detected, by checking if the MPP was reached and both the voltage and current are increased. Therefore, a permitted error is accepted (22) to detect that the MPP is reached.
The proposed algorithm is presented in Figure 21. So as shown, the addition is the check if the MPP was reached by using (22), then set Var to one. After that, when (22) is not met and Var is one, the proposed technique checks if both voltage and current are increased; in this case, the duty cycle is increased instead of decreased as made by the conventional algorithm. Hence, the INC algorithm is modified to overcome the incorrect decision made by the conventional algorithm when the irradiance is increased.
4. Results and Discussions
The test was made for the conventional and the proposed techniques. At first, the solar irradiance is suddenly increased from 500 W/m2 to 1000 W/m2 at t = 0.11 s, and then it is decreased from 1000 W/m2 to 500 W/m2 at t = 0.23 s. Figure 22 shows the test result of the conventional technique, and Figure 23 shows the test result of the proposed technique. Therefore, as shown in these figures, the steady-state oscillations are minimized by using the proposed algorithm and admit an error equals to 0.07. Contrary to the conventional algorithm, the power oscillates between (28.5–29.8 W), and this can generate a loss in PV energy.
In addition, as presented in Figure 22, by using the conventional algorithm, when solar irradiation is suddenly increased, the power diverges with a value greater than 62 W and after that, the conventional INC reverses the direction and the power diverges with a value lower than 56 W; and like that, the system takes a long time to converge around MPP that is due to the wrong decision made by the conventional technique. Also, even if the MPP is reached, the power oscillates between (59–61 W). On the other side, as presented in Figure 23, the proposed technique detects the fast increase of irradiance and makes a correct decision in duty cycle. As a result, the power converges to the new MPP from the first step and it is maintained at it (60 W). In addition, it only needs 0.001 s to reach the MPP. Hence, by using the proposed INC, the power converges faster compared with the conventional algorithm which needs more time to reach the MPP.
Table 3 summarizes a comparison of the proposed technique to other improved INC techniques proposed in the scientific literature in term of the oscillation level, tracking efficiency, the response time during sudden increase in irradiation, and if the technique makes an incorrect decision under sudden increase of irradiation. As presented, the proposed technique shows a very fast tracking speed, a higher efficiency, and neglected oscillations around the MPP compared to other techniques. Thus, only the proposed algorithm and that proposed in  make a correct decision under sudden increase of irradiation, contrary to the conventional technique and those proposed in [30, 31] which make an incorrect decision.
In this paper, PV panel’s parameters are found using MathWorks tool (PV array); hence, by using these parameters, a PV panel and a PV array are modeled, and the results show that the model is in accordance with experimental data of the used panel (MSX-60). In addition, a modified INC algorithm which can overcome the confusion faced by the conventional INC technique is proposed in this paper. As a result, the tests show that the modified technique detects the fast increase of irradiation and makes a correct decision, contrary to the conventional technique. Moreover, by using the modified algorithm, steady-state oscillations are almost neglected. Hence, the loss of energy is minimized; consequently, the efficiency is equal to 98.8% instead of 96% obtained by the conventional technique.
As a perspective, the modified INC algorithm can be more improved and then implemented in an embedded hardware device.
|a:||Diode’s ideality factor|
|I:||Output current of the panel (A)|
|Is:||Diode saturation current (A)|
|Iph:||Panel photocurrent (A)|
|G:||Solar irradiation (W/m2)|
|K:||Boltzmann constant (J⋅K−1)|
|q:||Electron charge (C)|
|R:||The load (Ω)|
|Req:||The resistance seen by the panel (Ω)|
|Rs:||Series resistance (Ω)|
|Rsh:||Shunt resistance (Ω)|
|T:||Junction temperature (K)|
|V:||Output voltage of the panel (V)|
|VO:||Output voltage of the boost converter (V)|
|IO:||Output current of the boost converter (A)|
|F:||Switching frequency (Hz)|
|∆V:||Input voltage ripple of the boost converter (V)|
|∆VO:||Output voltage ripple of the boost converter (V)|
|∆IL:||Inductor current ripple (A).|
|CCM:||Continuous conduction mode|
|FSCC:||Fractional short-circuit current|
|FOCV:||Fractional open-circuit voltage|
|MPP:||Maximum power point|
|MPPT:||Maximum power point tracking|
|P&O:||Perturb and observe|
|STC:||Standard test conditions.|
Conflicts of Interest
The authors declare that they have no conflicts of interest.
- H. Patel and V. Agarwal, “MATLAB-based modeling to study the effects of partial shading on PV array characteristics,” IEEE Transactions on Energy Conversion, vol. 23, no. 1, pp. 302–310, 2008.
- S. Motahhir, A. Chalh, A. Ghzizal, S. Sebti, and A. Derouich, “Modeling of photovoltaic panel by using proteus,” Journal of Engineering Science and Technology Review, vol. 10, no. 2, pp. 8–13, 2017.
- S. Motahhir, A. El Ghzizal, S. Sebti, and A. Derouich, “Proposal and implementation of a novel perturb and observe algorithm using embedded software,” in 2015 3rd International Renewable and Sustainable Energy Conference (IRSEC), pp. 1–5, Marrakech, Morocco, December 2015.
- H. S. Rauschenbach, Solar Cell Array Design Handbook, Van Nostrand Reinhold, 1980.
- N. Barth, R. Jovanovic, S. Ahzi, and M. A. Khaleel, “PV panel single and double diode models: optimization of the parameters and temperature dependence,” Solar Energy Materials and Solar Cells, vol. 148, pp. 87–98, 2016.
- K. Nishioka, N. Sakitani, Y. Uraoka, and T. Fuyuki, “Analysis of multicrystalline silicon solar cells by modified 3-diode equivalent circuit model taking leakage current through periphery into consideration,” Solar Energy Materials and Solar Cells, vol. 91, no. 13, pp. 1222–1227, 2007.
- N. Yıldıran and E. Tacer, “Identification of photovoltaic cell single diode discrete model parameters based on datasheet values,” Solar Energy, vol. 127, pp. 175–183, 2016.
- K. Ishaque, Z. Salam, S. Mekhilef, and A. Shamsudin, “Parameter extraction of solar photovoltaic modules using penalty-based differential evolution,” Applied Energy, vol. 99, pp. 297–308, 2012.
- K. Ishaque and Z. Salam, “An improved modeling method to determine the model parameters of photovoltaic (PV) modules using differential evolution (DE),” Solar Energy, vol. 85, no. 9, pp. 2349–2359, 2011.
- M. F. AlHajri, K. M. El-Naggar, M. R. AlRashidi, and A. K. Al-Othman, “Optimal extraction of solar cell parameters using pattern search,” Renewable Energy, vol. 44, pp. 238–245, 2012.
- Mathworks, PV Array, 2015, Mai 2017, https://fr.mathworks.com/help/physmod/sps/powersys/ref/pvarray.html.
- C. Carrero, J. Amador, and S. Arnaltes, “A single procedure for helping PV designers to select silicon PV modules and evaluate the loss resistances,” Renewable Energy, vol. 32, no. 15, pp. 2579–2589, 2007.
- T. Radjai, L. Rahmani, S. Mekhilef, and J. P. Gaubert, “Implementation of a modified incremental conductance MPPT algorithm with direct control based on a fuzzy duty cycle change estimator using dSPACE,” Solar Energy, vol. 110, pp. 325–337, 2014.
- J. Ahmed and Z. Salam, “A modified P&O maximum power point tracking method with reduced steady-state oscillation and improved tracking efficiency,” IEEE Transactions on Sustainable Energy, vol. 7, no. 4, pp. 1506–1515, 2016.
- A. Amir, J. Selvaraj, and N. A. Rahim, “Study of the MPP tracking algorithms: focusing the numerical method techniques,” Renewable and Sustainable Energy Reviews, vol. 62, pp. 350–371, 2016.
- A. Gupta, Y. K. Chauhan, and R. K. Pachauri, “A comparative investigation of maximum power point tracking methods for solar PV system,” Solar Energy, vol. 136, pp. 236–253, 2016.
- D. Verma, S. Nema, A. M. Shandilya, and S. K. Dash, “Maximum power point tracking (MPPT) techniques: recapitulation in solar photovoltaic systems,” Renewable and Sustainable Energy Reviews, vol. 54, pp. 1018–1034, 2016.
- J. P. Ram, T. S. Babu, and N. Rajasekar, “A comprehensive review on solar PV maximum power point tracking techniques,” Renewable and Sustainable Energy Reviews, vol. 67, pp. 826–847, 2017.
- S. Motahhir, A. El Ghzizal, S. Sebti, and A. Derouich, “Shading effect to energy withdrawn from the photovoltaic panel and implementation of DMPPT using C language,” International Review of Automatic Control, vol. 9, no. 2, pp. 88–94, 2016.
- K. S. Tey and S. Mekhilef, “Modified incremental conductance MPPT algorithm to mitigate inaccurate responses under fast-changing solar irradiation level,” Solar Energy, vol. 101, pp. 333–342, 2014.
- Solarex MSX60 and MSX64 photovoltaic panel, datasheet, 1998, April 2017, https://www.solarelectricsupply.com/media/custom/upload/Solarex-MSX64.pdf.
- F. Bayrak, G. Ertürk, and H. F. Oztop, “Effects of partial shading on energy and exergy efficiencies for photovoltaic panels,” Journal of Cleaner Production, vol. 164, pp. 58–69, 2017.
- S. Motahhir, A. El Ghzizal, S. Sebti, and A. Derouich, “MIL and SIL and PIL tests for MPPT algorithm,” Cogent Engineering, vol. 4, p. 1378475, 2017.
- S. Sivakumar, M. J. Sathik, P. S. Manoj, and G. Sundararajan, “An assessment on performance of DC–DC converters for renewable energy applications,” Renewable and Sustainable Energy Reviews, vol. 58, pp. 1475–1485, 2016.
- N. Mohan and T. M. Undeland, Power Electronics: Converters, Applications, and Design, John Wiley & Sons, 2007.
- S. Maniktala, Switching Power Supplies A-Z, Elsevier, 2012.
- M. H. Rashid, Power Electronics: Circuits, Devices, and Applications, Pearson Education India, 2009.
- A. Loukriz, M. Haddadi, and S. Messalti, “Simulation and experimental design of a new advanced variable step size incremental conductance MPPT algorithm for PV systems,” ISA Transactions, vol. 62, pp. 30–38, 2016.
- A. Belkaid, I. Colak, and O. Isik, “Photovoltaic maximum power point tracking under fast varying of solar radiation,” Applied Energy, vol. 179, pp. 523–530, 2016.
- M. A. G. De Brito, L. Galotto, L. P. Sampaio, G. de Azevedo e Melo, and C. A. Canesin, “Evaluation of the main MPPT techniques for photovoltaic applications,” IEEE Transactions on Industrial Electronics, vol. 60, no. 3, pp. 1156–1167, 2013.
- P. C. Sekhar and S. Mishra, “Takagi–Sugeno fuzzy-based incremental conductance algorithm for maximum power point tracking of a photovoltaic generating system,” IET Renewable Power Generation, vol. 8, no. 8, pp. 900–914, 2014.
Copyright © 2018 Saad Motahhir et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.