Next Article in Journal
Methods for Corrosion Detection in Pipes Using Thermography: A Case Study on Synthetic Datasets
Previous Article in Journal
Theoretical Study of the Wear of a Reduced-Diameter Wheel for Freight Wagons, Based on Its Diameter
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solving Optimal Power Flow Using New Efficient Hybrid Jellyfish Search and Moth Flame Optimization Algorithms

by
Chiva Mayouf
1,2,
Ahmed Salhi
2,
Fanta Haidara
1,
Fatima Zahra Aroua
2,
Ragab A. El-Sehiemy
3,*,
Djemai Naimi
2,
Chouaib Aya
1 and
Cheikh Sidi Ethmane Kane
1,*
1
Unité de Recherche des Nouvelles Technologies des Energies et Systèmes Thermo-Fluides (nTEST), Département de Physique, Faculté des Sciences et Technologie, Université de Nouakchott, Nouakchott BP 888, Mauritania
2
Electrical Engineering Départment, LGEB Laboratory, University of Biskra, BP 145, Biskra 07000, Algeria
3
Faculty of Engineering, Kafrelsheikh University, Kafrelsheikh 33516, Egypt
*
Authors to whom correspondence should be addressed.
Algorithms 2024, 17(10), 438; https://doi.org/10.3390/a17100438
Submission received: 10 August 2024 / Revised: 22 September 2024 / Accepted: 26 September 2024 / Published: 1 October 2024

Abstract

:
This paper presents a new optimization technique based on the hybridization of two meta-heuristic methods, Jellyfish Search (JS) and Moth Flame Optimizer (MFO), to solve the Optimal Power Flow (OPF) problem. The JS algorithm offers good exploration capacity but lacks performance in its exploitation mechanism. To improve its efficiency, we combined it with the Moth Flame Optimizer, which has proven its ability to exploit good solutions in the search area. This hybrid algorithm combines the advantages of both algorithms. The performance and precision of the hybrid optimization approach (JS-MFO) were investigated by minimizing well-known mathematical benchmark functions and by solving the complex OPF problem. The OPF problem was solved by optimizing non-convex objective functions such as total fuel cost, total active transmission losses, total gas emission, total voltage deviation, and the voltage stability index. Two test systems, the IEEE 30-bus network and the Mauritanian RIM 27-bus transmission network, were considered for implementing the JS-MFO approach. Experimental tests of the JS, MFO, and JS-MFO algorithms on eight well-known benchmark functions, the IEEE 30-bus, and the Mauritanian RIM 27-bus system were conducted. For the IEEE 30-bus test system, the proposed hybrid approach provides a percent cost saving of 11.4028%, a percent gas emission reduction of 14.38%, and a percent loss saving of 50.60% with respect to the base case. For the RIM 27-bus system, JS-MFO achieved a loss percent saving of 50.67% and percent voltage reduction of 62.44% with reference to the base case. The simulation results using JS-MFO and obtained with the MATLAB 2009b software were compared with those of JS, MFO, and other well-known meta-heuristics cited in the literature. The comparison report proves the superiority of the JS-MFO method over JS, MFO, and other competing meta-heuristics in solving difficult OPF problems.

1. Introduction

Today, planners and operators in the field of electrical engineering are interested in solving the Optimal Power Flow (OPF) problem as an important tool to determine the optimal operating states of electrical power systems. This tool plays one of the most important roles in modern power system operation and has received a lot of interest over recent decades [1,2]. Furthermore, the problem of economic power generation can be evoked by considering the security of the power system through simultaneously minimizing the total fuel cost and the load bus voltage deviation as a mono-objective optimization problem [3].
In addition to economic and security issues, environmental safety is a major concern that has attracted the attention of many researchers in the OPF field, where the objective function to be minimized is given by the total quantity of pollutant gasses emitted from thermal generating units based on fossil fuel burning. In this case, the environmental dispatch problem is evoked by considering the complex function of gas emission [4].
One of the OPF problems cited previously requires the optimal adjustment of power system control variables, including active power from generating units, the voltage magnitude of the generating buses, the tap ratio of regulating transformers, and the generated reactive powers from VAR compensators. The active power and voltage magnitude of generating buses are continuous variables, while the tap settings of transformers and reactive power from VAR compensators are discrete control variables, which yield to a Mixed-Integer Nonlinear Programming (MINLP) optimization problem. Various classical optimization methods have been applied in the literature to solve OPF problems, such as the gradient-based method [5], linear programming, the Newton method, and the interior point method [6], among others. These methods based on gradient concepts present some drawbacks in dealing with MINLP and highly dimensional OPF problems, particularly their convergence to a local minimum failing in achieving the global optimal solution with highly nonlinear constraints and discontinuous objective functions. To overcome these drawbacks, meta-heuristic optimization methods based on intelligent soft computing and evolutionary algorithms have been developed, especially in recent decades, such as Genetic Algorithms (GAs) [7], Particle Swarm Optimization (PSO) [8], Ant Colony Optimization (ACO) [9], Artificial Bee Colony (ABC) [10], Gravitational Search Algorithm (GSA) [11], Grey Wolf Optimizer (GWO) [12], Dragonfly Algorithm (DA) [13], Teaching–Learning-based Optimization Algorithm (TLBO) [14], Firefly Algorithm (FFA) [15], Harris Hawks Optimization [16], Turbulent Flow of Water-Based Optimization (TFWO) [17], and so on. They have been successfully applied to solve MINLP-OPF problems without considering the derivability and continuity of the objective functions and constraints.
The two main phases of a meta-heuristic algorithm that determine its performance in achieving a global optimal solution are exploration and exploitation [18]. In the past decade, various strategies have been developed to enhance the effectiveness of meta-heuristic techniques by balancing these two phases, particularly when the technique shows weaknesses in one of them. The first strategy focuses on improving the basic meta-heuristic optimization algorithm by adjusting the parameters of the search equation that guide exploration, exploitation, or both. In [19], a Modified Sine Cosine Algorithm (MSCA) is introduced to address the issue of premature convergence in the basic SCA, with two key modifications highlighted. The first modification involves replacing the existing linear update mechanism in the standard SCA with a new nonlinear update mechanism. The second modification adjusts the updated design variable of each agent to include an average design variable between the current design variable of each agent and the best current design variable. In [20], the author proposes an improved Whale Optimization Algorithm (WOA) based on a Nonlinear Adaptive Weight and Golden Sine Operator (NGS-WOA). First, NGS-WOA introduces a nonlinear adaptive weight, allowing search agents to adaptively explore the search space and balance the exploration and exploitation phases. Second, the improved Golden Sine Operator is integrated into the WOA. The second strategy for improving the performance of meta-heuristics involves combining or hybridizing two meta-heuristic techniques to leverage their strengths in managing both the exploration and exploitation phases [18].
To improve the solution quality and speed up the convergence characteristics of meta-heuristics-based OPF problems, various strategies using hybrid meta-heuristic techniques have been proposed in the literature. In [21], a hybrid PSO and Pattern Search (PS) algorithm named PSO-PS is considered to improve the OPF solution by merging the advantages of both basic algorithms with the inclusion of FACTS devices. In [22], a combination of Differential Evolution and Harmony Search algorithms (DE-HS) is elaborated to improve the OPF solution by considering various scales of test systems. In [23], the authors proposed a hybrid meta-heuristic optimization approach based on the Dragonfly Algorithm (DA) and PSO to enhance the search capability to achieve the best global solution quality for solving the single- and multi-objective OPF problems. As reported in [24], a new hybrid approach using the Firefly and PSO algorithms (HFPSO) has been applied to the OPF solution to enhance the exploration and exploitation mechanisms with fast convergence to deal with various complexities of the OPF problem. The recently published work in [25] presents a novel hybrid optimization technique combining ABC and Salp Swarm Algorithms (SSAs) named ABC-SSA to solve the Optimal Reactive Power Dispatch (ORPD) as a sub-problem of OPF. The efficiency of ABC-SSA in balancing between exploration and exploitation mechanisms has been proven for the large-scale test system IEEE 30-bus. The authors of [26] developed a new hybrid fruit fly-based artificial bee colony (HFABC) algorithm to solve single- and multi-objective OPF problems using a fuzzy decision-based mechanism strategy to obtain the best compromise solution from the Pareto front. A Firefly and JAYA Algorithm based on the hybrid evolutionary technique (HFAJAYA) is suggested in [27] to solve the OPF problem by efficiently considering renewable energy sources. Various scenarios in the standard IEEE 30-bus test system network have been implemented to validate the simulation results of the HFAJAYA.
The Jellyfish Search optimizer (JS) algorithm is a modern meta-heuristic optimization technique recently developed by Jui-Sheng Chou and Dinh-Nhat Truong [28] in 2020; it is inspired by the behavior of jellyfish when looking for food in the ocean. The JSO reaffirms its effectiveness by using a set of unimodal and multimodal benchmark functions and various engineering problems. In [29], the JS optimizer algorithm was applied to solve the economic load dispatch (ELD) in the electric power system.
In order to obtain the best solution quality for the OPF problem, this paper proposes a new hybrid approach based on the JS and Moth Flame Optimization (MFO) algorithms. The MFO algorithm was developed by Mirjalili in 2015 [30] and is inspired by the flying behavior of moths with respect to the moon, as they can become trapped in spiral paths around artificial lights. Later, a spiral function was proposed to enhance the MFO algorithm’s performance, as in [31]. The JS optimizer shows an adequate search in the exploration stage, but it suffers during the exploitation phase because it lacks a proficient exploitation operator of the best solutions in the search space. Inversely, the MFO algorithm, by imitating the flying spiral path performed by moths around flame lights, can be useful in the exploitation phase by exploiting the best solutions; nevertheless, it can become trapped in local optima. The main contribution issues can be summarized as follows:
  • Providing a hybrid JS-MFO algorithm to associate the advantages of both the JS optimizer in exploration and MFO in exploitation.
  • The effectiveness of the proposed hybrid approach has been investigated by solving various cases related to the OPF problem.
  • Applications were carried out on two test systems, the IEEE 30-bus and the Mauritanian RIM 27-bus electric power systems, to verify the JS-MFO’s performance.
This is an assessment study based on comparing our simulation results with those of other meta-heuristic methods reported in the literature.
The rest of the structure of this work is ordered as follows: Section 2 presents the formulation of the OPF problem. The mathematical representation of the proposed optimization problem is given in Section 3. Section 4 provides the simulation results and discussions. Finally, a brief conclusion is given in Section 5.

2. Optimal Power Flow Problem

A. Formulation of OPF problem
The OPF problem aims to determine the optimal adjustments of the control equipment by assessing their optimal control variables in the power system network to obtain an optimized objective function, while simultaneously guaranteeing various operating restrictions. These restrictions are given by power flow equations as equality constraints and the functional limits of control equipment are given as inequality constraints. The most recognized objective function to be minimized in the OPF area is that of the total generation cost, with various complexities, including the quadratic, non-quadratic, and piecewise quadratic of the fuel cost function besides the valve point effects, affecting the fuel cost curves of the generating units [2,5,6]. In general, the conventional OPF problem can be mathematically formulated as follows:
min f x , u s u b j e c t   t o   : g x , u = 0 h x , u 0
where f is the objective function to be optimized, g is the set of equality constraints represented by the nonlinear power flow equations, h is the set of inequality constraints reflecting the operating limits of control equipment, u is the vector of control variables, and x is the vector of state variables. The vector u can be expressed as
u = P G 2 . P G N g ,   V G 1 . V G N g ,   T 1 . T N t ,   Q C 1 . Q C N c T
Vector x is given as follows:
x = P G 1 ,   V L 1 . V L N p q ,     Q G 1 . Q G N g ,   S L 1 . S L N l T
The equality constraints represent the balance equations of the active and reactive powers assigned to each bus, and the expression of g(x,u) can be formulated in terms of control and state variables as follows:
P G i P D i V i j = 1 N B V j G i j c o s δ i j + B i j s i n δ i j = 0 Q G i Q D i V i j = 1 N B V j G i j s i n δ i j B i j c o s δ i j = 0
The inequality constraints expressed as h(x,u) are as follows:
Limits of Generators:
The capability constraints for each generation unit are presented for limiting the active and reactive power outputs within the operational limits as follows:
P G i , m i n P G i P G i , m a x Q G i , m i n Q G i Q G i , m a x         i = 1 , .... N g
Similar to the power capability constraints, the voltage at the generator terminals is kept within the minimum and maximum limits as follows:
V G i , m i n V G i V G i , m a x     i = 1 , .... N g
Limits of tap setting transformers:
In Equation (7), the transformer tap settings are kept between the accepted boundaries that are able to control the voltage at the transformer terminals as follows:
T i , m i n T i T i , m a x     i = 1 , .... N t
Limits of reactive power compensators:
Reactive power compensators are bounded in Equation (8) as follows:
Q C i , m i n Q C i Q C i , m a x     i = 1 , .... N C
Limits of voltage magnitude for load buses
V L i , m i n V L i V L i , m a x     i = 1 , .... N P Q
Power flow limits of transmission lines:
S L i , m i n S L i S L i , m a x     i = 1 , .... N l
To handle an inequality constraint of a state variable, a penalty function is introduced in the augmented objective function as follows [6]:
F a u g = f x , u + k = 1 N s λ p k ( x k x k l i m ) 2
where xk is the variable of the k-th violation of the state, xlimk is the limit of the variable of the k-th violation of the state, f(x,u) is the objective function, and λPk is the penalty factor for the penalty function of the variable of the k-th violation of the state.
B. Objective functions:
The mono-objective optimization problem is considered for eight cases of objective function f(x,u) in Equation (1):
Case 1: minimization of total fuel cost
The total cost of electrical power generation is expressed as an objective function (quadratic function) [9] by:
f = F C = i = 1 N g F i P G i = i = 1 N g ( a i + b i P G i + c i P G i 2 )
where ai, bi, and ci are the cost coefficients of the i-th generating unit.
Case 2: minimizing total gas emissions
The gasses emitted by each generating unit can be expressed as a combination of the quadratic and exponential functions of the generated active power [10], and the total amount of gas emissions produced by all generating units:
f = F E m = i = 1 N g ( α i + β i P G i + y i P G i 2 + δ i e x p ε i P G i )     ton / h
where FEm is the total gas emissions (ton/h) and αi, βi, γi, δi, and εi are the emission coefficients of the i-th unit.
Case 3: minimization of total fuel cost considering the valve point effect
A sine component is added to the expression of the objective function in (12) to consider the effect of the valve point to evaluate the total fuel cost as follows [11]:
f = F C = i = 1 N g F i P G i = i = 1 N g ( a i + b i P G i + C i P G i 2 + f v p i )
where fvpi = |di · sin [ei(PGi,minPGi)]|, while ei and di are the fuel cost coefficients of the since component.
Case 4: minimization of total active losses
The total active transmission losses for the power system can be expressed as an objective function:
f = F P l o s s e s = K = 1 N l g k V i 2 + V j 2 2 V i V j cos ( δ i δ j )
where gk is the conductance of the k-th transmission line between buses i and j.
Case 5: minimization of the total load bus voltage deviation
This optimization case is devoted to the minimization of total voltage magnitude deviation related to load buses with the goal of enhancing power system security; the total voltage deviation function for all load buses is given by:
f x , u = F V D = k = 1 N P Q V L k 1
Case 6: improvement of voltage stability index
The stability of a power system is its capacity to maintain all bus voltages within normal limits. The power system reaches a voltage instability state due to severe perturbations, rising electric load, or network configuration changes, which occasions a progressive and unexpected voltage drop. The improvement of the system voltage stability for each bus is a fundamental measure to ensure a safe and stable power system using voltage stability (L-index). Its value varies between 0 and 1 (0 corresponds to a zero-load situation, while 1 means a voltage drop). If the power system has N P Q load buses (PQ buses) and NPV generator buses (PV buses), the L-index is evaluated by [32,33]:
L j = 1 i = 1 N P V F j i V G i V L j           j = 1 ,   2 , ,   N P Q
where VGi is the complex voltage of the generation bus i and VLj is the complex voltage of the load bus j.
F j i = Y L L 1 · Y L G
where the sub-matrices Y L L and Y L G have been obtained from the YBUS matrix after rearranging nodal current injections with respect to nodal voltages, as defined as:
I L I G = Y L L Y L G Y G L Y G G V L V G
The maximum value Lmax of Lj among all load buses is a system stability indicator that can be considered as an objective function presented in the following expression:
F L m a x = max ( L j )       w h e r e           j = 1 , 2 , , N P Q
Case 7: quadratic total fuel cost function and load bus voltage deviation
Minimizing both the total cost of fuel in (6) and the total voltage deviation in (14) to achieve the economical and secure state of the power system, these two functions are introduced in one objective function, as stated in the following expression:
f = F C V D = i = 1 N g a i + b i P G i + c i P G i 2 + Υ V D × F T V D
where   Υ V D is a weighting factor, which is selected as 100.
Case 8: quadratic total fuel cost function and voltage stability index
In order to reduce the total fuel cost and improve the voltage stability of the system simultaneously, the two objective functions in (6) and (13) are implemented into a single-objective function, formulated as follows:
f = F L m a x = i = 1 N g a i + b i P G i + c i P G i 2 + Υ L × F L m a x
The selected value of the weighting factor Υ L is 6000.

3. Mathematical Representation of the Optimization Technique

A. Jellyfish Optimization Algorithm:
The main inspiration of the JS algorithm is the search behavior of jellyfish for food sources in the ocean, considering their intelligent movement behavior when detecting the best location of aliments [28]. Three rules are adopted in JS modeling:
  • Two types of jellyfish movement can be handled by the time control mechanism: moving in the jellyfish swarm or following the ocean current.
  • When searching for food sources in the ocean, jellyfish move in an attractive manner towards the best locations of nutrients.
  • The quantity of food found by jellyfish and its location correspond to the objective function of the problem and the corresponding solution, respectively.
(a) Ocean current
The ocean current includes a considerable amount of food sources, and jellyfish are drawn to its flow. In such a case, the jellyfish update their positions based on the following equation:
X i t + 1 = X i t + r a n d 0 , 1 μ X * β × r a n d 0 , 1 ,   i = 1 ,   2 ,   3 , N pop
where X * : the current best global position of jellyfish in the swarm, r a n d 0 , 1   :   the uniform random number in the interval [0, 1], μ :   the mean location of all jellyfish, and β : the distribution coefficient.
The determination of β and the ocean current direction are given in [28], with β = 3 referring to numerical tests. X i t and   X i t + 1 are the jellyfish positions at the iterations t and (t + 1), respectively, and N p o p : the number of jellyfish.
(b) Jellyfish swarm
In swarms, jellyfish show passive (type A) and active (type B) motions. In the beginning when the swarm is created, jellyfish exhibit the passive motion type (A). Over time, they increasingly exhibit a type B motion. Type A motion is the motion of jellyfish around their own locations, and the consistent updated location of each jellyfish is given by:
X i t + 1 = X i t + γ × r a n d 0 , 1 × U b L b
In order to simulate type B motion, a jellyfish (j) is selected randomly, other than the current jellyfish (i), and a distance vector is formed between jellyfish (i) and jellyfish (j) to establish the step and the direction of the movement. When the quantity of food at the location of the selected jellyfish (j) exceeds that at the location of the jellyfish (i) of interest, the latter moves toward the former; it moves directly away from it if the quantity of food available to the selected jellyfish (j) is lower than that available to the jellyfish of interest (i). So, each jellyfish moves toward a better direction to find the food in a swarm [28]. This behavior can be simulated by the following equations:
S t e p = X i t + 1 X i
where
S t e p = r a n d   0 ,   1 × D i r e c t i o n
D i r e c t i o n = X j t X i t             i f   f X i   f X j X i t X j t         i f   f X i < f   X j    
Here, f is an objective function of location X, and hence the active motions (type B) are accomplished using the following expression:
X i t + 1 = X i t + S t e p
A control parameter, time control mechanism, is employed to switch between the passive and active types of motion. It can also be used to manage the movement of jellyfish in the direction of the ocean current.
(c) Time control mechanism
At each level of changes in wind or temperature, the jellyfish change their positions to another ocean current and a new jellyfish swarm is created. Inside this swarm, the jellyfish favor the passive motion type A in the beginning of the search phase. Over time, the jellyfish begin to prefer active motion (type B). In order to simulate this behavior, a time control mechanism is introduced by regulating the movement, switching between jellyfish moving in a swarm and following the ocean current. This time control mechanism is modeled as a time control function, which randomly fluctuates over the time between 0 and 1. Equation (29) formulates the time control function as below:
C t = 1 t M a x I t × 2 × rand 0 ,   1 1
where MaxIt is the maximal number of iterations and rand is a random number in the interval [0, 1].
(d) Initialization of population
The initial population of jellyfish is generated randomly. To obtain more diversification of the population, a logistic map is developed, which provides more diverse initial populations than random selection and provides a lower probability of premature convergence [28]. Consequently, the logistic map can be used in the initial population based on the following:
X i + 1 = η × X i × 1 X i   w i t h   0 X o 1
where Xi is the chaotic logistic value of the location of the i-th jellyfish, Xo is used to generate the initial population of jellyfish, and the parameter η is set to 4.
(e) Restriction of boundaries
During the swarm movement, jellyfish can move outside of the boundaries of the search space; in this case, the location is updated by returning to the opposite bound, referring to Equation (31):
X i , d       = X i , d         U b , d + L b d         i f   X i , d > U b , d X i , d       = X i , d         L b , d + U b d         i f   X i , d < L b , d
where X i , d is the location of the i-th jellyfish in the d-th dimension; X i , d is the updated location after checking the boundary constraints; and U b , d and X b , d are the upper and lower bounds of the d-th dimension in search spaces, respectively.
B. Moth Flame Optimization Algorithm
As insects, moths share important characteristics with butterflies. The navigation strategy used by moths is a particularly intriguing aspect of their behavior and deserves closer examination. To cover great distances in a straight line, they fly at a fixed angle relative to the moon. This efficient technique is known as transverse orientation [30]. The distance from the light source substantially affects the efficacy of transverse orientation. When the light source is close, the moth spirals around it. This spiral flight path guides the moth to light [31]. The Moth Flame Optimization method was introduced by Mirjalili [30] based on observed behaviors and mathematical models. This method employs moths as potential solutions and their locations are presented as a decision variable vector, allowing each moth to move freely within the problem’s solution space.
M O = m o 1 , 1 m o 1 , 2 . . . . m o 1 , D 1 m o 1 , D m o 2 , 1 m o 2 , 2 . . . . . . m o 2 , D 1 m o 2 , D . . . m o N P O P 1 , 1 m o N P O P 1 , 2 . . m o N P O P 1 , D 1 m o N P O P 1 , D m o N P O P , 1 m o N P O P , 2 . . m o N P O P , D 1 m o N P O P , D
where NPOP is the initial population size of moths and D denotes the number of decision variables. The flame matrix is another fundamental element of the MFO algorithm. Given that each moth exhibits a tendency to navigate in close proximity to its individual flame, it may be inferred that the sizes of the flame matrix are the same as those of the moth’s matrix.
F l = F l 1 , 1 F l 1 , 2 . . . . F l 1 , D 1 F l 1 , D F l 2 , 1 F l 2 , 2 . . . . . . F l 2 , D 1 F l 2 , D . . . F l N P O P 1 , 1 F l N P O P 1 , 2 . . F l N P O P 1 , D 1 F l N P O P 1 , D F l N P O P , 1 F l N P O P , 2 . . F l N P O P , D 1 F l N P O P , D
The distinction between a moth and a flame is based on their respective roles within a problem-solving context. The moth actively engages in a search process, using its flight capabilities to explore various options in order to identify more optimal solutions. In contrast, the flame represents the most optimal solution that has been discovered so far by the moth. To illustrate the helical flight trajectory of moths in proximity to a certain light source, a logarithmic spiral function is used [31].
  MO i x + 1 = MO i x F l i · e b v · cos 2 π v + F l i
where v represents a uniformly distributed random variable ranging from −1 to 1. This parameter quantifies the proximity of the moth’s subsequent location to the associated flame. To improve the efficiency of exploring the solution space during the first iterations and exploiting the solution space during the final iterations, an adaptive approach is suggested to decrease the parameter t throughout the iterations [31].
a = 1 + t 1 M a x I t
v = a 1 × r a n d + 1
The constant, a, exhibits a linear reduction from −1 to −2 during the iterative process. This well-constructed temporal change ensures both exploration and exploitation in the MFO algorithm.
C. Hybrid Jellyfish Search–Moth Flame Optimizer (JS-MFO)
In the artificial JS optimizer, the movement toward an ocean current represents the exploration phase, while the movement within a jellyfish swarm depicts the exploitation phase. A time control mechanism alternates between these two processes [28]. For passive motion (type A) within the jellyfish swarm, each jellyfish updates its position (new solution) relative to its current location based on (24), exploiting the current solution without focusing on the best solutions found nearby. However, exploiting the best solutions in the vicinity of the current solution can help to discover more promising solutions. For active motion (type B), referring to the initial position of jellyfish (i) (current solution (i)), a jellyfish (j) is randomly selected to determine the step size and direction of movement. This updates the position of jellyfish (i) (solution (i)) using (28). This random selection of jellyfish (j) neglects the focus on the best positions (i.e., best solutions corresponding to the highest food availability), causing the algorithm to explore newly discovered solutions randomly without exploiting the best solutions or the best global solution found so far, potentially limiting its ability to find promising solutions. These two shortcomings of passive (A) and active (B) motions reduce the efficiency of the JS algorithm in the exploitation phase.
On the other hand, and based on the literature review, the MFO algorithm signals a useful exploitation capability by carrying out a good search in the proximity of the best solutions. This performance is achieved by using a spiral path around each exploited solution. However, the MFO algorithm presents a drawback; it is unable to share the best global solution obtained so far between agents during the search stage. These characteristics lead to a powerful exploitation mechanism and weak exploration process of the MFO algorithm.
Achieving a new efficient optimization approach is the reason behind our hybridization of the JS and MFO algorithms, which benefits from both advantages of the two algorithms, as shown in Figure 1. In this work, the exploitation capability of the basic JS algorithm is improved by introducing the spiral movement of MFO around the best solutions during the search process. In the early iterations of JS and just after the swarm is formed, the search space is not sufficiently explored, so in this case we try to exploit the good solutions found so far (in order to obtain promising solutions) in the first half-course of iterations (t varies from 0 to MaxIt/2); the elite solutions (which represent 30% of the sorted population of solutions, from the best to the worst) are selected to be exploited, and hence they may be imposed on the spiral movement using (37) and (38). After the first half-course of iterations, the search space is relatively well explored, so in the second half-course of iterations (t changes from MaxIt/2 to MaxIt) we concentrate on the exploitation of the best global solution found so far, which must be imposed on the spiral movement using (39) and (40).
In such a manner, the elite solutions are exploited to achieve promising solutions in the early iterations, and inversely the best global solution is more exploited in the later iterations to improve the quality of the best solution.
D i = X i t X e l i t e
X i t + 1 = D i . e b t . c o s 2 π t + X e l i t e
D i = X i t X *
X i t + 1 = D i   . e b t . c o s 2 π t + X *

4. Simulation Results and Discussion

4.1. Performance Evaluation of JS-MFO for Benchmark Functions

The numerical efficiency of JS-FO was analyzed on eight popular benchmark functions to extensively evaluate its performance. These functions are divided into four unimodal functions and four multimodal functions, which are listed in Table 1 [28]. The unimodal functions were used to verify the exploitation of algorithms. The multimodal functions were used to check the exploration of algorithms. Experimental tests were performed for JS-MFO, MFO, and JS separately in order to verify the effectiveness of the hybrid approach. This task was accomplished for 30 trials, with the maximum number of iterations (MaxIt) = 500 and the population size equals 50 (NPOP = 50). Statistical indices such as the average (mean), standard deviation (SD), and minimum (best) of 30 trials were determined for each benchmark function and are shown in Table 2. The findings in the previous table distinctly exhibit the capacity of the proposed JS-MFO in overcoming the basic algorithms JS and MFO in terms of solution quality (in reaching the near-global solution with slightly more execution time).

4.2. Implementation of JS-MFO for Solving OPF Problems

Two test systems were used to investigate the aptitude of the proposed JS-MFO algorithm to reach the near-global OPF solution for different cases of mono-objective optimization problems: the IEEE 30-bus, a real part of the USA network, and the RIM 27-bus. The modeling of the tested networks was carried out in the MATLAB 2009b environment using a personal computer running Windows XP Professional, Pentium P-IV CPU 3 GHz processor, and RAM of 1 GB.
A. IEEE 30-bus test system
The IEEE 30-bus model has 6 generators, 41 branches (37 lines and 4 transformers with off-nominal tap ratios), and 24 load buses as shown in Figure 2. Shunt Volt Ampere Reactive (VAR) compensators are installed in buses 10, 12, 15, 17, 20, 21, 23, 24, and 29, where the reactive power injection is controlled between 0 and 5 MVAR as the lower and upper limits, respectively [11]. The total system demand was (2.834 + j1.262) p.u. for the apparent power at the 100 MVA base. Bus 1 was taken as the slack bus. The upper and lower limits of active power generation, the reactive power limits, and the generator cost coefficients are taken from [34].
In the base case, the system presents the following data: Fc (quadratic fuel cost) = 901.935 USD/h, FEm = 0.23908 ton/h, FPlosses = 5.8332 MW, FVD = 1.1663 p.u., and Lmaxr = 0.1731. The parameters of the JS, MFO. and JS-MFO algorithms are as follows: the population number (population size) NPOP is 50, the maximum number of iterations MaxIt = 100, and the number of control variables is D = 25. For the present test system, six cases (six objective functions) were studied to evaluate the performance of the proposed JS-MFO compared to that of JS and MFO in order to reach the optimal solution. In order to validate the proposed JS-MFO in dealing with and solving complex OPF problems, 30 independent trial runs were performed for each case study by implementing JS, MFO, and JS-MFO. Three complex cases, particularly 1, 3, and 4, were considered to provide the minimum, average, maximum, and standard deviation (SD) for each objective function as shown in Table 3. A comparative report of each index between the JS, MFO, and JS-MFO algorithms, shown in the previous table, can clearly prove the superiority of JS-MFO over JS and MFO in providing the best solution quality and stability over the trial runs. The optimal control variable settings and the corresponding objective function value are depicted for each case in Table 4 and Table 5. Executing 30 independent runs for each case led to extracting the solution for the run that provided the best solution.
Case 1: Minimization of the quadratic total fuel cost function
The objective function in Equation (12) was evaluated for the optimal settings of the control variables by running the JS, MFO, and JS-MFO algorithms, and the simulation results are indicated in Table 4 for case 1. The optimal total fuel costs obtained were 799.085 USD/h, 801.02 USD/h, and 800.79 USD/h using JS-MFO, JS, and MFO, respectively. The best solution of the total fuel cost for JS-MFO was compared with that of the JS and MFO algorithms and other optimization techniques in the literature shown in Table 6, such as the Novel Improved Social Spider Optimization algorithm (NISSO) [35], Firefly Algorithm (FFA) [15], Moth Swarm Algorithm (MSA) [15], Efficient Sine Cosine Algorithm (ESCA) [36], Symbiotic Organism Search algorithm (SOS) [35], Harris Hawks Optimization (HHO) [16], Dragonfly Algorithm (DA) [13], Adaptive Gaussian Teaching–Learning-based Optimization (AGLBO) [14], PSO-PS [21], DA-PSO [23], HF-PSO [24], Artificial Bee Colony based on hybrid fruit fly (HF-ABC) [25] and HFA-JAYA [27]. It is very clear from Table 6 that the JS-MFO approach is able to reduce the fuel cost to the lowest value of 799.085 USD/h. The percent cost savings with respect to the base cases of JS, MFO and JS-MFO are 11.188%, 11.2137%, and 11.4028 %, respectively, which lead to annual cost savings for JS-MFO of 16,951 USD/year compared to JS and 14,936 USD/year compared to MFO. The evolution of the total fuel cost during the simulation is indicated in Figure 3. It can be seen that the optimal total fuel cost converges towards the best optimal solution using JS-MFO.
Case 2: Minimization of Total Gas Emissions
The total gas emission function in (13) was minimized using the proposed JS-MFO, MFO, and JS algorithms. The gas emission coefficients are given in [34]. The optimal solutions achieved using JS-MFO, MFO, and JS are 0.2047, 0.2048, and 0.2049 ton/h, respectively, based on the optimal setting of the control variables shown in Table 4 for case 2. Table 7 gives a statistical comparison report considering other optimization methods applied to the same test system and the same number of control variables as the Harris Hawks optimization algorithm (HHO) [16], Turbulent Flow Of Water-Based Optimization algorithm (TFWO) [17], hybrid DA-PSO [23], Adaptive Gaussian AGTLBO [14], and Social Spider Optimization algorithm (SSO) [37]. By examining the previous table, the minimum total gas emissions obtained by the proposed JS-MFO algorithm was shown to be better than that of the JS and MFO algorithms and the comparison optimization techniques reported in the literature. The percent gas emission reductions were 14.30%, 14.34%, and 14.38% for JS, MFO, and JS-MFO, respectively, referring to the base case. It can be seen that the proposed hybrid technique JS-MFO achieves total gas emission reductions of 1.752 tons/year and 0.876 tons/year compared to JS and MFO, respectively. The convergence of the total gas emission function with respect to the number of iterations is depicted in Figure 4, showing that the proposed JS-MFO algorithm gives faster convergence to the optimal solution with a better solution quality than the JS and MFO algorithms.
Case 3: Minimization of total fuel cost considering valve point effect
The objective function given by (14) is provided for the non-smooth curve of the total fuel cost function and that for the two generators in buses 1 and 2, where their cost coefficients are taken from [10]. The fuel cost curves of the remaining generators keep the same characteristics as in case 1. The minimum total fuel costs obtained using the proposed JS-MFO, JS, and MFO algorithms are 918.073 USD/h, 922.48 USD/h, and 924.02 USD/h, respectively, considering the valve point effect. The optimal solution reflects the optimal control variables presented in Table 4 for case 3. The percent cost savings considering the valve point effect for JS, MFO, and JS-MFO are 15.01%, 14.87%, and 15.42%, respectively, compared to the base case. It is clearly shown that JS-MFO gives the highest cost saving. The annual cost savings of JS-MFO correspond to 38,632 USD/year and 52,122 USD/year with respect to those of JS and MFO, respectively. The JS-MFO algorithm shows a better solution than that of the JS and MFO algorithms, and the other methods displayed in Table 8 such as the Improved Salp Swarm Algorithm (ISSA) [38], Salp Swarm Algorithm (SSA) [38], Teaching–Learning-based Optimization Algorithm (TLBO) [39], Gbest-guided ABC (Gbest-ABC) [40], and ABC-based Grenade Explosion Method (GABC) [41]. The progress of the non-quadratic total fuel cost with respect to the iteration evolutions is given in Figure 5 for the JS, MFO, and JS-MFO algorithms. These curves confirm the superiority of JS-MFO over JS and MFO in attaining the best global solution.
Case 4: Minimization of Total Active Losses
Active transmission losses in Equation (15) are adopted as an objective function in this case and are minimized by carrying out the proposed JS-MFO, JS, and MFO techniques. The optimal total active losses resulting from the simulation are given in Table 5 for case 4, along with their corresponding control and state variables. It is reported that JS-MFO is able to reduce the active transmission losses to the lowest value of 2.8810 MW, while JS and MFO lead to 3.232 MW and 3.049 MW, respectively. The percent loss savings of the presented techniques JS, MFO, and JS-MFO are, respectively, 44.59%, 47.72%, and 50.60% with reference to the base case.
For comparison purposes, Table 9 depicts the optimal active transmission losses for JS-MFO, JS, MFO, and other techniques inspired from previous works in the literature such as AGTLBO [14], FFA [15], MSA [15], DE-HS [22], DA [23], DA-PSO [23], HFA-JAYA [27], SS0 [35], NISSO [35], ESCA [36], and Adaptive Multiple-Team Perturbation-Guiding JAYA [42]. From this table, it can be observed that the results obtained using JS-MFO are better than those of the other methods presented in the current literature. Figure 6 shows the evolution of total active losses with respect to the progression of iterations during the simulation, depicting the convergence characteristics of the JS and JS-MFO algorithms for this case.
Case 7: Minimization of the total fuel cost along with the deviation of the total load bus voltage
This case study is devoted to the simultaneous minimization of the quadratic total fuel cost and the total load bus voltage deviation as the objective function depicted in Equation (21). The problem is solved by running the JS, MFO, and JS-MFO algorithms separately. The simulation results for case 7 are listed in Table 5 by showing the numerical optimal settings of the control and state variables. Table 10 is provided to compare the optimal total fuel cost and the total load bus voltage deviation of the JS-MFO method to those of the JS and MFO algorithms and the other techniques given by recently published works in the literature. It is clearly shown that JS-MFO dominates all other comparison techniques in terms of the optimal fuel cost with a slightly higher optimal voltage deviation compared to AGTLBO [14], HFA-JAYA [27], ESCA [36], PSO-SSA [37], MVO [43] and ECHT-DE [44]. The ESCA is reported to present the best optimal voltage deviation among all other techniques, but it provides the highest optimal voltage deviation.
It is observed from the results that JS-MFO provides simultaneous reductions in quadratic total fuel cost and total voltage deviation of 10.90% and 91.48%, respectively, with reference to the base case, while the other existing techniques, JS, and MFO show less important results. By examining the convergence curves of the JS and JS-MFO algorithms given by Figure 7, it can be seen that the JS-MFO algorithm gives the best final global optimal solution compared to the JS and MFO algorithms.
Case 8: Minimization of Quadratic Total Fuel Cost Along With Vulnerability Stability Enhancement
The present case study examines the minimization of the total fuel cost along with the enhancement of voltage stability by optimizing the objective function mentioned in (22). The optimal settings obtained for the control and state variables after the JS, MFO, and JS-MFO algorithm running processes are illustrated in Table 5 for case 8. The optimal solution presented by JS-MFO gives the best simultaneous reductions in total fuel cost and voltage stability index than that of JS and MFO. By referring to the base case, JS-MFO presents, simultaneously, a percent cost saving of 11.37% and an improvement of the stability margin of 34.08%. A comparison in Table 11 with other meta-heuristics-based optimization techniques such as ESCA [36], PSO-SSA [37], AMTPG-JAYA [42], MVO [43], and ECHT-DE [44] shows that the proposed JS-MFO presents the best global solution among others.
B. Mauritanian electric power system RIM 27-bus:
B.1. Simulation of various cases for RIM 27-bus
The RIM 27-bus power grid is controlled by the SOMELEC society, involving three voltage levels: 225 kV, 90 kV, and 33 kV. This system model consists of six generation buses: 1 (slack bus), 8, 10, 23, 24, and 25, feeding a total load of (1.786 + j0.736) p.u. at a base power of 100 MVA, and three shunt compensators at buses 12, 14, and 16. Two solar power plants are installed on buses 8 and 23 with active power capacities of 50 MW and 15 MW, respectively, and a wind power plant on bus 10 with a capacity of 30 MW, as shown in Figure 8. Also, the RIM 27-bus system has 21 load buses and 36 branches. The voltage magnitude of the generation buses is controlled between 0.92 p.u. and 1.08 p.u. The number of control variables for this system is 15. In the base case, the RIM 27-bus system presents the following data: FPlosses = 6.6015 MW, FVD = 1.6358 p.u., and Lmaxr = 0. 522.
Case 4: Minimization of Total Active Losses for RIM 27-bus system:
In this case study, the aim is to minimize the transmission active losses of the RIM 27-bus system, adopting (15) as the objective function, and the simulation results are tabulated in Table 12. The optimal active losses affected by JS, MFO, and JS-MFO are 3.35 MW, 3.56 MW, and 3.25 MW, which represent loss percent savings of 49.19%, 46.04%, and 50.67% compared to the base case, respectively. The best solution quality was found with the JS-MFO technique, which yielded the lowest total active power losses and the highest loss percent saving. The optimal total active transmission losses attributed to JS-MFO are 2.94% and 8.58% lower than those obtained by the JS and MFO algorithms. It is clearly illustrated in Figure 9 that the JS-MFO algorithm dominates the JS and MFO algorithms in terms of solution quality and convergence characteristics.
Case 5: Minimizing of the total voltage deviation for the RIM 27-bus system:
The objective function presented in (16) is minimized in case 5 and depicts an optimal solution for the control and state variables shown in Table 12. The RIM 27-bus power system security is enhanced via the minimization of the total voltage deviation FVD using JS, MFO, and JS-MFO, signaling reductions of 59.97%, 55.25%, and 62.44%, respectively, compared to the base case. The JS-MFO pointed out the best optimal FVD compared to that of JS and MFO. The optimal total voltage deviation achieved by JS-MFO is 0.6144 p.u., which is better than that of the JS and MFO algorithms, which correspond to 0.6547 p.u. and 0.732 p.u., respectively. The best performance of the JS-MFO algorithm is verified in Figure 10, illustrating the convergence proprieties of JS-MFO overcoming those of the JS and MFO algorithms. Case 6: Voltage stability enhancement for the RIM 27-bus system:
The FLmax voltage stability index is optimized for case 6 in order to enhance voltage stability using the JS, MFO, and JS-MFO algorithms. The stability margin of the RIM 27-bus system is improved using JS, MFO, and JS-MFO, where FLmax reaches the optimal values of 0.3861, 0.4201, and 0.3389, respectively. The JS-MFO approach denoted a percent reduction of FLmax of 12.22% compared to JSO and 19.33% compared to MFO, which proves the capacity of JS-MFO to overcome the JS and MFO algorithms.
It appears from Table 12 that the JS-MFO algorithm outperforms the JS and MFOs in solving the OPF problem. The JS-MFO simulation result provides 12.22% and 19.33% improvements in the stability margin compared to the JS and MFO algorithms, respectively.
B.2. Sizing of reactive power compensators for renewable sources
Table 13 points out the reactive power generation for each optimal state of cases 4, 5, and 6. The reactive power compensator sizing for each renewable energy source in generation buses 8, 10, and 23 is, respectively, −17.5/17.5 MVAR, −16/16 MVAR, and −15/15 MVAR. Figure 11 indicates the voltage magnitude profile for each optimal state for cases 4–6, describing that the voltage magnitude is within the predefined limits between 0.92 p.u. and 1.08 p.u.

5. Conclusions

A new meta-heuristic optimization technique has been developed in this paper, based on the hybridization of two methods, Jellyfish Search (JS) and Moth Flame Optimizer (MFO), called JS-MFO. This new hybrid approach is designed to improve the solution quality when dealing with complex Optimal Power Flow (OPF) problems. To improve the performance of the JS method in the exploitation mechanism, a spiral path of the MFO method on the best solutions found after the exploration phase is introduced. In this way, we can benefit from the advantages of the JS and MFO algorithms. The efficiency and precision of JS-MFO have been confirmed by applying it to several complex OPF optimization problems. These problems deal with the minimization of total power generation cost, gas emissions from power plants, total voltage deviation, and total transmission losses. Two test grids were employed to solve these problems: the popular IEEE 30-bus grid and the Mauritanian RIM 27-bus grid. For the IEEE 30-bus test system, JS-MFO presented a percent cost saving of 11.40% (optimization of the quadratic total fuel cost) with respect to the base case. This benefit is the highest when it is compared to those of the JS and MFO algorithms. It results in an annual cost saving of 16,951 USD/year with respect to the optimal value given by the JS algorithm. The percent gas emission reduction when dealing with the environmental optimization problem is about 14.38% for the JS-MFO algorithm (the highest value compared to those of JS and MFO) with reference to the base case. The JS-MFO algorithm shows a gas emission reduction of 1.752 tons/year with respect to the value given by the JS algorithm, and of 0.876 tons/year with respect to the amount given by the MFO algorithm. By optimizing the active transmission losses, a significant percent loss saving of 50.60% is assigned to the JS-MFO algorithm (greater than those of JS and MFO) with reference to the base case. For the RIM 27-bus power system, the loss percent saving using JS-MFO is 50.67% compared to the base case. The optimization of the total voltage deviation using JS-MFO showed a reduction of 62.44% with respect to the base case. The results obtained prove the capacity of the proposed approach to dominate the basic meta-heuristic techniques, JS, MFO, and others recently developed in the literature based on the comparison report for the same test network and the same objective functions. The aptitude of this hybrid JS-MFO technique to deal with complex OPF problems makes it competitive for (i) treating more problems in the power system, such as a variety of cases with complex objective functions, security constraints, prohibited zones, and different test systems, and (ii) offering prospects for improving the performance of the Mauritania RIM 27-bus grid by optimizing the location and sizing of new renewable sources. In view of the potential and superior characteristics of the proposed JS-MFO algorithm, a multi-objective algorithm based on JS-MFO is recommended to be developed and applied to solve multi-objective OPF problems.

Author Contributions

Conceptualization: C.M. and A.S.; methodology: A.S., D.N., and C.S.E.K.; software: C.M., A.S., and F.Z.A.; validation: R.A.E.-S., C.A., and D.N.; formal analysis: C.M., R.A.E.-S., F.H., and F.Z.A.; investigation: A.S., R.A.E.-S., D.N., and C.S.E.K.; resources: C.A. and C.S.E.K.; data curation: C.M., R.A.E.-S., and C.A.; writing—original draft preparation: C.M., F.Z.A., and F.H.; writing—review and editing: C.M., F.Z.A., F.H., and R.A.E.-S.; visualization: F.H., D.N., and C.S.E.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

All required data are included in the text.

Conflicts of Interest

The authors declare no conflicts of interest.

Nomenclature

PG1Generated active power of slack bus
PGiGenerated active power of i-th bus
VLiVoltage magnitude of load bus i
QGiReactive power output of generator i
PGi,minLower active generation capacity in i-th bus
PGi,maxUpper active generation capacity in i-th bus
PDiActive power of load demand on bus i
QDiReactive power of load demand on bus i
NgNumber of generators (power plants)
NPQNumber of load buses (PQ buses)
SLiApparent power in transmission line i
Vi, VjVoltage magnitudes at node i and node j, respectively
SLi,minLower limit of apparent power in transmission line i
SLi,maxUpper limit of apparent power in transmission line i
VLi,minLower limit of voltage magnitude in load bus i
VLi,maxUpper limit of voltage magnitude in load bus i
QGi,minLower reactive generation capacity in i-th bus
QGi,maxUpper reactive generation capacity in i-th bus
VGiVoltage magnitude of generator connected to bus i
QciReactive power injected by shunt VAR compensator in i-th bus
VGi,minMinimum voltage magnitude of generator connected to bus i
VGi,maxMaximum voltage magnitude of generator connected to bus i
δi, δjVoltage angles of buses i and j, respectively
δijDifference angle between voltage angles δi and δj, respectively
BijSusceptance of the transmission line between buses i and j (imaginary part of the admittance)
GijConductance of the transmission line between buses i and j (real part of the admittance)
NtNumber of regulating transformers
NBNumber of buses
NlNumber of the transmission lines
NCNumber of shunt compensators
NPVNumber of generator buses (PV buses)
TiTap setting of i-th transformer
Ti,minLower tap setting of i-th transformer
Ti,maxUpper tap setting of i-th transformer
PG1Generated active power of slack bus
NgNumber of generator buses
VLjVoltage magnitude of load bus j
NPQNumber of load buses (PQ buses)
QGiReactive power output of generator i
QDiActive power of load demand on bus i
SLiApparent power in transmission line i
Vi, VjVoltage magnitudes at node i and node j, respectively
PGiActive power output of generator i
PDiActive power of load demand on bus i
VGiVoltage magnitude of generator connected to bus i
QciReactive power injected by shunt VAR compensator in i-th bus
BijSusceptance of the transmission line between buses i and j (imaginary part of the admittance)
GijConductance of the transmission line between buses i and j (real part of the admittance)
NTNumber of regulating transformers
NBNumber of buses
TiTap setting of transformer i
NCNumber of shunt compensators
NPVNumber of generator buses (PV buses)
NlNumber of the transmission lines
δijDifference angle between voltage anglesδi and δj, respectively
δi,δjVoltage angles of buses i and j, respectively

References

  1. Frank, S.; Steponavice, I.; Rebennack, S. Optimal power flow: A bibliographic survey I. Energy Syst. 2012, 3, 221–258. [Google Scholar] [CrossRef]
  2. Niknam, T.; Azizipanah-Abarghooee, R.; Narimani, M.R. Reserve constrained dynamic optimal power flow subject to valve-point effects, prohibited zones and multi-fuel constraints. Energy 2012, 47, 451–464. [Google Scholar] [CrossRef]
  3. El Ela, A.A.A.; El-Sehiemy, R.A.; Shaheen, A.M.; Shalaby, A.S. Application of the crow search algorithm for economic environmental dispatch. In Proceedings of the 2017 IEEE Nineteenth International Middle East Power Systems Conference (MEPCON), Cairo, Egypt, 19–21 December 2017; pp. 78–83. [Google Scholar]
  4. Chaib, A.E.; Bouchekara, H.; Mehasni, R.; Abido, M.A. Optimal power flow with emission and non-smooth cost functions using backtracking search optimization algorithm. Int. J. Electr. Power Energy Syst. 2016, 81, 64–77. [Google Scholar] [CrossRef]
  5. Dommel, H.W.; Tinney, W.F. Optimal power flow solutions. IEEE Trans. Power Appar. Syst. 1968, 10, 1866–1876. [Google Scholar] [CrossRef]
  6. Momoh, J.A.; El-Hawary, M.E.; Adapa, R. A review of selected optimal power flow literature to 1993. II. Newton, linear programming and interior point methods. IEEE Trans. Power Syst. 1999, 14, 105–111. [Google Scholar] [CrossRef]
  7. Osman, M.S.; Abo-Sinna, M.A.; Mousa, A.A. A solution to the optimal power flow using genetic algorithm. Appl. Math. Comput. 2004, 155, 391–405. [Google Scholar] [CrossRef]
  8. Abido, M.A. Optimal power flow using particle swarm optimization. Int. J. Electr. Power Energy Syst. 2002, 24, 563–571. [Google Scholar] [CrossRef]
  9. Bouktir, T.; Slimani, L. Optimal power flow of the Algerian electrical network using an ant colony optimization method. Leonardo J. Sci. 2005, 6, 43–57. [Google Scholar]
  10. Adaryani, M.R.; Karami, A. Artificial bee colony algorithm for solving multi-objective optimal power flow problem. Int. J. Electr. Power Energy Syst. 2013, 53, 219–230. [Google Scholar] [CrossRef]
  11. Duman, S.; Güvenç, U.; Sönmez, Y.; Yörükeren, N. Optimal power flow using gravitational search algorithm. Energy Convers. Manag. 2012, 59, 86–95. [Google Scholar] [CrossRef]
  12. El-Fergany, A.A.; Hasanien, H.M. Single and multi-objective optimal power flow using grey wolf optimizer and differential evolution algorithms. Electr. Power Compon. Syst. 2015, 43, 1548–1559. [Google Scholar] [CrossRef]
  13. Ouafa, H.; Linda, S.; Tarek, B. Multi-objective optimal power flow considering the fuel cost, emission, voltage deviation and power losses using Multi-Objective Dragonfly algorithm. In Proceedings of the International Conference on Recent Advances in Electrical Systems. In Proceedings of the International Conference on Recent Advances in Electrical Systems, Hammamet, Tunusia, 22–24 December 2017. [Google Scholar]
  14. Alanazi, A.; Alanazi, M.; Memon, Z.A.; Mosavi, A. Determining Optimal Power Flow Solutions Using New Adaptive Gaussian TLBO Method. Appl. Sci. 2022, 12, 7959. [Google Scholar] [CrossRef]
  15. Vijaya Bhaskar, K.; Ramesh, S.; Abudhahir, A. Modern swarm intelligence-based algorithms for solving optimal power flow problem in a regulated power system framework. Turk. J. Comput. Math. Educ. (TURCOMAT) 2021, 12, 1786–1793. [Google Scholar] [CrossRef]
  16. Islam, M.Z.; Wahab, N.I.; Veerasamy, V.; Hizam, H.; Mailah, N.F.; Guerrero, J.M.; Nasir, M.N.M. A Harris Hawks optimization based single-and multi-objective optimal power flow considering environmental emission. Sustainability 2020, 12, 5248. [Google Scholar] [CrossRef]
  17. Sarhan, S.; El-Sehiemy, R.; Abaza, A.; Gafar, M. Turbulent flow of water-based optimization for solving multi-objective technical and economic aspects of optimal power flow problems. Mathematics 2022, 10, 2106. [Google Scholar] [CrossRef]
  18. Mirjalili, S. (Ed.) Handbook of Whale Optimization Algorithm: Variants, Hybrids, Improvements, and Applications; Academic Press Elsevier: Amsterdam, The Netherlands, 2023. [Google Scholar]
  19. Shaheen, A.M.; Elattar, E.E.; El-Sehiemy, R.A.; Elsayed, A.M. An improved sunflower optimization algorithm-based Monte Carlo simulation for efficiency improvement of radial distribution systems considering wind power uncertainty. IEEE Access 2020, 9, 2332–2344. [Google Scholar] [CrossRef]
  20. Zhang, J.; Wang, J.S. Improved whale optimization algorithm based on nonlinear adaptive weight and golden sine operator. IEEE Access 2020, 8, 77013–77048. [Google Scholar] [CrossRef]
  21. Berrouk, F.; Bounaya, K. Optimal power flow for multi-FACTS power system using hybrid PSO-PS algorithms. J. Control. Autom. Electr. Syst. 2018, 29, 177–191. [Google Scholar] [CrossRef]
  22. Reddy, S.S. Optimal power flow using hybrid differential evolution and harmony search algorithm. Int. J. Mach. Learn. Cybern. 2019, 10, 1077–1091. [Google Scholar] [CrossRef]
  23. Khunkitti, S.; Siritaratiwat, A.; Premrudeepreechacharn, S.; Chatthaworn, R.; Watson, N.R. A hybrid DA-PSO optimization algorithm for multi-objective optimal power flow problems. Energies 2018, 11, 2270. [Google Scholar] [CrossRef]
  24. Khan, A.; Hizam, H.; Wahab, N.I.B.A.; Othman, M.L. Optimal power flow using hybrid firefly and particle swarm optimization algorithm. PLoS ONE 2020, 15, e0235668. [Google Scholar] [CrossRef] [PubMed]
  25. Salhi, S.; Naimi, D.; Salhi, A.; Abujarad, S.; Necira, A. A novel hybrid approach based artificial bee colony and salp swarm algorithms for solving ORPD problem. Indones. J. Electr. Eng. Comput. Sci. 2021, 23, 1825–1837. [Google Scholar] [CrossRef]
  26. Mallala, B.; Papana, V.P.; Sangu, R.; Palle, K.; Chinthalacheruvu, V.K. multi-objective optimal power flow solution using a non-dominated sorting hybrid fruit fly-based artificial bee colony. Energies 2022, 15, 4063. [Google Scholar] [CrossRef]
  27. Alghamdi, A.S. A Hybrid Firefly–JAYA Algorithm for the Optimal Power Flow Problem Considering Wind and Solar Power Generations. Appl. Sci. 2022, 12, 7193. [Google Scholar] [CrossRef]
  28. Chou, J.-S.; Truong, D.-N. A novel metaheuristic optimizer inspired by behavior of jellyfish in ocean. Appl. Math. Comput. 2021, 389, 125535. [Google Scholar] [CrossRef]
  29. Rai, H.; Verma, H.K. Economic load dispatch using jellyfish search optimizer. In Proceedings of the 2021 10th IEEE International Conference on Communication Systems and Network Technologies (CSNT), Bhopal, India, 18–19 June 2021. [Google Scholar] [CrossRef]
  30. Mirjalili, S. Moth-flame optimization algorithm: A novel nature-inspired heuristic paradigm. Knowl.-Based Syst. 2015, 89, 228–249. [Google Scholar] [CrossRef]
  31. Khalilpourazari, S.; Pasandideh, S.H. Multi-item EOQ model with nonlinear unit holding cost and partial backordering: Moth-flame optimization algorithm. J. Ind. Prod. Eng. 2016, 34, 42–51. [Google Scholar] [CrossRef]
  32. Kessel, P.; Glavitsh, H. Estimating the voltage stability of power systems. IEEE Trans. Power Syst. 1986, 1, 346–354. [Google Scholar] [CrossRef]
  33. Shaheen, A.M.; El-Sehiemy, R.A.; Farrag, S.M. A novel adequate bi-level reactive power planning strategy. Int. J. Electr. Power Energy Syst. 2016, 78, 897–909. [Google Scholar] [CrossRef]
  34. Bouktir, T.; Slimani, L. Object-Oriented Economic Power Dispatch of Electrical Power System with minimum pollution using a Genetic Algorithm. J. Electr. Syst. 2005, 1, 19–34. [Google Scholar]
  35. Almabsout, E.A.; El-Sehiemy, R.A.; An, O.N.U.; Bayat, O. A Hybrid Local Search-Genetic Algorithm for Simultaneous Placement of DG Units and Shunt Capacitors in Radial Distribution Systems. IEEE Access 2020, 8, 54465–54481. [Google Scholar] [CrossRef]
  36. Messaoudi, A.; Belkacemi, M. Optimal Power Flow Solution using Efficient Sine Cosine Optimization Algorithm. J. Intell. Syst. Appl. 2020, 2, 34–43. [Google Scholar] [CrossRef]
  37. El Sehiemy, R.A.; Selim, F.; Bentouati, B.; Abido, M.A. A novel multi-objective hybrid particle swarm and salp optimization algorithm for technical-economical-environmental operation in power systems. Energy 2020, 193, 116817. [Google Scholar] [CrossRef]
  38. Abd el-sattar, S.; Kamel, S.; Ebeed, M.; Jurado, F. An improved version of salp swarm algorithm for solving optimal power flow problem. Soft Comput. 2021, 25, 4027–4052. [Google Scholar] [CrossRef]
  39. Xing, H.; Song, F.; Yan, L.; Pan, W. A modified artificial bee colony algorithm for load balancing in network-coding-based multicast. Soft Comput. 2019, 23, 6287–6305. [Google Scholar] [CrossRef]
  40. Roy, R.; Jadhav, H.T. Optimal power flow solution of power system incorporating stochastic wind power using Gbest guided artificial bee colony algorithm. Int. J. Electr. Power Energy Syst. 2015, 64, 562–578. [Google Scholar] [CrossRef]
  41. Salhi, A.; Naimi, D.; Bouktir, T. Optimal power flow resolution using artificial bee colony algorithm-based grenade explosion method. J. Electr. Syst. 2016, 12, 734–756. [Google Scholar]
  42. Warid, W. Optimal power flow using the AMTPG-Jaya algorithm. Appl. Soft Comput. 2020, 91, 106252. [Google Scholar] [CrossRef]
  43. Bentouati, B.; Chettih, S.; Jangir, P.; Trivedi, I.N. A solution to the optimal power flow using multi-verse optimizer. Electr. Syst. 2016, 12, 716–733. [Google Scholar]
  44. Biswas, P.P.; Suganthan, P.N.; Mallipeddi, R.; Amaratunga, G.A. Optimal power flow solutions using differential evolution algorithm integrated with effective constraint handling techniques. Eng. Appl. Artif. Intell. 2018, 68, 81–100. [Google Scholar] [CrossRef]
Figure 1. Flowchart of the JS-MFO algorithm.
Figure 1. Flowchart of the JS-MFO algorithm.
Algorithms 17 00438 g001
Figure 2. IEEE 30-bus single-line diagram.
Figure 2. IEEE 30-bus single-line diagram.
Algorithms 17 00438 g002
Figure 3. Convergence curves of JS, MFO and JS-MFO algorithms for case 1.
Figure 3. Convergence curves of JS, MFO and JS-MFO algorithms for case 1.
Algorithms 17 00438 g003
Figure 4. Convergence rates of JS, MFO, and JS-MFO algorithms for case 2.
Figure 4. Convergence rates of JS, MFO, and JS-MFO algorithms for case 2.
Algorithms 17 00438 g004
Figure 5. Convergence rates of JS, MFO, and JS-MFO algorithms for case 3.
Figure 5. Convergence rates of JS, MFO, and JS-MFO algorithms for case 3.
Algorithms 17 00438 g005
Figure 6. Convergence rates of JS, MFO, and JS-MFO for Case 4.
Figure 6. Convergence rates of JS, MFO, and JS-MFO for Case 4.
Algorithms 17 00438 g006
Figure 7. Convergence curves of JS, MFO, and JS-MFO algorithms for case 7.
Figure 7. Convergence curves of JS, MFO, and JS-MFO algorithms for case 7.
Algorithms 17 00438 g007
Figure 8. Single-line diagram of the Mauritanian electric power system 27 bus.
Figure 8. Single-line diagram of the Mauritanian electric power system 27 bus.
Algorithms 17 00438 g008
Figure 9. Convergence characteristics of JS, MFO, and JS-MFO algorithms for case 4.
Figure 9. Convergence characteristics of JS, MFO, and JS-MFO algorithms for case 4.
Algorithms 17 00438 g009
Figure 10. Convergence characteristics of JS, MFO, and JS-MFO algorithms for case 5.
Figure 10. Convergence characteristics of JS, MFO, and JS-MFO algorithms for case 5.
Algorithms 17 00438 g010
Figure 11. Voltage profile for each RIM 27-bus case using the JS-MFO algorithm.
Figure 11. Voltage profile for each RIM 27-bus case using the JS-MFO algorithm.
Algorithms 17 00438 g011
Table 1. Benchmark function characteristics.
Table 1. Benchmark function characteristics.
FunEquationBMinDT
f 1 f x = i = 1 n x i 2 [−100, 100]030Unimodal
f 2 f x = i = 1 n i x i 2 [−10, 10]030Unimodal
f 3 f x = i = 1 n 1 100 x i + 1 x i 2 2 + x i 1 2 [−30, 30]030Unimodal
f 4 f x = i = 1 n x i + i = 1 n x i [−10, 10]030Unimodal
f 5 f x = 20 e x p 0.2 1 n i = 1 n x i 2 e x p 1 n i = 1 n cos 2 π x i + 20 + e [−32, 32]030Multimodal
f 6 f x = i = 1 n x i 2 10 cos 2 π x i + 10 2 [−5.12, 5.12]030Multimodal
f 7 f x = π n 10 s i n 2 π y 1 + i = 1 n 1 y i 1 2 1 + 10 sin 2 π y i + 1 + y n 1 2 + i = 1 n u x i , 10 , 100 , 4 ,       y i = 1 + x i + 1 4 [−50, 50]030Multimodal
f 8 f x = 0.1 sin 2 3 π x 1 + i = 1 n x i 1 2 1 + sin 2 3 π x i + 1 + x n 1 2 1 + sin 2 2 π x n + i = 1 n u x i ,   5 ,   100 ,   4 [−50, 50]030Multimodal
Table 2. Comparison results between JS, MFO, and JS-MFO for f1f8.
Table 2. Comparison results between JS, MFO, and JS-MFO for f1f8.
AlgorithmIndexUnimodal FunctionsMultimodal Functions
f 1 f 2 f 3 f 4 f 5 f 6 f 7 f 8
JSMean2.00 × 10−182.03 × 10−190.71002.67 × 10−103.49 × 10−1012.5858.68 × 10−50.0098
SD2.84 × 10−182.68 × 10−190.88352.65 × 10−101.49 × 10−109.08851.94 × 10−40.0368
Best9.14 × 10−201.57 × 10−200.01793.97 × 10−111.50 × 10−103.05741.36 × 10−61.57 × 10−5
Time (s)1.24881.23321.07411.33841.47801.42032.49552.7288
MFOMean2.67 × 103383.5271.88 × 10432.093611.1412139.564.50981.36 × 107
SD5.20 × 103505.0983.64 × 10420.39438.738727.522.88507.48 × 107
Best0.31040.0524178.54370.10660.132377.670.47060.2722
Time (s)0.72180.71270.82120.82021.23630.91982.04992.2922
JS-MFOMean6.53 × 10−834.94 × 10−600.05202.59 × 10−393.13 × 10−153.89 × 10−135.08 × 10−292.41 × 10−25
SD3.43 × 10−821.79 × 10−590.12961.40 × 10−381.74 × 10−151.98 × 10−121.88 × 10−281.31 × 10−24
Best2.38 × 10−891.19 × 10−721.00 × 10−81.25 × 10−468.88 × 10−1603.63 × 10−321.62 × 10−32
Time (s)1.78421.75611.94011.94342.63282.06384.16974.6286
Table 3. Comparative statistical study between JSO, MFO, and JS-MFO for case 1, case 3, and case 4 of OPF problem.
Table 3. Comparative statistical study between JSO, MFO, and JS-MFO for case 1, case 3, and case 4 of OPF problem.
Study CaseIndexMethod
JSOMFOJS-MFO
Case 1
Quadratic total fuel cost (USD/h)
Minimum801.023800.79799.08
Average801.84801.92799.09
Maximum803.41802.87799.14
SD0.7811.0400.0133
Case 3
Total fuel cost with valve point effect (USD/h)
Minimum929.48924.02918.07
Average943.60929.06925.61
Maximum959.14984.18952.54
SD12.0713.5912.05
Case 4
Total active losses (MW)
Minimum3.2323.0492.881
Average3.1893.4432.909
Maximum3.2663.753.095
SD0.0530.21470.033
Table 4. Optimal tuning of control variables and objectives for Cases 1–3—IEEE 30-bus system.
Table 4. Optimal tuning of control variables and objectives for Cases 1–3—IEEE 30-bus system.
Control Variables/ObjectivesCase 1Case 2Case 3
JSMFOJS-MFOJSMFOJS-MFOJSMFOJS-MFO
PG1(MW)177.68178.13177.1564.52963.80463.921199.60199.60199.589
PG2(MW)48.68949.14648.72767.39367.7967.49020.14420.00020.002
PG5(MW)21.44821.82421.32949.95850.0049.99919.62522.29022.027
PG8(MW)19.43619.30720.89734.97434.9934.99924.92027.90423.029
PG11(MW)13.06911.74011.92429.98330.0029.99915.03810.00314.982
PG13(MW)12.27612.42612.00539.99340.0039.99914.52814.57613.192
VG1(p.u)1.08041.10001.10001.06001.0591.0991.06601.0231.097
VG2(p.u)1.06071.08111.08761.05511.0541.0951.03520.9991.076
VG5(p.u)1.02681.04911.06051.01981.0351.0760.98670.9501.043
VG8(p.u)1.03531.05031.06841.04061.0471.0831.00320.9661.058
VG11(p.u)1.06691.03631.10001.06131.1001.0971.06061.0991.096
VG13(p.u)1.06330.95001.09911.04181.1001.0990.99371.0871.098
T11(6_9)1.00921.09191.06660.98501.0191.0501.03990.9210.988
T12(6_10)0.96001.04050.90611.02590.9010.9140.98710.9851.022
T15(4_12)1.02361.06001.00991.03530.9720.9911.00161.0821.031
T36(28_27)1.00801.07450.97480.99180.9500.9720.99320.9310.968
QC10(MVar)3.60545.00004.99141.83903.6754.9841.645202.657
QC12(MVar)3.81403.43134.93413.616804.7391.93520.0173.497
QC15(MVar)2.05834.97444.92633.55155.00004.7642.157001.655
QC17(MVar)2.7145 04.95441.88044.7164.9001.56540.1372.442
QC20(MVar)2.87584.99884.85322.78212.8394.9063.18754.9314.492
QC21(MVar)2.59084.74934.99833.84694.9994.9632.88675.001.828
QC23(MVar)2.37904.94604.06533.89021.7104.3694.21004.1933.439
QC24(MVar)2.91895.00004.98392.93505.0004.79453.41873.2242.184
QC29(MVar)1.78365.00003.21823.78813.7142.90244.00900.0433.073
F C ($/h)801.02800.79799.085944.30944.6943.690922.48924.02918.07
F E m (ton/h)0.36730.36890.36630.20490.20480.20470.43950.44030.439
F P l o s s (MW)9.20269.1798.63683.43243.1953.009510.46610.9759.424
F V D (p.u)0.5320.98891.64860.39941.4861.92570.47060.3291.195
F L m a x 0.13830.11450.11820.13390.1200.11580.14330.14310.124
Table 5. Optimum tuning of control variables and objectives for Cases 4, 7, and 8—IEEE 30-bus system.
Table 5. Optimum tuning of control variables and objectives for Cases 4, 7, and 8—IEEE 30-bus system.
Control Variables and Objectives Case 4Case 7Case 8
JSMFOJS-MFOJSMFOJS-MFOJSMFOJS-MFO
PG1(MW)63.59151.48351.349178.73177.89176.174174.09179.01177.466
PG2(MW)68.15080.0079.98848.71148.7648.82648.90449.48348.746
PG5(MW)49.99350.0049.99722.00621.0321.97021.63222.10221.785
PG8(MW)34.96235.0034.95018.21321.4921.35622.19516.76620.639
PG11(MW)29.96929.96529.99812.43312.23912.72213.18211.95911.458
PG13(MW)39.96539.99939.99713.24612.00812.13312.35113.31712.015
VG1(p.u)1.0871.0641.0991.04971.0331.0441.0861.0921.100
VG2(p.u)1.0761.0561.0961.04291.0091.0261.0661.071.084
VG5(p.u)1.0551.0391.0771.01340.9961.0101.0231.0241.053
VG8(p.u)1.0701.0451.0841.00001.0041.00651.0451.0441.058
VG11(p.u)1.0691.09931.1001.00381.0731.05751.0831.1001.099
VG13(p.u)1.0811.09941.0991.04461.01870.98991.0781.0931.097
T11(6_9)0.9760.97351.0721.07801.10001.07791.0780.9400.990
T12(6_10)1.0250.93230.9001.06900.91100.90101.0690.9000.933
T15(4_12)0.9960.96930.9991.03200.99870.94161.0320.9530.969
T36(28_27)0.9880.94480.9781.06800.96150.96711.0680.9210.948
QC10(MVar)3.3883.70294.8702.93445.0004.17394.1194.9894.979
QC12(MVar)3.1125.00004.9840.62090.4680.38903.6514.9954.999
QC15(MVar)3.7135.00004.9464.77535.0004.90324.4945.0004.983
QC17(MVar)3.0725.00004.9642.85385.0000.08443.3504.6604.648
QC20(MVar)2.8383.92154.8154.45684.7814.98592.2085.0004.552
QC21(MVar)2.5604.07824.8102.52505.0004.97133.4171.1534.985
QC23(MVar)2.8603.4353.8933.83495.0004.99704.1284.8994.800
QC24(MVar)3.1144.7424.9774.00134.9014.98363.7874.9644.897
QC29(MVar)2.9153.9152.8652.72911.7562.37404.33404.680
F C ($/h)945.10967.46967.02803.57803.93803.554801.16800.96799.318
F E m (ton/h)0.20480.2070.2070.37000.36810.36310.3570.37120.367
F P l o s s (MW)3.2323.0492.8819.949110.0439.78458.9589.2448.711
F V D (p.u)1.1741.6391.88060.14800.11360.09931.2821.8721.8981
F L m a x 0.1260.1180.11660.13710.13610.13700.1170.11820.1141
Table 6. Comparisons of the results obtained for case 1 of the IEEE 30-bus system.
Table 6. Comparisons of the results obtained for case 1 of the IEEE 30-bus system.
Technique Total Fuel Cost (USD/h)Technique Total Fuel Cost (USD/h)
JS-MFO799.085DA [13]802.316
JS801.02AGTLBO [14]800.481
MFO800.79PSO-PS [21]799.8723
NISSO [35]799.762DA-PSO [23]802.124
FFA [15]802.130HF-PSO [24]799.123
MSA [15]802.223HF-ABC [26]800.212
ESCA [36]800.219HFA-JAYA [27]800.480
SOS [35]801.5733
HHO [16]801.829
Table 7. Comparison report provided for case 2 of the IEEE 30-bus system.
Table 7. Comparison report provided for case 2 of the IEEE 30-bus system.
Technique Total Gas Emissions (ton/h)Technique Total Gas Emissions (ton/h)
JS-MFO0.2047AGTLBO [14]0.2048
JS0.2049SSO [35]0.2315
MFO0.2048NISSO [35]0.2048
HHO [16]0.2850SSA [37]0.205
TFWO [17]0.2050PSO-SSA [37]0.205
DA-PSO [23]0.2048
Table 8. Comparison of the results obtained for case 3 of the IEEE 30-bus system.
Table 8. Comparison of the results obtained for case 3 of the IEEE 30-bus system.
Technique Total Fuel Cost (USD/h)Technique Total Fuel Cost (USD/h)
JS-MFO918.093TLBO [39]919.394
JS922.48Gbest-ABC [40]931.745
MFO924.02GABC1 [41]919.597
ISSA [38]919.191GABC2 [41]918.435
SSA [38]920.706
IHS [38]919.843
Table 9. Comparison of the results for Case 4 of the IEEE 30-bus system.
Table 9. Comparison of the results for Case 4 of the IEEE 30-bus system.
TechniqueTotal Transmission Losses (MW)TechniqueTotal Transmission Losses (MW)
JS-MFO2.881DA [23]3.198
JS3.232DA-PSO [23]3.189
MFO3.049HFA-JAYA [27]4.529
AGTLBO [14]3.090SS0 [35]3.911
FFA [15]3.643NISSO [35]2.945
MSA [15]3.649ESCA [36]3.021
DE-HS [22]3.054AMTPG-JAYA [42]3.080
Table 10. Comparisons of the results obtained for case 7 of the IEEE 30-bus system.
Table 10. Comparisons of the results obtained for case 7 of the IEEE 30-bus system.
TechniqueTotal Fuel Cost FC (USD/h)Total Voltage Deviation FVD (p.u)
JS-MFO803.55490.0993
JS803.57270.1480
MFO803.930.1136
AGTLBO [14]803.7380.0947
DA-PSO [23]803.660.1163
HFA-JAYA [27]803.70360.0948
ESCA [36]804.99680.09163
PSO-SSA [37]803.9890.0940
MVO [43]803.9080.1056
ECHT-DE [44]803.7190.0945
Table 11. Comparison of the results obtained for case 8 of the IEEE 30-bus system.
Table 11. Comparison of the results obtained for case 8 of the IEEE 30-bus system.
TechniqueTotal Fuel Cost FC (USD/h)Voltage Stability Index FLmax
JS-MFO799.31870.1141
JS801.16210.1171
MFO800.9630.1182
ESCA [36] 800.4100.1224
PSO-SSA [38] 830.3520.1250
AMTPG-JAYA [42] 840.9010.1240
MVO [43]802.4660.1146
ECHT-DE [44]800.4200.1374
Table 12. Settings of control variables and objectives of RIM 27-bus system for cases 4–6.
Table 12. Settings of control variables and objectives of RIM 27-bus system for cases 4–6.
Control Variables and ObjectiveCase 4Case 5Case 6
JSMFOJS-MFOJSMFOJS-MFOJSMFOJS-MFO
PG1(MW)85.6782.04483.627990.66749.7536.953077.5546108.26109.0601
PG8(MW)20.4921.97121.857020.00550.0049.802718.1765.5536.5210
PG10(MW)19.2719.32717.815111.78914.4227.615025.608129.9510.4073
PG23(MW)9.3458.5819.02374.68299.9363.015212.612014.663.5084
PG24(MW)33.4335.93935.992630.70036.0035.959519.638510.7123.6076
PG25(MW)13.4814.07313.312324.3332.74529.999229.096212.8629.9499
VG1(p.u)0.9340.96180.92000.95700.9540.95590.95180.93790.9445
VG8(p.u)0.9410.97050.92710.97440.9920.99360.96470.93980.9392
VG10(p.u)0.9380.96630.92350.96480.9660.97240.96460.95000.9492
VG23(p.u)0.9350.96410.92190.95740.9590.95740.95110.94220.9411
VG24(p.u)0.9290.95040.92000.93740.9410.94150.94210.93130.9303
VG25(p.u)0.9700.96910.94310.97190.9200.99930.92400.92000.9673
QC12(MVar)−19.67−19.948−19.994−19.73−20−19.998−16.248−19.34−19.903
QC14(MVar)−9.774−9.9932−9.9996−9.907−10−9.971−7.9491−9.760−7.072
QC16(MVar)−24.97−25.000−24.998−24.944−25−24.999−24.216−25.00−15.7419
F P l o s s e s (MW)3.3543.56213.25623.93574.4704.92354.35563.88664.9284
F V D (p.u)0.9160.60451.24020.65470.7320.61440.80050.93060.9650
F L m a x 0.4260.37960.45970.38690.3910.38760.38610.42010.3389
Table 13. Reactive power generation in the optimal state for each case—RIM 27-bus power system.
Table 13. Reactive power generation in the optimal state for each case—RIM 27-bus power system.
State VariablesCase 4Case 5Case 6
Generated Reactive PowerJS-MFOJS-MFOJS-MFO
QG1(MVAR)−79.6379−104.3660−99.8009
QG8(MVAR)−13.7852−6.2700−17.1219
QG10(MVAR)−15.8118−6.1590−8.9077
QG23(MVAR)−7.5879−5.8648−11.5592
QG24(MVAR)−5.9939−25.0052−17.6271
QG25(MVAR)0.32725.10914.5488
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

Mayouf, C.; Salhi, A.; Haidara, F.; Aroua, F.Z.; El-Sehiemy, R.A.; Naimi, D.; Aya, C.; Kane, C.S.E. Solving Optimal Power Flow Using New Efficient Hybrid Jellyfish Search and Moth Flame Optimization Algorithms. Algorithms 2024, 17, 438. https://doi.org/10.3390/a17100438

AMA Style

Mayouf C, Salhi A, Haidara F, Aroua FZ, El-Sehiemy RA, Naimi D, Aya C, Kane CSE. Solving Optimal Power Flow Using New Efficient Hybrid Jellyfish Search and Moth Flame Optimization Algorithms. Algorithms. 2024; 17(10):438. https://doi.org/10.3390/a17100438

Chicago/Turabian Style

Mayouf, Chiva, Ahmed Salhi, Fanta Haidara, Fatima Zahra Aroua, Ragab A. El-Sehiemy, Djemai Naimi, Chouaib Aya, and Cheikh Sidi Ethmane Kane. 2024. "Solving Optimal Power Flow Using New Efficient Hybrid Jellyfish Search and Moth Flame Optimization Algorithms" Algorithms 17, no. 10: 438. https://doi.org/10.3390/a17100438

APA Style

Mayouf, C., Salhi, A., Haidara, F., Aroua, F. Z., El-Sehiemy, R. A., Naimi, D., Aya, C., & Kane, C. S. E. (2024). Solving Optimal Power Flow Using New Efficient Hybrid Jellyfish Search and Moth Flame Optimization Algorithms. Algorithms, 17(10), 438. https://doi.org/10.3390/a17100438

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop