Application of the SSA for Optimal Reactive Power Compensation in Radial and Meshed Distribution Using D-STATCOMs

: This paper deals with the problem regarding the optimal placement and sizing of distribution static compensators (D-STATCOMs) in radial and meshed distribution networks. These grids consider industrial, residential, and commercial loads within a daily operation scenario. The optimal reactive power ﬂow compensation problem is formulated through a mixed-integer nonlinear programming (MINLP) model. The objective function is associated with the minimization of the expected energy losses costs for a year of operation by considering the investment costs of D-STATCOMs. To solve the MINLP model, the application of a master–slave optimization approach is proposed, which combines the salp swarm algorithm (SSA) in the master stage and the matricial backward/forward power ﬂow method in the slave stage. The master stage is entrusted with deﬁning the optimal nodal location and sizes of the D-STATCOMs, while the slave stage deals with the power ﬂow solution to determine the expected annual energy losses costs for each combination of nodes and sizes for the D-STATCOMs as provided by the SSA. To validate the effectiveness of the proposed master–slave optimizer, the IEEE 33-bus grid was selected as a test feeder. Numerical comparisons were made against the exact solution of the MINLP model with different solvers in the general algebraic modeling system (GAMS) software. All the simulations of the master–slave approach were implemented in the MATLAB programming environment (version 2021b). Numerical results showed that the SSA can provide multiple possible solutions for the studied problem, with small variations in the ﬁnal objective function, which makes the proposed approach an efﬁcient tool for decision-making in distribution companies.


Introduction
Distribution networks play a fundamental role in electric power systems, which is why reducing their losses and minimizing investment costs are among the greatest challenges faced by current systems [1]. When the voltage level increases, the current circulating through the network decreases, therefore reducing the losses caused by the Joule effect. In the particular case of distribution networks, the voltage is reduced in such a way that the current circulating through the network is higher. This effect generates greater losses, despite being necessary, given that reducing the voltage levels has a significant impact on distribution networks' design and investment costs. In the same way, these losses have an incidence on the power factor, given that the generators are responsible for covering the total demand (i.e., including these losses). This causes greater uncertainty in energy planning, which increases the design difficulty of a distribution network. The above has an Algorithms 2022, 15, 345 2 of 16 effect on the overruns that must be assumed by the network operator and, therefore, by the users [2].
Distributed generation (DG) may be regarded as the best option for reducing power losses. However, its initial installation costs can be very high, especially when compared to strategies such as the reconfiguration and connection of branch capacitors [3]. The main issue with capacitor banks is that they inject power in fixed reactive power steps, which does not contemplate a dynamic daily active and reactive power demand scenario along electrical distribution networks [4].
The use of power compensation techniques by means of devices that implement power electronics has allowed for the development of strategies that consider system dynamics in order to present increasingly efficient solutions in the search for loss reduction [5]. In this particular case, a focus on power electronics devices was chosen, specifically on D-STATCOMs (Distribution Static Compensators), which are connected in parallel to the network [6]. Among the main features of these compensators, the reactive power absorption or injection by varying the phase amplitude and angle with respect to the connected network can be highlighted. All this, as a function of the demand behavior [7]. These devices may be located at any node of the network. Nevertheless, it is necessary to determine their location and sizing, so that a considerable decrease in network losses may be achieved. The implementation of D-STATCOMs has the following advantages: (i) high reliability, (ii) low operating costs, and (iii) durability (typically between 5 and 15 years) [5].
The specialized literature presents different reports on the optimal siting and sizing of D-STATCOMs in distribution systems. Some of these works are described below. The authors of [8] proposed a multiobjective particle swarm optimizer that considers the reconfiguration of the electrical network in order to locate and size D-STATCOMs while analyzing the minimization of power losses (objective function), voltage stability, and the loadabilty factor of distribution lines. The main aspect of the aforementioned article is that the optimization process is only carried out under peak load conditions, which could lead to oversizing in the D-STATCOMs.The authors of [9] present a combination between a genetic algorithm and a mathematical programming method based on second-order cones in order to place and size D-STATCOMs in distribution networks. The main contribution of this work is the utilization of industrial, residential, and commercial demand curves. The results obtained in the IEEE 33-and 69-node systems show the effectiveness of the proposed methodology in comparison with the solution of the exact model in the GAMS software.
The authors of [7] present a complete review of the specialized literature with regard to the optimal placement and sizing of D-STATCOMs in distribution systems. The main contribution of this work is the classification of the different existing methodologies for solving the problem. Among these methodologies, the following stand out: analytical methods, neural networks-based methods, combinatorial optimization algorithms, sensitivity algorithms, and the combination between metaheuristic methods and sensitivity indicators. The authors of [10] present the combination of a Chu and Beasley genetic algorithm (CBGA) and a specialized power flow algorithm for distribution networks, with the purpose of siting and sizing D-STATCOMs in distribution networks. The main contribution of this methodology is the possibility of working with radial and meshed distribution networks. The numerical results of the proposed CBGA in the IEEE 33-and 69-node systems showed greater reductions in the network's annual operating costs in comparison with the solution of the exact model in the GAMS software. Other methodologies presented in the specialized literature for the optimal integration and sizing of D-STATCOMs in distribution systems include the particle swarm optimization algorithm with a discrete-continuous coding [11], bat-based optimization [12], methods based on exact mathematical programming from quadratic next-integer programming models [13], and mixed-integer second-order cone programming [5]. Table 1 summarizes the main contributions made by the literature regarding the problem of the optimal placement and sizing of D-STATCOMs in distribution networks. Mixed-integer second-order cone programming Power losses minimization and investment and operating costs reduction [5] 2022 The main characteristics of the optimization approaches in Table 1 are as follows: (i) the most common objective function is power losses minimization; (ii) there are two tendencies regarding the solution of the D-STATCOM optimal location and sizing problem: the first involves metaheuristic optimization techniques, and the second one proposes convex formulations or combines convex and combinatorial methods; and (iii) during the last decade, the problem regarding reactive power compensation in electrical systems continues to be a research area of interest for academia and the industry.
The main contributions of this research are listed below: i. The optimal integration of D-STATCOMs in distribution systems by applying the SSA through a discrete-continuous coding. ii. The combination of the SSA and the specialized power flow algorithm for distribution networks known as the generalized backward/forward method, which allows working with radial and meshed topologies with no modifications to its iterative formula.
It is important to mention that the combination between the SSA and generalized iterative sweeping constitutes a new master-slave optimization strategy that is applicable to the placement of parallel compensation devices in radial and meshed distribution systems, which has not been previously proposed in the specialized literature.
Note that the selection of the salp swarm algorithm (SSA) to determine the best nodal locations for the D-STATCOMs in distribution networks was based on the following facts: (i) the SSA is a recently developed optimization methodology with excellent numerical performance in both nonlinear and mixed-integer optimization problems [18], which has not previously been applied to locate and size reactive power compensators in distribution networks; and (ii) the computational implementation of the SSA is simple and only requires some abilities in mathematical programming. In addition, it has a few evolution rules that simplify its numerical validation in any programming language.
The rest of this work is structured as follows. Section 2 presents the mathematical formulation of the problem regarding the optimal integration of D-STATCOMs in distribution systems by means of a mixed-integer nonlinear programming model. Section 3 describes the solution methodology based on a master-slave optimization algorithm. In the master stage, the SSA is implemented, which defines the nodes for placing the D-STATCOMs, as well as their optimal size, by means of a discrete-continuous coding. In the slave stage, a technical evaluation of the solution is carried out by solving a multi-period power flow via the generalized backward/forward method. Section 4 presents the main characteristics of the test system, namely the IEEE 33-node system, which includes residential, industrial, and commercial demand scenarios. Section 5 presents the simulation results along with their analysis and discussion. Finally, Section 6 describes the main conclusions derived from this research.

Mathematical Modeling
This section of the article presents the model for the optimal placement and sizing of D-STATCOMs for radial and meshed distribution networks. It must be clarified that the problem is of the MINLP type, which means that, in order to represent the notes where the D-STATCOMs will be located, discrete variables will be employed, as well as continues variables for their size. Moreover, these variables have a nonlinear relationship due to the power balance equations.

Objective Function
The objective function comprises the sum of two functions: one of them represents the operating costs due to energy losses ( f 1 ), and the other one represents the investment costs of the D-STATCOMs. Next, the way in which Equations (1) and (2) define each objective function is shown: Equation (1) defines the annual operating costs due to power losses in the distribution network, where, starting by the left C kWh is the average cost of energy in kilowatts per hour; T is a constant associated with the period under study (365 days in this case); Y km is the magnitude of the admittance matrix that relates k nodes to m nodes with an angle θ km ; V kh and V mh are the voltages that relate nodes k and m in period h with voltage angles denoted as δ kh and δ mh , respectively; ∆ h represents the time associated with the power flow (1 h); and H and N refer to the set of time periods and nodes in the network, respectively. On the other hand, Equation (2) represents the objective function corresponding to the investment costs of the D-STATCOMs, where k 1 and k 2 greater than zero are constant. The first represents the annual investment entailed by the D-STATCOM and the second refers to the device's useful life. The parameters α, β, and γ are scalars greater than zero that correspond to the installation costs of the D-STATCOMs with nominal power and generation capacity Q DS k . Thus, the general objective function is proposed as defined in Equation (3).
Note that A cost represents the expected annual cost of the distribution network, which corresponds to the sum of the power losses and installation costs of the D-STATCOMs.
The main feature of the objective function defined in (3) and its components f 1 and f 2 in (1) and (2) is that both components (i.e., f 1 and f 2 ) are nonlinear, non-convex functions of the main decision variables. The former is nonlinear and non-convex due to the presence of products among voltages and trigonometric functions [19], and The latter is nonlinear and non-convex with a mixed-integer structure; its non-convexity is governed by the presence of a cubic function in its formulation [20].
Note that, due to the complexities of the objective function in Equations (1)-(3) and the large dimensions of the solution space regarding the discrete variables, most of authors deal with their solution by applying metaheuristic optimization methods [7], which is the case of this research (SSA), as will be detailed later in this document.

Set of Constraints
The constraints for the problem regarding the optimal location and size of D-STATCOMs in distribution networks are given by their maximum number and size, the voltage limits, Algorithms 2022, 15, 345 5 of 16 and the active and reactive power balance. This set is presented below. It is important to mention that Equations (4)-(7) are defined for ∀k ∈ N and ∀h ∈ H: Equations (4) and (5) represent the active and reactive power balance for each node and period. The most relevant aspect of Equation (5) is the term Q DS k , which represents the reactive injection provided by the D-STATCOM in the nodes where it will be allocated. The terms P g kh and Q g kh are the active and reactive power injected by the generation connected to node k in a period of time h. Given that this article uses a demand curve that varies over time, it is highlighted that the terms Q d kh and P d kh are the dynamical loads associated with node k at time h. On the other hand, inequality (6) defines the maximum upper and lower voltage limits at each node. This constraint is closely related to that defined by the relevant authorities of the electrical sector. Equation (7) defines the size limit of a D-STATCOM. Z k is a binary variable (i.e., it takes values of 0 or 1) that indicates the presence or not of a D-STATCOM at node k. The last inequality is the theoretical limit of D-STATCOMs to be installed in the distribution network, where N DS A corresponds to the maximum number of D-STATCOMs that will be able to be installed in the distribution network.

Solution Methodology
In order to solve the MINLP problem regarding the optimal placement and sizing of D-STATCOMs in distribution networks, dividing the problem in two stages is proposed, i.e., by means of a master-slave programming strategy [21,22]. In the master stage, the SSA defines the nodes and locations of the D-STATCOMs, while the slaves stage, via generalized iterative sweeping, evaluates the multi-period power flow and determines the network's operating costs, i.e., the costs of the energy losses. Next, each step of the proposed optimization methodology are detailed.

Slave Stage: Generalized Backward/Forward Method
The sweeping method is a numerical method based on graph theory that uses the incidence matrix as a means to solve the power flow problem [23]. The incidence matrix is a representation of the relationship between the nodes and branches of the distribution system [24]. To create an incidence matrix (node-branch), the following structure is employed: If the current leaves node i −1 If the current reaches node i 0 If the branch j is not connected to node i    To develop the problem, this matrix is divided into two submatrices representing generation (slack node) and demand.
With the above matrices, it is possible to obtain the relationship between the voltage drops in branches E and the voltages at the nodes, as follows: where V s is the voltage at the slack node, and V d is the vector comprising the voltages at the nodes where there is demand. This is thus formulated because the power flow problem revolves around finding the variable V d . Now, Ohm's law is applied in order to relate the branch current J to voltages E: with Z being the matrix that contains all the primitive impedances of the lines. Now the net injected currents I can be separated as follows as a function of the slack node current I s and the current at the demand nodes I d : Now, Kirchoff's current law is applied in order to find the relationship between the net injected currents at the system nodes I and the branch currents J, with which the following structure is obtained in matrix form: Note that, by considering Equation (14) and substituting Equations (10) and (11), a general expression for the demanded current is obtained: with Y being the primitive admittances matrix, i.e., Y = Z −1 .
On the other hand, it is known that the relationship between the demanded currents and the voltages has a nonlinear hyperbolic form [25], which is defined below.
Now, by equalizing Equations (15) and (16), a general equation is obtained which relates the demand node voltages between each other in a nonlinear way.
finally, Equation (17) is adjusted in order to make it iterative and achieve that the expected voltage value is reached with the predefined convergence error via loop programming logic. The iterative sweeping formula is presented below.
Note that the convergence error is defined as ε = 1 × 10 −10 and the iterative process starts with t = 0 and V 0 d = 1∠0 • , which ends if and only if max ||V t+1 d | − |V t d || ≤ ε.

Master Stage: Salp Swarm Algorithm
The master stage is entrusted with modifying the network by implementing the D-STATCOMs. This paper suggests solving this stage of the problem by using the salp swarm algorithm, which is explained below.
Salps are organisms that move by contraction as they pump water through their gelatinous body. The pumped water used for moving also allows them to obtain their food, which consists of phytoplankton filtered from water [26]. When salps look for food, they do so in swarms, which are also called chains. These chains significantly improve Algorithms 2022, 15, 345 7 of 16 efficiency when looking for food. This model was imitated and proposed as a combinatorial optimization algorithm by [27] in 2017.

Generating the Initial Population
Equation (19) is used to generate the initial population of the salp chain S k , where each salp in the population is a possible solution to the problem, i.e., it defines the node and size of each D-STATCOM available for installation in the distribution system.
where sub-index k represents the salp chain k, and sub-index j represents element j in the salp chain (variable of interest -node-or -size-), with d being the number of variables in the problem. The values generated by Equation (19) within the solution space are defined by X min and X max , which were mentioned in the problem's set of constraints. This initial population can be represented by a matrix of size nxd, where n is the number of salps as possible solutions to the problem. Equation (20) defines the initial population employed in the SSA.

Calculating the Objective Function
Each individual contained in M S (matrix row) is evaluated in the objective function with the help of the slave stage of the algorithm's solution. This evaluation is in turn organized in a vector of size nx1, which is denoted as MO S . It is important to remember that, given the problem constraints, it is guaranteed that the salps are within the solution space, (i.e., each candidate solution S k is checked in order to ensure that the nodes are integers and the sizes of the D-STATCOMs are within the permitted range [28]).
. . . f n ([S n,1 , S n,2 , · · · , S n,d ]) Once the entire population has been evaluated and built, the vector MO S organizes the evaluated results from lowest to highest, as the main proposal of this article is solving a minimization problem. Note that this organization is also applied to the initial population, with the purpose of preserving the correspondence between the candidate solution S k and its objective function f k (S k ).

Salp Chain Movement
The proposed salp swarm algorithm employs an iterative process, i.e., with each iteration, a new chain is generated and the salps are divided into two: the leader and the followers [27]. The leader is the head of the chain, while the others follow it. The leader is responsible for leading the swarm throughout the solution space [29]. The salp with the highest impact on the network (i.e., the lowest objective function) is chosen as the leader S l . with S 1 being the first row of the matrix M S , which has been organized from lowest to highest with respect to the objective function value. The leader is updated in each iteration outside the iterative algorithm. Moreover, in the first cycle, the leader is chosen as food (phytoplankton). This is due to the fact that the leader is the one that finds the area with the best food, which translates to the solution space, and it is within this same space that advancing with the chain is decided. The phytoplankton is represented by a vector of dimension 1xd, as it stores the information of each chain leader. Note that there will be as many leaders as variables d in the problem under study.

Modifying the Population by Means of the Leader
Within this first strategy, the initial population is sectioned in two equal parts. The first half is updated using the information of the leader, as expressed in Equation (23).
Equation (23) is made up of the actual position of the salp in the variable S (i,j) and a set of control constants that allow the salp chain to be exclusively within the solution space. The equation is the movement mechanism used by the first part of the salp chain, where the position is updated in S i,j . C 2 and C 3 are constants with a uniform distribution within a range of 0 to 1 which limit the movement of the chain in order to guarantee that it is in the best problem solution space. In addition to these constants is C 1 , which is entrusted with controlling the exploration and exploitation of the solution within the movement, as given by the following equation.
where the variables l are the program's current iteration, and L is the maximum number of iterations.

Modifying the Population Using NEWTON's Movement Laws
For this part of the algorithm, the chain moves using Newton's third law, which is represented via the following equation: This movement mechanism involves the salps sharing information with each other in order to move through the solution space and generate new potential solutions to the nonlinear programming problem. However, with this movement, it must be verified that the variable limits are still fulfilled. If not, the necessary corrections for the new possible candidate solution.

Remark 1.
Note that the SSA approach evolves through the solution space by using the evolution rules proposed in [27], which were designed for continuous optimization problems. However, this research applied the same rules, and, after each iteration, the component of the coding regarding the nodes where the D-STATCOMs must be placed is rounded as recommended in [30] for optimization problems with discrete variables.

Summary of the Proposed Solution Methodology
With the purpose of illustrating the application of the SSA to general optimization problems, the Algorithm 1 is employed. Define the values of the constants C 2 and C 3 ; Check the feasibility of each individual; 17 for k = n 2 + 1 : n do 18 for j = 1 : d do 19 S i,j = Check the feasibility of each individual; 21 Report the best solution S best Additionally, in order to clarify the application of the SSA approach to locate and size D-STATCOMs in radial and meshed distribution networks, the flow chart of the proposed methodology is presented in Figure 1.

Is new leader better than after phytoplankton?
Define the size of the initial population NO Aprende más sobre esta plan Usa este diagrama de flujo básico pa -Definir y visualizar claramente lo pasos de tu proceso. -Analizar y mejorar tu proceso. -Acceder a la biblioteca de figura de diagramas de flujo de Lucidchart.

Test System
To evaluate the effectiveness of the proposed solution methodology for placing and sizing D-STATCOMs in medium-voltage distribution systems, the IEEE 33-node system is employed. The electrical configuration of this system is presented in Figure 2. As for the parametric information of the network, the peak consumption data at each demand node are presented in Table 2. Table 2. Impedance and demand data for the IEEE 33-node system.

Node i
Node j R ij (Ω) X ij (Ω) P j (kW) Q i (kvar) Since this project is focused on selecting reactive compensators for radial and meshed distribution systems with varying demand according to the type of consumers across the distribution network, in order to evaluate these operating conditions, the demand data are considered according to that reported in [10], as shown in Table 3. Table 3. Load consumption profiles per user type: residential (Res.), industrial (Ind.), and commercial (Com.).

Hour
Ind Note that the selection of these power consumption profiles was based on the typical load variation in a distribution circuit in the metropolitan area of Bogotá, Colombia. These values were estimated after multiple measurements taken during work days in a distribution substation operated with 11.4 kV. No geographical information is provided due to confidentiality commitments made with the energy distribution company.

Results and Discussion
To validate the proposed solution methodology, the solution of the optimization model for placing and sizing D-STATCOM in electrical distribution systems by using the AESDC was implemented in MATLAB's R2021b version on a computer with an Intel Core i7-10750H @2.6 GHz and 16.0 GB of DDR4 RAM at 2300 MHz running the 64-bits version of Microsoft Windows 10 Home.

Radial Configuration Results
To verify the efficiency of the proposed solution methodology, 20 evaluations were performed for the radial configuration. The best result obtained had a cost of 108,249.36. Furthermore, this solution had a repeatability higher than 50% of the total evaluations. The best solution found for this test system corresponds to nodes 13, 25, and 30, with D-STATCOMs sizes of 258.51, 105.47, and 538.02 kvar, respectively. Note that, in this simulation, 20 individuals and 1000 iterations were set to evaluate the performance of the SSA.
The best five solutions reached by the methodology proposed for the IEEE 33-node system with a radial topology are presented in Table 4. The criteria used to select the five best solutions is the lowest operating cost using SSA without considering the execution time. Note that the fifth result in Table 4 corresponds to a value of 108,428.91 dollars per year and defines the D-STATCOM locations at nodes 14, 10, and 30, with sizes of 180.37, 82.17, and 605.97 kvar, respectively. In this way, it can be stated that, for both the best and the worst solution, there is an optimal location at node 30. Furthermore, the difference between the two is 179.55 dollars, which means an approximate error of 0.166%. Even though this is a low error, it can be said that the gap between the two demonstrates that the solution space is large, despite the fact that they converge to a similar result.
The processing times for each algorithm evaluation regarding the best and worst result are 2.4317 and 2.0876 min, respectively. Which is an adequate computational effort due to the complexity of the studied problem. It is also important to mention that the nodes at which the D-STATCOMs are optimally placed are numbers 30 (commercial), 25 (industrial), and 13 (industrial).

Meshed Configuration Results
As with the radial topology, a population of 20 individuals (salps) was used to test the developed algorithm, as well as 20 consecutive evaluations of the whole methodology, considering 1000 iterations per evaluation. Table 5 summarizes the best five solutions obtained for the meshed IEEE 33-node system.  Table 5). The difference between the best and worst result is 26.50 dollars. Note that this difference is very small in comparison with the objective function values, which demonstrates that the proposed methodology can find very high-quality solutions in each evaluation.
Note that the processing times for each algorithm evaluation regarding the best and worst result are 1.9360 and 1.9392 min, respectively, i.e., the algorithm has a tendency towards finding the space and solution of the problem in similar periods of time.
The next section only compares the proposed algorithm within the meshed configurations, given that, in the specialized literature, the use of industrial, commercial, and residential curves in the studied problem has not yet been presented for radial configurations.

Comparative Analysis: Meshed Configuration
For the particular case of the meshed network, the baseline case presented in [10] was employed, as well as the XPRESS, SBB, DICOPT, and LINDO optimization algorithms available in the General Algebraic Modeling System software. The percentage differences for the aforementioned algorithms, in comparison with the proposed solution methodology, are favorable by 11.57, 2.13, and 1.90%, respectively. Table 6 shows the figures associated with the percentages describing the differences between the implemented solutions. Finally, it is important to mention that the number of salps in the population was set as 20 in the radial and meshed configuration cases because of the compromise between the total processing time and the quality of the optimal solutions reached. This number of individuals was reached after constructing a meshed grid with sizes of populations between 5 and 30, evaluating the total processing time and the variability of the final solution after 100 consecutive iterations per population size.

Conclusions and Future Works
This paper studied the problem regarding the optimal placement and sizing of D-STATCOMs for electrical distribution networks with radial and meshed configurations. The optimization problem was formulated using an MINLP model, where the binary variables were associated with the nodes where the D-STATCOMs must be placed. The continuous variables were related to the power flow solution, i.e., the voltages, currents, and power variables. The objective of the MINLP model was to reduce the expected annual energy loss costs after installing D-STATCOMs across the distribution network. The main contribution of this research is the consideration of three different load profiles to classify the distribution grid as a function of energy usage (residential, industrial, and commercial). To solve the MINLP model, a master-slave optimization strategy was proposed. The master stage was entrusted with defining the nodes where the D-STATCOMs will be located and their corresponding sizes. In the slave stage, an efficient power flow formulation was implemented, which allowed calculating the expected annual grid operating costs. Numerical results in the IEEE 33-bus grid with radial and meshed configurations confirmed the effectiveness of the proposed solution methodology in comparison with the MINLP solvers available in the GAMS software.
An important result when comparing radial and meshed distribution topologies is that, even if the proposed optimization methodology is independent of the grid configuration, the expected annual operating costs are drastically affected by the grid configuration. The final objective function value was about USD/year 108,249.36, whereas, for the meshed topology, it was about USD/year 77,870.17, i.e., a difference of 30,379.19 dollars per year of operation. These differences are clearly explained by the nature of the meshed topol-ogy, where the voltage profiles are improved and a better distribution of the load flows is presented.
The main limitations of the proposed master-slave optimization method are its dependence on the experience of the programmer with the whole algorithm, including the selection of the algorithm parameters, in addition the fact that a statistical analysis is required to determine the average performance of the methodology regarding the objective function values and the processing times. Another important limitation of our proposal corresponds to its extension to power systems: even if the SSA approach can define the location and size of the reactive power compensators, the power flow employed in this research is only applicable to distribution networks without voltage-controlled nodes, which is not the case of power systems, where multiple power generates control voltage profiles in different nodes of the network. This implies the need for replacing the power flow solver with one that can deal with multiple voltage-controlled sources.
As future work, it will be possible to conduct the following studies: (i) combining D-STATCOMS and renewable energy resources to reduce the expected grid operating cost for a planning horizon of 10 or 20 years; (ii) considering the combination of D-STATCOMs and batteries in the same MINLP model and the proposed master-slave optimization methodologies using recently developed metaheuristic optimizers such as the generalized normal distribution optimizer and the arithmetic optimization algorithm; (iii) improving the mathematical formulation used to represent the D-STATCOMs in order to allow for simultaneous dynamic and steady-state analyses of the distribution network; and (iv) implementing new metaheuristic optimizers such as the Farmland Fertility algorithm, the African Vultures Optimization algorithm, and the Artificial Gorilla Troops algorithm to deal with the problem studied in this research.