INTRODUCTION
Shunt capacitors installed in radial distribution systems plays a vital role to reduce power loss, to increase available capacity of the feeder and for feeder voltage profile improvement. Harmony Search Algorithm (HSA) optimization technique which is heuristic in nature based on analogies with natural phenomena is integrated with harmonic power flow algorithm and Backward/Forward sweep for optimum shunt capacitor location and perfect capacitor sizing is proposed (Muthukumar and Jayalalitha, 2012). Hybrid Particle Swarm Optimization (HPSO) technique integrated with harmonic power flow technique for optimum capacitor allocation and proper sizing to bring down the total real power loss subjected to operating constraint and power quality constraints has been formulated by Eajal and ElHawary (2010). Positioning of capacitor at proper locations using Harmony search algorithm in balanced radial distribution system has been suggested and validated on the IEEE 9 and IEEE 34 bus Radial Distribution System (Sirjani et al., 2010).
Network reconfiguration and capacitor sizing problems are solved simultaneously using Harmony Search Algorithm (HSA) has been proposed by Srinivasa Rao (2010) and the placement of shunt capacitor at proper location is identified based on voltage stability index values. A three phase load flow method was formulated with mathematical modeling to include the mutual coupling between the phases (Subrahmanyam, 2009). A technique for three phase load flow suitable for radial distribution networks based on networktopology has been discussed by Teng (2000).
Evolutionary Programming (EP) approach is proposed in this study to determine the proper sizing of shunt capacitor in Unbalanced Radial Distribution System (URDS) to compensate reactive power demand at appropriate location based on computation of voltage stability index value.
MATHEMATICAL FORMULATION OF OPTIMIZATION PROBLEM
The choice of optimum capacitor size and appropriate location is framed as a problem of nonlinear integer optimization.
• 
objective function: The objective function is mathematically formulated such as the minimization of cost function related to real power loss along with investment cost of shunt capacitors as shown below (Eajal and ElHawary, 2010): 
Where:
K_{p}: 
The yearly cost ($) of the RPL kW^{1} 
K_{ci}: 
Yearly cost in ($) kVAR^{1} of the reactive power supplied at bus i by the capacitors 
Q_{ci}: 
Injection of reactive power (kVAR) at bus i by the capacitor 
nc: 
Number of shunt capacitors considered for installation 
Ploss: 
The total amount of active power loss 
Finally the sum of active power loss is attained by adding all the branch power losses as:
where, nb is the count of total number of branches.
The real power loss is estimated by radial distribution power flow (RDPF) (Subrahmanyam, 2009) algorithm based on Backward/ Forward sweep technique. Size of shunt capacitors chosen for installation has impact on the cost of reactive power injection.
To avoid solution diverge from constraint limits, it is necessary to add quadratic penalty factors with the fitness function to get a modified fitness function F _{c} as,
Where:
V: 
Bus voltage obtained from backward/forward sweep based load flow after the installation of capacitors 
λ_{1,} λ_{2,} λ_{3}: 
Penalty factors corresponding to minimum voltage, maximum voltage and total size of capacitor placed which can be adjusted during the optimization procedure 
For reactive power compensation, shunt capacitors can be placed at buses except substation bus since it act as slack bus. The different sizes of capacitor that can be put at these locations are considered as the control variables.
•  Constraints: Equality constraints are closely related with nonlinear load flow equation. The vector representation is as: 
Where:
x: 
State variable. 
y: 
Control variable 
Inequality constraints includes the rms values of bus voltage magnitudes and the total count of shunt capacitors considered for installation.
Throughout the optimization process, the rms value of the bus voltage magnitudes are maintained within specified operating limits.The rms value voltage magnitudes of each nodes should be between the lower and upper limits.
where, n is the total count of buses.
Number of capacitors to be placed and the available standard capacitor sizes (discrete) are also considered as the constraints:
Where:
Q_{0}: 
Least capacitor size (commercial available) 
Qci: 
The total injection of reactive power 
nc: 
The total number of shunt capacitors to be considered for installation 
The injection of reactive power from the shunt capacitors installed in distribution network should be kept lower than the reactive demand of distribution network:
where, Q_{T} is the total reactive power demand.
EVOLUTIONARY PROGRAMMING
Evolutionary Programming is one of the AI based search algorithms works on replicated procedure of natural selections, mutation and selection. EP uses population of potential solutions which can move over the search space to discover a global optimum solution. It is initialized with random number generation which represents the parameters responsible for the optimization of the fitness value. EP involves statistics, calculation of fitness value, mutation and finally the selection process to get new population (Lee and Yang, 1998).
Step 1: 
Initialization: Evolutionary Programming works with a random set of initial data (capacitor sizes) which forms the parent population. The size of capacitor connected to any phase of a bus is limited by the upper bound (Qmax) of the systems total reactive power demand and the minimum size of commercially available capacitor is considered as the lower bound (Qmin). The control variable (size of capacitors) is initialized by randomly selecting “m” values ranging over [Qmin, Qmax], where “m” is the population size. The generated populations are considered as the parent population and corresponding fitness functions (Fc) are evaluated by running Backward/ Forward sweep based loadflow for every individuals of parent population 
Step 2: 
Fitness calculation and statistics: The fitness function considered here is cost of real power loss which needs to be minimised. In this study, the fitness is calculated by using Backward/Forward sweep load flow algorithm. Subsequently, the maximum value of fitness, minimum value of fitness, sum and average values of fitnesses will be estimated and it will be utilized during the mutation process 
Step 3: 
Mutation: This process is carried out based on the random numbers (populations) to produce siblings (child). The procedure is implemented by using Eq. 9 and 10. Component k of individual i is represented as (i, k): 
μ = 0
Where:
C_{i, k}: 
In child population 
P_{i, k }: 
In parent population 
N (μ, σ^{2}): 
Gaussian random variables with mean μ=0, variance σ2 
f_{i }: 
Fitness of i^{th }individual in parent population 
f_{maximum}: 
Maximum fitness value for parent population 
x_{k, maximum}: 
Maximum limit of jth element 
x_{k, minimu}: 
Minimum limit of jth element 
β: 
Mutation scale, 0< β≤1 
n: 
Total number of control variables 
In Evolutionary Programming, a small fixed value of mutation factor can lead to premature convergence and thereby getting trapped into local optimum while searching with large value of mutation factor, prevents EP from converging. By using an adaptive mutation scale, this problem can be solved by choosing mutation scale as follows:
where, “k” is generation number.
β_{init}, β_{final}, β_{step} are fixed values and its values would be β_{init} = 1, β_{final} = 0.005, β_{step} is between 0.001 to 0.01, depends on the maximum generation number (Ma and Lai, 1996). As the search process continues, the mutation factor will decrease. Mutation scale decrease faster, when fitness value is low. Such type of adaptive mutation factor prevents premature convergence.
Step 4: 
Selection: The offspring’s generated as the result of mutations are pooled along with their parents for selection process to find out the candidates that have the probability to be moved into the next new generation. Initially parent populations (size of capacitors) and its corresponding fitness values (cost of real power loss) will be estimated. After mutation, offsprings are created and merged with parent populations and arranged in the ascending order. The first half of the population ranked in the ascending order will be considered as the new population and choosen individuals are considered as the the new parent population for next generation process. The mutation and selection process is carried out repetitively until the convergence criterion is satisfied 
Step 5: 
convergence test: Finally the convergence condition is achieved, when the variation between maximum and minimum fitness value is less than 0.0001 which is the resolution required. The flow chart as illustrated in Fig. 1 represents the overall EP procedure 

Fig. 1:  Flowchart of proposed EP algorithm 
METHODS
The optimization problem is formulated to identify the best possible capacitor ratings for an unbalanced RDS using Evolutionary Programming to achieve minimum value of real power losses with improved voltage profiles. The best location for placing shunt capacitor is identified by computi the voltage stability index values. Load flow based on backward/forward sweep technique is utilized (Subrahmanyam, 2009) to obtain active and reactive power flows, bus voltage magnitudes, phase angles, line currents, active and reactive power losses.
The suggested approach is examined on modified IEEE 13 bus unbalanced RDS network without voltage regulator (Kersting, 1991) to assess the performance of the algorithms.

Fig. 2:  One line diagram of two bus distribution network 
If a distribution network with “P” number of buses and “Q” number of available size of capacitors, then the possible combination of solutions are (Q+1)^{P}. To minimize the computational burden and size of the solution vector, it is necessary to decide the most receptive nodes that require reactive power compensation which can be identified by computation of voltage stability index.
VOLTAGE STABILITY INDEX
In RDS for stable operation, the stability index value should be near to unity. Minimum value of VSI is the indication of buses in distribution systems that are prone to voltage instability and voltage collapse (Chakravorty and Das, 2001). Such nodes in radial distribution system are the exact places for placement of shunt capacitors.
Derivation f or Voltage Stability Index (VSI).
Computing voltage stability index in all the buses as follows.
From Fig. 2:
Where:
i, j: 
Transmitting and receiving end 
P_{2}, Q_{2}: 
Total Real and Reactive Power load connected in bus j 
V_{1}, V_{2}: 
Voltage magnitude at bus i and j 
Using the above equations, Voltage stability index of bus “j” is computed as:
Table 1:  Voltage stability index values before capacitive compensation 

Table 2:  Standard capacitor sizes with associated costs/KVAR (KiloVolt Ampere Reactive) 

For radial distribution system to be in stable conditions, VSI (j) = 0, for j = 2, 3...n.
In this proposed method, branch currents and bus voltages, active and reactive power flow of all the lines are estimated using the backward/forward sweep based load flow technique and from Eq. 13, the VSI of each bus is computed and tabulated in Table 1. From the tabulated results, for the proposed test system the bus number 671, 3rd phase is chosen for the perfect location for placing shunt capacitor which is having lower voltage stability index of 0.7032. To generate parent population in Evolutionary Programming, the capacitor values are chosen randomly from 150 to 1500 KVAR as listed from Table 2 (Eajal and ElHawary, 2010).
RESULTS AND DISCUSSION
The proposed Evolutionary Programming approach is utilized to identify the optimal shunt capacitor sizing in modified IEEE 13 bus unbalanced Radial distribution system is deployed in MATLAB computing environment on Intel Duo Core of 1.60 GHz and RAM of 1 GB.

Fig. 3: 
Modified IEEE 13bus unbalanced radial distribution system (URDS) 
The modified IEEE 13 busURDS consists of three, two and single phase lines and loads. The total active power loads on the system is 3466 kW and the total reactive power load is 2102 kVAR, respectively. The system loads considered as distributed loads and spot loads. Bus 1 is considered as substation and act as slack bus. The remaining buses (213) except the slack bus are modeled as constant PQ buses. 10 MVA is taken as the base value and the line to line base voltage is the feeder nominal voltage 4.16 kV. Throughout the optimization process, the limits imposed on the rms values of bus voltages are selected between Vmin = 0.9 pu and Vmax = 1.1pu. The cost of real power losses are taken as 168 U.S.$ kW^{1} year^{1} (Eajal and ElHawary, 2010). The pictographic view of the modified IEEE 13 bus unbalanced radial distribution system as shown in Fig. 3. A detailed description of IEEE 13 bus distribution network can be found by Kersting (1991).
Installation of shunt capacitor in bus number 671 at phase “3” with the size of 450 KVAR (case 2) results the total active power loss reduction from 116.4267 to 99.76 kW and net savings achieved is 2686.15 $ year^{1}. The total cost of RPL before the capacitive compensation is 19559.68 $ year^{1} (case 1) and after capacitive compensation, it is reduced to 16873.53 $ year^{1} (case 2). Comparatively, the reactive power compensation by shunt capacitor installed at optimal place has reduced the overall cost.
Solutions obtained from the proposed method are tabulated in Table 3 and 4. which indicates the improvement in voltage profiles after the installation of shunt capacitor at appropriate node. The best location and suitable size of shunt capacitor is identified to achieve power loss reduction and summary of results are tabulated in Table 5.
Table 3:  Bus voltage profile (magnitude and phase angles) before capacitor installation 

Table 4:  Bus voltage profile (magnitude and phase angles) after capacitor installation 

Table 5:  Simulation result before and after capacitor placement 

CONCLUSION
The proposed Evolutionary Programming approach is used to identify the optimum capacitor sizing in modified IEEE 13 bus URDS. This results in reduction of total active power loss with voltage profile improvement. The backward/forward sweep based load flow technique is used to estimate the total real power loss of the test system and bus voltage profiles. The suggested EP technique to find the optimum size of shunt capacitors has been examined on the modified IEEE 13 bus unbalanced radial distribution system and simulation results are presented. The simulations result shows that EP is well suited for multi objective reactive power planning, particularly in noncontinuous, nonsmooth situations. This approach can be extended to IEEE standard bus systems to test the accuracy of the proposed algorithm.