1 Introduction

In 2020s, energy consumption has been growing at an exponential rate thanks to the higher energy efficiency and strong demand. As early as 1973, Lincoln et al. [28] pointed out that energy conservation should be paid attention to and the sustainable development of environment and resources should be guaranteed through technological progress and social science [55]. With the significant improvement of computational capacity, an increasing number of industrial equipment began to use computer-aided design and optimization to improve the energy efficiency and reduce economic expenditure [49, 51]. Plate-fin heat exchanger (PFHE) is a common type of heat exchanger, also known as brazed aluminum heat exchanger in the thermal engineering. Compared with other types of heat exchangers, PFHE has compact structure, relatively small size, relatively high thermal load and lightweight, due to which it has been widely used in aerospace, railway, energy and chemical industries. In real-world cases, the majority of the design for heat exchanger is based on the experience of designers, without any unified design method and scheme [45]. Bio-optimization, which is of high compatible capacity in solving various optimization problems, is promising to effectively solve PFHE optimal design problem. To solve the design problem of plate-fin heat exchanger, the geometric limitation of heat exchanger, such as shape and size, should be determined based on the actual engineering needs. Then the objective function is derived based on the design objective. The fluid data and constraints are brought into the operation, and the appropriate algorithm is selected for them. Finally, the optimal design data are obtained.

Bio-inspired optimization is based on a variety of physical and biology phenomena in the nature. Several typical examples include genetic algorithm (GA) based on evolutionary biology [46], particle swarm optimization (PSO) from simplified social models [25], simulated annealing algorithm from metallurgy for specific rate cooling operations [42], ant colony optimization (ACO) derived from the behavior of ants in food searching and paths discovering [13]. Moreover, a teaching learning-based optimization (TLBO) algorithm has also been proposed to simulate the teaching process of teachers and students [39].

In the heat exchanger optimal design optimization, many scholars have utilized bio-inspired optimization methods for solving the problems. The optimization of heat exchanger can be divided into two types: single objective and multi-objective. Featured optimization targets of optimal design include the amount of entropy produced, total heat exchange, maintenance costs and production costs. Port diameter, port horizontal distance, port vertical distance, plate thickness, fin length and number of plates are often used as geometric design variables [36]. Some recent proposed methods, such as Gaussian process regression and adaptive resampling, are adopted to minimize the cost [7]. Some scholars have combined traditional and novel methods, such as computational fluid dynamics (CFD) combination with GA, to optimize heat exchangers with Colburn factor j and the friction factor f as targets [2, 29]. Moreover, the optimal design problem is mostly solved by simple and improved optimization algorithm. Jonh and Krishnakumar [22] used GA to optimize the perforated plate matrix heat exchanger surfaces, where the optimal objectives were Colburn factor j and the friction factor f. Rao and Saroj [38] used Jaya algorithm to optimize the shell and tube heat exchanger (STHE) economically, taking into account the consistency and maintainability caused by scaling. In the same year, Rao and Saroj [37] improved the Jaya algorithm and proposed the elitist-Jaya algorithm to optimize the setup and operation costs of the STHE. Dhavle et al. [12] used cohort intelligence (CI) to optimize STHE design, and the total cost of the heat exchanger was taken as the target, and the exterior scene of the tube and the spacing of baffles were taken as the design variables. Vasconcelos et al. [43] proposed a falcon optimization algorithm (FOA) based on the predation behavior of falcons. The total cost of STHE was used as the objective function. Iver et al. [21] combined GA and CI and proposed an adaptive range genetic algorithm (ARGA) to solve the costs of the economic optimization problems and extended them to solve other complicated core mechanical engineering application problems.

Sine cosine algorithm (SCA) is a new meta-heuristic algorithm first proposed by Mirjalili in 2016 [31]. This optimization method is based on sine and cosine functions in mathematics and has been immediately widely applied in various optimization problems. Ekiz et al. [14] expounded the feasibility of SCA in solving constraint optimization problems and concluded that the successful feature of SCA lies in the smooth balance of exploration and development. Banerjee and Nabi [4] used SCA to optimize the shuttle’s re-entry trajectory back to Earth to find the best trajectory during the return. Attia et al. [3] used the modified sine cosine algorithm (MSCA) to solve the problem of optimal power flow (OPF), which accelerated the speed of SCA search and avoided the local optimal selection of the original algorithm as much as possible. Hekimoglu [19] used SCA to solve the parameters of automatic voltage regulator and proportional integral derivative (PID). In the optimization problem of automatic voltage regulator system, SCA has strong robustness. Mahadad and Srairi [30] use SCA to optimize load margin stability to improve power system security. Wang et al. [44] employed the multi-objective sine cosine algorithm (MOSCA) combined with wavelet neutral network (WNN) to accurately predict the wind behavior of wind power plants. Das et al. used SCA to optimize the hydrothermal scheduling problem in the economics of power system to minimize the total cost of power generation and solve the short-term hydrothermal scheduling problem [11]. Chandrasekaran et al. [8] solved the partial shading detection (PSD) and MPPT (maximum power point tracking) problems using SCA.

Since the proposal of canonical SCA, many scholars have proposed various algorithm improvement schemes or combinations with other algorithms to solve the complex problems. Elaziz et al. [1] combined the opposition-based learning (OBL) idea with SCA to improve the search accuracy of SCA. Nenavath and Jatoth [32] combined SCA with differential evolution (DE) and demonstrated that the hybrid SCA-DE algorithm solves the problems effectively. Sindhu et al. [40] combined elite strategy with SCA to make simple improvements to the algorithm and improve the search efficiency of SCA. Bureerat and Pholdee [6] proposed adaptive sine cosine algorithm integrated with differential evolution (ASCA-DE) to solve the problem of structural damage detection. Gupta and Deep [17] proposed an improved SCA with crossover scheme and proved the validity using CEC2014 benchmark test. In addition, they also solved the multilevel threshold processing of image segmentation using the new method and obtained competitive result. Chegini et al. [9] combined the PSO and Levy flight methods and proposed the PSOSCALF with stronger jumping ability and higher searching ability. Gupta and Deep [16] added self-adaptive to SCA to enhance the exploitation ability. Pasandideh and Khalilpourazari [33] combined crow and SCA and proposed the sine cosine crow search algorithm to balance the exploration and exploitation. Li et al. [26] proposed the enhanced brain storm SCA (EBS-SCA) to solve the problem of premature convergence of SCA on complex optimization problem, which significantly enhanced the algorithm performance. Gupta and Deep [18] proposed SCA with global optimization based on their previous studies, referred to simulated quenching algorithm search mechanism for improved algorithm, and used it to train multilayer perceptrons. Though a number of SCA variants have proposed, the optimization capacity of the algorithm is to be further improved, in particular toward the convergence speed aspect.

This paper tries to use dynamic opposite learning (DOL) to significantly enable the diversity of solutions in the SCA, where the elitist strategy is also used to accelerate the convergence. The DOL strategy was first proposed by Xu et al. [50], which significantly improved the comprehensive performance of TLBO . Combining DOL and elite policies with SCAs, a new EDOLSCA is proposed, taking the advantages of the adopted strategies to speed up SCA search and improve precision for specific problems. The ideas in this paper are evident in the CEC2014 benchmark test and PFHE design issues.

The major contribution of this article is summarized as follows:

Firstly, the objective function of PFHE optimization design is derived by using the method of number of transfer units (NTU), and two design objectives and fluid parameters are given.

Secondly, elite strategy and DOL are adopted to strengthen SCA, generating the EDOLSCA. The improvement makes the convergence speed of SCA greatly increase.

Finally, EDOLSCA was tested and compared using the CEC2014 benchmark and then adopted for solving an engineering example of PFHE optimal design.

The structure of this paper is as follows:

Section 2 demonstrates the thermodynamic formula of PFHE, establishes the optimization objective and evaluation function, as well as the constraint conditions. Sections 3 and 4 elaborate the preliminaries of the elite and DOL strategies and illustrate the principle and specific algorithm scheme of EDOLSCA. Section 5 provides a comprehensive test of EDOLSCA, including tests based on the CEC2014 benchmark and PFHE design, with comprehensive analysis of the results. Section 6 summarizes the whole paper and outlooks the future directions.

2 Previous work on plate-fin heat exchanger optimization

A standard two-flow straight finned plate-fin heat exchanger is shown in Fig.1. The plate-fin part of the plate-fin heat exchanger is shown in Fig.(a), in which the design variables mentioned below are marked. The fins are generally welded from thin aluminum plates. The picture shows a relatively simple form of fin. Figure(b) is the working principle diagram of the plate-fin heat exchanger. The two fluids flow alternately in different fins, exchanging heat with maximum efficiency. The optimization work in this paper is based on this featured heat exchanger. In this section, the formula is derived first; then, the optimization constraints and objectives are determined.

Fig. 1
figure 1

Plate-fin heat exchanger and rectangular offset strip fin

2.1 Formula derivation for plate-fin heat exchanger

Heat duty Q is defined first, then Q defines the overall energy exchange level of the heat exchanger, and the formula is given by Yousefi et al. [53]:

$$\begin{aligned} Q= & {} \epsilon C_{min} \left( T_{Ain}-T_{Bin} \right) \end{aligned}$$
(1)
$$\begin{aligned} C_{min}= & {} m*C_{p} \end{aligned}$$
(2)

where C is the heat capacity rate, \(T_{Ain},T_{Bin}\) refer to fluid A (usually hot fluid) and fluid B (usually cold fluid), in refers to the fluid flowing into the heat exchanger, and min refers to the less numeric of the two fluids, max is the larger of the two fluids.

The effectiveness \(\epsilon\) is given by Incropera et al. [20],

$$\begin{aligned} \epsilon =1-e^{\left( \left( \frac{1}{C_{r}} \right) NTU^{0.22} \left( e^{\left( -C_{r} NTU ^{0.78}\right) }-1\right) \right) } \end{aligned}$$
(3)

where e is the Eulers number,

$$\begin{aligned} C_{r} = C_{min}/ C_{max} \end{aligned}$$
(4)

Heat transfer coefficient method is a common method of thermodynamic solution, which is usually used to solve problems without knowing outlet temperature. The heat transfer coefficient is expressed as:

$$\begin{aligned} \begin{aligned}&\frac{1}{NTU}=\\&C_{min}\left( \frac{A_{ffA}}{j_{A}Cp_{A}Pr_{A}^{-0.667}m_{h}A_{A}} + \frac{A_{ffB}}{j_{B}Cp_{B}Pr_{B}^{-0.667}m_{B}A_{B}}\right) \end{aligned} \end{aligned}$$
(5)

where j is Colburn factor, \(A, A_{HT}\) is heat transfer area (m\(^{2}\)), Cp is specific heat of fluid (W\(\cdot\) kg\(^{-1}\) \(\cdot\)K\(^{-1}\) ), Pr is Prandtl number, and \(N_{s}\) is number of entropy generation units.

The free flow area of two fluids is defined as:

$$\begin{aligned} Aff_{A}= & {} \left( H_{t}-t_{A} \right) \left( 1-n_{A} t_{A} \right) L_{B} N_{A} \end{aligned}$$
(6)
$$\begin{aligned} Aff_{B}= & {} \left( H_{B}-t_{B} \right) \left( 1-n_{B} t_{B} \right) L_{A} N_{B} \end{aligned}$$
(7)

In the formula, H is height of the fin, H is height of the fin, t is fin thickness, n is fin frequency, L denotes heat exchanger length, N is number of fin layers, \(\mu\)is dynamic viscosity (N \(\cdot\) s \(^{-1}\) \(\cdot\) m \(^{-2}\)), and \(\rho\) is density (kg \(\cdot\) m \(^{-3}\)).

Fluid B has one more layer than fluid B:

$$\begin{aligned} N_{B}=N_{A} + 1 \end{aligned}$$
(8)

The heat transfer area of the fluid comes from Yousefi [52]:

$$\begin{aligned} A_{A}= & {} L_{A} L_{B} N_{A} \left( 1 + \left( 2n_{A}\left( H_{A} -t_{A}\right) \right) \right) \end{aligned}$$
(9)
$$\begin{aligned} A_{B}= & {} L_{A} L_{B} N_{A} \left( 1 + \left( 2n_{B}\left( H_{B} -t_{B}\right) \right) \right) \end{aligned}$$
(10)

j is the Colburn factor and f is fanning friction factor [23].

For laminar flow \(\left( Re \le 1500 \right)\),

$$\begin{aligned} j= & {} 0.53(Re)^{-0.5}\left(\frac{l}{D_{h}}\right)^{-0.15}\left(\frac{s}{H}-t\right)^{-0.14} \end{aligned}$$
(11)
$$\begin{aligned} f= & {} 8.12(Re)^{-0.74}\left(\frac{l}{D_{h}}\right)^{-0.41}\left(\frac{s}{H}-t\right)^{-0.02} \end{aligned}$$
(12)

For turbulent flow\(\left( Re > 1500 \right)\),

$$\begin{aligned} j= & {} 0.21(Re)^{-0.4}\left(\frac{l}{D_{h}}\right)^{-0.24}\left(\frac{s}{H}-t\right)^{0.02} \end{aligned}$$
(13)
$$\begin{aligned} f= & {} 1.12(Re)^{-0.36}\left(\frac{l}{D_{h}}\right)^{-0.65}\left(\frac{s}{H}-t\right)^{0.17} \end{aligned}$$
(14)

where Re is Reynolds number, s is fin spacing.

The fin spacing of the PFHE:

$$\begin{aligned} s=(1/n-t) \end{aligned}$$
(15)

where

$$\begin{aligned} Re= \frac{m\cdot D_{h}}{A_{ff}\mu } \end{aligned}$$
(16)

For the fin shape in this paper, the hydraulic diameter \(D_{h}\) of the heat exchanger can be expressed as:

$$\begin{aligned} D_{h}=\frac{4s\left( H-t \right) l}{2\left( sl + \left( H-t \right) l + t\left( H-t \right) \right) + ts} \end{aligned}$$
(17)

where \(D_{h}\) is hydraulic diameter and l is lance length of the fin.

The pressure drop due to friction between the fluid in the channel and the fin can be calculated as:

$$\begin{aligned} \Delta P_{A}= & {} \frac{2 f_{A} L_{A}\left( \frac{m_{A}}{Aff_{A}} \right) ^{2}}{\rho _{A} D_{h,A}} \end{aligned}$$
(18)
$$\begin{aligned} \Delta P_{B}= & {} \frac{2 f_{B} L_{B}\left( \frac{m_{B}}{Aff_{B}} \right) ^{2}}{\rho _{B} D_{h,B}} \end{aligned}$$
(19)

where \(\Delta P\) is pressure drop (N \(\cdot\)m\(^{-2}\))

2.2 Objective function for plate-fin heat exchanger

For most heat exchangers, the flow of fluid in the channel is accompanied by friction, which leads to the process of entropy increase. The design of the heat exchanger determines the amount of entropy increase. Therefore, entropy increase is considered as an effective evaluation index in the design of heat exchanger. Under the same fluid condition, the smaller the number of entropy increases, the better the performance of the heat exchanger.

In this article, one of the target functions is defined as entropy generation units [5]:

$$\begin{aligned} \begin{aligned} Ns =&\left( 1-\epsilon \right) \left[ \frac{\left( T_{Bin}-T_{Ain} \right) ^{2}}{T_{Bin}T_{Ain}} \right] + \left( \frac{Rcte_{A}}{Cp_{A}} \right) \left( \frac{\Delta P_{A}}{P_{Ain}} \right) \\& + \left( \frac{Rcte_{B}}{Cp_{B}} \right) \left( \frac{\Delta P_{B}}{P_{Bin}} \right) \end{aligned} \end{aligned}$$
(20)

where Rcte is specific gas constant (J\(\cdot\) kg\(^{-1}\) \(\cdot\)K\(^{-1}\))

In addition to entropy generation, the effective heat transfer area of heat exchanger is also one of the important indexes to evaluate heat exchanger. The heat transfer area is also closely related to the Reynolds number and efficiency of the heat exchanger. Therefore, the second objective function in this paper is the heat transfer area of the heat exchanger

The total heat transfer area is:

$$\begin{aligned} A_{HT}=A_{A} + A_{B} \end{aligned}$$
(21)

2.3 Application example of plate-fin heat exchanger

For the two PFHE design examples prepared in this paper, the difference lies in the difference of fluid parameters and size constraints. The specific constraints and fluid data are shown in Table 1.

In order to meet the heat duty requirement of the heat exchanger, effective constraints are carried out in the calculation process of the algorithm, and the penalty function is used to achieve this requirement. The formula is as follows:

$$\begin{aligned} g\left( X \right) \Rightarrow \xi \left( X \right) -Q=0 \end{aligned}$$
(22)

In the design case, the overall size of the heat exchanger, the density of fins, the length, width and height of fins, the number of fluid layers of the two fluids and the thickness of fins are all the constraints of the design. The volume of the heat exchanger is limited to 1m*1m, and heat duty is 160kW of Case 1 and 1050kW of Case 2, respectively. In the calculations, fluids are considered ideal gases.

Table 1 Design constraints and fluid data

3 Algorithm preliminaries

Heuristic optimization algorithms are inspired by biology as well as natural sciences, such as physics and mathematics. Different ideas of particle updating determine the different characteristics of each algorithm. For different algorithm ideas, the parameters that need to be adjusted are also different. Some algorithms do not need to adjust parameters, while others may need to adjust 4–6 parameters to ensure the best results.

3.1 Sine cosine algorithm

For SCA, the main idea of the algorithm is the sine and cosine functions in mathematics. Four random numbers determine the motion direction and motion distance of particles, and there are no parameters to be adjusted in the algorithm. Therefore, SCA has relatively strong adaptability, which provides the possibility for wide application.

The core of the algorithm is as follows:

$$\begin{aligned} X_{i}^{t + 1} = \left\{ \begin{matrix} X_{i}^{t} + r_{1}\times sin(r_{2}) \times \left| r_{3} P_{t}^{i}-X_{i}^{t}\right| , &{} r_{4} < 0.5 \\ X_{i}^{t} + r_{1}\times cos(r_{2}) \times \left| r_{3} P_{t}^{i}-X_{i_{i}}^{t}\right| , &{} r_{4} \ge 0.5\\ \end{matrix}\right. \end{aligned}$$
(23)
$$\begin{aligned} r_{1}=a-t\frac{a}{T} \end{aligned}$$
(24)

where X is location of position, t is the \(t-th\) iteration , P is final solution, and T is maximum number of iterations. Since the main updating method in the formula involves sine and cosine function, it is named sine and cosine algorithm. The rules of the algorithm are shown more intuitively in Figs.2, 3. Figure 2 shows how the parameters in the equation affect the search space for the solution, which is a two-dimensional representation. Figures 2 and 3 together show the influence of sines and cosines on the next position in the calculation. The search space can be achieved by changing the range of values of sines and cosines. And the random position in space or out of space is also realized by equation parameters.

Fig. 2
figure 2

Rules for the direction of particle motion in SCA

Fig. 3
figure 3

Relationship between particle motion and sine and cosine functions in SCA

In SCA, \(r_{1}\) determines the direction of movement, \(r_{2}\) determines the distance of movement, ranging in [0,2\(\pi\)], \(r_{3}\) determines the direction whether to move inward or outward, bringing randomness to the position, ranging in [-2,2], \(r_{4}\) determines whether to move by sine or cosine manner, ranging in [0,1].

3.2 Dynamic-opposite learning

In some featured cases, optimization algorithm can be improved effectively by modifying certain steps, such as convergence speed and search accuracy. Oppositional-based learning (OBL) [35] is a popular modification method in recent years, from which quasi-opposite-based learning (QOBL) [34] and quasi-reflection-based learning (QRBL) [15] are evolved.

OBL was first proposed to accelerate the convergence of DE. A strategy for finding the opposite number is proposed:

$$\begin{aligned} X^O=a + b-X \end{aligned}$$
(25)

where \(X^{O}\) is the opposite X, \(X\in [a,b]\).

In practical application, if D becomes the dimension of the problem, Eq.25 can be rewritten as:

$$\begin{aligned} X_j^O=a_j + b_j-X_j, j=1:D \end{aligned}$$
(26)

In practical application, OBL, QOBL and QRBL are all likely to fall into the local optimal position in the search space. Xu et al. [50] proposed the dynamic opposite learning (DOL) strategy to solve this problem:

$$\begin{aligned} X^{DO}=X + w*rand*(rand*X^O-X) \end{aligned}$$
(27)

where ’rand’ is a random number in [0, 1] , w is the weight that determines the combination of DOL and the original algorithm, that is, the weight of balance, and the ’\(rand*X^O\)’ makes the search space become asymmetric, so as to solve the problem of falling into local optimal, and make the search become dynamic and constantly changing.

In the same way, DOL is applied to the real world and dimension D is added, which can be rewritten as:

$$\begin{aligned} X_j^{DO}=X_j + w*rand*(rand*X_j^O-X_j), j=1:D \end{aligned}$$
(28)

Applying the above DOL strategy to the initial population and iteration makes the algorithm search better and faster.

4 Elite and DOL-based SCA

In this section, we describe in details how EDOLSCA works. EDOLSCA combines elite and DOL policies to bring greater speed and precision to SCA.

4.1 DOL population initialization

After the population is initialized, the first step will be calculated:

$$\begin{aligned} oP_{ij}^{DO}=oP_{ij} + r1_i*(r2_i*(a_j + b_j-oP_{ij})-oP_{ij}) \end{aligned}$$
(29)

where oP is a random population generated within the upper and lower bound by using random numbers, and i and j are the \(j^{th}\) variables of the \(i^{th}\) individual. a and b are the boundaries of the variables, and \(r_{1}\) and \(r_{2}\) are random numbers. This step can improve the adaptability of the improved initial population through initial processing of the population. After modifying the population, the particle should be checked to see whether it exceeds the set boundary.

4.2 DOL generation jumping

After the DOL processes the initial population, the algorithm enters the SCA processing phase and is updated according to the SCA rules. In order to cooperate with different algorithms, jumping rate (Jr) is introduced into DOL as the basis for whether the algorithm enters into DOL step. Jr is a number in the range of [0, 1] , and the specific value should be determined by experiments. In the algorithm step, the program randomly generates a random number in the range of [0, 1] . If the random number is less than the set Jr, then the algorithm will enter the DOL step:

$$\begin{aligned} P_{ij}^{DO}=P_{ij} + w*r3_i*(r4_i*(a_j + b_j-P_{ij})-P_{ij}) \end{aligned}$$
(30)

The boundary between A and B will be dynamically updated:

$$\begin{aligned} \begin{aligned}&a_j=min(P_{ij})\\&b_j=max(P_{ij}) \end{aligned} \end{aligned}$$
(31)

As with the other steps, the updated population needs to check the boundaries.

4.3 EDOLSCA steps

In addition to the DOL strategy, the introduction of elites is intended to address the exploration speed. The concept of elite is to select a part of the best individuals before each iteration, skipping the intermediate steps and keeping until the end of each iteration to form a new population, so as to ensure the overall superiority of the population.

At this point, the specific process of EDOLSCA has been introduced, and the flowchart of the algorithm is shown in Fig. 4, and steps of EDOLSCA are shown in Algorithm.1.

Fig. 4
figure 4

The flowchart of EDOLSCA

figure a

5 Experimental results and discussion

5.1 Benchmark test

In this test, the benchmark proposed by CEC2014 [27] was selected as the test function. The 14 alternative benchmarks are shown in Table 2, the first four benchmarks are unimodal functions, and the remaining twelve are multi-modal functions. These benchmarks are effective to test the exploration and development capabilities of algorithms. In order to make the comparison of algorithms more intuitive, this paper selects four standard and classical algorithms with improvements, namely, PSO [25], cfPSO [10], cfwPSO and ETLBO [24]. Algorithm parameters are as follows: for PSO, \(C_{1}=C_{2}=2\); for cfPSO and cfwPSO, \(C_{1}=C_{2}=2.05\), \(K=0.729\). For the parameters of EDOLSCA, w range is 1–15; Jr range is 0.1–1. The algorithms run 10 times for independent testing, and the results are shown in Tables 3, 4. The average with the smallest index is chosen, that is \(w=10\) and \(Jr=1\). The population size in the test was set to 50, and the number of iterations was 20,000. Results of the mean and standard deviation of each algorithm at each benchmark are recorded to compare the algorithms. Average index is considered as one of the criteria of comprehensive evaluation algorithm, which is calculated from the ranking of each algorithm. The calculated results are shown in Table 5, and the convergence diagram is shown in Fig.5.

Table 2 Unimodal and multimodal benchmark
Table 3 The sensitivity analysis of w
Table 4 The sensitivity analysis of Jr
Fig. 5
figure 5

The convergence trends of all algorithms on all unimodal and simple multimodal functions

Table 5 The mean and standard value of unimodal/multi-modal test functions

First of all, it can be seen from the calculation results in the table that, for F1-F4, three of the four baselines EDOLSCA reach the optimal value, and in F5-F14, most EDOLSCA also get the minimum value of all algorithms. ’Ave index’ shows that EDOLSCA has the best comprehensive result among all comparison algorithms, and only the improved SCA in one benchmark performs less well than the original algorithm. It is observed in the figure that EDOLSCA converges much faster than SCA in most benchmarks, which proves that the improvement is successful. In addition to using the average index method to detect EDOLSCA, the Wilcoxon rank-sum test [47] method is also introduced to analyze the results of the algorithm. The significance is set to 0.05, two-tailed. The test results are shown in Table 6.

Table 6 Wilcoxon rank-sum test on unimodal and multi-modal functions

In the result table, ’\(\uparrow\)’ denotes the EDOLSCA is better than other counterpart, ’\(\leftarrow\)’ means other algorithms (the algorithm pointed by the arrow) is better, and ’=’ means there is no significant difference between the two algorithms after the test. The results show that the convergence accuracy has a small improvement compared with SCA, but the overall improvement ratio is not high. It shows that EDOLSCA has an absolute advantage in convergence speed, but it can still be improved in terms of accuracy.

5.2 Application to PFHE design problems

In this section, the optimization object is extended to the design problem of PFHE. In fact, the design optimization problem of PFHE is different from the repeatability optimization. The design optimization problem only needs an optimal design, and the stability of the repeatability optimization is more important [54]. Therefore, the stability of the algorithm is not considered in the PFHE problem, where only the accuracy of the algorithm is considered. For the heat transfer area objective function, both cases 1 and 2 are applicable. In order to optimize the environment realistically, the number of iterations is set to 500,000, and the population is still 50. For the objective function of entropy increase, only test case 1, the number of iterations is set to 4,000, and the population is 50. In order to compare the fairness of the study, we add two varieties of algorithms in this section: elite sine cosine algorithm (ESCA), dynamic opposite learning sine cosine algorithm (DOLSCA) and two common algorithms: differential evolution (DE) [41] and neighborhood field for cooperative optimization (NFO) [48]. The results are shown in Figs.6,  7.

Fig. 6
figure 6

The solution result of the AHT objective, cases 1 and 2

Fig. 7
figure 7

The solution result of the Ns objective, cases 1

It can be seen from the result image that SCA and EDOLSCA have the same advantage in the AHT objective function. The EDOLSCA continued to converge under a high number of iterations, and the improved algorithm has higher accuracy than the original algorithm. On the Ns objective function, EDOLSCA converges faster than SCA, and the accuracy difference is larger. At the same time, in the three cases, the accuracy of EDOLSCA is the highest, which proves its potentials in engineering design optimization problems. The results of benchmark and engineering tests show that EDOLSCA has excellent dynamic search capabilities and has unique advantages in speed and accuracy, as well as demonstrating excellent comprehensive performance.

6 Conclusion

This paper proposed a new EDOLSCA adopted elite and DOL strategies to strengthen the canonical SCA, which greatly improves the performance of the original algorithm in terms of convergence accuracy and convergence speed. The proposed method was compared with the state-of-the-art algorithms counterparts in the CEC2014 benchmark test proved to perform well. The algorithm is further adopted in solving the optimal design problem of PFHE, where the dynamic and random search space of the DOL strategy has achieved excellent results. The proposed algorithm again demonstrated superiority in solving the complex optimal design problem. Future work lies on the wide application of the proposed EDOLSCA in various industrial and science optimization scenarios. The well robustness of the algorithm is promising to provide a useful computational tool for engineering optimization.