International Journal of Electrical and Computer Engineering (IJECE) Vol. 8, No. 2, April 2018, pp. 780~792 ISSN: 2088-8708, DOI: 10.11591/ijece.v8i2.pp780-792  780 Gravitational-Search Algorithm for Optimal Controllers Design of Doubly-fed Induction Generator Saher Albatran1, Muwaffaq I. Alomoush2, Ahmed M. Koran3 1 Department of Electrical Engineering, Faculty of Engineering, Jordan University of Science and Technology, Jordan 2,3 Department of Electrical Power Engineering, Hijjawi Faculty for Engineering Technology, Yarmouk University, Jordan Article Info ABSTRACT Article history: Recently, the Gravitational-Search Algorithm (GSA) has been presented as a promising physics-inspired stochastic global optimization technique. It takes its derivation and features from laws of gravitation. This paper applies the GSA to design optimal controllers of a nonlinear system consisting of a doubly-fed induction generator (DFIG) driven by a wind turbine. Both the active and the reactive power are controlled and processed through a back-toback converter. The active power control loop consists of two cascaded proportional integral (PI) controllers. Another PI controller is used to set the q-component of the rotor voltage by compensating the generated reactive power. The GSA is used to simultaneously tune the parameters of the three PI controllers. A time-weighted absolute error (ITAE) is used in the objective function to stabilize the system and increase its damping when subjected to different disturbances. Simulation results will demonstrate that the optimal GSA-based coordinated controllers can efficiently damp system oscillations under severe disturbances. Moreover, simulation results will show that the designed optimal controllers obtained using the GSA perform better than the optimal controllers obtained using two commonly used global optimization techniques, which are the Genetic Algorithm (GA) and Particle Swarm Optimization (PSO). Received Oct 22, 2017 Revised Dec 25, 2017 Accepted Jan 7, 2018 Keyword: DFIG Global optimization Gravitational search algorithm Optimal controllers Wind energy Copyright © 2018 Institute of Advanced Engineering and Science. All rights reserved. Corresponding Author: Saher Albatran, Department of Electrical Engineering, Faculty of Engineering, Jordan University of Science and Technology, Irbid 221-10, Jordan. Email: saalbatran@just.edu.jo 1. INTRODUCTION Generated electric power from wind is growing at an increasing annual rate worldwide in the last decade as its cost has decreased considerably. Widespread use of wind power can largely reduce production costs, reduce carbon emissions, contribute to resolving global warming calls, and reduce the dependency on fossil fuels [1]. With the recent developments in power electronics technology, DFIG driven by a wind turbine [2] are becoming very popular and their worldwide market share is increasing quickly. DFIG has been widely used in grid-connected variable-speed constant-frequency wind energy conversion systems (WECS) due to their advantageous features such as higher efficiency, lower cost, wide-range operations, capability of producing reactive power, better energy capture, smoother operation, and flexible control. Therefore, with the increasing utilization of wind power in the interconnected power network, improving DFIG transient performance is becoming more significant in power system operations. Tuning the parameters of a classical PI controller of a grid-connected DFIG-based WECS is complicated due to the high nonlinearity and complexity of the system. Optimization techniques are applied Journal homepage: http://iaescore.com/journals/index.php/IJECE Int J Elec & Comp Eng ISSN: 2088-8708  781 to determine the optimal parameters of a classical PI controller to enhance the performance of the WECS as its response change certainly with time. In [3] an adaptive fuzzy gain scheduling of PI controller is presented. The PI controller parameters are optimized online using an adaptive mechanism based on fuzzy logic. In [4] a non-dominated sorting GA alongside discrete Fourier transform is implemented to quantify the damping characteristics of the rotor-side current due to system disturbance. In recent years, different evolution algorithms [5], [6] including mean variance optimization (MVO) algorithm [7], GA [8], PSO [9-12], bacteria foraging optimization (PFO) algorithm [13], hybrid artificial neural network [14] and Gbest guided artificial bee colony algorithm [15], [16] are applied to optimize the parameters of the PI controller of the rotor side converter (RSC) and hence improve the damping of oscillatory modes in DFIG-based WECS. In [17] a PSO is presented to generate an on-off controller to ensure robustness and quality of the energy produced, wherein [18] a nonlinear rotor-side controller is designed based on H2 optimal control theory to demonstrate the synthesis of a maximum power point tracking (MPPT) algorithm. Authors in [19] presented an optimal tuning of the MPPT algorithm and the active disturbances rejection controller to reduce the impact of the change on the characteristics of DFIG-based WECS. In [20], authors presented a computational intelligence control method to optimize the performance of the back-to-back converter and to regulate the low voltage ride-through capability of DFIG-based WECS. In [21], the DFIG-based WECS small signal stability was analyzed and their impacts on the eigenvalues of a DFIG in a single machine infinite bus system was investigated. In [22], the states of DFIG are used to obtain the optimal feedback control with utilizing the rotor currents beside the rotor voltages. Fuzzy logic nonlinear controller and particle swarm optimization are used in [23] to handle the nonlinearities in DFIG and to optimize the controllers’ gains. In [24], the authors optimized the performance of the grid side inverter (GSI) and RSC of grid-connected DFIG to minimize the tracking error based on discrete-time PSO inverse optimal controller. Differential geometry theory is used in [25] to control linearized DFIG-wind turbine model based on linear quadratic regulator. To enhance the DFIG during and after power system transients by optimizing the PI-controllers of the RSC, PSO is used in [26] and GA in [27]. In [28], the transients of the DFIG are optimized by minimizing the sum squared error deviation of the dc link voltage and the generator speed deviation. GSA is one of the recent heuristic algorithms, which has been introduced as a physics-inspired method in 2009 [29], [30]. It takes derivation and the main characteristics of heuristics from “Newton’s law of gravitation”; “Every particle in the universe attracts every other particle with a force that is directly proportional to the product of their masses and inversely proportional to the square of the distance between them” [31]. The GSA updates gravitational and inertia masses with the help of heavy masses and finds the optimum [29-31]. In the last few years, GSA has been successfully applied to solve world engineering optimization problems and demonstrate high quality performance in solving these problems [32-35]. In this paper, GSA is used to design optimal controllers for improving transient response of a nonlinear DFIG-based WECS. The optimal controllers are tested on different common types of disturbances. The results are compared to the results obtained using the two well-recognized global optimization algorithms, which are GA and PSO, to show the effectiveness of using GSA to attain a global optimal solution of the design problem under study. 2. MODEL OF THE DFIG-BASED WECS Figure 1 depicts a grid-connected DFIG-based WECS. The system shown includes a wind turbine, an induction generator, a back-to-back converter (RSC and GSI) [36], and the utility-grid. The ac/dc and the dc/ac conversions are used to control the bidirectional power delivered from/to the rotor circuit and to/from the grid, respectively. The dynamic electrical and mechanical equations are summarized below. The induction generator considered in this paper is a wound rotor induction machine. If the magnetic branch is assumed linear, the dynamic modeling of DFIG in arbitrary reference frame can be written as seen from the stator side as follows: vqs  rsiqs  ds  pqs (1) vds  rsids  qs  pds (2) v'qr  r 'r i'qr (  r ) 'dr  p 'qr (3) (   r ) ' qr  p ' (4) v' dr  r ' r i' dr dr Gravitational-Search Algorithm for Optimal Controllers Design of Doubly-fed …. (Saher Albatran)  782 ISSN: 2088-8708 qs  Llsiqs  M (iqs  i'qr ) GEAR BOX (5) DFIG AC GRID ar br DC LINK vdc cr GENERATOR SIDE Zg as Zg bs Zg cs GRID SIDE Figure 1. General configuration of grid-connected DFIG-based WECS ds  Llsids  M (ids  i'dr ) (6)  'qr  L'lr i'qr  M (iqs  i'qr ) (7)  'dr  L'lr i'dr  M (ids  i'dr ) (8) are the stator voltages, currents, and flux linkages in qd-reference where v qs , v , i qs , i , qs and  ds ds ds ' ' , v ' , i ' , i ' , ' frame, respectively. v qr dr qr dr qr and dr are the rotor voltages, currents, and flux linkages in qd reference frame, respectively.  is the arbitrary angular speed,  r is the electrical angular speed of the rotor, p is the derivative operator, rs is the stator windings resistance, r' r is the rotor windings resistance, M  3 2 Lms , Lms is the stator magnetizing inductance, L'lr is the rotor leakage inductance, and Lls is the stator leakage inductance. If indicates number of poles, the electromechanical developed torque Td can be expressed as: Td  3P ( iqs  qsids ) 2 2 ds (9) The dynamics of the mechanical shaft can be expressed based on the torque applied on the shaft Tsh from both the generator end and the turbine end as follow: Tsh  T  J d m d  Tm  J m dt dt (10) where T is the applied mechanical torque at the turbine side, Tm is the electrical torque seen at the generator side, J  and J m are the inertia of the turbine and the generator, respectively, and  and m are the angular speed of the turbine and the generator, respectively. The existence of the speed difference on the shaft ends and the compliance coefficient of the shaft K sh can determine the shaft-torsion (  g   ) as in (11): T  g     m     dt  sh K sh Int J Elec & Comp Eng, Vol. 8, No. 2, April 2018 : 780 – 792 (11) Int J Elec & Comp Eng  ISSN: 2088-8708 783 As shown in Figure 2, both the active and the reactive power are controlled through the back-toback converter. The active power control loop, which sets the d-component of the rotor voltage v' , consists dr of two cascaded PI-controllers; the first one is used to regulate the rotor speed and to set the reference active power to the second PI-controller which is used to compensate the generated active power Ps . Another PI controller is used to set the q-component v'qr by compensating the generated reactive power Q s . The mathematical description of the PI controllers is given by: GENERATOR SIDE ar AC GRID as ROTOR SIDE STATOR SIDE bs br DFIG LPF is,abc θm LPF is,dq ir,abc ir,dq vs,dq Speed Measurement ωm Ps SPEED CONTROL ωm Qs,ref ωg(t) Ki(ω) Kp(Ps) ∑ ∫ωg(t)dt y(t) ∑ x(t) Ps ∑ Modulator Ki(Ps) Kp(Qs) w(t) v'dr vr,dq v'qr Power Calculation Qr Kp(ω) ∑ LPF vr,abc abc/dq LPF Power Calculation Qs LPF dθ/dt vs,abc abc/dq ωm,ref vdc cr cs Ki(Qs) Pr ACTIVE POWER CONTROL ∑ ∫x(t)dt REACTIVE POWER CONTROL ∑ ∫w(t)dt Qs Figure 2. Block diagram of the grid-connected DFIG-based WECS with its associated control system ' (t )  K vdr pPS x(t )  KiPS y (t )  K p g (t )  Ki t  x( )d 0 (12) t   g ( )d 0 ' (t )  K vqr pQS w(t )  KiQS (13) t  w( )d 0 (14) The constants K p , K i , K pP , K iP , K pQ , KiQ are the proportional gains and integral S S S S gains of the PI controller as shown in Figure 2. Each gain of them is responsible for doing a specific task made on the input deviation signal; proportional gains ( K p , K pP , K pQ ) are concerned with system S S stability; the integral gains ( K i , K iP , KiQ ) are used of the steady-state error. In order to restrict the S S level of output of the PI controlled speed fluctuation during transient conditions, limits are imposed as shown in Figure 2. The parameters of the control system to be optimally tuned are the six PI controller constants ( K p , K i , K pP , K pQ , KiQ ). In control system analysis and design, performance indices are S S S used to assess a system behavior. This behavior is judged to be optimal when the performance induces reach minimum values and all system constraints are satisfied. In this paper, the integral of the time-weighted absolute error (ITAE) measure is used, which is given as in (15). Gravitational-Search Algorithm for Optimal Controllers Design of Doubly-fed …. (Saher Albatran) 784  ISSN: 2088-8708 ts ITAE   t | e(t ) | dt 0 (15) where t s is the simulation time, and e(t ) is the error signal which represent the difference between a reference value and a measured variable. The error signal is selected as the speed error signal ( e(t )  m,ref  m ). The optimal values of controller gains are optimally tuned solving the following generalized optimization problem: Minimize ITAE ( K p , Ki , K pQ , KiQ , K pP , KiP ) S S S S Subject to (16) K min  K p  K max p p K min  K i  K max i i K min  K pQS  K max pQS pQS K min  KiQS  K max iQS iQ S K min  K pPS  K max pPS pPS K min  KiPS  K max iPS iPS 3. GRAVITATIONAL SEARCH ALGORITHM GSA is a physics-inspired optimization method as it is derived and takes derivation and the main features of heuristics from Newton's laws of gravity. It depends on the gravity and mass interactions [29-36]. Based on Newtonian theory, three masses are assigned to a body, which describe its quite different properties; inertial mass ( M I ), passive gravitational mass ( M P ), and active gravitational mass ( M A ) [29-31]. The mass M I is the quantity presented in Newton's second law, which is a measure of a body's inertia. It is a measure of body's resistance to change in motion when a force is applied on the body, i.e., the body opposes that change through its inertial mass. The larger M I , the harder it is to accelerate for a given strength of force. M P is a measure of body's reaction to a gravitational field, i.e., it measures a body's response to being placed in a gravitational field, which implies that if M P is placed at some point that has gravitational potential, then it will experience a force on it. M A is a measure of body's source strength for producing a gravitational field by itself [29-32]. Newton’s law of gravity indicates that the force between two particles is directly proportional to the product of their masses and inversely proportional to the square of the distance between them, which can be written as [29]: F  G M1M 2 / R 2 (17) where F is the magnitude of the gravitational force, G is the gravitational constant, M1 and M 2 are the masses of the first and second particles, respectively, and R is the distance between them. Newton’s second law indicates that a particle's acceleration ( a ) depends on the force applied to the particle and the inertial mass ( M I ) of the particle, which is expressed as: a  F / MI Int J Elec & Comp Eng, Vol. 8, No. 2, April 2018 : 780 – 792 (18) Int J Elec & Comp Eng  ISSN: 2088-8708 785 As in the case of physics where performance of objects (particles) is measured by their masses, the agents in this method are attracted to each other based on gravitational force, and their tendency is towards heavier masses, and the heavy masses move in a manner much slower than the lighter ones, and this represents the exploitation behavior. The GSA algorithm depends on the fact that every mass (agent) has four specifications: position, inertial mass, active gravitational mass, and passive gravitational mass. The algorithm updates gravitational and inertia masses with the help of heavy masses and finds the optimum. This artificial world of masses obeys Newton’s law of gravity and motion. To present the GSA algorithm, consider a system with N agents (masses). The position of the ith agent is defined by [29]: d n i  1, 2, ..., N X i  ( x1 i ... xi ...xi ) ; (19) where xid is the position of ith agent in the dth dimension. The force acting on mass specific time t , is defined as follows [29]: Fijd  G(t ) M Pi (t ) M Aj (t ) Rij (t )   agent j , G(t ) is the gravitational constant at time i from mass j at a ( x dj (t )  x dj (t )) (20) where M Pi (t ) is the passive gravitational mass of agent distance between two agents i t ,  i , M Aj (t ) is the active gravitational mass of is a small constant, and Rij (t ) is the Euclidian and j [29]: Rij (t )  X i (t ), X j (t ) (21) The total force that acts on agent i in a dimension d is a randomly weighted sum of the dth components of the forces exerted by other agents, which is expressed as [29]: Fid (t )  N  j 1, j i rand j Fijd (t ) (22) where rand j is a random number in the interval [0, 1] , which is introduced the equation to give a stochastic characteristic to the algorithm. If M Ii (t ) is the inertial mass of ith-agent, then by the law of motion, the acceleration in direction d of the agent i at time t , is given by [29]: aid (t )  Fid (t ) / M Ii (t ) (23) The next velocity of an agent is considered as a fraction of its current velocity added to its acceleration. Therefore, its position and its velocity could be calculated as follows [29]: vid (t  1)  ri  vid (t )  aid (t ) (24) xid (t  1)  xid (t )  vid (t  1) (25) where ri is a uniform random variable in the interval [0, 1] , which is employed in the equation to give a randomized characteristic to the search. The gravitational constant G is a function of the initial value ( Go ) and time ( t ), i.e., it is initialized at the start and will be decreased as time passes to control the search accuracy. This constant is expressed as [29]: Gravitational-Search Algorithm for Optimal Controllers Design of Doubly-fed …. (Saher Albatran) 786  ISSN: 2088-8708 G(t )  G(Go , t )  Go e ( / T )t (26) where  is a constant specified by a user, t is current iteration and T is the total number of iterations. Fitness evaluation is used in this algorithm to calculate gravitational and inertia masses. A more efficient agent has a heavier mass, which indicates that better agents have higher attractions with slower walks. The values of gravitational and inertia masses, which are assumed equal, are calculated using the map of fitness. The gravitational and inertial masses are updated by the following Equations [29]: N M i (t )  mi (t ) /  m j (t ) j 1 (27) mi (t )  ( Fi (t )  W (t )) /(B(t )  W (t )) (28) M Ai  M Pi  M Ii  M i ; i  1, 2, ..., N (29) where Fi (t ) represents the fitness value of the agent i at time B(t ) are defined as follows [29]: B(t )  W (t )  t . For a minimization problem, W (t ) and min F j (t ) j{1,...,N} max F j (t ) j{1,...,N} (30) (31) For a maximization problem, W (t ) and B(t ) are defined as follows [29]: B(t )  W (t )  max F j (t ) j{1,...,m} min F j (t ) j{1,...,m} (32) (33) Equation (25) can be modified by reducing the number of agents with lapse of time to perform a good compromise between exploration and exploitation, where only a set of agents with bigger mass apply their forces to the others. Let K B is the set of the first K agents with the best fitness values and biggest masses, (25) can be written as [29]: Fid (t )  N  jKb , j i rand j Fijd (t ) (34) The function K b is a function of time, which has an initial value of K o at the beginning and decreases as time passes, which implies that all agents apply the force at the beginning, and as time passes Kb is decreased linearly. Finally, just one agent with the best fitness values and biggest masses apply force to the others. 4. SIMULATION RESULTS This section presents the simulation results of the DFIG-based WECS using the GSA-based coordinated optimal tuning of PI controllers. Per-unit system parameters are summarized in Table 1. In order Int J Elec & Comp Eng, Vol. 8, No. 2, April 2018 : 780 – 792 Int J Elec & Comp Eng  ISSN: 2088-8708 787 to indicate the effectiveness of the GSA to tune the parameters of PI controllers, the results will be compared with the results obtained using GA and PSO as shown Table 2. The GA and PSO are selected in the paper for comparisons as they are becoming very common for solving the optimization problems in different disciplines, essentially in the design of the controllers because of being robust in finding an optimal solution and capable to attain a near-optimal solution close to a global minimum. Table 1. Parameters of the Simulated System rs 0.05 [p.u.] Stator windings resistance Rotor windings resistance r'r 0.02 [p.u.] Lms 5.00 [p.u.] Rotor leakage inductance L'lr 0.10 [p.u.] Stator leakage inductance Lls K sh 0.10 [p.u.] Inertia of the turbine H 2.50 [sec] Inertia of the generator Hg 0.50 [sec] Stator magnetizing inductance Compliance coefficient of the shaft 0.40 [p.u./rad] Two disturbances are discussed in this paper to show the advantages of using the GSA to tune the PI controller over the other algorithms (GA and PSO). Two performance indices are used to compare between the results. The first one is the objective function (ITAE) and the second one is the settling time for 5% criterion. The comparisons between the three algorithms are listed in Table 2 and Table 3. Table 2. Optimal Controller Parameters of GSA, GA, and PSO Fault Disturbance Reference Speed Disturbance GSA 0.10585 GA 0.28642 PSO 0.61399 GSA 1.3265 GA 1.38187 PSO 4.63480 20.00 00 16.881 35 19.146 40 20.000 0 18.954 94 17.629 59 K i 20.00 00 16.384 74 14.289 63 20.000 0 19.036 54 19.027 92 KpQ 0.659 10 0.5154 7 0.4850 1 10.131 7 0.0826 3 0.8635 0 KiQ 0.050 00 3.9751 7 17.598 22 20.000 0 3.7749 6 6.4840 3 KpPs 0.050 00 0.2964 6 0.0500 16.396 43 9.2915 5 14.761 53 15.920 4 14.143 75 18.090 68 ITAE Kp KiPs 0.641 30 0.1587 1 3.1956 2 Table 3. Calculated Settling time [in seconds] Fault Disturbance Reference Speed Disturbance GSA GA PSO GSA GA PSO 0.82 1.19 1.38 2.40 2.71 2.82 0.10 0.23 0.85 1.56 1.79 1.85 ωm 0.04 0.06 0.11 1.04 1.42 1.54 Ps 0.83 0.96 1.13 2.34 2.6 2.68 Qs 0.22 0.29 0.94 0.47 0.72 0.72 v’qr v’dr By solving the optimization problem using GSA, GA and PSO, the optimal controller parameters of the system under study are obtained and summarized in Table 2. As can be seen, the values of ITAE Gravitational-Search Algorithm for Optimal Controllers Design of Doubly-fed …. (Saher Albatran) 788  ISSN: 2088-8708 associated with the GSA-based optimal controllers are better than those obtained using the GA-based and PSO-based optimal controllers. This fact will be reflected in the following response subsections, Figures and Tables. Only five variables ( v'qr , v' , Ps , Qs and  m ) are presented in this paper while the response of dr the other variables can be estimated from the presented results. 4.1. Three-phase to Ground Fault Three-phase to ground fault was applied after the system reaches to the steady state (t=10 s) and cleared after 30 ms. From Table 2, 60% and 80% reduction of ITAE is achieved using GSA relative to GA and PSO, respectively. Figure 3 and Table 3 show that v' and v'qr can reach the new steady state much dr faster when GSA is adopted. Figure 4 and Table 3 manifest that using GSA helps the rotor speed to damp in 40 ms which is two thirds of the time required when GA is used and about one third of the settling time needed when PSO is used. Figure 5 and Figure 6 illustrates the advantages of using GSA to compensate Ps and Qs , respectively. The settling time for the active power is at least 13 ms faster when GSA is considered. For the reactive power, GSA reduced the settling time by about 25% at least. (a) (b) Figure 3. The rotor voltages when three-phase to ground fault is applied using GA, PSO and GSA: (a) v’dr and (b) v’qr Figure 4. The rotor angular speed when three-phase to ground fault is applied using GA, PSO and GSA Int J Elec & Comp Eng, Vol. 8, No. 2, April 2018 : 780 – 792 Int J Elec & Comp Eng  ISSN: 2088-8708 (a) 789 (b) Figure 5. The stator active and reactive power when three-phase to ground fault is applied using GA, PSO and GSA: (a) Ps and (b) Qs (a) (b) Figure 6. The rotor voltages when the reference speed is changed using GA, PSO and GSA: (a) vdr (b) vqr 4.2. Variable Reference Speed In this disturbance and as shown in Figure 6, the reference speed m,ref was initially set to 0.7 p.u., after that, the speed is reduced to 0.35 p.u. in one second. Next, the reference speed stays constant for one second and finally increases linearly to the initial reference value. Figure 7. The rotor speed with changing the reference speed using GA, PSO and GSA Gravitational-Search Algorithm for Optimal Controllers Design of Doubly-fed …. (Saher Albatran) 790  ISSN: 2088-8708 From Table 2, a smaller ITAE is associated with the GSA as seen in Figure 6. The settling time of the output variables listed in Table 3 is damped faster by at least 0.38 s when GSA is used. The rotor voltages shown in Figure 7 clarifies that both GA and PSO give almost a similar response while GSA can reduce the maximum over shoot during the acceleration or the deceleration periods. The stator active and reactive power shown in Figure 8 clearly show the advantage of using GSA to damp the response much faster. Using GSA, the system requires 0.47 s to damp the reactive power while it takes 150% of this time when GA or PSO are adopted. (a) (b) Figure 8. The stator active and reactive power under variable reference speed using GA, PSO and GSA: (a) Ps (b) Qs 5. CONCLUSIONS In this paper, GSA is used to design optimal controllers for improving the transient response of a nonlinear DFIG-based WECS. The designed optimal controllers are tested under different common types of disturbances. The results are compared to the results obtained using the two well-recognized global optimization algorithms, which are GA and PSO, to show the effectiveness of using GSA to attain a global optimal solution of the design problem under study. A time-weighted absolute error (ITAE) is used in the objective function to stabilize the system and increase its damping when subjected to different types of disturbances. The comparative simulations confirm the fruitiness of the GSA and its effectiveness to damp the system faster than the GA and the PSO. Mathematically, the considered measured performance indices (ITAE and the settling time) show clearly the advantages of applying GSA over the other two algorithms for different disturbances. REFERENCES [1] [2] [3] [4] [5] [6] Boulahia A, Nabti K, Benalla H.“Direct Power Control for AC/DC/ACConverters in Doubly Fed Induction Generators Based Wind Turbine”, International Journal of Electrical and Computer Engineering (IJECE). 2012; 2(3); 425-432 Anju M, Rajasekaran R.“Power System Stability Enhancement and Improvement of LVRT Capability of a DFIG Based Wind Power System by Using SMES and SFCL”, International Journal of Electrical and Computer Engineering (IJECE). 2013; 3(5); 618-628. Hamane B, Benghanemm M, Bouzid A, Belabbes A, Bouhamida M, Draou A. “Control for Variable Speed Wind Turbine Driving a Doubly Fed Induction Generator using Fuzzy-PI Control”. Energy Procedia. 2012; 18; 476-485. Zamanifar M, Fani B, Golshan M, Karshenas H. “Dynamic modeling and optimal control of DFIG wind energy systems using DFT and NSGA-II”. Electric Power Systems Research; 2014; 108; 50-58. Chakravarty P, Venayagamoorthy G. “Development of optimal controllers for a DFIG based wind farm in a smart grid under variable wind speed conditions”. 2011 IEEE International Electric Machines & Drives Conference (IEMDC). 2011. Niagara Falls, ON, Canada Yang L, Yang G.Y, Xu Z, Dong Z.Y, Wong K.P, Ma X. “Optimal controller design of a doubly-fed induction generator wind turbine system for small signal stability enhancement”, IET Generation, Transmission & Distribution 2010; 4 (1); 579–597. Int J Elec & Comp Eng, Vol. 8, No. 2, April 2018 : 780 – 792 Int J Elec & Comp Eng [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] [18] [19] [20] [21] [22] [23] [24] [25] [26] [27] [28] [29] [30] [31] [32] [33] [34] ISSN: 2088-8708  791 Suryoatmojo H, Zakariya A.M.B, Anam S, Musthofa A, Robandi I. "Optimal controller for doubly fed induction generator (DFIG) using Differential Evolutionary Algorithm (DE)," in Intelligent Technology and Its Applications (ISITIA); 2015;159-164 Vieira J.P.A, Nunes M.V.A, Bezerra U.H, Nascimento A.C. “Designing optimal controllers for DFIGs using a genetic algorithm”, IET Generation, Transmission& Distribution. 2009; 3 (1); 472–484. Sitharthan R, Geethanjali M. “An adaptive Elman neural network with C-PSO learning algorithm based pitch angle controller for DFIG based WECS”. Journal of Vibration and Control, 2016; 23(5), pp.716-730. Hoseinpoor A, Gardoon Piraei A, Asgar M, Mohammadzaki M. “DFIG WIND TURBINES CONTROLL BY PSO ALGORITHM”, Journal of Fundamental and Applied Sciences. 2016; 8 (2); 3638-3646. Mukhopadhyay B, Mandal R. “Application of Artificial Intelligence-Based Techniques in Controlling the STATCOM Used for Compensation for Voltage Dips in DFIG-Based Grid-Connected Wind Power System’. Lecture Notes in Electrical Engineering, 2017; pp.177-189. Bharti, Om Prakash, R. K. Saket, and S. K. Nagar. "Controller Design of DFIG Based Wind Turbine by Using Evolutionary Soft Computational Techniques." Engineering, Technology & Applied Science Research 7.3 (2017): pp-1732. Panda and M. Tripathy, “Optimal power flow solution of wind integrated power system using modified bacteria foraging algorithm”, International Journal of Electrical Power & Energy Systems. 2014; 54; 306-314. Venu Madhav G, Obulesu Y. P. “A New Hybrid Artificial Neural Network Based Control of Doubly Fed Induction Generator”, International Journal of Electrical and Computer Engineering (IJECE). 2015; 5 (3); 379-390. Dahiya P, Sharma V, Sharma R. “Optimal Generation Control of Interconnected Power System Including DFIGBased Wind Turbine”, IETE Journal of Research. 2015; 61 (3). 285-299 Roy R, Jadhav H. “Optimal power flow solution of power system incorporating stochastic wind power using Gbest guided artificial bee colony algorithm”, International Journal of Electrical Power & Energy Systems. 2015; 64 (1); 562-578. Kahla S, Soufi Y, Sedraoui M, Bechouat M. “On-Off control based particle swarm optimization for maximum power point tracking of wind turbine equipped by DFIG connected to the grid”, International journal of hydrogen energy. 2015; 40 (39); 13749-13758. Vepa R. “Nonlinear, optimal control of a wind turbine generator”, IEEE Transactions on Energy Conversion, 2011; 26 (2); 468-478. Yun-tao S, Qi K, Yan-jiao H, Yang Z, De-hui S, Zheng-xi L. "A Hybrid optimal control structure of DFIG Wind Energy Convert System," in Intelligent Control and Automation (WCICA), 2014 11th World Congress on , vol., no., pp.4982-4987, June 29 2014-July 4 2014 Vrionis T, Koutiva X, Vovos N. “A Genetic Algorithm-Based Low Voltage Ride-Through Control Strategy for Grid Connected Doubly Fed Induction Wind Generators”, IEEE Trans. Power Syst., 2014; 29 (3); 1325-1334. Xiang D, Ran L, Tavner P, Yang S. “Control of a Doubly Fed Induction Generator in a Wind Turbine During Grid Fault Ride-Through”, IEEE Transactions on Energy Conversion, 2006; 21 (3); 652-662. Barros L, Mota W, da Silva J, Barros C. "An optimal control strategy for DFIG, " 2010 IEEE International Conference. march 2010; 1727 -1732. Sheikhan M, Shahnazi R, Nooshad Yousefi A. “An optimal fuzzy PI controller to capture the maximum power for variable-speed wind turbines”, Neural Comput & Applic. 2012; 23 (5); 1359-1368. Ruiz-Cruz R, Sanchez E, Ornelas-Tellez F, Loukianov A, Harley R. “Particle Swarm Optimization for DiscreteTime Inverse Optimal Control of a Doubly Fed Induction Generator,” IEEE Trans. on Cybernetics, 2013; 43 (6); 1-9. Feng W, Xiao-Ping Z, Ping J, Sterling M. “Decentralized Nonlinear Control of Wind Turbine With Doubly Fed Induction Generator”, IEEE Trans. Power Syst.2008; 23 (2); 613-621. Qiao W. Vanayagamorthy G.K, Harley R.G. "Design of optimal PI controllers for doubly fed induction generators driven by wind turbines using particle swarm optimization", Intl. Joint Conf. on Neural Networks 2016, pp.1982-1987 Vieira J, Bezerra U, Nunes M, Nascimento A. “Designing optimal controllers for doubly fed induction generators using a genetic algorithm”, IET Generation, Transmission & Distribution, 2009; 3 (5); 472-484. Jeba, V.E.; Ravichandran, S.; Devi, R.P.K., "Optimal control of grid connected variable speed wind energy conversion system," in Energy Efficient Technologies for Sustainability (ICEETS); 393-399, 10-12 April 2013 Rashedi E, Nezamabadi-pour H, Saryazdi S. “GSA: A Gravitational Search Algorithm”, Information Sciences, Vol. 179, No. 13, pp. 2232-2248, 2009. Rashedi E, Nezamabadi-pour H, Saryazdi S. “BGSA: binary gravitational search algorithm”, Natural Computing, 2009; 9 (3);727–745. Balachandar S, Kannan K. “Newton’s Law of Gravity-Based Search Algorithms,” Indian Journal of Science and Technology. 2013; 6 (2); 2013. Ibrahim T.N.T, et al.,“A Brief Analysis of Gravitational Search Algorithm (GSA) Publication from 2009 to May 2013, ” Int. Conf. on Recent treads in Engineering & Technology (ICRET’2014), Feb 13-14, 2014, Batam, Indonesia. Ajami and A. J. Shendi, “Application of UPFC Tuned Based on Gravitational Search Algorithm to Damping of Low Frequency Oscillations,” Int. J. of Advanced Research in Computer Science and Software Engineering. 2013; 3 (7); 2013. Tosun S, Oztürk A, Ertay M. M, Yalçın M.A, Zengin A.“An Approach For Designing Passive Power Filters For Industrial Power Systems By Using Gravitational Search Algorithm, ” Tehnički vjesnik. 2015; 22 (2); 343-349. Gravitational-Search Algorithm for Optimal Controllers Design of Doubly-fed …. (Saher Albatran) 792  ISSN: 2088-8708 [35] Duman S, Guvenc U, Sonmez Y, Yorukeren N. “Optimal power flow using gravitational search algorithm,” Energy Conversion and Management, 2012; 59; 86–95. [36] Boulahia, Adel M, Hocine B. “Predictive Power Control of Grid and Rotor Side converters in Doubly Fed Induction Generators Based Wind Turbine,” International Journal of Electrical and Computer Engineering (IJECE). 2013; 3 (3); 287-293. BIOGRAPHIES OF AUTHORS Saher Albatran received the B.Sc. degree in electric power engineering from Yarmouk University, Irbid, Jordan, in 2005. He received his M.Sc. degree in electric power and control engineering from the Jordan University of Science and Technology, Irbid, in 2008 and his Ph.D. degree in electrical engineering from Mississippi State University, Mississippi State, MS, USA in 2013. He is currently an Assistant Professor with the Department of Electrical Engineering, Jordan University of Science and Technology, Irbid, Jordan. His research interests include integration distributed generation, optimization of power system and motor drives. Muwaffaq I. Alomoush was born in Jordan in 1967. He received the B.Sc. (1990) and the M.Sc. (1994) degrees in electrical engineering from the Jordan University of Science and Technology (JUST), Irbid, Jordan, and the Ph.D. degree (1999) in electrical engineering from Illinois Institute of Technology (IIT), Chicago, U.S. He is currently a professor of the Electrical Power Engineering Department, Hijjawi Faculty for Engineering Technology, Yarmouk University, Irbid, Jordan. His areas of interest are power system control, power system restructuring, economic operation of power system, FACTS applications in restructured systems, and optimization of power system. Ahmed M. Koran received the B.Sc. (2004) in Electrical Power Engineering from Yarmouk University, Irbid, Jordan. He received the M.Sc. (2007) and the Ph.D. (2013) degrees in Electrical Engineering from Illinois Institute of Technology (IIT), Chicago, U.S. and Virginia Polytechnic Institute and State University (VT), Blacksburg, USA, respectively. He is currently an assistant professor and the chairman of the Electrical Power Engineering Department, Hijjawi Faculty for Engineering Technology, Yarmouk University, Irbid, Jordan. His research interests include the modeling, design, and optimization of Power Electronics systems for renewable energy applications. Int J Elec & Comp Eng, Vol. 8, No. 2, April 2018 : 780 – 792