Improving Profile Parameters of the Power System Network Using Krill Herd Algorithm with Facts Device: UPFC # Sunil Kumar Jilledi Abstract-In the electric power system the optimal power flow plays a key role in optimal location, steady state operation point, minimization of generation cost, losses etc. Large Power transmission networks are facing lot of problems to improve the power profiles, improvement of losses, security constraints, and stability issues. The advent of Flexible AC Transmission System (FACTS) devices has led the controlling capabilities to improve the power profiles. Out of these Unified Power Flow Controller (UPFC) has the capability to improve all the three system variables namely line reactance, magnitude and phase angle difference of voltage across the line simultaneously or individually. Many natural inspired algorithms have implemented to solve the optimal power flow with incorporation of UPFC. The main scope of this paper is to implement the KRILL HERD ALGORITHM (KHA) a novel meta-heuristic approach which is influenced from the herding actions of the krill swarms search for food or communication with each other. Considering the generator buses play down the real power losses, generator bus voltages and reactive power injection. The proposed KRILL HERD ALGORITHM based optimal power flow has been tested for the IEEE 14-30 bus systems. The results have been presented clearly for the proposed algorithm with and without incorporation of UPFC. Proposed algorithm results have been compared with GA, PSO, FF and ABC algorithms. Keywords: krill herd algorithm (kha), UPFC, factsoptimal power flow. Nomenclature # Introduction he optimal power flow problem is becoming a peculiar topic in the power systems. Due to increasing of load demand the power systems are becoming large by interconnecting with different regional systems. Interconnected systems are facing more failures [1]. It is becoming a tedious task for the power system engineers to utilize the existing transmission lines efficiently. Optimal power solution is the best process to get better output with the existing systems, by generation relocation. For efficient utilization of the existing system the shunt capacitors and shunt reactors are incorporating to improve the voltage profile and transmission line reactance as well as power transfer capability. To improve the phase shift between receiving and sending voltages phase shifting transformers are using. Moreover the faster expansion and interconnection of the regional systems voltage stability, power system securities are facing in the deregulated market. In the literature different authors described about the voltage collapse [2,3]. The power electronic devices are playing a key role in the recent era. The advanced development in the power electronics controllers leads to develop the Flexible AC Transmission System (FACTS) to supply flexible power in the system. Optimize the utilization of the existing system by incorporating the FACTS devices. The FACTS devices technology was presented by Electric Power Research Institute (EPRI) in the year 1980s. These devices has the capability to control the different parameters of the transmission line such as shunt/ series impedance, phase angle, real and reactive power compensation, etc The FACTS family include number of devices such as Stastic VAR Compensator(SVC), thyristor controlled reactor (TCRs), are Shunt FACTS devices, later the series FACTS devices [4].UPFC powerful FACTS device, combination of Static Synchronous Compensator(STATCOM) and Static Synchronous Series Compensator(SSSC) coupled by DC link [5]. Optimal power flow problem is solved by adjusting several variables in the objective function considering generations cost, loss function.etc. Over the decades many researchers presented different solutions to optimal power flow by using different methods Newton Method, Genetic algorithm [6], Differential Evolution and Evolutionary programming, BAT Algorithm [7,8]. Researchers are showing interest on meta-heuristic techniques which includes Genetic Algorithm (GA), Practical Swarm optimization(PSO),Ant colony Algorithm. In the literature authors [9] proposed optimal power flow using GA other [10], gravitational search algorithm (GSA) [11], artificial bee colony (ABC) optimization [12] using swarm intelligence for the optimal power flow. Researchers proposed these algorithms to overcome the failures of the conventional implemented to solve the optimal power flow problem. In this paper the Krill Herd (KH) a Meta heuristic algorithm is proposed it is one of the bio-based swarm intelligence algorithms. The Krill Herd is developed based on the behaviour of Krill Swarms [13] i.e., distance between food and highest density of swarms simulates the objective function of individual krill. Comparing with other optimization techniques in the KH the controlling variables are very few. The Krill Herd already using in some research areas like optimization problems [14]. This paper solves the optimal power flow without and with UPFC using the Krill Herd algorithm for different IEEE standard bus systems. The main objective function considered is minimization of the real power losses, voltage deviation, incorporation of UPFC is considered based on the real power losses. The results obtained are presented clearly. KH algorithm optimal power flow results with UPFC is compared with GA and BAT algorithm. The paper organization is follows in the coming section about the KHA, Formulation of UPFC model, optimal power flow using conventional method and proposed method using KHA. Problem solving using the Matlab simulation results and discussion finally the conclusion of the paper and the future work II. # Power Flow Model of UPFC Gyugyi proposed the Unified Power Flow Controller (UPFC), for real time control and dynamic compensation of AC transmission systems. UPFC consists of Static Synchronous series Compensator (SSSC) and STATCOM connected by a DC link capacitor. UPFC is capable to control the active and reactive power and voltage magnitudes simultaneously at the terminals of UPFC [15]. UPFC consists of two converters, Converter 2 controls the power flow of the device by infuse of an AC voltage V_pq in controllable magnitude and phase angle in series to the transmission line. Similarly the converter 1 can absorb or supply the real power demand by the converter 2 at the DC link. Each converter can supply or absorb the real and reactive power demanded by the system independently [16]. Finding the load flows of any power system is the initial stage to evaluate the power system. Many iterative solutions are there for finding the load flow like Gauss, Newton Raphson method, decouple, fast decouple, Ranga-Kutta methods are available. In this paper the load flows are performed by Newton Rahson Method by using the polar coordinates.Fig. 01 shows the clear model of UPFC connected between the bus i and j. and power flow directions of real and reactive power at the shunt and series elements where UPFC is connected. For each bus the real and reactive powers are computed by Eqs ( 1) and ( 2) ?? ?? = ? ?? ?? ?? ?? (?? ???? ??ð??"ð??"?? ?? ???? + ?? ???? ?????? ?? ???? ) ?? ?? =1 (1) Q i = ? V i V j (G ij sin ? ij ? B ij cos ? ij ) N j=1 (2) After finding the load flows by the conventional method, the UPFC is incorporated in the system, compute the power flows with UPFC. The UPFC voltage sources are given in Eqs. ( 3) and ( 4) ?? ??? = ?? ??? (??ð??"ð??"?? ?? ??? + ?????? ?? ???(3)?? ???? = ?? ???? (cos ?? ???? + sin ?? ???? )(4) The active and reactive power equations are given in Eqs. ( 5) -( 8 At bus -j ?? ?? = [?? ?? ?? ?? ?? ???? sin??? ?? ? ?? ?? ? + ?? ?? ?? ???? ?? ???? sin??? ?? ? ?? ???? ?](7)?? ?? = ??? ?? 2 ?? ???? ? ??? ?? ?? ?? ?? ???? cos??? ?? ? ?? ?? ?? ? [ ?? ?? ?? ???? ?? ???? cos??? ?? ? ?? ???? ?](8) Power flow equations at the converter terminals of UPFC Eqs. ( 9)-( 12) At the series converter ?? ???? = [?? ???? ?? ?? ?? ???? sin(?? ???? ? ?? ?? )] + [?? ?? ?? ??? ?? ???? sin??? ???? ? ?? ?? ?](9)?? ???? = ??? ???? 2 ?? ???? ? ??? ?? ?? ???? ?? ???? cos??? ?? ? ?? ?? ?? ? [ ?? ?? ?? ???? ?? ???? cos??? ?? ? ?? ???? ?] (10) At the shunt converter ?? ??? = ?? ??? ?? ?? ?? ??? sin(?? ???? ? ?? ?? )(11)?? ??? = ?? ??? 2 ?? ??? ? [?? ?? ?? ??? ?? ??? cos(?? ???? ? ?? ?? )](12) For the analysis in this paper the source reactance are considered as X_sr=X_sh=0.1 p.u. The UPFC Source voltage and phase angles are considered as V_sr=0.02 p.u,V_sh=1 p.u,?_sr=85 ?_sh=0. When UPFC is connected between bus-i and j in the power system ? ?? ?? ?? ?? ? = ? ?? ???? + ?? ??? ??? ???? ??? ???? ??? ??? ??? ???? ?? ???? ?? ???? 0 ? ? ?? ?? ?? ?? ?? ???? ?? ??? ? (13) Where ?? ???? = The imaginary distances between the krill herd and food give the best fitness value .The main two characteristics considered in the engineering optimization problems are exploration and random search are needed for better performance. The main objective function of the KHA is from the Lagrangian model [17][18][19]. In the two dimensional problems the above mentioned factors are sufficient, for the ndimensional problem analysis for the i^th krill individual is given by ???? ?? ???? = ?? ?? + ?? ?? + ?? ??(14) Motion Induced by other Krill Individuals: The fitness function mainly depends on the concreteness of the krill's in the searching space. The main significant and essential thing to obtain the optimum solution is to maintain the krill density. The motion of the individual krill is promptly dependent on the adjacent individual krill and the effects between them. The direction of individual krill movement is designed on different swarm densities [19]. a) Local effect provided by local krill density b) Target effect provided by target krill density c) Repulsive effect provided by repulsive swarm density ?? ?? ?????? = ?? ?????? ?? ?? + ð??"ð??" ?? ?? ?? ð??"ð??"???? (15) Where ?? ?? = ?? ?? ??ð??"ð??"?????? + ?? ?? # ??????ð??"ð??"???? The effect of krill individual on the nearest krill is calculated by ?? ?? ??ð??"ð??"?????? = ? ?? ???? ?? ???? ???? ?? =1(16) Where ?? ???? = ?? ?? ??? ?? ???? ?? ??? ?? ??+?? ; ?? ???? = ?? ?? ??? ?? ?? ??ð??"ð??"?????? ??? ???????? ; To know the distance between each individual is given by ?? ???? = ? ||?? ?? ? ?? ?? || ???? ?? =1(17) Foraging activity: Foraging activity is computed based on two main factors, First factor is current food location, second is prior food location information. The foraging velocity of the ith krill individual is given by the formula (18) [13] ?? ?? = ?? ð??"ð??" ?? ?? + ð??"ð??" ð??"ð??" ?? ?? ð??"ð??"???? (18) Where ?? ?? = ?? ?? ð??"ð??"ð??"ð??"ð??"ð??"?? + ?? ?? ???????? Food attraction is calculated by Eqs. ( 19) ?? ?? ð??"ð??"ð??"ð??"ð??"ð??"?? = ?? ð??"ð??"ð??"ð??"ð??"ð??"?? ?? ??,ð??"ð??"ð??"ð??"ð??"ð??"?? ?? ??,ð??"ð??"ð??"ð??"ð??"ð??"??(19) Where ?? ð??"ð??"ð??"ð??"ð??"ð??"?? = 2(1 ? ?? ?? ?????? ) # Physical Diffusion: In the diffusion process mainly considered to increase density of population. This motion is a based on maximum diffusion speed and random directional vector, It is given by Eqs. (20) [13] ?? ?? = ?? ?????? ?? (20) Motion process in KHA Depending up on the local effect, global effect, presence of food, best fitness position, the presence of the i^th krill stays in the time interval [t,t+Î?"t] given by Eqs. ( ???? ?? = ??(??) + ??(??) + ??(??)(22) The scaling factor Î?"t is formulated in Eqs.( 23) Î?"?? = ?? ?? ? (???? ?? ? ???? ?? ) ???? ?? =1(23) Step In equality constrains are Generator active and reactive powers, voltage magnitudes, Transformer tap settings, UPFC settings. The limits of the generator real and reactive powers limits at P^th bus should lie between maximum and minimum limits Eqs. (26)(27). The voltage magnitude at the each load bus is given in Eqs. (28). Transformer tap setting minimum and maximum conditions is given in Eqs. (29). Transmission line loading should not violate the loading limits is Eqs (30) [20][21][22][23][24][25] ?? ???? ?????? ? ?? ???? ? ?? ???? ?????? ????????? ?? = 1,2,3, ? ? ?? ????(26) ?? ???? ?????? ? ?? ???? ? ?? ???? ?????? ????????? ?? = 1,2,3, ? ? ?? ???? (27) where ?? ???? ?????? , ?? ???? ?????? , ?? ???? ?????? , ?? ???? ?????? are minimum and maximum limits of real and reactive powers at ?? ??? bus. ?? ???? ?????? ? ?? ???? ? ?? ???? ?????? ????????? ?? = 1,2,3, ? ? ?? ???? (28) where ?? ???? ?????? , ?? ???? ?????? are minimum and maximum limits of voltage at ?? ??? bus . ?? ?? ?????? ? ?? ?? ? ?? ?? ?????? ????????? ?? = 1,2,3, ? ? ?? ??(29) where ?? ?? ?????? , ?? ?? ?????? are minimum and maximum tap setting limits of transformer . ?? ???? ? ?? ???? ?????? ????????? ?? = 1,2,3, ? ? ?? ??(30) Where ?? ???? , ?? ???? ?????? are the total power flow in the ?? ??? ???????????. # Objective function The objective functions considered in this article are based on the fuel cost [25] ???????? ??ð??"ð??"???? (?? ?? ) = ? (?? ?? + ?? ?? ?? ???? + ?? ?? ?? ???? 2 ) ?? ???? ??=1(31) For minimization of transmission losses, the mathematical formula is given as Min P Loss = ? ?? ?? ?? ?? ??=1 [?? ?? 2 + ?? ?? 2 ? 2|?? ?? ||?? ?? |cos(?? ?? ? ?? ?? )](32) Line identification is very essential to locate the UPFC in the proposed system. Optimal location of UPFC is calculated by using the Performance Index (PI) given Eqs.( 33) By satisfying the equality and inequality constants minimize the objective function is the main objective of optimal power flow (OPF). OPF is used for the incorporation of UPFC in the system, considering four different objective functions by fulfilling equality and inequality constraints. PI = ?? ?? 2?? ( ?? ?? ??? ?? ?? ??? ?????? ) 2??(33) The general optimization problem constraints are as follows Objective function to be minimised is Min(u, v), and subjected to g(u, v)=0; h(u, v)?0, the g(u, v) is equality constraints, h(u, v) is inequality constraints, u is dependent variable , v is independent variable. The dependent variables considered in the problem formulation are generator active power (slack bus) P_G1, load voltages (V_(L1,) V_L2, ??.V_(?LN_PQ )), reactive power at the generators (Q_(G1,) Q_G2,?? Q_(GN?_PV)), Line loading of the transmission system (S_L1,S_L2,??.S_(LN?_L)). The independent variables are active power of the generators apart from slack bus (P_G1,P_G2??.P_(?GN?_PV )), generator voltages (V_(G1,) V_G2, ??.V_(?GN?_PV)), Transformer tap settings (T_(1,) T_2,??.. T_(N_T)), active power injections (P_c1,P_c2,??.P_(CN_U)), reactive power injection (Q_c1,Q_c2,??.Q_(CN_U)) [20][21][22][23][24][25]. # Equality and Inequality constraints: Mentioned above g is the set of equality constraint and h is inequality constraint. With the help of load flow equations the equality constraints are represented by Eqs. (24)(25). The inequality constraint h is the operating limits represented by Eqs. (26)(27)(28)(29)(30) Voltage Deviation should be very minimum at all the bus formulated as Eqs.( 34) ?? ?????? = min(??????) = min ? |?? ?? ? ?? ?? ????ð??"ð??" | ?? ??=1(34) # V. Simulation Results and Discussion For better understanding analysis of the proposed KHA is simulated by using IEEE 14 and 30 bus standard systems. At the initial state IEEE 14 and 30 bus system load flows are run by Newton Raphson Method using the polar coordinates in the MATLAB environment. IEEE 14 bus system is included by 5 generation units which are located at the Bus No. 1, 2, 3, 6, 8 and 20 transmission lines are used to interconnect the system and tap changing transformers are connected between the buses(4-7,4-9 and 5-6) and for the Bus-9 and 14 shunt VAR compensators are connected . The total demand by the system is 2.98p.u. at 100MVA base. Control variables and line data is considered [26] The data of the Modified IEEE 30 bus system is having six generators located at the buses -1, 2, 5, 8, 11, 13 and remaining 24 are the load buses, 41 transmission lines are used to interconnect the system. The slack bus is considered as bus -1. Total demand by IEEE 30 bus system is 2.83 p.u at 100 MVA base. In the system load bus, voltages are considered in the range of 0.95 to 1.1p.u. IEEE-14 bus system minimum and maximum constraints is shown in Table .01 [26] [28], GA [30], ABC [31], PSO [33]. The total power generated by using KHA is decreased by 1.4% as compared with ABC with the incorporation of UPFC between bus 24 and 25. The results of ABC are mentioned in [31]. Similarly there is a decrease of 1.06% for GA is reported in [30], and 1.248% decrease for FF reported in [28]. By using the fuel cost optimization for the KHA method, the power losses have reduced to 4.6986% as compared with the other optimization techniques. # Conclusion and Future Scope A novel Meta heuristic algorithm KHA is used to solve the Optimal power flow problem of the proposed power system networks IEEE-14 and 30 bus systems. Two main objective functions has considered (i) cost function (ii) Active power transmission losses due to high impact of equality and inequality constraint each objective function is studied individually. For the analysis of the KHA, FACTS device UPFC is incorporated in the system. Results obtained using KHA is compared with Genetic Algorithm, Practical Swarm Algorithm, Fire Fly and ABC algorithms and compared with the other popular optimization techniques for the optimal power problem, The results obtained from the KHA are better and robustness, stability and the convergence rate is faster than the other methods. By this article the new algorithm KHA may be extended for other optimization methods for the further research. In future the KHA can be extended to OKHA, and can be implemented for the other FACTS devices like IPFC, UPQC etc,. for the better analysis. ![?? ???? = Suspectance between ?? ??? ?????? ?????? ?? ??? ?????? ?? ð??"ð??" = Empirical constant [0, 2] ?? ð??"ð??"ð??"ð??"ð??"ð??"?? = Coefficient of Effective food ?? ?????? =Maximum diffusion speed ?? ?? = Random Diffusion ?? ?? ð??"ð??"???? = Last foraging motion ?? ?? =Foraging motion ?? ???? = Conductance between ?? ??? ?????? ?????? ?? ??? ?????? ?? ?? = current at the ?? ??? ?????? ?? ?????? = Total number of Iterations ?? ???????? and ?? ????ð??"ð??"???? = Best and worst fitness of each individual ?? ?? ?????? ?? ?? = Fitness value of the ?? ??? ?????? ?? ??? krill individual NN = Total number of Neighbours NV = Number of Variables ?? ???? = Number of generator buses ?? ???? = Number of Load buses ?? ?? = Number of Transmission lines ?? ?? = Number of Tap setting Transformers n = exponent taken as"1" N = Number of buses ?? ???????? = bus where UPFC is connected ?? ?? ð??"ð??"???? = Motion induced by the previous Krill ?? ?? =Motion induced on the ?? ??? krill individual depending on the other krill individual ?? ?????? =Maximum induced speed ?? ???? = Total power generated at ?? ??? ?????? ?? ?? = Total power demand ?? ?? = Total power Losses ?? ???? ?????? ?????? ?? ???? ?????? = Minimum and maximum real power generated at the ?? ??? ð??"ð??"????????????ð??"ð??"?? ?? ???? ?????? ?? ???? = real and reactive powers at ?? ??? bus ?? ???? ?????? ?? ???? = real and reactive powers injected by the UPFC at ?? ??? bus ?? ????? = Apparent power in the line connected between buses P and Q bus ???? ?? ?????? ???? ?? = Upper and Lower boundaries of the variables ?? ?? = ?? ??? bus voltage ?? ??? = Controllable voltage at the shunt converter ?? ???? = Controllable voltage at the series converter ?? ???? and ?? ??? = Admittance at the series and shunt converter ?? ???? and ?? ??? = Impedance at the series and shunt converter X= Relative position of each Krill ?? ??? = Phase angle of voltage source at the shunt converter ?? ???? = Phase angle of the voltage source at the series converter ð??"ð??" ?? =Weight of Inertia ?? ?? ??ð??"ð??"?????? = Local effect provided by neighbouring krill ?? ?? ??????ð??"ð??"???? = Target effect provided by individual the best krill individual ? = Small positive number ?? ð??"ð??" = Foraging speed ð??"ð??" ð??"ð??" = Inertia motion of foraging speed ?? ?? ð??"ð??"ð??"ð??"ð??"ð??"?? ?????? ?? ?? ???????? = Effect due to presence of food and Effect due to current Krill's best fitness value recorded ?? = Random directional vector ?? ???? = Admittance angle of the transmission line connected between p-bus and Q -bus ð??"ð??" ð??"ð??" = real non negative weighing coefficient](image-2.png "") 01![Fig. 01: UPFC voltage source model connected between i^(th ) and j^(th ) bus](image-3.png "Fig. 01 :") ![Krill Herd Algorithm (KHA) proposed by the researchers Gandomi and Alavi in 2012. KHA is a metaheuristic algorithm enthused by bio-based swarm intelligence algorithm. KHA is simulated based on the behaviour of the Krill Swarms. Mostly based on the food of the highest density of the Swarms forming the objective function of each Krill folk. The position of each Krill folk is dependent on following factors [17]: a) Movement induced by other Krill folk b) Foraging Activity c) Random Diffusion](image-4.png "") ) At bus-i ???? ?? ????(21) ? ??? ?? ??=1 ?? ?? ??=1?? ???? ? ?? ???? ?? ???? ? ?? ????+ ? + ? ?? ???????? ??=1 ?? ???????? ??=1?? ???? ?? ????= ? = ? ? ?? ?? ??=1 ?? ?? ? ??=1 ?? ?? ??=1 ? ?? ?? |?? ?? | |?? ?? | ??=1(24)( ) Volume XVII Issue III Version IJournal of Researches in EngineeringGlobal??? ??? Generating Unit?? ???????? (p.u)?? ???????? (p.u)Pg10.000.1Pg2-0.40.5Pg30.000.4Pg6-0.060.24Pg8-0.060.24Voltage Limits Transformer tap changer Line voltage?? ?? ?????? = 0.95 ?? ?????? = 0.9 ?? ?? ?????? = 0.95?????? = 1.05 ?? ?? ?? ?????? = 1.1 ?? ?? ?????? = 1.05 Fig. 04:Global Journal of Researches in Engineering ( ) Volume XVII Issue III Version ITable-02: Cost constraints and maximum and minimum power limits of the generator units [27] Generating unit P(min) P (Max) ?? ?? ?? ?? * ???? ??? ?? ?? * ???? ??? MW MW $/? $/????? $/???? 2 ? Pg1 50 200 0.00 200 Pg2 20 80 0.00 175 175.0 Pg5 15 50 0.00 100 625.0 Pg8 10 35 0.00 325 83.0 Pg11 10 30 0.00 300 250.0 Pg13 12 40 0.00 300 250.0 Case study-i: smoothly reduced as compared with the GA and PSO. From the convergence results, it is clearly observed that reduced by 80% with respective to other algorithms. 12 13 14 15 16 17 18 19 20 Active power Transmission Losses(p.u) Genetic Algorithm PSO Krill Herd Algorithm 37.5 Total APTL are presented clearly in the Fig.03, from which it is clearly observed that the APTL are APTL are reduced by 0.8p.u in contrast with GA and PSO. By implementing the KHA almost APTL areGlobal Journal of Researches in Engineering ( ) Volume XVII Issue III Version I1104080Number of Fitness FunctionFig. 03: Comparative Active power transmission losses (APTL) of GA,PSO and Krill Herd for IEEE 14 bus sytem. © 2017 Global Journals Inc. (US) Generator ABCGAFF OPFPSOKHA OPFPG1180.5218176.7307176.7311174.26174.16PG248.784548.848848.845449.7748.754PG521.259821.494121.82PG818.646921.688121.692321.420.61PG1111.814512.153012.153511.9311.95PG1312.101112.000912.00012.0012.01TOTAL293.3805292.3805292.9154290.41289.304COST802.1649802.717802.3646802.36795.41?? ????????9.72869.51569.51559.30649.292 Generator PG1PG2PG5PG8PG11PG13TOTALCOST?? ????????KHA OPF74.35663.50949.99933.43128.7734.251284.316952.562.753 Year 2017 F © 2017 Global Journals Inc. (US) © 2017 Global Journals Inc. (US) Year 2017 F Improving Profile parameters of the Power System Network using Krill Herd Algorithm with FACTS device: UPFC Year 2017 F © 2017 Global Journals Inc. (US)The IEEE 30 bus system active power generating limits and unit cost of generators are presented in table-2.Improving Profile parameters of the Power System Network using Krill Herd Algorithm with FACTS device: UPFC Year 2017 F Improving Profile parameters of the Power System Network using Krill Herd Algorithm with FACTS device: UPFC Case study-ii IEEE 30 bus system is considered for the enhanced analysis. Voltage profiles, real power at generating units, APTL, Cost analysis is evaluated in the MATLAB environment. The bus system is simulated with and without incorporation of UPFC. Based on TVD the UPFC is incorporated. The UPFC is installed between bus 24 and 25. The bus data, line data, generation data are considered from [28]. The voltage profiles of the IEEE 30 bus system is obtained by simulating in MATLAB,and the results are presented in Table .05 is incorporating the UPFC in line 33 between and 25 buses. As compared with NR, GA, FF, ABC [28][29][30][31][32] with KHA the voltage profiles are smoothly and drastically increased in the system. Voltage profile is improved almost 0.06% compared with the conventional and remaining algorithms. * On computational issues of market based optimal power flow WHongye CEMurillo Sanchez RZimmerman Thomasr IEEE Trans.Power systems 3 22 2007 * Power system stability and control PKundur 1994 McGraw-Hill * Power system voltage stability CWTaylor 1994 McGraw-Hill Inc * Power flow control and solutions with multiple and multi-type FACTS devices MANarayana Prasad Padhy Moamen Electric Power Systems Research 74 2005 * Improvement of transient stability using unified power flow controller power delivery RMihalic PZunko DPovh IEEE Transactions 11 1996 * Optimal power flow for large-scale power system with shunt FACTS using efficient parallel GA BMahdad KSrairi TBouktir Int J Electr Power Energy Syst 32 5 2010 * Transmission loss minimization using SVC based on particle swarm optimization JumaatSiti Amely MusirinIsmail MokhlisHazlie 2011 IEEE symposium on industrial electronics and applications (ISIEA2011) Langkawi, Malaysia * Optimal VAR control for real power loss minimization using differential evolution algorithm AmraneYoucef BoudourMohamed LadjiciAhmed Amine ElmaouhabAli Int J Electr Power Energy Syst 66 2015 * Genetic algorithm for solving optimal power flow problem with UPFC VKrishnasamy Int. J. Softw. Eng. Appl 5 1 2011 * Loss minimization using optimal power flow based on swarm intelligences NSinsuphun ULeeton UKwannetr DUthitsunthorn TKul-Worawanichpong ECTI Trans. Electr. Eng. Electr. Commun 9 1 2011 * Optimal reactive power dispatchusing a gravitational search algorithm SDuman YSonmez UGuvenc NYorukeren IET Gener Transm Distrib 6 6 2012 * Artificial bee colony algorithm solution for optimal reactive power flow KAyan UKilic Appl Soft Comput 12 5 2012 * Krill Herd: A new bio inspired optimization algorithm AHGandomi AHAlavi 2012 17 * Application of the bio-inspired krill herd optimization technique to phase equilibrium calculations KMoodley JRarey DRamjugernath Comput Chem Eng 74 2015 * Steady-state and Dynamic Models of Unified Power Flow Controller (UPFC) for Power System Studies ANabavi-Niaki MRIravani IEEE Trans. Power Systems 11 4 1996 * NGHingorani & Naren Understanding FACTS * Lagrangian modelling studies of Antarctic krill (Euphasia superba) swarm formation EEHofmann AgeHaskell JMKlinck CMLascara 2004 * Krill herd algorithm for optimum design of truss structures AHGandomi STalatahari FTadbiri AHAlavi Int J Bio-Inspired Comput 5 5 2013 * Fuzzy krill herd optimization algorithm EFattahi MBidar HRKanan Proc.of IEEE 1st int. conf. on networks & soft computing (ICNSC) .of IEEE 1st int. conf. on networks & soft computing (ICNSC)Guntur Aug. 19-20. 2014 * Multi-objective optimal power flow with FACTS devices MBasu Energy Convers Manage 52 2 2011 * A hybrid GA approach for OPF with consideration of FACTS devices TSChung YKLi IEEE Power Eng Rev 20 8 2000 * Optimal power flow with FACTS devices by hybrid TS/SA approach WOngsakul PBhasaputra Int J Electr Power Energy Syst 24 10 2002 * Dynamic strategy based fast decomposed GA coordinated with FACTS devices to enhance the optimal power flow MBelkacem TBouktir KSrairi Energy Convers Manage 51 7 2010 * Optimal power flow with TCSC and TCPS modeling using craziness and turbulent crazy particle swarm optimization PKRoy SPGhoshal SSThakur Int J Swarm Intel Res, IGI Global Pub 1 3 2010 * Particle swarm optimization with an aging leader and challengers algorithm for optimal power flow problem with FACTS devices RPSingh VMukherjee SPGhoshal Int J Electr Power Energy Syst 64 2015 * Optimal power flow for steady state secu-rity enhancement using enhanced genetic algorithm with FACTS devices DDevaraj RNBanu Asian Power Electr. J 4 3 2010 * Optimal reactive power dispatch using self-adaptive real coded genetic algorithm PSubbaraj PNRajnarayanan Electr Power Syst Res 79 2 2009 * Obj ect-Oriented Optimal Power Flow TBouktir MBelkacemi Electric Power Components and systems 2003 31 * Optimal Power Flow Using Firefly Algorithm with Consideration of FACTS Devices "UPFC OuafaHerbadji TarekBouktir International Journal on Electrical Engineering and Informatics 7 1 March 2015 * Biogeography Based Optimization Approach for Solving Optimal Power Flow Problem OHerbadji LSlimani TBouktir International Journal of Hybrid Information Technology 6 5 2013 * A Genetic Algorithm for Solving the Optimal Power Flow Problem TBouktir LSlimani MBelkacemi Leonardo Journal of Sciences 2004 * Optimal Power Flow with Emission Controlled using Artificial Bee Colony Algorithm LSlimani TBouktir 12th International confer ence on Sciences and Techniques of Automatic control & computer engineering Sousse, Tunisia 2011 * Milkias Berhanu Tuka 2 , Dr. Sultan F. Meko 3 , Shalini J 4 , Dawit Leykuen 5, Line Losses in the 14-Bus Power System Network using UPFC SunilKumar J1 ACEEE Int. J. on Electrical and Power Engineering 5 1 February 2014 * Optimal Location of FACTS for Voltage stability using Modified Particle Swarm Optimization Whei-MinLin Member Kai-HungIeee Chih-MingLu Chia-ShengHong Tu Proceedings of the International Multi conference of Engineers and Computer II science32012 IMECS