Optimal Reactive Power Generation for Radial Distribution Systems Using a Highly Effective Proposed Algorithm

Faculty of Electrical and Electronics Engineering, Ho Chi Minh City University of Technology and Education, Ho Chi Minh City 700000, Vietnam Faculty of Electrical Engineering Technology, Industrial University of Ho Chi Minh City, Ho Chi Minh City, Vietnam Power System Optimization Research Group, Faculty of Electrical and Electronics Engineering, Ton Duc ,ang University, Ho Chi Minh City 700000, Vietnam


Introduction
Electric distribution networks have a very important role in receiving electricity from transmission power network and supplying the electricity to loads. e main difference between the distribution networks and transmission networks is voltage level, leading to another difference, which is total active power loss due to the impact of resistance of conductors. e active power loss is dependent on the result of RI 2 [1] (where R is the resistance of conductor and I is current flowing the conductor). Current value is a main factor to result in a high active power loss in distribution networks while R is a constant in the networks. e smaller the voltage is, the higher the active power loss is. Hence, active power loss is a significant issue when distribution networks are working for supplying power energy to loads.
In order to reduce the high active power loss in distribution networks, experts have proposed two basic methods including network reconfiguration [2,3] and shunt capacitor installation [4,5]. e network reconfiguration method is to change status of switches, either open or close to change direction of current. Basically, distribution networks are supplied at one point, which is called slack node, and it is obvious that all loads in the networks are being supplied by the slack node via distribution lines. us, the method cannot reduce power supplied by the slack node and just mainly reduces high current in lines with high resistance or long length and increases lower current in other lines. By using the method, power loss is effectively reduced. In addition, voltage profile is also improved, but the improvement is not certain or insignificant. On the contrary, the second method using shunt capacitors can reduce reactive power that is supplied by the distribution lines. Loads can consume reactive power from both the added capacitors and the distribution lines or only consume reactive power from the added capacitors. As a result, current in distribution lines can be reduced considerably and power factor is increased effectively. In addition, another benefit from the capacitor installation is the reduction of voltage drop in the line. In fact, as current is smaller, the voltage drop is decreased accordingly. In addition to the two basic methods, other methods can be applied such as (1) placement of distributed generators [6,7], (2) the combination of reconfiguration and capacitor placement [8,9], (3) the combination of reconfiguration and distributed generator placement [10,11], and (4) the combination of capacitor placement and distributed generator placement [12,13]. In this paper, we focus on the second basic method of optimally installing capacitors with the task of determining the best location and the best rated reactive power. e best location and the best rated power of these added capacitors are for reaching two single-objective functions in which the first objective function is to minimize the total active power loss on all distribution lines  and the second objective function is to minimize the total cost of energy loss and capacitor investment [34,[46][47][48][49]. In addition, operation limit of conductor and operating voltage of loads are always supervised seriously via the consideration of maximum current of lines [50] and the consideration of upper and lower voltage [51]. e problem of capacitor placement has attracted a huge number of researchers in proposing optimization tools and capacitor placement strategies based on configuration and practical analysis. Approximately all the applied methods are different; however, the common study of all the methods is the active power loss reduction. In [1], a proposed method with two stages was applied for two systems with 15 and 33 nodes. In the first stage, a sequence of compensated nodes is first determined by using an iterative algorithm with the placement of one capacitor for minimizing power loss. en, the optimal size of capacitors at the determined nodes was found by minimizing a loss saving equation, which is a function of capacitors' current. is method could reach lower power loss than original networks without capacitor placement. However, the method had to suffer from the limits of application for large-scale problem with a high number of load nodes because each capacitor is tried to be placed at all nodes excluding slack node in the first stage. So, it will be time consuming for trying one by one node in a large-scale system with too many nodes. For example, it must try fourteen times for 15-node network and 32 times for 33-node network. us, the higher the number of nodes is, the longer the simulation time is. e method can solve the high power loss issue, but it is not a good choice for the radial distribution networks because there was no comparison between the method and other ones in the study. Another similar method was proposed in [4] for maximizing saving energy loss as compared to original radial distribution networks. e study replicated the first stage of determining compensated nodes where reactive power is necessary to reduce current flowing in distribution lines. en, the second stage is to determine the most appropriate size for each shunt capacitor by maximizing the saving power loss compared to original network. e method can solve the problem easily and successfully, but its applications for large-scale problem also suffer from the same restriction as the two-stage method [1] because of the first stage. In fact, the method was only applied for 15-node and 33-node networks. e method was only superior to the two-stage method [1]. In 2013, another two-stage method (TSM) [14] was applied for the same problem but the application was wider and more successful thanks to the modifications on the first stage. e first stage for finding the most suitable locations is performed by using cross check fuzzy expert system and loss reduction index. So, the two-stage method could avoid the significant restriction of the methods [1,4]. e large-scale problem with 69 nodes was successfully solved and the method could reach better loss reduction than other previous methods; however, the simulation time was still the major disadvantage of the method. In [15], the two-stage method proposed in [4] was applied to determine distributed generator location and size in the radial distribution network. e method could find location and size of the distributed generator successfully and effectively as the obtained power loss was less than that in capacitor location and size determination problem. However, the method one more time shows its disadvantage since the most complicated study case was 33-node network. Clearly, the two-stage methods could not reach the highest performance for the problem of determining location and size of capacitors and distributed generators. Due to major disadvantages such as not applicable for large-scale network and time consuming, the two-stage methods were not applied widely and they must be replaced with more potential metaheuristics such as genetic algorithm variants, particle swarm optimization (PSO) variants, and other recent ones. PSO based on inertia weight and constriction factor (IWC-PSO) was applied for finding reactive power generation of capacitors while the loss sensitivity factor method was proposed to determine candidate nodes, where capacitor placement is necessary [16]. e loss sensitivity factor method was used to identify lowvoltage nodes or capacitor location that can improve voltage of other low-voltage nodes, where capacitors are not installed. Single and multiple capacitors were installed in five distribution networks with 10, 15, 34, 69, and 85 nodes, and voltage profile was significantly improved as compared to voltage profile in original networks and results from [1]. It should be noted that objective function of the study [1] was loss reduction, whereas that in [16] was voltage enhancement. So, the comparison between the two-stage method [1] and IWC-PSO [16] was not suitable. e IWC-PSO continued to be applied for the problem with two single-objective functions including power loss and voltage profile [17]. e fuzzy method was used to identify candidate nodes, and then the PSO method determined the most suitable size 2 Complexity for capacitors at the nodes. e largest study case was the 69node network, and results were compared to original networks. Another study [18] also applied the fuzzy method to find the most suitable locations to place capacitors, and then multiagent particle swarm optimization (MAPSO) was proposed to determine the size of capacitors. MAPSO was demonstrated to be superior to only the conventional PSO for the 69-node radial distribution network with the consideration of active and reactive power losses and voltage profile. A set of different PSO methods with different distributions (including Gaussian, Cauchy, and chaotic distributions) and different equations for calculating velocity (including weight inertia factor and constriction factor) was applied for identifying location and size of capacitors [19]. e study is different from other applications of PSO methods above since capacitors' location was selected to be control variables. Due to the selection of control variables, the study [19] could skip the first stage of determining location of capacitors by using loss reduction index as the studies [1,4], loss sensitivity factor as the study [16], and fuzzy method as the study [17,18]. ere were fifteen PSO methods to be applied for determining the best one for the problem of finding both location and size of capacitors. e comparisons from two study cases in the 9-feeder radial distribution network showed that the PSO method with uniform distribution and chaotic distribution was the best one for the smallest power loss. In addition to the comparison among the PSO methods, the best PSO method was compared to genetic algorithm (GA) and tabu search algorithm (TSA). In general, the PSO method was the best one among fifteen PSO methods and superior to two other lowly effective methods such as GA and TSA for only a small-scale system with 9 feeders and 10 nodes. Hence, the real performance of this method was still a question for the problem. Different GA variants including conventional GA [20][21][22][23], micro GA (MGA) [24], real coded genetic algorithm (RCGA) [25], and the combination of fuzzy and GA (FGA) [26] were the applied solution methods to optimally place capacitors in the radial distribution networks. e applications of conventional GA did not demonstrate the high performance of GA because the study cases were simple and comparisons were mainly between the original networks and networks with capacitor placement. In fact, Taiwan network and Iran network were, respectively, studied in [20,23] while 23-node network and 33-node network were, respectively, studied in [21,22]. ese studies were poor in comparisons and study cases. In [24], MGA was applied for Italian network and compared with GA for comparison. In [25], capacitors were installed in three networks with 15, 34, and 69 nodes by using RCGA. e power loss reduction of the cases with and without capacitor placement was compared. Clearly, all the studies have the same shortcoming of poor study cases and comparisons. In addition to the application of GA for single-objective problem with only power loss reduction, a multiple-objective problem with voltage profile improvement and total cost reduction was solved by the implementation of GA for getting a set of solutions. en, the fuzzy method was employed to determine the most appropriate compromise solution. e paper only executed the comparison of networks with and without capacitor placement rather than showing the real performance of GA as compared to other methods. So, GA was not a real effective method for the problem [27].
In addition to these method groups, other smaller groups were also applied for the same problem of capacitor placement such as mixed integer nonlinear programming-based method (MINPM) [27], gravitational search algorithm (GSA) [28], the combination of GSA and weight inertia factor-based PSO (WIFPSO-GSA) [29], bacterial foraging optimization algorithm (BFOA) [30,31], flower pollination algorithm (FPA) [32,33], teaching-learning algorithm (TLA) [34], whale optimization algorithm (WOA) [35], power loss index-based improved harmony algorithm (PLI-IHA) [36], cuckoo search algorithm (CSA) [37], improved mutation technique-based differential evolution (IMT-DE) [38], moth swarm algorithm (MSA) [39], ant colony algorithm based on loss sensitivity factor (LSF-ACA) [40], heuristic method based on network configuration (NCB-HM) [41,42], combined practical method (CPM) [43,44], hybrid method (HM) [45], direct search optimization algorithm (DSOA) [46], penalty free method-based heuristic algorithm (PFHA) [47], inclusion and variable interchange algorithm (IVIA) [48], water cycle algorithm [49], and grey wolf algorithm (GWA) [49]. Among the methods, MINPM, NCB-HM, CPM, IVIA, and DSOA are not metaheuristic algorithms based on population and they are mainly dependent on the real configuration of networks. So, the application of the methods is not performed for arbitrary systems without the analysis on the power loss and voltage drop. Other metaheuristic algorithms can reach better results than PSO and GA method groups; however, the real performance of these methods was not demonstrated clearly. In fact, these methods have been run by setting different values to population and iterations without comparisons. It is noted that metaheuristic algorithms can result in good solutions if they spend high computation time due to high value of population and iterations. In terms of considered objective functions, almost all previous studies focused on the purpose of reducing power loss of all branches and neglecting the total compensation capacity. Observing the results from BFOA [30,31], FPA [33], NCB-HM [42], CPM [43,44], and HM [45], it could be seen that only the obtained power loss was compared, methods with smaller power loss were concluded to be more effective, and total compensation of all capacitors was not discussed. For some cases, methods with higher capacity could reach less power loss, but for other cases, lowperformance methods even with higher compensation capacity still obtained higher power loss. e shortcoming has been pointed out, and it was noted that higher compensation capacity will use higher capacitor investment purchase cost [46]. For tackling the issue, power loss and compensation capacity were converted into cost by calculating energy loss cost and considering capacitor purchase cost. e sum of energy loss cost and capacitor purchase cost was then considered as an objective for performance comparison.
In summary, the previous studies have two main shortcomings in which the former is not to further investigate the convergence speed of compared methods and the latter is not to consider compensation capacity. In this paper, Complexity the two shortcomings are solved by collecting population, iteration, and the simulation time in addition to the consideration of sum of energy loss cost and capacitor purchase cost as the second single-objective function. A proposed method, called modified stochastic fractal search optimization algorithm (MSFS), is applied to reach two single-objective functions including total power loss and the total cost of loss and capacitor of distribution systems. e performance of the proposed method is investigated by comparing total loss, total cost, and convergence speed. Furthermore, the proposed method is also demonstrated to be more effective than its original method, stochastic fractal search optimization algorithm (SFS). SFS was also an effective optimization tool developed in 2015 for reaching results from a 23function set [52]. e method could attract the concern from a huge number of researchers thanks to its superiority to eight popular and effective algorithms such as PSO, differential evolution algorithm (DE), gravitational search optimization algorithm (GSOA), artificial bee colony optimization algorithm (ABCOA), cuckoo search optimization algorithm (CSOA), improved version of CSOA (ICSOA), animal migration optimization algorithm (AMOA), and backtracking search optimization algorithm (BSOA). SFS was applied and reached promising results for different optimization problems in electrical engineering such as distributed generator placement in distribution systems [53], environmental and economic dispatch [54], PID controller design [55], economic load dispatch [56], and reconfiguration of distribution networks [11]. SFS is characterized by three mechanisms of generating new candidate solutions including diffusion mechanism and two update mechanisms. e diffusion mechanism uses Gaussian distribution mean while the two update mechanisms use mutation of DE with modifications. e original structure of SFS was considered to be ineffective for complicated problems, and a number of improved versions were suggested for different problems such as multiple area economic load dispatch [57], optimal dispatch of reactive power [58], photovoltaic system parameter estimation [59], optimal transmission network power flow [60], and optimal generation control of interconnected power systems [61]. In the proposed MSFS method, we use modifications on the three mechanism to improve the search performance. In the diffusion mechanism, we employ another new formula and keep one out of two old formulas of SFS to balance the global search and local search. In the two update mechanisms, a new algorithm is proposed for the two update mechanisms and only performed for the first half of the whole population with better quality. e new algorithm can balance the global and local search effectively while the number of fitness evaluations is decreased to half of SFS. Each proposed modification and all the proposed modifications are tested on three radial distribution networks (RDNs) including 33, 69, and 85 nodes with two different single-objective functions such as reduction of power loss and reduction of the sum of energy loss cost and capacitor investment cost. As a result, the core contributions of paper are summarized as follows: (i) Core shortcomings of SFS are found. (ii) New formula and new algorithms are proposed for tackling shortcomings of SFS. (iii) Each proposed modification and the whole modifications are tested and proved (iv) e proposed MSFS is superior to its conventional method in terms of using smaller number of computation iterations, spending shorter simulation time, and reaching better results.
In terms of the organization of the paper, the remaining parts of the paper are as follows. Two single-objective functions and considered constraints are presented in detail in Section 2. e original SFS method and the proposed MSFS method are clearly explained in Section 3. Section 4 explains the implementation of MSFS for the problem. In Section 5, three systems with two different objective functions are employed to run SFS and MSFS for comparison and discussion. Finally, the summary and conclusions are presented in Section 6.

Objective Functions.
e conventional problem of adding capacitors in the distribution system aimed to reduce the active power loss [14] as much as possible while the total compensation capacity was not considered as long as the total compensation capacity was not higher than reactive power demand of loads. However, some studies have pointed out the major issue of the problem was the cost of buying capacitors and energy loss. So, energy loss cost and capacitor investment cost were considered as an objective. In this paper, we consider two single-objective functions in which the first objective is to minimize the total active power loss (TPL) while the second objective is to minimize the sum of energy loss cost and capacitor investment cost. e two objectives are presented as follows.

TPL Reduction.
Distribution network is comprised of a high number of distribution lines, and total active power loss on the lines is significant as a result. Consequently, capacitors are added at suitable nodes to reduce current of lines and reduce the total active power loss on all lines. e first objective of the problem is to minimize the total active power loss as shown in the following formula [30]: where I i is the current magnitude in the ith branch and R i is the ith branch resistance.

Energy Loss and Investment Cost Reduction.
When placing capacitors in distribution networks, power loss and energy loss can be reduced or even highly reduced if the compensation capacity is very high. e high compensation capacity leads to a high cost for capacitor investment, and the energy loss reduction cannot compensate the high cost of capacitor investment. As the sum of energy loss cost and the investment cost is low, the capacitor placement is effective. us, the energy loss cost and the capacitor investment cost are the main factors of the second objective, and the objective is mathematically formed as follows [46]: where TC is the total cost of energy loss cost and investment cost; TPL m is the total active power loss of the mth load level; Pr e and Pr c are the price of each kWh and each kVAr; T m is the number of hours for the mth load level; and M is the number of load levels.

Constraints.
e distribution power network can work stably and effectively if voltage magnitude of all nodes is within working range and current of branches is not higher than the maximum limit of conductor. In addition, the total reactive power of all capacitors must not be higher than the total reactive power of all loads. e whole constraints are as follows: (i) Voltage limit constraint: voltage of each node in the system must be within the working range of the network. e working range is the same for all nodes while the voltage magnitude of nodes can be different. e constraint is mathematically formulated as follows [34]: where U min and U max are the minimum and maximum working voltage of all nodes in the radial distribution networks and U i is the voltage of the ith node. e constraint plays a very important role in stabilizing the operation of loads and the whole distribution system as well as achieving the high benefit for both customers and power companies. In case that voltage of loads is smaller than U min , load cannot work at the rated power and loads can be damaged if the voltage of loads is higher than U max . Limit constraint of total capacity of capacitors: capacitors are very useful in improving voltage profile and reducing power loss of the distribution networks. However, as the size of capacitors is high enough, the power loss becomes higher, whereas voltage profile is always improved and even voltage of loads can be equal to rated value or higher than rated value. Consequently, the total compensated reactive power of capacitors must be limited as the following constraints [30]: where Q cj is the reactive power generation of the jth capacitor and Q c,max is the maximum reactive power generated by all capacitors. As demonstrated in [62], the maximum compensated power of capacitors could not be higher than the total reactive power of all loads because the redundant reactive power will continue to flow into other nodes where reactive power is necessary. e active power flows will cause power loss unintentionally. So, the compensation capacity must satisfy the constraint below [43]: where Q i is the ith load reactive power.
Limit constraint of branch current: one main parameter of each distribution line is the maximum current magnitude that conductor can work stably in operating time. So, active and reactive power flows that are distributed in all branches should be appropriate and not higher than the capacity of distribution lines. In order to satisfy the requirement, active and reactive power flows are converted into current magnitude and the current magnitude must not be higher than the maximum current of conductor as shown in the model below: where I max i is the ith branch capacity.

The Proposed Method
In the equation, normrnd (Gbest, sigma) produces random numbers around Gbest with a standard deviation of sigma, which is obtained by

e First Update Technique.
In the first update technique, SFS can update the whole population at most, but normally the number of newly updated solutions is less than the population size. In the first step, the whole population must be sorted based on the fitness function. Better solutions with smaller fitness function are placed at the end of the population, whereas worse solutions with higher fitness function are put at the top of the population. e order number of each solution after sorting is stored and set to R s . Each solution is evaluated by calculating a quality index as the following equation: e quality index of each solution cannot be lower than 0 and higher than 1. e best solution with the last position will get the highest index with QI s � 1, and the worst solution with the first position will get the lowest index with QI s � 1/ N Po .
In the second step, each solution will produce a random number rd s within zero and 1. If rd s is equal to or higher than QI s , the considered solution will be updated by the following formula: For another case (i.e., rd s is less than QI s ), the sth solution remains unchanged. Clearly, solutions with high quality do not have a high possibility to be newly updated.

e Second Update Technique.
In the second update technique, all solutions have the same update possibility, but the methods for updating the solution can be different. Each solution also produces a random number rd s , and this value is compared to 0.5 for determining which method is employed for the considered solution. e second update technique is formulated as follows:

Discussion on Diffusion Technique and New Diffusion
Technique. e MSFS is proposed to tackle shortcomings existing in its original version. In the diffusion technique, two models in equation (7) can be used, but the new update using the second model, normrnd (X s , sigma), cannot be effective because the result from normrnd (X s , sigma) is a random method without a certain base while the first model, [normrnd (Gbest, sigma) + ε (Gbest − X s )], is much more effective by using the best solution and a step size based on a certain condition using [ε × (Gbest − X s )]. Clearly, the first model uses local search by searching around the so-far best solution Gbest with a certain step size, whereas the second model uses global search by searching around all current solutions with randomization. So, the first model is suitable for problem with many near global solutions while the second model may not be effective because it is a random search without certain condition.
In the new diffusion technique, shortcoming of equation (7) can be avoided and improved by using the following equation.
In equation (12), the above way can search around the best solution with large distance where the below way can exploit the search around the sth considered solution with a smaller distance. e proposed equation can support to exploit global search (i.e., nearby each considered solution) and local search (i.e., nearby the best solution) effectively. In addition, the two ways are based on certain determinations using a small step size with two solutions and a larger step size with four solutions, whereas the below way of SFS in equation (7) uses a random step size without certain condition.

Discussion on the Update Techniques.
e first update technique shown in Section 3.1.2 has two key shortcomings including the condition for update and the method for update. e condition for update is the comparison between random number rd s and the quality index QI s . Only solutions with rd s > QI s are updated while high quality solutions with high value of QI s do not have many chances to be updated. So, the condition limits the search around promising spaces but highly allows to exploit search space where low-quality solutions are existing. e second 6 Complexity shortcoming of the technique is equation (10), which uses the subtraction, i.e., "−" from the current solution instead of addition, i.e., "+" like other methods such as PSO, bat algorithm, and cuckoo search algorithm.
In the second update technique, observing two ways shown in equation (11), we can see the major disadvantage of SFS in producing new solutions by using the second update technique. e above way has the same feature as mutation technique of differential evolution algorithm (DE), whereas the below way has two differences from the mutation technique of DE. e first difference is to use "−" to subtract a step size from the current solutions to produce new solutions, whereas the second difference is to use Gbest instead of a random solution to calculate a jumping step. As indicated in [63] for improving the performance of spider optimization algorithm and indicated in [58,60] for improving the performance of SFS, the use of subtraction to decrease a step size from the current solutions could lead to worse results and slow the convergence to the global optimum solutions. e main disadvantage of SFS was shown clearly in [58] for reactive power dispatch problem and in [60] for optimal power flow problem. e two problems are complicated in terms of a high number of decision variables and complex constraints. Clearly, the two update techniques need to be improved for getting more effective result from complicated problems.

Modification on the Two Update Techniques.
For reaching good solutions in the two update mechanisms, we propose a new technique for producing new solutions. In the first update, only a half of population with better quality is newly updated. en, the new solutions are evaluated and compared to old solution to keep better one and abandon worse ones. e whole population is sorted to arrange good solutions with lower fitness at the top of the population, and worse ones are arranged at the end of the population. As a result, the first half of population is newly updated based on the following steps: Step 1: sort all solutions in ascending order of fitness function (Fit s ).
Step 2: calculate average fitness function of the population, Fit mean .
Step 3: update new solutions for the first half of population based on the following equation.
Formula (13) shows two different ways for updating new solutions X new s based on comparison result between the considered old solution and the mean solution of the whole population. In the first method, old solutions with better quality than the mean solution of the whole population (i.e., Fit s ≤ Fit mean ) will be searched around with a distance [ε ( On the contrary, old solutions with worse quality than the mean solution (i.e., Fit s > Fit mean ) are abandoned and the current best solution Gbest will be searched around with a step size of [ε (X 1 − X 2 + X 3 − X 4 )] in the second method. e two methods can exploit promising search zones effectively and reduce search time around ineffective zones. Furthermore, the application of the larger step with four solutions X 1 , X 2 , X 3 and X 4 can form a large enough jumping step to escape local zones and move fast to other global zones. e effectiveness of the proposed update techniques will be demonstrated in numerical results.

The Implementation of the Proposed
Method for the Considered Problem 4.1. Initialization. In the first step of implementation, selection of control variables should be accomplished and then boundaries of the variables should be determined. In the problem, shunt capacitors are placed in order to reduce total power loss. e total loss can be reduced in maximum if the location and the capacity are the most appropriate. In a large distribution power network with a high number of nodes, the determination of the most appropriate nodes for placing capacitors is not easy if optimization algorithms are not applied. Furthermore, it is supposed to determine exact nodes to place capacitors, and the determination of capacity is also a difficulty. So, both location and capacity are control variables and must be determined by using the proposed method in the study. e location of capacitors is from node 2 to node Nb while the total capacity of all capacitors should not be higher than the total reactive power of loads as shown in equation (5) and demonstrated in Section 2.1. Each solution is represented as follows: where L j,s and Q cj,s are the location and size of the jth capacitor corresponding to the sth solution.
Each solution X s is randomized satisfying the constraint below: where X min and X max are called minimum solution and maximum solution and are mathematically expressed as follows: where L min j,s and L max j,s are the minimum location and the maximum location in the considered distribution network and Q min cj,s and Q max cj,s are the minimum and maximum reactive power generation of the jth capacitor.

Determination of Dependent
In the equation above, control variables that are smaller than minimum values are set to the minimum values while other ones with higher value than the maximum values are set to the maximum values. For remaining control variables within the minimum and maximum values, the correction will not be executed. Because there is no violation of control variables, control variables are always accepted. However, valid control variables are not called optimal solutions. Optimal solutions must have valid control variables, valid dependent variables, and suitable objective function value.

Fitness Function.
Fitness function must be calculated to reflect the quality of solutions. In the problem, fitness function is the sum of objective function and penalty terms, in which objective function (obj) can be either TPL in equation (1) or TC in equation (2), whereas the penalty terms are always the penalty of the violation of the branch current and penalty of the violation of node voltage [62].
In the fitness function, ΔI i,s and ΔU i,s are the violation of branch current in the ith branch and the violation of voltage at the ith node corresponding to the sth solution. e two factors are calculated as follows: Equations (19) and (20) indicate that penalty terms can be equal to 0 if current is not higher than the maximum value of conductor and voltage of nodes is within the operating range of loads. In this case, the fitness function and the objective function are the same and the considered solution can be called valid solution. Nevertheless, the selection of an optimal solution needs a predetermined solution set and the best solution is then decided by using the lowest fitness value of all solutions. In the paper, fifty valid solutions are found for each study case and then the best solution is taken from the fifty available solutions.

e Procedure of Implementing MSFS for the Problem.
e whole procedure of applying the proposed MSFS for finding the best site and size of capacitors can be plotted in Figure 1 and described in detail as follows: Step 1: set values to N Po , N Iter , N Df , and DF.
Step 3: solve power flow by running the BFS technique to obtain I i,s and U i,s .
Step 6: perform the diffusion technique using equation (12) to produce new solutions.
Step 7: correct violation of new solutions by using equation (17).
Step 8: solve power flow by running the BFS technique to obtain I i,s and U i,s .
Step 10: compare new and old solutions to keep higher quality ones.
Step 11: sort the population in ascending order of fitness value.
Step 12: calculate average fitness function of the population Fit mean .
Step 13: update new solutions for the first half of population using equation (13).
Step 14: correct violation of new solutions by using equation (17).
Step 15: solve power flow by running the BFS technique to obtain I i,s and U i,s .
Step 17: at each position in the population, compare solutions to keep better one.
Step 18: sort all solutions in ascending order of fitness function.
Step 19: calculate average fitness function of the population Fit mean .
Step 20: update new solutions for the first half of population using equation (13). 8 Complexity Step 21: correct violation of new solutions by using equation (17).
Step 22: solve power flow by running the BFS technique to obtain I i,s and U i,s .
Step 24: at each position in the population, compare solutions to keep better one.
Step 26: if C Iter � N Iter , stop the iterative algorithm.
Otherwise, set C Iter � C Iter + 1 and go back to step 6.

Numerical Results
In this section, the performance of the proposed MSFS is investigated by implementing the capacitor bank placement in three distribution systems with the increase of bus number, For demonstration of the high performance of the proposed MSFS method, we focus on the comparison criteria as follows: (1) e best solution of 50 trial runs: the four implemented methods are based on randomization, so the best solution should be collected for indicating the strong search ability of compared method. As compared to other previous methods, the best solution comparison can be performed by using the total power loss and total cost. Update new solutions for the first half of population using equation (13) Correct violation of new solutions by using equation (18) Calculate fitness function by using equations (19)-(21) Compare new and old solutions to keep higher quality ones Sort all solutions in ascending order of fitness function Calculate average fitness function of the population Fit mean Update new solutions for the first half of population using equation (13) Correct violation of new solutions by using equation (18) Compare new and old solutions to keep higher quality ones

Complexity
(2) e fluctuations of 50 trial runs: the fluctuations can be reflected by observing the mean loss of 50 runs, the standard deviation of 50 runs, and the detail of 50 runs. Methods with lower mean and lower standard deviation can reach more stable search. e low fluctuations of 50 runs can confirm the lower mean and standard deviation obviously. As compared to other previous methods, mean power loss, maximum power loss, and standard deviation can be used for comparison; however, approximately all previous studies have ignored these values for comparisons. For the objective of total cost, the comparisons cannot be done because all previous methods neglected mean and maximum cost.
(3) e convergence speed: the four implemented methods are run by setting the same control parameters (i.e., the same diffusion number, the same population, and the same number of iterations). So, the method that can reach the best solution with better quality is the faster one. In addition, calculation time (CT) is also employed to evaluate the convergence speed. So, the criterion is also applied for comparing with other previous methods. Methods with lower control parameters and/or shorter CT can be faster.

Result Comparison for the 33-Bus Network with Power Loss
Objective. In this section, the first distribution system with 33 buses is employed for placing capacitors with the aim to reduce total power loss. e system is plotted in Figure 2 [38]. e load of the system is comprised of 3715 kW and 2300 kVAr. e whole data of the system are taken from the study [38] and also reported in Table 1. In order to obtain simulation results, four methods consisting of SFS, MSFS1, MSFS2, and the proposed MSFS are run by setting N Df , N Po , and N Iter to 2, 10, and 15, respectively. e two study cases for the system corresponding to the different number of capacitor banks are shown in the following sections.

Case 1: Two Capacitor Banks (Nc � 2).
e summary of the 50 runs including minimum and standard deviation and the mean and maximum loss is, respectively, reported in Figures 3 and 4 while power loss from 50 runs found by the four executed methods is shown in Figure 5. Figure 3 shows that three modified methods can reduce the same total power loss to 141.8433 kW while SFS cannot reach the value although the loss is insignificantly higher. However, the standard deviation of the proposed method is the lowest with 0.78 and that of SFS is the highest with 2.19. As shown in Figure 4, the mean loss of the proposed method is the lowest and equal to 142.0859 kW, whereas that of SFS is 143.86 kW and the highest. e proposed method can reach the second-best maximum loss with 145.9195 kW while SFS is still the worst method with the highest maximum loss with 151.23 kW. Namely, the proposed method can reach less mean loss and less maximum loss than SFS by 1.77 and 5.31 kW, respectively. Clearly, the performance of SFS is the lowest among four implemented methods. Figure 5 can indicate the lowest stability of SFS since points in green have the highest fluctuations. On the contrary, the proposed method with points in blue has the lowest fluctuations and there are many points with the same loss as the best point.
In summary, the comparison of results can indicate the advantages as follows: (1) e newly proposed diffusion technique and the two newly proposed update techniques are really effective. In fact, MSFS1 and MSFS2 can find less minimum loss, less mean loss, and less maximum loss than SFS. In addition, the standard deviation and the fluctuations of MSFS1 and MSFS2 are lower than those of SFS. (2) e proposed method is much effective than its conventional method, SFS. In fact, the three comparison criteria are applied to indicate the statement exactly. e proposed method reaches less minimum loss, so it reaches superiority as shown in the first comparison criterion. e mean, maximum, and standard deviation of the proposed method are less than those of SFS. In addition, the proposed method can find many solutions with the same quality as the best solution. So, the second comparison criterion is confirmed. e proposed method and SFS use the same control parameters, but the proposed method can reach much better results. Hence, the proposed method is faster than SFS as said in the third comparison criterion.
In addition, the proposed method is also compared to previous methods as shown in Table 2. e results indicate that the proposed method can reach the smallest power loss with 141.8433 kW while two other methods reach higher loss with 141.94 and 141.9 kW. As compared to the base network without capacitors, the loss reduction from HM [45], NCB-HM [42], and the proposed method is, respectively, 32.72%, 32.74%, and 32.77%. e calculation time of the proposed method is only 0.382 seconds, whereas the time of the two compared methods was not reported for comparison. It is noted that the time of 0.382 seconds is very small while the two methods have applied analytical methods. Furthermore, HM has used the power loss sensitivity factor to find the most appropriate location and then it has used PSO to search the size of each selected capacitor bank. So, the two methods are hardly ever faster than thẽ  and 7 shows the same characteristic that the height of bars tends to be decreased from SFS to MSFS. is means the improvement is gradually increased from MSFS1, MSFS2 to MSFS, whereas SFS is the worst method with the highest bars. e exact comparisons indicate that the proposed can find less minimum loss, less mean loss, and less maximum

Complexity
proposed method can find eight solutions with the same quality as the best solution, but SFS cannot find even one run with the same loss as the best solution. ere is a main difference of the case with 3 capacitor banks and above case with 2 capacitor banks in finding the best solution. Only MSFS can find the best solution for this case while other ones are searching solutions. But the four methods could find the same solutions for the above case. e fact that this case is more complicated is due to the higher number of control variables (i.e., 3 capacitor banks are corresponding to six control variables, three locations, and three sizes). So, SFS as well as MSFS1 and MSFS2 needs more iterations for finding the same solution as the proposed MSFS method.
For giving more evidences of the high performance of the proposed method, Table 3 is established for comparisons with previous methods. Looking through the minimum loss in the table, we can see that the lowest minimum loss is from MSFS with 138.31 kW while the second-lowest minimum loss is from HM [45] and CPM [44] with 138.37 kW and the highest minimum loss is from BFOA [30] with 144.04 kW. e exact calculation indicates the proposed method can reduce the loss from 0.06 kW to 5.72 kW as compared to other ones. As compared to the base network without capacitor banks, the proposed method can reach the total loss reduction of 34.45% while other ones can reach the reduction from 31.73% to 34.42%. Clearly, the improvement of the proposed method is higher than that of other ones from 0.03% to 2.72%. e comparison of convergence speed is not successful excluding the comparison with BFOA [30] and FPA [33]. Population size and the number of iterations from the proposed method are less than those from the two methods. ey are 10 for the proposed method, but they are 50 for BFOA [30] and 20 and 200 for FPA [33]. e simulation time of FPA was 7.75 seconds, but it is only 0.36 seconds for the proposed MSFS method. In summary, the proposed method can reach better solution than other methods and its speed can be as fast as them or much faster than them. Consequently, the proposed method is an effective one for the study case with 3 capacitor banks in the 33-node distribution system. e voltage profile of the 33-node distribution system for three cases, the first case without capacitor placement (called Base), the second case with 2 capacitors, and the third case with 3 capacitors is plotted in Figure 9. e three curves show the significant improvement of voltage for the cases with capacitor bank placement. e lowest voltage is much increased from less than 0.91 pu of the base network to higher than 0.93 pu of the two cases with capacitors. Furthermore, voltage of approximately all nodes of the cases with capacitors is higher than that of the base network. Clearly, the capacitor placement is necessary for the voltage profile improvement purpose. However, the improvement of the case with 3 capacitors over the case with 2 capacitors is not much since the same nodes from the two cases have approximately the same voltage excluding nodes 23, 24, and 25 where the case with 3 capacitors has better voltage.

Result Comparison for the 69-Bus Network with Power
Loss Objective. In this section, the four executed methods are tested on a larger scale system with 69 nodes as shown in Figure 10 [38]. e total active power and reactive power of all loads are, respectively 3801 kW and 2695 kVAr. e data of the system are taken from the study [38] and also reported in Table 4. In order to solve system with two cases, two capacitor banks and three capacitor banks, the four executed methods are run by setting N Df , N Po , and N Iter to 2, 10, and 20, respectively. e simulation result comparisons are presented in the following sections.

Case 1: Two Capacitor Banks (Nc � 2).
e minimum loss and standard deviation and the mean loss and maximum loss from the four methods are given in Figures 11 and  12. e minimum loss of the study case is similar to the case with 2 capacitor banks of the 33-node distribution system since MSFS1, MSFS2, and the proposed MSFS have approximately the same loss about 146.44 kW while that of SFS is much higher and equal to 146.6221 kW. For other results, the proposed method continues to show outstanding performance because it gets the lowest values of STD, mean loss, and maximum loss. SFS is still the worst method with all highest values. e exact comparison between SFS and the proposed method shows that the proposed method can reach smaller minimum, mean, and maximum loss than SFS by 0.18, 3.45, and 7.33 kW, respectively. SFS is less stable in searching solution since its standard deviation is higher by 1.67. In addition, Figure 13 is also a clear evidence for the improvement of the proposed method over SFS. e fluctuations of the proposed method are much smaller than those of SFS, and MSFS can find more than thirty best solutions, but SFS cannot find the best solution even one time. Consequently, the proposed method outperforms SFS in effectively finding solution and converging faster to the best solution. e proposed method is compared to other ones in Table 5. Because mean loss, maximum loss, STD, and CT for all methods were not reported, the comparison for further investigation failed. However, the minimum loss can be seen from the best performance of the proposed method since it reaches the lowest loss with 146.44 kW while that of the base network, the second-best method, and the worst method is 225, 146.5, and 152.05 kW. Obviously, the proposed method can reduce loss by 78.56, 0.06, and 5.6141 kW, respectively. e proposed method can reduce power loss up to 34.9% compared to the base network without capacitor placement, and its improvement level over other ones is from 0.04 to 3.69%. ese values can show the high performance of the proposed method. NCB-HM [42] and HM [45] have reached results as good as the proposed method, but these methods have been based on the features of networks and convergence speed cannot be measured. PSO and GA are much slower than the proposed method. e simulation time of the two methods was not reported, but the population size and the number of iterations are very high, 30 and 50, while those are only 10 and 20 for the proposed method. So, the proposed method is very useful for the study case.    � 3). For the case, the minimum loss, standard deviation, and the mean loss have the same shape as shown in Figures 14 and 15 e comparison with SFS indicates that the proposed method can reach the reduction of minimum loss, mean loss, and maximum loss by 0.46, 1.53, and 1.12 kW. Figure 16 shows the strong search of the proposed method with many points having smaller loss than other ones, especially compared to SFS. e proposed method can find three best solutions with the lowest loss, whereas SFS cannot find the same quality solution over 50 runs. Consequently, the proposed method still reaches the best performance for the study case because it can find approximately all better results excluding the maximum loss. Table 6 summarizes comparisons of the proposed method with other ones for the study case. e optimal location, size, and the minimum loss of all methods have been reported, but other results such as mean loss, maximum loss, standard deviation, and CT have been reported for only some methods. In addition, control parameters such as population and iterations have not been reported for all methods. e minimum loss of the compared methods is from 145.2 kW to 156.62 kW while the best loss of the proposed method is 145.13 kW. Clearly, the proposed method can reduce the loss by from 0.07 kW to 11.49 kW as compared to other ones and by 79.87 kW as compared to the base network without capacitors. From the power loss reduction, the improvement level of the proposed MSFS method over the worst compared method and the best compared method can be calculated to be 7.34% and 0.05%. Furthermore, the power loss improvement can be up to 35.5% as compared to the original network without the placement of capacitors. Only FPA [33] reported better mean loss, maximum loss, and standard deviation than the proposed method, but it used 20 for population size and 200 for iterations and took much longer calculation time, 18.36  seconds. TLA [34] reported better maximum loss and standard deviation but worse minimum loss, worse mean loss, and longer simulation time together with higher population size and more iterations. e population size and iterations of the method are 50 and 100, but those of the proposed method are only 10 and 20. e values resulted in very high simulation time of TLA with 15.76 seconds, which is much higher than 0.96 seconds of the proposed method. Clearly, the proposed method can be higher than five times faster than TLA through the comparison of control parameters and higher than sixteen times faster than TLA via the comparison of the simulation time. Other methods with the presentation of population and iterations are FGA [26] (60 and 100), GSA [28] (2500 and 750), and MSA [39] (50 and 100), but only MSA [39] reported simulation time of 11.42 seconds. rough the comparison of control parameter settings, the proposed method is, respectively, about 10, 3125, and 8 times faster than FGA [26], GSA [28], and MSA [39]. ese methods are much slower, but they find lower quality of optimal solutions. e result comparison analysis can confirm the following advantages of the proposed method.

Case 2: ree Capacitor Banks (Nc
(1) Reaching better minimum loss than all compared methods. (2) Using smaller population size and smaller number of iterations than metaheuristics, whose control parameters were reported. (3) Spending smaller calculation time than other methods.
As a result, it can be concluded that the proposed method is more suitable than other ones when applied for determining optimal location and size of capacitors in the 69-node distribution system. Figure 17 shows the improvement of node voltage when installing capacitors in the 69-bus distribution system. e blue curve shows the voltage of all nodes before installing capacitors while the green and red curves are voltage of all nodes for the cases with two and three capacitors. e lowest voltage can be seen for the blue curve which is about less than 0.91 pu but that of the green curve and red curve is much higher, about higher than 0.93 pu. At other nodes, the improvement is not as high as at the lowest node voltage, but the improvement is significant. e deviation of voltage is high from bus 10 to bus 29 and from bus 57 to bus 66. e deviation can confirm the high contribution of the capacitors in enhancing the voltage profile. e comparison between the case with 2 capacitors and 3 capacitors reveals that the voltage profile is not much improved when increasing the number of capacitors from 2 to 3. Consequently, if the voltage profile is an objective, 3 capacitors cannot be the best option. is maybe dependent on other issues such as investment capital and operation cost.

Result Comparison for the 85-Bus Network with Power
Loss Objective. In the section, the four methods are run on the largest distribution system with 85 buses and the number of capacitors is increased from 2, 3, 4 to 14. e total active power and the total reactive power of all loads are 2570.28 kW and 2622.08 kVAr.
e system is plotted in Figure 18 [35]. e entire data are taken from the study [35] and also given in Table 7. Due to the large system and the different number of control variables, N Df , N Po , and N Iter are set to different values for different cases. ey are shown in detail in the following sections. � 2). In order to solve the study case with 2 capacitors, the four executed methods are run by setting N Df , N Po , and N Iter to 2, 10, and 25, respectively. All results are reported in  e minimum loss in Figure 19 shows that all the methods have approximately the same best solution; however, the proposed method can reach the best solution with the lowest loss, which is slightly less than that of other ones. e standard deviation, 0.5669 for MSFS and 1.4205 for SFS, from the figure also shows the most stable search of MSFS and the worst stable search of SFS. Figure 20 continues to give the better result of MSFS since its mean loss and maximum loss are the smallest values, and SFS still reaches the highest value for the mean loss and the second-highest value for the maximum loss. e proposed method can reach less minimum loss, mean loss, and maximum loss than SFS  16 Complexity e outstanding performance of MSFS over SFS can be clearly seen through Figure 21. e fluctuations of SFS are higher than MSFS1 and MSF2 and much higher than MSFS. MSFS not only has very smaller oscillations but also reaches many solutions with the same quality as the best solution. e manner can    Minimum loss STD Figure 14: Minimum power loss and standard deviation of 50 trial runs obtained by SFS methods for case 2 of the 69-bus system.

Case 1: Two Capacitor Banks (Nc
confirm the real performance of the proposed method as compared to its conventional version. However, the proposed method's performance cannot be further investigated for the case because previous studies have not considered the case of 2 capacitor banks for the system with 85 nodes. So, the performance investigation can be implemented for the next cases with 3, 4, and 14 capacitor banks. (Nc � 3). In this section, the proposed methods and other SFS methods are applied for the case with three capacitors by setting N Df , N Po , and N Iter to 2, 10, and 25, respectively. Result comparisons for the case are shown in Figures 22-24. e first two figures have the same characteristics since SFS is the worst method with the highest values of loss and standard deviation while MSFS is the best method with the lowest values. e proposed method can reach less minimum loss, mean loss, and maximum loss than SFS by 1.88, 4.62, and 7.42 kW, respectively. Figure 24 shows that all the fluctuations of the proposed method are much lower than those of SFS and approximately all points of the proposed method have lower loss than those of SFS. In the figure, it is hard to indicate a better method between MSFS1 and MSFS2, but the two methods are more effective than SFS. e proposed method finds three best solutions and many local optimal solutions with slightly less loss than these best solutions while SFS cannot converge to nearby global optimal solutions or the same local optimal solutions. Clearly, the new diffusion      Table 8. e loss of the base network is 316.12 kW and that of the proposed method is 150.8555 kW while that of BFOA [31] and FPA [33] is 152.25 and 151.807 kW. e loss indicates that the proposed method can reach less loss than the base network and these methods by 165.3445, 1.3945, and 0.9515 kW, respectively. e loss reduction is, respectively, equivalent to 52.3%, 0,92%, and 6.3%, and these values are also better effectiveness levels of the proposed method. FPA reported better mean loss, maximum loss, and standard deviation, but it used much higher values for control parameters and took longer simulation time. e method has used 20 and 200 for population and iterations and spent 18.36 seconds for simulation time while these parameters are only 10 and 25 and simulation time is 2.5 seconds for the proposed method. Clearly, FPA could reach better stability thanks to the higher population size and more iterations. However, its best solution is worse than that of the proposed method and its search speed is slower. So, FPA cannot be a competitive method of the proposed method. BFOA [31] did not report mean loss, maximum loss, and calculation time for comparison. But the comparison with the method can easily to lead to a conclusion because it used higher population size and much more iterations and reached higher loss. So, the proposed method is superior to BFOA and FPA for the study case. � 4). In this section, the number of capacitors is increased to 4, leading to the increase of control variables to 8, 4 control variables for location, and 4 control variables for size of capacitors. So, the four methods are executed by setting N Df , N Po , and N Iter to 2, 20, and 30, respectively. e simulation results of the case are shown in Figures 25-27. e results of minimum loss, standard deviation, and mean loss have the same characteristic with the effective order of SFS, MSFS1, MSFS2, and MSFS. Figure 27 shows the highest loss of points from SFS and many points with lower loss from the proposed method. For the large study case, the proposed method can reach less loss than SFS by 2.03 kW equivalent to 1.35%, and approximately all runs of the proposed method have less loss than those of SFS. Clearly, MSFS is still the best method with much better performance than SFS. e comparison with MINPM [27] is shown in Table 9. e method only reported the minimum loss with 159.87 kW, which is higher than 148.7077 kW of the proposed method by 11.1623 kW. e less loss is corresponding to the improvement level of 6.98%. Clearly, the proposed method is much more effective than MINPM [27]. e proposed method's simulation time is 7.7 seconds while that of MINPM [27] was not reported for comparison. e minimum loss can confirm the outstanding performance of the proposed method as compared to MINPM [27] for the study case. (Nc � 14). In this section, the most complicated case with fourteen capacitors is considered for testing the real performance of the proposed method and other ones. e number of control variables is 28, which is much higher than all above study cases. So, the setting of N Df , N Po , and N Iter is 2, 20, and 75, respectively. Figures 28-30 show that the proposed method is the best method and SFS method is the worst method. Namely, the proposed method finds less loss than SFS bỹ   Table 10.

Complexity
e proposed method can reach less minimum loss and shorter simulation time than TLA [34]. e two factors are 142.4999 kW and 15.72 seconds for the proposed method and 143.18 kW and 18.38 seconds for TLA [34]. e proposed method can reduce power loss more effectively than TLA by 0.68 kW corresponding to 0.47%. For other results such as mean loss, maximum loss, and standard deviation, TLA is more effective; however, it should be noted that TLA was run by setting 50 and 100 to population and iteration while the values of the proposed method are 20 and 75. e stability of search from TLA is better because it used higher control parameters, but the best solution of TLA is worse than that of the proposed method. So, it can be concluded that the proposed method is very effective for the study case. Figure 31 shows the improvement of voltage profile after installing different number of capacitors in the 85-bus distribution system. e voltage profile of the base system without capacitors in blue is the worst among five cases of voltage. Voltage profile of four other cases is much better than that of the base network. e lowest voltage of base network is about less than 0.88 pu but that of other cases is higher than 0.92 pu. Furthermore, the voltage profile from node 25 to node 85 in the base network without capacitors is under 0.9 pu, but it is higher than 0.92 pu in the modified network with capacitors. e comparisons between the case without capacitors and other cases with capacitors reveal that the voltage profile is not much improved when installing capacitors. However, the voltage of the case with 3 and 4 capacitors is not much better than that of the case with 2 capacitors. Similarly, voltage profile of network with 14 capacitors is not much better than that of the case with 3 and 4 capacitors. is manner can be understood easily because the purpose of the capacitor placement in the paper is for power loss reduction rather than for voltage profile improvement.

Result Comparison for the 69-Bus Distribution Network with the Objective Function of Total Cost.
In this section, the total cost of energy loss over a year and capacitor investment is considered as objective function as shown in equation (2). e number of load levels, M, is three in which the number of hours in the first load level (called light load), in the second load level (called normal load), and the third load level (called peak load) is, respectively, 2000, 5260, and 1500 hours. e light load, normal load, and peak load are corresponding to 50%, 100%, and 160% of the current loads in the 69-bus distribution network [46]. e main difference of the total cost objective function and the power loss objective function is the size of capacitors. e size is a continuous variable in the previous section, but it is a discrete variable in this section. e size of capacitors is selected from 150 to 4050 with a step size of 150 kVAr. e price of each kWh and each kVAr is $0.06 and $3.0, respectively [46]. e data of the system are also taken from Table 4. For reaching the total cost of the 69-bus system with the placement of three capacitors, SFS methods are run 50 independent trials by setting the number of diffusions, the size of population, and the number of iterations to 2, 10, and 30. e results obtained by the four SFS methods are reported in Table 11. e table can show the best performance of the proposed method since it finds the most optimal solution for placing three capacitors in the system. It gets the lowest minimum cost, but SFS suffers the highest cost. MSFS has achieved a lower cost than SFS about $25.9, corresponding to 0.03 % of the cost of SFS. e results from the proposed method are compared to those of TLA [34] and DSOA [46], which are shown in Table 12. In addition, location and size of each capacitor, power loss of the whole system corresponding to three load statuses including light load (LL), normal load (NL), and peak load (PL), and total compensation reactive power are also reported in the table. As shown in the total compensation capacity row, MSFS and TLA [34] use the same 2100 kVAr and much less than 3600 kVAr of DOSA. So, the proposed method and TLA [34] expend the same capacitor investment cost of $6300 while DSOA uses much higher capacitor cost, which is $10,800 and higher than that of MSFS and DSOA [46]   power loss, price of each kWh, and the number of hours, the energy loss cost can be obtained for the three methods, and it is equal to $89,016.38 for TLA, $89,112.60 for DSOA [46] and $87,855.44 for MSFS. By employing equation (2), total cost is calculated and reported for evaluating the real performance of the proposed method. As a result, the proposed MSFS is the best method with the lowest total cost of $94,155.43 while that of TLA and DSOA is much higher and equal to $95,316.38 and $99,912.60. Clearly, the selection of capacitor can support the proposed method save $1,160.95 and $5,757.17 as compared to TLA and DSOA. e saving cost is corresponding to 1.22% and 5.76% of the total cost from TLA and DSOA. e saving cost and the improvement level are significant when implementing the proposed MSFS for the purpose of placing capacitors in the 69-bus distribution network instead of using TLA and DSOA. Furthermore, the proposed method can find more effective solutions much faster than TLA. In fact, N Po and N Iter are set to 50 and 100 for TLA but only 10 and 30 for the proposed MSFS method.   Table 7 while other information such as price of energy, price of capacitor, time of each load status, and load levels is the same as that of the 69-bus distribution network. Summary of total cost from 50 independent trial runs by implementing four SFS methods is reported in Table 13. MSFS is still the best method, and SFS is always the worst method with the highest minimum, mean, and maximum cost. e proposed method can find less cost than SFS by $2,286.24, which is corresponding to 2.31% of the total cost from SFS. e improvement level of the system is 2.31% while that of the 69-bus distribution network is only 0.03%. Obviously, the outstanding performance of the proposed method is much clearer for larger scale system.

Result Comparison for the 85-Bus Distribution
Results from the proposed method and two compared methods such as TLA [34] and DSOA [46] are presented in Table 14. Location and size of each capacitor corresponding to LL, NL, and PL are optimal solutions, and they can support to calculate power loss for the three load statuses. e power loss is (34.11, 143.2493, 411.6) kW for TLA [34], (34.76 144.01, 410.69) kW for DSOA [46], and (34.943, 144.871, 415.389) kW for the proposed method. By using the power loss, the time of each load level, and the price of energy, the energy loss cost of TLA, DSOA, and the proposed can be obtained, and it is equal to $86,346.679, $86,586, and $87,299.544, respectively. Although the energy loss cost of the proposed method is higher than that of

Conclusion
In this paper, the determination of the most appropriate site and size of capacitor banks has been implemented by the proposed MSFS method for reaching two single-objective functions which were the reduction of power loss and the reduction of total cost of energy loss and capacitor investment. e proposed method has been developed by performing the new diffusion mechanism and two new update mechanisms on the conventional SFS method. ree systems with 33, 69, and 85 buses have been solved by the proposed method, conventional SFS, and two other improved versions of SFS. e results from the proposed method, SFS, and other previous methods were compared for evaluation of performance. For the comparison with SFS, the advantages of the proposed method over SFS can be summarized as follows: (1) For the first system with 33 buses, SFS and the proposed method could find the same best solution, but for the second study case of the system, SFS could not find the same best solution as the proposed method. e proposed method could find less loss than SFS by 0.263 kW corresponding to 0.002% of power loss. (2) For the second system with 69 buses, the proposed method could find less power loss than SFS by 0.182 kW and 0.457 kW, which were corresponding to 0.001% and 0.003% of loss from SFS. For the more complicated objective with the sum of energy loss cost and capacitor purchase cost, the proposed method could reduce the total cost to $25.9 corresponding to 0.03% of cost from SFS. (3) For the third system with 85 buses, the improvement of the proposed method over SFS was more significant for four study cases. For the first case with two capacitors, the loss reduction was 0.0294 kW corresponding to 0.0002%, but for the following cases with higher number of capacitors, the reduction was higher. In fact, the loss reduction was, respectively, 1.8781, 2.0274, and 2.6125 kW corresponding to 0.0123%, 0.01345%, and 0.018% for the cases with three, four, and fourteen capacitors. When considering the sum of energy loss cost and capacitor purchase cost, the proposed method could reach much higher performance than SFS by reducing the total cost to $2,286.24, which was corresponding to 2.31%.
ese values of improvement level, which was represented by percent, were not high because SFS was also an effective method. Here, the performance improvement of the proposed method was about the ability to find the best solution and speed up the search process. e proposed method could find many solutions with the same quality as the best solution, whereas SFS could not reach the best solution even one time over fifty trial runs. When competing against other previous methods in solving the capacitor placement, the proposed method could find either the same good solutions or much better solutions. From the simulation results, the superiority of MSFS over other methods is as follows: (1) For the 33-bus system, the power loss of the proposed method could be equal to that of a few methods but much smaller than many ones. e power loss reduction of the proposed method could be 0.07% and 3.98% of compared methods.
(2) For the 69-bus system, MSFS could achieve the improvement level over other ones from 3.7% to 7.3% for power loss and from 1.22% to 5.76% for the sum of energy loss cost and the capacitor purchase cost. (3) For the largest system with 85 buses, MSFS could reach the improvement level from 0.92% to 6.98% for power loss and from 0.87% to 3.77% for the sum of energy loss cost and the capacitor purchase cost.
In terms of convergence speed, the proposed method was always much faster than other methods by setting smaller population size and smaller number of iterations. e parameter settings and simulation time can be summarized as follows: (1) For the 33-bus system, the population size and iteration were set to 10 and 15 for the proposed method but 20 and 200 and 50 and 50 for other methods.
(2) For the 69-bus system, these parameters were 10 and 20 and 10 and 30 for the proposed method but 30 and 1000, 60 and 100, 50 and 100, and 2500 and 750 for other methods. (3) For the largest system with 85 buses, the proposed method was run by setting 10 and 25, 20 and 30, 20 and 75, and 30 and 100 for different cases, but they were much higher and equal to 50 and 100, 100 and 50, and 20 and 200 for other methods. (4) e simulation time of the proposed method was shorter than 1 second for the two first systems and from two to fifteen seconds for the last system, whereas the time of other methods was about ten seconds for the first two systems and about from ten to twenty seconds for some cases of the last system.
In summary, the proposed method could reach much better solutions than previous method, use less control parameter, and spend shorter simulation time. So, it can be concluded that the proposed method is an effective method for application of determining site and size of capacitors in distribution networks.

Nomenclature
Nbr: Number of branches of the considered distribution network Nb: Number of buses in the considered distribution network Nc: Number of capacitor banks

Data Availability
Data of the 33-and 69-node distribution systems are taken from [38], and data of the 85-node distribution system are taken from [35].

Conflicts of Interest
e authors declare that they have no conflicts of interest.