Design and Implementation of Frequency Controller for Wind Energy-Based Hybrid Power System Using Quasi-Oppositional Harmonic Search Algorithm

: An innovative union of fuzzy controller and proportional-integral-derivative (PID) controller under the environment of fractional order (FO) calculus is described in the present study for an isolated hybrid power system (IHPS) in the context of load frequency control. The proposed controller is designated as FO-fuzzy PID (FO-F-PID) controller. The undertaken model of IHPS presented here involves different independent power-producing units, a wind energy-based generator, a diesel engine-based generator and a device for energy storage (such as a superconducting magnetic energy storage system). The selection of the system and controller gains was achieved through a unique quasi-oppositional harmony search (QOHS) algorithm. The QOHS algorithm is based on the basic harmony search (HS) algorithm, in which the combined concept of quasi-opposition initialization and HS algorithm fastens the proﬁle of convergence for the algorithm. The competency and potency of the intended FO-F-PID controller were veriﬁed by comparing its performance with three different controllers (integer-order (IO)-fuzzy-PID (IO-F-PID) controller, FO-PID and IO-PID controller) in terms of deviation in frequency and power under distinct perturbations in load demand conditions. The obtained simulation results validate the cutting-edge functioning of the projected FO-F-PID controller over the IO-F-PID, FO-PID and IO-PID controllers under non-linear and linear functioning conditions. In addition, the intended FO-F-PID controller, considered a hybrid model, proved to be more robust against the mismatches in loading and the non-linearity in the form of rate constraint under the deviation in frequency and power front.


Introduction
A number of studies have confirmed that the availability of electricity to communities can deliver several socio-economic gains such as a superior education system, effective business models and better healthcare opportunities [1,2]. The energy crisis in developing countries affects nearly 2 billion people who have no access to electric energy due to the absence of an electric grid in their regions. Moreover, electricity in underdeveloped regions lacks reliability-this not only means that the lights sometimes flicker, but also that the complete region can be left in darkness for hours and even days on end. The solution lies in the existence of the grid; however, the establishment of grids for many of these locations is impractical as most of the population is scattered, living in a harsh landscape or both. Hence, the population in these regions generally have to depend on diesel engine generators (DEGs) for electricity [3]. However, this is a very costly and complex affair as it involves huge costs such as transportation costs for the fuel and other associated complications related to energy such as oil price escalation, global warming, global environmental challenges, etc. These challenges along with the energy requirement demand make small stand-alone renewable energy systems (RESs) essential [4].
RESs (such as wind and solar) are the only substitute for fossil fuels. At present, solar photovoltaic and wind energy systems can provide electricity to a majority of individual systems. However, when compared, the energy conversion proficiency of wind energy systems is higher than that of solar photovoltaic systems. Consequently, wind energy systems have an advantage over solar photovoltaic systems [5].
Wind energy-based generators (WTGs) produce power with shifting frequencies and magnitudes due to irregular attributes from wind energy [4] that interrupt the reliability and stability of the local network [6]. However, these problems related to WTGs may be prevented by hybridizing WTGs with DEGs. The local network undertaken in the current study is categorized as an isolated hybrid power system (IHPS) [7] since it is not connected to the grid.
In the case of WTG failure, the DEG takes control whenever there is a mismatch in the power generation and load demand for IHPS. In a DEG, the fuel ejection system may be fueled by the operational load and the switching frequency. The operation of DEG is not satisfactory when there is frequent switching action and lower load demand. The solution lies in the utilization of energy storage devices (ESDs) along with WTGs and DEGs in IHPS [6]. In conditions with a lower speed of wind or load peak demand, the unconsumed active energy stored by the ESDs (flywheel energy storage system, battery energy storage system, compressed air energy storage, supercapacitors, superconducting magnetic energy storage (SMES), etc.) is made available [8]. However, from the available ESDs catalogue, SMES have the quickest and extremely adaptable control functioning. In [9], SMES is employed for suppressing the power and frequency due to load perturbation in IHPS. A potent design of SMES for an IHPS with a wind-diesel generation unit has been suggested with the standard lag/lead compensator of predefined configuration [10]. Several studies have been carried out that prove the potential of SMES units for IHPS application and hence, a controlling strategy that will synchronize with distinct energy sources for accurate functioning of the system is necessary.
Currently, the primary problem of the IHPS lies in the difficulty associated with the process of controlling the wind speed deviations of a wind turbine (as its power output holds the cubic relation to the wind velocity) and the demanded load. Largely, classical controller tuning is engaged for IHPS controlling, such as proportional-integral (PI) controller, proportional-derivative (PD) controller, proportional-integral-derivative (PID) controllers, etc. However, the sensitivity of the classical controllers concerning load variations, system parametric variations, etc., is high. Thus, in order to resolve this problem associated with controlling the action, wind energy conversion is generally carried out utilizing a fuzzy logic controller (FLC) [11], as proposed by Zadeh [12]. However, the FLC has design deficiencies and still the standards are not well-defined for (a) membership functions (MFs) design, (b) rule base pattern, (c) number of considered linguistic divisions, (d) defuzzification methodology, and (e) appropriate inference mechanism. Generally, in designers' experiences, the method of hit and trial is utilized for a selection of these parameters. Therefore, new practices for resolving these complexities of FLC design are required. In [13], the PID controller is hybridized with FLC for improving its utilization capabilities and controlling the action. Carvajal et al. [14] discussed a controller structure combining the PID controller with FLC and specified its capability of handling systems with non-linear, higher-order uncertainty, as well as complex systems. Furthermore, the Energies 2021, 14, 6459 3 of 23 differentiation and integration (differ-integral) operators provide additional flexibility in design with the combination of fractional order (FO), proving its competency as an improved estimator for the control signals compared to the integer order (IO). The primary idea for the insertion of the input and output FO operators of the FLCs is based on the heuristic understanding of a confident rate of change in error. This degree of change in error as perceived through the operator and a comparable activity can be carried out and cannot be invariable in nature [15].
In recent years, the demand for fractional calculus-based control systems amongst researchers has increased manifold because of its added features of being flexible and high-functioning. Podlubny [16] expanded the logic of the classical PID controller by utilizing the concept of FO differ-integral to the FO-PID controller by proposing two new scales of freedom: the orders of the integrator (λ) and differentiator (µ). Owing to the features of high flexibility and greater adjustment property, the FO-PID controllers find their application within the realm of various control schemes for controlling energy in wind-driven generators [17], nuclear power plants [18] and water-driven generators [19]. Furthermore, the FO-based fuzzy-PID (F-PID) (FO-F-PID) controller is broadly applied for controlling systems with non-linear, time-varying and vague characteristics. To handle the oscillatory FO processes, the FO-F-PID control with the tuning process based on the optimal time domain is projected in [20,21].
From the literature survey, it can be stated that among all the controllers, from the perspective of superior performance and robustness, the FO-F-PID controller has an edge. An FO-based control structure for the IHPS model is explored in the present study, inspired by the diversified applications of FO control systems based on computational intelligence for various power systems with the automatic regulation of voltage [22], frequency deviation suppression on the two-area system under load perturbation [23] and control of frequency in microgrid configurations [24].
In [25], the HS algorithm's amended version is proposed for upgrading the convergence shape by introducing a reasonable optimal global solution in real-time. The swarm intelligence perception has been exploited to realize a global best HS algorithm [26]. Mahdavi et al. proposed amending its crucial variables dynamically to achieve a better HS algorithm [27]. In [28], Banerjee et al. applied the opposition-reflected learning concept for the opposition-based HS (OHS) algorithm in the optimizing problem for the compensation of reactive power in IHPS. Similarly, by applying the concept of quasi-oppositional learning to the fundamental algorithm of HS, the developed algorithm is the quasi-oppositional harmony search (QOHS) algorithm, which is utilized in IHPS to stabilize the deviation of frequency in [29]. Therefore, in the current study, the newest and effective QOHS algorithm [29] is implemented with the view of controlling the deviation in power flow alongside the variation by optimal tunning of the vital optimizable variables of the considered IHPS on its installed controllers and ESD.
Given the aforementioned literature survey, the core motivation of the current study evolves with the subsequent considerations.
(a) For the end-users, an alternative source of electricity may be the modelling of an IHPS. (b) For delivering the uninterrupted power supply for the ultimate consumer with quality, the appropriateness of SMES as an ESD may be investigated. (c) In designing the controller, the effectiveness of both FO calculus and IO calculus may be studied and the most suitable one may be adopted. (d) A new modified HS algorithm considers the concept of quasi-oppositional learning, i.e., the QOHS algorithm, recognized for efficient optimization process in regulating the deviations of frequency and power. This may be achieved by appropriately optimizing the vital optimizable variables of the IHPS system under study.
Therefore, the major contributions of the study are as follows. The present study makes an effort to: (a) Improve the stability and reliability of IHPS (consists of WTG, DEG and SMES) for rural as well as urban areas with nil effects on the environment. (b) Enhance the quality of governor-regulated DEG supply power and the pitch-controlled WTG power output together with power obtained from SMES. (c) Implement the theory of fractional calculus for controller design to regulate its working principle. (d) Enhance the objective function-based convergence profile while employing the QOHS algorithm (a new HS algorithm variant, inspired by the concept of quasi-oppositional learning) for optimizing the tuning purpose of the optimizable variables of the controllers with additional vital variables of the system under study. (e) Evaluating the efficacy of the system while dealing with ambiguities, random load demand deviation and stochastic variation in step. (f) Analyze the QOHS algorithm and controllers' robustness. Furthermore, the system robustness is also verified with load mismatch variation and the insertion of rate constraint-type non-linearity.
The proposed work is presented in a flow chart manner in Figure 1, representing the whole idea and work contained in the manuscript in a step-by-step manner. The manuscript is organized as follows. Section 2 describes important elements of the IHPS. In Section 3, architecture and the concept of FO and IO controllers are explained. The details of the optimization task and objective function are illuminated in Section 4. In Section 5, the base variant and adaptation incorporated in the basic HS algorithm to form the QOHS algorithm are demonstrated. The analysis of the simulation results is presented in Section 6. Finally, in Section 7, the present work is concluded.

System Model
In the studied IHPS, during its regular operation, undergoes a minor deviation in the load demand. Hence for the analysis of the IHPS model linear model has been undertaken for dynamic illustration [10]. Consequently, for the studied IHPS (a wind turbine generator (WTG), a DEG and an ESD (a SMES device)), the schematic block diagram of the considered model based on transfer functions (TF/TFs) is depicted in Figure 2 with the installed controllers. The rated capacity of 150 kW is chosen for both the WTG and the DEG respectively and the nominal system parameters are furnished from [28][29][30].
The variations in frequency and power in this studied IHPS, are reduced by the collective work of the SMES, the DEG speed governor and the devoted support of the WTG pitch controller. For the considered IHPS arrangement (Figure 2), the overall variation in the entire power output (∆P TOTAL ) is evaluated using Equation (1): In Equation (1), ∆P TOTAL , ∆P DEG , ∆P WTG , ∆P SMES and ∆P LD represent the change in total power, DEG power, WTG power output, the power output of the SMES and input load demand, respectively.

Wind Turbine Generator (WTG)
The WTG-generated power is governed by the wind speed available to the wind turbine. The wind speed is subjected to variables involving maximum, rated, or zero output power [31], depicted in (2): where P rat is the rated output power, V c_in is the wind cut-in speed, V r is the wind-rated speed and V c_out is the wind cut-off speed. When the wind speed changes, a standard WTG generates power that varies from V c_in to V c_out that leads to a linear increase in the generated power [31]. The pitching action creates non-linearity to support the WTG speed. The WTG is rationalized to a first-order system and denoted as the TF from (3)- (7): where K P , K FC , K PC , K PA and K PB are the gains of the data fit pitch response, the hydraulic pitch actuator, the pitch control, the blade characteristic and fluid coupling, T PA (pitch controller time constant), T PB (hydraulic pitch actuator time constant) and T PC (data fit pitch time constant), respectively, while the ∆F represents the frequency deviation of the system. Hence, the value of ∆P DEG can be estimated as in Equation (1).

Diesel Engine Generators (DEG)
For the DEG model, the speed of the generator along with the mechanical output power from the engine are described as the function of fuel consumption rate (FCR) [32]. The FCR of the DEG can be indicated as (8): where α GED (per kWh) and β GED (per kWh) represent the coefficients from the fuel consumption curve and Q(t) is the fuel consumption of DEG on an hourly basis, made available from the maker. Furthermore, P(t) GED.gen (kW) is the DEG-generated power and P(t) GED.rat (kW) is the DEG-rated power. Here, it should be noted that non-linearity in a DEG exists by virtue of the dead time, which is the time-varying nature between the mechanical torque production and fuel injection. The relation between engine mechanical power and fuel consumption is illustrated in [33] where a simple first-order TF model for DEG has been considered, as in Figure 2, and may be described as: where K DA is the gain of the speed governor, T DA , T DB and T DC are the speed governor time constants and T DT is the turbine time constant. Hence, the calculation of ∆T GD in Equation (1) can be carried out.

SMES Configuration
Throughout the regular functioning, charged SMES establishes the value through the system. In case of an unexpected surge in the power requirement, SMES instantly releases the energy stored to accomplish the power demand change. For this duration, it begins the charging instantly, absorbing a certain percentage of the energy which is a surplus to the system [34]. The firing angle (α) for the converter is furnished by the DC voltage (E L ) developed at the inductor. The firing angle α changes constantly in a specific range of negative and positive values. When the losses from the converter and the transformer are neglected, the expression for DC voltage may be stated as: where I d denotes the inductor current (kA), V d0 denotes the circuit bridge maximum voltage (kV) and R C denotes commutating resistance (Ω). E L is constantly monitored with the SMES control loop input signal [34]. Figure 3 shows the SMES unit block diagram. Hence, the variation of E L and I d may be expressed as: where ∆E L represents incremental converter voltage deviation (kV), ∆I D represents SMES current incremental change (kA), K ID represents I D feedback gain (kV/kA), T DC represents the time delay of converter (sec), K F denotes constant of gain (in kV/unit) and L denotes inductance of coil (H). In the time domain, the change in the SMES unit's real power is represented as: In the time domain, the stored energy at any time in the SMES unit is expressed as: where K ID , K F and T DC are considered as adjustable variables for the definite limits of maximum and minimum.

Designing of the System Controller
For the controllers employed in this study, the most vital objectives are to (i) acquire the highest possible power output from the WTG, (ii) control the frequency and power deviation and (iii) manage power among the consumption and production of IHPS. This section depicts the design characterization and execution procedure of the controllers (i.e., conventional or IO-PID, FO-PID, IO-F-PID and FO-F-PID).

Essentials of FO Calculus
Fractional calculus is the expansion to the nth-order sequential differentiation and integration (differ-integral) of any random function possessing any real value as its order. The operator for differ-integral is symbolized as a D α b , where α (upper limit) R represents the order of the function a and b (lower limit) are widely employed in fractional calculus.
The differ-integral operator is a single expression to symbolize both the fractional derivative (µ) and integral (λ) [35], expressed as: Among the most popular definitions of the differ-integral operator (such as definitions by Riemann and Liouville, Grunwald and Letnikov, and Caputo) [15], Riemann and Liouville's definition for the operator is generally employed [35] and expressed as: where m is the integer quantity of α, m − 1 < α < m, m N, Γ(x) is the noted Euler's gamma function of x and f (τ) stands for the related function.

IO and FO Based PID Controllers
The TF for the simplified PID controller (GC(s)) [16] can be written as: where G C (s), E(s) and U(s) denote the TF, error signal and the output signal for the controller, respectively; K P represents proportional controller gain; K I represents integral controller gain; K D represents derivative controller gain; and µ and λ represent the differintegral order operator's two added variables for designing controllers. In order to obtain the integer-order (IO) controller, µ and λ should be the same or have a value of zero. By considering variables for the order of differ-integral operator as non-integer, the IO controllers are extended to an FO controller. Hence, the conventional structure of the FO controller holds two independent extra knobs for tuning, i.e., µ and λ, and its IO counterpart. Therefore, some feasible FO controllers are termed as PI λ D µ , PI λ , PD µ and PID µ controllers.
The schematic view for FO-PID along with the IO-PID controller on the λ-µ plane and justification for the differ-integral operator's order is presented in Figure 4. The order of the differ-integral operator may vary depending upon its vertical and horizontal axis positioning [36].

IO and FO Fuzzy-PID Controllers
Three error signals for the proportional controller (P), integral controller (I) and derivative controller (D) are required for the implementation of the fuzzy-PID controller. These error signals first require fuzzification and later need to be defuzzied in the output signal. More membership functions (MFs) are required for expressively including the entire response curve. However, this will result in a large rule base with regard to its three numbers of input and one number of output-based fuzzy-PID controller. Hence, it becomes difficult to handle such a large rule base. Consequently, the traditional fuzzy-PID controller combines the fuzzy-PI and fuzzy-PD controllers. In this arrangement, the two input scaling factors (SFs) are K a and K b and the output SFs are the PI and PD controller gains ( Figure 5) [37,38].
In the "IO-F-PID controller", for FLC, at the input, the order of the rate of change in error as well as at the output order of the integral are integer values (µ = 1 and λ = 1 ( Figure 5), respectively). However, this statement does not stand true in the case of the "FO-F-PID controller"; rather, at the input, the order for the rate of change in error as well as at the output order are not integer values but substituted by counterpart fractional numbers (µ, λ). In FLCs, the input-output SFs of a fuzzy controller [37,39] has more effect Energies 2021, 14, 6459 9 of 23 than tuning the MFs shape for the closed-loop performance. Therefore, in the present study, to attain optimal characteristics with the time domain, the rule base and MFs for the fuzzy controller are assumed as fixed for both "IO-F-PID" and "FO-F-PID" controllers [39]. In addition, more emphasis is put on the input-output SFs (i.e., K a , K b , K PI and K PD ) for tuning, and in the F-FO-PID controller, µ and λ are tuned. The considered MFs shape and rule base for fuzzy controllers are portrayed in Figure 6 and Table 1. The fuzzy linguistic variables depicted in Table 1 such as LNV, SNV, ZEV, SPV and LPV denote a large-negative variable, small-negative variable, zero variable, small-positive variable and large-positive variable, respectively.

Performance Index
For enhanced accomplishment of the controller in controlling frequency and power for the Figure 2 system, an optimization method is required. As FLC has non-linear characteristics, there may be no straightforward approach to optimization. The solution lies with the performance indices (such as integral of absolute error (IAE), integral of squared where t represents considered sampling time.
The tuned FLC MFs along with the input-output SFs using a genetic algorithm (GA) to minimize a weighted summation of ISE normalized by maximum error, the maximum percentage of overshoot and the settling time normalized by simulation time. The researchers improved electromagnetism with GA to minimize the ISE while searching for optimal PI λ D µ parameters. The input-output SFs and differ-integrals of the FO-F-PID controller are tuned while minimizing the weighted sum of various error indices and control signals. The parameter-tuning method of the FO-PID controller is based on GA to find the minimum value of the objective function. To reduce unit step response overshoot, ISE is a more ideal performance criterion; thus, ISE performance criteria were used as the objective function in the design of the controller considering each variable's limits. Therefore, the problem of designing the system is framed as a constraint problem for optimization, expressed in (24)- (30): Exposed to the subsequent limits: (a) For classical controllers: , j = 'W' f or WTG and 'D' f or DEG unit (25) (b) For SMES unit: (c) For SFs of FLC: X λmin λ λmax X = 'W' f or WTG and 'D' f or DEG unit (30) where K P represents proportional, K I represents integral and K D represents derivative gains of classical controllers; the tunable SMES parameters are K SMES , K ID and T DC ; K a and K b are the input of FLC and SF, respectively; and the output of SF and FLC are K PI j and K PDj , respectively. From Equations (24)-(30), the minimum and maximum variable values are limited by the min and max range. Here, the QOHS algorithm is utilized for optimization purposes. The output of the optimization program is ISE (i.e., J); at the same time, values are calculated for IAE, ITAE and ITSE performance indices by the optimization program and these obtained values are contended for exhibiting the efficacy of different considered controllers therein [40].

HS Algorithm
In [25], the authors proposed a novel amendment to the metaheuristic algorithm HS which imitates systematic and natural occurrences. The algorithm is motivated by musicians, analogous to the different combinations of harmony for creating a superior state of harmony and the manner of pitching the instruments. In this algorithm, the obtained solution vector of optimization practice is similar to that of search patterns for the local and global solutions. The distinguishing qualities, such as fewer numerical computations, stochastic random search participation and design to improve vector solution by investigating each prevalent vector solution, contributes to increased reliability and obtaining better solutions for the HS algorithm [41].
The method of optimization for the HS algorithm can ensure better enlightenment through three major measures. i.
Initialization: The variables of the algorithm along with the objective function is specified. The initial condition for the harmony memory (HM) may also be defined. ii.
Harmony improvisation: The novel vector solution is generated by adjusting the pitch and manner of randomization. Additionally, they are associated with the HM vector solution. iii. Selection: Until terminating norms are met, the HM continues to find the best solution vector (harmony).
The pseudocode of HS algorithm and the symbols depicted are detailed in [25,28].

Quasi-Oppositional Learning: A Concept
In [42], to catalyze the convergence rate for an intelligent computational method, the impression of opposition-based learning (OBL) is practiced. Many researchers in their work have reported that the chance of attaining the global optimal vector solution with oppositional vector solution rather than any random vector solution has strong prospects [43]. However, later on, it was demonstrated by Rahnamayan et al. [44] that there is a higher chance of reaching the global best solution with the quasi-opposite points as opposed to opposite points. The definition of the opposite and quasi-opposite numbers along with their respective points as considered for OBL practice and the concept of quasi-opposition is stated in [29,45].

Quasi-Oppositional Population Initialization
The QOHS algorithm initializes its population arbitrarily like various other populationbased optimization algorithms. Even deprived of prior understanding of solution vector(s), with quasi-oppositional learning concepts, better fitting solution vector(s) may be achieved. The structure in this context, as proposed in [29], is as follows: (a) initialization of the evenly spread random population; (b) establishment of the population with the quasi-opposite learning concept; (c) evaluation of the objective function for each variable; and (d) best fitting population selection from the early population of a set.

Quasi-Oppositional Generation Jumping
The present vector of solution compels us to develop a new, healthier solution vector with the application of the evolutionary mechanism [29]. With the assistance of randomization and the pitch adjustment process, a new population is generated, also using the quasi-opposite learning concept, and from the combined population, based on the size of the HM, individuals with higher fitness are chosen [29] considering jumping probability (called jumping rate (J r )). For estimation of the generation jumping of quasi-opposite population, each variable's opposite and their respective middle points are computed, as explained in [45]. Therefore, with the progress of the search process, there is a decline in the search space range of the different points from the initial limit [29].

QOHS Algorithm
The pseudo-code of the QOHS algorithm is presented in [29]. Steps 2 and 3 of the algorithm describe HM initialization randomly and with quasi-oppositional concepts, respectively. Quasi-oppositional generation jumping for the QOHS algorithm is presented in Step 6 [29].

Simulation and Result Analysis
In this study, to analyze the frequency and power deviations, the IHPS model TF (Figure 2) involving wind and diesel generators along with ESD (i.e., SMES) is taken into consideration. Controlling the actions of the DEG unit through the governor and the WTG unit with the pitch angle is performed by Controllers 1 and 2, respectively, in the considered IHPS TF model (Figure 2). A comparative study based on simulations of the power and frequency deviations of four controllers for the various IHPS model configurations are presented below: (a) Controllers 1 and 2 are considered as "PID controller" and the studied model is labelled as "PID"; (b) Controllers 1 and 2 are weighted by "FO-PID controller" and the studied model is labelled as "FO-PID"; (c) Controllers 1 and 2 are chosen as "IO-F-PID controller" and the studied model is labelled as "IO-F-PID"; and (d) Controllers 1 and 2 are realized as "FO-F-PID controller" and the model is labelled as "FO-F-PID".
In the current study, due to change in load demand, deviations occur in the frequency and power for the IHPS model considered. The QOHS algorithm is utilized for minimizing the objective function value (i.e., J, as expressed in Equation (24)) of various optimizable variables from the ESD and the undertaken controllers, thus restricting the frequency and power deviations. The considered model is analyzed for its performance and robustness.  Figure 7 for (a) step and (b) random load demand for the present simulation work.
Step load demand is operated at t = 1 s In the corresponding tables, the findings with concern are marked as bold.

Performance Analysis of Different Controller-Based Configurations
In this segment, the performances of different configurations for the controllers ("PID", "FO-PID", "IO-F-PID" and "FO-F-PID") in the time domains of the IHPS model are analyzed. Two different conditions are considered for the analytical performance study, and the assessments of the relative performance under different cases for the considered IHPS model for four types of controllers are presented. The QOHS algorithm is employed for each case and every controller type for the optimal controller parameters.

Case Study I: Performance Analysis for 1% Step Input Load Demand Condition
For performance analysis of the four distinct configurations of the IHPS model ("PID", "FO-PID", "IO-F-PID" and "FO-F-PID"), a 1% step abrupt increment in demand load is operated at t = 1 s (Figure 7a). The related results of frequency and power deviation are exhibited in Figure 8. The optimal gains are represented in Tables 2 and 3 of the controllers and the SMES unit for diverse IHPS model configuration under a 1% rise in the demand load. It can be observed from Figure 8 that the IHPS configurations based on the FO controller (i.e., "FO-PID" and "FO-F-PID") have an advantage over the controller configurations based on IO (i.e., "PID" and "IO-F-PID") in terms of maximum overshoot (Mp) (i.e., less Mp) and less settling time (ts). Through the inclusive assessment, the "FO-F-PID" configuration is confirmed as the most effective in terms of t s and M P and is established as superior to the other three configurations. Additionally, the report regarding the frequency deviation under the transient state (from the plots obtained for different configurations) is presented in Table 4, representing the individual controller-based IHPS models. Table 4 confirms the "FO-F-PID" configuration's improved results in terms of controlling the variation in power and frequency over the other considered configurations. In addition, in Table 4, the objective function value (J = ISE) and the generated values for the remaining three performance indices are depicted. It can be observed that the "FO-F-PID" configuration of the IHPS has the lowest values for the performance indices as well as for the objective function. This endorses the character of responses depicted in Figure 8 (i.e., nature of response for various undertaken IHPS model configurations). Table 2. "IO-F-PID" and "FO-F-PID" optimizable variables of IHPS for DEG and WTG controllers along with the SMES block.

Gains of Different Controllers
Step Random Step  Table 3. "IO-F-PID" and "FO-F-PID" optimizable variables of IHPS for DEG and WTG controllers along with the SMES block.

Gains of Different Controllers
Step Random Step   To verify the characteristics of various configurations of the undertaken IHPS model under the practical ecosystem, the different configurations were exposed to a randomly

Case Study II: Performance Analysis under Random Load Demand Condition
To verify the characteristics of various configurations of the undertaken IHPS model under the practical ecosystem, the different configurations were exposed to a randomly varying load demand (Figure 7b). Here, to obtain the optimized parameters against diverse controller configurations and SMES parameters, the QOHS algorithm is utilized. Figure 9a,b depict the results of the acquired deviation for the frequency and power of each IHPS model configuration for the variation of random load. At the same time, control over the system and dominating the power and frequency deviation at the earliest after every variation in load can be observed (Figure 9), in which the IHPS model with "FO-F-PID" configuration has an edge over the other considered configurations (i.e., "PID", "FO-PID" and "IO-F-PID").
The controller gain parameters optimized with the QOHS algorithm for the considered configurations and tunable SMES parameters with the studied case are shown in Tables 2 and 3. In addition, a zoomed-in view of the deviation in power and frequency is displayed in Figure 9, depicting, with additional transparency, the nature of every IHPS configuration in terms of ts, steady-state error (Ess) and Mp (i.e., less Mp).

Robustness Analysis of Different Controller-Based Configurations
In this study, to validate the robustness of the QOHS algorithm and FO controller, two categories of investigation were carried out with the undertaken configurations of the IHPS model. In the first investigation, significant nonlinearity was introduced [46] under the practice of incorporating rate constraint nature in the DEG unit or ESD of the model. During the second investigation, a mismatch for the demand load was instigated with the optimally adjusted variables, designed to study the robustness of the considered IHPS. The first method (i.e., non-linearity by rate constraint) limits specific source energy from accumulating or discharging energy immediately by establishing an upper and lower boundary, which can be regarded as strongly resembling a more practical situation. The application of the non-linearity of rate constraint with the first-order TF is shown in Figure 10.

Case Study III: Rate Constraint Type Non-Linear ESD Operation
For Case Study III, all the configurations undertaken for the study of the IHPS model have been taken into consideration for a robustness study. In the same vein, the nonlinearity block is formulated under the SMES block in Figure 2. In the SMES block, a non-linearity component is inserted considering rate constraint ( Figure 10) with a value of |P SMES | < 0.005. The controllers gain variables and the tunable SMES parameters, which are optimized by the QOHS algorithm for the linear operation.
The performance of the IHPS model's every configuration and the SMES configuration in conjunction with non-linearity in the form of rate constraint under an abrupt increase in load demand in step manner is presented in Figure 11 along with the change in the frequency and power. It can be observed that FO controllers have an edge over the IO controllers. The IHPS model with the FO-F-PID controller proved its competency over the other controllers (PID, FO-PID and IO-F-PID) in controlling the abrupt shift in load ( Figure 11). The transient performance analyses through simulation ( Figure 11) are presented in Table 5.

Case Study IV: Rate Constraint Type Non-Linear DEG Operation
In Case Study IV, the DEG unit non-linearity in the form of rate constraint ( Figure 10) is presented, as depicted in Figure 2. Its variation in power and frequency for the various considered configurations of the IHPS model ("PID", "FO-PID", "IO-F-PID" and "FO-F-PID") are presented in Figure 12 for a load variation with an abrupt increase (Figure 7a). The value employed regarding the rate constraint type is |P DEG | < 0.001. Considering the various sub-plots associated with Figure 12, the situation revealed that an improved controlling nature may be obtained with FO controllers over IO controllers under the non-linearity effect consideration. The HIPS model with the "FO-F-PID" configuration has an edge over the remaining three configurations ("PID", "FO-PID" and "IO-F-PID") during unexpected deviations in load demand with non-linearity at the DEG unit in the form of rate constraint. Figure 12 shows the nature of transient performance and Table 6 depicts the performance indices conceded.  For the analysis of the robustness of the system with loading mismatch, optimization is achieved for the controllers' gains (positioned on the side of the DEG and the WTG) along with the SMES optimizable variables with the QOHS algorithm for 90% of the loading ability. In this case, if there is a loading mismatch, the system is loaded only to 85% loading ability. Therefore, robustness evaluation for the present case is carried out in consideration of the two scenarios, as depicted in Table 7.
The optimized gains for Controller 1 and Controller 2 are offered at Table 3 along with the SMES block variables for the IHPS model under the "FO-F-PID" configuration as achieved with the QOHS algorithm. The comparative profiles for the deviation in the frequency and power of the analyzed FO-F-PID controller in accordance with the IHPS model under Scenario I and Scenario II are presented in Figure 13a,b. Figure 13 demonstrates that the deviation profile for the frequency and the power in the two considered scenarios are overlapping, which offers robust optimization of the FO-F-PID controller through the QOHS algorithm. Information related to transience and the distinguished impact of the current work for both the considered scenarios is exhibited in Table 8.  Case Study VI is characterized by 95% system loading capacity, and the associated FO-F-PID controller gains for DEG, WTG along with SMES tunable variables are optimized for the capacity of 90% loading. Hence, in Table 9, the robustness assessment on this case is shown by considering both situations.
The QOHS algorithm as well as the FO-F-PID controller's robustness can be easily verified with assistance from Figure 14a,b under both considered scenarios for the deviation in frequency and strength with the FO-F-PID controller-based IHPS model. Table 3 depicts the FO-F-PID controller's optimal values and the adjustable variables from SMES for a sudden 10% load demand increase in step manner (same as Case V). Figure 14 depicts the two different scenarios as having identical responses, confirming the QOHS algorithm's robustness for optimizing the adjustable variables. The dissimilarity between the situations is represented in Table 8 in the manner of transient data for the generated characteristics of transient nature.

Conclusions
In the present study, challenges such as (a) the overemployment of DEG systems for electrification, (b) fossil fuel rapid exhaustion, (c) countryside electrification and (d) global warming are foreseen and solved with a trustworthy and controlled inexhaustible energy source based on IHPS. The following outcomes are achieved with various system configurations: (a) For the studied IHPS, the hybridized classical fuzzy controller using the concept of FO of calculus is projected and the discrete system configurations framed are labelled as "PID", "FO-PID", "IO-F-PID" and "FO-F-PID". (b) The FO concept applied in the controllers upholds its capability in governing the dynamics of the system. (c) For the IHPS application, the priority is to set up a static configuration for MF and the rule base of FLCs. (d) In the context of the considered associated power engineering problems, the utilization of the QOHS scheme proves to be effective in generating close proximity to the global optimal solution. (e) For the studied IHPS model, the "FO-PID" configuration's performance and robustness are no less than that of the "IO-F-PID" arrangement. (f) In governing the deviations for both power and frequency in rate constraint-centered robustness analysis, the FO controllers proved their effectiveness and were found to be more competitive compared to the IO controller. (g) The complete performance and robustness analysis conclude that under the changing load conditions, the "FO-F-PID" controller is the most effective in limiting the changes for both power and frequency, followed by "IO-F-PID", "FO-PID" and "PID" controllers.