Next Article in Journal
Investigating the Effects of a Fractional Operator on the Evolution of the ENSO Model: Bifurcations, Stability and Numerical Analysis
Next Article in Special Issue
First-of-Its-Kind Frequency Enhancement Methodology Based on an Optimized Combination of FLC and TFOIDFF Controllers Evaluated on EVs, SMES, and UPFC-Integrated Smart Grid
Previous Article in Journal
Numerical Analysis of Direct and Inverse Problems for a Fractional Parabolic Integro-Differential Equation
Previous Article in Special Issue
Fractional-Order Model-Free Predictive Control for Voltage Source Inverters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Design Optimization of Improved Fractional-Order Cascaded Frequency Controllers for Electric Vehicles and Electrical Power Grids Utilizing Renewable Energy Sources

by
Fayez F. M. El-Sousy
1,*,
Mohammed H. Alqahtani
1,
Ali S. Aljumah
1,
Mokhtar Aly
2,
Sulaiman Z. Almutairi
1 and
Emad A. Mohamed
1,3
1
Department of Electrical Engineering, College of Engineering, Prince Sattam bin Abdulaziz University, Al Kharj 16278, Saudi Arabia
2
Facultad de Ingeniería, Arquitectura y Diseño, Universidad San Sebastián, Bellavista 7, Santiago 8420000, Chile
3
Department of Electrical Engineering, Faculty of Engineering, Aswan University, Aswan 81542, Egypt
*
Author to whom correspondence should be addressed.
Fractal Fract. 2023, 7(8), 603; https://doi.org/10.3390/fractalfract7080603
Submission received: 26 May 2023 / Revised: 17 July 2023 / Accepted: 18 July 2023 / Published: 4 August 2023

Abstract

:
Recent developments in electrical power grids have witnessed high utilization levels of renewable energy sources (RESs) and increased trends that benefit the batteries of electric vehicles (EVs). However, modern electrical power grids cause increased concerns due to their continuously reduced inertia resulting from RES characteristics. Therefore, this paper proposes an improved fractional-order frequency controller with a design optimization methodology. The proposed controller is represented by two cascaded control loops using the one-plus-proportional derivative (1 + PD) in the outer loop and a fractional-order proportional integral derivative (FOPID) in the inner loop, which form the proposed improved 1 + PD/FOPID. The main superior performance characteristics of the proposed 1 + PD/FOPID fractional-order frequency controller over existing methods include a faster response time with minimized overshoot/undershoot peaks, an ability for mitigating both high- and low-frequency disturbances, and coordination of EV participation in regulating electrical power grid frequency. Moreover, simultaneous determination of the proposed fractional-order frequency controller parameters is proposed using the recent manta ray foraging optimization (MRFO) algorithm. Performance comparisons of the proposed 1 + PD/FOPID fractional-order frequency controller with existing PID, FOPID, and PD/FOPID controllers are presented in the paper. The results show an improved response, and the disturbance mitigation is also obtained using the proposed MRFO-based 1 + PD/FOPID control and design optimization methodology.

1. Introduction

Climate change and its serious effects on different environmental conditions has motivated the urgent transition to new renewable and clean sources in energy-related sectors [1,2]. In the energy generation sector, renewable energy sources (RESs) have dominated the recent installation shares. Specifically, photovoltaic (PV) and wind energy sources are the most commonly used in modern electrical power grids [3]. In addition, the increased use of electric vehicles (EVs) in the transportation sector helps to combat climate change and utilize clean energy sources. In [4], a study on the benefits and impacts of adding more functionalities of EVs into electrical grids with RESs was presented. The study proved the ability of EVs to actively share the required tasks of RES inverters, and hence there was no need for high-capacity energy storage devices. Therefore, cleaner and sustainable energy systems characterize modern electrical power grid systems [5].
Connecting RESs with electrical power grids requires power electronic-based power converters for grid integration of different RESs and for maximizing the extracted power from RESs. This makes RESs different from traditional synchronous-generator-based power generation; therefore, RESs impact electrical power grids with their lowered inertia [6]. The reduced inertial characteristics of modern electrical power grids with high participation rates of RESs make them highly susceptible to voltage and frequency stability problems. Therefore, proper frequency regulation is mandatory to mitigate the reduced inertia of RESs in modern electrical power grids. The load frequency control (LFC) has been shown to be a more suitable method for regulating frequency in modern electrical power grids with RESs [7].
The literature includes several research works on the topic of developing proper frequency regulation controllers [8]. Data-driven-based neural network (NN) controllers were recently presented for several electrical power grid case studies. However, they require huge computational burdens for data processing and designing NNs. From another side, model-based predictive controllers (MPCs) were introduced in the literature with different MPC structures. However, they are sensitive to the modelling of the process and system parameters. Additional LFC proposals have been proposed in the literature, such as the sliding mode LFC, machine and deep learning-based LFC, linear matrix inequalities (LQR), etc. [9]. Additionally, type-2 fuzzy modelling and control methods have been shown to have better performance in the literature. For instance, the type-2 fuzzy modelling has provided improved modelling for non-linear systems in [10,11]. Interval type-2 fuzzy modelling with fractional-order (FO)-based LFC methods has been shown in [12,13]. These methods merge the benefits of fuzzy type-2 with FO control methods.
From another perspective, integer-order (IO) frequency regulators and fractional-order (FO) frequency regulators have been widely discussed in the literature [14]. Different combinations have been presented using basic terms such as proportional (P), derivative (D), integral (I), tilt (T), filtering (F), fractional filtering (FF), and FO operators. Furthermore, various cascaded solutions of IO and FO have been developed in the literature to provide a better disturbance rejection capability. Metaheuristic-based optimized determination methods for control parameters have been developed in the literature to ensure the best control parameters without the need for complex control modelling and design tools [15].
An I-based controller was designed for LFC in interconnected electrical power grids in the literature [16]. Another adaptive I controller with an optimized design using Jaya balloon optimizers (JBO) was discussed in [17]. A PI controller was introduced in [18] with a binary-moth flame optimizer (MFO), whereas the hybrid gravitation searching with the firefly optimizer algorithm (hGFA) was introduced in [19]. Moreover, a PI controller with Harris–Hawks optimization (HHO) was provided in [20], and the grey wolf optimizers (GWO) were provided in [21]. PI-based LFC methods have succeeded in improving system dynamics; however, they fail to deal with system non-linearities and parameter uncertainties. Additionally, a PID controller with the imperialist competitive optimization algorithm (ICA) was proposed in [22]. Moreover, a PID virtual inertia controller was proposed in [23] in order to improve the inertial response for a real electrical power grid case study. IO-based frequency controllers showed simple design requirements, easy implementation, and lower costs for implementation. However, IO-based frequency regulators cannot fully mitigate uncertainty and fluctuations resulting from the electrical power grid parameters. Moreover, they showed a poor response to any uncertainties in the power system’s parameters and at low-inertia operation.
Different control structures have been provided in the literature using cascaded control loops [24]. In cascaded LFC methods, two loops are used to construct the frequency regulation controller. The area control error (ACE) is used as an input for the outer loop, and the frequency deviation signal is employed for the inner loop [25]. This leads to a higher degree of freedom and better rejection of existing disturbances. In some cascaded controllers, the tie-line power is also included in the inner loop. The ACE signal represents the slow dynamics loop, and the frequency deviation represents a faster loop of system frequency disturbances.
The cascaded IO-based frequency regulator using the PD-PI scheme was introduced in [26], applying the enhanced slime mould optimization algorithm (ESMOA) for design optimization. The PI-PDF controller was provided in [27] with the driver training-based optimizer (DTBO). The results showed improved disturbance mitigation of the power grids. Another fuzzy logic control (FLC) cascaded PI-PDF controller with scaled factors and the modified Dragonfly optimizer algorithm was proposed in [28]. Additional IO-based frequency regulators have been proposed in the literature, such as PIDF [29], 2DoF-PID [24], PD-PID [30], PI-(1 + DD) [25], PID2D [31], IPD-(1 + I) [32], FLC-PID [33], and the neuro-fuzzy LFC [34]. The associated optimizers include the slap-swarm-based algorithm (SSA) [35], and the flower pollination-based algorithm (FPA) [36]. Some local modelling LFC was provided in [37,38]. The cascaded IO-based frequency regulators have shown better mitigation for existing electrical power grid disturbances. However, the use of IO control methods has a lower number of tunable parameters compared with FO-based control methods.
In the literature, FO-based frequency controllers have also found widespread use in regulating frequency in electrical power grids with different structures. Some metaheuristic optimizers have been presented for the design optimization of FO-based frequency regulators, such as the sine cosine-based algorithm (SCA) [39], genetics algorithm (GA) [22], and movable-damped-wave-based algorithm (MDWA) [40]. In [41], a review of the possible cascaded and multiple input-based LFC methods has been presented. Several two- and three-input schemes exist in the literature. A cascaded FOPID FO-based frequency regulator was proposed in [42] for stand-alone electrical power grids, whereas the FOPIDF was provided in [43]. A higher degree-of-freedom cascaded 3DOF TID-FOPID was provided in [44] to enhance electrical power grid stability. Other cascaded FO-based frequency controllers have been provided based on the FOID-FOPIDF in [45], FO-IDF in [46], and PI-TDF in [47]. Examples of associated optimizers include the pathfinder-based algorithm (PFA) in [48], the artificial-bee-colony-based algorithm (ABC) in [49], differential evolution-based algorithm (DE) in [50], and the SSA optimization in [47]. A cascaded FOPD-PI controller, considering plug-in EVs (PEVs), was presented in [51]. Another ID-T cascaded controller was proposed in [52] with an Archimedes optimizer algorithm (AOA). The inclusion of FO operators in cascaded LFC methods increase the flexibility and the number of tunable parameters. This can be reflected as a better optimization of the system’s response to disturbances. A common difficulty in applying FO control systems is their implementation complexity. For instance, using Oustaloup’s recursive approximation (ORA) representation of the fifth-order leads to eleven-order IO equivalent representations. However, this difficulty can be solved with the recent powerful microcontroller systems.

Paper Contribution

It has become obvious that several types and structures exist of IO- and FO-based frequency regulators for LFC in interconnected electrical power grids. Moreover, several techniques of metaheuristic optimization algorithms have been associated with the presented controllers for design optimization and reaching the best parameters. Proper and optimum design and selection of LFC and optimization method are crucial when facing expected reduced inertia with high participation levels of RESs. Additionally, probable local minimum settling represents another issue for several metaheuristic optimizers. Therefore, this paper presents an improved 1 + PD/FOPID FO-based frequency regulator, using the recent manta ray foraging optimization (MRFO) for design optimization.
Based on authors’ knowledge, this is the first time the (1 + PD) has been combined with the FOPID controller in a cascaded way to provide a hybrid high-flexibility frequency regulation controller. Additionally, the integration of the MRFO algorithm with the proposed 1 + PD/FOPID controller leads to providing joint optimum behaviour and a better parameter determination process. This confirms the aforementioned findings that the power grid frequency regulation performance is determined jointly by the used control methodology and the applied optimization algorithm. Table 1 provides a summary of the existing controllers and the proposed controller. The main contributions in this paper are as follows:
1.
An improved controller and design optimization method is proposed for frequency regulation in interconnected electrical power grids with high participation levels of RESs in addition to active participation of EVs in regulating frequency. The proposed controller and design methodology can effectively lead to mitigating various existing frequency fluctuations in electrical power grids. The proposed method can be generalized and applied to various electrical power grid systems and components.
2.
The proposed frequency regulation control methodology is formed using a cascaded 2DoF 1 + PD/FOPID control method, which utilizes two input signals (namely the frequency deviation in each area, and control error in each area (ACE)). The utilization of two different signals is beneficial for the mitigation process of low- and high-frequency existing disturbances.
3.
The proposed frequency regulation methodology using a 2DoF 1 + PD/FOPID controller provides better frequency regulation responses compared with the widely utilized PID, FOPID, and PD/FOPID LFCs, providing better disturbances rejections capabilities. The proposed 2DoF 1 + PD/FOPID structure is capable of mitigating various deviations in area frequency and electrical power grid tie-line power as a direct result of employing two cascaded loops with frequency and ACE signals.
4.
Benefiting from EVs’ batteries in the effective participation in frequency regulation is coordinated through the proposed 2DoF 1 + PD/FOPID structure. Therefore, the proposed 2DoF 1 + PD/FOPID structure reduces the frequency regulation complexities due to employing the centralized frequency regulation structure that coordinates the connected EVs’ batteries and LFC regulator.
5.
An improved design optimization methodology using the recent manta ray foraging optimization (MRFO) to determine the best parameters for the proposed 2DoF 1 + PD/FOPID frequency regulation. The optimized values of LFCs in different electrical power grids are simultaneously searched using the MRFO optimizer, thus minimizing the desired objectives.
The remainder of the paper is divided as follows: Section 2 provides the mathematical model representations for the interconnected electrical power grids. The proposed 2Dof 1 + PD/FOPID frequency regulation is detailed in Section 3. The proposed design optimization of the 1 + PD/FOPID controller is described in Section 4. Section 5 provides the obtained simulation results of the interconnected two-area electrical power grids with EVs participations and RESs. Finally, the paper’s conclusions are presented in Section 6.
Table 1. Summarized comparison of existing LFC methods and the paper’s contribution.
Table 1. Summarized comparison of existing LFC methods and the paper’s contribution.
Ref.CategoryControl SchemesCharacteristics
 [16,17,18,19,20,22]Conventional IO LFC (single input)I, PI, PID, PIDF
  • Simple structure and implementation;
  • Low ability to mitigate disturbances;
  • Lower robustness against parameters uncertainty.
 [22,39,40,42,43]Conventional FO LFC (single input)FOPI, FOPID, FOPIDF
  • Increased flexibility and number of parameters compared to IO LFC methods;
  • Limited rejection capability of existing disturbances;
  • Lower mitigation of high-frequency deviations.
 [24,25,27,28,29,30,32,33]Cascaded IO LFC (multiple inputs)PD-PI, PI-PDF, 2DoF-PID, PD-PID, PI-(1 + DD), IPD-(1 + I), FLC-PID
  • Easy to be implemented;
  • Improved ability to mitigate disturbances;
  • limited number of tunable parameters (only gains can be tuned).
 [44,45,46,47,48,49,51,52]Cascaded FO LFC (multiple inputs)Cascaded FOPID, 3DOF TID-FOPID, FOID-FOPIDF, FO-IDF, PI-TDF, ID-T, FOPD-PI
  • Higher number of tunable parameters compared to IO LFC methods;
  • Increased design flexibility compared to other IO-based LFC methods;
  • Enhanced disturbance rejection capability compared to single-input FO LFC methods;
  • Improved mitigation ability of high-frequency deviations.
ProposedProposed cascaded FO LFC (multiple inputs)Proposed cascaded 2DoF 1 + PD/FOPID LFC method
  • New hybrid (1 + PD) cascaded with FOPID controller for LFC application;
  • Merging the characteristics of IO-based control with the added flexibility of FO control;
  • Simultaneous design optimization process of all tunable parameters, considering connected devices using a powerful MRFO algorithm;
  • Active contribution of EVs in regulating frequency, and this functionality is considered during the design process;
  • High disturbance rejection capability;
  • Better mitigation of high-frequency fluctuations due to using two loops and employing the frequency deviation signal.

2. Modelling of Interconnected Electrical Power Grids

2.1. Electrical Power Grid Description

The interconnected electrical power grid using two areas is selected to verify the proposed 1 + PD/FOPID and design optimization method. Figure 1 shows two electrical power grid areas (area a and area b) connected using an AC bus as tie-line between the areas. The electrical power grid area a includes a thermal power plant, local loads, EV batteries, and a wind RES, whereas electrical power grid area b includes a hydraulic power plant, local loads, EV batteries, and a PV RES. It is assumed in the analysis that EVs are equally shared by the two electrical power grids. Each electrical power grid has a frequency regulation controller that controls the power injection of each power/storage devices in the area. The parameters for the studied electrical power grids in this work are taken from [14] and listed in Table 2. Figure 2 presents the complete implemented model of the studied electrical power grid elements with the EV batteries and RESs.
The controlled system consists of the aforementioned two-area interconnected power grids with the connected elements in each area. The input error to be controlled is the ACE signal of each area, and it has to be maintained at a zero reference value. The output of the control method adjusts the contribution of the connected generation and/or energy storage devices to mitigate the frequency changes. When the system has unbalanced generation and loading, this is reflected as an increase/decrease in the system frequency. Therefore, the LFC method mitigates the frequency fluctuations and preserves the frequency deviations at a zero value.

2.2. RES Behaviour Models

In the studied electrical power grids, participation of PV and wind is considered with their intermittency proprieties. They rely on environmental conditions, such as wind speed, solar irradiance levels, and temperatures. For these reasons, continuous tracking is mandatory for the maximum power operating point (MPPT controller), in which, power electronics conversion blocks play a vital role. In addition, they are responsible for grid integration and synchronization of RESs. The outputted powers from RESs are continuously fluctuating due to variations in the weather conditions.
For PV generation, PV outputted power is unpredictable due to its intermittency. This, in turn, results in high-frequency power fluctuations that can lead to severe stability problems in the electrical power grids. In this work, the PV model from [53] is employed for PV power as follows:
P P V = η Φ s o l a r A P V [ 1 0.005 ( T a + 25 ) ]
where  η  stands for PV panel’s conversion efficiency (in %),  Φ s o l a r  stands for solar insolation (W/m 2 ),  A P V  stands for the area of the PV unit (m 2 ), and  T a  stands for the ambient temperature (°C). The implemented PV configuration model represents a realistic PV power generation model such as in the presented model in [54].
From another side, the outputted mechanical wind power from the turbine possesses high fluctuations resulting from its intermittency characteristics. The wind speed continuously varies, and hence the outputted power differs from instant to another. The calculation of the mechanical power is as follows [55]:
P w i n d = 1 2 ρ A r C p V w 3
where  ρ  stands for the density of air (in kg/m 3 ),  A r  stands for the swept area (in m 2 ),  C p  stands for the power coefficient, and  V w  stands for the speed of wind (in m/s). The implemented configuration of the wind power model is based on the realistic representation using the presented modelling in [54].

2.3. EV Behaviour Model

In this study, participation of the installed PV batteries is also considered for the frequency regulation control of the electrical power grids. This, in turn, leads to eliminating the need for additional energy storage capacities in the electrical power grids. The connected EV batteries are charged/discharges based on control signals coming from the electrical power grids. Accordingly, a better performance of the electrical power grid’s reliability, stability, efficiency, and transient response is achieved. In frequency regulation in the electrical power grids, the EV batteries regulate frequencies against the fluctuations resulting from the RESs and connected load proprieties. The used EV battery model is included in Figure 2 for frequency regulation studies in electrical power grids as in [56], in which, Nernst’s equation-based EV battery model determines the dependency of the open circuit battery voltage ( V o c ) on their state of charge ( S O C ) as follows:
V o c ( S O C ) = V n o m + S p R T F l n ( S O C C n o m S O C )
where  V o c ( S O C )  stands for the  V o c  dependency on  S O C V n o m  stands for the nominal battery voltage, and  C n o m  stands for the nominal battery capacities of the EVs (in Ah). Furthermore,  S p  denotes the sensitivity parameter of the  V o c  and  S O C  of the batteries. R and F are the gas and Faraday constants, respectively, and T is the temperature.

2.4. System State Space Model

The studied system in Figure 2 is represented mathematically in Appendix A. The state space-based linear representation is employed for the proposed analysis. The models in Appendix A are collected in the state space model. The general representation of the state space model is modelled as:
x ˙ = A x + B 1 ω + B 2 u
y = C x
where x stands for a vector including state variables, y stands for a vector including output states,  ω  stands for a vector including the disturbance of a system, and u stands for a vector including the control output. The generation system is modelled using Laplace transform, and this model is employed to define the system state variables vector x. Whereas the load and RESs are considered as disturbances in this model for defining the vector  ω . In (4), vectors x and  ω  are expressed as:
x = Δ f a Δ P g a Δ P g a 1 Δ P W T Δ f b Δ P g b Δ P g b 1 Δ P g b 2 Δ P P V Δ P t i e T
ω = Δ P l a P W T Δ P l b P P V T
where  Δ P g a  and  Δ P g a 1  are the governor and turbine outputs of the thermal unit in area a, respectively. Whereas  Δ P g b Δ P g b 1 , and  Δ P g b 2  are the governor, droop compensation, and penstock turbines outputs of the hydraulic generation in area b, respectively. In (4), vector u includes frequency regulation signals of each electrical power grid  A C E o a , and  A C E o b  in addition to the participated power by EVs ( Δ P E V a  and  Δ P E V b ) as follows:
u = A C E o a Δ P E V a A C E o b Δ P E V b T
The representative matrices in state space modelling in (4), (A B 1 B 2 , and C) are obtained from the electrical power grid model in Figure 2 and its parameters. They are represented as:
A = D a 2 H a 1 2 H a 0 1 2 H a 0 0 0 0 0 1 2 H a 0 1 T t 1 T t 0 0 0 0 0 0 0 1 R a T g 0 1 T g 0 0 0 0 0 0 0 0 0 0 1 T WT 0 0 0 0 0 0 0 0 0 0 D b 2 H b 1 2 H b 0 0 1 2 H b 1 2 H b 0 0 0 0 2 T R R b T 1 T 2 2 T W 2 T 2 + 2 T W T 2 T W 2 T R 2 T 1 T 1 T 2 0 0 0 0 0 0 T R R b T 1 T 2 0 1 T 2 T 1 T R T 1 T 2 0 0 0 0 0 0 1 R b T 1 0 0 1 T 1 0 0 0 0 0 0 0 0 0 0 1 T PV 0 2 π T tie , eq 0 0 0 2 π T tie , eq 0 0 0 0 0
B 1 = 1 2 H a 0 0 0 0 0 0 0 0 0 0 0 0 K WT T WT 0 0 0 0 1 2 H b 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 K PV T PV 0 0 0 0 ,   and B 2 = 0 1 2 H a 0 0 0 0 0 0 1 T g 0 0 0 0 0 0 0 0 0 0 1 2 H b 0 0 2 T R T 1 T 2 0 0 0 T R T 1 T 2 0 0 0 1 T 1 0 0 0 0 0 0 0 0 0
C = 1 0 0 0 0 0 0 0 0 0 B a 0 0 0 0 0 0 0 0 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 B b 0 0 0 0 1
Table 2. The electrical power grid parameters for the modelled case study (with  x { a , b } ), [14].
Table 2. The electrical power grid parameters for the modelled case study (with  x { a , b } ), [14].
SymbolsValue
Area aArea b
P r x  (MW)12001200
R x  (Hz/MW)2.42.4
B x  (MW/Hz)0.42490.4249
Valve min. limit  V v l x  (p.u.MW)−0.5−0.5
Valve max. limit  V v u x  (p.u.MW)0.50.5
T g  (s)0.08-
T t  (s)0.3-
T 1  (s)-41.6
T 2  (s)-0.513
T R  (s)-5
T w  (s)-1
H x  (p.u.s)0.08330.0833
D x  (p.u./Hz)0.008330.00833
T P V  (s)-1.3
K P V  (s)-1
T W T  (s)1.5-
K W T  (s)1-
EVs Modelling
Penetration level5–10%5–10%
V n o m  (V)364.8364.8
C n o m  (Ah)66.266.2
R s  (ohms)0.0740.074
R t  (ohms)0.0470.047
C t  (farad)703.6703.6
  R T / F 0.026120.02612
Minimum EVs SOC %1010
Maximum EVs SOC %9595
Minimum capacity EVs limit (p.u.MW)−0.1−0.1
Maximum capacity EVs limit (p.u.MW)+0.1+0.1
C b a t t (kWh)24.1524.15
Figure 2. Modelling of the various components of the studied electrical power grids.
Figure 2. Modelling of the various components of the studied electrical power grids.
Fractalfract 07 00603 g002

3. The Proposed 2Dof 1 + PD/FOPID Frequency Regulation

3.1. FO-Based Frequency Regulator Representation

The FO-based frequency regulator systems are based on FO calculus and representations of non-integer systems. There are several schemes provided in the literature, such as Grunwald–Letnikov, Riemann–Liouville, and Caputo [57]. The  α t h  FO-based derivative for function f between a and t is defined using the Grunwald–Letnikov method:
D α | a t = lim h 0 1 h α r = 0 t a h ( 1 ) r n r f ( t r h )
where h is the employed step time, and  [ · ]  represents the used operator based on integer terms for the Grunwald–Letnikov method. Whereas n should satisfy the condition ( n 1   < α < n ). The binomial-based coefficients are represented using:
n r = Γ ( n + 1 ) Γ ( r + 1 ) Γ ( n r + 1 )
where the used gamma function within (13) is usually defined using:
Γ ( n + 1 ) = 0 t x 1 e t d t
From another side, the Riemann–Liouville method defines FO-based derivatives by avoiding the utilization of the sum and limits, while it uses the integer derivative in addition to the integral representations as follows:
D α | a t = 1 Γ ( n α ) ( d d t ) n a t f ( τ ) ( t τ ) α n + 1 d τ
The Caputo representation for the FO derivative is as follows:
D α | a t = 1 Γ ( n α ) a t f ( n ) ( τ ) ( t τ ) α n + 1 d τ
Generally, the FO non-integer operator  D α | a t  can be denoted as follows:
D α | a t = α > 0 d α d t α FO derivative α < 0 t 0 t f d t α FO integral α = 0 1
Implementing FO-based frequency regulators using Oustaloup’s recursive approximation (ORA) has been shown as a suitable digital representation scheme in real-time implementations [57]. The ORA is selected in this paper for FO-based frequency regulator implementations. In which the  α t h  derivatives ( s α ) are represented as [57]:
s α ω h α k = N N s + ω k z s + ω k p
where  ω k p  and  ω k z  stand for the pole/zero locations within the  ω h  sequence, and their calculations are as follows:
ω k z = ω b ( ω h ω b ) k + N + 1 α 2 2 N + 1
ω k p = ω b ( ω h ω b ) k + N + 1 + α 2 2 N + 1
ω h α = ( ω h ω b ) α 2 k = N N ω k p ω k z
where the approximated FO-based frequency regulator operator function possesses a  ( 2 N + 1 )  pole/zero number. Therefore, N is related to ORA’s order of representation and equals  ( 2 N + 1 ) . In this work, ORA is employed for FO-based frequency regulator representation with  ( M = 5 )  within the frequency range ( ω [ ω b , ω h ] ) between  [ 0.001 , 1000 ]  rad/s.

3.2. Controllers from the Literature

The literature has several proposals of frequency regulation control schemes, including IO- and FO-based frequency controllers. Some existing IO-based frequency regulators include the I, PI, PID, and PIDF LFCs, and their transfer functions (TFs)  C ( s )  are as follows:
C I ( s ) = Y ( s ) E ( s ) = K i s C P I ( s ) = Y ( s ) E ( s ) = K p + K i s C P I D ( s ) = Y ( s ) E ( s ) = K p + K i s + K d s C P I D F ( s ) = Y ( s ) E ( s ) = K p + K i s + K d s N f s + N f
Whereas FO-based frequency regulators in the literature include the following:
C F O I ( s ) = Y ( s ) E ( s ) = K i s λ C F O P I ( s ) = Y ( s ) E ( s ) = K p + K i s λ C F O P I D ( s ) = Y ( s ) E ( s ) = K p + K i s λ + K d s μ C F O P I D F ( s ) = Y ( s ) E ( s ) = K p + K i s λ + K d s μ N f s + N f C T I D ( s ) = Y ( s ) E ( s ) = K t s ( 1 n ) + K i s + K d s C F O T I D ( s ) = Y ( s ) E ( s ) = K t s ( 1 n ) + K i s λ + K d s μ C P F O T I D ( s ) = Y ( s ) E ( s ) = K p + K t s ( 1 n ) + K i s λ + K d s μ
It has become evident that each frequency regulator includes a number of tunable parameters to optimize its performance. The number of tunable parameters relies on the employed control scheme.

3.3. Proposed 1 + PD/FOPID Controllers

Figure 3 shows the structure of the proposed 1 + PD/FOPID controller for interconnected electrical power grids. The proposed 1 + PD/FOPID controller uses two cascaded loops using the ACE signal of each electrical power grid in the area in the outer loop and the frequency deviation signal of the area in the inner loop. The 1 + PD is employed for the outer loop and FOPID is employed for the inner loop. Therefore, the proposed 1 + PD/FOPID controller can benefit from both the characteristics of 1 + PDF in the outer loop and the FOPID FO controller in the inner loop. Moreover, the proposed 1 + PD/FOPID controller uses two different inputs signals with different characteristics. That is, the ACE can mitigate low-frequency-related disturbances, while frequency deviation can mitigate high-frequency-related disturbances. Therefore, the proposed 1 + PD/FOPID achieves a fast and robust response, low values of overshoot/undershoot peaks, and a high rate for rejecting disturbances with various frequency scales. The modelling of inputs into the ACE loop ( ( A C E a )  and  ( A C E b ) ) and into the controller inputs ( E a 1 ( s )  and  E b 1 ( s ) ) is provided as:
E a 1 ( s ) = A C E a = Δ P t i e + B a Δ f a E b 1 ( s ) = A C E b = A a b Δ P t i e + B b Δ f b
where  ( A a b )  represents the ratio among the capacities of the electrical power grids’ areas a and b, whereas the outputted signals from  Y a 1 ( s )  and  Y b 1 ( s )  of the 1 + PD controller’s loop are expressed as:
Y a 1 ( s ) = [ 1 + K p 1 + K d 1 s ] . E a 1 ( s ) Y b 1 ( s ) = [ 1 + K p 2 + K d 2 s ] . E b 1 ( s )
From (25), each electrical power grid has two parameters for tuning in the outer loop. The electrical power grid in area a has two tunable parameters ( K p 1  and  K d 1 ), whereas the electrical power grid in area b has  K p 2  and  K d 2  for the tuning process. The output of the outer loop is fed into the inner FOPID loop. The representations of the error signals  E a 2 ( s )  and  E b 2 ( s )  are as follows:
E a 2 ( s ) = Y a 1 ( s ) Δ P t i e Δ f a E b 2 ( s ) = Y b 1 ( s ) Δ P t i e Δ f b
The representations of the FOPID loops are as follows:
Y a 2 ( s ) = [ K p 3 + K i 1 s λ 1 + K d 3 s μ 1 ] . E a 2 ( s ) Y b 2 ( s ) = [ K p 4 + K i 2 s λ 2 + K d 4 s μ 2 ] . E b 2 ( s )
From (27), the electric power grid in area a has five tunable control parameters ( K p 3 K i 1 K d 3 λ 1  and  μ 1 ), and the electrical power grid in area b has  K p 4 K i 2 K d 4 λ 2  and  μ 2  as tunable parameters.
Figure 3. Proposed 1 + PD/FOPID controller.
Figure 3. Proposed 1 + PD/FOPID controller.
Fractalfract 07 00603 g003

4. The Proposed Design Optimization

4.1. MRFO Optimizer

The MRFO is a metaheuristic optimizer that belong to foraging strategies followed by manta rays during the catching process of their prey [58]. It is mainly composed of three foraging processes (chain, cyclone and somersault foraging). During chain foraging, manta rays consider highly concentrated plankton, which represent the desired optimization objectives and tracker. This, in turn, makes them align with the foraging chain. In which, everyone is directed towards food by the manta rays within its front. Then, an updated process of individuals is obtained from the best solution in each iterations. The chain process of foraging is mathematically expressed as [58]:
x i t + 1 = x i t + r . ( x b e s t t x i t ) + ω 1 ( x b e s t t x i t ) , i = 1 x i t + r . ( x i 1 t x i t ) + ω 1 ( x b e s t t x i t ) , i = 2 : N
where  x i t  stands for the  i t h  position for the current iteration t, r represents a random vector,  x b e s t t  stands for best solution in the  t t h  iteration, N is the number of manta rays, and  ω 1  is the weighting coefficient, expressed as:
ω 1 = 2 × r × | log ( r ) |
From (28), the individuals’ positions are calculated for all individuals except the first  ( i 1 ) t h  individual and the best individual  x b e s t t . After plankton patch position determinations using manta rays, a chain is formed by their combination and they swim in a spiral shape towards their prey. Additionally, individuals swim towards the front-sided manta ray. The cyclone foraging process is expressed as:
x i t + 1 = x b e s t + r . ( x i 1 t x i t ) + e b ω . cos ( 2 π ω ) . ( x b e s t x i t ) y i t + 1 = y b e s t + r . ( y i 1 t y i t ) + e b ω . cos ( 2 π ω ) . ( y b e s t y i t )
where  ω  stands for a random number. Then, the cyclone process foraging is as follows:
x i t + 1 = x b e s t + r . ( x b e s t t x i t ) + ω 2 ( x b e s t t x i t ) , i = 1 x b e s t + r . ( x i 1 t x i t ) + ω 2 ( x b e s t t x i t ) , i = 2 : N
where  ω 2  represents the weighting factor of cyclone foraging and is determined as:
ω 2 = 2 e r 1 ( T t + 1 T ) . sin ( 2 π r 1 )
where t stands for the current iteration, T stands for the maximum iterations number, and  r 1  stands for a random number. The improved exploitation is obtained through the cyclone foraging process to determine the best solution region. This is because all existing manta rays contribute to the food search processes based on their reference positions. Additionally, the exploitation process is also enhanced by forcing the individuals to search for new positions located away from the current best position. The random position in the search space is determined as:
x r a n d = L b + r . ( U b L b )
x i t + 1 = x r a n d + r . ( x r a n d t x i t ) + β ( x r a n d t x i t ) , i = 1 x r a n d + r . ( x i 1 t x i t ) + β ( x r a n d t x i t ) , i = 2 : N
where  U b  and  L b  are the upper/lower limits, respectively, of the desired variables,  x r a n d  is an assigned random position in the search space. In somersault foraging, the food is recognized in this stage as a hinge. Wherein, each manta ray swims backwards and forwards around the food hinge, and then they tumble to the new position. The somersault foraging process is performed as follows:
x i t + 1 = x i t + S . ( r 2 . x b e s t t r 3 . x i t ) , i = 1 , 2 , . . . , N
where S is the somersault factor employed to determine the somersault range for the manta rays, and  r 2  and  r 3  are random numbers. A flowchart representation of the MRFO stages is shown in Figure 4.

4.2. Design Optimization

The MRFO is proposed in this paper to determine the best parameters for the proposed 1 + PD/FOPID for both electrical power grids. The main driving objectives for the optimum parameters include minimizing the existing fluctuations in the frequency in both areas as well as the tie-line power among the electrical power grids. A measure of frequency deviation ( Δ f a  of area a and  Δ f b  of area b) and tie-line power ( Δ P t i e  between electrical power grids) is performed and fed into the optimization process. Then, they are combined in a single-objective function to drive the optimization process for simulation time  t s  while taking problem constraints into consideration in the process. The employed objective functions can be expressed as follows:
I S E = integral squared error = 0 t s ( ( Δ f a ) 2 + ( Δ f b ) 2 + ( Δ P t i e ) 2 ) d t
I T S E = integral time squared error = 0 t s t . ( ( Δ f a ) 2 + ( Δ f b ) 2 + ( Δ P t i e ) 2 ) d t
I A E = integral absolute error = 0 t s ( a b s ( Δ f a ) + a b s ( Δ f b ) + a b s ( Δ P t i e ) ) d t
I T A E = integral time absolute error = 0 t s t . ( a b s ( Δ f a ) + a b s ( Δ f b ) + a b s ( Δ P t i e ) ) d t
The ISE and IAE are based on using the integration of the square and absolute error values, respectively, within the simulation time  t s . The ISE provides better consideration of the large error values due to the square operation (when errors are more than 1). In addition, IAE provides equal consideration for large and low error values. Whereas ITSE and ITAE consider the time during the integration compared to the ISE and IAE objectives, which leads to lower/zero steady state error compared to ISE and IAE. The four objectives are considered in this paper to provide a comprehensive comparison of the studied controllers.
The proposed design optimization process based on MRFO is summarized in Figure 5. The employed optimization constraints in our proposed design process are:
k p m i n k p 1 , k p 2 , k p 3 , k p 4 k p m a x k i m i n k i 1 , k i 2 k i m a x k d m i n k d 1 , k d 2 , k d 3 , k d 4 k d m a x λ m i n λ 1 , λ 2 λ m a x μ m i n μ 1 , μ 2 μ m a x
where the lower/upper constraints are represented by  m i n , and  m a x , respectively, for the proposed 1 + PD/FOPID LFC. The used minimum constraints for  k p m i n k i m i n , and  k d m i n  were set at zero, and  k p m a x k i m a x , and  k d m a x  were set at five during the proposed optimization stages. The set minimum values for  λ m i n  and  μ m i n  were zero and for  λ m a x  and  μ m a x  were 1.

5. Simulation Results and Performance Verification

The proposed system and proposed design optimization were implemented using the MATLAB R2021a joint m-file and Simulink platforms. The objective function and optimizers were programmed using m-file and linked with Simulink platform. The proposed design optimization process is based on using 20 populations with a maximum of 100 iterations for all the studied optimizers. The same process was used for the design of all the compared controllers for a fair comparison. The two-area electrical power grid system was tested, and the performance of the proposed 1 + PD/FOPID controller was evaluated and compared with the PID, FOPID, and PD/FOPID controllers. Moreover, the convergence performance of the MRFO optimizer was compared with some metaheuristic optimization algorithms from the literature.
The considered optimizers are the GA, PSO, and MPA. The design optimization was made using a personal computer with an Intel Core i7, CPU of 2.9 GHz, and a 64-bit version. Figure 6 compares the ISE and IAE convergence curves of the studied optimizers, whereas Figure 7 shows the ITSE and ITAE comparisons. It has become evident that the MRFO-based design optimization possesses the best convergence with the lowest objective function for all the studied ISE, IAE, ITSE, and ITAE objective functions. In addition, MRFO achieves a very fast conversion with better determination of the control parameters compared with the other studied optimization algorithms. Table 3 summarizes the obtained controllers’ parameters using the proposed design optimization process. The considered test scenarios are organized as follows:
  • Scenario (1): Impacts of the stepped load perturbations (SLP);
  • Scenario (2): Impacts of multiple SLPs on the two interconnected electrical power grids;
  • Scenario (3): Impacts of multiple connection/disconnection of RESs;
  • Scenario (4): Impacts of randomly varying loads;
  • Scenario (5): Joint impacts of RES fluctuations with various load-type variations.
Table 3. The optimum controllers’ parameters using MRFO design optimization.
Table 3. The optimum controllers’ parameters using MRFO design optimization.
ControllerAreaParameters
  k p 1   k p 2   k i 1   k d 1   k d 2   λ 1   μ 1
PIDArea a1.90621.85471.8637
Area b0.88080.28230.4233
FOPIDArea a1.81841.5670.99690.830.56
Area b1.98091.1891.94970.890.73
PD/FOPIDArea a4.37494.98371.92313.11521.64030.910.76
Area b2.58394.77020.95440.70113.31580.620.93
1 + PD/FOPIDArea a4.52813.27513.40074.22124.94970.970.82
Area b3.71130.63611.63414.31582.94660.770.91
Figure 6. Convergence curves for the proposed design optimization; (a) ISE; (b) IAE.
Figure 6. Convergence curves for the proposed design optimization; (a) ISE; (b) IAE.
Fractalfract 07 00603 g006
Figure 7. Convergence curves for the proposed design optimization; (a) ITSE; (b) ITAE.
Figure 7. Convergence curves for the proposed design optimization; (a) ITSE; (b) ITAE.
Fractalfract 07 00603 g007

5.1. Results of Scenario (1)

Figure 8 shows the obtained results during Scenario (1) with an SLP of 2%. The proposed 1 + PD/FOPID LFC achieves the best transient response compared with the studied controllers. The proposed 1 + PD/FOPID LFC has a maximum undershoot (MU) in  Δ f a  of 0.0018 compared with 0.0101, 0.0061, and 0.0044 under the PID, FOPID, and PD/FOPID, respectively. Moreover, the MU in  Δ f b  was 0.0002 under the proposed 1 + PD/FOPID controller compared with 0.0071, 0.0023, and 0.0016 under the PID, FOPID, and PD/FOPID, respectively. Furthermore, the proposed 1 + PD/FOPID had the lowest settling time (ST) compared with the studied controllers. In addition, the effect of the proposed controller on the thermal, hydraulic, and EV performances is depicted on Figure 9. It can be noted from this figure that the output powers from the thermal power unit and EVs in area a do not exceed their maximum bounds based on different control signals from the frequency variations or area control error. This is reflected as an improvement in the stability behaviour and response of the studied electrical power grids with various expected sudden load changes. For this, the proposed controller succeeded at preserving a better response with subjected load changes compared with the studied controllers.

5.2. Results of Scenario (2)

From another side, Scenario (2) was made using multiple SLPs in different areas. Figure 10 shows the applied load powers in different electrical power grids as a test scenario. Figure 11a–c shows the obtained performance response in this scenario. The response of the frequency and tie-line power transients shows the proposed controller with better transients in all the tested SLP changes in this scenario. From the measured response, the proposed 1 + PD/FOPID has the best frequency deviation response in areas a and b. The proposed 1 + PD/FOPID has an MU in area a of 0.0005 and in area b of 0.0011. Whereas the PID has values of 0.0078 and 0.0103 in areas a and b, respectively; FOPID has values of 0.0051 and 0.0081 in areas a and b, respectively; and PD/FOPID has values of 0.0037 and 0.00581 in areas a and b, respectively. Therefore, the proposed 1 + PD/FOPID achieves the lowest peaks during this scenario. Furthermore, the superior impact of the proposed controller on the performance of the hydraulic generation unit and EVs outputs and its static of charge can be seen in Figure 12, which shows that both of them can regulate the system frequency without exceeding their maximum limits. In addition, it can be observed that there is a cross-coupling between the areas during the transient state, and hence each area produces its own power at a steady state.

5.3. Results at Scenario (3)

In this scenario, multiple RES connections/disconnections have been made to test the proposed controller. Figure 13 shows the PV and wind powers in this scenario. The wind is connected at time 30 s and disconnected at time 80 s. Whereas the PV power is connected at time 50 s. In addition, an SLP is made at the start of the scenario. The obtained results are shown in Figure 14a–c for this scenario. It can be seen that connecting/disconnecting wind/PV affects the response of the system due to its participation level. The proposed 1 + PD/FOPID has the best performance metrics in this scenario, and the PIF has the worst response. For instance, the obtained MO values in area a at 30 s are 0.1202, 0.0596, 0.0456, and 0.0085 for the PID, FOPID, PPD/FOPID, and proposed 1 + PD/FOPID LFC, respectively. The measured performance metrics for all the studied scenarios are shown in Table 4. In which, a performance enhancement is observed when using the proposed control and design optimization method.

5.4. Results at Scenario (4)

An important factor to be considered is the characteristics of the connected electrical load in the electrical power grids. The load varies all day, and hence so too do the expected different demands of power in each moment. These variations are often reflected as fluctuations in the operating frequency of the electrical power grids. Therefore, in this tested scenario, a randomly changing electrical loading is considered as shown in Figure 15. The associated results for this scenario are shown in Figure 16. It can be seen that the proposed 1 + PD/FOPID has the lowest peak fluctuations with varying load profiles. Whereas the PID-based LFC has the highest level of fluctuations during this level. The PD/FOPID comes in second place in terms of improving the electrical power grid response; this followed by the response of the FOPID control.

5.5. Results at Scenario (5)

The interconnected electrical power grids are subjected to joint intermittency of the connected RESs with the connected electrical loading. In this scenario, the domestic and industrial load profiles are considered and studied as shown in Figure 17. Additionally, the fluctuated RESs are included in this scenario with their participation levels and connection/disconnection events in the scenario, as shown in Figure 18. At time 0 s, all the renewable sources and loads are connected, which represents the worst-case scenario to test all the studied controllers and the proposed design optimization method. Figure 19 presents the obtained results in this scenario. The PID at the start of scenario with all RESs and loading step has 0.0927 MO in  Δ f a  of area a and 0.1463 MO in  Δ f b  of area b. Furthermore, regarding the deviations in  Δ P t i e , the PID has an MU of 0.0461 in this scenario. Therefore, it has the lowest performance of the four studied controllers. From another side, the proposed controller has an MU of 0.0334 and 0.0879 in  Δ f a , and  Δ f b , respectively. In addition, the proposed control has deviations in  Δ P t i e  of 0.0026 in this scenario. Therefore, the best performance is achieved through the proposed 1 + PD/FOPID LFC method. A full measure of system performance during the five scenarios through various metrics is detailed in Table 4.
Table 4. Measurements of the settling time (ST), peak undershoot (PU), and peak overshoot (PO) for the studied scenarios (where FU stands for a fluctuated condition.
Table 4. Measurements of the settling time (ST), peak undershoot (PU), and peak overshoot (PO) for the studied scenarios (where FU stands for a fluctuated condition.
ScenarioController Δ f 1 Δ f 2 Δ P t i e
POPUST (s)POPUST (s)POPUST (s)
No. 1 at 0 sPID0.00080.0101130.00110.007190.00060.002719
FOPID0.00150.0061110.00140.0023110.00010.002316
PD/FOPID0.00020.00448-0.0016100.00040.001410
1 + PD/FOPID0.00010.00184-0.00023- 9.3 × 10 5 3
No. 2 at 30 sPID0.00030.0078>20 s0.00050.0103>20 s0.00350.0001>20 s
FOPID0.00060.0051190.00060.0081>20 s0.00310.0009>20 s
PD/FOPID-0.0037220.00120.0058190.00290.000520
1 + PD/FOPID-0.00057-0.001150.0003-6
No. 3 at 30 sPID0.12020.0111FU0.07390.0209FU0.02640.0039FU
FOPID0.05960.0155FU0.02360.0132FU0.02430.0058FU
PD/FOPID0.04560.0174130.02050.0043110.01940.0012FU
1 + PD/FOPID0.00850.002170.0031-50.0015-5
No. 4PID0.05690.0534FU0.05370.0608FU0.01790.0192FU
FOPID0.02390.0272FU0.02190.0209FU0.01390.0124FU
PD/FOPID0.00940.0088FU0.00760.0092FU0.00940.0089FU
1 + PD/FOPID0.00170.0015FU0.00120.0013FU0.00070.0006FU
No. 5 at 0 sPID0.09270.0073FU0.14630.0127FU0.00770.0461FU
FOPID0.05270.0004FU0.08790.0092FU0.00190.0285FU
PD/FOPID0.04710.0141FU0.07110.0065FU0.00220.0221FU
1 + PD/FOPID0.0334-150.0191-19-0.002623

6. Conclusions

An improved fractional-order controller based on a cascaded 1 + PD/FOPID control was proposed in this paper with MRFO-based design optimization to regulate the frequency in interconnected electrical grids. The proposed controller is advantageous at mitigating disturbances at a wide range of frequencies due to employing two cascaded control loops. Additionally, the employment of the MRFO and design optimization process leads to simultaneous design and determination of the best control parameters. The proposed 1 + PD/FOPID and MRFO-based design optimization were implemented and simulated in MATLAB. Various scenarios of load power changes and renewable power connection/interconnection scenarios were considered in the performance investigation. In addition, the PID, FOPID and PD/FOPID LFCs were compared with the proposed controller. Results verify the reduced peak overshoot and settling times under the proposed 1 + PD/FOPID LFC compared with the studied controllers. Future research includes the use of more practical modelling of electric power systems (integer and fractional-order modelling), stability analysis with system non-linearities, and comprehensive comparisons of existing LFC methods.

Author Contributions

Conceptualization, F.F.M.E.-S., M.H.A., E.A.M. and M.A.; data curation, A.S.A., S.Z.A., E.A.M. and M.H.A.; formal analysis, F.F.M.E.-S., S.Z.A., A.S.A. and M.A.; funding acquisition, M.H.A., A.S.A., F.F.M.E.-S. and S.Z.A.; investigation, F.F.M.E.-S., M.H.A., E.A.M. and M.A.; methodology, F.F.M.E.-S., M.H.A., E.A.M., A.S.A. and S.Z.A.; project administration, S.Z.A., E.A.M., F.F.M.E.-S. and M.H.A.; resources, F.F.M.E.-S., A.S.A., E.A.M. and M.A.; software, F.F.M.E.-S., M.H.A., A.S.A. and M.A.; supervision, A.S.A., S.Z.A., E.A.M. and M.A.; validation, F.F.M.E.-S., M.H.A., E.A.M., M.A. and S.Z.A.; visualization, F.F.M.E.-S., E.A.M., M.H.A. and M.A.; writing—original draft, F.F.M.E.-S., M.H.A., A.S.A., E.A.M. and M.A.; writing—review and editing, F.F.M.E.-S., M.H.A. and M.A. All authors have read and agreed to the published version of the manuscript.

Funding

The authors extend their appreciation to the Deputyship for Research and Innovation, Ministry of Education in Saudi Arabia for funding this research work through the Project number (IF2/PSAU/2022/01/23085).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Appendix A.1

For power system representation in areas a and b, the representation of  Δ f a  and  Δ f b  based on the system representation in Figure 2 are represented as:
Δ f a = ( Δ P g a + Δ P W T Δ P E V a Δ P l a Δ P t i e ) 1 2 H a s + D a
Δ f b = ( Δ P g b + Δ P P V Δ P E V b Δ P l b + Δ P t i e ) 1 2 H b s + D b
By multiplying both parts in (A1) by  2 H a s + D a  and also both parts in (A2) by  2 H b s + D b , we obtain:
2 H a s Δ f a + D a Δ f a = Δ P g a + Δ P W T Δ P E V a Δ P l a Δ P t i e
2 H b s Δ f b + D b Δ f b = Δ P g b + Δ P P V Δ P E V b Δ P l b + Δ P t i e
where the term  s Δ f a  represents  Δ f ˙ a  and  s Δ f b  represents  Δ f ˙ b . Thus, (A3) and (A4) become:
Δ f ˙ a = D a 2 H a Δ f a + 1 2 H a Δ P g a + 1 2 H a Δ P W T 1 2 H a Δ P E V a 1 2 H a Δ P l a 1 2 H a Δ P t i e
Δ f ˙ b = D b 2 H b Δ f b + 1 2 H b Δ P g b + 1 2 H b Δ P P V 1 2 H b Δ P E V b 1 2 H b Δ P l b + 1 2 H b Δ P t i e

Appendix A.2

For thermal generation, it is represented by:
Δ P g a = 1 T t s + 1 Δ P g a 1
Δ P g a 1 = 1 T g s + 1 ( A C E o a 1 R a Δ f a )
By simplifying (A7) and (A8), and  s Δ P g a  is replaced with  Δ P ˙ g a  and  s Δ P g a 1  is replaced with  Δ P ˙ g a 1 , we obtain:
Δ P ˙ g a = 1 T t Δ P g a + 1 T t Δ P g a 1
Δ P ˙ g a 1 = 1 T g Δ P g a 1 + 1 T g A C E o a 1 T g R a Δ f a
For hydraulic generation, it is represented by:
Δ P g b = T w s + 1 0.5 T w s + 1 Δ P g b 1
Δ P g b 1 = T R s + 1 T 2 s + 1 Δ P g b 2
Δ P g b 2 = 1 T 1 s + 1 ( A C E o b 1 R b Δ f b )
By simplifying (A11)–(A13), and  s Δ P g b  is replaced with  Δ P ˙ g b s Δ P g b 1  is replaced with  Δ Δ ˙ P g b 1  and  s Δ P g b 2  is replaced with  Δ P ˙ g b 2 , we obtain:
Δ P ˙ g b = 2 T w Δ P g b + 2 T 2 + 2 T w T 2 T w Δ P g b 1 + 2 T R 2 T 1 T 1 T 2 Δ P g b 2 + 2 T R T 1 T 2 R b Δ f b + 2 T R T 1 T 2 A C E o b
Δ P ˙ g b 1 = 1 T 2 Δ P g b 1 + T 1 T R T 1 T 2 Δ P g b 2 T R T 1 T 2 R b Δ f b T R T 1 T 2 A C E o b
Δ P ˙ g b 2 = 1 T 1 Δ P g b 2 + 1 T 1 A C E o b 1 T 1 R b Δ f b

Appendix A.3

For wind generation, it is represented by:
Δ P W T = K W T T W T s + 1 P W T
By simplifying (A17), and  s Δ P W T  is replaced with  Δ P ˙ W T , we obtain:
Δ P ˙ W T = 1 T W T Δ P W T + K W T T W T P W T
For PV generation, it is represented by:
Δ P P V = K P V T P V s + 1 P W T
By simplifying (A19), and  s Δ P P V  is replaced with  Δ P ˙ P V , we obtain:
Δ P ˙ P V = 1 T P V Δ P P V + K P V T P V P P V

References

  1. Blaabjerg, F.; Yang, Y.; Kim, K.A.; Rodriguez, J. Power Electronics Technology for Large-Scale Renewable Energy Generation. Proc. IEEE 2023, 111, 335–355. [Google Scholar] [CrossRef]
  2. Biswas, S.; Mahata, S.; Roy, P.K.; Chatterjee, K. Application of Empirical Bode Analysis for Delay-Margin Evaluation of Fractional-Order PI Controller in a Renewable Distributed Hybrid System. Fractal Fract. 2023, 7, 119. [Google Scholar] [CrossRef]
  3. Sami, I.; Ullah, S.; Khan, L.; Al-Durra, A.; Ro, J.S. Integer and Fractional-Order Sliding Mode Control Schemes in Wind Energy Conversion Systems: Comprehensive Review, Comparison, and Technical Insight. Fractal Fract. 2022, 6, 447. [Google Scholar] [CrossRef]
  4. Salama, H.S.; Said, S.M.; Aly, M.; Vokony, I.; Hartmann, B. Studying Impacts of Electric Vehicle Functionalities in Wind Energy-Powered Utility Grids With Energy Storage Device. IEEE Access 2021, 9, 45754–45769. [Google Scholar] [CrossRef]
  5. Alilou, M.; Azami, H.; Oshnoei, A.; Mohammadi-Ivatloo, B.; Teodorescu, R. Fractional-Order Control Techniques for Renewable Energy and Energy-Storage-Integrated Power Systems: A Review. Fractal Fract. 2023, 7, 391. [Google Scholar] [CrossRef]
  6. Daraz, A.; Malik, S.A.; Basit, A.; Aslam, S.; Zhang, G. Modified FOPID Controller for Frequency Regulation of a Hybrid Interconnected System of Conventional and Renewable Energy Sources. Fractal Fract. 2023, 7, 89. [Google Scholar] [CrossRef]
  7. Fayek, H.H. 5G Poor and Rich Novel Control Scheme Based Load Frequency Regulation of a Two-Area System with 100% Renewables in Africa. Fractal Fract. 2020, 5, 2. [Google Scholar] [CrossRef]
  8. Ahmed, E.M.; Mohamed, E.A.; Selim, A.; Aly, M.; Alsadi, A.; Alhosaini, W.; Alnuman, H.; Ramadan, H.A. Improving load frequency control performance in interconnected power systems with a new optimal high degree of freedom cascaded FOTPID-TIDF controller. Ain Shams Eng. J. 2023, 14, 102207. [Google Scholar] [CrossRef]
  9. Mohamed, E.A.; Aly, M.; Elmelegi, A.; Ahmed, E.M.; Watanabe, M.; Said, S.M. Enhancement the Frequency Stability and Protection of Interconnected Microgrid Systems Using Advanced Hybrid Fractional Order Controller. IEEE Access 2022, 10, 111936–111961. [Google Scholar] [CrossRef]
  10. Kavikumar, R.; Kwon, O.M.; Lee, S.H.; Sakthivel, R. Input-output finite-time IT2 fuzzy dynamic sliding mode control for fractional-order nonlinear systems. Nonlinear Dyn. 2022, 108, 3745–3760. [Google Scholar] [CrossRef]
  11. Kavikumar, R.; Sakthivel, R.; Kwon, O.M.; Selvaraj, P. Robust tracking control design for fractional-order interval type-2 fuzzy systems. Nonlinear Dyn. 2022, 107, 3611–3628. [Google Scholar] [CrossRef]
  12. Gheisarnejad, M.; Khooban, M.H.; Dragicevic, T. The Future 5G Network-Based Secondary Load Frequency Control in Shipboard Microgrids. IEEE J. Emerg. Sel. Top. Power Electron. 2020, 8, 836–844. [Google Scholar] [CrossRef]
  13. Gheisarnejad, M.; Karimaghaie, P.; Boudjadar, J.; Khooban, M.H. Real-Time Cellular Wireless Sensor Testbed for Frequency Regulation in Smart Grids. IEEE Sens. J. 2019, 19, 11656–11665. [Google Scholar] [CrossRef]
  14. Ahmed, E.M.; Selim, A.; Alnuman, H.; Alhosaini, W.; Aly, M.; Mohamed, E.A. Modified Frequency Regulator Based on TIλ-TDμFF Controller for Interconnected Microgrids with Incorporating Hybrid Renewable Energy Sources. Mathematics 2022, 11, 28. [Google Scholar] [CrossRef]
  15. Aly, M.; Mohamed, E.A.; Noman, A.M.; Ahmed, E.M.; El-Sousy, F.F.M.; Watanabe, M. Optimized Non-Integer Load Frequency Control Scheme for Interconnected Microgrids in Remote Areas with High Renewable Energy and Electric Vehicle Penetrations. Mathematics 2023, 11, 2080. [Google Scholar] [CrossRef]
  16. Dahab, Y.A.; Abubakr, H.; Mohamed, T.H. Adaptive Load Frequency Control of Power Systems Using Electro-Search Optimization Supported by the Balloon Effect. IEEE Access 2020, 8, 7408–7422. [Google Scholar] [CrossRef]
  17. Ewais, A.M.; Elnoby, A.M.; Mohamed, T.H.; Mahmoud, M.M.; Qudaih, Y.; Hassan, A.M. Adaptive frequency control in smart microgrid using controlled loads supported by real-time implementation. PLoS ONE 2023, 18, e0283561. [Google Scholar] [CrossRef]
  18. Arora, K.; Kumar, A.; Kamboj, V.K.; Prashar, D.; Shrestha, B.; Joshi, G.P. Impact of Renewable Energy Sources into Multi Area Multi-Source Load Frequency Control of Interrelated Power System. Mathematics 2021, 9, 186. [Google Scholar] [CrossRef]
  19. Gupta, D.K.; Soni, A.K.; Jha, A.V.; Mishra, S.K.; Appasani, B.; Srinivasulu, A.; Bizon, N.; Thounthong, P. Hybrid Gravitational–Firefly Algorithm-Based Load Frequency Control for Hydrothermal Two-Area System. Mathematics 2021, 9, 712. [Google Scholar] [CrossRef]
  20. Yousri, D.; Babu, T.S.; Fathy, A. Recent methodology based Harris Hawks optimizer for designing load frequency control incorporated in multi-interconnected renewable energy plants. Sustain. Energy Grids Netw. 2020, 22, 100352. [Google Scholar] [CrossRef]
  21. Salama, H.S.; Magdy, G.; Bakeer, A.; Vokony, I. Adaptive coordination control strategy of renewable energy sources, hydrogen production unit, and fuel cell for frequency regulation of a hybrid distributed power system. Prot. Control Mod. Power Syst. 2022, 7, 34. [Google Scholar] [CrossRef]
  22. Shabani, H.; Vahidi, B.; Ebrahimpour, M. A robust PID controller based on imperialist competitive algorithm for load-frequency control of power systems. ISA Trans. 2013, 52, 88–95. [Google Scholar] [CrossRef] [PubMed]
  23. Youssef, A.R.; Mallah, M.; Ali, A.; Shaaban, M.F.; Mohamed, E.E.M. Enhancement of Microgrid Frequency Stability Based on the Combined Power-to-Hydrogen-to-Power Technology under High Penetration Renewable Units. Energies 2023, 16, 3377. [Google Scholar] [CrossRef]
  24. Sahu, R.K.; Panda, S.; Rout, U.K.; Sahoo, D.K. Teaching learning based optimization algorithm for automatic generation control of power system using 2-DOF PID controller. Int. J. Electr. Power Energy Syst. 2016, 77, 287–301. [Google Scholar] [CrossRef]
  25. Kouba, N.E.Y.; Menaa, M.; Hasni, M.; Boudour, M. Optimal load frequency control based on artificial bee colony optimization applied to single, two and multi-area interconnected power systems. In Proceedings of the 2015 3rd International Conference on Control, Engineering Information Technology (CEIT), Tlemcen, Algeria, 25–27 May 2015. [Google Scholar] [CrossRef]
  26. Abid, S.; El-Rifaie, A.M.; Elshahed, M.; Ginidi, A.R.; Shaheen, A.M.; Moustafa, G.; Tolba, M.A. Development of Slime Mold Optimizer with Application for Tuning Cascaded PD-PI Controller to Enhance Frequency Stability in Power Systems. Mathematics 2023, 11, 1796. [Google Scholar] [CrossRef]
  27. Zhang, G.; Daraz, A.; Khan, I.A.; Basit, A.; Khan, M.I.; Ullah, M. Driver Training Based Optimized Fractional Order PI-PDF Controller for Frequency Stabilization of Diverse Hybrid Power System. Fractal Fract. 2023, 7, 315. [Google Scholar] [CrossRef]
  28. Singh, B.; Slowik, A.; Bishnoi, S.K.; Sharma, M. Frequency Regulation Strategy of Two-Area Microgrid System with Electric Vehicle Support Using Novel Fuzzy-Based Dual-Stage Controller and Modified Dragonfly Algorithm. Energies 2023, 16, 3407. [Google Scholar] [CrossRef]
  29. Mohanty, B.; Panda, S.; Hota, P. Controller parameters tuning of differential evolution algorithm and its application to load frequency control of multi-source power system. Int. J. Electr. Power Energy Syst. 2014, 54, 77–85. [Google Scholar] [CrossRef]
  30. Dash, P.; Saikia, L.C.; Sinha, N. Automatic generation control of multi area thermal system using Bat algorithm optimized PD–PID cascade controller. Int. J. Electr. Power Energy Syst. 2015, 68, 364–372. [Google Scholar] [CrossRef]
  31. Raju, M.; Saikia, L.C.; Sinha, N. Automatic generation control of a multi-area system using ant lion optimizer algorithm based PID plus second order derivative controller. Int. J. Electr. Power Energy Syst. 2016, 80, 52–63. [Google Scholar] [CrossRef]
  32. Latif, A.; Hussain, S.M.S.; Das, D.C.; Ustun, T.S. Optimization of Two-Stage IPD-(1I) Controllers for Frequency Regulation of Sustainable Energy Based Hybrid Microgrid Network. Electronics 2021, 10, 919. [Google Scholar] [CrossRef]
  33. Gheisarnejad, M. An effective hybrid harmony search and cuckoo optimization algorithm based fuzzy PID controller for load frequency control. Appl. Soft Comput. 2018, 65, 121–138. [Google Scholar] [CrossRef]
  34. Prakash, S.; Sinha, S. Simulation based neuro-fuzzy hybrid intelligent PI control approach in four-area load frequency control of interconnected power system. Appl. Soft Comput. 2014, 23, 152–164. [Google Scholar] [CrossRef]
  35. Latif, A.; Paul, M.; Das, D.C.; Hussain, S.M.S.; Ustun, T.S. Price Based Demand Response for Optimal Frequency Stabilization in ORC Solar Thermal Based Isolated Hybrid Microgrid under Salp Swarm Technique. Electronics 2020, 9, 2209. [Google Scholar] [CrossRef]
  36. Hussain, I.; Das, D.C.; Latif, A.; Sinha, N.; Hussain, S.S.; Ustun, T.S. Active power control of autonomous hybrid power system using two degree of freedom PID controller. Energy Rep. 2022, 8, 973–981. [Google Scholar] [CrossRef]
  37. Bakeer, A.; Magdy, G.; Chub, A.; Jurado, F.; Rihan, M. Optimal Ultra-Local Model Control Integrated with Load Frequency Control of Renewable Energy Sources Based Microgrids. Energies 2022, 15, 9177. [Google Scholar] [CrossRef]
  38. Yakout, A.H.; AboRas, K.M.; Kotb, H.; Alharbi, M.; Shouran, M.; Samad, B.A. A Novel Ultra Local Based-Fuzzy PIDF Controller for Frequency Regulation of a Hybrid Microgrid System with High Renewable Energy Penetration and Storage Devices. Processes 2023, 11, 1093. [Google Scholar] [CrossRef]
  39. Ayas, M.S.; Sahin, E. FOPID controller with fractional filter for an automatic voltage regulator. Comput. Electr. Eng. 2021, 90, 106895. [Google Scholar] [CrossRef]
  40. Fathy, A.; Alharbi, A.G. Recent Approach Based Movable Damped Wave Algorithm for Designing Fractional-Order PID Load Frequency Control Installed in Multi-Interconnected Plants With Renewable Energy. IEEE Access 2021, 9, 71072–71089. [Google Scholar] [CrossRef]
  41. Latif, A.; Hussain, S.S.; Das, D.C.; Ustun, T.S.; Iqbal, A. A review on fractional order (FO) controllers’ optimization for load frequency stabilization in power networks. Energy Rep. 2021, 7, 4009–4021. [Google Scholar] [CrossRef]
  42. Singh, K.; Amir, M.; Ahmad, F.; Refaat, S.S. Enhancement of Frequency Control for Stand-Alone Multi-Microgrids. IEEE Access 2021, 9, 79128–79142. [Google Scholar] [CrossRef]
  43. Zaheeruddin; Singh, K. Load frequency regulation by de-loaded tidal turbine power plant units using fractional fuzzy based PID droop controller. Appl. Soft Comput. 2020, 92, 106338. [Google Scholar] [CrossRef]
  44. Oshnoei, S.; Aghamohammadi, M.; Oshnoei, S.; Oshnoei, A.; Mohammadi-Ivatloo, B. Provision of Frequency Stability of an Islanded Microgrid Using a Novel Virtual Inertia Control and a Fractional Order Cascade Controller. Energies 2021, 14, 4152. [Google Scholar] [CrossRef]
  45. Peddakapu, K.; Srinivasarao, P.; Mohamed, M.; Arya, Y.; Kishore, D.K. Stabilization of frequency in Multi-Microgrid system using barnacle mating Optimizer-based cascade controllers. Sustain. Energy Technol. Assess. 2022, 54, 102823. [Google Scholar] [CrossRef]
  46. Arya, Y.; Kumar, N.; Dahiya, P.; Sharma, G.; Çelik, E.; Dhundhara, S.; Sharma, M. Cascade-I ΛDμN controller design for AGC of thermal and hydro-thermal power systems integrated with renewable energy sources. IET Renew. Power Gener. 2021, 15, 504–520. [Google Scholar] [CrossRef]
  47. Malik, S.; Suhag, S. A Novel SSA Tuned PI-TDF Control Scheme for Mitigation of Frequency Excursions in Hybrid Power System. Smart Sci. 2020, 8, 202–218. [Google Scholar] [CrossRef]
  48. Priyadarshani, S.; Subhashini, K.R.; Satapathy, J.K. Pathfinder algorithm optimized fractional order tilt-integral-derivative (FOTID) controller for automatic generation control of multi-source power system. Microsyst. Technol. 2020, 27, 23–35. [Google Scholar] [CrossRef]
  49. Oshnoei, A.; Khezri, R.; Muyeen, S.M.; Oshnoei, S.; Blaabjerg, F. Automatic Generation Control Incorporating Electric Vehicles. Electr. Power Compon. Syst. 2019, 47, 720–732. [Google Scholar] [CrossRef]
  50. Sahu, R.K.; Panda, S.; Biswal, A.; Sekhar, G.C. Design and analysis of tilt integral derivative controller with filter for load frequency control of multi-area interconnected power systems. ISA Trans. 2016, 61, 251–264. [Google Scholar] [CrossRef]
  51. Elkasem, A.H.; Khamies, M.; Hassan, M.H.; Nasrat, L.; Kamel, S. Utilizing controlled plug-in electric vehicles to improve hybrid power grid frequency regulation considering high renewable energy penetration. Int. J. Electr. Power Energy Syst. 2023, 152, 109251. [Google Scholar] [CrossRef]
  52. Ahmed, M.; Magdy, G.; Khamies, M.; Kamel, S. Modified TID controller for load frequency control of a two-area interconnected diverse-unit power system. Int. J. Electr. Power Energy Syst. 2022, 135, 107528. [Google Scholar] [CrossRef]
  53. Khokhar, B.; Dahiya, S.; Parmar, K.P.S. A Robust Cascade Controller for Load Frequency Control of a Standalone Microgrid Incorporating Electric Vehicles. Electr. Power Compon. Syst. 2020, 48, 711–726. [Google Scholar] [CrossRef]
  54. Das, D.C.; Roy, A.; Sinha, N. GA based frequency controller for solar thermal–diesel–wind hybrid energy generation/energy storage system. Int. J. Electr. Power Energy Syst. 2012, 43, 262–279. [Google Scholar] [CrossRef]
  55. Ray, P.K.; Mohanty, S.R.; Kishor, N. Proportional–integral controller based small-signal analysis of hybrid distributed generation systems. Energy Convers. Manag. 2011, 52, 1943–1954. [Google Scholar] [CrossRef]
  56. Ahmed, E.M.; Mohamed, E.A.; Elmelegi, A.; Aly, M.; Elbaksawi, O. Optimum Modified Fractional Order Controller for Future Electric Vehicles and Renewable Energy-Based Interconnected Power Systems. IEEE Access 2021, 9, 29993–30010. [Google Scholar] [CrossRef]
  57. Micev, M.; Ćalasan, M.; Oliva, D. Fractional Order PID Controller Design for an AVR System Using Chaotic Yellow Saddle Goatfish Algorithm. Mathematics 2020, 8, 1182. [Google Scholar] [CrossRef]
  58. Zhao, W.; Zhang, Z.; Wang, L. Manta ray foraging optimization: An effective bio-inspired optimizer for engineering applications. Eng. Appl. Artif. Intell. 2020, 87, 103300. [Google Scholar] [CrossRef]
Figure 1. Structure of studied electrical power grids and the proposed 1 + PD/FOPID control.
Figure 1. Structure of studied electrical power grids and the proposed 1 + PD/FOPID control.
Fractalfract 07 00603 g001
Figure 4. Flowchart representation of the MRFO stages.
Figure 4. Flowchart representation of the MRFO stages.
Fractalfract 07 00603 g004
Figure 5. Proposed design optimization of the 1 + PD/FOPID controller.
Figure 5. Proposed design optimization of the 1 + PD/FOPID controller.
Fractalfract 07 00603 g005
Figure 8. Performance evaluations with an SLP of 2%. Scenario (1): (a Δ f a ; (b Δ f b ; (c Δ P t i e .
Figure 8. Performance evaluations with an SLP of 2%. Scenario (1): (a Δ f a ; (b Δ f b ; (c Δ P t i e .
Fractalfract 07 00603 g008
Figure 9. EV participation in Scenario (1): (a) EV output power; (b) EV battery SOC; (c) generator output power.
Figure 9. EV participation in Scenario (1): (a) EV output power; (b) EV battery SOC; (c) generator output power.
Fractalfract 07 00603 g009
Figure 10. Load profiles at multiple SLPs in Scenario (2).
Figure 10. Load profiles at multiple SLPs in Scenario (2).
Fractalfract 07 00603 g010
Figure 11. Performance evaluations at multiple SLPs Scenario (2): (a Δ f a ; (b Δ f b ; (c Δ P t i e .
Figure 11. Performance evaluations at multiple SLPs Scenario (2): (a Δ f a ; (b Δ f b ; (c Δ P t i e .
Fractalfract 07 00603 g011aFractalfract 07 00603 g011b
Figure 12. EV participation in Scenario (2): (a) EV output power; (b) EV battery SOC; (c) generator output power.
Figure 12. EV participation in Scenario (2): (a) EV output power; (b) EV battery SOC; (c) generator output power.
Fractalfract 07 00603 g012aFractalfract 07 00603 g012b
Figure 13. Load and RES generation profiles in Scenario (3).
Figure 13. Load and RES generation profiles in Scenario (3).
Fractalfract 07 00603 g013
Figure 14. Performance evaluations for RES changes in Scenario (3): (a Δ f a ; (b Δ f b ; (c Δ P t i e .
Figure 14. Performance evaluations for RES changes in Scenario (3): (a Δ f a ; (b Δ f b ; (c Δ P t i e .
Fractalfract 07 00603 g014aFractalfract 07 00603 g014b
Figure 15. Random load change profiles of Scenario (4).
Figure 15. Random load change profiles of Scenario (4).
Fractalfract 07 00603 g015
Figure 16. Performance evaluations at random load changes in Scenario (4): (a Δ f a ; (b Δ f b ; (c Δ P t i e .
Figure 16. Performance evaluations at random load changes in Scenario (4): (a Δ f a ; (b Δ f b ; (c Δ P t i e .
Fractalfract 07 00603 g016
Figure 17. Load profiles at high RES participation in Scenario (5).
Figure 17. Load profiles at high RES participation in Scenario (5).
Fractalfract 07 00603 g017
Figure 18. PV and wind generation profiles at high RES participation in Scenario (5).
Figure 18. PV and wind generation profiles at high RES participation in Scenario (5).
Fractalfract 07 00603 g018
Figure 19. Performance evaluations at high RES participation in Scenario (5): (a Δ f a ; (b Δ f b ; (c Δ P t i e .
Figure 19. Performance evaluations at high RES participation in Scenario (5): (a Δ f a ; (b Δ f b ; (c Δ P t i e .
Fractalfract 07 00603 g019
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

El-Sousy, F.F.M.; Alqahtani, M.H.; Aljumah, A.S.; Aly, M.; Almutairi, S.Z.; Mohamed, E.A. Design Optimization of Improved Fractional-Order Cascaded Frequency Controllers for Electric Vehicles and Electrical Power Grids Utilizing Renewable Energy Sources. Fractal Fract. 2023, 7, 603. https://doi.org/10.3390/fractalfract7080603

AMA Style

El-Sousy FFM, Alqahtani MH, Aljumah AS, Aly M, Almutairi SZ, Mohamed EA. Design Optimization of Improved Fractional-Order Cascaded Frequency Controllers for Electric Vehicles and Electrical Power Grids Utilizing Renewable Energy Sources. Fractal and Fractional. 2023; 7(8):603. https://doi.org/10.3390/fractalfract7080603

Chicago/Turabian Style

El-Sousy, Fayez F. M., Mohammed H. Alqahtani, Ali S. Aljumah, Mokhtar Aly, Sulaiman Z. Almutairi, and Emad A. Mohamed. 2023. "Design Optimization of Improved Fractional-Order Cascaded Frequency Controllers for Electric Vehicles and Electrical Power Grids Utilizing Renewable Energy Sources" Fractal and Fractional 7, no. 8: 603. https://doi.org/10.3390/fractalfract7080603

APA Style

El-Sousy, F. F. M., Alqahtani, M. H., Aljumah, A. S., Aly, M., Almutairi, S. Z., & Mohamed, E. A. (2023). Design Optimization of Improved Fractional-Order Cascaded Frequency Controllers for Electric Vehicles and Electrical Power Grids Utilizing Renewable Energy Sources. Fractal and Fractional, 7(8), 603. https://doi.org/10.3390/fractalfract7080603

Article Metrics

Back to TopTop