Next Article in Journal
Collaborative Gold Mining Algorithm: An Optimization Algorithm Based on the Natural Gold Mining Process
Previous Article in Journal
Simultaneous Beam Forming and Focusing Using a Checkerboard Anisotropic Surface
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Novel MOGNDO Algorithm for Security-Constrained Optimal Power Flow Problems

1
Department of Electrical Engineering, Shri K.J. Polytechnic, Bharuch 392 001, India
2
Department of Computer Science and Engineering, Vel Tech Rangarajan Dr Sangunthala R&D Institute of Science and Technology, Chennai 600 062, India
3
Department of Control Systems and Instrumentation, Faculty of Mechanical Engineering, VSB-Technical University of Ostrava, 17. Listopadu 2172/15, 708 00 Ostrava, Czech Republic
4
School of Mechanical Engineering, Vellore Institute of Technology, Chennai 600 127, India
5
Rajasthan Rajya Vidyut Prasaran Nigam, Losal, Sikar 332 025, India
*
Author to whom correspondence should be addressed.
Electronics 2022, 11(22), 3825; https://doi.org/10.3390/electronics11223825
Submission received: 26 September 2022 / Revised: 15 November 2022 / Accepted: 16 November 2022 / Published: 21 November 2022

Abstract

:
The current research investigates a new and unique Multi-Objective Generalized Normal Distribution Optimization (MOGNDO) algorithm for solving large-scale Optimal Power Flow (OPF) problems of complex power systems, including renewable energy sources and Flexible AC Transmission Systems (FACTS). A recently reported single-objective generalized normal distribution optimization algorithm is transformed into the MOGNDO algorithm using the nondominated sorting and crowding distancing mechanisms. The OPF problem gets even more challenging when sources of renewable energy are integrated into the grid system, which are unreliable and fluctuating. FACTS devices are also being used more frequently in contemporary power networks to assist in reducing network demand and congestion. In this study, a stochastic wind power source was used with different FACTS devices, including a static VAR compensator, a thyristor- driven series compensator, and a thyristor—driven phase shifter, together with an IEEE-30 bus system. Positions and ratings of the FACTS devices can be intended to reduce the system’s overall fuel cost. Weibull probability density curves were used to highlight the stochastic character of the wind energy source. The best compromise solutions were obtained using a fuzzy decision-making approach. The results obtained on a modified IEEE-30 bus system were compared with other well-known optimization algorithms, and the obtained results proved that MOGNDO has improved convergence, diversity, and spread behavior across PFs.

1. Introduction

Constraint-based optimization problems with multiple objectives are the most prevalent type. In contrast to single-objective optimization problems, multi-objective optimization problems have a wide variety of optimal solutions. The PF is an assortment of perfect responses [1,2]. A multi-objective optimization approach must be able to locate solutions that are uniform in the generated PFs and are workable optimum solutions to address multi-objective problems [3]. Multi-objective optimization approaches are challenged by the simultaneous achievement of these many objectives [4]. MH algorithms are typically tested on simpler, well-known optimization scenarios. However, unlike classic search problems, engineering design tasks can have different specifications. Modifying and developing the algorithm for them is the most effective way to optimize for them. The realm of application of multi-objective optimization algorithms is quite vast, ranging from machining processes [5,6], to vehicle routing [7], to optimizing AI systems [8].
Power systems researchers have been seeking solutions to the OPF challenges for many decades. One issue with managing power systems and making plans for modern electrical energy networks is working with systems that use non-conventional sources of energy. Ullah et al. [9] developed a hybrid phasor particle swarm optimization and gravitational search algorithm to address the OPF problem in the wind and solar energy systems that are connected to electrical power grids, while accounting for the control variables (PPSOGSA). For the OPF problem with wind and solar systems, the developed PPSOGSA algorithm produced outstanding and helpful results. In Elattar’s research, the OPF problem was principally modelled mathematically using a combined heat and power system with stochastic wind energy. On an IEEE 30-bus test system under various test situations, the suggested method was assessed. The formulation of the OPF problem covering energy sources and the suggested approach to solve it produced effective answers in contrast to pre-existing algorithms [10] that were used to address similar problems. Anongpun et al. [11] used IEEE 30- and 118-bus test systems to study the use of enhanced particle swarm optimization (PSO) to solve a multi-objective OPF problem using a wind energy system that combined chaotic mutation and stochastic weights. When compared to the other algorithms included in their study, the suggested method produced better results. Salkuti [12] used the glowworm swarm optimization method to offer a solution to a multi-objective OPF problem requiring a modern electrical energy system that utilized wind energy. On IEEE 30- and 300-bus test systems, the methodology was examined in many operational scenarios. The simulation findings indicated that the proposed methodology might offer an alternative. A flower pollination algorithm was used by Kathiravan et al. [13] to address the OPF problem using coal-based, wind, and solar energy systems. In a variety of test situations, the authors used their method to test systems for the IEEE 30-bus and Indian utility 30-bus. Duman et al. [14] used differential evolutionary particle swarm optimization to address the OPF problem (DEEPSO) with manageable wind and solar (PV) energy sources. IEEE 30-, 57-, and 118-bus test systems were used to evaluate the DEEPSO technique to explore the issue under various objective functions. DEEPSO produced better simulation results when compared to the other optimization approaches that were looked at [14]. They recommended using FACTS devices such as a thyristor-controlled phase shifter (TCPS) and a thyristor-controlled series capacitor (TCSC) to solve the OPF problem. To account for the uncertainties associated with wind energy installation, they used chaotic maps and a modified version of the PSOGSA (particle swarm optimization and gravitational search algorithm). The method presented [15] appears to be a potential approach for a solution based on the findings of simulations. Biswas et al. [16] solved the OPF challenge, which incorporated coal-based, wind, solar, and small-hydro energy sources coupled to IEEE 30-bus test systems, by running multiple rounds of the multi-objective evolutionary algorithm. The constrained multi-objective population extremal optimization (CMOPEO) technique was used to handle the wind and solar-integrated OPF problem by Chen et al. [17], who also tested the method on an IEEE 30-bus for different scenarios. Additionally, research has been done on the evolutionary particle swarm optimization (EPSO) [18], the hybrid differential evolution and symbiotic organisms search algorithm (HMICA-SQP) [19], the success-history-based adaptation of differential evolution with superiority of feasible solutions (SHADE-SF) [20], and the hybrid modified imperialist competitive algorithm and sequential quadratic programming algorithm (HMICA-SQP) [21]. Pandya and Jariwala [22] addressed single and multi-objective OPF issues by integrating with various sustainable energy sources using recently developed metaheuristics algorithms. Biswas et al. [23] analyzed the integration of three FACTS devices, wind turbines, and coal-fired power plants. The success history-based adaptive differential evolution (SHADE) method was used to conduct the investigation. According to the “No Free Lunch (NFL)” theorem [24], no metaheuristic can solve every issue that occurs in real-world situations. This theorem has opened the door to the creation of both novel metaheuristic techniques and improvements to existing ones.
When using the Generalized Normal Distribution Optimization method [25] in multi-objective optimization scenarios, several things need to be taken into account. The initial problem in multi-objective generalized normal distribution optimization is balancing convergence and divergence archives. The algorithm’s generated solution set is filtered according to a certain quality metric, and the non-dominating solution set is kept in separate, external archives. In the literature, there are numerous recommendations for constructing archives that could be employed in the multi-objective Generalized Normal Distribution Optimization method. The pre-defined maximum size archives are widely employed, since more non-dominating solutions can emerge quickly. The algorithm’s ability to begin with the straightforward determination of the required population size and the terminal condition is the most noticeable feature of the GNDO. The location of the person is automatically changed by the generalized normal distribution function (GNDO), which has a simple construction. The benefits and drawbacks of the GNDO algorithm are as follows:
  • It offers a faster and smoother convergence, especially for difficult problems, and it strikes the perfect balance between exploration and exploitation.
  • Local minima are less likely to become entangled in relaxed convergence.
  • Effortlessly simple, adaptable, and simple to use
  • The traditional GNDO may have issues with convergence trends or become stuck in narrow, deceptive optima for challenging optimization tasks, such as high-dimensional and multimodal problems.
Currently, both conventional and non-conventional energy sources require more studies. The current body of research recommends using coal-based plus wind and FACTS devices, combined with single and multi-objective optimum power flow (MOOPF) problems. The conventional IEEE 30-bus network has been altered to include non-conventional sources for research purposes. Using Weibull PDF, non-conventional units’ stochastic behaviors are calculated. The generating cost is suitably adjusted to account for reserve cost if these stochastic units are over-estimated and adjusted for penalty cost in the case that they are underestimated. Using the Generalized Normal Distribution Optimization method, Pareto solution clusters are discovered for the multi-objective problem. The following is a list of the contributions made by this study:
  • This work focuses on the mathematical modelling of the single and multiple-objective OPF issue modelled, which takes into account both conventional units and non-conventional sources of energy units, as well as FACTS devices.
  • The appropriate probability density functions (PDFs) are modeled in the second stage to describe the wind power plants’ random behavior.
  • Stochastic non-conventional sources of energy sources are among the single and multiple-objective OPF issues for which the Non-Dominated Sorting Generalized Normal Distribution Optimization (NSGNDO) technique is used to develop solutions.
  • Studies and performance evaluations of the MOGNDO algorithm using empirical comparisons are conducted.
The notion of the mathematical models for coal-based power, wind power, and FACTS devices is presented in Section 2 of the study. An explanation of the objectives that need to be optimized is included in Section 3. Section 4 provides an explanation and illustrations of the multi-objective GNDO technique. Section 5 presents numerical results and discussion, and Section 6 provides concluding remarks.

2. Mathematical Representations

The case studies presented here restructure the original IEEE 30-bus test apparatus. The modified approach incorporates wind turbines and FACTS devices, and is listed in Table 1. The equipment used for the analysis is depicted in Figure 1. The placement and ratings of FACTS devices are depicted in the diagram with dotted lines because they have been optimized. The section below provides information on the costs of traditional coal-based production facilities and plants using non-conventional sources of energy.

2.1. Cost of Coal-Based Power Units

The generalized quadratic equation for the calculation of generation cost is expressed in (1) in $/h [23]:
C T 0 P T G = i = 1 N T G a i + b i P T G i + c i P T G i 2
For a more practical case, the valve point effect included:
C T P T G = i = 1 N T G a i + b i P T G i + c i P T G i 2 + d i × sin e i × P T G i m i n P T G i
The values of both coal-based price constants and emanation constants with various scenarios are shown in [23].

2.2. Toxic Gas Emanation

Polluted gases are released by using coal-based plants. So, toxic gas emanations in tons per hour can be determined as (in ton/h):
F 2 = t , E = i = 1 N T G α i + β i P T G i + γ i P T G i 2 × 0.01 + ω i e μ i P T G i
The toxic gas emanation constants of coal-based power plants are taken from [22].

2.3. Direct Cost of Stochastic Non-Conventional Sources Plants

It is particularly challenging to integrate non-conventional energy sources into the power grid since they are stochastic. The independent system operator (ISO) is responsible for managing these non-conventional energy sources. Due to this, the private operator must contract with the grid or ISO for a specific quantity of planned power. The scheduled electricity must be maintained by the ISO scheduled power. If these non-conventional sources are unable to maintain the planned power, the ISO is liable for the absence of power. So, if a need arises, there are spinning reserve requirements. This spinning reserve increases costs for the ISO, and this circumstance is known as an overestimation of non-conventional sources. Conversely, if non-conventional sources formed more energy than was planned, it might go to waste, due to underuse. Therefore, the ISO must accept the penalty charge. The scheduled power cost, the overestimation cost caused by the spinning reserve, and the penalized cost caused by the underestimation are the three costs related to electricity. The direct cost linked to wind farms is demonstrated with the P w s scheduled power from the same sources as:
C w P w s = g w P w s

2.4. Indeterminate Non-Conventional Sources of Wind Power Cost

Due to the erratic nature of wind, the wind farm occasionally produces less energy than expected. This means that if demand increases, it needs the spinning reserve to maintain the agreed-upon amount of scheduled power. It is sometimes feasible that the real power generated by wind farms won’t be enough to meet demand and will have lower values. Such power is referred to as exaggerated power by an ambiguous resource. To control this kind of uncertainty and provide end users with a reliable power source, the network ISO operates spinning reserves. The price of hiring a backup generator to supply the overestimated power is known as the reserve cost.
Reserve cost for the wind unit is formulated by:
C R w P w s P w a v = K R w P w s P w a v = K R w 0 P w s P w s p w f w p w d p w
The possibility exists that the wind farm will generate more power than is required, which is the opposite of the overestimation scenario. Underestimated power is the term used to describe such a situation. If there is no provision for managing the output power from coal-based units, the excess power will be lost. Regarding the extra power, the ISO needs to be penalized.
The penalty charge for the wind unit is given by:
C P w P w a v P w s = K P w P w a v P w s = K P w P w s P w r p w P w s f w p w d p w

2.5. Uncertainty Models of Stochastic Wind Units

In the redesigned IEEE-30, the wind power generating units installed at buses 5 and 11, which were originally thermal generators, were replaced. It should be noted that, as for a comparison point of view with the published reference article [23], in this paper, the thermal DGs were also replaced with the wind turbines. This will ensure compatibility of the results obtained by the proposed algorithm to the already published research article [23]. The scale (c) and shape (k) constants for the proposed Weibull model are detailed in Table 2.
The Weibull curve and wind frequency distributions in Figure 2 (for the bus 5 wind plant) and Figure 3 (for the bus 11 wind plant) were produced using 8000 Monte-Carlo settings. The standard provided explains the need for wind turbine design and specifies the maximum turbulence class IA that is confirmed to operate at the highest yearly average wind velocity of 10 m/s at hub height. The formed shape k and scale ( c ) parameters of wind farms are given particular attention, because their highest Weibull mean value is fixed at around 10. It is commonly known that the wind speed distribution follows the Weibull PDF curve.
The following formula can be used to calculate the probability of wind velocity v, in m/s, pursuing the Weibull PDF with shape factor (k) and scale factor (c) [23]:
f v v = k c v c k 1 e v c k f o r       0 < v <
The Weibull distribution’s mean is given as follows [23]:
M w b l = c Γ 1 + k 1
and the gamma function Γ x is expressed in Equation (9):
Γ x = 0 e t t x 1 d t

2.6. Average Power Calculation for Wind Plants

The combined outputs of the 25 turbines in the farm are taken as the wind unit connected at bus 5. Every turbine has a 3 MW output rating. The wind velocity affects the wind turbine’s precise output, which varies. We used the following equations to express turbine output power in terms of wind velocity (v) [23]:
p w v = 0 ,                               for   v v i n   a n d   v v o u t p w r v v i n v r v i n                           for   v i n v v r p w r                                                       for   v r < v v o u t
The Enercon E82-E4 design specification is referred to for the 3 MW wind turbine. The various speeds are v i n = 3 m/s, v r = 16 m/s, and v o u t = 25 m/s.

2.7. Wind Power Probabilities Calculation

In certain ranges of wind speeds, uncertain wind generation is noticeable. The generated power would be 0 if the wind speed was greater than or less than the cut-out speed or cut-in speed. The turbine thereby produces the specified amount of power within the range of the rated and cut-out wind speeds. These are possible ways to describe the likelihood of these areas [23]:
f w p w p w = 0 = 1 exp v i n c k + exp v o u t c k
f w p w p w = p w r = exp v r c k exp v o u t c k
Between the cut-in velocity and the rated velocity of the wind, the wind production remains constant. The following can be used to express the likelihood of the continuous zone [23]:
f w p w = β v r v i n α β p w r v i n + p w p w r v r v i n β 1 exp v i n + p w p w r v r v i n α β

2.8. Thyristor-Controlled Series Compensator (TCSC) Modeling

The basic circuitry of the TCSC is depicted in Figure 4. It consists of a fixed series capacitor (XC) and a reactor (XL) operated by a thyristor. For the TCSC to function as a variable capacitive reactance, reactance X C < X L is taken into consideration. By varying the firing angle ( α ) of the thyristors, the inductive reactance is changed, and for high values of inductive reactance, the least corresponding capacitive reactance is produced (Open circuit Inductive branch). As a result, the TCSC’s effective reactance with constant capacitive reactance X C and variable inductive reactance X L α can be written as [23]:
X T C S C α = X C X L α X L α X C = j X C
The TCSC static model, which is situated in the path between buses m and n, is shown in Figure 4. Following the TCSC’s integration (described as a variable capacitive reactance mode), the transmission line’s adjusted reactance ( X e q ) is given by [23]:
X e q = X m n X T C S C = 1 τ X m n
where
τ = X T C S C X m n
The power flow equations of the line incorporating the TCSC are written as [23]:
P m n = V m 2 g m n V m V n g m n cos δ m δ n . V m V n b m n sin δ m δ n
Q m n = V m 2 b m n V m V n g m n sin δ m δ n . + V m V n b m n cos δ m δ n
P n m = V n 2 g m n V m V n g m n cos δ m δ n . + V m V n b m n sin δ m δ n
Q n m = V n 2 b m n + V m V n g m n sin δ m δ n . + V m V n b m n cos δ m δ n
where
g m n = r m n r m n 2 + ( x m n x c ) 2
b m n = x m n x c r m n 2 + ( x m n x c ) 2

2.9. Model of Thyristor-Controlled Phase Shifter (TCPS)

Figure 5 displays the model of the TCPS placed between the line that connects buses m and n . The power flow equations of the line can be expressed as below, assuming that is the phase shift angle ϕ is introduced by the TCPS:
P m n = V m 2 g m n cos 2 ϕ V m V n cos ϕ g m n cos δ m δ n + ϕ . + b m n sin δ m δ n + ϕ
Q m n = V m 2 b m n cos 2 ϕ V m V n cos ϕ [ g m n sin δ m δ n + ϕ . b m n cos δ m δ n + ϕ ]
P n m = V n 2 g m n V m V n cos ϕ [ g m n cos δ m δ n + ϕ . b m n sin δ m δ n + φ ]
Q n m = V n 2 b m n + V m V n cos ϕ [ g m n sin δ m δ n + ϕ . + b m n cos δ m δ n + φ ]
The inserted actual and reactive power of the TCPS at bus m and n is [23]:
P m s = g m n V m 2 tan 2 ϕ V m V n tan ϕ [ g m n sin δ m δ n . b m n cos δ m δ n ]
Q m s = b m n V m 2 tan 2 ϕ + V m V n tan ϕ [ g m n cos δ m δ n . + b m n sin δ m δ n ]
P n s = V m V n tan ϕ g m n sin δ m δ n + b m n cos δ m δ n
Q n s = V m V n tan ϕ g m n cos δ m δ n b m n sin δ m δ n

2.10. Model of Static VAR Compensator (SVC)

The basic circuit architecture and the SVC model are depicted in Figure 6. It is made up of a thyristor-controlled reactor ( X L = ω L ) and a fixed capacitor ( X C = 1 / ω C ) . By changing the thyristor firing angle α , the reactance can be changed. The equivalent susceptibility is computed as:
B e q = B L α + B C
where
  B L α = 1 ω L 1 2 α π , B c = ω × C
The reactive power offered by the SVC can be expressed in terms within the context of power flow:
Q S V C = V m 2 · B S V C

3. Objectives of Optimization

The best active power allocation and the best VAR power allocation are goals in the OPF. The following are some examples of how this section incorporates the goals of wind power flow optimization.

3.1. Reducing Overall Costs While Using Non-Conventional Energy Sources

Objective 1:
Minimizing the whole cost is the first goal. Direct, reserve, and penalty charges for non-conventional resources are added to the coal-based unit cost to determine the overall generation cost. Therefore, the comprehensive cost for coal-based and wind power plants is denoted as:
Minimize—
F 1 = C T o t = C T P T G + C w P w s + C R w P w s P w a v + C P w P w a v P w s

3.2. Reduction of Voltage Variation with the Use of Non-Conventional Energy Sources

One of the most crucial safety and administrative superiority lists is the bus voltage. By restricting the voltage deviations of the PQ bus from 1.0 for each unit, the improving voltage profile will be acquired. The objective function is going to come from:
Objective 2: Minimize—
F 2 = i = 1 N p q v i 1.0

3.3. Minimization of APL Including Non-Conventional Energy Sources

The optimization of actual power losses P L O S S (MW) maybe calculated by:
Objective 3: Minimize—
F 3 = P L O S S = i = 1 N B P G i i = 1 N B P D i

3.4. Enhancement of VSI Including Non-Conventional Energy Sources

The L m a x index is the most important indicator for evaluating each bus’s voltage constancy margin, since it keeps the voltage constant within a reasonable range during typical operation. For every PQ bus, the L m a x index offers a scalar number. Between ‘0’ (no load) and ‘1’ is where the L m a x index is located (voltage collapse). The following formula is used to get the jth bus’s voltage collapse indicator amount:
L j = 1 i = 1 N g F j i V i V j       j = 1 , 2 , , N L
F j i = Y 1 1 Y 2
The objective function of voltage stability enhancement is written by:
F 4 = L = max L j         j = 1 , 2 , , N L

3.5. Minimization of Entire Gross Cost Including Non-Conventional Energy Resources

The generating cost is significantly higher in the latter scenario, whereas the loss is greater in the former, as shown by Objectives 1 and 3. The requirement for an aim that includes both the cost and the loss is increased by this very circumstance. Making a cost model that converts the loss into an equivalent energy cost is a straightforward way to take into consideration both goals. The price of energy taken into account in this analysis is $0.10 per kWh. The goal of gross cost in dollars per hour might be stated as follows:
F 5 = C T o t + P L O S S 10000 0.10

3.6. Equality Constraints

Power flow equations provide equality boundaries, and demonstrate that both real and fictitious power generated in a system should fulfill the load demand and system losses:
P G i P D i V i j = 1 N B V j G i j cos δ i j + B i j sin δ i j = 0       i N B
Q G i Q D i V i j = 1 N B V j G i j sin δ i j B i j cos δ i j = 0       i N B

3.7. Inequality Constraints

Inequality bounds are the operational boundaries of devices and the security bounds of lines and PQ buses.
  • Generator bounds:
P T G i m i n P T G i P T G i m a x   i N T G
P ws m i n P w s P ws m a x
Q T G i m i n Q T G i Q T G i m a x   i N T G
Q ws m i n Q w s Q ws m a x
V G i m i n V G i V G i m a x , i = 1 , , N G
  • Security bounds:
V L p m i n V L p V L p m a x , p = 1 , , N L
S l q S l q m a x , q = 1 , , n l
  • FACTS devices bounds:
τ TCSCm m i n τ T C S C m τ TCSCm m a x m N T C S C
ϕ TCPSn m i n τ T C P S n τ TCPSn m a x n N T C P S
Q S V C j m i n Q S V C j Q S V C j m a x j N S V C
The real power output limits of coal-based and wind units are shown in Equations (43) and (44), respectively. Then, Equations (45) and (46), which show the VAR power size of producing units, are used. The overall voltage regulator buses are shown in NG. PV bus voltage restrictions are shown in Equation (47), while PQ bus voltage restrictions are shown in Equation (48), where NL is the number of PQ buses. When NL is the total number of lines in a system, Equation (49) can be used to calculate line loading limitations. Equations (50)–(52) show the limits of the TCSC, TCPS, and SVC devices, respectively.

4. Generalized Normal Distribution Optimization Algorithm

4.1. Inspiration

The normal distribution theory serves as the foundation for GNDO. The normal distribution, commonly referred to as the Gaussian distribution, is a crucial tool for describing natural phenomena. The following is a definition of a normal distribution. Assume that random variable x follows a probability distribution with location μ and scale δ parameters, and that its probability density function may be written as:
f x = 1 2 π δ e x p x μ 2 2 δ 2
Following that, x can be referred to as a normal random variable, and this distribution can be referred to as a normal distribution. Two variables, the location parameter and scale parameter, are part of a normal distribution, according to Equation (53). It is possible to describe the mean value and standard deviation of random variables using the location parameter and scale parameter, respectively. In general, population-based optimization approaches’ search procedures consist of the three stages listed below. The scattered distribution contains all initialized people to start. Following that, everyone begins to move in the direction of the global optimal solution, and are guided by the designed exploration and exploitation tactics. The optimal answer is attained, and everyone congregates around it. Multiple normal distributions can adequately characterize this search process. To put it more precisely, individuals’ positions can be thought of as random variables with a normal distribution. The ideal position and the mean position are farther apart in the initial stage. The positions of all people exhibit a significantly high standard deviation. The gap between the average and ideal positions gradually narrows in the second stage. With each individual’s position, the standard variance decreases. The standard deviation of each individual’s location can be as low as possible in the final stage, which also sees the shortest distance between the mean position and the ideal position.

4.2. Local Exploitation

The suggested GNDO framework is depicted in Figure 3. As can be seen, GNDO has a fairly straightforward structure, and its local exploitation and global exploration information exchange mechanisms are built specifically for GNDO. The generalized normal distribution model that has been constructed—which is based on the current mean position and the present optimal position—is the foundation for local exploitation. Three people that were chosen at random were tied to global exploration. The following provides a thorough overview of the two learning strategies. Local exploitation is the process of locating better solutions within a search space made up of everyone’s present placements. A generalized normal distribution model for optimization can be constructed based on the correlation between the population’s distribution of people and the normal distribution:
v i t = μ i + δ i × η , i = 1 , 2 , 3 , , N
where v i t is the trailing vector of the ith individual at time t , μ i is the generalized mean position of the ith individual, δ i is generalized standard variance, and η is the penalty factor. Moreover, μ i ,   δ i , and η can be defined as
μ i = 1 3 x i t + x Best t + M
δ i = 1 3 x i t μ 2 + x Best t μ 2 + M μ 2
η = l o g λ 1 × c o s 2 π λ 2 ,   if   a < = b l o g λ 1 × c o s 2 π λ 2 + π ,   otherwise
where a , b, λ 1 , and λ 2 are random numbers between 0 and 1, x Best t is the present fitness location, and M is the mean location of the present population. In addition, M can be calculated by:
M = i = 1 N x i t N

4.3. Global Exploration

A speech space is searched globally to identify promising regions. The worldwide exploration of GNDO is based on three individuals who were chosen at random, as shown in Figure 7.
v i t = x i t + β × λ 3 × v 1 Local   information   sharing + 1 β × λ 4 × v 2 Global   information   sharing
where λ 3 and λ 4 are two random numbers subject to the standard normal distribution, β is called the adjust limit and is a random number between 0 and 1, and v 1 and v 2 are two trail vectors. Moreover, v 1 and v 2 can be calculated by:
v 1 = x i t x p 1 t , if f x i t < f x p 1 t x p 1 t x i t ,   otherwise
v 2 = x p 2 t x p 3 t , if f x p 2 t < f x p 3 t x p 3 t x p 2 t ,   otherwise
where p 1 , p 2 , and p 3 are three random integers selected from 1 to N, which meets p1 ≠ p2 ≠ p3 ≠ i. In the context of Equations (60) and (61), the second term to the right of Equation (59) can be referred to as the local learning term, indicating that solution p1 shares information with solution i , and the third term to the right of Equation (59) can be referred to as the global information sharing term, indicating that the information is provided to the individual i by the individuals p2 and p3. To balance the two information-sharing options, z utilizes the adjusted parameter β . Furthermore, because λ 3 and λ 4 are random numbers with a typical normal distribution, the search space for the GNDO can be expanded when conducting a global search. The absolute symbol in Equation (59) is used to stay steady with the screening mechanism in Equations (60) and (61).

4.4. The Implementation of the Proposed Method for Optimization

This section presents the GNDO implementation. The defined local exploitation and global exploration tactics serve as the foundation for the proposed GNDO. The two tactics are equally important and equally likely to be chosen for the GNDO. Additionally, similar to other population-based optimization methods, GNDO initializes its population by;
x i , j t = l j + u j l j × λ 5 , i = 1 , 2 , 3 , , N , j = 1 , 2 , 3 , , D
where D is the total number of design variables, l j is the jth design variable’s lower boundary, u j is its upper boundary, and λ 5 is a random number between 0 and 1. Note that neither a local exploitation approach nor a global exploration strategy will guarantee that the ith individual will find a better solution. A screening system is created to ensure that the population of the future generation receives the best solution, and it may be described as:
x i t + 1 = v i t ,                             if f ( v i t )   < f ( x i t ) x i t ,                                             otherwise
Figure 8 provides the pseudocode of the GNDO.

4.5. Basic Definitions of Multi-Objective Optimization

One optimization technique or tool that allows for more than one objective function to be used for any type of issue is multi-objective optimization. The following is a formulation of the fundamental elements of any multi-objective optimization:
Minimize : F x = f 1 x , f 2 x , , f o x Subjected   to : y i x = 0 ,       i = 1 , 2 , , n z i x 0 ,       i = 1 , 2 , , p l b i x i u b i ,       i = 1 , 2 , m
where y i signifies the ith equality constraint, z i denotes the ith inequality constraint, n signifies the number of equality constraints, p denotes the number of inequality constraints, x = x 1 , x 2 , , x j T denotes the decision/optimization variables, the minor and upper limits of the decision variables are represented by lb and ub, the number of decision variables is denoted by m, and the number of objective functions is denoted by o. In any multi-objective optimization problem, the relational operators are no longer valid to compare the search space solutions. The fundamental descriptions of multi-objective optimization problems are stated as follows, and a new operator known as Pareto optimality may be utilized to compare the solution.
Definition 1.
Pareto Optimality.
The solution x X is called Pareto optimum if and only if:
  y X   |   F y F x
Definition 2.
Pareto Dominance.
Let two different vectors be represented as x = x 1 , x 2 , , x j and y = y 1 , y 2 , , y j . The vector y is said to dominate the vector x (symbolized as y x ), if and only if:
i 1 , 2 , , j :   f i y f i x i 1 , 2 , , j : f i y < f i x
Definition 3.
Pareto Optimal Set.
All Pareto optimal solution sets are called the Pareto set, and are expressed as follows:
P s = { x , y X   |   F y F x }
Definition 4.
Pareto Optimal Front.
The Pareto optimal front is a collective of Pareto optimal solutions in the Pareto optimal set, as shown in (68):
P f = F x | x P s
Any multi-objective optimization issue must be solved using the Pareto optimum set during the multi-objective optimization process. The search space (set of dominated solutions) and objective space (set of non-dominated solutions) are depicted in Figure 9. The Pareto optimum front describes the interaction between the objective space and search space.

4.6. Multi-Objective Generalized Normal Distribution Optimization (MOGNDO)

The proposed MOGNDO algorithm optimizer uses both the crowding distance (CD) mechanism and the elitist non-dominated sorting (NDS) method. The NDS consists of the following stages:
  • Locating the non-dominated solution is the first step.
  • The second step is the use of the NDS strategy.
  • Performing non-dominated ranking (NDR) calculations on all non-dominated solutions.
Between two fronts, the NDR process takes place. The first front’s solutions offer a “0” index because no solutions are dominated by them, but at least one solution from the first front dominates the second front’s solutions. A solution’s NDR is equal to the number of solutions that predominate it. The CD process is used to keep the created solutions diverse. The following is a definition of the CD mechanism:
C D j i = f o b j j i + 1 f o b j j i 1 f o b j j m a x f o b j j m i n  
where f o b j j m a x and f o b j j m i n are the maximum and minimum values of j th objective function. The diagrammatic illustration of an NDS-based approach is illustrated in Figure 10.
The MOGNDO algorithm’s pseudocode is displayed in Algorithm 1. The MOGNDO method begins by specifying the necessary inputs, such as population size (Np), termination criteria, the maximum number of generations, and the maximum number of iterations (Maxit). Then, each objective function in the objective space vector F for Po is evaluated using a randomly generated parent population Po in the feasible search space region S. Thirdly, Po is subjected to the elitist-based CD and NDS. Fourthly, Po is merged with a fresh population of Pj to create a population, Pi. This Pi is sorted using the CD and NDR data, as well as elitist non-dominance. To establish a new parent population, the best Np options are evaluated. The process is then repeated until the termination criteria are met. MOGNDO’s flowchart is displayed in Figure 11.
Algorithm 1: Pseudocode of Multi-objective Generalized Normal Distribution Optimization (MOGNDO).
Step 1: Initially Generate population (Po) randomly in solution space (S)
Step 2: Evaluate objective space (F) for the generated population (Po)
Step 3: Sort the based on the elitist non-dominated sort method and find the non-dominated rank (NDR) and fronts
Step 4: Compute crowding distance (CD) for each front
Step 5: Update solutions (Pj)
Step 6: Merge Po and Pj to create Pi = Po U Pj
Step 7: For Pi perform Step 2
Step 8: Based on NDR and CD sort Pi
Step 9: Replace Po with Pi for Np first members of Pi

4.7. Constraint Handling Approach

The majority of engineering design issues in the actual world are multi-objective and highly nonlinearly constrained. To solve constrained MOPs, managing all constraints within their bounds is crucial. A static penalty technique is used in the MOGNDO algorithm because it transforms a constrained problem into an unconstrained problem, despite the literature survey giving various constrained handling approaches. This approach adds a significant penalty to the relevant goal function if a constraint is broken. The following is a presentation of the static penalty system:
f j X = f j X + i = 1 p P i max g i X , 0 + i = p N C P i max h i X δ , 0
where f j X , j = 1 , 2 n is the objective function to be optimized (here minimized), X = x 1 , x 2 , x m are design variables, g i X 0 , i = 1 , 2 p are inequality constraints, h i X = 0 , i = p + 1 N C are equality constraints, and δ is the tolerance in equality constraints.

4.8. Fuzzy Approach for the Multi-Objective Problem

The fuzzy membership approach can be used in multi-objective functions to identify the best compromising outcome out of all the non-inferior results. The fuzzy membership function μ f i uses a fuzzy membership function to keep track of the minimum f i m i n and maximum f i m a x values for each objective aim. Now, the membership function of the ith the objective is given as:
μ f i = 1                                                                                                 f i   f i m i n f i m a x f i f i m a x f i m i n                                           f i m i n < f i < f i m a x 0                                                                                                   f i   f i m a x
The standards of membership functions lie on the measure of (0–1) and display in the way that satisfies the function f i . Later, the decision-making function μ k should be calculated as follows:
μ k = i = 1 N μ f i k k = 1 M i = 1 N μ f i k
For non-inferior findings, the decision-making function can also be thought of as the normalized membership function, which displays the ordering of the undominated results. The end outcome is regarded as the best attainable compromise among all PFs, with a maximum value of   m a x i m u m   { μ k :   k = 1 , 2 , 3 M } .

5. Simulation Results, Analysis, and Comparative Study

This section discusses the outcomes of the MOGNDO algorithm, which optimized the optimal power flow with non-conventional and FACTS device problems with control variables. The initialization of the algorithm’s population size, archive size, the maximum number of iterations, and boundary condition for optimal power flow problems all came first. To identify the best optimal tradeoff points between multiple objective functions, the MOGNDO algorithm was then used to obtain the initial position and objective function values. Optimal power flow with non-conventional sources and FACTS devices were used to apply the MOGNDO algorithm’s performance, which was initially verified on eight unconstrained multi-objective problems. On a computer with 4 GB of RAM and a 3.20 GHz clock speed, the simulation was run using the MATLAB program. The benchmark functions for each unconstrained test were solved using 10 separate runs. The population size was set to 30, the maximum number of iterations was set to 100, and the archive size was set to 30 when the control parameters for the proposed MOGWO algorithm were first set. The performance measures for the MOGNDO algorithm, including Generational Distance (GD), Inversion Generational Distance (IGD), Spacing Metrics (SP), Diversity Metrics (DM), and Spread Metrics (SD), are covered in this section.

5.1. MOGNDO Results for Test Benchmark Problems

Before tackling real-world issues, the MOGNDO was used to evaluate the performance of the benchmark unconstraint test function provided in [26]. Eight benchmark unconstrained test functions—ZDT1, ZDT2, ZDT3, ZDT4, ZDT6, KURSAVE, SCHAFFER-1, and SCHAFFER-2 (Figure 12) were taken into account, and a thorough simulation was performed using the MOGNDO technique. Any algorithm’s control parameters are crucial to the resolution of the optimization problem. As a result, the number of populations was decided after conducting a comparative analysis that took into account various population sizes, while holding all other variables constant. Following careful consideration, the population size, maximum iterations, and archive size were chosen as 30, 100, and 30, respectively, for the unconstrained test benchmark functions. The MOGNDO algorithm’s performance was evaluated using performance metrics, such as Generational Distance (GD), Inversion Generational Distance (IGD), Spacing Metrics (SP), Diversity Metrics (DM), and Spread Metrics (SD), for convergence measurement. Table 3, Table 4, Table 5, Table 6 and Table 7 demonstrate that MOGNDO could achieve the best outcomes for all performance metrics, including Generational Distance (GD), Inversion Generational Distance (IGD), Spacing Metrics (SP), Diversity Metrics (DM), and Spread Metrics (SD), which cover convergence and solution accuracy. It follows that the suggested MOGNDO can provide the best convergence on all benchmark functions. The outcomes (archive solutions) of all eight test benchmark issues are displayed in Figure 1, Figure 2, Figure 3, Figure 4 and Figure 5. As can be shown, the MOGNDO method was capable of approximating the PF. By comparing the PF estimations, it can also be seen that the suggested MOGNDO could provide acceptable performance. Thus, it was determined that the MOGNDO algorithm is more suitable for the stochastic OPF problem with three FACTS devices and wind power plants.

5.2. Multi-Objectives OPF Problem with Wind Power Plants and Three FACTS Devices

The GNDO algorithm was used to address the stochastic OPF problem with wind power plants and three FACTS devices in this study. The solution to the optimum power flow problem was evaluated in parallel using newly created algorithms, such as the Multi-Verse Optimization (MVO), the Sine-Cosine Algorithm (SCA) [27], the Grey Wolf Optimization (GWO), the Moth Fame Optimization (MFO), the Ant Lion Optimization (ALO) [28], and Ion Motion Algorithms (IMA) [29]. The proposed approach was demonstrated using a modified IEEE-30 bus infrastructure with wind power plants and FACTS devices. Table 1 lists the major characteristics of the customized IEEE-30 bus framework. The following are two scenarios:
  • Scenario-1 (Solo objective OPF with wind power plants and FACTS devices)
  • Scenario-2 (Multi-objective OPF with wind power plants and FACTS devices)
As shown in Table 8, there were a total of thirteen different test scenarios to evaluate. In this section, the results of case studies using various metaheuristics methodologies are tabulated and presented. The first six case studies are for single-objective optimization, while the latter seven are multi-objective optimization problems that include non-conventional sources of energy resources, as well as optimal FACTS device sizes and locations. The search agent value was set to 40, and each algorithm underwent 500 iterations of analysis. Please refer to the original research for a detailed discussion of those procedures. Table 4 shows the parameter settings for these methods.

5.3. Scenario-1 (Single Objective OPF with Wind Power Plants and FACTS Devices)

With the use of GNDO, MVO, ALO, SCA, and IMO methods, all of the objective goals indicated in the mathematical formulation were simultaneously handled as solo objective optimization issues. The limitations of all control variables, as well as proper FACTS device locations and sizing, are listed below. From case 1 to case 6, the outcomes of objective functions are tabulated in Table 9, Table 10 and Table 11, with the best minimum values containing five different recent techniques.
The overall fuel cost with GNDO, which included the two non-conventional sources of power plants and optimal placement of FACTS devices, was 806.999 $/h, which was the best in comparison with the other cited algorithm shown in Table 9. The reductions in TFC in comparison with MVO, ALO, SCA, IMO, SHADE-SF, DE-SF, ABC-SF, PSO-SF, FPA-SF, and MSA-SF were 1.031 $/h, 1.031 $/h, 2.45 $/h, 11.655 $/h 7.866 $/h, 0.0176 $/h, 0.4917 $/h, 0.4 $/h, 1.2553 $/h, 3.398 $/h, and 1.0403 $/h, respectively. This demonstrated the GNDO algorithm’s superiority over other cited metaheuristics algorithms.
Figure 13 illustrates the convergence traits of the TFC minimization. Similar convergence traits for APL, voltage deviations, and VSI are depicted in Figure 14, Figure 15 and Figure 16. Figure 17 also displays a comparison of the fuel cost decrease with various algorithms. In example 2, the GNDO method resulted in a pollutant gas emission of 0.138 tons per hour. In instance 3, the APL of the various transmission lines using the GNDO approach was 1.647 MW. The APL was 0.088 MW, 0.039 MW, 0.835 MW, 0.233 MW, 0.0997 MW, 0.0997 MW, 0.2598 MW, 0.2494 MW, 0.6127 MW, and 0.4972 MW less compared to MVO, ALO, SCA, IMO, SHADE-SF, DE-SF, ABC-SF, PSO-SF, and MSA-SF, respectively. A crucial factor for the grid’s ability to operate reliably was the voltage divergence of each bus from 1.0 per unit. Therefore, in scenario 4, the moth flame algorithm produced the lowest voltage variation (0.124 p.u), making it the best of the five optimization methods. The VSI, sometimes referred to as the L max index, varied between zero (no load) and one (voltage collapse). Therefore, in case 5, 0.096 was the lowest value for the L max index. In scenario 6, the overall gross fuel cost using the GNDO method was 1120.996 dollars per hour. It is interesting to note here that in Table 12, the total gross fuel cost of the proposed GNDO was more than the SHADE-SF, which further enforces the narrative of the “No free lunch theorem,” which states that no algorithm gives the best result in every problem.
Figure 18 and Figure 19 provide comparison charts with the minimizing of APL and overall fuel cost. The results of the simulations were compared to those of the most recent algorithms, including MVO, ALO, SCA, IMO, and other mentioned optimization approaches. It was found that the proposed method of Generalized Normal Distribution Optimization methodology produced superior results.

5.4. Scenario-2 (Multi-Objective OPF with Non-Conventional Sources Energy Resources)

In this case, two, three, and four objectives were simultaneously optimized using the Multi-Objective Generalized Normal Distribution Optimization (MOGNDO) algorithm technique. In multi-objective optimization, the non-dominated sorting optimization technique is used to simultaneously find solutions for numerous objectives. To find the PF for the modified IEEE 30-bus architecture, 30 non-dominate solutions are retained. The scenarios in cases 7 through 10 are thought of as two-objective optimization cases. Three-objective optimization problems are what are known as cases 11 and 12. Contrarily, case 13 is referred to as a set of four problems involving objective optimization. Among all the Pareto archives, the best compromising solution was found using the fuzzy decision-making method. For cases 7 through 13, the most optimal compromise solutions using the proposed MOGNDO algorithm and other cited metaheuristics techniques are shown in boldface and stated in Table 13, Table 14, Table 15 and Table 16. The best PFs of TFC and pollution minimization for case 7 are shown in Figure 20, utilizing various metaheuristics techniques. Likewise with case 7, case 8’s PF with two goal optimizations resulted in APL and TFC, which are depicted in Figure 21. Figure 22 shows the PF of the TFC with the carbon tax and voltage deviation minimization. In scenario 11, Figure 23 shows the three objectives PFs for minimizing APL, TFC, and toxic gas emanations. Figure 24 depicts the PF for the minimization of voltage variation, APL, and TFC, using various algorithms.
The MOGNDO methodology was one of the finest methods for finding the best solutions to the multi-objective OPF problem that integrated with wind power plants and the appropriate placement of FACTS devices, according to the tabulated data.

6. Conclusions

The optimal location and size of FACTS devices in this study, as well as single- and multi-objective optimal power flow (MOOPF) concerned coal-based and wind power plants, which were all addressed by the solution technique. Different probability density functions were used to express inconsistencies in unconventional resource availability. The method for integrating each unit was described in detail. When utilizing non-conventional sources of energy sources and FACTS devices, single objectives were optimized, such as generation cost, toxic gas emanation, voltage deviation, active loss, and VSI. A multi-objective form of the OPF problem was looked into in light of the current situation of the electric network. The outcomes were contrasted with a recently created optimization strategy. Based on the results, it can be said that the proposed MOGNDO outperformed existing algorithms in terms of convergence, and delivered higher quality and more usable solutions for each situation involving optimal power flow. All of the results point to the suggested technique’s significant advantage in obtaining the best solutions to OPF issues with one or more objectives. Finally, it was shown that by integrating wind farms with FACTS devices utilizing a non-dominated sorting technique, MOGNDO could be successfully employed to address small and large optimal power flow challenges. Based on the extensive analysis of the proposed MOGNDO, the following can be summarized as its advantages—
  • Randomization in MOGNDO includes the diversity of the Pareto front being enhanced, since all solutions in the first dominated front will have an equal chance of being selected, and multi-objectives are made uniformly significant while performing local exploration.
  • MOGNDO can deal with large-scale search spaces and is less dependent on problem characteristics. Moreover, these algorithms are capable of estimating multiple points in the search domain simultaneously, due to their population-based nature.
  • MOGNDO strikes a good balance between exploitation and exploration, providing powerful searchability for finding the optimum solution
  • MOGNDO is superior in terms of the balance of diversity and convergence, the distribution of PF, and better convergence.

Author Contributions

Conceptualization, S.B.P., J.V. and M.M.; Methodology, S.B.P., J.V., M.M. and T.K.M.; Software, T.K.M. and P.J.; Validation, S.B.P., M.M. and T.K.M.; Formal analysis, J.V., M.M. and T.K.M.; Investigation, S.B.P.; Resources, P.J.; Data curation, J.V. and P.J.; Writing—original draft, S.B.P., J.V., M.M., T.K.M. and P.J.; Writing—review & editing, S.B.P.; Visualization, T.K.M. and P.J.; Supervision, S.B.P.; Funding acquisition, M.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the project SP2022/60 Applied Research in the Area of Machines and Process Control, which is supported by the Ministry of Education, Youth and Sports, in the Czech Republic.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in the article.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

Acronyms
OPFOptimal Power Flow
MasMeta Heuristics Algorithms
MOGNDOMulti-Objective Generalized Normal Distribution Optimization
TG Thermal Generating unit
WGWind Generation
ISOIndependent System Operator
PDFProbability Density Function
BCS Best Compromise Solution
MOMFO Multi-Objective Moth Flame Optimization
MOOPFMulti-Objective Optimal Power Flow
SHADE-SFSuccess History-based Adaptive Differential Evolution using Superiority of Feasible solutions method
DE-SFDifferential Evolution using Superiority of Feasible solutions method
ABC-SFArtificial Bee Colony using Superiority of Feasible solutions method
PSO-SFParticle Swarm Optimization using Superiority of Feasible solutions method
FPA-SFFlower Pollination Algorithm using Superiority of Feasible solutions method
MSA-SFMoth Swarm Algorithm using Superiority of Feasible solutions method
TFCTotal Fuel Cost
APLActive power Loss
VSIVoltage Stability Index
Nomenclature
a i ,   b i , c i , e i   and   d i Price constants for ith coal-based power plants.
α i , β i ,   γ i , ω i   and   μ i Toxic gas emanation constants concerning the ith coal-based units.
g w Direct cost constant
P w s Scheduled power of the wind unit.
K R w Reserve cost coefficient regarding wind unit
K P w Penalty cost coefficient of wind unit
P w s Accessible power from the wind unit
P w r Specified output power from the wind unit
f w p w Wind energy probability density function for the wind unit.
v i n , v r and v o u t Cut-in, rated, and cut-out wind velocity of the turbine respectively
p w r Rated value of the generated output of the wind turbine
τ Degree of series compensation
X m n Line inductive reactance linking buses m and n
R m n Resistance of the line linking buses m and n
V m and V n Bus voltage magnitudes linking buses m and n .
δ m and δ n Phase angles of the linking buses m and n
g m n and b m n Conductance and susceptance of the line linking buses m and n .
N p q Number of load (PQ) buses
v i pu voltage level of i th bus.
P G i and P D i Generation and dispatch at i th bus
Number of buses
Y 1 and Y 2 Sub-matrices of
δ i j = δ i δ j Variance in phase angles of voltage among bus i and bus
P D i and Q D i Real and VAR power demand respectively at i th bus
P G i and Q G i Real and VAR outputs respectively of i th bus by either unit (coal-based or non-conventional) as applicable
G i j and B i j Conductance and susceptance between bus i and bus j

References

  1. Joshi, M.; Ghadai, R.K.; Madhu, S.; Kalita, K.; Gao, X.-Z. Comparison of NSGA-II, MOALO and MODA for Multi-Objective Optimization of Micro-Machining Processes. Materials 2021, 14, 5109. [Google Scholar] [CrossRef] [PubMed]
  2. Cao, B.; Li, M.; Liu, X.; Zhao, J.; Cao, W.; Lv, Z. Many-objective deployment optimization for a drone-assisted camera network. IEEE Trans. Netw. Sci. Eng. 2021, 8, 2756–2764. [Google Scholar] [CrossRef]
  3. Liang, X.; Luo, L.; Hu, S.; Li, Y. Mapping the knowledge frontiers and evolution of decision making based on agent-based modeling. Knowl.-Based Syst. 2022, 250, 108982. [Google Scholar] [CrossRef]
  4. Zhang, K.; Wang, Z.; Chen, G.; Zhang, L.; Yang, Y.; Yao, C.; Wang, J.; Yao, J. Training effective deep reinforcement learning agents for real-time life-cycle production optimization. J. Pet. Sci. Eng. 2022, 208, 109766. [Google Scholar] [CrossRef]
  5. Ganesh, N.; Ghadai, R.K.; Bhoi, A.K.; Kalita, K.; Gao, X.-Z. An intelligent predictive model-based multi-response optimization of EDM process. Comput. Model. Eng. Sci. 2020, 124, 459–476. [Google Scholar] [CrossRef]
  6. Ghadai, R.K.; Kalita, K.; Gao, X.-Z. Symbolic regression metamodel based multi-response optimization of EDM process. FME Trans. 2020, 48, 404–410. [Google Scholar] [CrossRef]
  7. Cao, B.; Zhang, W.; Wang, X.; Zhao, J.; Gu, Y.; Zhang, Y. A memetic algorithm based on two_Arch2 for multi-depot heterogeneous-vehicle capacitated arc routing problem. Swarm Evol. Comput. 2021, 63, 100864. [Google Scholar] [CrossRef]
  8. Liu, Y.; Zhang, Z.; Liu, X.; Wang, L.; Xia, X. Ore image classification based on small deep learning model: Evaluation and optimization of model depth, model structure and data size. Miner. Eng. 2021, 172, 107020. [Google Scholar] [CrossRef]
  9. Ullah, Z.; Wang, S.; Radosavljevic, J.; Lai, J. A solution to the optimal power flow problem considering WT and PV generation. IEEE Access 2019, 7, 46763–46772. [Google Scholar] [CrossRef]
  10. Elattar, E.E. Optimal power flow of a power system incorporating stochastic wind power based on modified moth swarm algorithm. IEEE Access 2019, 7, 89581–89593. [Google Scholar] [CrossRef]
  11. Man-Im, A.; Ongsakul, W.; Singh, J.G.; Madhu, M.N. Multi-objective optimal power flow considering wind power cost functions using enhanced PSO with chaotic mutation and stochastic weights. Electr. Eng. 2019, 101, 699–718. [Google Scholar] [CrossRef]
  12. Salkuti, S.R. Optimal power flow using multi-objective glowworm swarm optimization algorithm in a wind energy integrated power system. Int. J. Green Energy 2019, 16, 1547–1561. [Google Scholar] [CrossRef]
  13. Kathiravan, R.; Kumudini Devi, R.P. Optimal power flow model incorporating wind, solar, and bundled solar-thermal power in the restructured Indian power system. Int. J. Green Energy 2017, 14, 934–950. [Google Scholar] [CrossRef]
  14. Duman, S.; Rivera, S.; Li, J.; Wu, L. Optimal power flow of power systems with controllable wind-photovoltaic energy systems via differential evolutionary particle swarm optimization. Int. Trans. Electr. Energy Syst. 2020, 30. [Google Scholar] [CrossRef]
  15. Duman, S.; Li, J.; Wu, L.; Guvenc, U. Optimal power flow with stochastic wind power and FACTS devices: A modified hybrid PSOGSA with chaotic maps approach. Neural Comput. Appl. 2020, 32, 8463–8492. [Google Scholar] [CrossRef]
  16. Biswas, P.P.; Suganthan, P.N.; Qu, B.Y.; Amaratunga, G.A.J. Multiobjective economic-environmental power dispatch with stochastic wind-solar-small hydro power. Energy 2018, 150, 1039–1057. [Google Scholar] [CrossRef]
  17. Chen, M.-R.; Zeng, G.-Q.; Lu, K.-D. Constrained multi-objective population extremal optimization based economic-emission dispatch incorporating renewable energy resources. Renew. Energy 2019, 143, 277–294. [Google Scholar] [CrossRef]
  18. Chang, Y.-C.; Lee, T.-Y.; Chen, C.-L.; Jan, R.-M. Optimal power flow of a wind-thermal generation system. Int. J. Electr. Power Energy Syst. 2014, 55, 312–320. [Google Scholar] [CrossRef]
  19. Saha, A.; Bhattacharya, A.; Das, P.; Chakraborty, A.K. A novel approach towards uncertainty modeling in multiobjective optimal power flow with renewable integration. Int. Trans. Electr. Energy Syst. 2019, 29. [Google Scholar] [CrossRef]
  20. Biswas, P.P.; Suganthan, P.N.; Amaratunga, G.A.J. Optimal power flow solutions incorporating stochastic wind and solar power. Energy Convers. Manag. 2017, 148, 1194–1207. [Google Scholar] [CrossRef]
  21. Ben Hmida, J.; Chambers, T.; Lee, J. Solving constrained optimal power flow with renewables using hybrid modified imperialist competitive algorithm and sequential quadratic programming. Electr. Power Syst. Res. 2019, 177, 105989. [Google Scholar] [CrossRef]
  22. Pandya, S.; Jariwala, H.R. Single- and multiobjective optimal power flow with stochastic wind and solar power plants using moth flame optimization algorithm. Smart Sci. 2022, 10, 77–117. [Google Scholar] [CrossRef]
  23. Biswas, P.P.; Arora, P.; Mallipeddi, R.; Suganthan, P.N.; Panigrahi, B.K. Optimal placement and sizing of FACTS devices for optimal power flow in a wind power integrated electrical network. Neural Comput. Appl. 2021, 33, 6753–6774. [Google Scholar] [CrossRef]
  24. Wolpert, D.H.; Macready, W.G. No free lunch theorems for optimization. IEEE Trans. Evol. Comput. 1997, 1, 67–82. [Google Scholar] [CrossRef] [Green Version]
  25. Zhang, Y.; Jin, Z.; Mirjalili, S. Generalized normal distribution optimization and its applications in parameter extraction of photovoltaic models. Energy Convers. Manag. 2020, 224, 113301. [Google Scholar] [CrossRef]
  26. Kumar, S.; Jangir, P.; Tejani, G.G.; Premkumar, M.; Alhelou, H.H. MOPGO: A new physics-based multi-objective plasma generation optimizer for solving structural optimization problems. IEEE Access 2021, 9, 84982–85016. [Google Scholar] [CrossRef]
  27. Mirjalili, S. SCA: A Sine Cosine Algorithm for solving optimization problems. Knowl. Based Syst. 2016, 96, 120–133. [Google Scholar] [CrossRef]
  28. Mirjalili, S. The ant lion optimizer. Adv. Eng. Softw. 2015, 83, 80–98. [Google Scholar] [CrossRef]
  29. Javidy, B.; Hatamlou, A.; Mirjalili, S. Ions motion algorithm for solving optimization problems. Appl. Soft Comput. 2015, 32, 72–79. [Google Scholar] [CrossRef]
Figure 1. Adapted IEEE 30-bus scheme with power units and FACTS policies [23].
Figure 1. Adapted IEEE 30-bus scheme with power units and FACTS policies [23].
Electronics 11 03825 g001
Figure 2. Weibull PDF (bus 5).
Figure 2. Weibull PDF (bus 5).
Electronics 11 03825 g002
Figure 3. Weibull PDF (bus 11).
Figure 3. Weibull PDF (bus 11).
Electronics 11 03825 g003
Figure 4. Basic Structure and model of TCSC [23].
Figure 4. Basic Structure and model of TCSC [23].
Electronics 11 03825 g004
Figure 5. Model of TCP [23].
Figure 5. Model of TCP [23].
Electronics 11 03825 g005
Figure 6. Basic structure and model of SVC [23].
Figure 6. Basic structure and model of SVC [23].
Electronics 11 03825 g006
Figure 7. The framework of the proposed GNDO [25].
Figure 7. The framework of the proposed GNDO [25].
Electronics 11 03825 g007
Figure 8. The pseudocode of GNDO [25].
Figure 8. The pseudocode of GNDO [25].
Electronics 11 03825 g008
Figure 9. Objective space and search space in the multi-objective optimization problem.
Figure 9. Objective space and search space in the multi-objective optimization problem.
Electronics 11 03825 g009
Figure 10. The procedure of the non-dominated sorting approach.
Figure 10. The procedure of the non-dominated sorting approach.
Electronics 11 03825 g010
Figure 11. Flowchart of MOGNDO algorithm.
Figure 11. Flowchart of MOGNDO algorithm.
Electronics 11 03825 g011
Figure 12. Best Pareto optimal front obtained of unconstrained test functions by MOGNDO algorithm.
Figure 12. Best Pareto optimal front obtained of unconstrained test functions by MOGNDO algorithm.
Electronics 11 03825 g012
Figure 13. Convergence characteristics of TFC minimization.
Figure 13. Convergence characteristics of TFC minimization.
Electronics 11 03825 g013
Figure 14. Convergence characteristics of APL minimization.
Figure 14. Convergence characteristics of APL minimization.
Electronics 11 03825 g014
Figure 15. Convergence characteristics of voltage deviation minimization.
Figure 15. Convergence characteristics of voltage deviation minimization.
Electronics 11 03825 g015
Figure 16. Convergence characteristics of VSI minimization.
Figure 16. Convergence characteristics of VSI minimization.
Electronics 11 03825 g016
Figure 17. Comparison of the TFC reduction ($/h) with other algorithms.
Figure 17. Comparison of the TFC reduction ($/h) with other algorithms.
Electronics 11 03825 g017
Figure 18. Comparison chart for TFC minimization with different algorithms.
Figure 18. Comparison chart for TFC minimization with different algorithms.
Electronics 11 03825 g018
Figure 19. Comparison chart for APL minimization with different algorithms.
Figure 19. Comparison chart for APL minimization with different algorithms.
Electronics 11 03825 g019
Figure 20. PF of TFC and emission minimization.
Figure 20. PF of TFC and emission minimization.
Electronics 11 03825 g020
Figure 21. PF of TFC and APL minimization.
Figure 21. PF of TFC and APL minimization.
Electronics 11 03825 g021
Figure 22. PF of TFC with a carbon tax and voltage deviation minimization.
Figure 22. PF of TFC with a carbon tax and voltage deviation minimization.
Electronics 11 03825 g022
Figure 23. PF of TFC, emission, and APL minimization with different algorithms.
Figure 23. PF of TFC, emission, and APL minimization with different algorithms.
Electronics 11 03825 g023
Figure 24. PF of TFC, APL, and voltage deviation minimization with different algorithms.
Figure 24. PF of TFC, APL, and voltage deviation minimization with different algorithms.
Electronics 11 03825 g024
Table 1. Test system key features for analysis.
Table 1. Test system key features for analysis.
ParticularsQuantityDetails
Total buses30[23]
Total branches 41[23]
Coal-based generators (TG1; TG2; TG3; TG4)4Buses: 1 (swing), 2, 8 and 13
Wind generators (WG1; WG2)2Bus-5 and Bus-11
Tap changing transformers4Branches: 11, 12, 15 and 36
SVC2Optimal bus and rating derived
TCSC2Optimal branch position and rating derived
TCPS2Optimal placement and rating derived
Demand-283.4 MW, 126.2 MVAr
Table 2. PDF constants of wind power plants [23].
Table 2. PDF constants of wind power plants [23].
WindfarmNo. of TurbinesRated PowerWeibull PDF ParametersCost Constants
DirectReservePenalty
WG5 (bus 5)2575 c = 9, k = 21.603.01.50
WG11 (bus 11)2060c = 10, k = 21.753.01.50
Table 3. Results of G D M E T R I C S on test functions.
Table 3. Results of G D M E T R I C S on test functions.
TEST FUNCTIONSMinimumAverageMedianMaximumStd Dev
ZDT-10.000147950.000250510.000250430.000331336.658 × 10−5
ZDT-20.000152930.000166650.000169150.000182718.9928 × 10−6
ZDT-30.000483140.000596010.000579430.000773549.9396 × 10−5
ZDT-40.000121040.000228080.000208340.00047589.4401 × 10−5
ZDT-69.084 × 10−50.0469630.000113160.227770.084189
KURSAVE0.000797410.00144520.00129050.00252110.00051722
SCHAFFER-10.000336070.000409920.000422030.000478974.714 × 10−5
SCHAFFER-24.3236 × 10−57.2824 × 10−55.2463 × 10−50.000169394.0743 × 10−5
Table 4. Results of I G D   M E T R I C S on test functions.
Table 4. Results of I G D   M E T R I C S on test functions.
TEST FUNCTIONSMinimumAverageMedianMaximumStd Dev
ZDT-10.000866230.000988620.00103480.00107188.6005 × 10−5
ZDT-20.000865780.000988480.000948370.00123580.00012678
ZDT-30.0012110.00207460.00133860.00772050.0020015
ZDT-40.000780660.000895180.000870080.00122760.00012621
ZDT-60.00042070.000455250.000446880.000497482.8268 × 10−5
KURSAVE0.000544360.000623990.000578450.000854089.8135 × 10−5
SCHAFFER-10.00131540.00152460.00150270.00170640.0001302
SCHAFFER-20.000388080.000425940.000425860.000449331.8173 × 10−5
Table 5. Results of S P A C I N G M E T R I C S on test functions.
Table 5. Results of S P A C I N G M E T R I C S on test functions.
TEST FUNCTIONSMinimumAverageMedianMaximumStd Dev
ZDT-10.0555130.0689220.0687730.0780860.0087192
ZDT-20.0552690.0669820.062250.0827060.010406
ZDT-30.237080.301990.304890.371950.038816
ZDT-40.0554610.0796730.0827980.09840.012829
ZDT-60.0602750.364670.0826721.29840.50689
KURSAVE1.71232.08172.05412.38070.20644
SCHAFFER-10.477340.64440.657220.784760.086121
SCHAFFER-24.31846.2286.14588.17970.99173
Table 6. Results of D I V E R S I T Y M E T R I C S on test functions.
Table 6. Results of D I V E R S I T Y M E T R I C S on test functions.
TEST FUNCTIONSMinimumAverageMedianMaximumStd Dev
ZDT-10.397740.47260.469590.537720.060856
ZDT-20.346660.43690.442150.571390.069382
ZDT-30.424020.521220.506250.648130.07401
ZDT-40.30350.375730.368220.430950.043003
ZDT-60.361480.69070.452361.34550.42847
KURSAVE0.287740.38390.390090.457630.050257
SCHAFFER-10.287740.38390.390090.457630.050257
SCHAFFER-20.91610.956070.959921.00330.030768
Table 7. Results of S P R E A D   M E T R I C S on test functions.
Table 7. Results of S P R E A D   M E T R I C S on test functions.
TEST FUNCTIONSMinimumAverageMedianMaximumStd Dev
ZDT-10.386490.457940.450490.52510.059019
ZDT-20.335110.424590.428420.567090.070091
ZDT-30.56040.645030.647490.738140.066426
ZDT-40.303640.367410.362040.42130.040601
ZDT-60.352910.679630.44791.32930.42341
KURSAVE0.377290.434910.434480.48720.030263
SCHAFFER-10.271130.373750.378320.444410.051198
SCHAFFER-20.564310.612930.615610.650220.02504
Table 8. Summary of case studies for adapted IEEE-30 bus test system.
Table 8. Summary of case studies for adapted IEEE-30 bus test system.
Test SystemCase #Single and Multi-Objectives Functions
Adapted
IEEE 30-bus test system
Case # 1TFC includes FACTS devices, wind farms, and coal-based plants
Case # 2Reduction of total toxic gas emanations with the use of coal-based, wind, and FACTS technologies.
Case # 3Minimization of APL in FACTS devices, wind farms, and coal-based plants.
Case # 4Minimization of the total voltage variation using coal-based, wind, and FACTS devices.
Case # 5Voltage stability improvement in coal-based, wind, and FACTS equipment.
Case # 6Total Gross Generation Cost, includes FACTS devices, wind farms, and coal-based plants.
Case # 7Minimizing TFCs and toxic gas emanations while using non-conventional sources of energy and FACTS devices
Case # 8TFC and APL Minimization Including non-conventional sources and FACTS Devices
Case # 9TFC and VSI minimization including non-conventional sources and FACTS devices
Case # 10Total Gross Generation Cost and voltage deviation minimization with non-conventional sources and FACTS devices
Case # 11TFC, Toxic gas emanation, and APL minimization together with non-conventional sources and FACTS devices
Case # 12TFC, APL, and VSI minimization including non-conventional sources and FACTS devices
Case # 13TFC, Toxic gas emanation, APL, and voltage deviation minimization including non-conventional sources and FACTS devices
Table 9. Single objectives simulation results for case 1 and case 2.
Table 9. Single objectives simulation results for case 1 and case 2.
Control &
State Variables
MinMaxCase-1Case-2
GNDOMVOALOSCAIMOGNDOMVOALOSCAIMO
PTG2208041.42740.31140.96035.45830.88146.63446.63946.63448.35746.726
PWG507549.81549.45949.07741.71954.22674.81874.93471.36275.00074.507
PTG8103510.00010.35213.03815.17214.36135.00035.00035.00035.00035.000
PWG1106040.79942.13539.36247.56136.32952.36551.21554.90560.00048.761
PTG13124012.00212.00012.00014.32116.64840.00040.00040.00040.00040.000
V10.951.11.0911.1001.1001.0191.1001.0900.9971.1000.9971.100
V20.951.11.0751.0901.0910.9921.1001.0801.0131.0780.9501.100
V50.951.11.0491.0711.0720.9781.1001.0671.0370.9751.1001.100
V80.951.11.0461.0761.0790.9501.1000.9621.1001.1001.1001.100
V110.951.11.1001.1001.0711.0471.1001.0731.0821.1000.9501.100
V130.951.11.0691.0621.0380.9501.1001.0870.9701.0471.1001.100
T110.91.11.0371.0071.0590.9591.0900.9770.9381.0201.1001.100
T120.91.10.9971.0871.0840.9001.0900.9431.0551.0730.9651.100
T150.91.11.0271.0991.0950.9381.0901.0231.0631.0800.9001.100
T360.91.10.9571.0241.0870.9001.0901.0171.0241.0540.9001.100
SVC1 Location--2427619222711212529
SVC2 Location--72710915620302530
SVC1 Rating−101010.0009.893−6.5413.8821.740−3.3809.2748.78610.00010.000
SVC2 Rating−10105.764−6.2771.2095.482−3.094−9.5590.6398.550−10.0004.480
TCSC1 Location--1539103143738323433
TCSC2 Location--1229204323624374139
TCSC1 Rating00.50.4910.2190.2180.0000.4520.1100.1960.4970.0000.500
TCSC2 Rating00.50.4960.2390.4400.0000.4920.4540.0600.4850.1930.500
TCPS1 Location--14163415143814264040
TCPS2 Location--16221913035532141
TCPS1 Rating−552.7144.503−3.9903.9171.2354.787−0.1750.9265.0005.000
TCPS2 Rating−551.7053.761−4.2811.1500.1134.0951.1484.860−3.9634.566
TFC ($/h) 806.999808.030809.449818.654814.865-----
Emission (Ton/h) -----0.1380.1380.1380.1380.138
Table 10. Single objectives simulation results for case 3 and case 4.
Table 10. Single objectives simulation results for case 3 and case 4.
Control &
State Variables
MinMaxCase-3Case-4
GNDOMVOALOSCAIMOGNDOMVOALOSCAIMO
PTG2208069.00575.90079.74580.00079.64479.72178.27128.57749.60634.817
PWG507575.00074.96675.00073.82574.66747.50411.23613.4190.00047.506
PTG8103534.99934.65635.00030.72634.84425.35732.73616.32835.00034.441
PWG1106059.99958.42260.00050.63559.74826.23420.2924.6020.00055.391
PTG13124039.97631.23939.82331.70639.83226.91632.99624.66924.06039.576
V10.951.11.0361.1001.0981.0981.0891.0080.9580.9670.9500.955
V20.951.11.0371.1001.1001.0911.0891.0321.0541.0511.0721.052
V50.951.11.0271.0901.0911.1001.0891.0101.0161.0170.9600.955
V80.951.11.0301.0921.0961.1001.0901.0250.9921.0131.0241.015
V110.951.11.1001.0991.1001.1001.0890.9501.0691.0091.1001.004
V130.951.11.1001.1001.0770.9901.0891.0041.0671.0841.0431.038
T110.91.11.0121.0471.0461.0301.0900.9381.0630.9491.0170.905
T120.91.10.9030.9031.0451.1001.0900.9070.9030.9040.9120.986
T150.91.10.9981.0421.1001.0181.0900.9451.0481.0651.0441.069
T360.91.10.9350.9831.0681.1001.0900.9340.9300.9360.9570.936
SVC1 Location--18122711271924281919
SVC2 Location--24153016301014162123
SVC1 Rating−10104.9016.8713.9516.5359.9568.0828.769−1.8984.9699.593
SVC2 Rating−101010.0005.4485.238−0.9364.2195.421−2.342−5.291−1.4659.547
TCSC1 Location--3454034014132240
TCSC2 Location--11404124118259529
TCSC1 Rating00.50.4940.1070.5000.0000.4710.3530.3000.0750.0270.500
TCSC2 Rating00.50.5000.1980.4970.0000.4980.4990.3910.0870.0000.486
TCPS1 Location--16124043319383137
TCPS2 Location--19154113415415536
TCPS1 Rating−551.744−0.898−3.3780.4334.979−4.9973.754−1.6385.0002.660
TCPS2 Rating−550.2574.748−3.4411.2490.9460.7670.550−3.8740.128−1.703
APL (MW) 1.6471.7351.6862.4821.880-----
Voltage Deviation (p.u) -----0.1240.1500.1770.2270.165
Table 11. Single objectives simulation results for case 5 and case 6.
Table 11. Single objectives simulation results for case 5 and case 6.
Control &
State Variables
MinMaxCase-5Case-6
GNDOMVOALOSCAIMOGNDOMVOALOSCAIMO
PTG2208078.16128.76576.48720.00074.36344.89447.80255.20320.00056.164
PWG507575.00016.02074.2400.00067.67774.99874.56471.25675.00068.555
PTG8103535.00034.85534.40510.00032.81535.00032.89333.27135.00030.826
PWG1106054.3550.00055.98349.37914.17359.00658.03050.10260.00058.732
PTG13124012.00128.84038.08016.65935.96121.58322.22232.27419.54023.229
V10.951.11.1001.1001.1001.1001.0971.0461.1001.1001.1001.100
V20.951.11.1001.1001.1001.1001.0971.0421.0971.0981.1001.100
V50.951.11.1001.1001.1001.1001.0971.0321.0871.0871.1001.100
V80.951.11.1001.1001.1001.1001.0971.0351.0921.0911.1001.098
V110.951.11.1001.1001.1001.1001.0971.0981.1001.1001.1001.100
V130.951.11.1001.1001.1001.1001.0971.0351.1001.0831.1001.098
T110.91.10.9050.9100.9901.1001.0891.0771.0541.0071.1001.084
T120.91.10.9040.9090.9900.9001.0890.9010.9031.0891.1001.084
T150.91.10.9010.9000.9300.9001.0891.0651.0531.0811.1001.100
T360.91.10.9010.9090.9100.9000.9100.9851.0011.0331.1001.087
SVC1 Location--10292929262428261524
SVC2 Location--2930303026131426310
SVC1 Rating−10109.9994.2607.32210.0009.3129.999−3.6393.88610.0008.490
SVC2 Rating−10109.9992.9829.6089.7479.6223.844−1.4873.411−0.0657.591
TCSC1 Location--383838243616739118
TCSC2 Location--153640138192934331
TCSC1 Rating00.50.5000.4900.5000.0020.4990.5000.4680.4700.0020.500
TCSC2 Rating00.50.5000.3250.4750.0130.4990.0130.4900.3550.0000.147
TCPS1 Location--36433338144353130
TCPS2 Location--41173813942251133
TCPS1 Rating−55−4.9993.7884.884−5.0004.5403.1730.4864.5912.7611.494
TCPS2 Rating−55−4.9981.9934.772−1.1264.697−0.508−1.527−0.9015.0002.191
VSI 0.0960.1000.0960.1080.102-----
Total Gross Fuel Cost ($/h) -----1120.9961125.9701138.3571187.2871148.359
Table 12. Comparison of the simulation results for single objectives.
Table 12. Comparison of the simulation results for single objectives.
Objectives FunctionsGNDOMVOALOSCAIMOSHADE-SFDE-SFABC-SFPSO-SFFPA-SFMSA-SF
Total F.C ($/h)806.999808.030809.449818.654814.865807.0166807.4907807.399808.2543810.397808.0393
Emission (T/h)0.1380.1380.1380.1380.138------
Ploss (MW)1.6471.7351.6862.4821.8801.74671.74671.90681.89642.25972.1442
V.D (p.u)0.1240.1500.1770.2270.165------
Lmax0.0960.1000.0960.1080.102------
Total Gross F.C ($/h)1120.9961125.9701138.3571187.2871148.3591104.0771113.6761116.3651118.6011164.7191122.331
Table 13. Multi-objective simulation results for case 7 and case 8.
Table 13. Multi-objective simulation results for case 7 and case 8.
Control & State VariablesMinMaxCase-7Case-8
GNDOGWOMFOMVOGNDOGWOMFOMVO
PTG2208045.89654.59247.91547.56831.21244.40640.20037.027
PWG507554.82858.92253.22052.05564.77665.06960.99063.027
PTG8103523.32220.82125.40324.40117.91323.01525.01924.049
PWG1106046.87842.84841.48446.11648.38151.17151.42943.893
PTG13124019.17117.91422.28522.09414.45323.60419.01514.758
V10.951.11.0311.0431.0671.0361.0851.0991.0871.065
V20.951.11.0271.0191.0551.0251.0751.0901.0821.049
V50.951.11.0221.0001.0291.0101.0571.0771.0691.027
V80.951.11.0231.0011.0321.0201.0541.0711.0741.041
V110.951.11.0311.0351.0391.0221.0871.0771.0731.069
V130.951.11.0261.0201.0451.0291.0751.0791.0521.037
T110.91.10.9860.9561.0070.9870.9721.0031.0341.016
T120.91.11.0060.9661.0021.0011.0201.0371.0450.966
T150.91.11.0410.9531.0511.0381.0491.0711.0571.017
T360.91.11.0030.9561.0120.9810.9900.9961.0331.004
SVC1 Location----1511151412221619
SVC2 Location----2225191017122313
SVC1 Rating−10100.6960.846−1.6332.629−0.0188.596−3.9986.891
SVC2 Rating−10104.6115.523−0.5491.4223.179−3.6513.2312.030
TCSC1 Location----2413151812192020
TCSC2 Location----1617261427142517
TCSC1 Rating00.50.2940.4800.2210.2610.4370.4030.2790.204
TCSC2 Rating00.50.2130.2020.1600.2430.1440.3570.3220.184
TCPS1 Location----1916232328252120
TCPS2 Location----2226332314153119
TCPS1 Rating−55−1.239−2.128−1.398−0.7470.3040.797−0.1870.372
TCPS2 Rating−551.618−0.2980.3730.0180.988−4.1131.4730.005
TFC ($/h)----861.489865.902863.797865.736845.768883.242870.299852.611
Emission (Ton/h)----0.1660.1640.1650.163--------
APL (MW)------------4.0103.0243.4584.066
Table 14. Multi-objective simulation results for case 9 and case 10.
Table 14. Multi-objective simulation results for case 9 and case 10.
Control & State VariablesMinMaxCase-9Case-10
GNDOGWOMFOMVOGNDOGWOMFOMVO
PTG2208041.85440.05635.40641.92952.33950.96047.71157.834
PWG507546.21344.29045.65943.00073.72959.37867.05272.012
PTG8103510.53917.20612.90517.79427.54833.07925.00134.062
PWG1106040.50737.04742.06737.37455.54750.93047.00445.086
PTG13124014.44114.58317.10715.54325.04628.26033.66229.799
V10.951.11.0951.0851.0771.0771.0271.0451.0391.019
V20.951.11.0881.0681.0701.0621.0261.0361.0371.030
V50.951.11.0821.0541.0521.0351.0211.0241.0261.017
V80.951.11.0891.0811.0811.0781.0111.0161.0091.004
V110.951.11.0811.0671.0821.0811.0421.0721.0481.033
V130.951.11.0931.0821.0781.0801.0311.0101.0361.045
T110.91.11.0220.9790.9711.0010.9841.0101.0320.979
T120.91.10.9570.9450.9591.0120.9630.9560.9490.929
T150.91.11.0051.0170.9790.9810.9640.9640.9850.980
T360.91.10.9090.9060.9270.9100.9470.9330.9430.946
SVC1 Location--17211521823229
SVC2 Location--2414222719101919
SVC1 Rating−10107.2701.6291.3953.6040.7036.4960.374−0.378
SVC2 Rating−10108.3033.6206.5838.4186.197−2.4546.1934.988
TCSC1 Location--383821151471126
TCSC2 Location--1615381330242827
TCSC1 Rating00.50.4700.4460.2730.4730.2810.0970.3570.309
TCSC2 Rating00.50.4870.4100.4030.4380.1860.4330.4120.384
TCPS1 Location--3013111423152531
TCPS2 Location--317142329142815
TCPS1 Rating−551.5270.6930.0010.2860.8881.0630.6352.025
TCPS2 Rating−552.980−0.065−1.394−0.2720.2402.6460.126−0.864
TFC ($/h)--810.049814.201813.755817.163----
VSI--0.1050.1100.1160.115----
Total Gross F.C ($/h)------1173.3221208.4161212.9591199.459
Voltage Deviation (p.u)------0.1760.2300.2090.169
Table 15. Multi-objective simulation results for case 11 and case 12.
Table 15. Multi-objective simulation results for case 11 and case 12.
Control & State variablesMinMaxCase-11Case-12
GNDOGWOMFOMVOGNDOGWOMFOMVO
PTG2208045.29344.97550.29845.85631.58960.14733.79936.035
PWG507560.66668.79562.29457.15562.04261.71766.09452.766
PTG8103524.27124.48527.17926.03124.65127.55025.23225.207
PWG1106046.55147.89953.20647.44040.37557.22045.11241.884
PTG13124019.14617.28720.73414.41918.30715.56916.09019.407
V10.951.11.0331.0931.0601.0151.0661.0461.0771.045
V20.951.11.0191.0861.0601.0081.0581.0401.0631.032
V50.951.10.9971.0691.0551.0011.0400.9961.0411.017
V80.951.11.0151.0811.0511.0061.0461.0391.0391.009
V110.951.11.0201.0401.0551.0511.0511.0751.0341.037
V130.951.11.0351.0511.0211.0191.0281.0291.0521.050
T110.91.11.0000.9720.9990.9660.9870.9630.9630.964
T120.91.11.0031.0701.0551.0071.0200.9761.0020.973
T150.91.11.0611.0991.0751.0431.0351.0871.0511.015
T360.91.10.9521.0761.0240.9850.9640.9740.9640.965
SVC1 Location--1121131814151818
SVC2 Location--131215820231918
SVC1 Rating−10100.4956.230−2.056−0.1090.0573.863−0.308−2.293
SVC2 Rating−10103.8614.6655.1598.1320.5993.889−0.8187.033
TCSC1 Location--201621426392524
TCSC2 Location--1924242217221823
TCSC1 Rating00.50.1730.2690.3060.3310.2170.1010.1330.278
TCSC2 Rating00.50.2050.2880.3360.3650.3120.1940.2590.237
TCPS1 Location--152430820191725
TCPS2 Location--143420361620225
TCPS1 Rating−550.083−2.2131.397−0.2220.555−0.8910.3840.657
TCPS2 Rating−552.118−0.2410.8630.454−0.207−1.208−0.6520.936
TFC ($/h)--868.806877.252887.402860.794849.334893.297860.198845.253
Emission (Ton/h)--0.1620.1570.1500.167----
APL (MW)--4.0343.3703.2544.3474.3283.4533.9224.825
VSI------0.4160.3560.3780.269
Table 16. Multi-objective simulation results for case 13.
Table 16. Multi-objective simulation results for case 13.
Control & State VariablesMinMaxCase-13
GNDOGWOMFOMVO
PTG2208049.76352.44359.63445.029
PWG507547.29162.62261.40453.782
PTG8103526.71713.34015.91117.846
PWG1106029.20658.71947.38047.174
PTG13124026.88819.44318.67218.177
V10.951.11.0481.0341.0411.053
V20.951.11.0381.0121.0271.051
V50.951.11.0280.9800.9901.009
V80.951.11.0400.9981.0261.027
V110.951.11.0531.0691.0471.024
V130.951.11.0291.0541.0511.021
T110.91.11.0120.9310.9740.978
T120.91.10.9681.0460.9941.005
T150.91.11.0230.9751.0321.055
T360.91.10.9690.9431.0060.967
SVC1 Location--10152420
SVC2 Location--18251520
SVC1 Rating−1010−4.1766.2331.5601.932
SVC2 Rating−10103.7456.2594.7660.737
TCSC1 Location--2263119
TCSC2 Location--19271722
TCSC1 Rating00.50.4230.2080.1680.228
TCSC2 Rating00.50.2780.3000.1250.114
TCPS1 Location--20282720
TCPS2 Location--27263220
TCPS1 Rating−55−0.164−3.317−0.896−0.585
TCPS2 Rating−551.413−2.786−0.9760.658
TFC ($/h) 863.417883.480874.361850.096
Emission (Ton/h) 0.1730.1570.1600.175
APL (MW) 5.1894.6284.2654.809
Voltage Deviation (p.u) 0.3590.2250.4500.529
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Pandya, S.B.; Visumathi, J.; Mahdal, M.; Mahanta, T.K.; Jangir, P. A Novel MOGNDO Algorithm for Security-Constrained Optimal Power Flow Problems. Electronics 2022, 11, 3825. https://doi.org/10.3390/electronics11223825

AMA Style

Pandya SB, Visumathi J, Mahdal M, Mahanta TK, Jangir P. A Novel MOGNDO Algorithm for Security-Constrained Optimal Power Flow Problems. Electronics. 2022; 11(22):3825. https://doi.org/10.3390/electronics11223825

Chicago/Turabian Style

Pandya, Sundaram B., James Visumathi, Miroslav Mahdal, Tapan K. Mahanta, and Pradeep Jangir. 2022. "A Novel MOGNDO Algorithm for Security-Constrained Optimal Power Flow Problems" Electronics 11, no. 22: 3825. https://doi.org/10.3390/electronics11223825

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