Optimal Scheduling of Non-Convex Cogeneration Units Using Exponentially Varying Whale Optimization Algorithm

: This paper proposes an Exponentially Varying Whale Optimization Algorithm (EVWOA) to solve the single-objective non-convex Cogeneration Units problem. This problem seeks to evaluate the optimal output of the generator unit to minimize a CHP system’s fuel costs. The nonlinear and non-convex characteristics of the objective function demands a powerful optimization technique. The traditional Whale Optimization Algorithm (WOA) is improved by incorporating four different acceleration functions to ﬁne-tune its performance during exploration and exploitation phases. Among the four variants of the proposed WOA, the emphasis is laid on the EVWOA which uses the exponentially varying acceleration function (EVAF). The proposed EVWOA is tested on six different small-scale to large-scale systems. The results obtained for these six test systems, followed by a statistical study highlight the supremacy of EVWOA for ﬁnding the best optimal solution and the convergence traits.


Introduction
Cogeneration units are a single entity with the ability to produce thermal and electrical energy. The ever-growing need for electrical and heat energy has contributed to an increase in demand for these systems. There are manifold advantages of these schemes. In comparison to coal-based thermal power plants and boilers, some of them have higher productivity, economic benefits, lower environmental emissions, etc. These CHP systems have a performance of about 80 to 85% [1], in comparison to the lower values of around 30 to 40% [2] of coal-based power plants. In CHP systems, by using waste heat that is supplied to different regions as per their needs, they boost the aggregate efficacy of the system. Figure 1 shows a sample sketch of the Heat Recovery Steam Generator, which is one of the configurations of operating a CHP unit. Holding all these advantages of cogeneration systems in mind, having the system run at the highest possible efficiency is of utmost importance. This is made possible by knowing the system's best dispatch schedule. It leads to a CHPED conundrum in which, considering all the operational constraints, electrical and heat energy outputs of each generator unit are computed, such that the fuel cost is the lowest. This CHPED conundrum becomes complicated and cumbersome due to the nonconvex and nonlinear structure of cost function and considering various operational difficulties. Therefore, to compute the best optimal solution, a strong optimization technique is required. The technique should be able to prevent local trapping in the optimal solution to achieve this best optimal solution and therefore pursue a global best solution by having a better equilibrium between the exploration and exploitation phases. To solve the above-mentioned conundrum, several various methodologies have been used in numerous studies. Figure 2 shows a summary of the various optimization techniques. Some of the conventional mathematical approaches used in numerous studies include Lagrangian Relaxation [3], Dual and Quadratic Programming [4], and Bender's decomposition [5], etc. However, these methods are unable to find the required solution adequately because they are unable to work successfully with the non-convexities and hard constraints. To deal with the problem of not being able to handle the complexities, researchers have now started using meta-heuristic methods and have found a rational solution. In earlier literature, some of the methods used were: Ant Colony Optimization (ACO) [6], Gravitational Search Algorithm (GSA) [7], Krill Herd Algorithm (KHA) [8], Particle Swarm Optimization (PSO) [9] and its variants such as Selective PSO [10], Time-Varying Acceleration Coefficients with WOA. Reference [25] suggested WOA focused on the Levy Flight Trajectory to solve the issue of global optimization. To ensure a better negotiation between the exploration and exploitation of the WOA, the suggested approach is very useful. Diego et al., in [32], introduced enhanced Chaotic WOA (CWOA) for solving solar photovoltaic cell parameter estimation. Chaotic maps are used for calculation in this CWOA method and the internal parameters of the method are automatically chosen. Aziz et. Mohamed Abd El Aziz et al. to achieve the optimum multi-level thresholds for image segmentation in [33], suggested a hybrid approach consisting of WOA and Moth-Flame Optimization (MFO). A new MLP training method based on the recently developed WOA for the optimization of relation weights in neural networks is proposed in Reference [34]. Mohamed Abdel-Basset et al. presented a framework for hybrid whale optimization with a local search approach to address the issue of permutation flow shop scheduling in [35]. Kaur et al., in [26], by considering multiple chaotic maps, the chaotic WOA (CWOA) approach was suggested. For fine-tuning the key parameters of the WOA system, these chaotic maps are very useful, for maintaining the exploration and exploitation of the proposed CWOA method. A hybrid approach consisting of Support Vector Machines and WOA for solving the challenge of detecting spammers in the issue of online social networks is proposed in Reference [36]. An improved WOA (IWOA) for the correct extraction of the parameters of various solar photovoltaic model problems is proposed in Reference [27]. Two prey-seeking techniques have been implemented by this IWOA approach to effectively preserve the proper balancing between exploration and exploitation to boost WOA efficiency. In addition, Jadhav et al., recommended a hybrid solution consisting of a grey wolf optimizer with whale optimization to solve the optimization problem of data clustering [37]. To solve the 0-1 knapsack problem with various scales, Mohamed Abdel-Basset et al. suggested an improved WOA (IWOA) [38].
In this paper, an Exponentially Varying WOA (EVWOA) is proposed to improve and enhance the performance of the basic WOA in terms of convergence, maintaining the balancing between global exploration and local exploitation throughout the search process and improving its efficiency. The performance of basic WOA is improved by introducing an exponentially varying acceleration function. This helps the proposed EVWOA method for providing proper balancing of both exploration and exploitation phases throughout the search process. The solution in the proposed EVWOA converges faster, is robust and efficient, and obtains a better optimal solution without premature convergence. The proposed EVWOA technique is implemented in this paper for giving a better optimal solution than the traditional WOA along with different variants of WOA. The proposed WOA and its different variants have been tested on six different test systems ranging from a small-scale system of four units to a large-scale system comprising of 96 units with several related constraints. The results obtained by the proposed EVWOA method after 100 different unbiased trials are compared with the basic WOA, its variants, and different published recent algorithms. It can be observed that the proposed WOA method provides good quality solutions in terms of cost, robustness, feasible and obtain effective convergence characteristic.
The main highlights of the paper are as follows: • Four different variants of WOA are considered for improving and enhancing the performance of the basic WOA. In the first variant of WOA, the acceleration function is randomly selected, and it is written as Randomly Varying WOA (RVWOA The paper is organized as: the problem formulation is explained in Section 2; a description of the basic WOA and its shortcomings is given in Section 3. This section also explains, in brief, the proposed Exponentially Varying WOA (EVWOA); Section 4 presents the simulation results and discussion; the conclusion is presented in Section 5.

Problem Formulation
The objective of the CHPED conundrum is to obtain the best output values for the interconnected electrical energy, heat energy, and CHP plants in a manner that the total fuel cost is minimum, and all the related constraints are satisfied, and energy demands are met.
This can be stated as [3]: The equation of the electrical energy generators taking into account valve-point loading is given as [7]: where e l , f l and g l are the cost coefficients for the lth electrical energy generator unit and h l and i l are the coefficients for modeling the valve-point. The non-convex nature of the CHPED problem can be attributed to this sinusoidal term. Figure 3 shows graphically the effect of valve-point loading. The fuel cost for the CHP plants is given by [6]: where e m , f m , g m , h m , i m and j m are the coefficients of the mth CHP unit. The fuel cost for the heat energy generators is given by [7]: where e n , f n and g n are the cost coefficients of the nth heat energy generator.

Constraints
The equality and inequality constraints, which form the energy balance equations are stated as [18]:

Equality Constraints
The electrical energy balance equation is stated as: The above Equation (5) signifies that the total electrical energy output from the electrical energy generators and the CHP units should be sufficient to satisfy the electrical energy demand, as well as certain energy losses that may arise in the system in due course of transmission from the plant to the consumer. This electrical energy loss is modeled using Kron's loss formula and is stated as below: In addition, the heat energy balance equation is given as: Equation (7) signifies that the total heat energy output from the CHP units and heat generator units should be sufficient to meet the thermal energy requirement. Equations (8) to (11) form the heat-power Feasible Operating Region (FOR) of the CHP units. The FOR, as the name suggests is the zone in which CHP units should have their heat and electrical energy outputs to operate. A sample sketch of this FOR is shown in Figure 4.

Proposed Exponentially Varying Whale Optimization Algorithm (EVWOA)
This section discusses the equations which form the WOA and some limitations which hampers its performance. Later, an acceleration function term [ζ(itr)] is proposed, which is then appended in the characteristic equations of the WOA. The WOA is an experience-based algorithm which is having a hunting strategy adopted by humpback whales. The peculiar bubble-net hunting pattern of these whales is the inspiration for the exploitation strategy. The characteristic equations of the WOA are given as: [23][24][25][26][27].

Encircling Prey
The location of the prey is traced by the leader whale and is then encircled by the other whales. Initially, a leader whale is selected among the group and the others follow it. This phenomenon is given by [23]: Here, α and β are the coefficient vectors, itr is the present iteration, while Z and Zbest * are the position and best position vectors respectively, which are updated in each cycle. The coefficient vectors are given as: where τ varies from 2 to 0 in a decreasing manner, while rand is a random value in [0,1]. In basic WOA, as per Equation (13), the updated position is regulated by multiplication of coefficient vector ( α) and distance of the ith whale to the prey ( Dist) which is subtracted from current best position, which is further dependent upon coefficient vector ( α), where α is calculated from Equation (14). From Equation (12), Dist is dependent upon coefficient vector ( β) and from Equation (15), it can be said that β is dependent upon rand which is a random vector. It can be observed that α is dependent upon τ which is decreasing linearly with respect from 2 to 0 as per the iterations. This may cause uncontrolled updated positions during the iterations. It may also create an unbalance between local exploration and global exploitation in the search space and further it may result in poor convergence. For improving this, ref. [39,40] introduced an inertia weight into WOA similar to the PSO to fine-tune the influence of the current best solution. According to references [39,40], in the proposed WOA, an acceleration function (ζ) is introduced.

Exploitation Phase
As stated earlier, the exploitation phase is based on the bubble-net hunting pattern of the whales. Here, two different techniques are used to mimic this behavior.

Shrinking Encircling
Shrinking encircling is obtained when τ reduces to 0 from 2 in Equation (14). Hence as τ decreases, a lesser variation is seen in the value of α, which lies in [− τ, τ].

Spiral Updating Position
In this method, it is assumed that the whale and the prey are present at (Z,W) and (Z best , W best ). The following equation defines the helical movement of the whale towards the prey [24]: Z(itr + 1) = Dist times γη × cos (2πη)+ Zbest * (itr) (18) Dist is the distance between the ith whale and the prey. It is written as: γ lies between [−1,1] and is fixed valued and helps to model the spiral shape of the path.
By assuming that around the prey, whales either swim in a spiral path or a shrinking encircling path, an equal probability is assumed for selecting either of the paths. This is stated as [25]: where prob is the chance of selecting either path.
In the proposed EVWOA, the above equations are modified as:

Exploration Phase
This phase comprises searching for the prey, in which the value of α can be changed for locating the prey. Whenever |α| > 1, the current search agent is to be moved away from the leader agent. The location of an agent in this phase is improved by erratically selecting a particle. This phase is described by the following [26]: In the proposed EVWOA, these equations are varied as [28,29]:

Choice of Acceleration Function
Four different variants of WOA are considered by incorporating different acceleration functions for improving and enhancing the WOA performance. These different variants are as follows: is linearly varied and decreasing from 0.9 to 0.1 as in References [29,39,40]. The expression is written below: 3. Sinusoidally Varying WOA (SVWOA): In this variant, the acceleration function (ζ) is varied sinusoidally from 0.9 to 0.1 as in [39]. The mathematical expression is given below: where θ = X * itr + Y and the coefficients X and Y are calculated by Equations (30) and (31) and iteration (itr) is varied from itrmin to itrmax.
4. Exponentially Varying WOA (EVWOA): In this proposed variant, the acceleration function (ζ) is varied exponentially as in [40]. The exponential variation of ζ is calculated by the following relation: where l = itr itr max , itr min ≤ itr ≤ itr max and k w is the ratio of maximum minimum bounds of the acceleration function.
The linear, sinusoidal, and exponential variations of acceleration function with iterations are shown in Figure 5. It can be analyzed from the figure that ζ-linear is decreasing linearly. By this, there may be a chance of local trapping or premature convergence during the iterations. Moreover, in the case of ζ-sinusoidal, during the first half (50%) of iterations, ζ is maintained higher for wider exploration, and in the second half, ζ is lower for better exploitation during the computation process. However, when ζ-exponential is considered, initially exploration is wider and exploitation capability is improved by a sharp fall in ζ during the iterations. This will improve and efficiently maintain the exploration and exploitation phase of the proposed method. Hence, it is favorable that ζ exponential is selected for the best performance for the WOA. Therefore, further in this paper, the emphasis is laid on the EVWOA.     For the above Figures 7 and 8, the population of whales considered was 100 whales and number of generators considered were 60.

Testing of the WOA Variants on Well-Known Benchmark Functions
The proposed EVWOA is tested on ten different well-known benchmark functions. The benchmark functions are defined in Table 1 and their optimum and bound limits are known to be a priori. The functions used are Ackley, Generalized Griewank, Generalized Rastrigin, Generalized Rosenbrock, Schwefel 1.2, Generalized Schwefel 2.26, Sphere, Sixhump Camel-back, and Goldstein-Price. The data for the benchmark functions are referred from [41].
The parameters for EVWOA are: Population of whales = 100 and the Maximum number of iterations (itrmax) = 1000. The results obtained using proposed EVWOA after 100 different trials are compared with some of the well-known optimization algorithms in Table 2. MATLAB ver. 2020a is used to simulate these benchmark functions. The computer specifications are Processor-Intel i5 7th Generation @ 2.50 GHz, RAM -8 GB, and Storage Capacity-1TB.  From Table 2 it can be observed that for the non-shifted functions F1, F2, F3, F5, and F7, the proposed EVWOA can obtain the local minima of an absolute zero, compared to some of the other algorithms considered for comparison. For the shifted functions F4, F6, F8, and F9, the performance of the proposed EVWOA is found to be satisfactory compared to the other algorithms, where except in F9, EVWOA obtains the minima of 0, especially when compared with FWA variants [41], PSO [42] and HS [42].

Results and Discussions
In this paper, six different case studies have been considered and solved using WOA and all the different variants of the proposed WOA, i.e., RVWOA, LVWOA, SVWOA, and EVWOA. A summary of the six different case studies is presented in Table 3.  1  1  2  1  4  200  115  2  1  3  1  5  300  150  3  4  2  1  7  600  150  4  13  6  5  24  2350  1250  5  26  12  10  48  4700  2500  6  52  24  20  96  9400  5000 For solving the proposed problem, WOA and different proposed variants of WOA are implemented in MATLAB ver. 2020a on PC with 8 GB RAM, 1 TB storing capacity, and an Intel i5 processor. Each case study is simulated 100 different times, and the schedule obtained for the best cost using EVWOA is presented in this paper. For showing the algorithm's robustness, the best, mean and worst cost, along with the standard deviation (STD) of cost and computation time is also presented. These factors would establish the algorithm's superiority compared to the basic WOA and other latest techniques in the literature. For each of the case studies, the population of whales and iterations, both are considered to be 100.

Case Study 1
This case study data has been referred from [3] and is given in Appendix A. The electrical and thermal energy requirement of the whole system is 200MW and 115MWth. This system is considered to be a small-scale system. The optimal schedule obtained by the proposed EVWOA for the best cost is presented in Table 4. The results obtained by the proposed EVWOA after 100 different trials are presented in Table 5. It can be observed from Table 5 that the results obtained by WOA, and other variants of WOA are giving much better results as compared with recently published methods such as Firefly Algorithm (FA) [43], Genetic Algorithm (GA) [44], Differential Evolution (DE) [15] and Grey Wolf Algorithm (GWO) [22]. It can also be analyzed from Table 5 that among WOA and different variants of WOA, the basic WOA is improved by adding different acceleration functions, the quality of solution in terms of best cost (minimum cost), mean cost (average cost) and Standard Deviation (STD) is also improved. From Table 5, proposed EVWOA gives minimum best cost, minimum average cost, and minimum STD than WOA and other different variants of WOA. This establishes the robustness of the proposed EVWOA. The average CPU time taken by the proposed EVWOA is less than DE [15] and GWO [22]. It shows that the proposed EVWOA converges faster with less CPU time. Figure 9 compares the best cost of all the above-mentioned algorithms in a graphical format. Among recently published methods i.e., DE [15], GWO [22], FA [43] and GA [44], DE [15] has the lowest cost of $9,236.14/h. By comparing the minimum cost obtained by DE [15] and the proposed EVWOA, EVWOA gives an hourly saving of $146.72, which means that $1,285,267 could be saved annually. This leads to the fact that the proposed EVWOA is capable of obtaining better quality solutions compared to other methods in terms of the best and average fuel cost and is also an efficient optimization algorithm to solve dispatch problems.

Case Study 2
This case study data has been taken from [45] and is given in Appendix B. The optimal generator output schedule obtained by the proposed EVWOA is given in Table 6. The best optimal solution generated by WOA and different variants of WOA after independence 100 trials are compared with the cost obtained with Firefly Algorithm (FA) [43], Harmony Search (HS) [45], Invasive Weed Optimization (IWO) [20], Cuckoo Optimization Algorithm (COA) [46] and Exchange Market Algorithm (EMA) [21] in Table 7. It is observed from this Table 7 that WOA and other different proposed variants of WOA provide best minimum best cost as compared with latest published methods i.e., FA [43], HS [45], IWO [20], COA [46] and EMA [21]. It can also be observed from this Table 7 that the proposed EVWOA again obtained the best quality solution with almost 0 STD than WOA and RVWOA, LVWOA, and SVWOA. This shows that EVWOA is efficient, robust, and feasible. Figure 10 compares the best cost obtained by all the above-mentioned algorithms. From Table 7, COA [46] has the lowest hourly cost of $13,672.83/h. Therefore, by implementing EVWOA for this test system, hourly savings of $0.04 are possible, which amounts to a yearly saving of $350.40. It still proves that EVWOA can work more efficiently than other latest techniques.

Case Study 3
The case study parameters have been taken from [47] and are also provided in Appendix C. Power system transmission losses are taken into account for this system. The optimal generator schedule by EVWOA is provided in Table 8. The quality solution obtained by WOA and its proposed variants after 100 trials are compared with the best cost obtained using Teaching-Learning-Based Optimization (TLBO) [13], Gravitational Search Algorithm (GSA) [7] and Differential Evolution (DE) [48] in Table 9. This shows that the minimum cost of WOA and its proposed variants are much better than TLBO [13], GSA [7] and DE [48] with less CPU time. From this comparison, it can be concluded that EVWOA provides minimum best cost and minimum average cost than the WOA and other variants of WOA i.e., RVWOA, LVWOA, and SVWOA. The standard deviation of EVWOA is also quite good and is about 0.005, meaning that in almost each trial best cost is achieved. Figure 11 compares the best cost achieved by all the above algorithms in a graphical format. Among the other published algorithms TLBO [13], GSA [7] and DE [48], GSA [7] provides the lowest cost of $9,912.69/h. However, by using the proposed EVWOA an hourly saving of $173.18 is achieved, which leads to a yearly saving of $1,517,057, which is a huge amount. This high figure of saving re-establishes the superiority of the proposed EVWOA as it can provide a better solution with less CPU time for a system where system transmission losses and other constraints are considered.

Case Study 4
The data for this test system is referred from [24] and is given in Appendix D. The schedule for the best cost obtained by EVWOA is provided in Table 10. The solutions obtained by WOA and its variants after 100 trials are presented in Table 11. It can be observed from this table that the quality solution of WOA and different variants of WOA are much better than Oppositional Teaching-Learning-Based Optimization (OTLBO) [13], Exchange Market Algorithm (EMA) [21], Improved Group Search Optimization (IGSO) [49] and Civilized Swarm Optimization (CSO) [48]. It can also be analyzed from this table that again proposed EVWOA gives minimum cost and minimum average cost with a lesser standard deviation of around 0.003 than WOA, RVWOA, LVWOA, and SVWOA. Figure 12 compares the best cost of all the above algorithms graphically.
EMA [21] among other published algorithms OTLBO [13], IGSO [49] and CSO [48] has the lowest best and average costs. When comparing the average cost obtained by EMA [21] and proposed EVWOA, it can be said that the proposed EVWOA provides an hourly saving of $4459.5447 and this leads to a yearly saving of $39,065,637.852 which is indeed a huge sum. This high figure is due to better handling of the related constraints while using EVWOA to solve the dispatch problem.

Case Study 5
In this case study, 48 units large test system is considered. The data is referred from [9] and is provided in Appendix E. The optimal generator outputs obtained by EVWOA are as shown in Table 12. The cost obtained by WOA and its variants are compared with the costs obtained using Cuckoo Optimization Algorithm (COA) [46], Crisscross Optimization Algorithm (CSO) [50], Modified Particle Swarm Optimization (MPSO) [51] and Group Search Optimization (GSO) [52] in Table 13. It can be said that the quality solution, i.e., minimum cost, average cost, and STD obtained by WOA, RVWOA, LVWOA, SVWOA, and EVWOA are much better than the solution presented by COA [46], CSO [50], MPSO [51] and GSO [52]. It shows that WOA and its variants perform superior to other established methods with less CPU time. It can also be observed from Table 13 that minimum cost, average cost, and STD obtained by the proposed EVWOA outperforms WOA, RVWOA, LVWOA, and SVWOA. It implies that the proposed EVWOA dominates all other techniques. Figure 13 compares the best cost of all the above algorithms graphically. GSO [52] has better best minimum cost and average costs among COA [46], CSO [50] and MPSO [51]. When comparing the average cost of $116,463.65/h obtained by GSO [52] and EVWOA's average cost of $96,657.72/h, it can be observed that EVWOA gives an hourly saving of $19,805.93 and the yearly savings are $173,499,947 which is a very huge amount. This again is attributed to better constraint handling employed while solving with EVWOA.

Case Study 6
This is also considered to be a large system consisting of 96 units. This case study is included in this work to highlight the superiority of EVWOA in efficiently handling large test systems. The data is taken from [24] and is included in Appendix F. The optimal generator outputs using EVWOA are presented in Table 14. The quality solution obtained by WOA and its different variants after 100 trials are presented in Table 15 and is compared with Real Coded Genetic Algorithm with improved Muhlenbein mutation (RCGA-IMM) [53]and TVAC-PSO [9]. It can be observed that WOA and its variants perform superior to RCGA-IMM [53] and TVAC-PSO [9] with less CPU time. Again, it can be said by analyzing the results in Table 15 that proposed EVWOA performs better than WOA, RVWOA, LVWOA, and SVWOA. This shows that by introducing exponentially varying acceleration function in basic WOA, the proposed method performs well and during the computation process, the exploration and exploitation are properly balanced for getting the best optimal quality solution. Figure 14 compares the best cost of all the above algorithms graphically. The best minimum cost of EVWOA $164,691.71/h is compared to the best minimum cost of TVAC-PSO [9], i.e., $239,139.5018/h. The hourly savings are calculated and found by EVWOA to be $74,447.7918, which when seen yearly, amounts to $652,162,656. These high figures prove the worthiness of EVWOA for obtaining a good solution, with good standard deviation and a good computation time compared to other different published algorithms. The main objective of including this case study in this paper is to highlight the supremacy of the EVWOA algorithm compared to others. Lesser minimum cost and average cost obtained by EVWOA prove the point. The standard deviation is also very low, i.e., 0.009.

Discussion About Convergence Characteristics
The convergence characteristics obtained by WOA, RVWOA, LVWOA, SVWOA, and EVWOA for best cost vs no. of iterations for this large test system during sample trial are shown in Figure 15. It can be observed from this figure that initially basic WOA explores the search area very well but in the latter part of the iterations, it traps in local optima. However, by incorporating randomly varying acceleration function in basic WOA, during initial iterations, exploration of RVWOA is not proper but in the latter half exploitation is improved but it is not able to get an optimal solution. Next, LVWOA is not performing better because there is an unbalance between local exploration and global exploitation. It can be observed from Figure 15 that during some initial iterations it explores well, but after that, it remains constant till the end of iterations. In SVWOA, a sinusoidally varying acceleration function is used. The basic tendency of the sinusoidal function is that for the initial half, it provides sufficient search space for exploration, and later half this search space is reduced for exploitation. It can be seen from this figure that the exploration of SVWOA is good but still, the exploitation potential remains weak. For improving the efficiency and maintaining the balance between exploration and exploitation, an exponentially varying acceleration function is used in the proposed EVWOA. It can be observed from Figure 15. that by adding exponentially varying acceleration function, EVWOA well explores the search area in a better way during the initial half and in the remaining half, the proposed method getting the best optimal solution because of better exploitation by avoiding the local optima. This proves that the proposed EVWOA obtains better quality solutions by properly managing the balance between exploration and exploitation throughout the search process.

Conclusions
The optimal scheduling of non-convex CHP units is a highly complex, non-convex, and hard constrained nonlinear optimization problem. In this paper, an Exponentially Varying WOA (EVWOA) method is proposed and implemented to enhance and improve the performance of basic WOA in terms of convergence, maintaining the balance between global exploration and exploitation throughout the search process, and improving its efficiency. This variation is done by the introduction of an exponentially varying acceleration function in the characteristic equations of the original WOA.
This proposition, as well as the other variations of WOA, namely RVWOA, LVWOA, and SVWOA, are tested on several well-known standard benchmark functions to test their ability and effectiveness in finding the solution. The obtained results are compared to some of the existing optimization algorithms in the literature. It is found that all the variations perform satisfactorily and can get the desired results.
Next, the proposed EVWOA is tested on six small to large different test systems of CHPED problem consisting of several related operational constraints. The simulation results provided by the proposed EVWOA after 100 independent trials are compared with basic WOA, other remaining variants of WOA, and recently published results obtained by different methods. The comparison results show that the proposed EVWOA performs much better than other latest published methods and can obtain a much better optimal quality solution in terms of best minimum cost, mean of cost, and STD of cost in less CPU time. It can also be said that the proposed EVWOA is providing promising results with an efficient constraint handling method for solving the CHPED problem.
Apart from the conventional single-objective CHPED problem, the proposed EVWOA can be used for solving various other power system optimization problems. Some of them to be named are the multi-objective CHP economic emission dispatch (CHPEED), combining economic load dispatch consisting of conventional electrical generator units with renewable energy sources, multi-area economic dispatch, dynamic load dispatch, etc.

Appendix A
Case Study 1 details are referred from Reference [3]. Electrical energy generator units data is as per Table A1.  Cogeneration unit data is as per Table A2. Heat energy generator unit data is as per Table A3.

Appendix B
Case Study 2 details are taken from [47].
Electrical energy generator data is as per Table A4.  Cogeneration unit data is as per Table A5. Heat energy generator data is as per Table A6.  Appendix C Details of Case Study 3 are referred from ref. [45]. Electrical energy generator data is as per Table A7.  Cogeneration unit data is as per Table A8. Heat energy generator data is as per Table A9.

Appendix D
Data of Case Study 4 are referred from ref. [24]. Electrical energy generator data is as per Table A10.  Cogeneration unit data is as per Table A11.   Cogeneration unit data is as per Table A14. Appendix F Case Study 6 data are referred from ref. [24]. Electrical energy generator data is as per Table A16.  Cogeneration unit data is as per Table A17.