# Introduction t the present time, The Optimal power flow (OPF) is a very significant problem and most focused objective for power system planning and operation [1]. The OPF is the elementary tool which permits the utilities to identify the economic operational and secure states in the system [2]. The OPF problem is one of the utmost operating desires of the electrical power system [3]. The prior function of OPF problem is to evaluate the optimum operational state for Bus system by minimizing each objective function within the limits of the operational constraints like equality constraints and inequality constraints [4]. Hence, the optimal power flow problem can be defined as an extremely non-linear and non-convex multimodal optimization problem [5]. From the past few years too many optimization techniques were used for the solution of the Optimal Power Flow (OPF) problem. Some traditional methods used to solve the proposed problem have some limitations like converging at local optima and so they are not suitable for binary or integer problems or to deal with the lack of convexity, differentiability, and continuity [6]. Hence, these techniques are not suitable for the actual OPF situation. All these limitations are overcome by metaheuristic optimization methods. Some of these methods are [7][8][9][10]: genetic algorithm (GA) [11], hybrid genetic algorithm (HGA) [12], enhanced genetic algorithm (EGA) [13][14], differential evolution algorithm (DEA) [15][16], artificial neural network (ANN) [17], particle swarm optimization algorithm (PSO) [18], tabu search algorithm (TSA) [19], gravitational search algorithm (GSA) [20], biogeography based optimization (BBO) [21], harmony search algorithm (HSA) [22], krill herd algorithm (KHA) [23], cuckoo search algorithm (CSA) [24], ant colony algorithm (ACO) [25], bat optimization algorithm (BOA) [26], Ant-lion optimizer (ALO) [27][28] and Multi-Verse optimizer (MVO) [29]. In the present work, a newly introduced hybrid meta-heuristic optimization technique named Hybrid Particle Swarm Optimization-Moth Flame Optimizer (HPSO-MFO) is applied to solve the Optimal Power Flow problem. HPSO-MFO comprises of best characteristic of both Particle Swarm Optimization [30] and Moth-Flame Optimizer [31][32] algorithm. The capabilities of HPSO-MFO are finding the global solution, fast convergence rate due to the use of roulette wheel selection, can handle continuous and discrete optimization problems. According to No Free Lunch Theorem [27,29,30], particular meta-heuristic algorithm is not best for every problem. So, we considered HPSO-MFO for continues optimal power flow problem based on No Free Lunch Theorem. In this work, the HPSO-MFO is presented to standard IEEE-30 bus test system [33] to solve the OPF [34-37] problem. There are five objective cases considered in this paper that have to be optimize using HPSO-MFO technique are Fuel Cost Reduction, Voltage Stability Improvement, Voltage Deviation Minimization, Active Power Loss Minimization and Reactive Power Loss Minimization. The results show the optimal adjustments of control variables in accordance with their limits. The results obtained using HPSO-MFO technique has been compared with Particle Swarm Optimisation (PSO) and Moth Flame Optimizer (MFO) techniques. The results show that HPSO-MFO gives better optimization values as compared to other methods which prove the effectiveness of the proposed algorithm. This paper is summarized as follow: After the first section of the introduction, the second section concentrates on concepts and key steps of standard PSO and MFO techniques and the proposed Hybrid PSO-MFO technique. The third section presents the formulation of Optimal Power Flow problem. Next, we apply HPSO-MFO to solve OPF problem on IEEE-30 bus system in order to optimize the operating conditions of the power system. Finally, the results and conclusion are drawn in the last section. # II. Standard PSO and Standard MFO a) Particle Swarm Optimization The particle swarm optimization algorithm (PSO) was discovered by James Kennedy and Russell C. Eberhart in 1995 [30]. This algorithm is inspired by the simulation of social psychological expression of birds and fishes. PSO includes two terms ?? ???????? and ?? ???????? . Position and velocity are updated over the course of iteration from these mathematical equations: 1 1 1 2 2 ( ) ( ) t t t t t t ij ij v wv c R Pbest X c R Gbest X + = + ? + ? (1) 1 1 t t t X X v + + = + ( ) i 1, 2...NP = And ( ) j 1, 2...NG =(2) Where max min max ( )* max w w iteration w w iteration ? = ? ,(3) w max =0.4 and w min =0.9. t ij v , 1 t ij v + is the velocity of aj th member of ani th particle at iteration number (t) and (t+1). (Usually C 1 =C 2 =2), r 1 and r 2 Random number (0, 1). # b) Moth-Flame Optimizer A novel nature-inspired Moth-Flame optimization algorithm [31] based on the transverse orientation of Moths in space. Transverse orientation for navigation uses a constant angle by Moths with respect to Moon to fly in straight direction in night. In MFO algorithm that Moths fly around flames in a Logarithmic spiral way and finally converges towards the flame. Spiral way expresses the exploration area and it guarantees to exploit the optimum solution [31]: Moth-Flame optimizer is first introduced by Seyedali Mirjalili in 2015 [31]. MFO is a populationbased algorithm; we represent the set of moths in a matrix: 1,1 1,2 1, 2,1, 2,2 2, ,1 ,2 , , , , , , , , d dn n n d m m m m m m M m m m ? ? ? ? ? ? ? ? ? ? = ? ? ? ? ? ? ? ? ? ? ? ? ?(4) Where n represents a number of moths and d represents a number of variables (dimension). For all the moths, we also assume that there is an array for storing the corresponding fitness values as follows: 1 2 . . n OM OM OM OM ? ? ? ? ? ? ? ? = ? ? ? ? ? ? ? ? ? ? (5) Where n is the number of moths. Note that the fitness value is the return value of the fitness (objective) function for each moth. The position vector (first row in the matrix M for instance) of each moth is passed to the fitness function and the output of the fitness function is assigned to the corresponding moth as its fitness function (OM 1 in the matrix OM for instance). Other key components in the proposed algorithm are flames. We consider a matrix similar to the moth matrix [31]: 1, 1 1,2 1, 2, 1 2, 2 2, , 1 , 2 , . . . . . . . . . .. . . . . . d d n n n d FL FL FL FL FL FL F FL FL FL ? ? ? ? ? ? ? ? = ? ? ? ? ? ? ? ?(6) Where n shows a number of moths and d represents a number of variables (dimension). We know that the dimension of M and F arrays are equal. For the flames, we also assume that there is an array for storing the corresponding fitness values [31]: 1 2 . . n OFL OFL OF OFL ? ? ? ? ? ? ? ? = ? ? ? ? ? ? ? ? (7) Where n is the number of moths. Here, it must be noted that moths and flames both are solutions. The variance among them is the manner we treat and update them, in the iteration. The moths are genuine search agents that move all over the search space while flames are the finest location of moths that achieves so far. Therefore, every moth searches around a flame and updates it in the case of discovering an enhanced solution. With this mechanism, a moth never loses its best solution. The MFO algorithm is three rows that approximate the global solution of the problems defined like as follows [31]: ( ) MFO I, P, T = (8) I is the function that yields an uncertain population of moths and corresponding fitness values. The methodical model of this function is as follows: { } : , I M OM ? ? (9) The P function, which is the main function, expresses the moths all over the search space. This function receives the matrix of M and takes back its updated one at every time with each iteration. : P M M ? (10) The T returns true and false according to the termination Criterion satisfaction: { } : , T M true false ? (11) In order to mathematical model this behavior, we change the location of each Moth regarding a flame with the following equation: ( , ) i i j M S M F = (12) Where indicate the moth, indicates the flame and S is the spiral function. In this equation flame FL n,d (search agent * dimension) of equation ( 6) modifies the moth matrix of equation (12). Considering these points, we define a log (logarithmic scale) spiral for the MFO algorithm as follows [31]: ( ) ( ) , * cos 2 bt i j i j S M F D e t F ? = +(13) Where: i D expresses the distance of the moth for thej th flame, b is a constant for expressing the shape of the log (logarithmic) spiral, and t is a random value in [-1, 1]. i j i D F M = ?(14) Where: i M represent the i th moth, F j represents the j th flame, and where i D expresses the path length of the i th moth for the j th flame. The no. of flames are adaptively reduced with the iterations. We use the following formulation: 1 * N flame no round N l T ? ? ? = ? ? ? ? ? (15) Where l is the present number of iteration, N is the maximum number of flames, and Tshows the maximum number of iterations. O MFO O t O Quick sort O position update = + 2 2 ( ) ( ( * )) ( ) O MFO O t n n d O tn tnd = + = +( ) ( ) ( ) ( ) ( ) Where n shows a number of moths, t represents maximum no. of iterations, and d represents no. of variables. # c) The Hybrid PSO-MFO Algorithm The drawback of PSO is the limitation to cover small search space while solving higher order or complex design problem due to constant inertia weight. This problem can be tackled with Hybrid PSO-MFO as it extracts the quality characteristics of both PSO and MFO. Moth-Flame Optimizer is used for exploration phase as it uses logarithmic spiral function so it covers a broader area in uncertain search space. Because both of the algorithms are randomization techniques so we use term uncertain search space during the computation over the course of iteration from starting to maximum iteration limit. Exploration phase means the capability of an algorithm to try out a large number of possible solutions. The position of particle that is responsible for finding the optimum solution to the complex non-linear problem is replaced with the position of Moths that is equivalent to the position of the particle but highly efficient to move solution towards optimal one. MFO directs the particles faster towards optimal value, reduces computational time. As we know that that PSO is a well-known algorithm that exploits the best possible solution for its unknown search space. So the combination of best characteristic (exploration with MFO and exploitation with PSO) guarantees to obtain the best possible optimal solution of the problem that also avoids local stagnation or local optima of the problem. A set of Hybrid PSO-MFO is the combination of separate PSO and MFO. Hybrid PSO-MFO merges the best strength of both PSO in exploitation and MFO in exploration phase towards the targeted optimum solution. 1 1 1 2 2 ( _ ) ( ) t t t t t t ij ij v wv c R Moth Pos X c R Gbest X + = + ? + ?(17) III. # Optimal Power Flow Problem Formulation As specified before, OPF is the optimized problem of power flow that provides the optimum values of independent variables by optimizing a predefined objective function with respect to the operating bounds of the system [1]. The OPF problem can be mathematically expressed as a non-linear constrained optimization problem as follows [1]: Minimize f(a,b) (18) Subject to s(a,b)=0(19) # i. Control variables The control variables should be adjusted to fulfill the power flow equations. For the OPF problem, the set for control variables can be formulated as [1], [4]: 2 1 1 1 [ ] , , , NTr NGen NGen NCom T G G G G C C P P V V Q b Q T T = ? ? ? ?(21) Where, G P = Real power output at the PV(Generator) buses excluding at the slack (Reference) bus. G V = Magnitude of Voltage at PV (Generator) buses. C Q = shunt VAR compensation. T = tap settings of the transformer. NGen, NTr, NCom= No. of generator units, No. of tap changing transformers and No. of shunt VAR compensation devices, respectively. The control variables are the decision variables of the power system which could be adjusted as per the requirement. # ii. State variables There is a need of variables for all OPF formulations for the characterization of the Electrical Power Engineering state of the system. So, the state variables can be formulated as [1], [4]: 1 1 1 1 ] [ , , , NLB NGen T l l L L G G G Nline P V V Q Q S S a = ? ? ?(22) Where, # b) Constraints There are two OPF constraints named inequality and equality constraints. These constraints are explained in the sections given below. # i. Equality constraints The physical condition of the power system is described by the equality constraints of the system. These equality constraints are basically the power flow equations which can be explained as follows [1], [4]. # a. Real power constraints The real power constraints can be formulated as follows: [ ( ) ( )] 0 NB i j ij ij ij Di Gi J i ij P P V V G Cos B Sin ? ? = ? ? + = ? (23) b. Reactive power constraints [ ( ) ( )] 0 NB i j ij ij ij ij Di Gi J i Q Q V V G Cos B Sin ? ? = ? ? + = ?(24) Where, j ij i ? ? ? = ? is the phase angle of voltage between buses i and j.NB= total No. of buses, G P = real power output, G Q = reactive power output, ij ij ij Y G jB = + shows the susceptance and conductance between bus i and j, respectively, ij Y is the mutual admittance between buses I and j. # ii. Inequality constraints The boundaries of power system devices together with the bounds created to surety system security are given by inequality constraints of the OPF [4], [5]. # a. Generator constraints For all generating units including the reference bus: voltage magnitude, real power and reactive power outputs should be constrained within its minimum and maximum bounds as given below [27]: , i i i upper l er G o G G w V V V ? ? i=1,?, NGen (25) i i i upper lower G G G P P P ? ? , i=1,?, NGen(26)i i i upper lower G G G Q Q Q ? ? , i=1,?, NGen (27) b. Transformer constraints Tap settings of transformer should be constrained inside their stated minimum and maximum bounds as follows [27]: i i i upper lower G G G T T T ? ? , i=1,?,NGen(28) c. Shunt VAR compensator constraints Shunt VAR compensation devices need to be constrained within its minimum and maximum bounds as given below [27]: i i i upper lower C GC C Q Q Q ? ? , i=1,?,NGen(29) d. Security constraints These comprise the limits of a magnitude of the voltage at PQ buses and loadings on the transmission line. Voltage for every PQ bus should be limited by their minimum and maximum operational bounds. Line flow over each line should not exceed its maximum loading limit. So, these limitations can be mathematically expressed as follows [27]: i i i lower upper L L L V V V ? ? , i=1,?,NGen (30) i i upper l l S S ? , i=1,?,Nline(31) The control variables are self-constraint. The inequality constrained of state variables comprises the magnitude of PQ bus voltage, active power production at reference bus, reactive power production and loadings on line may be encompassed into an objective function in terms of quadratic penalty terms. In which, the penalty factor is multiplied by the square of the indifference value of state variables and is included in the objective function and any impractical result achieved is declined [27]. Penalty function may be mathematically formulated as follows: ( ) 1 1 2 2 2 1 1 0 ( ) ( ) i i i i NLB NGen Nline aug P V L L G G Q S l l i i i lim lim max J J P P V V S S = = = = + ? ? + ? ? + ? + ? ? ? ? ?(32) Where, , , , The reactive power constraints can be formulated as follows: minimum limit lim U takings the value of that limit. This can be shown as follows [27]: P V Q S = ? ? ? ? ; ; upper upper lim lower lower U U U U U U U = ? > ? < ? (33) IV. # Application and Results The PSO-MFO technique has been implemented for the OPF solution for standard IEEE 30bus test system and for a number of cases with dissimilar objective functions. The used software program is written in MATLAB R2014b computing surroundings and used on a 2.60 GHz i5 PC with 4 GB RAM. In this work the HPSO-MFO population size is selected to be 40. # a) IEEE 30-bus test system With the purpose of elucidating the strength of the suggested HPSO-MFO technique, it has been verified on the standard IEEE 30-bus test system as displays in fig. 2. The standard IEEE 30-bus test system selected in this work has the following features [6] In addition, generator cost coefficient data, the line data, bus data, and the upper and lower bounds for the control variables are specified in [33]. In given test system, five diverse cases have been considered for various purposes and all the acquired outcomes are given in Tables 3, 5, 7, 9, 11. The very first column of this tables denotes the optimal values of control variables found where: -P G1 through P G6 and V G1 through V G6 signifies the power and voltages of generator 1 to generator 6. -T 4-12 , T 6-9 , T 6-10 and T 28-27 are the transformer tap settings comprised between buses 4-12, 6-9, 6-10 and 28-27. -Q C10 , Q C12 , Q C15 , Q C17 , Q C20 , Q C21 , Q C23 , Q C24 and Q C29 denote the shunt VAR compensators coupled at buses 10, 12, 15, 17, 20, 21, 23, 24 and 29. Further, fuel cost ($/hr), real power losses (MW), reactive power losses (MVAR), voltage deviation and Lmax represent the total generation fuel cost of the system, the total real power losses, the total reactive power losses, the load voltages deviation from 1 and the stability index, respectively. Other particulars for these outcomes will be specified in the next sections. The control parameters for HPSO-MFO, MFO, PSO used in this problem are given in table 1. In table 1, no. of variables (dim) shows the six no. of generators used in the 30 bus system. It gives the optimization values for different cases as they depends on the decision variables. In all 5 cases, results are the average value obtained after 10 number of runs. The very common OPF objective that is generation fuel cost reduction is considered in the case 1. Therefore, the objective function Y indicates the complete fuel cost of total generating units and it is calculated by following equation [1]: 1 ($ / ) NGen i i Y f hr = = ? (34) Where, i f is the total fuel cost of th i generator. i f , may be formulated as follow: 2 which displays that the results obtained by PSO-MFO are better than the other methods. The optimal values of control variables obtained by different algorithms for case 1 are specified in Table 3. By means of the same settings i.e. control variables boundaries, initial conditions and system data, the results achieved in case 1 with the PSO-MFO technique are compared to some other methods and it display that the total fuel cost is greatly reduced compared to the initial case [6]. Quantitatively, it is reduced from 901.951$/hr to 799.056$/hr. Case 2: Voltage profile improvement Bus voltage is considered as most essential and important security and service excellence indices [6]. Here the goal is to reduce the fuel cost and increase voltage profile simultaneously by reducing the voltage deviation of PQ (load) buses from the unity 1.0 p.u. # Fig. 3: Fuel cost variations with different algorithms Hence, the objective function may be formulated by following equation [4]: cost voltage deviation Y Y wY ? = + (36) Where, w is an appropriate weighting factor, to be chosen by the user to offer a weight or importance to each one of the two terms of the objective function. cos 1 NGen t i i Y f = = ? (37) _ 1 | 1.0 | NGen voltage deviation i i Y V = = ? ? (38) The variation of voltage deviation with different algorithms over iterations is sketched in fig. 3. It demonstrates that the suggested method has good convergence characteristics. The statistical values of voltage deviation obtained with different methods are shown in table 4 which display that the results obtained by PSO-MFO are better than the other methods excluding GSA method. The optimal values of control variables obtained by different algorithms for case 2 are specified in Table 5. By means of the same settings the results achieved in case 2 with the PSO-MFO technique are compared to some other methods and it display that the voltage deviation is greatly reduced compared to the initial case [6]. It has been made known that the voltage deviation is reduced from 1.1496 p.u. to 0.1056p.u. using PSO-MFO technique.GSA [2] gives better result than the HPSO-MFO method only in case of voltage deviation among five cases. Due to No Free Lunch (NFL) theorem proves that no one can propose an algorithm for solving all optimization problems. This means that the success of an algorithm in solving a specific set of problems does not guarantee solving all optimization problems with different type and nature. NFL makes this field of study highly active which results in enhancing current approaches and proposing new meta-heuristics every year. This also motivates our attempts to develop a new Hybrid meta-heuristic for solving OPF Problem. Case 3: Voltage stability enhancement Presently, the transmission systems are enforced to work nearby their safety bounds, because of cost-effective and environmental causes. One of the significant characteristics of the system is its capability to retain continuously tolerable bus voltages to each node beneath standard operational environments, next to the rise in load, as soon as the system is being affected by disturbance. The unoptimized control variables may cause increasing and unmanageable voltage drop causing a tremendous voltage collapse [6]. Hence, voltage stability is inviting ever more attention. By using various techniques to evaluate the margin of voltage stability, Glavitch and Kessel have introduced a voltage stability index called L-index depends on the viability of load flow equations for every node [34]. The L-index of a bus shows the probability of voltage collapse circumstance for that particular bus. It differs between 0 and 1 equivalent to zero load and voltage collapse, respectively. For the given system with NB, N Gen and NLB buses signifying the total no. of buses, the total no. of generator buses and the total no. of load buses, respectively. The buses can be distinct as PV (generator) buses at the head and PQ (load) buses at the tail as follows [4]: [ ] L L LL LG L bus G G GL GG G I V Y Y V Y I V Y Y V = = ? ? ? ? ? ?? ? ? ? ? ? ? ?? ? ? ? ? ? ? ?? ?(39) Where, LL Y , LG Y , GL Y and GG Y are co-matrix of bus Y . The subsequent hybrid system of equations can be expressed as: [ ] L L LL LG L G G GL GG G H I V V V I H H I H H = = ? ? ? ? ? ?? ? ? ? ? ? ? ?? ? ? ? ? ? ? ?? ?(40) Where matrix H is produced by the partially inverting of bus The matrix H is given by: [ ] 1 LL LL LG LL LL GL LL GG GL LL LG Z Z Y H Z Y Y Z Y Y Z Y ? ? ? ? = = ? ? ? ? ?(41) Hence, the L-index denoted by j L of bus j is represented as follows: 1 1 i j LG i j ji NGen v L H v = = ? ? j=1,2?,NL(42) Hence, the stability of the whole system is described by a global indicator max L which is given by [6], max max( ) j L L = j=1,2?,NL(43) The system is more stable as the value of max L is lower. The voltage stability can be enhanced by reducing the value of voltage stability indicator L-index at every bus of the system. [6]. Thus, the objective function may be given as follows: cos _ _ t voltage Stability Enhancement Y Y wY = +(44) Where, cos 1 NGen t i i Y f = = ? (45) max _ _ voltage stability enhancement Y L = (46) The variation of the Lmax index with different algorithms over iterations is presented in fig. 4. The statistical results obtained with different methods are shown in table 6 which display that PSO-MFO method gives better results than the other methods. The optimal values of control variables obtained by different algorithms for case 3 are given in Table 7. After applying the PSO-MFO technique, it appears from Table 7 that the value of Lmax is considerably decreased in this case compared to initial [6] from 0.1723 to 0.1126. Thus, the distance from breakdown point is improved. # Case 4: Minimization of active power transmission losses In the case 4 the Optimal Power Flow objective is to reduce the active power transmission losses, which can be represented by power balance equation as follows [6]: 9. By means of the same settings the results achieved in case 4 with the PSO-MFO technique are compared to some other methods and it display that the real power transmission losses are greatly reduced compared to the initial case [6] # Case 5: Minimization of reactive power transmission losses The accessibility of reactive power is the main point for static system voltage stability margin to support the transmission of active power from the source to sinks [6]. Thus, the minimization of VAR losses are given by the following expression: 1 1 1 i Gi Di NGen NGen NGen i i i J Q Q Q = = = = = ? ? ? ?(48) It is notable that the reactive power losses are not essentially positive. The variation of reactive power losses with different methods shown in fig. 6. It demonstrates that the suggested method has good convergence characteristics. The statistical values of reactive power losses obtained with different methods are shown in table 10 which display that the results obtained by hybrid PSO-MFO method are better than the other methods. The optimal values of control variables obtained by different algorithms for case 5 are given in Table 11. It is shown that the reactive power losses are greatly reduced compared to the initial case [6] from -4.6066 MVAR to -25.335MVAR using hybrid PSO-MFO method. Table 12 show the comparison of elapsed time taken by the different methods to optimize the different objective cases. The comparison shows that the time taken by all three algorithms is not same which indicates the different evaluation strategy of different methods. V. # Robustness Test In order to check the robustness of the HPSO-MFO for solving continues Optimal Power Flow problems, 10 times trials with various search agents for cases Case 1, Case 2, Case 3, Case 4 and Case 5. Table 2, Table 3, Table 4, Table 5, Table 6, Table 7, Table 8, Table 9, Table 10 and Table 11 presents the statistical results achieved by the HPSO-MFO, MFO and PSO algorithms for OPF problems for various cases. From these tables, it is clear that the optimum objective function values obtained by HPSO-MFO are near to every trial and minimum compare to MFO and PSO algorithms. It proves the robustness of hybrid PSO-MFO algorithm (HPSO-MFO) to solve OPF problem. # VI. # Conclusion Particle Swarm Optimization-Moth Flame Optimizer (PSO-MFO), Moth Flame Optimizer and Particle Swarm Optimization Algorithm are successfully applied to standard IEEE 30-bus test systems to solve the optimal power flow problem for the various types of cases. The results give the optimal settings of control variables with different methods which demonstrate the effectiveness of the different techniques. The solutions obtained from the hybrid PSO-MFO method approach has good convergence characteristics and gives the better results compared to MFO and PSO methods which confirm the effectiveness of proposed algorithm. VII. 1![Fig. 1: A conceptual model of position updating of a moth around a flame We utilize Quick sort algorithm, the sorting is in the ( log ) o n n best and](image-2.png "Fig. 1 :") ![a=vector of state variables, b=vector of control variables, f(a,b)=objective function, s(a,b)=different equality constraints set, h(a,b)=different inequality constraints set. The evaluation function for the OPF problem is given as follows: Evaluation Function = Search Agents * Maximum Iterations = 40*500= 20000 a) Variables](image-3.png "") ![Magnitude of Voltage at Loadbuses.G Q = Reactive power generation of all generators. l S = Transmission line loading. NLB, Nline= No. of PQ buses and the No. of transmission lines, respectively.](image-4.png "LV=") ![= active power load demand, D Q = reactive power load demand, ij B and ij G = elements of the admittance matrix ( )](image-5.png "DP") ![, [33]: NGen = No. of generators = 6 at buses 1,2,5,8,11 and 13, NTr = No. of regulating transformers having offnominal tap ratio = 4 between buses 4-12, 6-9, 6-10 and 28-27, NCom = No. of shunt VAR Compensators = 9 at buses 10,12,15,17,20,21,23,24 and 29 and NLB = No. of load buses = 24.](image-6.png "") 2![Fig. 2: Single line diagram of IEEE 30-bus test system Case 1: Minimization of generation fuel cost.The very common OPF objective that is generation fuel cost reduction is considered in the case 1. Therefore, the objective function Y indicates the complete fuel cost of total generating units and it is calculated by following equation[1]:](image-7.png "Fig. 2 :") ![and voltage deviation Y ? are specified as follows[4]:](image-8.png "costY") Year 2017V Version Iof Researches in Engineering ( ) Volume XVII Issue FGlobal Journalpenalty factors lim U = Boundary value of the state variable U.If U is greater than the maximum limit, lim U takings the value of this one, if U is lesser than the 1Sr. No.ParametersValue1Population (No. of Search agents) (N)402Maximum iterations count (t)5003No. of Variables (dim)254Random Number[0,1]5source acceleration coefficient (??_1, ??_2)2 2MethodFuel Cost ($/hr)Method DescriptionHPSO-MFO799.056Hybrid Particle Swarm Optimization-Moth Flame OptimizerMFO799.072Moth Flame OptimizerPSO799.704Particle Swarm OptimizationDE799.289Differential Evolution [15]BHBO799.921Black Hole-Based Optimization [6] 3Control VariableMinMaxInitialHPSO-MFOMFOPSOP G15020099.2230178.133177.055177.105P G220808048.95648.69848.748P G515505021.38521.30421.318P G810352021.70621.08420.986P G1110302010.00011.88312.049P G1312402012.00012.00012.000V G10.951.11.051.1001.1001.100V G20.951.11.041.0881.0881.088V G50.951.11.011.0621.0621.061V G80.951.11.011.0701.0691.070V G110.951.11.051.1001.1001.100V G130.951.11.051.1001.1001.100T 4-1201.11.0780.9391.0440.976T 6-901.11.0691.1000.9000.975 4MethodVoltage Deviation (p.u)Method DescriptionHPSO-MFO0.1056Hybrid Particle Swarm Optimization-Moth Flame OptimizerMFO0.1065Moth Flame OptimizerPSO0.1506Particle Swarm OptimizationGSA0.0932Gravitational Search Algorithm [2]DE0.1357Differential Evolution [15]BHBO0.1262Black Hole-Based Optimization [6] 5Control VariableMinMaxInitialHPSO-MFOMFOPSOP G15020099.2230177.650180.212175.922P G220808049.09249.58446.389P G515505015.00015.00021.597P G810352010.00024.34919.396P G1110302030.00012.65717.656P G1312402012.00012.00012.000V G10.951.11.051.0331.0331.047V G20.951.11.041.0171.0171.034V G50.951.11.011.0151.0050.999V G80.951.11.010.9970.9991.005V G110.951.11.051.0471.0710.999V G130.951.11.051.0161.0521.018T 4-1201.11.0781.0651.1000.954T 6-901.11.0690.9140.9000.969T 6-1001.11.0320.9731.0720.989 6MethodL maxMethod DescriptionHPSO-MFO0.1126Hybrid Particle Swarm Optimization-Moth Flame OptimizerMFO0.1138Moth Flame OptimizerPSO0.1180Particle Swarm OptimizationGSA0.1162Gravitational Search Algorithm [2]DE0.1219Differential Evolution [15]BHBO0.1167Black Hole-Based Optimization [6] 7Control VariableMinMaxInitialHPSO-MFOMFOPSOP G15020099.2230182.308177.299158.331P G220808045.36048.79249.050P G515505021.10921.31618.956P G810352021.55720.35131.224P G1110302010.00012.37015.906P G1312402012.00012.01217.801V G10.951.11.051.1001.1001.098 8MethodActive Power Loss (MW)Method DescriptionHPSO-MFO2.831Hybrid Particle Swarm Optimization-Moth Flame OptimizerMFO2.853Moth Flame OptimizerPSO3.026Particle Swarm OptimizationBHBO3.503Black Hole-Based Optimization [6] 9Control VariableMinMaxInitialHPSO-MFOMFOPSOP G15020099.223051.26951.25351.427P G220808080.00080.00080.000P G515505050.00050.00050.000P G810352035.00035.00035.000P G1110302030.00030.00030.000P G1312402040.00040.00040.000V G10.951.11.051.1001.1001.100V G20.951.11.041.1001.0981.100V G50.951.11.011.0821.0801.083V G80.951.11.011.0861.0871.090V G110.951.11.051.1001.1001.100V G130.951.11.051.1001.1001.100T 4-1201.11.0781.0441.0560.977T 6-901.11.0690.9010.9001.100T 6-1001.11.0320.9930.9821.100T 28-2701.11.0680.9870.9730.998QC 100505.0005.0004.065QC 120504.5705.0000.000QC 150504.9693.0705.000QC 170504.9425.0005.000QC 200504.3375.0000.000QC 210505.0005.0005.000QC 230505.0005.0005.000QC 240505.0005.0000.000QC 290502.4122.5080.000PLoss (MW)--5.82192.8312.8533.026 10MethodReactive Power Loss (MVAR)Method DescriptionHPSO-MFO-25.335Hybrid Particle Swarm Optimization-Moth Flame OptimizerMFO-25.204Moth Flame OptimizerPSO-23.407Particle Swarm OptimizationBHBO-20.152Black Hole-Based Optimization [6] 12Case No.Elapsed Time (Seconds)MFOPSOHPSO-MFO1166.2097250.2674211.79152191.8238266.5375229.68733196.6275270.3358243.29194161.6395248.8739259.97315173.5987253.3971209.4387 Year 2017 F Optimal Power Flow using a hybrid Particle Swarm Optimizer with Moth Flame Optimizer G P = Real power generation at reference bus.© 2017 Global Journals Inc. (US)Global Journal of Researches in Engineering © 2017 Global Journals Inc. (US) weighting function (w) 0.65© 2017 Global Journals Inc. (US)Global Journal of Researches in Engineering Year 2017 F Optimal Power Flow using a hybrid Particle Swarm Optimizer with Moth Flame Optimizer © 2017 Global Journals Inc. (US)Global Journal of Researches in Engineering Year 2017 F Optimal Power Flow using a hybrid Particle Swarm Optimizer with Moth Flame Optimizer Year 2017 F Optimal Power Flow using a hybrid Particle Swarm Optimizer with Moth Flame Optimizer Year 2017 F Optimal Power Flow using a hybrid Particle Swarm Optimizer with Moth Flame Optimizer Year 2017 F Optimal Power Flow using a hybrid Particle Swarm Optimizer with Moth Flame Optimizer ## Acknowledgment The authors would like to thank Professor Seyedali Mirjalili for his valuable comments and support. http://www.alimirjalili.com/MFO.html. * Optimal power flow using the league championship algorithm: a case study of the Algerian power system HR E HBouchekara MAAbido AEChaib RMehasni Energy Convers. Manag 87 2014 * Optimal power flow using gravitational search algorithm SDuman UGüvenç YSönmez NYörükeren Energy Convers 59 Manag2012 * A modified shuffle frogleaping algorithm for multiobjective optimal power flow TNiknam MNarimani MJabbari ARMalekpour Energy 36 11 2011 * Contribution à l'étude du Dispatching Economique JCarpentier Bull. Soc. Fran-caise Electriciens 1962 * Optimal power flow using Teaching-Learning-Based Optimization technique HR E HBouchekara MAAbido MBoucherma Electr. Power Syst. Res 114 2014 * Optimal power flow using black-hole-based optimization approach HR E HBouchekara Appl. Soft Comput 2013 under review * Optimal power flow: a bibliographicsurvey I, formulations and deterministic methods SFrank ISteponavice SRebennack Energy Syst 3 3 2012 * MRAlrashidi MEEl-Hawary Applications of computational intelligence techniques for solving the revived optimal power flow problem 2009 79 * Optimal power flow: a bibliographic survey II, nondeterministic and hybrid methods SFrank ISteponavice SRebennack Energy Syst 3 3 2012 * A comparative study of populationbased optimization algorithms for turning operations ARYildiz Inf. Sci. (Ny) 2012 * Optimal choice and allocation of FACTS devices in deregulated electricity market using genetic algorithms ILj Cai GErlich Stamtsis 2004 IEEE * A hybrid GA approach for OPF with consideration of FACTS devices TSChung Li IEEE Power Engineering Review 2001 * Enhanced Genetic Algorithm for Optimal Electric Power Flow using TCSC and TCPS KKalaiselvi KKumar Chandrasekar Proc. World 2010 II * Optimal power flow by enhanced genetic algorithm. Power Syst AGBakirtzis Biskas VZoumas Petridis IEEE Trans 17 2 2002 * Optimal power flow using differential evolution algorithm AAAbou El Ela MAAbido Electr. Power Syst. Res 80 7 2010 * Optimal Power Flow Solution of the Algerian Electrical Network using Differential Evolution Algorithm LSlimani TBouktir TELKOMNIKA Indonesian Journal of Electrical Engineering 10 2 2012 * Towards the concept of integrated security: optimal dispatch under static and dynamic security constraints. Electric Power Syst HBadrul Chowdhury Research 25 1992 * Application of particle swarm optimization technique for optimal location of FACTS devices considering cost of installation and system loadability MSaravanan PSlochanal Venkatesh Abraham Electric Power Syst 77 2007 Research * Optimal power flow using tabu search algorithm MAAbido Electric Power Compon Syst 30 2002 * Optimal power flow using gravitational search algorithm. Energy Convers US Duman YGüvenç NSönmez Yörükeren Manag 59 2012 * Application of biogeography-based optimisation to solve different optimal power flow problems ABhattacharya Pchattopadhyay IET Gener, Transm. Distrib 5 1 70 2011 * Application of Harmony Search to Optimal Power Flow Problems NSinsupan ULeeton TKulworawanichpong 2010 * Solution of optimal power flow using chaotic krill herd algorithm AMukherjee Mukherjee Chaos, Solitons & Fractals 78 2015 * Optimal Power Flow Using Cuckoo Optimization Algorithm. IJAREEIE MRama MohanaRao AVNareshbabu 2013 * Ant colony search algorithm for the optimal power flow problem JSoares Sousa Vale Morais Faria Gen, Meet 2011 IEEE Power Energy Soc * Voltage stability enhancement and voltage deviation minimization using BAT optimization algorithm NIndrajit MotilalTrivedi PradeepBhoye Jangir ASiddharth NarottamParmar ArvindJangir Kumar 10.1109/ICEES.2016.7510626 3rd International Conference on Electrical Energy Systems (ICEES) Chennai, India 2016. 2016 * Optimal power flow with enhancement of voltage stability and reduction of power loss using ant-lion optimizer Cogent Engineering Indrajit N. Trivedi, Pradeep Jangir & Siddharth A. Parmar 3 1208942 2016 * Voltage stability enhancement and voltage deviation minimization using Ant-Lion optimizer algorithm NIndrajit Trivedi ASiddharth RHParmar PradeepBhesdadiya Jangir 10.1109/AEEICB.2016.7538286 nd International Conference on Advances in Electrical, Electronics, Information, Communication and Bio-Informatics (AEEICB) Chennai, India 2016 2. 2016 * Voltage stability enhancement and voltage deviation minimization using Multi-Verse optimizer algorithm NIndrajit PradeepTrivedi NarottamJangir Jangir ASiddharth MotilalParmar ArvindBhoye Kumar 10.1109/ICCPCT.2016.7530136 2016 International Conference on Circuit, Power and Computing Technologies (ICCPCT) Nagercoil, India 2016 * Particle swarm optimization JKennedy REberhart Proceedings of the IEEE International Conference on Neural Networks the IEEE International Conference on Neural NetworksPerth, Australia 1995 * Moth-flame optimization algorithm: A novel nature-inspired heuristic paradigm SeyedaliMirjalili Knowledge-Based System 89 2015 * Moth-Flame optimization Algorithm for solving real challenging constrained engineering optimization problems NarottamJangir MaheshHPandya NIndrajit RHTrivedi PradeepBhesdadiya ArvindJangir Kumar 2016 IEEE Students' Conference on Electrical, Electronics and Computer