Optimal Power Flow Using Particle Swarm Optimization of Renewable Hybrid Distributed Generation

The problem of voltage collapse in power systems due to increased loads can be solved by adding renewable energy sources like wind and photovoltaic (PV) to some bus-bars. This option can reduce the cost of the generated energy and increase the system efficiency and reliability. In this paper, a modified smart technique using particle swarm optimization (PSO) has been introduced to select the hourly optimal load flow with renewable distributed generation (DG) integration under different operating conditions in the 30-bus IEEE system. Solar PV and wind power plants have been introduced to selected buses to evaluate theirs benefits as DG. Different solar radiation and wind speeds for the Dammam site in Saudi Arabia have been used as an example to study the feasibility of renewable energy integration and its effect on power system operation. Sensitivity analysis to the load and the other input data has been carried out to predict the sensitivity of the results to any deviation in the input data of the system. The obtained results from the proposed system prove that using of renewable energy sources as a DG reduces the generation and operation cost of the overall power system.


Introduction
Optimal power flow (OPF) studies aim to optimize specific objectives by adjusting some power system variables, provided that all equality and inequality constraints of the system are satisfied. Based on its energy management capability, OPF became one of the most important areas of study in the electric power field. OPF is a highly constrained and nonlinear problem with continuous and discrete variables. In 1962, Huneault et al. [1] published the first paper on OPF. Since then, researchers started studying the OPF problem extensively. The methods or techniques that have been implemented to optimize the OPF are divided into two categories: deterministic methods and evolutionary methods. Deterministic methods include linear [2] and nonlinear programming (LP) [3], quadratic programming (QP) [4], and interior point method (IM) [5]. These methods have a problem in handling many local minima due to the non-convexity of OPF problems. Gradient-based methods overcome the convergence problem, but sometime fail to meet inequality constraints [6]. Due to the limitations of deterministic methods, evolutionary methods were introduced to remedy these

Particles Swarm Optimization Technique
PSO is smart optimization technique that has been used in many applications. This was first introduced to scientific community by Eberhart and Kennedy [32]. As shown in this paper, this optimization technique is an adaptive algorithm based on a social-psychological metaphor; a population of individuals (referred to as particles) adapts by returning stochastically toward previously successful regions. In this technique a velocity and position should be updated continually. During each iteration, each particle is accelerated toward the particle's previous best position and the global best position. Under each iteration, the particle position will be updated using the new velocity which is counted as the distance from its previous best position, and the distance from the global best position will be calculated. The new velocity of each particle will be used to calculate the next position of the particle in the search space. This process will be repeated again for number of iterations until a minimum error is achieved [33].
The velocity and position of each particle during the PSO optimization steps are shown in Equations (1) and (2) [32]: v ij = ω.v ij + c 1 .r 1 .(p ij − x ij ) + c 2 .r 2 .(p g − x ij ) (1) Energies 2017, 10, 1013 3 of 14 (2) where, c 1 and c 2 are constants in the range from 0 to 4 and has been used as 2 in this paper, ω is the inertia, r 1 , r 2 are random variables which are uniformly distributed between 0 to 1, p g is the global best position, x ij is current position, p ij the best position for the current particle and v ij is the velocity of particles. j is the counter for iteration number, and i is the counter for particle number. The pseudo code of the procedure is as follows [34]: Begin; Generate random population of N solutions (particles); For each individual i N: calculate fitness (i); where N is the total number of particles. Initialize the value of the weight factor, ω; For each particle; Set p Best as the best position of particle i; If fitness (i) is better than p Best ; p Best = fitness (i); End; Set g Best as the best fitness of all particles; For each particle; Calculate particle velocity according to Equation (1); Update particle position according to Equation (2); End; Update the value of the weight factor, ω; Check if termination = true;

Optimal Power Flow Formulation
OPF is generally formulated mathematically using compact notation as follows [35]: with nonlinear equality constraints: G (x, u) = 0 (4) and inequality constraints: where, F(u, x) is the objective function which should be minimized (or maximized); G (x, u) represents the power flow nonlinear equations; H(u, x) represents various limits of the system control variables like generators power limits, bus voltages limits and lines limits; u is a vector of the system controlled variables like generators voltage, generators active power, transformers tap settings and shunt VAR compensations, load M and Mega Volt-Ampere Reactive MVAR (load shedding), dc power transmission; x is a vector of dependent variables that includes slack bus power, load bus voltages magnitude and phase angles, generator reactive power outputs for bus voltage control, reference bus voltage magnitude and angle (fixed) and transmission lines parameters.

Fuel Cost Objective Function
Total generation fuel cost function which is the objective function in this study is expressed as shown in the following equation [36]: where, γ, β and α are the cost function coefficients; and NG is number of generators. For a particular load the total generation cost should be minimum. The equality nonlinear constraints are given by the following equations: where, the active and reactive power generation at bus-bar i (P i and Q i ) can be obtained from the following equations: where, V i and V j are the buses voltages of busbars i and j, respectively. δ i , δ j are the angles of i th , and j th bus-bars, respectively, Y ij is the admittance between buses i and j, ϕ ij is the admittance angle of Y ij , and NB is number of buses. At any time, the total generation must meet the total load plus the transmission losses as shown in the following equation: where, ND is the number of load buses, NG is the number of generator buses, P Gi and P Di are the generation power and load power at bus-bar i, and P L is the total transmission losses. The above equations must be solved with the following generators constraints: where, P min and P max are the minimum and maximum active power constraints for each generator, respectively. Q min and Q max are the minimum and maximum reactive power constraints for each generator, respectively. V min and V max are the minimum and maximum voltage magnitude constraints for each bus. δ min and δ max are the minimum and maximum voltage angle constraints for each bus, respectively.

Wind and Solar Energy Modeling
The wind power generated by wind turbines is calculated using the following equation [37]: where, ρ is the air density in kg/m 3 , A is the rotor swept area in m 2 , C p is the turbine coefficient of performance (The theoretical value of C p is 59.3%), w is the wind speed at hub position m/s, N t is number of turbines and η is the net efficiency of wind turbine. However, the modern turbines reached a value of 50% [38]. In this paper the coefficient of performance has been used as 45% except in sensitivity analysis showing its effect on the total cost. The solar energy is calculated using the following equation: where, I is the solar radiation in W/m 2 , A p is the panel area, N p is the number of panels, and η p is the panels' efficiency.

System Description
This paper studies the OPF problem in the case of integration of renewable energy generations with the power system as DG. The study was applied on the IEEE 30-bus system as shown in Figure 1. A full description of this system's components can be found in [39]. Variable loads at the system buses were assumed. Four renewable energy power plants are inserted into the power system, mainly, two wind energy power plants and two solar energy power plants. Two wind energy power plants are inserted at buses 7 and 12 with 250 wind turbines at each wind farm. Two solar PV power plants are inserted at buses 8 and 18. The number of the PV panels is 40,000 and 380,000, respectively. The total capacity of the renewable power was chosen so that it represents 10% of the conventional generation to avoid technical problems [40,41]. The wind speed and solar radiation are assumed variable and different at each location of renewable generation power plant to be flexible when using this strategy in any real power system. where, I is the solar radiation in W/m 2 , Ap is the panel area, Np is the number of panels, and ηp is the panels' efficiency.

System Description
This paper studies the OPF problem in the case of integration of renewable energy generations with the power system as DG. The study was applied on the IEEE 30-bus system as shown in Figure 1. A full description of this system's components can be found in [39]. Variable loads at the system buses were assumed. Four renewable energy power plants are inserted into the power system, mainly, two wind energy power plants and two solar energy power plants. Two wind energy power plants are inserted at buses 7 and 12 with 250 wind turbines at each wind farm. Two solar PV power plants are inserted at buses 8 and 18. The number of the PV panels is 40,000 and 380,000, respectively. The total capacity of the renewable power was chosen so that it represents 10% of the conventional generation to avoid technical problems [40,41]. The wind speed and solar radiation are assumed variable and different at each location of renewable generation power plant to be flexible when using this strategy in any real power system. The system consists of six generators at buses 1, 2, 5, 8, 11 and 13. The Newton-Raphson method has been used in calculating the load flow equations of the power system at each time point. The PSO algorithm was implemented to find the optimal generation scheduling under variable load and renewable generation for certain period of time (24 h to represent a day has been used in this study).
The selection of these buses for integration of renewable energy generation is based on selecting a bus that is connected to other loaded buses. The idea behind this is to locate the renewable energy power plants near load centers to minimize the transmission losses. Both wind speed and solar radiation are assumed different at each bus to provide more flexibility to the proposed program. The wind speed, solar radiation, and the loads of the system buses are varied over a 24 h period to simulate the actual real time scenarios for one day. The data of wind speed and solar radiation were obtained from the Dammam site at Saudi Arabia as an example. This enables us to predict the dynamic performance of the power system in terms of economics and technical point of views. The operation and maintenance cost (O&M) of wind energy has been taken as 1.655 $/MWh for PV systems [42]. For wind the O&M cost is taken as 2.250 $/MWh [43]. The commercially available SM100 PV array (Siemens, Munich, Germany) has been considered in the design of 10 kW PV system with 15% efficiency. The electric characteristics of this PV module are rated power of 100 W, rated voltage of 34 V, rated current of 2.95 A, nominal temperature of 20 °C, and rated saturation current of 0.00405 A. The PV array consists of several modules in series (Ns) and parallel rows (Np).
To check the suitability of the selected buses for integration of renewable energy sources, the buses' voltage magnitude can be investigated. Figure 2 shows the buses' voltage magnitude for the base case without renewable sources and for the system with integration of renewable energy The system consists of six generators at buses 1, 2, 5, 8, 11 and 13. The Newton-Raphson method has been used in calculating the load flow equations of the power system at each time point. The PSO algorithm was implemented to find the optimal generation scheduling under variable load and renewable generation for certain period of time (24 h to represent a day has been used in this study).
The selection of these buses for integration of renewable energy generation is based on selecting a bus that is connected to other loaded buses. The idea behind this is to locate the renewable energy power plants near load centers to minimize the transmission losses. Both wind speed and solar radiation are assumed different at each bus to provide more flexibility to the proposed program. The wind speed, solar radiation, and the loads of the system buses are varied over a 24 h period to simulate the actual real time scenarios for one day. The data of wind speed and solar radiation were obtained from the Dammam site at Saudi Arabia as an example. This enables us to predict the dynamic performance of the power system in terms of economics and technical point of views. The operation and maintenance cost (O&M) of wind energy has been taken as 1.655 $/MWh for PV systems [42]. For wind the O&M cost is taken as 2.250 $/MWh [43]. The commercially available SM100 PV array (Siemens, Munich, Germany) has been considered in the design of 10 kW PV system with To check the suitability of the selected buses for integration of renewable energy sources, the buses' voltage magnitude can be investigated. Figure 2 shows the buses' voltage magnitude for the base case without renewable sources and for the system with integration of renewable energy sources. It's clear from this figure that the deviation of buses' voltage magnitude from 1 p.u. decreased after the integration of the renewable energy sources at buses 7, 8, 12 and 18. Therefore, the selection of the buses 7, 8, 12 and 18 for integration of renewable energy sources shows an improvement in voltage stability.

Validation
To validate this work, the results of the base case which is the base load without integrating renewable energy power plants will be compared with results reported in the literature. The total load of the original 30 bus IEEE system is 283.4 MW. The total generation cost at that load is 801.8 $/h according to our model. Table 1 [19] shows different values found in the literature. One notes that our result is quite consistent with those results which give a validation of the work done in this paper.

Results and Discussion
The total load of the system is shown in Figure 3 for 24 h. The peak load occurs at midday and the maximum load is 283.1 MW at hour 13. The renewable energy generation over the 24 h period is shown in Figure 4. The total generated power from solar panels, the total wind power generation, and the total renewable generation which is the sum of solar PV power and wind power generation are illustrated in

Validation
To validate this work, the results of the base case which is the base load without integrating renewable energy power plants will be compared with results reported in the literature. The total load of the original 30 bus IEEE system is 283.4 MW. The total generation cost at that load is 801.8 $/h according to our model. Table 1 [19] shows different values found in the literature. One notes that our result is quite consistent with those results which give a validation of the work done in this paper.

Results and Discussion
The total load of the system is shown in Figure 3 for 24 h. The peak load occurs at midday and the maximum load is 283.1 MW at hour 13.
The renewable energy generation over the 24 h period is shown in Figure 4. The total generated power from solar panels, the total wind power generation, and the total renewable generation which is the sum of solar PV power and wind power generation are illustrated in Figure 4. Maximum PV power occurs at hour 11. Usually, solar radiation is maximum at the hour 12, but sometimes cloud and dust reduce the radiation reaching the PV panels level at a particular time. This explains the shape of solar PV power generation curve shown in Figure 4. The wind power generation curve could have any shape according to the wind speed. It may be like constant with little or large fluctuations or it may be intermittent with very low values at some hours and high values at others. Here in this paper, the wind power generation has its peak amplitude between the hours 12 and 13. A good feature in the total renewable energy production is the correlation between the generated power from renewable energy sources with the electric loads as shown in Figures 3 and 4. This helps in extending the total power generation capacity of the system and meeting the extra demand without installing new conventional fuel power plants. The maximum total renewable power generation is 31.8 MW at time 12, which represents 10.3% of the conventional fuel power generation. according to our model. Table 1 [19] shows different values found in the literature. One notes that our result is quite consistent with those results which give a validation of the work done in this paper.

Results and Discussion
The total load of the system is shown in Figure 3 for 24 h. The peak load occurs at midday and the maximum load is 283.1 MW at hour 13. The renewable energy generation over the 24 h period is shown in Figure 4. The total generated power from solar panels, the total wind power generation, and the total renewable generation which is the sum of solar PV power and wind power generation are illustrated in Figure 4. Maximum PV power occurs at hour 11. Usually, solar radiation is maximum at the hour 12, but sometimes cloud and dust reduce the radiation reaching the PV panels level at a particular time. This explains the shape of solar PV power generation curve shown in Figure 4. The wind power generation curve could have any shape according to the wind speed. It may be like constant with little or large fluctuations or it may be intermittent with very low values at some hours and high values at others. Here in this paper, the wind power generation has its peak amplitude between the hours 12 and 13. A good feature in the total renewable energy production is the correlation between the generated power from renewable energy sources with the electric loads as shown in Figures 3 and 4. This helps in extending the total power generation capacity of the system and meeting the extra demand without installing new conventional fuel power plants. The maximum total renewable power generation is 31.8 MW at time 12, which represents 10.3% of the conventional fuel power generation.

Total Conventional Generation
The total conventional generation of the system with and without renewable energy generation integration is shown in Figure 5. This figure illustrates that the renewable energy system provides the power to the load at peak time which explains the correlation between the generation and the load from renewable energy system which presents its superiority to be used in the power system as a DG. This figure shows a significant reduction in conventional generation and hence, fossil fuel saved when renewable energy generation is integrated. The availability of renewable energy in the peak load period decreases the generators loading in these times. The peak of generators loading is shifted from midday hours, between hours 16-19, in this case. The availability of solar power during peak load period will definitely decrease of generators' loading, but wind speed may contribute at any time during the day. Therefore, different generators' loading curves or shapes will be produced. During the early and late hours, the two curves are close to each other. This is because of low contribution from renewable energy during these periods.

Total Renewable Generation,MW
Total conventional geneation, MW

Total Conventional Generation
The total conventional generation of the system with and without renewable energy generation integration is shown in Figure 5. This figure illustrates that the renewable energy system provides the power to the load at peak time which explains the correlation between the generation and the load from renewable energy system which presents its superiority to be used in the power system as a DG. This figure shows a significant reduction in conventional generation and hence, fossil fuel saved when renewable energy generation is integrated. The availability of renewable energy in the peak load period decreases the generators loading in these times. The peak of generators loading is shifted from midday hours, between hours 16-19, in this case. The availability of solar power during peak load period will definitely decrease of generators' loading, but wind speed may contribute at any time during the day. Therefore, different generators' loading curves or shapes will be produced. During the early and late hours, the two curves are close to each other. This is because of low contribution from renewable energy during these periods.
Furthermore, there are different incremental rates in conventional generation over the 24 h. During the hours between 1 and 10, the incremental rate is linear while it is closer to nonlinear between 10 and 19. This change in the increment rate is due to the simultaneous change in load demand and renewable energy generation. In the hours between 16 and 19, there is a slight decrease in conventional generation followed by obvious increase. These trends in few hours may occur in power systems integrated with renewable generation. However, this represents undesirable behavior in power systems and introduces big challenge in the operation and control of such power systems. load from renewable energy system which presents its superiority to be used in the power system as a DG. This figure shows a significant reduction in conventional generation and hence, fossil fuel saved when renewable energy generation is integrated. The availability of renewable energy in the peak load period decreases the generators loading in these times. The peak of generators loading is shifted from midday hours, between hours 16-19, in this case. The availability of solar power during peak load period will definitely decrease of generators' loading, but wind speed may contribute at any time during the day. Therefore, different generators' loading curves or shapes will be produced. During the early and late hours, the two curves are close to each other. This is because of low contribution from renewable energy during these periods. Further studies should focus on studying this phenomenon to overcome such a technical problem and others associated with integrating renewable energy generation.

Total Generation Cost
The cost of the total generation for the cases of conventional generation with and without adoption of renewable energy is shown in Figure 6. During the first and last day hours, the difference in the cost of the two cases is low due to the absence of PV output which makes a limited contribution from renewable energy during these periods. demand and renewable energy generation. In the hours between 16 and 19, there is a slight decrease in conventional generation followed by obvious increase. These trends in few hours may occur in power systems integrated with renewable generation. However, this represents undesirable behavior in power systems and introduces big challenge in the operation and control of such power systems. Further studies should focus on studying this phenomenon to overcome such a technical problem and others associated with integrating renewable energy generation.

Total Generation Cost
The cost of the total generation for the cases of conventional generation with and without adoption of renewable energy is shown in Figure 6. During the first and last day hours, the difference in the cost of the two cases is low due to the absence of PV output which makes a limited contribution from renewable energy during these periods. It can be observed that the cost of generated power is decreased significantly when renewable energy is integrated into the system. Before integration of renewable energy distribution generators, the optimal cost at highest load was 855.4 $/h. However, when renewable energy was integrated, the maximum cost decreased to 809.4 $/h (a 5.4% reduction in cost) under the same conditions. In addition, the occurrence of the maximum cost shifted from hour 13 to hour 16. This has the advantage of shifting the maximum generation cost to low energy cost periods in the case a load management strategy is applied. Moreover, the cost increase trend with load changed from a quadratic form (which is the case of conventional generation cost) to a piecewise linear curve. This means that, the rate of incremental fuel cost is reduced. It was concluded from Figure 6 that the cost of power generation is decreased significantly when renewable energy is integrated into the power system. Figure 7 shows the transmission power losses for both cases with and without renewable energy. The transmission losses are decreased when distributed renewable generation is integrated into the system. The loads that connected to the buses, where the renewable energy power plants are inserted, are supplied by the renewable energy power plants. Therefore, the power coming from the It can be observed that the cost of generated power is decreased significantly when renewable energy is integrated into the system. Before integration of renewable energy distribution generators, the optimal cost at highest load was 855.4 $/h. However, when renewable energy was integrated, the maximum cost decreased to 809.4 $/h (a 5.4% reduction in cost) under the same conditions. In addition, the occurrence of the maximum cost shifted from hour 13 to hour 16. This has the advantage of shifting the maximum generation cost to low energy cost periods in the case a load management strategy is applied. Moreover, the cost increase trend with load changed from a quadratic form (which is the case of conventional generation cost) to a piecewise linear curve. This means that, the rate of incremental fuel cost is reduced. It was concluded from Figure 6 that the cost of power generation is decreased significantly when renewable energy is integrated into the power system.  Figure 7 shows the transmission power losses for both cases with and without renewable energy. The transmission losses are decreased when distributed renewable generation is integrated into the system. The loads that connected to the buses, where the renewable energy power plants are inserted, are supplied by the renewable energy power plants. Therefore, the power coming from the conventional fuel power plants; far away from these buses; is decreased. Consequently, the power flow in the corresponding transmission lines was low and the transmission lines losses were reduced. If the generated power from renewable energy sources is more than the connected loads to the relevant buses, the extra power is transferred to the buses in the neighborhood depending on OPF methodology. It can be observed that transmission losses decreased from 10.14 MW; when the contribution of renewable energy sources is zero; to 9.15 MW with the adoption of renewable sources.

Sensitivity Analysis
This section is studying effect of some parameters on the power system performance. Since the renewable energy resources are variable and intermittent in nature, the power system performance will definitely be affected. Sensitivity analysis is additionally conducted for considering the error in measurement. Moreover, the load is changing in nature which makes the sensitivity analysis is vital to predict the results with any deviations in the input data.

Effect of Renewable Energy and Load Changes on Generation Cost
The effect of change in renewable energy and change in load demand on total generation cost is shown in Figure 8. A small change of ±5% in both renewable energy generation and load was assumed. This figure shows that the worst case occurs when the load is increased by 5% and the renewable energy contribution is zero. The total cost reached more than 900 $/h. On the other hand, the lowest cost was achieved when there is a renewable energy contribution while the load is decreased by 5%, the total cost becomes 730 $/h. It can be observed that transmission losses decreased from 10.14 MW; when the contribution of renewable energy sources is zero; to 9.15 MW with the adoption of renewable sources.

Sensitivity Analysis
This section is studying effect of some parameters on the power system performance. Since the renewable energy resources are variable and intermittent in nature, the power system performance will definitely be affected. Sensitivity analysis is additionally conducted for considering the error in measurement. Moreover, the load is changing in nature which makes the sensitivity analysis is vital to predict the results with any deviations in the input data.

Effect of Renewable Energy and Load Changes on Generation Cost
The effect of change in renewable energy and change in load demand on total generation cost is shown in Figure 8. A small change of ±5% in both renewable energy generation and load was assumed. This figure shows that the worst case occurs when the load is increased by 5% and the renewable energy contribution is zero. The total cost reached more than 900 $/h. On the other hand, the lowest cost was achieved when there is a renewable energy contribution while the load is decreased by 5%, the total cost becomes 730 $/h. Figure 9 shows the effect of renewable energy capacity variation and the effect of load change on transmission losses of the power system. The load and renewable energy generation are assumed to be changed in the range of ±5%. The solid curve with square marker represents the base case of base load without renewable energy generation contribution and the maximum losses exceeded 10 MW at hour 13. It's obvious that, the load change has a greater effect than the renewable energy generation change. When load is increased by 5% without using renewable energy sources, the transmission losses reached up to 11 MW, so with a 5% increase in load, one has a 10% increase in transmission line losses in the case the renewable energy sources are not adoped. However, when the load decreased by 5% with zero renewable energy contribution; the transmission losses were less than 9.36 MW. The integration of renewable energy sources reduced these losses to 9.95 MW when the load is increased by 5% (instead of 11 MW), so with a 5% increase in load, there will be more than a 10% reduction in transmission line losses in the case of adoption of the renewable energy sources. Also, the transmission line losses became 8.4 MW when the load is decreased by 5% instead of 9.36 MW, so there is about a 16% reduction in transmission line losses when the load is reduced by 5% when adopting the renewable energy sources. Furthermore, the integration of renewable energy generation transferred the peak value of these losses from hour 13 to hour 16. Therefore, the probability of transmission lines to be congested is shifted from the time when the conventional peak load occurs. renewable energy sources is zero; to 9.15 MW with the adoption of renewable sources.

Sensitivity Analysis
This section is studying effect of some parameters on the power system performance. Since the renewable energy resources are variable and intermittent in nature, the power system performance will definitely be affected. Sensitivity analysis is additionally conducted for considering the error in measurement. Moreover, the load is changing in nature which makes the sensitivity analysis is vital to predict the results with any deviations in the input data.

Effect of Renewable Energy and Load Changes on Generation Cost
The effect of change in renewable energy and change in load demand on total generation cost is shown in Figure 8. A small change of ±5% in both renewable energy generation and load was assumed. This figure shows that the worst case occurs when the load is increased by 5% and the renewable energy contribution is zero. The total cost reached more than 900 $/h. On the other hand, the lowest cost was achieved when there is a renewable energy contribution while the load is decreased by 5%, the total cost becomes 730 $/h.  Figure 9 shows the effect of renewable energy capacity variation and the effect of load change on transmission losses of the power system. The load and renewable energy generation are assumed to be changed in the range of ±5%. The solid curve with square marker represents the base case of base load without renewable energy generation contribution and the maximum losses exceeded 10 MW at hour 13. It's obvious that, the load change has a greater effect than the renewable energy generation change. When load is increased by 5% without using renewable energy sources, the transmission losses reached up to 11 MW, so with a 5% increase in load, one has a 10% increase in transmission line losses in the case the renewable energy sources are not adoped. However, when the load decreased by 5% with zero renewable energy contribution; the transmission losses were less than 9.36 MW. The integration of renewable energy sources reduced these losses to 9.95 MW when the load is increased by 5% (instead of 11 MW), so with a 5% increase in load, there will be more than a 10% reduction in transmission line losses in the case of adoption of the renewable energy sources. Also, the transmission line losses became 8.4 MW when the load is decreased by 5% instead of 9.36 MW, so there is about a 16% reduction in transmission line losses when the load is reduced by 5% when adopting the renewable energy sources. Furthermore, the integration of renewable energy generation transferred the peak value of these losses from hour 13 to hour 16. Therefore, the probability of transmission lines to be congested is shifted from the time when the conventional peak load occurs.

Effect of Changing Solar PV Efficiency on Total Cost
The efficiency of PV system has a great influence on the total cost of the PV system. The efficiency of PV system used in all above calculations was 15%. The sanctity analysis of PV system efficiency has been carried out by varying its value to 10%, 15% and 20% as shown in Figure 10. It is clear from this figure that the lowering the PV efficiency will increase the cost during the hours having PV contributions and vice versa. During the hours that don't have PV generation contribution the total cost will not change.

Effect of Changing Solar PV Efficiency on Total Cost
The efficiency of PV system has a great influence on the total cost of the PV system. The efficiency of PV system used in all above calculations was 15%. The sanctity analysis of PV system efficiency has been carried out by varying its value to 10%, 15% and 20% as shown in Figure 10. It is clear from this figure that the lowering the PV efficiency will increase the cost during the hours having PV contributions and vice versa. During the hours that don't have PV generation contribution the total cost will not change.
The efficiency of PV system has a great influence on the total cost of the PV system. The efficiency of PV system used in all above calculations was 15%. The sanctity analysis of PV system efficiency has been carried out by varying its value to 10%, 15% and 20% as shown in Figure 10. It is clear from this figure that the lowering the PV efficiency will increase the cost during the hours having PV contributions and vice versa. During the hours that don't have PV generation contribution the total cost will not change.

Effect of Changing Coefficient of Performance on Total Cost
As has been shown above for the sensitivity analysis on changing the PV system efficiency, the same study for changing the coefficient of performance of wind energy system has been performed. Three different values, 40%, 45% and 50%, have been used in this study. Due to the limited contribution the effect of changing the coefficient of performance on the cost does not show a considerable influence on the total cost of the power system as seen in Figure 11. It is clear from this figure that the cost is reduced by 1-2% when the coefficient of performance increases from 45% to 50%. The cost also increased by almost the same percentage upon reducing the coefficient of performance from 45% to 50%.

Effect of Changing Coefficient of Performance on Total Cost
As has been shown above for the sensitivity analysis on changing the PV system efficiency, the same study for changing the coefficient of performance of wind energy system has been performed. Three different values, 40%, 45% and 50%, have been used in this study. Due to the limited contribution the effect of changing the coefficient of performance on the cost does not show a considerable influence on the total cost of the power system as seen in Figure 11. It is clear from this figure that the cost is reduced by 1-2% when the coefficient of performance increases from 45% to 50%. The cost also increased by almost the same percentage upon reducing the coefficient of performance from 45% to 50%.

Conclusions
The integration of renewable energy in existing power systems in the form of DG was addressed in this study. A 30-bus IEEE power system was used as a case study to assess the advantages and disadvantages of integration of different renewable energy sources in the form of DG. Wind energy and solar energy were integrated into the proposed power system. The results showed an enhancement in the economical operation of fossil fuel generators. The total cost of operation was decreased from 855.4 $/h to 809.4 $/h during peak load periods. The total transmission losses were decreased by 10% in case of renewable energy source use at the selected buses. A sensitivity analysis was performed to investigate the performance of the system under varying uncertain data conditions. In this analysis, the effect of changing the loads by ±5% with and without renewable energy sources on the cost and transmission losses have been examined. The results from the sensitivity analysis

Conclusions
The integration of renewable energy in existing power systems in the form of DG was addressed in this study. A 30-bus IEEE power system was used as a case study to assess the advantages and disadvantages of integration of different renewable energy sources in the form of DG. Wind energy and solar energy were integrated into the proposed power system. The results showed an enhancement in the economical operation of fossil fuel generators. The total cost of operation was decreased from 855.4 $/h to 809.4 $/h during peak load periods. The total transmission losses were decreased by 10% in case of renewable energy source use at the selected buses. A sensitivity analysis was performed to investigate the performance of the system under varying uncertain data conditions. In this analysis, the effect of changing the loads by ±5% with and without renewable energy sources on the cost and transmission losses have been examined. The results from the sensitivity analysis show a 5.4% reduction in cost, 16% reduction in transmission line losses in the case of renewable energy source adoption, and a 5% reduction in load compared with the case of no renewable energy source adoption, which prove the superiority of the proposed analysis.