Solution for Optimal Power Flow Problem Using WDO Algorithm

Wind driven optimization (WDO) algorithm is a best optimization method based on atmospherically motion, global optimization nature inspired method. The method is based on population iterative analytical global optimization for multifaceted and multi prototype in the search domain for constraints to implement. In this paper, WDO algorithm is accustomed to find optimal power flow solution. To find the efficacy of the technique, it is applied to IEEE 30 bus systems to find fuel cost for generation of power as a main objective. Obtained results were compared with other techniques shows the better solution for optimal power flow problem.


Introduction
For power system control and functioning, optimal power flow solution is back bone (R.B., 1961; Bhasha, 2020). Optimizing and optimal operating in a satisfying operating constrains is main objective of optimal power flow problem. Over a half of s century optimal power flow problem studied. There were many deterministically optimization method has successfully implemented. Some are gradient based technique sequential quadratic type algorithm; Newton based methods and interior point method. Optimality and grow with come theoretical suppositions like convexity, complexity and other are not suitable for optimal power flow condition. A survey for commonly used methods is given in (Huneault , 1991;ChinnamahammadBhasha, 2020). Many of the optimization methods are mostly used in industrial applications and have some short falls (Balamurugan, 2020). There are so many drawbacks with these methods. Some of them are; they are not able to handle integer variables or binary they cannot global (ChinnamahammadBhasha, 2020; Deepthi, 2019; Aroulanandam, 2020). Nowadays, rapid growth in computational tools has boosted investigation in area of nondeterministic that is 'heuristic' methods used to work out the optimization problems in last two decades are Biogeography based optimization (BBO),bacterial foraging algorithm (BFA), Artificial neutral network(ANN), Ant colony algorithm (ACA), differential evaluation (DE), particle swarm optimization (PSO),algorithm of chaos, evolutionary algorithm, evolutionary programming methods, evolution strategies (ES), gravitational search algorithm (GSA) (ChinnamahammadBhasha, 2020; Balamurugan, 2020), Tabu search (TS), for finding optimal solutions and avoided to be confine with local one these methods are known for capabilities, quick search of large solutions and their capacity to uncertainty in area of power systems (Balamurugan, 2020). A brief review on optimal power flow solution using optimization methods is given in (Frank, 2012; Karimulla, 2019; Garikapati, 2020; Latchoumi, 2020; Pavan, 2020).
WDO is the recently developed optimization method, and nature inspired optimization method. It has been proved that WDO is uncomplicated and successful in performing multifaceted optimization problems. The aim of this article is to find solution for optimal power flow problem using WDO method. Under the objective of generation fuel cost minimization, performance of WDO is tested and sought on IEEE-30 bus stranded test system.
In this paper first OPF is mathematically evaluated after that WDO is presented. To solve optimization issue the power operating constraints WDO method is used at last concluded with important points.

Optimal Power Flow
The non linear mathematical representation can be given the form of formulae as Minimize f(x,u) (1) Subjected to g(x,u) =0 (2) h(x,u)≤ 0 (3)

Fuel cost objective
The simple quadratic cost equation for i th unit for active power is given as F(PGi) = AiPGi 2 + BiPGi + Ci (4) Where Ai ,Bi and Ci are generators coefficients for fuel cost . Fuel cost generation Ft of all generators (Ng) can calculate using the expression A1= min(Ft) = ∑ f i (P Gi ) N g =1 (5) It must be satisfy the inequality as well as equality constrains while minimizing fuel cost and constraints can be given in equation form as:

Equality constraints are:
The constraints are load flow expressions which are solved in conventional load flow procedure. Balance expressions for active, reactive powers in load flow are: Where P Gi , PDi QG and QDi i are active and reactive power generation and demand at i th bus .
Nbus is total number of buses , θij is phase angle difference, Yij is the admittance magnitude connected between bus i th and j th

Generator constraints:
The other constraints are, tap setting transformer tap setting ≤ ≤ ∀ i NG Compensation limit of reactive power by shunt compensator:

Security constraints are:
Voltage value constraint at load buses loadings at transmission lines are comes under security constrains. Every load bus voltage can be constricted within its upper and lower operating limits. The constraints can be given in equation form as: Here Ti, VG, PG and Qshu are self restricted inequality constraints and can be satisfied with in load flow problem, remaining three constrains along with active power are the non-self-restricted inequality constraints can operated by penalty approach (Ahrens, 2007). By this generalized load flow problem can be expressed as x is value of PG, QGm, and Vm

Wind Driven Optimization (Wdo) Algorithm:
WDO algorithm is best optimization method based on atmospherically motion, and is global optimization nature inspired method. It is observed that wind driven optimization method is uncomplicated and successful in solving multifaceted optimization problems. The method is best suitable for population build heuristic iterative global optimization technique for multimodel and multifaceted search domain to implement different constraints. When compared with other iterative methods in WDO method, updated velocity equations, carioles and gravitational forces used to provide degree of freedom and robustness for fine tune. With the help of algorithm we can easily explain the WDO as: Step 1: initialize max no. of iterations and population size.
Step 2: Assign the boundaries, pressure function and coefficients Step 3: Assign velocity and random position Step 4: Calculate pressure for all air parcels Step 5: update velocity Step 6: check velocity limits.
Step 7: update position and check boundaries Step 8: repeat the steps i.e.4,5,6 and 7 until number iterations reached.
Step 8: display the result.

Mathematical approach for WDO:
Wind blows in atmosphere to make air pressure uniform, we know that air always moves to low pressure from a high pressure at some velocity, and proportional to pressure gradient (Bayraktar, 2010). Further in WDO some assumptions and simplified formulae which are derived from basic fundamentals are used. The starting step of WDO is based on Newton's 2 nd law, issued to find accurate result for examination of atmospheric motion (Robinson, 2004;.
Here Fl is density of air, ⃗ acceleration vector Fi is the forces acting on mass.
Ideal gas law gives the relation between air density, acceleration vector and temperature is given as below: P= Rt Here P, R and t are pressure, universal gas constant and temperature.
In equation 9 there are fore forces, which cause te wind to move in particular direction and deflect its path. Noticeable force give rise to air in pressure gradient force (Fpg), frictional force (Ff) acts to opposite motion as described in equation 14. Gravitational force (Fg) is acts in a vertical direction in three dimensional atmosphere when a N dimensional space mapped with it. Coriolis force (Fc) caused by movement of earth, and averts path of wind. In this algorithm a motion in one direction that effect in another is implemented (Baraktar et al. 2010).
The equations for the above mentioned forces is given by Here Δ⍴ : Represents pressure gradient, δV: is a inappreciable air volume, Ω : is movement of earth, G:

Represents gravitational acceleration and u is wind velocity vector. Addition of all forces is given in equation 15
⍴ ⃗⃗⃗⃗ Δ = ⍴ ⃗ + (−2 ⃗⃗⃗⃗ * Ω) + ( − ⃗⃗⃗⃗⍴ ) + (−Δ⍴δV) (15) Assume an inappreciable air parcel moving with the velocity update expression which is given by equation number 16. From ideal gas equation from (10), ⍴ in pressure, and a step unit time (Δt = 01) assumed for simplicity. Updated equation for velocity id given as: From equation 16 we can observe that updated velocity value unew depends on the previous iteration value uold.
The fitness value in genetic algorithm (GA) is analogous to pressure term in wind driven optimization technique. If it is compared with particle swam optimization (PSO) same updated velocity equation can be perceive. Once velocity is updated site of air parcel also should update as: ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ = ⃗⃗⃗⃗⃗⃗⃗⃗ + ⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ * Δt (17) In equation 17, we know that air parcel continue to move from one position to another creates friction. Δt Is force of maximal pressure that is the global best site for optimization. Like this WDO is simple and successful to solve the multifaceted complex problems. WDO allows air parcel to move within the boundary of [-01, 01]. There are various boundary conditions presented in particle based optimization (Abido, 2002). Every air parcel influenced by gravitational force, are brought to boundaries back into search space. It is also being noted that per iteration for air parcels updated velocities are restricted. It avoids the air parcels to take number of steps and over looking in search space. Velocity dimension is restricted by its limits and is given by the equation 18 umax is adjusted velocity after , restricted to high speed.

Results
In order to check the WDO method, it has applied on IEEE 30 bus stranded test system. Program is developed in MATLAB software environment. IEEE 30 bus test system includes of six generator buses located at buses 1,11,2,8,5,and 13, four off nominal tap changing transformers with tap ratios at lines28-27,6-10, 6-9,4-12,and and it also consists of nine shunt VAR compensators located at bus numbers 10,12,15,17,20,21,23,24,and 29 (Lee , 1985;Abou, 2010) . Minimum and maximum boundaries, generator data, bus data and line data are available in (Bouchekara, 2013). By taking fuel cost for generator as an objective function proposed WDO algorithm is tested for OPF problem.
The WDO method is run for OPF and optimal solution obtained is tabulated in table 1. It has been observed that total fuel cost for generation is minimized as compared with other method. It is observed that generator fuel cost is reduced from 901.9416 $/h to 799.069 $/h which is almost 11.39%.
Under same constraints i.e. system data control variable limits results obtained from WDO is compared with other methods as given in literature. Table 2 gives clear comparison of WDO with other methods. It has proved that proposed method is best, uncomplicated and successful method to solve multifaceted problems. Particle swarm optimization [17] 800.96 Genetic algorithm [18] 805.94

Conclusion
In this paper, novel method is proposed to solve OPF issue in power systems that are Wind driven optimization method. : WDO algorithm is best optimization method based on atmospherically motion, and it is global optimization nature inspired method. On a IEEE 30 bus system WDO has been effectively and successfully applied, by considering minimization generation fuel cost as objective function. This paper is highly significant among the others for the following reasons i) High ranking when compared with other techniques ii) Efficacy and simple method Efficiency of the WDO is proved by carried out analysis that has strength by converging optimal value is very close to it.