Next Article in Journal
Simulation of Groundwater Dissolved Organic Carbon in Yufu River Basin during Artificial Recharge: Improving through the SWAT-MODFLOW-RT3D Reaction Module
Previous Article in Journal
A Protocol of Prevention and Protection Measures on New Occupational Risk Factors in Green Jobs in Italian Workplaces
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Development and Sensitivity Analysis of an Improved Harmony Search Algorithm with a Multiple Memory Structure for Large-Scale Optimization Problems in Water Distribution Networks

1
Water & Wastewater Research Center, K-Water Research Institute, Yuseong-gu, Daejeon 34045, Republic of Korea
2
Department of Mechanical Engineering, Faculty of Engineering, University of Science and Culture, Tehran 13145-871, Iran
3
Department of Civil and Infrastructure Engineering, Gyeongsang National University, Jinju 52725, Republic of Korea
4
Department of Civil and Environmental Engineering, Dongshin University, Naju 58245, Republic of Korea
5
Department of Civil Engineering, University of Suwon, Hwaseong-si 18323, Republic of Korea
*
Authors to whom correspondence should be addressed.
Sustainability 2024, 16(15), 6689; https://doi.org/10.3390/su16156689
Submission received: 29 May 2024 / Revised: 23 July 2024 / Accepted: 28 July 2024 / Published: 5 August 2024
(This article belongs to the Section Sustainable Engineering and Science)

Abstract

:
The continuous supply of drinking water for human life is essential to ensure the sustainability of cities, society, and the environment. At a time when water scarcity is worsening due to climate change, the construction of an optimized water supply infrastructure is necessary. In this study, an improved version of the Harmony Search Algorithm (HSA), named the Maisonette-type Harmony Search Algorithm (MTHSA), was developed. Unlike the HSA, the MTHSA has a two-floor structure, which increases the optimizing efficiency by employing multiple explorations on the first floor and additional exploitations of excellent solutions. Parallel explorations enhance the ability in terms of exploration (global search), which is the tendency to uniformly explore the entire search space. Additional exploitations among excellent solutions also enhance the ability of local searches (effective exploitation), which is the intensive exploration of solutions that seem to have high possibilities. Following the development of the improved algorithm, it was applied to water distribution networks in order to verify its efficiency, and the numerical results were analyzed. Through the considered applications, the improved algorithm is shown to be highly efficient when applied to large-scale optimization problems with large numbers of decision variables, as shown in comparison with the considered optimizers.

1. Introduction

Following the introduction of the Genetic Algorithm (GA) [1], one of the most widely known meta-heuristic algorithms, Tabu Search (TS) [2], Simulated Annealing (SA) [3], Ant Colony Optimization (ACO) [4], Particle Swarm Optimization (PSO) [5], and Differential Evolution (DE) [6] were developed, mimicking natural phenomena, human behaviors, and physical phenomena. Subsequently, the Harmony Search Algorithm (HSA) was developed in 2001 [7,8], and many meta-heuristic algorithms followed, such as the Honeybee Algorithm (HBA) [9], the Firefly Algorithm (FA) [10], Cuckoo Search (CS) [11], the Bat Algorithm (BA) [12], the Water Cycle Algorithm (WCA) [13], and the Mine Blast Algorithm (MBA) [14]). Individual algorithms have been continuously enhanced through improvements, and recently they have actively been improved through combinations of different meta-heuristic algorithms or combinations with algorithms with different mathematical bases. Although many meta-heuristic algorithms have been developed over the last 30 years through diverse approaches, the ultimate objective of meta-heuristic algorithms is to drive general methodologies that will obtain efficient results close to global solutions in all optimization problems. In optimization problems using meta-heuristic algorithms, striking the balance between explorations/global searches in unknown areas and searches in obtained information/effective exploitation is important. Exploration means to the tendency uniformly to explore the entire search space, and exploitation is the capacity for intensive explorations of points that seem to have high possibilities.
Following the development of the HSA in 2001, studies to revise or improve the algorithm were initially conducted, such as changing the parameters (Harmony Memory Consideration Rate; HMCR, Pitch-Adjusting Rate; PAR) of the algorithm (Improved/Revised/Modified HSA: [15,16]). Subsequently, some concepts were developed relating to Self-adaptive HSA/Parameter Setting Free HSA, which change the parameters of the algorithm automatically within the algorithm, in order to improve its practical usability [17,18,19]. In addition, some Hybrid HSAs [20,21,22,23,24,25] that combine other optimization algorithms with the existing algorithm to improve its efficiency have recently been developed. Furthermore, the concept of multi-objective HSAs, which can have two or more purposes instead of just one, has been presented [26]. In particular, studies regarding the improvement of the HSA itself have recently been actively conducted. Som researchers [27,28,29,30] have presented improved methods for the application of PAR and bandwidth (BW), which are parameters of the HSA. New method for determining Harmony Memory Size (HMS) and HMCR values, along with efficient results for optimization is presented [31], and an improvement to the HSA is suggested [32], using concepts termed Multiple HMCRs and Dynamic PARs to create new harmony solution vectors. In detail, a method for implementing the practice of individual players first and then optimizing final harmonies among players, named this algorithm Melody Search (MS) is presented [33]. However, these improved HSAs only focus on improving optimization performance by enhancing parameter settings or algorithm operators. Beyond this, it is necessary to improve the characteristics of the memory structure of the HSA, to improve its exploration and exploitation performance.
The efficiency of the development and improvement of meta-heuristic algorithms has been verified through the application of the algorithms not only to mathematical problems but also to a diverse range of engineering problems. Among them, the optimal design of water distribution networks (WDNs) is a highly non-linear and constrained problem to which such algorithms have been applied by many researchers, and it has been widely employed to prove the efficiency of optimization algorithms. In general, the optimal design of WDNs is aimed at finding a low-cost pipe diameter combination that satisfies the hydraulic conditions (e.g., minimum nodal pressure) of WDNs. In the early stages, many researchers focused on basic meta-heuristic algorithms, such as the GA [34,35,36], ACO [37,38], SA [39], the shuffled frog-leaping algorithm (SFLA) [40], and the HSA [7], to overcome the drawbacks of mathematical techniques and generate good results. However, even though there have been many theoretical studies on the optimal design of water distribution networks using optimization techniques, the application to real-world problems is currently challenged by inefficiency and time-consuming tasks of the algorithm. Recently, in order to enhance the performance of the optimization results, many researchers have developed hybrid/improved methods. For instance, Geem [41] utilized HSA coupled with the concepts of PSO in order to form a hybrid method known as Particle Swarm Harmony Search (PSHS). In addition, Sedki and Ouazar [42] combined PSO with differential evolution (DE), which is named PSO-DE. When using the hybrid PSO-DE compared against PSO and other optimizers, they found minimum values of construction costs were reached through lesser computational efforts. Sadollah et al. [43] applied an Improved Mine Blast Algorithm (IMBA) to various benchmark WDNs (i.e., Two loop, Hanoi, New York, and Balerma networks), and they presented considerably more efficient optimization results compared to the existing results. Lee et al. [44] proposed an optimization method considering engineering knowledge of the WDNs, and they applied this to the Pescara network and the Balerma network design problems. Schwartz et al. [45] proposed a least-cost design method for WDNs under multiple loading conditions, and El-Ghandour and Elbeltagi [46] applied meta-heuristic algorithms, including the GA, PSO, and ACO, to three WDNs (Two Loop, New York, and El-Mostakbal networks) and compared the performance of these algorithms. Lee et al. [47] proposed a meta-heuristic algorithm motivated by a vision correction procedure, and they applied this to the Balerma network design problem.
However, among the WDNs that were applied in previous studies, those other than the Balerma network had small numbers of pipes (numbers of decision variables). In addition, the Balerma network took the form of a typical agricultural irrigation network. Going forward, diverse types of water distribution networks should additionally be considered in order to compare the efficiency levels of algorithms. Yoo et al. [48] implemented pipe diameter optimization of Saemangeum networks, considering diverse supply types (branched, looped, and pumped type), and compared and analyzed the optimization results by type with existing designs.
In the present study, a Maisonette-type Harmony Search Algorithm (MTHSA) was developed, based on the Harmony Search Algorithm (HSA). Using this MTHSA, optimal pipe diameter design problems, that is, finding low-cost pipe diameter combinations within limited hydraulic conditions, were solved. In particular, the MTHSA was applied to large-scale problems, and the effects of this application were analyzed, with the results presented. The first benchmark problem selected in this study was the Balerma network minimum-cost design, which is a representative example of an existing engineering optimization problem and is continuously being studied by various researchers. In this study, a sensitivity analysis was performed to evaluate the performance of the parameter combination of the newly developed optimization algorithm for the optimal design of the Balerma network, a representative optimization problem. In addition, in order to verify whether the developed optimization algorithm can be universally applied to pipe network optimization problems, the Saemangeum networks were additionally selected from among the existing optimization benchmark problems that have different characteristics from the Balerma network.

2. Problem Descriptions

2.1. Formulation of Optimal Pipe Design

In general, the minimization of pipe construction costs is considered an objective function, and the pressure and velocity requirements are used as a constraint. Equations (1)–(3) show the objective function with penalty functions.
M i n . C o s t = i = 1 n p C ( D i ) L i + j = 1 n n P j + j = 1 n p V j
P j = α h j h min / max + β
V j = γ v j v min / max + δ
In the above, C D i denotes costs according to pipe diameters per unit length. Li, Di, hj, hmin/max, Pj, vj, vmin/max, Vj, np, and nn denote the pipe length, pipe diameter, nodal pressure at node j, minimum/maximum nodal pressure requirements, penalty function for checking the pressure constraint, link velocity at pipe j, minimum/maximum pipe velocity requirements, penalty function for checking the velocity constraint, number of pipes, and number of nodes, respectively. α, β, γ, and δ are constants in the penalty functions.
In addition, in order to obtain the nodal pressure, hydraulic simulations of the water distribution network should be performed. In this study, the proposed optimization method was coupled with a hydraulic simulator EPANET [49] for optimization purposes. This is a free-access and widely used water distribution network solver, which employs gradients to determine the mass and energy conservation equations [42]. A pipe network design problem was solved by choosing the least-cost pipe diameter among commercial candidate pipes satisfying hydraulic flow demands, energy conservation, and head restriction. The mass conservation of node, energy conservation of loop, minimum/maximum pressure at node, minimum/maximum velocity in pipe, and pipe diameter size availability constraints can be defined by the following Equations (4)–(11):
Q i n Q o u t = Q e
Here, Qin, Qout, and Qe denote flow into the node, flow out of the node, and the external outflow rate (demand) at the node, respectively.
K L o o p l Δ H k = 0 , l N L
Here, ∆Hk is the head loss in pipe k and NL is the total number of closed loops in the system. Equation (5) means that the sum of head losses considering the flow direction of individual pipes within the loop is “0”. Head loss always occurs depending on the flow of water in individual pipes within the loop, but the loss depending on the flow direction appears as a positive or negative value. The head loss for each pipe is the head difference between inter-connected nodes, which can be calculated using head loss equations such as the Hazen–Williams, Darcy–Weisbach, and Manning equations, as shown in Equations (6)–(8).
H k = 10.67   L D 4.8704 Q 1.852 C 1.852
Here, Q is the volumetric flow rate, m3/s (cubic meters per second), L is length, D is diameter, and C is the pipe roughness coefficient in the Hazen–Williams equation.
H k = f L D V 2 2 g
Here, Hk is the head loss, f is the friction factor, L is length, D is diameter, V is velocity, and g is the gravity coefficient in the Darcy–Weisbach equation.
V = k   n R 2 3 ( H k L ) 1 / 2
Here, V is the cross-sectional average velocity, n is the Gauckler–Manning coefficient, R is the hydraulic radius, and L is length in Manning equations.
hmin ≤ hj ≤ hmax j = 1, 2, 3, …, nn
Here, hj, hmin, hmax, and nn denote the pressure head at node j, the minimum and maximum required pressure heads at node j, and the number of nodes in the network, respectively.
vmin ≤ vj ≤ vmax j = 1, 2, 3, …, np
Here, vj, vmin, vmax, and np denote the velocity at pipe j, the minimum and maximum required velocities at node j, and the number of pipes in the network, respectively.
D i D i n p
Here, Di is the diameter of pipe I, and {D} denotes the set of commercially available pipe diameters. In fact, all design variables (i.e., pipe diameter) should be chosen from the suggested commercial sizes (discrete optimization problem).

2.2. Studied Networks

2.2.1. Balerma Network

The Balerma network has four reservoirs, eight loops, 454 pipes, and 443 demand nodes, as shown in Figure 1. It was firstly presented by Reca and Martinez [50]. The minimum nodal pressure requirement is 20 m at each node, and the pipe velocity requirements are not considered for the Balerma network.
The discrete sizes of commercial diameters, accompanied by their respective cost per unit length, can be found in Table 1. Table 1 shows that each pipe can be selected from the 10 candidate diameters, resulting in a total of 10,454 potential network configurations.

2.2.2. Saemangeum Networks

Saemangeum networks were initially introduced by Yoo et al. [48] as an irrigation water supply network in the Republic of Korea. While originally planned as a looped network, the Saemangeum network was also considered as branched and pumped types in this study. The branched network (Plan 1) consists of 335 pipes with a combined length of 37.88 km, shown in Figure 2a. On the other hand, the looped network (Plan 2) has 356 pipes; we converted some pipe configurations and generated 10 more loops from the branched network, as shown in Figure 2b. Furthermore, the pumped irrigation network is made up of 345 pipes that are supplied through four pumping stations within the network (Plan 3, Figure 2c). The length of the pipes is 41.39 km, and based on four pumping stations, the whole network is divided into four sections.
The data regarding the pipe cost per unit length for the pipes of different diameters used in this study are shown in Table 2. In total, 18 types of commercial pipes were considered in the optimization process.
These were obtained from the ‘Water Facilities Construction Cost Estimation Report’, from K-water [51] in the Republic of Korea, which provides cost data on the construction costs for different types of pipes [48]. In addition to nodal pressure, flux conditions of the pipe line were also considered as hydraulic conditions for the Saemangeum networks. The reason why these water pressure and water flux standards exist in pipe network optimization problems is because water supply stability can be problematic if the water pressure is low or too high, or if the water flux is slow or very fast. Therefore, the purpose of optimal water pipe network design is to determine the pipe diameter that supplies water at the minimum cost within the appropriate water pressure and water flux range. In the case of nodal pressure, the minimum value and the maximum value are set as 10 m and 35 m, respectively. In the case of flux conditions of the pipe line, values should be maintained at between 0.01 m/s and 2.5 m/s.

3. Maisonette-Type Harmony Search Algorithm

3.1. Description of the Proposed Algorithm

A maisonette is defined as an apartment that is usually part of a larger building with two levels, often having its own entrance from the outside. The multi-floor structure has the advantage of making efficient use of space, because two floors can be utilized, although they are independent residential units. That is, the first floor and the second floor have the characteristic of being independent from each other, while still being mutually linked. Accordingly, the MTHSA in the present study, based on this structure, was developed as an optimization technology. Figure 3 shows the structure and components of the MTHSA. That is, multiple rooms can be arranged on the first floor, and independent spaces different from those on the first floor can be separately organized on the second floor. Therefore, in the case of the MTHSA, multiple Sub Harmony Memories (SHMs) are organized on the first floor. On the second floor, separate spaces are formed. Here, excellent solutions created in individual rooms on the first floor, and additional HSAs that can be implemented through the solutions, can be separately stored. That is, unlike existing HSAs, this structure increases optimization efficiency through multiple (parallel) searches on the first floor and additional searches among excellent solutions on the second floor. In the MTHSA, the Harmony Memory is divided into separated Sub Harmony Memories, and an independent search is performed in each Sub Harmony Memory, thereby improving the local search performance (exploitation). In addition, based on the best solutions of Sub Harmony Memories on the first floor, the solution memory of the second floor is formed, and further operations are performed, thereby improving the global search performance (exploration) in the MTHSA.

3.2. Model Parameters

In the present study, five parameters were used. HMS, HMCR, PAR, and BW were used as with HSAs, and Sub Harmony Memory Size (SHMS) was additionally employed. HMS means the total number of Harmony Memories used to solve optimization problems. The HMCR is a user parameter that shows the probability of randomly generating new solutions, or selecting new solutions from existing Harmony Memories. If the value of the HMCR is 0.7, this means that the probability of generating new solutions from existing HMs is 70%, and the probability of random selection of new solutions is 30%. The PAR is the probability of adjusting the value of each decision variable of new solutions within the designated range (BW). If the value of PAR is 0.1, the probability of adjusting the value of each decision variable is 10%. The designated range used in this case was determined by the bandwidth. In the case of water distribution network problems, because discrete variables for selecting one of the available sizes of commercial pipes were used, the BW is defined as increasing or reducing pipe diameters by one level. SHMS means the size of rooms arranged on the first floor of the MTHSA. For instance, if the value of HMS is 60 and the value of SHMS is 20, three HMs of size 20 will be formed on the first floor.

3.3. Procedure of the Proposed Algorithm

A pseudo code that shows the process of implementing the MTHSA is shown in Algorithm 1. First, the values of the HMS, SHMS, HMCR, PAR, and BW parameters of the MTHSA are set.
According to the set sizes of HMS and SHMS, the number of rooms on the first floor and the number of excellent solutions to be organized on the second floor are determined.
That is, if the values of HMS and SHMS are 60 and 20, respectively, the number of rooms on the first floor and the excellent solutions to be organized on the second floor will both be three. The next stage is to perform parallel searches on the first floor. Using the three Harmony Memory sets organized on the first floor, HSAs are implemented separately. In this case, the HMCR and PAR are used to create new solutions. The new solutions and the objective function of the worst solution in each set are then compared with each other, and if the values of the new solutions are better, the location of the Harmony Memory will be changed. The next stage is to perform searches among excellent solutions on the second floor. The best solutions in individual sets obtained on the first floor become the components of the initial Harmony Memories on the second floor. The HSA is implemented using the Harmony Memories on the second floor, organized as such. The process of optimization in the MTHSA is iterated until the terminating condition (e.g., the number of iterations) is satisfied. In general, whereas one new solution is generated at each iteration in the HSA, the number of new solutions generated at each iteration varies with the values of HMS and SHMS in the case of the MTHSA. That is, when the values of HMS and SHMS are 60 and 20, respectively, the HSA is implemented three times on the first floor and one time on the second floor in each iteration, so that a total of four new solutions are generated.
Algorithm 1. Pseudo code of MTHSA.
Begin
Objective function f(x), x = (x1,x2,…,xd)T
Generate initial harmonies (Define HMS & SHMS)
Define harmony memory considering rate (HMCR)
Define pitch adjusting rate (PAR), pitch limits and band width (BW)
while (t ≤ Max number of iterations)
Generate new harmonies by accepting best harmonies in first floor
Adjust pitch to get new harmonies (solutions)
if (rand < HMCR), choose an existing harmony randomly in each of sub harmony memories
else if (rand < PAR), adjust the pitch randomly within limits
else generate new harmonies via randomization
end if
Accept the new harmonics (solutions) if better
Gathering best harmonics in each of sub harmony memories
Generate new harmony by accepting best harmony in second floor
Adjust pitch to get new harmony (solution)
if (rand < HMCR), choose an existing harmony randomly
else if (rand < PAR), adjust the pitch randomly within limits
else generate new harmony via randomization
end if
Accept the new harmonics (solutions) if better
end while
Find the current best solutions
End

4. Application Results

4.1. Sensitivity Analysis of Parameters

Before verifying the applicability of the MTHSA, an analysis of the sensitivity of each user parameter was carried out, and the results are presented here. The sensitivity analysis was applied to the Balerma network. In order to have a fair comparison, 30 runs were implemented for each parameter combination. The Number of Function Evaluations (NFE) of all simulations was set to 45,400. Here, the NFEs are general calculation termination conditions for evaluating the performance of each technique in the field of optimization under the same conditions.
In order to analyze the sensitivity of each parameter, among various parameter combinations, a combination of 60 HMSs, 20 SHMSs, a PAR of 0.02, and an HMCR (MTHSA (1)) of 0.90, which produces relatively good results, was set as a reference combination. The results of the sensitivity analyses are presented in Table 3, Table 4, Table 5 and Table 6. Table 3 and Table 4 show the sensitivity analysis results for HMS and SHMS.
It can be observed that the smaller HMS values and the larger SHMS values show better results. One thing that must be noted is that a ratio between HMS and SHMS (HMS/SHMS) that is closer to one shows better results. Since the maximum NFEs have been determined, the larger the value of the ratio between HMS and SHMS, the more times function evaluations are implemented per time of iteration. In other words, those types of the MTHSA that can be advantageous in terms of the number of implementations of function evaluation, by performing searches among excellent solutions on the second floor while reducing the number of solution searches on the first floor, show relatively higher optimization efficiency. However, it is judged that the ratio between HMS and SHMS (HMS/SHMS) that shows the highest efficiency may vary with the value of the Maximum NFEs. Therefore, in the case of the MTHSA, determining appropriate ratios between HMS and SHMS is judged to be an important criterion.
Table 5 and Table 6 show the results of HMCR and PAR sensitivity analysis. It can be observed that relatively larger HMCR values and relatively smaller PAR values show better results. In the case of HMCR, the larger the value, the higher the probability of generating new solutions in the existing Harmony Memory. Therefore, applying HMCR with a large value means generating a new optimal solution on the second floor by using the information of the best solutions of the first floor with a high probability, and therefore the optimal solution can be found faster. It can be seen that in the case of PAR, because the optimization problem involves a large number of decision variables (i.e., 454 decision variables), large PAR values act as an element that obstructs the convergence of solutions. In particular, it can be seen that deviations in optimal solutions determined by the value of PAR appear much larger compared to other parameters. Therefore, in the case of large-scale problems with large numbers of decision variables and alternative solutions, in order to obtain reasonable optimal results, using relatively lower PAR values is judged to be desirable.

4.2. Optimization Results

4.2.1. Results of Balerma Network

Based on the results of the sensitivity analyses conducted Section 4.1, the results for the optimal design of the Balerma network were organized and analyzed, as presented in the following. Table 7 shows the combinations of parameters and the NFE used to compare and analyze the results.
MTHSA (I) and MTHSA (II) denote those parameter combinations that showed excellent results among the sensitivity analysis results. MTHSA (III) shows a case with the same parameters as those used in MTHSA (I), but a different value for the NFEs.
Table 8 presents those results from the existing literature that are different from the results of the MTHSA. The optimal results of the GA, SA, MSATS [36], the HAS (1), HAS (II), PSHS [41], the genetic heritage evolution by stochastic transmission (GHEST) [52], the MBA, IMBA (I), IMBA (II) [43], the decomposition and multistage approach for optimizing the design of WDSs using the DE (Decompose-DE) [53], the GANOME [50], the HD-DDS-1, HD-DDS-2 [54], and other optimization techniques from the literature are presented. Table 8 presents the results, organized in ascending order of the Maximum NFE values used in individual pieces of the literature.
First, in the case of Best Cost, Decompose-DE and NLP-DE2 show the least cost, at a value of 1.923 (106 EUR). It can be observed that the MTHSA (III) developed in the present study shows almost the same value, as it yields a value of 1.932 (106 EUR). In the case of the Decompose-DE method, before applying DE, pre-processing is implemented based on the engineering characteristics of water distribution networks. The pre-processing process used here should be preceded by gaining a sufficient understanding of water distribution networks. Therefore, the MTHSA has an advantage in that it shows sufficiently good results even without pre-processing based on the engineering knowledge of water distribution networks. Unlike the MTHSA, the NLP-DE2 is a combination of non-linear programming and DE, and it confirms that the MTHSA is efficient, as good optimal solutions were obtained when only meta-heuristic techniques were used in the case of the MTHSA.
On reviewing the results, when the value for the Maximum NFEs was fixed at 45,400, it can be seen that the proposed method and the IMBA show excellent results. For the case of the Best Cost, IMBA (I) showed that it can design solutions at a cost lower by 0.021 (106 EUR) than for MTHSA (II). However, the Average and Worst Costs of the proposed algorithm were shown to be lower by approximately 0.030 (106 EUR) and 0.063 (106 EUR), respectively, compared to IMBA (I). This means that the proposed algorithm achieved better optimal solutions on average compared to the IMBA.
When reviewing the results of HSA (I) and HSA (II), which are based on the HSA and PSHS, it can be seen that when a Maximum NFE value of 45,400 was applied, the MTHSA achieved much better results. Whereas a value of 2.018 (106 EUR) could be obtained by applying an NFE of 10,000,000 in the case of HSA (II), a value of 1.932 (106 EUR) was obtained by choosing 1,000,000 as the NFE in the case of MTHSA (III).
Figure 4 shows the results (nodal pressure) of a numerical analysis of WDNs, when the optimal pipe diameter combination from MTHSA (I) is chosen. The node that shows the closest result to the minimum required pressure of 20 m is no. 55 (20.01 m), and it can be seen that this value satisfies the minimum required pressure. In particular, in the case of WDNs, an even pressure distribution is very important in terms of leakage management. Figure 4b shows the cumulative distribution of the pressure of each node, appearing as a result of the optimal design, from which an even distribution without any particular unequally distributed pressure section can be identified.

4.2.2. Results of Saemangeum Network

Yoo et al. [48] applied the HSA to the Saemangeum network for the first time, in order to present optimal design results. In the present study, the results of the application of the proposed MTHSA algorithm were compared with the results of Yoo et al. [48] and the results of that comparison are presented. The parameters used in the proposed algorithm were an HMS of 60, SHMS of 20, PAR of 0.03, and HMCR of 0.90. Note that the parameters were determined through a sensitivity analysis for the Saemangeum network. The Maximum NFE was set to 15.5 (106), in consideration of the scales and characteristics of the problems, although Yoo et al. [48] did not present their values.
Table 9 presents the results that Yoo et al. [48] obtained by applying the HSA, and the results of the MTHSA. According to the application, the results for the MTHSA were of a higher quality for all three design plans.
That is, the MTHSA is shown to be capable of designing optimal pipe diameters that satisfy the hydraulic conditions, at costs that are lower by approximately 1.77% to 2.79%. In the case of a real water distribution network design, which requires large construction costs, a small percentage of cost reduction can reduce the construction costs by a large amount. In particular, it is observed that the MTHSA can provide economic saving effects amounting to approximately 300 million KRW in the case of the Looped Type (Plan 2), in which relatively more diverse water supply paths can appear.
Table 10 present the arithmetic statistical values obtained when 30 independent runs were implemented for each plan. It can be observed that the MTHSA achieved even values for all plans, without any large deviation.
In particular, as presented in the results of Yoo et al. [48], it can be seen that in the case of Plan 2 (Looped Type), in which more water supply paths can appear, relatively larger deviations in the minimum and maximum costs appeared compared to Plans 1 and 3. However, the value of the maximum cost of the MTHSA was shown to be smaller than or similar to that of the HSA, and this result can be interpreted as showing the optimization efficiency of the MTHSA.

5. Conclusions

In the present study, an improved version of the HSA, reinforced further in terms of diversification and intensification, was proposed. The MTHSA configures a two-storied structure, and it increases the optimization efficiency by implementing multiple searches on each floor and additional searches among the better solutions. In order to verify the proposed algorithm, it was applied to pipe diameter optimization problems for WDNs, which constituted representative discrete optimization problems. In particular, the algorithm was applied to large-scale WDNs, with at least 300 decision variables to verify in actual problems. First, through analyses of the sensitivity of the MTHSA, PAR was identified as the most sensitive parameter. However, it was observed that optimization efficiency could be increased by adjusting the ratio between HMS and SHMS (HMS/SHMS). In general, in the case of small NFEs, HMS/SHMS values closer to one showed a faster and more efficient convergence. However, when HMS was set to three times the value of SHMS, efficient optimization results could be obtained, on average. In the results, for the Balerma network optimal design problem, the developed MTHSA algorithm was found to reduce the minimum cost by about 19.85% compared to the existing HSA under the same termination conditions (45,400 NFEs). In addition, when the optimal solution search calculation was continued, it was confirmed that the MTHSA showed better results even under the condition of NFEs equivalent to 10% of the existing HSA. It could be identified that in the optimal design of pipe diameters, the MTHSA was advantageous over existing algorithms, as it does not involve a pre-processing process that requires expertise, as well as not using complicated mathematical technology. In particular, when compared to the existing literature in which only meta-heuristic algorithms were used, the MTHSA showed considerably more efficient optimization results. In the Saemangeum network, it was also confirmed that the newly developed algorithm showed a superior performance compared to the existing HSA. The results showed that design costs were reduced by 1.77 to 2.29% in all three types of Saemangeum networks. Additionally, with the MTHSA, the optimization calculation structure can be easily changed by adding only one parameter (SHMS).
However, the MTHSA developed here should be applied to diverse optimization problems (e.g., numerical benchmark problems, real-world problems), in order to further verify the quality of the algorithm. Single-objective problems were considered for the optimal cost design in this study; however, the MTHSA can also applied to WDNs with multi-objective features, which are more complex than single-objective problems. Therefore, in future studies, the MTHSA will be applied to multi-objective design problems for WDNs, and also other engineering problems with various problem characteristics. In addition, if technologies to reduce or automatically adjust the number of parameters using self-adaptation techniques considering problem characteristics are developed based on the results of the sensitivity analysis, the algorithm will become more practical.

Author Contributions

H.-M.L. and A.S. surveyed the previous studies. H.-M.L., D.-G.Y. and J.-G.J. wrote the original manuscript. H.-M.L. and Y.-H.C. conducted simulations. J.-G.J. and D.-G.Y. conceived the original idea of the proposed method. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Korea Environmental Industry & Technology Institute (KEITI) through the Water Management Program for Drought, funded by the Korea Ministry of Environment (MOE) (RS-2022-KE002032).

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Holland, J.H. Adaptation in Natural and Artificial Systems; University of Michigan Press: Ann Arbor, MI, USA, 1975. [Google Scholar]
  2. Glover, F. Heuristics for Integer Programming Using Surrogate Constraints. Decis. Sci. 1977, 8, 156–166. [Google Scholar] [CrossRef]
  3. Kirkpatrick, S.; Vecchi, M.P. Optimization by Simulated Annealing. Science 1983, 220, 671–680. [Google Scholar] [CrossRef] [PubMed]
  4. Dorigo, M. Optimization, Learning, and Natural Algorithms. Ph.D. Thesis, Politecnico di Milano, Milan, Italy, 1992. [Google Scholar]
  5. Eberhart, R.C.; Kennedy, J. A New Optimizer Using Particle Swarm Theory. In Proceedings of the Sixth International Symposium on Micro Machine and Human Science, Nagoya, Japan, 4–6 October 1995; Volume 1, pp. 39–43. [Google Scholar]
  6. Storn, R.; Price, K.V. Minimizing the Real Functions of the ICEC’96 Contest by Differential Evolution. In Proceedings of the International Conference on Evolutionary Computation, Nagoya, Japan, 20–22 May 1996; pp. 842–844. [Google Scholar]
  7. Geem, Z.W.; Kim, J.H.; Loganathan, G.V. A New Heuristic Optimization Algorithm: Harmony Search. Simulation 2001, 76, 60–68. [Google Scholar] [CrossRef]
  8. Kim, J.H.; Geem, Z.W.; Kim, E.S. Parameter Estimation of the Nonlinear Muskingum Model Using Harmony Search. J. Am. Water Resour. Assoc. 2001, 37, 1131–1138. [Google Scholar] [CrossRef]
  9. Nakrani, S.; Tovey, C. On Honey Bees and Dynamic Server Allocation in Internet Hosting Centers. Adapt. Behav. 2004, 12, 223–240. [Google Scholar] [CrossRef]
  10. Yang, X.S. Firefly Algorithm, Levy Flights and Global Optimization. In Research and Development in Intelligent Systems XXVI; Springer: London, UK, 2010; pp. 209–218. [Google Scholar]
  11. Yang, X.S.; Deb, S. Engineering Optimisation by Cuckoo Search. Int. J. Math. Model. Numer. Optim. 2010, 1, 330–343. [Google Scholar] [CrossRef]
  12. Yang, X.S. A New Metaheuristic Bat-Inspired Algorithm. In Nature Inspired Cooperative Strategies for Optimization (NICSO 2010); Springer: Berlin/Heidelberg, Germany, 2010; pp. 65–74. [Google Scholar]
  13. Eskandar, H.; Sadollah, A.; Bahreininejad, A.; Hamdi, M. Water Cycle Algorithm—A Novel Metaheuristic Optimization Method for Solving Constrained Engineering Optimization Problems. Comput. Struct. 2012, 110, 151–166. [Google Scholar] [CrossRef]
  14. Sadollah, A.; Bahreininejad, A.; Eskandar, H.; Hamdi, M. Mine Blast Algorithm for Optimization of Truss Structures with Discrete Variables. Comput. Struct. 2012, 102, 49–63. [Google Scholar] [CrossRef]
  15. Paik, K.; Kim, J.H.; Kim, H.S.; Lee, D.R. A Conceptual Rainfall-Runoff Model Considering Seasonal Variation. Hydrol. Process. 2005, 19, 3837–3850. [Google Scholar]
  16. Geem, Z.W. Improved Harmony Search from Ensemble of Music Players. In Knowledge-Based Intelligent Information and Engineering Systems; Springer: Berlin/Heidelberg, Germany, 2006; pp. 86–93. [Google Scholar]
  17. Wang, C.M.; Huang, Y.F. Self-Adaptive Harmony Search Algorithm for Optimization. Expert Syst. Appl. 2010, 37, 2826–2837. [Google Scholar] [CrossRef]
  18. Pan, Q.K.; Suganthan, P.N.; Tasgetiren, M.F.; Liang, J.J. A Self-Adaptive Global Best Harmony Search Algorithm for Continuous Optimization Problems. Appl. Math. Comput. 2010, 216, 830–848. [Google Scholar] [CrossRef]
  19. Geem, Z.W.; Sim, K.B. Parameter-Setting-Free Harmony Search Algorithm. Appl. Math. Comput. 2010, 217, 3881–3889. [Google Scholar] [CrossRef]
  20. Yildiz, A.R. Hybrid Taguchi-Harmony Search Algorithm for Solving Engineering Optimization Problems. Int. J. Ind. Eng. Theory Appl. Pract. 2008, 15, 286–293. [Google Scholar]
  21. Fesanghary, M.; Mahdavi, M.; Minary-Jolandan, M.; Alizadeh, Y. Hybridizing Harmony Search Algorithm with Sequential Quadratic Programming for Engineering Optimization Problems. Comput. Methods Appl. Mech. Eng. 2008, 197, 3080–3091. [Google Scholar] [CrossRef]
  22. Kaveh, A.; Talatahari, S. Particle Swarm Optimizer, Ant Colony Strategy and Harmony Search Scheme Hybridized for Optimization of Truss Structures. Comput. Struct. 2009, 87, 267–283. [Google Scholar] [CrossRef]
  23. Wang, L.; Pan, Q.K.; Tasgetiren, M.F. Minimizing the Total Flow Time in a Flow Shop with Blocking by Using Hybrid Harmony Search Algorithms. Expert Syst. Appl. 2010, 37, 7929–7936. [Google Scholar] [CrossRef]
  24. Wang, L.; Yang, R.; Xu, Y.; Niu, Q.; Pardalos, P.M.; Fei, M. An Improved Adaptive Binary Harmony Search Algorithm. Inf. Sci. 2013, 232, 58–87. [Google Scholar] [CrossRef]
  25. Yuan, X.; Zhao, J.; Yang, Y.; Wang, Y. Hybrid Parallel Chaos Optimization Algorithm with Harmony Search Algorithm. Appl. Soft Comput. 2014, 17, 12–22. [Google Scholar] [CrossRef]
  26. Sivasubramani, S.; Swarup, K.S. Multi-Objective Harmony Search Algorithm for Optimal Power Flow Problem. Int. J. Electr. Power Energy Syst. 2011, 33, 745–752. [Google Scholar] [CrossRef]
  27. Mahdavi, M.; Fesanghary, M.; Damangir, E. An Improved Harmony Search Algorithm for Solving Optimization Problems. Appl. Math. Comput. 2007, 188, 1567–1579. [Google Scholar] [CrossRef]
  28. dos Santos Coelho, L.; Mariani, V.C. An Improved Harmony Search Algorithm for Power Economic Load Dispatch. Energy Convers. Manag. 2009, 50, 2522–2526. [Google Scholar] [CrossRef]
  29. Jaberipour, M.; Khorram, E. Two Improved Harmony Search Algorithms for Solving Engineering Optimization Problems. Commun. Nonlinear Sci. Numer. Simul. 2010, 15, 3316–3331. [Google Scholar] [CrossRef]
  30. Gupta, C.; Jain, S. New Approach for Function Optimization: Amended Harmony Search. In Advances in Intelligent Informatics; Springer International Publishing: Cham, Switzerland, 2015; pp. 643–653. [Google Scholar]
  31. Yadav, P.; Kumar, R.; Panda, S.K.; Chang, C.S. An Intelligent Tuned Harmony Search Algorithm for Optimisation. Inf. Sci. 2012, 196, 47–72. [Google Scholar] [CrossRef]
  32. Arul, R.; Ravi, G.; Velusami, S. An Improved Harmony Search Algorithm to Solve Economic Load Dispatch Problems with Generator Constraints. Electr. Eng. 2014, 96, 55–63. [Google Scholar] [CrossRef]
  33. Ashrafi, S.M.; Dariane, A.B. Performance Evaluation of an Improved Harmony Search Algorithm for Numerical Optimization: Melody Search (MS). Eng. Appl. Artif. Intell. 2013, 26, 1301–1321. [Google Scholar] [CrossRef]
  34. Savic, D.A.; Walters, G.A. Genetic Algorithms for Least-Cost Design of Water Distribution Networks. J. Water Resour. Plan. Manag. 1997, 123, 67–77. [Google Scholar] [CrossRef]
  35. Dandy, G.C.; Simpson, A.R.; Murphy, L.J. An Improved Genetic Algorithm for Pipe Network Optimization. Water Resour. Res. 1996, 32, 449–458. [Google Scholar] [CrossRef]
  36. Reca, J.; Martinez, J.; Gil, C.; Banos, R. Application of Several Meta-Heuristic Techniques to the Optimization of Real Looped Water Distribution Networks. Water Resour. Manag. 2008, 22, 1367–1379. [Google Scholar] [CrossRef]
  37. Simpson, A.R.; Maier, H.R.; Foong, W.K.; Phang, K.Y.; Seah, H.Y.; Tan, C.L. Selection ff Parameters for Ant Colony Optimization Applied to the Optimal Design of Water Distribution Systems. In Proceedings of the International Congress on Modeling and Simulation, Canberra, Australia, 10–13 December 2001; pp. 1931–1936. [Google Scholar]
  38. Maier, H.R.; Simpson, A.R.; Zecchin, A.C.; Foong, W.K.; Phang, K.Y.; Seah, H.Y.; Tan, C.L. Ant Colony Optimization for Design of Water Distribution Systems. J. Water Resour. Plan. Manag. 2003, 129, 200–209. [Google Scholar] [CrossRef]
  39. Cunha, M.; Sousa, J. Water Distribution Network Design Optimization: Simulated Annealing Approach. J. Water Resour. Plan. Manag. 1999, 125, 215–221. [Google Scholar] [CrossRef]
  40. Eusuff, M.M.; Lansey, E.K. Optimisation of Water Distribution Network Design Using the Shuffled Frog Leaping Algorithm. J. Water Resour. Plan. Manag. 2003, 129, 210–225. [Google Scholar] [CrossRef]
  41. Geem, Z.W. Particle-Swarm Harmony Search for Water Network Design. Eng. Optim. 2009, 41, 297–311. [Google Scholar] [CrossRef]
  42. Sedki, A.; Ouazar, D. Hybrid Particle Swarm Optimization and Differential Evolution for Optimal Design of Water Distribution Systems. Adv. Eng. Inform. 2012, 26, 582–591. [Google Scholar] [CrossRef]
  43. Sadollah, A.; Yoo, D.G.; Kim, J.H. Improved Mine Blast Algorithm for Optimal Cost Design of Water Distribution Systems. Eng. Optim. 2014, 47, 1602–1618. [Google Scholar] [CrossRef]
  44. Lee, H.M.; Yoo, D.G.; Sadollah, A.; Kim, J.H. Optimal cost design of water distribution networks using a decomposition approach. Eng. Optim. 2016, 48, 2141–2156. [Google Scholar] [CrossRef]
  45. Schwartz, R.; Housh, M.; Ostfeld, A. Least-cost robust design optimization of water distribution systems under multiple loading. J. Water Resour. Plan. Manag. 2016, 142, 04016031. [Google Scholar] [CrossRef]
  46. El-Ghandour, H.A.; Elbeltagi, E. Comparison of five evolutionary algorithms for optimization of water distribution networks. J. Comput. Civ. Eng. 2017, 32, 04017066. [Google Scholar] [CrossRef]
  47. Lee, E.H.; Lee, H.M.; Yoo, D.G.; Kim, J.H. Application of a meta-heuristic optimization algorithm motivated by a vision correction procedure for civil engineering problems. KSCE J. Civ. Eng. 2018, 22, 2623–2636. [Google Scholar] [CrossRef]
  48. Yoo, D.G.; Lee, H.M.; Sadollah, A.; Kim, J.H. Optimal Pipe Size Design for Looped Irrigation Water Supply System Using Harmony Search: Saemangeum Project Area. Sci. World J. 2015, 2015, 651763. [Google Scholar] [CrossRef] [PubMed]
  49. Rossman. EPANET 2.0 User’s Manual; EPA: Washington, DC, USA, 2000.
  50. Reca, J.; Martínez, J. Genetic algorithms for the design of looped irrigation water distribution networks. Water Resour. Res. 2006, 42, W05416. [Google Scholar] [CrossRef]
  51. K-Water. Water Facilities Construction Cost Estimation Report; K-water: Daejeon, Republic of Korea, 2010. [Google Scholar]
  52. Bolognesi, A.; Bragalli, C.; Marchi, A.; Artina, S. Genetic Heritage Evolution by Stochastic Transmission in the Optimal Design of Water Distribution Networks. Adv. Eng. Softw. 2010, 41, 792–801. [Google Scholar] [CrossRef]
  53. Zheng, F.; Simpson, A.R.; Zecchin, A.C. A Decomposition and Multistage Optimization Approach Applied to the Optimization of Water Distribution Systems with Multiple Supply Sources. Water Resour. Res. 2013, 49, 380–399. [Google Scholar] [CrossRef]
  54. Tolson, B.A.; Asadzadeh, M.; Maier, H.R.; Zecchin, A. Hybrid Discrete Dynamically Dimensioned Search (HD-DDS) Algorithm for Water Distribution System Design Optimization. Water Resour. Res. 2009, 45, W12416. [Google Scholar] [CrossRef]
Figure 1. Layout of the Balerma network.
Figure 1. Layout of the Balerma network.
Sustainability 16 06689 g001
Figure 2. Layout of the Saemangeum networks.
Figure 2. Layout of the Saemangeum networks.
Sustainability 16 06689 g002
Figure 3. Schematic structure and components of the MTHSA.
Figure 3. Schematic structure and components of the MTHSA.
Sustainability 16 06689 g003
Figure 4. Nodal pressure distribution of optimal results in MTHSA (I).
Figure 4. Nodal pressure distribution of optimal results in MTHSA (I).
Sustainability 16 06689 g004
Table 1. Candidate pipe diameters and costs for the Balerma network.
Table 1. Candidate pipe diameters and costs for the Balerma network.
Pipe No.Pipe Diameter (mm)Cost (EUR/m)
11137.22
2126.69.1
3144.611.92
4162.814.84
5180.818.38
6226.228.6
728545.39
8361.876.32
9452.2124.64
10581.8215.85
Table 2. Commercial pipe diameters and their costs per length for the Saemangeum networks.
Table 2. Commercial pipe diameters and their costs per length for the Saemangeum networks.
Pipe No.Pipe Diameter (mm)Cost (KRW/m)
Construction CostsMaterial
Costs
Maintenance
Costs
18065,00015,0006500
210065,99927,5836600
315076,41040,6867641
420086,02858,7168603
525096,13581,1609614
6300105,325103,23110,533
7350113,818125,10711,382
8400126,797148,83612,680
9450136,250155,52213,625
10500147,792181,82314,779
11600171,991211,39617,199
12700211,413273,52821,141
13800307,640339,74030,764
14900359,048384,61935,905
151000415,702451,93241,570
161100482,074547,22448,207
171200576,736606,96257,674
181350687,390716,07568,739
Table 3. Effects of the MTHSA parameter (HMS) for the Balerma network.
Table 3. Effects of the MTHSA parameter (HMS) for the Balerma network.
ResultsSHMS = 20
Maximum NFEs = 45,400
HMCR = 0.90
PAR = 0.02
HMS
20 (MTHSA (II))60 (MTHSA (I))80140180
Best cost2,084,8442,148,7402,173,7992,257,7712,258,333
Average cost2,172,2832,241,8562,296,8992,398,9902,357,524
Worst cost2,275,7102,337,8672,387,6762,511,5072,419,805
Table 4. Effects of the MTHSA parameter (SHMS) for the Balerma network.
Table 4. Effects of the MTHSA parameter (SHMS) for the Balerma network.
ResultsHMS = 60
Maximum NFEs = 45,400
HMCR = 0.90
PAR = 0.02
SHMS
1520 (MTHSA (I))60
Best cost2,189,6822,148,7402,092,939
Average cost2,256,9332,241,8562,182,203
Worst cost2,339,4952,337,8672,235,220
Table 5. Effects of the HMCR parameter used in the MTHSA for the Balerma network.
Table 5. Effects of the HMCR parameter used in the MTHSA for the Balerma network.
ResultsHMS = 60
SHMS = 20
Maximum NFEs = 45,400
PAR = 0.02
HMCR
0.600.700.800.90 (MTHSA (I))0.95
Best Cost2,254,1432,249,9412,152,4852,148,7402,149,149
Average Cost2,364,2832,310,5432,259,3542,241,8562,236,766
Worst Cost2,497,4472,421,5212,377,7912,337,8672,309,875
Table 6. Effects of the PAR user parameter used in the MTHSA for the Balerma network.
Table 6. Effects of the PAR user parameter used in the MTHSA for the Balerma network.
ResultsHMS = 60
SHMS = 20
Maximum NFEs = 45,400
HMCR = 0.90
PAR
0.02 (MTHSA (I))0.050.100.200.30
Best Cost2,148,7402,328,4602,577,7562,932,6693,305,492
Average Cost2,241,8562,415,8942,661,4363,097,6953,505,115
Worst Cost2,337,8672,516,3482,732,5183,229,6483,785,797
Table 7. Initial parameters used for the Balerma network.
Table 7. Initial parameters used for the Balerma network.
Control ParameterCase
MTHSA (I)MTHSA (II)MTHSA (III)
HMS602060
SHMS202020
HMCR0.900.900.90
PAR0.020.020.02
NFEs45,40045,4001,000,000
Table 8. Numerical optimization results for the Balerma network.
Table 8. Numerical optimization results for the Balerma network.
MethodBest CostAverage CostWorst CostMaximum NFEs
Decompose-DE1.9231.9311.935N/A
GA3.738N/A 1N/A45,400
SA3.475N/AN/A45,400
MSATS3.298N/AN/A45,400
HAS (I)2.601N/AN/A45,400
PSHS2.633N/AN/A45,400
GHEST2.178N/AN/A45,400
MBA2.2112.3842.43345,400
IMBA (I)2.0642.2022.33845,400
MTHSA (I)2.1492.2412.33445,400
MTHSA (II)2.0852.1722.27545,400
IMBA (II)2.013N/AN/A250,000
MTHSA (III)1.9322.0292.0781,000,000
NLP-DE11.9561.9571.9591,000,000
NLP-DE21.9231.9271.9342,000,000
HAS (II)2.018N/AN/A10,000,000
DE31.9821.9861.98910,000,000
GANOME2.3022.3342.35010,000,000
HD-DDS-21.956N/AN/A10,000,000
HD-DDS-11.941N/AN/A30,000,000
1 Not available.
Table 9. Optimal design results and cost comparison with the previous study for the Saemangeum networks.
Table 9. Optimal design results and cost comparison with the previous study for the Saemangeum networks.
Target NetworkOptimal Design Costs (KRW)
HSA
(Yoo et al. [48])
MTHSAImproved (%)
Branched Type (Plan 1)10,044,962,4059,823,580,100(−) 2.20
Looped Type (Plan 2)10,182,733,2959,898,369,770(−) 2.79
Pumped Type (Plan 3)11,586,379,38011,380,850,160(−) 1.77
Table 10. Numerical optimization results for the Saemangeum networks.
Table 10. Numerical optimization results for the Saemangeum networks.
Target NetworkMethodBest CostAverage CostWorst CostMaximum NFEs
Branched Type (Plan 1)HSA10,044,962,405N/AN/AN/A
MTHSA9,823,580,1009,867,502,03210,030,934,66515,500,000
Looped Type (Plan 2)HSA10,182,733,295N/AN/AN/A
MTHSA9,898,369,77010,080,620,69110,264,994,76015,500,000
Pumped Type (Plan 3)HSA11,586,379,380N/AN/AN/A
MTHSA11,380,850,16011,443,353,94411,559,108,07015,500,000
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Lee, H.-M.; Sadollah, A.; Choi, Y.-H.; Joo, J.-G.; Yoo, D.-G. Development and Sensitivity Analysis of an Improved Harmony Search Algorithm with a Multiple Memory Structure for Large-Scale Optimization Problems in Water Distribution Networks. Sustainability 2024, 16, 6689. https://doi.org/10.3390/su16156689

AMA Style

Lee H-M, Sadollah A, Choi Y-H, Joo J-G, Yoo D-G. Development and Sensitivity Analysis of an Improved Harmony Search Algorithm with a Multiple Memory Structure for Large-Scale Optimization Problems in Water Distribution Networks. Sustainability. 2024; 16(15):6689. https://doi.org/10.3390/su16156689

Chicago/Turabian Style

Lee, Ho-Min, Ali Sadollah, Young-Hwan Choi, Jin-Gul Joo, and Do-Guen Yoo. 2024. "Development and Sensitivity Analysis of an Improved Harmony Search Algorithm with a Multiple Memory Structure for Large-Scale Optimization Problems in Water Distribution Networks" Sustainability 16, no. 15: 6689. https://doi.org/10.3390/su16156689

APA Style

Lee, H. -M., Sadollah, A., Choi, Y. -H., Joo, J. -G., & Yoo, D. -G. (2024). Development and Sensitivity Analysis of an Improved Harmony Search Algorithm with a Multiple Memory Structure for Large-Scale Optimization Problems in Water Distribution Networks. Sustainability, 16(15), 6689. https://doi.org/10.3390/su16156689

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