<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">CS</journal-id><journal-title-group><journal-title>Circuits and Systems</journal-title></journal-title-group><issn pub-type="epub">2153-1285</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/cs.2016.711304</article-id><article-id pub-id-type="publisher-id">CS-70453</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Computer Science&amp;Communications</subject><subject> Engineering</subject><subject> Physics&amp;Mathematics</subject></subj-group></article-categories><title-group><article-title>
 
 
  Hybrid Swarm Algorithm for Multiobjective Optimal Power Flow Problem
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>K.</surname><given-names>Rajalashmi</given-names></name><xref ref-type="aff" rid="aff1"><sup>1</sup></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>S.</surname><given-names>U. Prabha</given-names></name><xref ref-type="aff" rid="aff2"><sup>2</sup></xref></contrib></contrib-group><aff id="aff2"><addr-line>Department of Electrical and Electronics Engineering, Sri Ramakrishna Engineering College, Coimbatore, India</addr-line></aff><aff id="aff1"><addr-line>Department of Electrical and Electronics Engineering, Bannari Amman Institute of Technology, Sathyamangalam, India</addr-line></aff><author-notes><corresp id="cor1">* E-mail:<email>rajalashmik@bitsathy.ac.in(KR)</email>;</corresp></author-notes><pub-date pub-type="epub"><day>06</day><month>09</month><year>2016</year></pub-date><volume>07</volume><issue>11</issue><fpage>3589</fpage><lpage>3603</lpage><history><date date-type="received"><day>May</day>	<month>17,</month>	<year>2016</year></date><date date-type="rev-recd"><day>Accepted:</day>	<month>May</month>	<year>28,</year>	</date><date date-type="accepted"><day>September</day>	<month>8,</month>	<year>2016</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  Optimal power flow problem plays a major role in the operation and planning of power systems. It assists in acquiring the optimized solution for the optimal power flow problem. It consist
  s
   of several objective functions and constraints. This paper solves the multiobjective optimal power flow problem using a new hybrid technique by combining the particle swarm optimization and ant colony optimization. This hybrid method
   
  overcome
  s
   the drawback in local search such as stagnation and premature convergence and also enhances the global search with chemical communication signal. The best results are extracted using fuzzy approach from the hybrid algorithm solution. These methods have been examined with the power flow objectives such as cost, loss and voltage stability index by individuals and multiobjective functions. The proposed algorithms applied to IEEE 30 and IEEE
   
  118
  -
  bus test system and the results are analyzed and validated. The proposed algorithm results record the best compromised solution with minimum execution time compared with the particle swarm optimization
  .
 
</p></abstract><kwd-group><kwd>Multiobjective</kwd><kwd> OPF</kwd><kwd> Optimization</kwd><kwd> Hybrid</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>The power system is an interconnected electric network, which has generating plant and loads that are connected through transmission and distribution networks. It is a complicated network and it has many objectives to be solved. The reliable result of the objectives is obtained using optimal power flow. Optimized solution for the objective functions is attained by satisfying the power flow equations and constraints of the entire power system network. Various control variables are also influenced to achieve an optimal solution. Many conventional techniques are applied to solve the optimal power flow problem. They are quadratic programming [<xref ref-type="bibr" rid="scirp.70453-ref1">1</xref>] , interior programming [<xref ref-type="bibr" rid="scirp.70453-ref2">2</xref>] , linear programming [<xref ref-type="bibr" rid="scirp.70453-ref3">3</xref>] , gradient method [<xref ref-type="bibr" rid="scirp.70453-ref4">4</xref>] and newton’s method [<xref ref-type="bibr" rid="scirp.70453-ref5">5</xref>] . The drawback of these methods is that it is a local solver only. The solutions from these classical methods are not the most reasonable for nonlinear, discrete and continuous functions with constraints.</p><p>In order to overcome the drawbacks of classical methods, many computational methods are used to solve the optimal power flow problem. Particle Swarm Optimization (PSO) [<xref ref-type="bibr" rid="scirp.70453-ref6">6</xref>] inspires the school of fishes and a flock of birds and trapped in local optima. Genetic algorithm [<xref ref-type="bibr" rid="scirp.70453-ref7">7</xref>] [<xref ref-type="bibr" rid="scirp.70453-ref8">8</xref>] solves the optimal power flow problem using a random selection of starting values of voltage angles, but it has no control over global and local investigation. Evolutionary programming [<xref ref-type="bibr" rid="scirp.70453-ref9">9</xref>] algorithm solves the problem by enhancing the gradient of the function, and it has the problem of long computational time. Tabu search [<xref ref-type="bibr" rid="scirp.70453-ref10">10</xref>] solves the optimal power flow problem which is more composite for classical method. These manuscripts give the optimal solution for individual objective. Conversely the real world problem needs a concurrent optimal solution for all the objectives.</p><p>Multiobjective optimal power flow problem is solved by differential evolution, for active and reactive power dispatch [<xref ref-type="bibr" rid="scirp.70453-ref11">11</xref>] and fuel cost, voltage profile, and voltage stability enhancement [<xref ref-type="bibr" rid="scirp.70453-ref12">12</xref>] . The drawback of this technique is that the convergence time is more. Enhanced genetic algorithm [<xref ref-type="bibr" rid="scirp.70453-ref13">13</xref>] solves the problem with the combination of new developed quadratic load flow and its demerit is that it gets trapped in local optima. Self adaptive evolutionary programming [<xref ref-type="bibr" rid="scirp.70453-ref14">14</xref>] combines the evolutionary programming and random search technique and it has the problem of long computational time. As particle swarm optimization is improved [<xref ref-type="bibr" rid="scirp.70453-ref15">15</xref>] , the particles avoid being trapped in local optima. Harmony search algorithm [<xref ref-type="bibr" rid="scirp.70453-ref16">16</xref>] inspires the analogy with music improvisation which has solved the problem with slow convergence. These techniques solve the problem with the long computation time and being trapped in local optima.</p><p>In the weighted sum method [<xref ref-type="bibr" rid="scirp.70453-ref17">17</xref>] , different objectives are formatted into a single function. In this method, the demerit is to run the problem with multiple time to desire the optimal solution. Voltage security costs [<xref ref-type="bibr" rid="scirp.70453-ref18">18</xref>] are optimized using a penalty method which has the drawback of deciding the suitable penalty. The importance of each objective is not preserved while merging all the objective function into a single objective function. In order to maintain the significance of each objective pareto optimal method is adopted.</p><p>In the PSO algorithm, if the particles are small, local minimum problem will occur and the outcome could be evaluated in multiple runs. If the particles are too large, the global solution is obtained in each run and it reduces the speed of the algorithm. In order to overcome these difficulties and to increase the efficiency of PSO in this paper, a new hybrid technique technique is proposed by combining the particle swarm optimization and ant colony optimization with chemical communication signal.</p><p>In this paper cost, loss and Voltage Stability Index (VSI) objectives are analyzed using enhanced PSO with the multiobjective optimal power flow by pareto optimal me- thod. The multiobjective optimization of these objectives is practiced through three cases. In the first case, all the objective functions are simulated individually. In the second case, two objectives are considered simultaneously for optimization such as cost-loss, cost-VSI, loss-VSI and in the third case three objectives are considered concurrently as cost-loss-VSI. These objectives are optimized using the proposed technique with the IEEE 30 and IEEE 118-bus test system. A short introduction of an optimal power flow system has been presented in this section. This paper is planned as a tag on the section. In Section 2, the optimal power flow problem is formulated and discussed. Section 3 explains the central concepts of PSO, Ant colony Optimization (ACO) and Hybrid Swarm Optimization. Section 4 presents the simulated results with discussions of PSO, ACO and hybrid PSO. The conclusion is dealt in the Section 5.</p></sec><sec id="s2"><title>2. Problem Formulation for OPF Solution</title><p>In this paper the cost, transmission loss, voltage stability index are optimized, which is satisfying equality and voltage, generator, shunt VAR and transformer constraints.</p><sec id="s2_1"><title>2.1. Problem Objectives</title><p>1) Minimization of Fuel Cost:</p><p>This objective is to minimize the cost. The total fuel cost can be expressed as</p><disp-formula id="scirp.70453-formula268"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/13-7601183x2.png"  xlink:type="simple"/></disp-formula><p><img src="http://html.scirp.org/file/13-7601183x4.png" /><img src="http://html.scirp.org/file/13-7601183x3.png" /> (2)</p><disp-formula id="scirp.70453-formula269"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/13-7601183x5.png"  xlink:type="simple"/></disp-formula><p>a<sub>i</sub>, b<sub>i</sub>, c<sub>i</sub> is the cost coefficients of the ith generator, P<sub>i</sub> is the real power output of the ith generator.</p><p>2) Minimization of Transmission Loss:</p><p>This objective is to minimize the real power transmission loss. The transmission loss can be stated as</p><disp-formula id="scirp.70453-formula270"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/13-7601183x6.png"  xlink:type="simple"/></disp-formula><p>V<sub>i</sub> is the voltage at the ith line, V<sub>j</sub> is the voltage at the jth line, g<sub>k</sub> conductance.</p><p>3) Minimization of Voltage Stability Index:</p><p>This objective is to minimize the voltage stability index. The voltage level ought to maintain below critical level. If the voltage level greater than the critical limit the entire system becomes unstable. The fast voltage stability index expressed as</p><disp-formula id="scirp.70453-formula271"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/13-7601183x7.png"  xlink:type="simple"/></disp-formula><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x8.png" xlink:type="simple"/></inline-formula>is the impedance of line between i and j, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x9.png" xlink:type="simple"/></inline-formula>is the reactance of line between i and j, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x10.png" xlink:type="simple"/></inline-formula> is the sending end voltage, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x11.png" xlink:type="simple"/></inline-formula>is the receiving end reactive power.</p></sec><sec id="s2_2"><title>2.2. Problem Constraints</title><p>1) Voltage Constraint:</p><p>The bus voltage maintained the maximum and minimum level. The voltage limit conveyed as</p><disp-formula id="scirp.70453-formula272"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/13-7601183x12.png"  xlink:type="simple"/></disp-formula><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x13.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x14.png" xlink:type="simple"/></inline-formula> Minimum, maximum values of voltage magnitudes at bus i.</p><p>2) Generation Constraints:</p><p>The real and reactive power of generation sustained the maximum and minimum limit. The generation level of real and reactive power is expressed as</p><disp-formula id="scirp.70453-formula273"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/13-7601183x15.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.70453-formula274"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/13-7601183x16.png"  xlink:type="simple"/></disp-formula><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x17.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x18.png" xlink:type="simple"/></inline-formula> Minimum, maximum values of real power allowed at bus i.</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x19.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x19.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x20.png" xlink:type="simple"/></inline-formula> Minimum, maximum values of reactive power allowed at bus i.</p><p>3) Transformer Constraints:</p><p>The transformer tap setting values are maintained within the maximum and minimum level. The transformer tap setting limit is expressed as</p><disp-formula id="scirp.70453-formula275"><label>(10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/13-7601183x21.png"  xlink:type="simple"/></disp-formula><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x22.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x23.png" xlink:type="simple"/></inline-formula> Minimum, maximum values of transformer tap settings allowed at bus i.</p><p>4) Shunt VAR Constraints:</p><p>The shunt VAR restricted the maximum and minimum level. The shunt VAR limit is stated as</p><disp-formula id="scirp.70453-formula276"><label>(11)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/13-7601183x24.png"  xlink:type="simple"/></disp-formula><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x25.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x25.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x26.png" xlink:type="simple"/></inline-formula>minimum, maximum values of Shunt VAR allowed at bus i.</p><p>5) Equality Constraints:</p><p>The equality constraints are expressed as</p><disp-formula id="scirp.70453-formula277"><label>(12)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/13-7601183x27.png"  xlink:type="simple"/></disp-formula><disp-formula id="scirp.70453-formula278"><label>(13)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/13-7601183x28.png"  xlink:type="simple"/></disp-formula><p>n: Number of buses.</p><p>P<sub>Gi</sub>, Q<sub>Gi</sub> Real, reactive power generation at bus i.</p><p>P<sub>Di</sub>, Q<sub>Di</sub> Real, reactive load demand at bus i.</p><p>V<sub>i</sub>, V<sub>j</sub> Voltage magnitudes at bus i and j respectively.</p><p>G<sub>ij</sub>, B<sub>ij</sub> Transfer conductance and susceptance between bus i and j respectively.</p><p>δ<sub>i</sub>, δ<sub>j</sub> V oltage angle at bus i and j respectively.</p></sec></sec><sec id="s3"><title>3. Concepts of Swarm Intelligence</title><sec id="s3_1"><title>3.1. Basics of PSO</title><p>Particle swarm optimization (PSO) is a population based optimization technique which is motivated by group manners of birds or fishes. Each particle has its own experience and neighbor particle experience. According to that particle adjusts its position, within the bounded area. Each particle has its local best (lbest) or a personal best. The best position brings into being by all neighbours particles in the specific solution is global best (gbest). The optimal solution is acquired in the specific solution in the course of its current velocity and the experience.</p><p>The formation of individual best is</p><disp-formula id="scirp.70453-formula279"><label>(14)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/13-7601183x29.png"  xlink:type="simple"/></disp-formula><p>The pattern of global best is</p><disp-formula id="scirp.70453-formula280"><label>(15)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/13-7601183x30.png"  xlink:type="simple"/></disp-formula><p>The velocity and the positions are updated in the each iteration. The updating of velocity is</p><disp-formula id="scirp.70453-formula281"><label>(16)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/13-7601183x31.png"  xlink:type="simple"/></disp-formula><p>The updating of position of the particle is</p><disp-formula id="scirp.70453-formula282"><label>(17)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/13-7601183x32.png"  xlink:type="simple"/></disp-formula><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x33.png" xlink:type="simple"/></inline-formula>current velocity of individual i at iteration k, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x34.png" xlink:type="simple"/></inline-formula>Modified velocity of individual i at iteration k + 1, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x35.png" xlink:type="simple"/></inline-formula>Current position of individual i at iteration, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x36.png" xlink:type="simple"/></inline-formula>Inertia weight parameter, c<sub>1</sub>, c<sub>2</sub> Acceleration factors, rand<sub>1</sub>, rand<sub>2</sub> Random numbers between 0 and 1, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x37.png" xlink:type="simple"/></inline-formula>Best position of individual i until iteration k, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x33.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x36.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x37.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x38.png" xlink:type="simple"/></inline-formula>Best position of the group until iteration k.</p></sec><sec id="s3_2"><title>3.2. Ant Colony Optimization</title><p>Ant colony optimization (ACO) is an evolutionary and adaptive algorithm inspired by the behavior of real ant colonies. The ant when searching food and find food sources, it deposits a chemical called pheromone as a trail during the return path. Based on the quantity and quality of the food available the quantity of pheromone deposited. The pheromone trails allow the ants to identify the shortest path between the ant nest and the food source. Ant colony optimization algorithm updates pheromone values in order to update solutions during run time [<xref ref-type="bibr" rid="scirp.70453-ref19">19</xref>] .</p><p>The pheromone trail updating equation is stated as</p><disp-formula id="scirp.70453-formula283"><label>(18)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/13-7601183x39.png"  xlink:type="simple"/></disp-formula><p>α is the pheromone trails evaporation co-efficient, which lies between 0 and 1, ΔT<sub>ij</sub> is the increment of edge for the period Δt.</p><p>To enhance the tempo of search in the optimal power flow problem, the chemical communication between the insects will help. It is the signal of wasp which produces alarm pheromone to the inspiration of anxiety in the colony. This increases the speed of search in the problem area.</p><p>The chemical communication signal is expressed as</p><disp-formula id="scirp.70453-formula284"><label>(19)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/13-7601183x40.png"  xlink:type="simple"/></disp-formula><p>A is the chemical signal, B is the emission rate, D is the diffusion coefficient, r is the radius of the active space, t is the time from the beginning of emission, Comperr is the complementary error function.</p></sec><sec id="s3_3"><title>3.3. Hybrid Swarm Intelligence</title><p>ACO not suitable for large optimization problems because it takes long time to search the result and premature convergence. To overcome of finding the best solution with large search space, combine with another algorithm (PSO) which has a better solution. In the hybrid PSO-ACO, the PSO algorithm is used for pheromone update of ACO. This improves the global exploration capabilities of ACO algorithm. This also improves convergence performance of ACO.</p><p>The velocity updating equation is expressed as</p><disp-formula id="scirp.70453-formula285"><label>(20)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/13-7601183x41.png"  xlink:type="simple"/></disp-formula><p>A is the chemical communication signal, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x42.png" xlink:type="simple"/></inline-formula> is the increment of pheromone trails laid on the edge by ant between t and t + Δt, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x43.png" xlink:type="simple"/></inline-formula>is the local update value of ant k, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x43.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x44.png" xlink:type="simple"/></inline-formula>is the global update value of edge.</p></sec><sec id="s3_4"><title>3.4. Fuzzy Approach</title><p>In the proposed method, the aim is to achieve a set of solution of the multiobjective optimal power flow problem. In the optimal power flow problem with conflicting objectives the single optimal solution does not determine the real solution. The Pareto optimal approach, it is feasible to find the set of solution instead of the single optimal solution. From the set of solution, the best solution is obtained from the fuzzy approach. The equation for the fuzzy approach for best solution is expressed as</p><disp-formula id="scirp.70453-formula286"><label>(21)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/13-7601183x45.png"  xlink:type="simple"/></disp-formula><p>x total number of non-dominated solutions.</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x46.png" xlink:type="simple"/></inline-formula>= 0 when the objective function is greater than the maximum value of ith objective function.</p><p><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/13-7601183x47.png" xlink:type="simple"/></inline-formula>= 1 when the objective function is greater than the minimum value of ith objective function.</p></sec><sec id="s3_5"><title>3.5. Proposed Hybrid Algorithm for the Optimal Power Flow Problem</title><p>Step 1: Read input data include equality and inequality constraints.</p><p>Step 2: Initialize the population.</p><p>Step 3: Calculate the objective function.</p><p>Step 4: Apply Pareto optimal method to get the set of solution.</p><p>Step 5: Update the velocity using hybrid velocity updating equation.</p><p>Step 6: Update the position using position updating equation.</p><p>Step 7: Check for stopping criteria. If the iteration reaches the maximum iteration stop else go to step 5.</p><p>Step 8: Apply Pareto method and determine the non-dominated solution.</p><p>Step 9: Find the best solution using fuzzy approach.</p></sec></sec><sec id="s4"><title>4. Simulation Results</title><p>The proposed system has been tested to IEEE 30-bus test system. The 30-bus IEEE test system has 41 transmission lines, six generators and four transformers (T6-9, T6-10, T4-12 and T27-28). The transformer taps are set between 0.9 to 1.1. The proposed algorithm has the particles of 30.</p><p>Case 1: Single objective optimization.</p><p>The proposed algorithms analyzed with the various values of acceleration factors (c<sub>1</sub> and c<sub>2</sub>) individually. This analysis assists to achieve the suitable value of the acceleration factor for the further process. This optimal power flow problem was simulated with the acceleration factor (c<sub>1</sub>) varied from 1 to 3, using hybrid algorithm. The best value is obtained at c<sub>1</sub> &gt; 2.5. <xref ref-type="table" rid="table1">Table 1</xref> shows the best values of the cost, loss and VSI for c<sub>1</sub> = 2.5 are $801.7/hr, 5.21 MW and 0.14 p.u. respectively. The hybrid algorithm also simulated with the various values of the acceleration factor (c<sub>2</sub>). The best value is obtained at c<sub>2</sub> &gt; 2. <xref ref-type="table" rid="table2">Table 2</xref> shows the best values of the cost, loss and VSI for c<sub>2</sub> = 2 are, $801.72/hr, 5.24 MW and 0.13 p.u. respectively.</p><p><xref ref-type="table" rid="table3">Table 3</xref> shows the single objective optimized values with the control variables using the algorithms of PSO, ACO and Hybrid algorithm respectively. It shows that, the reliable generation cost of PSO, ACO and Hybrid are $802.15/hr, $804.86/hr and $801.54/hr respectively. <xref ref-type="fig" rid="fig1">Figure 1</xref> displays the generation cost convergence characteristics of a single objective optimization of IEEE 30 test bus system using PSO, ACO and Hybrid algorithm. <xref ref-type="fig" rid="fig2">Figure 2</xref> indicates the reduction of transmission loss using PSO, ACO and Hybrid algorithms is 5.974 MW, 9.42 MW and 5.208 MW respectively. <xref ref-type="fig" rid="fig3">Figure 3</xref> indicates the voltage stability of PSO, ACO and Hybrid algorithms are 0.1307 p.u., 0.142 p.u. and 0.13 p.u. respectively. These results show that, hybrid algorithm results are distributed.</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> Best optimized values of different C<sub>1</sub> using hybrid algorithm</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Acceleration factor</th><th align="center" valign="middle" >1</th><th align="center" valign="middle" >1.5</th><th align="center" valign="middle" >2</th><th align="center" valign="middle" >2.5</th><th align="center" valign="middle" >3</th></tr></thead><tr><td align="center" valign="middle" >Cost ($/hr)</td><td align="center" valign="middle" >803.7</td><td align="center" valign="middle" >802.4</td><td align="center" valign="middle" >801.9</td><td align="center" valign="middle" >801.7</td><td align="center" valign="middle" >801.9</td></tr><tr><td align="center" valign="middle" >Loss (MW)</td><td align="center" valign="middle" >5.41</td><td align="center" valign="middle" >5.322</td><td align="center" valign="middle" >5.24</td><td align="center" valign="middle" >5.21</td><td align="center" valign="middle" >5.91</td></tr><tr><td align="center" valign="middle" >VSI (p.u.)</td><td align="center" valign="middle" >0.151</td><td align="center" valign="middle" >0.15</td><td align="center" valign="middle" >0.149</td><td align="center" valign="middle" >0.14</td><td align="center" valign="middle" >0.153</td></tr></tbody></table></table-wrap><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> Best optimized values of different C<sub>2</sub> using hybrid algorithm</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Acceleration factor</th><th align="center" valign="middle" >1</th><th align="center" valign="middle" >1.5</th><th align="center" valign="middle" >2</th><th align="center" valign="middle" >2.5</th><th align="center" valign="middle" >3</th></tr></thead><tr><td align="center" valign="middle" >Cost ($/hr)</td><td align="center" valign="middle" >803.89</td><td align="center" valign="middle" >802.75</td><td align="center" valign="middle" >801.72</td><td align="center" valign="middle" >801.76</td><td align="center" valign="middle" >801.78</td></tr><tr><td align="center" valign="middle" >Loss (MW)</td><td align="center" valign="middle" >5.319</td><td align="center" valign="middle" >5.27</td><td align="center" valign="middle" >5.24</td><td align="center" valign="middle" >5.25</td><td align="center" valign="middle" >5.254</td></tr><tr><td align="center" valign="middle" >VSI (p.u)</td><td align="center" valign="middle" >0.131</td><td align="center" valign="middle" >0.129</td><td align="center" valign="middle" >0.13</td><td align="center" valign="middle" >0.132</td><td align="center" valign="middle" >0.131</td></tr></tbody></table></table-wrap><table-wrap id="table3" ><label><xref ref-type="table" rid="table3">Table 3</xref></label><caption><title> Single objective optimized values of cost, loss and VSI using PSO, ACO and Hybrid algorithms</title></caption><table><tbody><thead><tr><th align="center" valign="middle" ></th><th align="center" valign="middle" >PSO</th><th align="center" valign="middle" >ACO</th><th align="center" valign="middle" >Hybrid</th></tr></thead><tr><td align="center" valign="middle" >P<sub>g1</sub> (MW)</td><td align="center" valign="middle" >177.132</td><td align="center" valign="middle" >176.901</td><td align="center" valign="middle" >176.7</td></tr><tr><td align="center" valign="middle" >P<sub>g2</sub> (MW)</td><td align="center" valign="middle" >48.12</td><td align="center" valign="middle" >49.932</td><td align="center" valign="middle" >48.89</td></tr><tr><td align="center" valign="middle" >P<sub>g3</sub> (MW)</td><td align="center" valign="middle" >21.356</td><td align="center" valign="middle" >21.125</td><td align="center" valign="middle" >21.4719</td></tr><tr><td align="center" valign="middle" >P<sub>g4</sub> (MW)</td><td align="center" valign="middle" >22.878</td><td align="center" valign="middle" >21.468</td><td align="center" valign="middle" >21.6423</td></tr><tr><td align="center" valign="middle" >P<sub>g5</sub> (MW)</td><td align="center" valign="middle" >10.121</td><td align="center" valign="middle" >10.116</td><td align="center" valign="middle" >12.0878</td></tr><tr><td align="center" valign="middle" >P<sub>g6</sub> (MW)</td><td align="center" valign="middle" >12</td><td align="center" valign="middle" >12</td><td align="center" valign="middle" >12.00</td></tr><tr><td align="center" valign="middle" >V<sub>g1</sub></td><td align="center" valign="middle" >1.05</td><td align="center" valign="middle" >1.045</td><td align="center" valign="middle" >1.045</td></tr><tr><td align="center" valign="middle" >V<sub>g2</sub></td><td align="center" valign="middle" >1.0442</td><td align="center" valign="middle" >1.043</td><td align="center" valign="middle" >1.043</td></tr><tr><td align="center" valign="middle" >V<sub>g3</sub></td><td align="center" valign="middle" >1.446</td><td align="center" valign="middle" >0.998</td><td align="center" valign="middle" >0.998</td></tr><tr><td align="center" valign="middle" >V<sub>g4</sub></td><td align="center" valign="middle" >1.0408</td><td align="center" valign="middle" >1.009</td><td align="center" valign="middle" >1.009</td></tr><tr><td align="center" valign="middle" >V<sub>g5</sub></td><td align="center" valign="middle" >0.9601</td><td align="center" valign="middle" >1.014</td><td align="center" valign="middle" >1.014</td></tr><tr><td align="center" valign="middle" >V<sub>g6</sub></td><td align="center" valign="middle" >1.05</td><td align="center" valign="middle" >1.047</td><td align="center" valign="middle" >1.047</td></tr><tr><td align="center" valign="middle" >T<sub>6-9</sub>, (p.u.)</td><td align="center" valign="middle" >1.01</td><td align="center" valign="middle" >1.012</td><td align="center" valign="middle" >1.012</td></tr><tr><td align="center" valign="middle" >T<sub>6-10</sub>, (p.u.)</td><td align="center" valign="middle" >0.99</td><td align="center" valign="middle" >0.971</td><td align="center" valign="middle" >0.971</td></tr><tr><td align="center" valign="middle" >T<sub>4-12</sub>, (p.u.)</td><td align="center" valign="middle" >1.01</td><td align="center" valign="middle" >1.023</td><td align="center" valign="middle" >1.023</td></tr><tr><td align="center" valign="middle" >T<sub>27-28</sub>, (p.u.)</td><td align="center" valign="middle" >1.02</td><td align="center" valign="middle" >1.014</td><td align="center" valign="middle" >1.014</td></tr><tr><td align="center" valign="middle" >Cost ($/hr)</td><td align="center" valign="middle" >802.15</td><td align="center" valign="middle" >804.86</td><td align="center" valign="middle" >801.54</td></tr><tr><td align="center" valign="middle" >Loss (MW)</td><td align="center" valign="middle" >5.974</td><td align="center" valign="middle" >9.42</td><td align="center" valign="middle" >5.208</td></tr><tr><td align="center" valign="middle" >VSI (p.u)</td><td align="center" valign="middle" >0.1307</td><td align="center" valign="middle" >0.142</td><td align="center" valign="middle" >0.13</td></tr></tbody></table></table-wrap><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> The generation cost convergence characteristics of a single objective optimization of IEEE 30 test bus system using PSO, ACO and hybrid algorithm</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/13-7601183x48.png"/></fig><p>The accuracy of the proposed system is measured using the mean and standard deviation. <xref ref-type="table" rid="table4">Table 4</xref> indicates the mean and standard values of the PSO, ACO and Hybrid algorithms. It is obvious that the proposed hybrid algorithm can achieve better result in the best solution search. The convergence time for the cost optimization using a hybrid technique is 1.05 s, which is a faster than the PSO and ACO methods (1.24 s and 1.12 s</p><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> The transmission loss convergence characteristics of a single objective optimization of IEEE 30 test bus system using PSO, ACO and Hybrid algorithm</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/13-7601183x49.png"/></fig><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> The voltage stability index convergence characteristics of a single objective optimization of IEEE 30 test bus system using PSO, ACO and hybrid algorithm</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/13-7601183x50.png"/></fig><table-wrap id="table4" ><label><xref ref-type="table" rid="table4">Table 4</xref></label><caption><title> Performance measure of PSO, ACO and hybrid algorithms</title></caption><table><tbody><thead><tr><th align="center" valign="middle" ></th><th align="center" valign="middle" >ACO</th><th align="center" valign="middle" >PSO</th><th align="center" valign="middle" >Hybrid</th></tr></thead><tr><td align="center" valign="middle" >Mean</td><td align="center" valign="middle" >0.212</td><td align="center" valign="middle" >0.234</td><td align="center" valign="middle" >0.259</td></tr><tr><td align="center" valign="middle" >STD</td><td align="center" valign="middle" >0.042</td><td align="center" valign="middle" >0.036</td><td align="center" valign="middle" >0.024</td></tr></tbody></table></table-wrap><p>respectively). The convergence time taken for the hybrid technique is less. The accuracy and speed of convergence of hybrid are more than the ACO and PSO. Because the hybrid algorithm overcome the drawback in local search such as stagnation and premature convergence by the local experience of the ant and it processed in parallel. This enhances the chance of locating the finest solution in rapidity. The global search increased by the divergence of the pheromone tracks and also with the chemical communication signal. This leads the particle move the quickly to reach the solution. The global utilization, speed up the search and it support to acquire the best solution.</p><p>This optimal power flow problem was worked out using P4, 3 GHz. To evaluate the difficulty of the proposed system, the IEEE 14 bus is chosen to solve this OPF problem. The convergence time for this optimization in IEEE 14 bus system of hybrid method is 0.63 s. The convergence time for the hybrid technique is for IEEE 30 bus test system is 1.05 s The convergence speed of the IEEE 30 bus system is twice the times of IEEE 14 bus test system. This will persuade to carry on the proposed technique in higher buses.</p><p>Case 2: Two objective optimization.</p><p>In a real world, the problems involve simultaneous optimization of several objective functions. These functions are non commensurable and conflicting objective functions. It gives rise to a set of optimal solutions, instead of one optimal solution. The reason for the optimality of many solutions is that no one can be considered to be better than any other with respect to all objective functions. <xref ref-type="table" rid="table5">Table 5</xref> shows the optimized values of two</p><table-wrap id="table5" ><label><xref ref-type="table" rid="table5">Table 5</xref></label><caption><title> Optimal solution for objective functions with different algorithms</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  rowspan="2"  ></th><th align="center" valign="middle"  colspan="3"  >Cost-Loss</th><th align="center" valign="middle"  colspan="3"  >Cost-VSI</th><th align="center" valign="middle"  colspan="3"  >VSI-Loss</th></tr></thead><tr><td align="center" valign="middle" >PSO</td><td align="center" valign="middle" >ACO</td><td align="center" valign="middle" >Hybrid</td><td align="center" valign="middle" >PSO</td><td align="center" valign="middle" >ACO</td><td align="center" valign="middle" >Hybrid</td><td align="center" valign="middle" >PSO</td><td align="center" valign="middle" >ACO</td><td align="center" valign="middle" >Hybrid</td></tr><tr><td align="center" valign="middle" >Pg1, MW</td><td align="center" valign="middle" >136.467</td><td align="center" valign="middle" >132.374</td><td align="center" valign="middle" >134.717</td><td align="center" valign="middle" >178.167</td><td align="center" valign="middle" >179.232</td><td align="center" valign="middle" >178.12</td><td align="center" valign="middle" >177.458</td><td align="center" valign="middle" >172.486</td><td align="center" valign="middle" >169.451</td></tr><tr><td align="center" valign="middle" >Pg2, MW</td><td align="center" valign="middle" >49.396</td><td align="center" valign="middle" >51.585</td><td align="center" valign="middle" >48.856</td><td align="center" valign="middle" >46.544</td><td align="center" valign="middle" >47.54</td><td align="center" valign="middle" >45.421</td><td align="center" valign="middle" >45.425</td><td align="center" valign="middle" >49.581</td><td align="center" valign="middle" >46.859</td></tr><tr><td align="center" valign="middle" >Pg3, MW</td><td align="center" valign="middle" >32.826</td><td align="center" valign="middle" >27.847</td><td align="center" valign="middle" >31.471</td><td align="center" valign="middle" >21.412</td><td align="center" valign="middle" >21.324</td><td align="center" valign="middle" >21.433</td><td align="center" valign="middle" >23.192</td><td align="center" valign="middle" >21.247</td><td align="center" valign="middle" >22.294</td></tr><tr><td align="center" valign="middle" >Pg4, MW</td><td align="center" valign="middle" >21.686</td><td align="center" valign="middle" >32.694</td><td align="center" valign="middle" >30.642</td><td align="center" valign="middle" >16.234</td><td align="center" valign="middle" >17.415</td><td align="center" valign="middle" >16.574</td><td align="center" valign="middle" >17.497</td><td align="center" valign="middle" >19.782</td><td align="center" valign="middle" >18.482</td></tr><tr><td align="center" valign="middle" >Pg5, MW</td><td align="center" valign="middle" >18.085</td><td align="center" valign="middle" >28.878</td><td align="center" valign="middle" >19.078</td><td align="center" valign="middle" >12.201</td><td align="center" valign="middle" >13.127</td><td align="center" valign="middle" >12.401</td><td align="center" valign="middle" >12.108</td><td align="center" valign="middle" >13.875</td><td align="center" valign="middle" >12.158</td></tr><tr><td align="center" valign="middle" >Pg6, MW</td><td align="center" valign="middle" >26.011</td><td align="center" valign="middle" >27.200</td><td align="center" valign="middle" >26.002</td><td align="center" valign="middle" >14.704</td><td align="center" valign="middle" >12.332</td><td align="center" valign="middle" >12.228</td><td align="center" valign="middle" >14.704</td><td align="center" valign="middle" >13.229</td><td align="center" valign="middle" >12.759</td></tr><tr><td align="center" valign="middle" >Vg1</td><td align="center" valign="middle" >1.035</td><td align="center" valign="middle" >1.0265</td><td align="center" valign="middle" >1.045</td><td align="center" valign="middle" >1.049</td><td align="center" valign="middle" >1.056</td><td align="center" valign="middle" >1.045</td><td align="center" valign="middle" >1.049</td><td align="center" valign="middle" >1.056</td><td align="center" valign="middle" >1.045</td></tr><tr><td align="center" valign="middle" >Vg2</td><td align="center" valign="middle" >1.034</td><td align="center" valign="middle" >1.023</td><td align="center" valign="middle" >1.043</td><td align="center" valign="middle" >1.04</td><td align="center" valign="middle" >1.034</td><td align="center" valign="middle" >1.043</td><td align="center" valign="middle" >1.082</td><td align="center" valign="middle" >1.034</td><td align="center" valign="middle" >1.043</td></tr><tr><td align="center" valign="middle" >Vg3</td><td align="center" valign="middle" >0.998</td><td align="center" valign="middle" >0.948</td><td align="center" valign="middle" >0.998</td><td align="center" valign="middle" >0.890</td><td align="center" valign="middle" >0.998</td><td align="center" valign="middle" >0.998</td><td align="center" valign="middle" >1.089</td><td align="center" valign="middle" >0.998</td><td align="center" valign="middle" >0.998</td></tr><tr><td align="center" valign="middle" >Vg4</td><td align="center" valign="middle" >1.090</td><td align="center" valign="middle" >1.089</td><td align="center" valign="middle" >1.009</td><td align="center" valign="middle" >1.000</td><td align="center" valign="middle" >1.005</td><td align="center" valign="middle" >1.009</td><td align="center" valign="middle" >1.002</td><td align="center" valign="middle" >1.005</td><td align="center" valign="middle" >1.009</td></tr><tr><td align="center" valign="middle" >Vg5</td><td align="center" valign="middle" >1.0154</td><td align="center" valign="middle" >1.014</td><td align="center" valign="middle" >1.014</td><td align="center" valign="middle" >1.022</td><td align="center" valign="middle" >1.034</td><td align="center" valign="middle" >1.014</td><td align="center" valign="middle" >1.212</td><td align="center" valign="middle" >1.034</td><td align="center" valign="middle" >1.014</td></tr><tr><td align="center" valign="middle" >Vg6</td><td align="center" valign="middle" >1.0347</td><td align="center" valign="middle" >1.0477</td><td align="center" valign="middle" >1.047</td><td align="center" valign="middle" >1.057</td><td align="center" valign="middle" >1.064</td><td align="center" valign="middle" >1.047</td><td align="center" valign="middle" >1.057</td><td align="center" valign="middle" >1.064</td><td align="center" valign="middle" >1.047</td></tr><tr><td align="center" valign="middle" >T6-9, pu</td><td align="center" valign="middle" >1.0152</td><td align="center" valign="middle" >1.016</td><td align="center" valign="middle" >1.012</td><td align="center" valign="middle" >1.019</td><td align="center" valign="middle" >1.023</td><td align="center" valign="middle" >1.012</td><td align="center" valign="middle" >1.019</td><td align="center" valign="middle" >1.023</td><td align="center" valign="middle" >1.012</td></tr><tr><td align="center" valign="middle" >T6-10, pu</td><td align="center" valign="middle" >0.9671</td><td align="center" valign="middle" >0.9598</td><td align="center" valign="middle" >0.971</td><td align="center" valign="middle" >0.980</td><td align="center" valign="middle" >0.992</td><td align="center" valign="middle" >0.971</td><td align="center" valign="middle" >0.980</td><td align="center" valign="middle" >0.992</td><td align="center" valign="middle" >0.971</td></tr><tr><td align="center" valign="middle" >T4-12, pu</td><td align="center" valign="middle" >1.0123</td><td align="center" valign="middle" >1.002</td><td align="center" valign="middle" >1.023</td><td align="center" valign="middle" >1.034</td><td align="center" valign="middle" >1.039</td><td align="center" valign="middle" >1.023</td><td align="center" valign="middle" >1.034</td><td align="center" valign="middle" >1.039</td><td align="center" valign="middle" >1.023</td></tr><tr><td align="center" valign="middle" >T27-28, pu</td><td align="center" valign="middle" >1.0124</td><td align="center" valign="middle" >1.0133</td><td align="center" valign="middle" >1.014</td><td align="center" valign="middle" >1.034</td><td align="center" valign="middle" >1.045</td><td align="center" valign="middle" >1.014</td><td align="center" valign="middle" >1.034</td><td align="center" valign="middle" >1.045</td><td align="center" valign="middle" >1.014</td></tr><tr><td align="center" valign="middle" >Cost ($/hr)</td><td align="center" valign="middle" >825.15</td><td align="center" valign="middle" >831.72</td><td align="center" valign="middle" >821.72</td><td align="center" valign="middle" >809.54</td><td align="center" valign="middle" >815.65</td><td align="center" valign="middle" >802.5</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >Loss (MW)</td><td align="center" valign="middle" >5.376</td><td align="center" valign="middle" >5.637</td><td align="center" valign="middle" >5.003</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >5.2574</td><td align="center" valign="middle" >5.117</td><td align="center" valign="middle" >5.113</td></tr><tr><td align="center" valign="middle" >VSI (p.u)</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >0.1094</td><td align="center" valign="middle" >0.1096</td><td align="center" valign="middle" >0.1092</td><td align="center" valign="middle" >0.1047</td><td align="center" valign="middle" >0.1052</td><td align="center" valign="middle" >0.1007</td></tr><tr><td align="center" valign="middle" >Execution time (s)</td><td align="center" valign="middle" >1.4815</td><td align="center" valign="middle" >1.675</td><td align="center" valign="middle" >1.4167</td><td align="center" valign="middle" >1.598</td><td align="center" valign="middle" >1.679</td><td align="center" valign="middle" >1.5297</td><td align="center" valign="middle" >1.4298</td><td align="center" valign="middle" >1.4972</td><td align="center" valign="middle" >1.4475</td></tr></tbody></table></table-wrap><p>objective optimization such as cost-loss, cost-VSI and loss-VSI using PSO, ACO and Hybrid algorithm. It is not easy to pinpoint the correct concert solution from a single run for the corresponding objective. So, to get the appropriate solution 30 numbers of pare to solution sets are obtained from 30 generations.</p><p>In the cost-loss minimization, objectives contrasts to each other. The cost-loss divergence optimization, the cost and loss for the PSO, ACO and Hybrid techniques are $825.15/hr, $831.72/hr, $821.72/hr and 5.376 MW, 5.375 MW, 5.043 MW respectively. The reduction of cost and loss in hybrid is 1.12% and 11%higher than PSO and ACO techniques. <xref ref-type="fig" rid="fig4">Figure 4</xref> indicates the compromised solution of cost and loss are $821.72/hr and 5.43 MW obtained using hybrid algorithm. It shows that, the cost-loss multiobjec- tive optimization, the hybrid solution is best in cost and less reduction.</p><p>To avoid the stability problem in the system, the VSI should be maintained below critical limit. In addition to the VSI, cost also to be reduced. The cost-VSI conflict optimization, the cost and VSI for the PSO, ACO and Hybrid techniques are $809.54/hr, $815.65/hr, $802.5/hr and 0.1094 p.u, 0.1096 p.u., 0.1092 p.u, respectively. The reduction of cost and VSI of hybrid technique is 1.63% and 0.36% respectively. <xref ref-type="fig" rid="fig5">Figure 5</xref> shows that the best compromised solution of cost and VSI are $802.5/hr and 0.1092 p.u, obtained using hybrid algorithm. The execution time of the hybrid technique for the cost-VSI divergence optimization is 1.5297 s, which is a faster than the PSO and ACO methods (1.598 s and 1.679 s respectively).</p><p>In the loss-VSI case, the loss is minimized with the voltage stability, which is maintained below critical limit. The loss-VSI contrast optimization, the loss and VSI for the PSO, ACO and Hybrid techniques are 5.2574 MW, 5.117 MW, 5.113 MW and 0.1047 p.u, 0.1052 p.u., 0.1007 p.u, respectively. It shows that the reduction of loss in the ACO technique is less than the PSO technique <xref ref-type="fig" rid="fig6">Figure 6</xref> indicates the divergence of hybrid is</p><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> Loss-Cost optimization using PSO, ACO and hybrid algorithms of IEEE 30 test bus system</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/13-7601183x51.png"/></fig><fig id="fig5"  position="float"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Loss-voltage stability index optimization using PSO, ACO and hybrid algorithm of IEEE 30 test bus system</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/13-7601183x52.png"/></fig><fig id="fig6"  position="float"><label><xref ref-type="fig" rid="fig6">Figure 6</xref></label><caption><title> Voltage stability index-cost optimization using PSO, ACO and hybrid algorithm of IEEE 30 test bus system</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/13-7601183x53.png"/></fig><p>more than the ACO and PSO. The best compromised solution of loss and VSI are 5.113 MW and 0.1007 p.u, obtained using SPPSO algorithm. This result indicates the hybrid optimized result from the loss-VSI multiobjective optimization is better compared to the PSO and ACO techniques.</p><p>From the two objective optimization the reduction of cost is more in cost-VSI than the cost-loss optimization. The loss is optimized more in cost-loss multiobjective than the loss-VSI. Optimizing two objectives will give the dissimilar result for the same objective function. It initiates the objectives which are considered for the optimal power flow should be optimized simultaneously.</p><p>Case 3: Three objective optimization.</p><p><xref ref-type="table" rid="table6">Table 6</xref> indicates the optimal values of the three objectives such as cost-loss-VSI, using PSO, ACO and Hybrid algorithms. The cost of the three objective optimization using PSO, ACO and Hybrid algorithms are $814.76/hr, $815.6/hr and $805.5/hr respectively. In this, the cost reduced in the hybrid 1% more than the ACO and PSO techniques. The loss of the three objective optimization of the PSO, ACO and Hybrid techniques are 9.924 MW, 9.98 MW and 6.321 MW respectively. The loss reduction of hybrid technique is 57% higher than the PSO and ACO. The VSI of the PSO, ACO and Hybrid for the multi objective optimization are 0.1407 p.u, 0.1482 p.u and 0.1127 p.u respectively. The hybrid algorithm results, 24% more reduction in VSI than the PSO and ACO techniques. <xref ref-type="table" rid="table6">Table 6</xref> shows, the hybrid technique results the better accuracy compare to the PSO and ACO algorithm for multiobjective optimization problem.</p><p>In order to confirm the performance of the proposed algorithm in large systems, hybrid algorithm applied in IEEE 118 test bus system. The particles are increased to 50. To get the appropriate solution, 30 numbers of pareto solution sets are obtained from 30 generations. <xref ref-type="table" rid="table7">Table 7</xref> indicates the optimal solutions for the cost, loss and VSI optimized with a single objective, two objectives and three objective function using a hybrid</p><table-wrap id="table6" ><label><xref ref-type="table" rid="table6">Table 6</xref></label><caption><title> Multiobjective optimal solution for Cost-Loss-Voltage Stability Index with different algorithms</title></caption><table><tbody><thead><tr><th align="center" valign="middle" >Variables</th><th align="center" valign="middle" >PSO</th><th align="center" valign="middle" >ACO</th><th align="center" valign="middle" >Hybrid</th></tr></thead><tr><td align="center" valign="middle" >P<sub>g1</sub>, MW</td><td align="center" valign="middle" >168.163</td><td align="center" valign="middle" >167.937</td><td align="center" valign="middle" >105.71</td></tr><tr><td align="center" valign="middle" >P<sub>g2</sub>, MW</td><td align="center" valign="middle" >46.232</td><td align="center" valign="middle" >45.827</td><td align="center" valign="middle" >74.445</td></tr><tr><td align="center" valign="middle" >P<sub>g3</sub>, MW</td><td align="center" valign="middle" >23.957</td><td align="center" valign="middle" >23.211</td><td align="center" valign="middle" >38.329</td></tr><tr><td align="center" valign="middle" >P<sub>g4</sub>, MW</td><td align="center" valign="middle" >15.328</td><td align="center" valign="middle" >18.185</td><td align="center" valign="middle" >35.517</td></tr><tr><td align="center" valign="middle" >P<sub>g5</sub>, MW</td><td align="center" valign="middle" >18.541</td><td align="center" valign="middle" >17.248</td><td align="center" valign="middle" >19.278</td></tr><tr><td align="center" valign="middle" >P<sub>g6</sub>, MW</td><td align="center" valign="middle" >23.789</td><td align="center" valign="middle" >25.157</td><td align="center" valign="middle" >12.248</td></tr><tr><td align="center" valign="middle" >V<sub>g1</sub></td><td align="center" valign="middle" >1.065</td><td align="center" valign="middle" >1.074</td><td align="center" valign="middle" >1.045</td></tr><tr><td align="center" valign="middle" >V<sub>g2</sub></td><td align="center" valign="middle" >1.054</td><td align="center" valign="middle" >1.063</td><td align="center" valign="middle" >1.043</td></tr><tr><td align="center" valign="middle" >V<sub>g3</sub></td><td align="center" valign="middle" >0.9002</td><td align="center" valign="middle" >0.9098</td><td align="center" valign="middle" >0.998</td></tr><tr><td align="center" valign="middle" >V<sub>g4</sub></td><td align="center" valign="middle" >1.100</td><td align="center" valign="middle" >1.104</td><td align="center" valign="middle" >1.009</td></tr><tr><td align="center" valign="middle" >V<sub>g5</sub></td><td align="center" valign="middle" >1.017</td><td align="center" valign="middle" >1.019</td><td align="center" valign="middle" >1.014</td></tr><tr><td align="center" valign="middle" >V<sub>g6</sub></td><td align="center" valign="middle" >1.066</td><td align="center" valign="middle" >1.078</td><td align="center" valign="middle" >1.047</td></tr><tr><td align="center" valign="middle" >T<sub>6-9</sub>, pu</td><td align="center" valign="middle" >1.015</td><td align="center" valign="middle" >1.019</td><td align="center" valign="middle" >1.012</td></tr><tr><td align="center" valign="middle" >T<sub>6-10</sub>, pu</td><td align="center" valign="middle" >0.980</td><td align="center" valign="middle" >0.998</td><td align="center" valign="middle" >0.971</td></tr><tr><td align="center" valign="middle" >T<sub>4-12</sub>, pu</td><td align="center" valign="middle" >1.034</td><td align="center" valign="middle" >1.045</td><td align="center" valign="middle" >1.023</td></tr><tr><td align="center" valign="middle" >T<sub>27-28</sub>, pu</td><td align="center" valign="middle" >1.016</td><td align="center" valign="middle" >1.017</td><td align="center" valign="middle" >1.014</td></tr><tr><td align="center" valign="middle" >Cost ($/hr)</td><td align="center" valign="middle" >814.76</td><td align="center" valign="middle" >815.6</td><td align="center" valign="middle" >805.4</td></tr><tr><td align="center" valign="middle" >Loss (MW)</td><td align="center" valign="middle" >9.924</td><td align="center" valign="middle" >9.98</td><td align="center" valign="middle" >6.321</td></tr><tr><td align="center" valign="middle" >VSI (p.u)</td><td align="center" valign="middle" >0.1407</td><td align="center" valign="middle" >0.1482</td><td align="center" valign="middle" >0.1127</td></tr><tr><td align="center" valign="middle" >Execution time (s)</td><td align="center" valign="middle" >1.4724</td><td align="center" valign="middle" >1.6245</td><td align="center" valign="middle" >1.5348</td></tr></tbody></table></table-wrap><table-wrap id="table7" ><label><xref ref-type="table" rid="table7">Table 7</xref></label><caption><title> The optimal solution of single objective and multi objective functions using hybrid algorithm for IEEE 118 test bus system</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  rowspan="2"  >Objectives</th><th align="center" valign="middle"  colspan="3"  >Single objective</th><th align="center" valign="middle"  colspan="3"  >Two objectives</th><th align="center" valign="middle" >Three objectives</th></tr></thead><tr><td align="center" valign="middle" >Cost</td><td align="center" valign="middle" >Loss</td><td align="center" valign="middle" >VSI</td><td align="center" valign="middle" >Cost-loss</td><td align="center" valign="middle" >Cost-VSI</td><td align="center" valign="middle" >Loss-VSI</td><td align="center" valign="middle" >Cost-loss-VSI</td></tr><tr><td align="center" valign="middle" >Cost ($/hr)</td><td align="center" valign="middle" >128945</td><td align="center" valign="middle" >152014</td><td align="center" valign="middle" >149156</td><td align="center" valign="middle" >131992</td><td align="center" valign="middle" >131129</td><td align="center" valign="middle" ></td><td align="center" valign="middle" >131119</td></tr><tr><td align="center" valign="middle" >Loss (MW)</td><td align="center" valign="middle" >139.5</td><td align="center" valign="middle" >123.3</td><td align="center" valign="middle" >138.6</td><td align="center" valign="middle" >115.4</td><td align="center" valign="middle" ></td><td align="center" valign="middle" >125.12</td><td align="center" valign="middle" >143.5</td></tr><tr><td align="center" valign="middle" >VSI (p.u.)</td><td align="center" valign="middle" >0.47</td><td align="center" valign="middle" >0.479</td><td align="center" valign="middle" >0.412</td><td align="center" valign="middle" ></td><td align="center" valign="middle" >0.431</td><td align="center" valign="middle" >0.429</td><td align="center" valign="middle" >0.442</td></tr><tr><td align="center" valign="middle" >Execution time (s)</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td><td align="center" valign="middle" >35.487</td><td align="center" valign="middle" >35.897</td><td align="center" valign="middle" >35.127</td><td align="center" valign="middle" >48.087</td></tr></tbody></table></table-wrap><p>algorithm for IEEE 118 test bus system. It indicates the values of cost, loss and VSI of the single objective optimization are less when compared to the multiobjective optimization. It states that while considering more objective function, the optimized value and execution time are increased concurrently. Because of nonlinearity in large system the front is less consistent comparing to the smaller system.</p></sec><sec id="s5"><title>5. Conclusion</title><p>In this paper optimal power flow problem has three objectives with constraints, and it works out in the multiobjective optimization approach. This multiobjective problem is solved by single objective, two-objective optimization such as cost-loss, cost-voltage stability index, loss-voltage stability index, and three-objective optimization of cost-loss- voltage stability index using PSO, ACO and hybrid techniques. Hybrid swarm intelligence is used to avoid the stagnation of global and local search and get the best solution by fuzzy approach. The simulation results show that the hybrid swarm optimization provides better results and faster convergence compared to the ant colony optimization and particle swarm optimization. The proposed technique gives better results with less convergence time and it will help to choose the sensible solution to the multiobjective optimization problem in IEEE 30 &amp; IEEE 118 test bus system.</p></sec><sec id="s6"><title>Cite this paper</title><p>Rajalashmi, K. and Prabha, S.U. (2016) Hybrid Swarm Algori- thm for Multiobjective Optimal Power Flow Problem. Circuits and Systems, 7, 3589- 3603. http://dx.doi.org/10.4236/cs.2016.711304</p></sec></body><back><ref-list><title>References</title><ref id="scirp.70453-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Burchet, R.C., Happ, H.H. and Vierath, D.R. (1984) Quadratically Convergent Optimal Power Flow. IEEE Transactions on Power Apparatus and Systems, 103, 3267-3276.http://dx.doi.org/10.1109/TPAS.1984.318568</mixed-citation></ref><ref id="scirp.70453-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Capitanescu, F., Glavic, M., Ernst, D. and Wehenkel, L. (2007) Interior-Point Based Algorithms for the Solution of Optimal Power Flow Problems. Electric Power Systems Research, 77, 508-517. http://dx.doi.org/10.1016/j.epsr.2006.05.003</mixed-citation></ref><ref id="scirp.70453-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Olofsson, M., Anderson, G. and Soder, L. (1995) Linear Programming Based Optimal Power Flow Using Second Order Sensitivities. IEEE Transactions on Power Systems, 10, 1691- 1697. http://dx.doi.org/10.1109/59.466472</mixed-citation></ref><ref id="scirp.70453-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Salgado, R., Brameller, A. and Aitchison, P. (1990) Optimal Power Flow Solutions Using the Gradient Projection Method Part 2: Modelling of the Power System Equations, Generation, Transmission and Distribution. IEE Proceedings, 137, 429-435.</mixed-citation></ref><ref id="scirp.70453-ref5"><label>5</label><mixed-citation publication-type="other" xlink:type="simple">Sun, D.I., Ashley, B., Brewer, B., Hughes, A. and Tinney, W.F. (1984) Optimal Power Flow by Newton Approach. IEEE Transaction on Power Apparatus and Systems, PAS-103, 2864- 2880. http://dx.doi.org/10.1109/TPAS.1984.318284</mixed-citation></ref><ref id="scirp.70453-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Abido, M.A. (2002) Optimal Power Flow Using Particle Swarm Optimization. Electrical Power &amp; Energy System, 24, 563-571. http://dx.doi.org/10.1016/S0142-0615(01)00067-9</mixed-citation></ref><ref id="scirp.70453-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Todorovski, M. and Rajicic, D. (2006) An Initialization Procedure in Solving Optimal Power Flow by Genetic Algorithm. IEEE Transactions on Power Systems, 21, 480-487.http://dx.doi.org/10.1109/TPWRS.2006.873120</mixed-citation></ref><ref id="scirp.70453-ref8"><label>8</label><mixed-citation publication-type="other" xlink:type="simple">Bakirtzis, A.G., Biskas, P.N., Zoumas, C.E. and Petridis, V. (2002) Optimal Power Flow by Enhanced Genetic Algorithm. IEEE Transactions on Power Systems, 17, 229-236.http://dx.doi.org/10.1109/TPWRS.2002.1007886</mixed-citation></ref><ref id="scirp.70453-ref9"><label>9</label><mixed-citation publication-type="other" xlink:type="simple">Yuryevich, J. and Wong, K.P. (1999) Evolutionary Programming Based Optimal Power Flow Algorithm. IEEE Transactions on Power Systems, 14, 1245-1250.http://dx.doi.org/10.1109/59.801880</mixed-citation></ref><ref id="scirp.70453-ref10"><label>10</label><mixed-citation publication-type="other" xlink:type="simple">Abido, M.A. (2002) Optimal Power Flow Using Tabu Search Algorithm. Electric Power Components and Systems, 30, 469-483. http://dx.doi.org/10.1080/15325000252888425</mixed-citation></ref><ref id="scirp.70453-ref11"><label>11</label><mixed-citation publication-type="other" xlink:type="simple">Varadarajan, M. and Swarup, K.S. (2008) Solving Multi-Objective Optimal Power Flow Using Differential Evolution. IET Generation, Transmission &amp; Distribution, 2, 720-730.http://dx.doi.org/10.1049/iet-gtd:20070457</mixed-citation></ref><ref id="scirp.70453-ref12"><label>12</label><mixed-citation publication-type="other" xlink:type="simple">Abou El Elaa, A.A., Abidob, M.A. and Speaa, S.R. (2010) Optimal Power Flow Using Differential Evolution Algorithm. Electric Power Systems Research, 80, 878-885.http://dx.doi.org/10.1016/j.epsr.2009.12.018</mixed-citation></ref><ref id="scirp.70453-ref13"><label>13</label><mixed-citation publication-type="other" xlink:type="simple">Kumari, M.S. and Maheswarapu, S. (2010) Enhanced Genetic Algorithm Based Computation Technique for Multi-Objective Optimal Power Flow Solution. Electrical Power and Energy Systems, 32, 736-742. http://dx.doi.org/10.1016/j.ijepes.2010.01.010</mixed-citation></ref><ref id="scirp.70453-ref14"><label>14</label><mixed-citation publication-type="other" xlink:type="simple">Shi, L.B. and Xu, G.Y. (2001) Self-Adaptive Evolutionary Programming and Its Application to Multi-Objective Optimal Operation of Power Systems. Electric Power Systems Research, 57, 181-187. http://dx.doi.org/10.1016/S0378-7796(01)00086-4</mixed-citation></ref><ref id="scirp.70453-ref15"><label>15</label><mixed-citation publication-type="other" xlink:type="simple">Niknam, T., Narimani, M.R., Aghaei, J. and Azizipanah-Abarghooee, R. (2012) Improved Particle Swarm Optimisation for Multi-Objective Optimal Power Flow Considering the Cost, Loss, Emission and Voltage Stability Index. IET Generation, Transmission &amp; Distribution, 6, 515-527. http://dx.doi.org/10.1049/iet-gtd.2011.0851</mixed-citation></ref><ref id="scirp.70453-ref16"><label>16</label><mixed-citation publication-type="other" xlink:type="simple">Sivasubramani, S. and Swarup, K.S. (2011) Multi-Objective Harmony Search Algorithm for Optimal Power Flow Problem. Electrical Power and Energy Systems, 33, 745-752.http://dx.doi.org/10.1016/j.ijepes.2010.12.031</mixed-citation></ref><ref id="scirp.70453-ref17"><label>17</label><mixed-citation publication-type="other" xlink:type="simple">Fonseca, C.M. and Fleming, P.J. (1995) An Overview of Evolutionary Algorithms in Multiobjective Optimization. Evolutionary Computation, 3, 1-16. http://dx.doi.org/10.1162/evco.1995.3.1.1</mixed-citation></ref><ref id="scirp.70453-ref18"><label>18</label><mixed-citation publication-type="other" xlink:type="simple">Rosehart, W.D., Ca?izares, C.A. and Quintana, V.H. (2003) Multiobjective Optimal Power Flows to Evaluate Voltage Security Costs in Power Networks. IEEE Transactions on Power Systems, 18, 578-587. http://dx.doi.org/10.1109/TPWRS.2003.810895 </mixed-citation></ref><ref id="scirp.70453-ref19"><label>19</label><mixed-citation publication-type="other" xlink:type="simple">Shuang, B., Chen, J.P. and Li, Z.B. (2011) Study on Hybrid PS-ACO Algorithm. Application Intelligence, 34, 64-73. http://dx.doi.org/10.1007/s10489-009-0179-6</mixed-citation></ref></ref-list></back></article>