Adaptive coordination control strategy of renewable energy sources, hydrogen production unit, and fuel cell for frequency regulation of a hybrid distributed power system

Owing to the significant number of hybrid generation systems (HGSs) containing various energy sources, coordination between these sources plays a vital role in preserving frequency stability. In this paper, an adaptive coordination control strategy for renewable energy sources (RESs), an aqua electrolyzer (AE) for hydrogen production, and a fuel cell (FC)-based energy storage system (ESS) is proposed to enhance the frequency stability of an HGS. In the proposed system, the excess energy from RESs is used to power electrolysis via an AE for hydrogen energy storage in FCs. The proposed method is based on a proportional-integral (PI) controller, which is optimally designed using a grey wolf optimization (GWO) algorithm to estimate the surplus energy from RESs (i.e., a proportion of total power generation of RESs: Kn). The studied HGS contains various types of generation systems including a diesel generator, wind turbines, photovoltaic (PV) systems, AE with FCs, and ESSs (e.g., battery and flywheel). The proposed method varies Kn with varying frequency deviation values to obtain the best benefits from RESs, while damping the frequency fluctuations. The proposed method is validated by considering different loading conditions and comparing with other existing studies that consider Kn as a constant value. The simulation results demonstrate that the proposed method, which changes Kn value and subsequently stores the power extracted from the RESs in hydrogen energy storage according to frequency deviation changes, performs better than those that use constant Kn. The statistical analysis for frequency deviation of HGS with the proposed method has the best values and achieves large improvements for minimum, maximum, difference between maximum and minimum, mean, and standard deviation compared to the existing method.


Introduction
Hybrid generation systems (HGSs) are often used in a centralized form. However, a centralized HGS delivers a continuous power source, and has various disadvantages, such as high losses because of transmission systems, low effectiveness, increased emission, insufficient power supply reliability, high changing and maintenance costs, unstable peaking power, low resilience, and low sustainability [1]. Therefore, it is an urgent requirement to transfer to a distributed HGS, which has many merits because of the utilization of renewable energy sources (RESs). Advantages are such as high effectiveness, low emission, generation and delivery of electricity being more reliable and having high quality, low changing and maintenance costs, high resilience, and sustainability [1]. However, the unpredictable changes in the weather, such as wind speed and irradiance, lead to intermittent generation from RESs (e.g. PV and wind systems) [2]. Thus, the uncertain characteristics of PV and wind systems, and the probabilistic behavior of the loads, lead to fluctuations in system frequency [3]. In order to overcome these challenges, load frequency control (LFC) is mandatory for distributed HGS. In the literature, several methods are proposed to mitigate system frequency fluctuations. A black hole optimization algorithm is applied to fine-tune the cascaded proportional derivative-proportional integral derivative (PD-PID) controller parameters for the frequency regulation of a distributed HGS [4], while an adaptive differential evolution technique based on PD with filter plus (1 + PI) is proposed in [5]. Reference [6] presents a modified equilibrium optimizer process that relies on an interval type-2 fuzzy-PID controller for LFC of a distributed HGS. A modified moth swarm algorithm (mMSA) with a hybrid fuzzy PD-PI controller to achieve the frequency stability of a distributed HGS is described in [7]. In addition, a feedforward fractional-order PID strategy based on a harmony search algorithm is applied in [8] to regulate system frequency considering the integral of time multiplied squared error (ITSE), while the chaotic crow search algorithm is used to optimize the hybrid fuzzy proportional derivative-tilt integral derivative controller parameters to control the frequency of a standalone microgrid [9].

Open Access
A bacterial foraging optimization algorithm-based PID controller with a derivative filter is used for frequency control of a distributed HGS [10], while the researchers in [11] propose a compound control strategy based on model predictive control and a distributed leaderfollowing consensus control strategy. In these control algorithms, the distributed energy resources and energy storage technologies are used to support the power system frequency [12,13]. In [14], the PI + clegg integrator (CI) reset control technique is proposed for the secondary frequency control in an islanded microgrid, considering RESs/load changes. A cascaded PI-PD controller is tuned by the salp swarm optimization technique for LFC of a standalone microgrid integrated with electric vehicles (EVs) [15].
In [16], a comparison is presented between the proposed fractional-order (FO)-fuzzy-PID controller method and different types of control techniques such as optimized FO-Fuzzy-PID, PID, FO-PID, and Fuzzy-PID. The droop controller and PID controller are used with superconducting magnetic energy storage (SMES) to inject power to the microgrid during sudden load changes, thus achieving the frequency stability aim [17]. The proposed control strategy in [18] depends on two main control actions. The first is injecting active power from the battery to reach a constant power between the load and generation where the generated power from RESs is intermittent. The second control action is related to the frequency deviation, in which the proposed method is used to damp the frequency fluctuation and to compensate for the low-inertia of the system. The issue of reducing the lifetime of a battery energy storage system (BESS) (which helps in damping frequency deviation through a large amount of charging and discharging) is addressed in [19], where an adaptive droop control method is proposed between the BESS and FC through the controlling of the power-sharing between FC and battery.
One of the newest optimization algorithms is the grey wolf optimization (GWO) algorithm based on the metaheuristic optimization technique [20]. The GWO algorithm is used in many applications of electrical power systems because of its simplicity and easy applicability [21]. Moreover, the GWO has several advantages such as flexibility, fewer algorithm parameters, comprehensiveness, and fast programmability features. The PID controller gains are used with the integral time absolute error (ITAE) as inputs to the control process that depends on the GWO to optimize the value of the PID gains, in order to achieve the secondary frequency control of an isolated multi-microgrid [22]. The particle swarm optimization (PSO) and GWO techniques are applied with adaptive fuzzy logic control based on the PI controller to obtain the optimal values of membership functions (MFs) in the regulation and control of microgrid frequency [23].
The coordination control strategy between RESs and FC-based energy storage systems (ESS) is presented in [24], where an adaptive fuzzy PID controller is proposed based on a simplified GWO algorithm for the LFC of a distributed power generation system considering generation rate constraints and time delays to include nonlinearity features. It also presents a comparison between several optimization algorithms including the simplified GWO algorithm to show its superiority over other methods. One of the most important applications of a hydrogen production unit is presented in [25] to overcome the large-scale power curtailment resulting from the connection between RESs and large-scale power systems. The hydrogen production unit is suitable for such an application, since it can have large-capacity and long-term absorption of electric energy.
The presented coordination control strategy has the following levels. The first derives energy conversion efficiency from electrical to hydrogen energy by a hydrogen production unit. The second takes into account the storage state of charge and the instantaneous power in the control process. The maximum power point tracking of a PV system is designed and included in the control process. The Bahr AL-Najaf area is used as a case study that utilizes a PV and wind system with fuel cells [26]. The chaotic PSO is used in [27] to optimize the parameters of FO fuzzy control. The proposed method is compared with PID and the integer-order fuzzy PID controller, considering the linear and nonlinear operating systems, to show its effectiveness over the mentioned control methods. The management strategy is to merge with the DC grid in [28] to optimize the bus voltage and coordinate the power-sharing between the energy sources during different loading conditions with variable solar irradiance, considering the state of FC. The presented method is based on DC bus regulation through supercapacitor energy, and a PI controller is used to reduce the fluctuations in the DC bus voltage. The energy management strategy relies on the optimal control for a battery-supercapacitor in HGS, and is targeted to define the adaptive gains to state the reference power for the battery and the supercapacitor by using the slap swarm algorithm and load following control technique [29].
The high penetration of integrated RES with EVs is aggregated as a virtual power plant to reduce the cost of power balancing between the generation and demand sides, and the proposed control framework focuses on the control of complementary generation based on lifelong learning [30]. In this control strategy, lifelong learning represents the primary layer, while the secondary layer is set in line with the ascending order of the regulation costs of the available reserve resources. In addition, the tertiary layer is represented by the coordination control for RESs, and an online optimization is applied considering the EV charging demands. In [31], an optimal dispatch for automatic generation control is proposed to distribute the AGC command between generation units in a hybrid system with RESs and conventional resources. To do this, an adaptive distributed auction-based algorithm (ADAA) is used instead of the multiple fixed swap sizes of the auction-based algorithm, where ADAA has a fast convergence speed (i.e., low communication latency) and model-free features.
The rest of this paper is structured as follows: the relevant literature is reviewed in Sect. 2, to highlight how the proposed method differs from existing ones. Section 3 presents the model of the HGS, while Sect. 4 discusses the proposed adaptive coordination control strategy. The simulation results are presented in Sect. 5, and finally, Sect. 6 concludes the study.

Literature survey
In [32], a super-twisting algorithm (STA) technique is proposed to improve the fluctuation of bus voltage, the current circulation, and the power distribution rationality. However, the impact of the wind system and high load variations are not investigated, and the fraction factor of the power-sharing between RESs and FC is also not considered. A two-time scale energy management system is presented in [33] to deal with the FC and electrolysis specifications, but the high load variation and the fraction factor consideration are not presented with the proposed method. The P-f droop control strategy is used in [34] for an energy management system between the AC microgrid components, considering their electrical characteristics, although the fraction factor which is used to define the power-sharing between RESs and FC is not studied. In [35], a linear quadratic regulator with an integral controller is applied to control the frequency of a microgrid (MG). However, the effectiveness of the proposed method is validated under small load perturbation, and the fraction factor is not taken into consideration. References [25,36] propose a decentralized coordination control strategy among generation, ESS, hydrogen production unit, and FC in island DC and AC microgrids, respectively. However, both studies fail to take the electrical characteristics of the FC and the fraction factor into consideration. A stochastic model is proposed for the coordinated scheduling of RESs in [37], but without including state-of-charge (SoC) protection and overpower protection studies. In addition, in the decentralization approach, the fraction factor is not considered. The most relevant studies are compared in Table 1.
Based on Table 1, several research gaps can be identified: (i) SoC protection and overpower protection; (ii) decentralization approach; (iii) electrical characteristics of the FC; (iv) high load variations; and (v) the fraction factor coordination to state the power-sharing between RESs and FC.
Based on the above discussion and the presented studies, all the studies have either used a constant value for the fraction of the total power of RESs (Kn) or not considered. Therefore, this study proposes an adaptive coordination control application to control the power used from the RESs for energizing the FC-based ESS. The control approach proposed here is based on controlling the value of Kn that depends on the GWO algorithm to obtain the optimal gains of the PI controller. These are used to define the value of Kn with different loading conditions and uncertainties of PV and wind power generation. Moreover, it defines the power extracted from the RESs injected into the FC to achieve the best benefits from the RESs and FC, while fulfilling the frequency control aim.
In order to clarify the disparities among the different studies, Table 2 provides a comparison among the investigated studies and the present work. The main contributions of this study are as follows: • Benefits from RESs and FC-based ESSs using an adaptive coordination control strategy are obtained. • The proposed adaptive coordination control method controls the value of Kn and consequently the power extracted from RESs to be used with FC-based ESSs. • The GWO algorithm is used with the proposed adaptive coordination control method to adjust the PI controller gains of fraction factor Kn with different loading conditions and uncertainties of RESs powers. • The advantages of the proposed coordination strategy are proved and established over existing studies that are based on the constant value of Kn. • A comparison between the proposed method and the existing methods is presented. • The effectiveness and robustness of the proposed method are proved using different loading conditions and uncertainties of RESs/loads power. • The frequency fluctuation is stabilized successfully by the newly proposed control method.

Hybrid generation system (HGS) model
The studied HGS is presented in Fig. 1. This shows different types of energy sources: wind generators (WGs), a photovoltaic (PV) system, an aqua electrolyzer (AE), fuel cells (FCs), a diesel generator (DG), electric vehicles (EVs), BESS, and a flywheel energy storage system (FESS) [38]. All these components are represented by transfer function (TF) models, where the values of the gain (K) and time constant (T) are listed in Table 3. The LFC-PID controller is designed using the mMSA, as discussed in [7].

Wind generator
The studied wind generation system consists of three wind turbines, and the power produced by each is P WG . The well-known equation describing the mechanical power of a wind generator (P m ) is presented in (1), where C p is the coefficient of performance and is related to β (blade pitch angle) and λ (tip speed ratio), ρ is air density, A is turbine swept area, and v wind is wind speed [38].
The nonlinearities and uncertainties are considered to be the main features of the wind generator system, e.g., the pitch control system. The wind generator system is represented by the TF model as in (2), where K WG and T WG represent the gain and time constant, respectively, and P m and P WG stand for the input and output power of the wind generator system, respectively. (1)

Photovoltaic model
A PV system is based on converting the sun's radiation into electricity, and the PV power (P PV ) can be calculated by (3), where A PV is the PV array area, φ is the solar radiation,η is the PV cell conversion efficiency.
The TF of the PV system is defined in (4) [38], where K PV is the gain and T PV is the time constant of the PV system, while �φ and P PV are the input and output of the PV system, respectively.

Aqua electrolyzer (AE)
The power generated by the PV and wind generation systems is used by an AE to generate hydrogen, while the generated hydrogen can be injected into an FC-based ESS to produce power. The TF of the AE is specified in (5) [38], where K AE and T AE are its gain and time constant, respectively. The term [(1 − K _n) × P_RESs] represents the input of the AE and P AE is its output.
(3)  Fig. 1 The studied HGS with the proposed coordination control strategy Table 3 Nominal parameters of the studied HGS

Fuel cells
The studied system has two identical FCs to provide high efficiency by sharing the power, and to increase the system reliability in case one fails. As the FC here is the main energy storage that is integrated to make coordination with the RESs, it is important to always keep them in the system. The TF used in this study is based on the first-order presented in (6) [38], where K FC and T FC are the gain and time constant of the FC system-based ESS, respectively, while P AE and P FC are the input and output of the fuel cell system, respectively.

Diesel generator
The diesel generator can be used as an independent power source to minimize the mismatch between generation and demand. The TF of the diesel generator is given in (7) [38], where K DG and T DG are the gain and time constant of the diesel generator, respectively, while u DG and P DG are its input and output, respectively.

Electric vehicles
The EVs can operate as supply in the discharging mode and as load in the charging mode, so their performance helps to mitigate system fluctuations. The TF of EVs in the first-order is specified in (8) [38], where K EV and T EV are the gain and time constant of EVs, respectively.

Battery energy storage system
The BESS plays an essential role in preserving system frequency stability, as it can operate either as a load or a power source. The TF of BESS is presented in (9) [38], where K B is the gain of BESS and T B is its time constant.

Flywheel energy storage system
The FESS aids in keeping the power system as stable as possible during different operating conditions. The TF of the FESS is specified in (10) [38] where K F is the gain of FESS and T F is its time constant.

Power system model
The dynamic power system model has one input and one output. The input is the change in the power (the difference between the load and generation), while the output changes the system frequency as presented in (11) [38], where D and J are the damping factor and the virtual inertia factor, respectively.

The proposed adaptive coordination control strategy
In this study, the algorithm of GWO is programmed using the MATLAB program. The PI controller of Kn is optimized, as shown in Fig. 2. The main role of the optimization is to achieve the adaptive control method of the RESs, the AE for hydrogen production and the FC-ESS, and enhance the frequency of the power system in different loading conditions in addition to the main frequency control. The different loading conditions will be presented in the results section. The objective function used in this study is the ITAE shown in Fig. 2. This integrates the multiplication of the absolute error (i.e., system frequency deviation ( f )) and the time as: where t sim is the simulation time.
Recently, one of the most recommended optimization techniques has been GWO, which is considered as a meta-heuristic technique [20]. Wolves have a hierarchy in strict dominance when they live in a group, which typically has 5 to 12 members. The leader of this group is called alpha (α), who has a strong aider called beta (β), which helps the main leader make the right decision. The omega (ω) is considered the lowest position wolf, which acts as the scapegoat. The rest of the wolves are called delta (δ), and all of the deltas follow alphas and betas. All of these roles are presented in Fig. 3. The GWO method has three main steps to reach a prey: (i) search for the prey and approach it, (ii) surround the prey and limit its motion, and (iii) attack and hunt down the prey. The surrounding of the prey can be presented in arithmetic form as in (13), where P i is the position of the GW, P pi is the position of prey, and S is the distance between them. A and C are vectors that are computed from Eqs. where where − → n 1 and − → n 2 express the random numbers between zero and one.
Equation (15) presents the factor a that reduces linearly from two to zero when increasing the number of iterations. The diverging technique is used for searching a prey location where |A| > 1 , while the convergence technique is used for getting a prey where |A| < 1 . The hunting process uses α which is then followed by β and δ as presented in (18)- (20). The A and C factors help the GWO method to keep away from the stagnation event.
The flowchart of the proposed adaptive coordination control strategy based on the GWO technique is shown in Fig. 4, to achieve the coordination between RESs, aqua electrolyzer for hydrogen production, and FC-ESS, thus enhancing the frequency stability of HGSs. The overall proposed adaptive coordination control strategy is organized as follows: Step 1: Measure and collect the various signals, and feed them into the adaptive coordination control process. These parameters include the generation active power, load active power, and the parameters of the GWO algorithm.
Step 2: The computation of the frequency deviation is presented in this step from the known difference between the load active and generation active power. The power mismatch between the load and generation plays a vital role in power system frequency stability and defines the value of the frequency deviation. The frequency deviation must be within permissible limits based on the signal of an adaptive coordination control strategy that defines the fraction of total power generation of RESs (Kn) and the power fed into the AE for hydrogen production and FC. When the frequency deviation is out of limit, the absolute value of the frequency deviation is taken and used in Step 3.
Step 3: The GWO algorithm is described in this step, where the best value of the Kn can be obtained by applying the GWO to get the optimal value of the PI controller. The value of Kn is responsible for coordinating the RESs and the AE of hydrogen production and the FC.
Step 4: The best value of Kn is used in the equation presented in Step 4 and also shown in Fig. 1. The power extracted from RESs is adaptively used with an AE, and the remaining power is fed into the HGS. By using this adaptive coordination control strategy, the power of FC and RESs can be optimally distributed to help preserve the system frequency within limits.
The fitness value and gain factors of K P and K I are listed in Table 4. Fig. 3 The hierarchy of grey wolves  (13) to (20) Kn is the best value Step 1 Step 2 Step 3 Fig. 4 Flowchart of the complete proposed method

Results and discussion
The studied system shown in Fig. 1 is simulated using the MATLAB/Simulink software. The simulation results comprise four cases of different loading conditions to indicate the performance of the proposed method compared to the existing work in [7] that used a fixed value for the coordination gain of Kn. At the same time, the frequency regulation control is based on mMSA in both studies to make a fair comparison between the two systems. Therefore, a better coordination between the RESs and the ESS-based fuel cell can be achieved, where the PID controller based on mMSA damps the frequency regulation in both studies. The PI controller based on GWO is proposed to determine the value of Kn that defines the injected power into the AE to power the FC. This case is called "Proposed PI-GWO of Kn" and the existing case is referred to as "Constant value of Kn". The power of the wind generator and PV systems used during all cases at different loading conditions is presented in Fig. 5. It is clear that the RESs profile has random output characteristics as a reflection of the change in environmental conditions which lead to fluctuation in system frequency. The fluctuations of output power from the wind and PV systems are overcome by a suitable control approach. Different cases are tested with variable load conditions to examine the proposed control method.

Case 1
In this case, the random load variation is a combination of two load types: one being an industrial load and the other a domestic load, as shown in Fig. 6, where the total load profile is the sum of the two. The total load changes from 0.045 pu to 0.307 pu, while the generated power from the wind varies between 0.025 pu to 0.225 pu (see Fig. 5a), and the PV between 0.001 pu to 0.1 pu (see Fig. 5b). The frequency deviation in this case study is presented in Fig. 7. It is evident that at 400 s, the proposed PI-GWO of Kn controller mitigates the frequency fluctuation better than that with the constant value of Kn. The improvement in the system frequency results from the impact of the Kn controller, where the power of the FC is changed, as presented in Fig. 8a. The power injected into the AE is controlled by following the value of Kn as shown in Fig. 8b. The response of the controller is evident in Kn as presented in Fig. 8c, where the proposed method varies the value of Kn according to the frequency state, whereas the existing work uses constant Kn regardless of the frequency deviation in the system. The variation of input   power from different power sources of HGS is presented in Fig. 9, in which the proposed work provides smoother and lower fluctuations of the power from the different sources than the existing work during the variation in the power of RESs and load. This leads to reduced variations of the generated power from the BESS, FESS, EV, and DG, as is evident in Fig. 9. The above discussion proves that the proposed control method is suitable for making the HGS more stable and robust in the case of large system disturbances.

Case 2
In this case, the load demand is varied as shown in Fig. 10, while the power profile of wind and PV systems is the same as in Fig. 5. The frequency deviation is given in Fig. 11. It is observed from the curve that the HGS frequency is damped using the proposed method despite the uncertainties of the power generated from RESs and the load variations. The responses of the FC power and the power injected to the AE, and the Kn response are shown in Fig. 12, which indicates the effect of the proposed control. The variations of power from the generation sources and ESSs are presented in Fig. 13. By applying the proposed control method, the power from the variable sources can be better smoothed than with the existing method.

Case 3
The proposed control method is examined with multisharp variation in the load demand curve, as shown in Fig. 14, considering the same profile of the wind and PV generation systems as in Fig. 5. The frequency fluctuation due to the variable load demand and intermittent RESs is minimized by applying the proposed coordination method using the adaptive gain Kn (Fig. 15). The proposed method has successfully minimized the overshoot, undershoot, and steady-state error of the system frequency during the overall simulation time. The impact of the proposed control technique on HGS is indicated through the performance of the FC power energized by the AE. This is controlled by the Kn value that is based on the PI-GOW technique, as displayed in Fig. 16. Also, the power variations from different generation sources are presented in Fig. 17 with the two control systems.

Case 4
In this case, the proposed control is tested by multichanges in the load demand profile as displayed in Fig. 18 to inspect the quality of the control. In the presence of the proposed control method, the frequency deviation is kept approximately constant (around zero) during the variations of the load demand and power generated from the RESs. In contrast, in existing work with a constant value of Kn, the frequency deviation varies between 0.041 pu and -0.03 pu, as presented in Fig. 19. This improvement is achieved by using a suitable control of Kn, which affects the FC operation as shown in Fig. 20a through controlling the power required by the AE as shown in Fig. 20b. The values of Kn with the proposed method and existing work are compared in Fig. 20c. These values are responsible for stating the different power used from the RESs to energize the FC. The BESS, FESS, EV, and DG achieve the balance between generation units and loads in the HGS as displayed in Fig. 21. The smaller fluctuations in the generated power are the results of the proposed control method. The proposed control method and existing work in [7] are further compared in Fig. 22 using the statistical analysis of the frequency deviation presenting computational results of minimum (Min), maximum (Max), the difference between maximum and minimum (Max-Min), mean, and standard deviation (Std). The maximum and minimum values are employed to define the overshoot and undershoot of the frequency response during the different loading conditions, respectively. The difference between them is used to express the range of changes, and if this indicator is near zero, it means that the frequency deviation is a near-constant and stable value.
The mean value reflects the frequency response during the whole simulation time, while the standard deviation expresses how much the frequency response differs from the mean value of the frequency response curve. A comparison between the proposed work and the existing work for the frequency deviation response is presented to evaluate the effectiveness of the proposed method.
The proposed method in Case 1 reduces the undershoot by 82.2% compared to the existing work, and the overshoot by 32.4%. The proposed method has the lowest value of the difference between the maximum and minimum values, mean, and standard deviation. In Case 2, the improvements reach 93.4% and 87.8% in the undershoot and overshoot, respectively, compared to the existing work. The mean, standard deviation values, and the difference between the maximum and minimum validate the proposed control method. These indicators have smaller values than those in the existing work.  In addition, the importance of the proposed control technique compared to the existing work is observed in Case 3, where it improves the lowest minimum value by 78.6% and lowest maximum value by 72.1%. Also, the other indicators, including the difference between the maximum and minimum, mean, and standard deviation achieve enhancements of 76%, 100%, and 91%, respectively. The proposed PI-GWO of Kn control method in Case 4 makes the system more stable and reliable through the improvements of 91.7%, 85.7%, 89.2%, 98.4%, and 95.4% for minimum, maximum, the difference between the maximum and minimum, mean, and standard deviation, respectively.

Conclusion
This paper presents an adaptive coordination control method for optimal operation of the RESs and fuel cell  system. The proposed method is based on the PI-GWO controller to set the power-sharing between RESs and FCs, and therefore to stabilize the power system frequency with different loading conditions. The studied system has multi-energy sources, including intermittent PV and wind generation systems, leading to more fluctuations in the power system frequency. The parameters of the PI controller are tuned using the GWO technique to obtain optimal controller values. The proposed control method is examined with four cases with different load profiles, and the effectiveness and robustness of the proposed method are validated by comparing it with existing work. The simulation results validate the effectiveness of the proposed method compared with the existing one. The key findings can be summarized as follows: • The proposed control approach for HGS constructively mitigates the frequency fluctuations and improves the power exchange between the genera-tion and load sides under different generation and load conditions. • The variable Kn effectively maintains the load/frequency at an approximately constant level during the assessment tests. • With the proposed coordination control strategy based on PI-GWO of Kn, the variability of Kn shows a superior performance over the constant Kn in terms of minimum, maximum, the difference between the maximum and minimum, mean, and standard deviation of the frequency deviation. • The computational results of the frequency deviation are listed and discussed, and prove the superiority of the proposed coordination control method.