Enhancement of Hosting Capacity with Soft Open Points and Distribution System Reconﬁguration: Multi-Objective Bilevel Stochastic Optimization

: Soft open points (SOPs) are power electronic devices that replace the normal open points in active distribution systems. They provide resiliency in terms of transferring electrical power between adjacent feeders and delivering the beneﬁts of meshed networks. In this work, a multi-objective bilevel optimization problem is formulated to maximize the hosting capacity (HC) of a real 59-node distribution system in Egypt and an 83-node distribution system in Taiwan, using distribution system reconﬁguration (DSR) and SOP placement. Furthermore, the uncertainty in the load is considered to step on the real beneﬁts of allocating SOPs along with DSR. The obtained results validate the e ﬀ ectiveness of DSR and SOP allocation in maximizing the HC of the studied distribution systems with low cost.


Introduction
Renewable energy integration has been crucial in recent decades to limit the effect of green-house gases on the environment [1][2][3].Various strategies have been adopted in the modern smart grids to increase their capabilities to accommodate the intermittent renewable resources [4][5][6].Hosting capacity (HC) is the mathematical expression that represents the ability of a distribution system to host distributed generation (DG) without violating its operational limits [5][6][7][8].Many methodologies for improving HC [5][6][7][8] have been proposed in the literature, including power quality (PQ) enhancement, network reinforcement, distribution system reconfiguration (DSR), static var compensators (SVCs), energy storage systems (ESS), and soft open points (SOPs), among others.The maximum DG penetration was determined in [9] for the 18-node and 33-node distribution systems, considering the IEEE 519 allowable voltage harmonic limits.A constrained harmonic distortion study was carried out in [10] to maximize HC using harmonic filters.In [11], network reinforcement was employed to maximize the HC of an existing Egyptian distribution network.The deterministic HC was assessed for a real-life grid in Jordan called 73-node distribution system having a high X/R ratio, and a hypothetical grid having a lower X/R ratio called 19-node distribution system designed based on the Jordanian Energies 2020, 13, 5446 2 of 20 standards [12].The HC was successfully maximized in [12] via two strategies, including network reinforcement and reactive power control.DSR was employed in [13] to maximize the HC of the 33-node distribution system using linear load flow formulation.HC was successfully enhanced in [14] via DSR and ESS for a real distribution system in central Italy, including 5 feeders and 287 busbars.A multi-period DSR was employed in [15,16] to enhance the HC of the IEEE 123-node and 1001-node distribution systems.Static and dynamic DSR were investigated in [17], and a static reconfiguration was found to be beneficial at the planning stage.In contrast, dynamic reconfiguration was found to be useful for active distribution networks, especially when a higher number of remotely controlled switches are available.DSR was employed in [18] to maximize HC by selecting the best configurations of the distribution network capable of maximizing the HC of a large distribution network in Japan including 235 switches.A stochastic optimization was developed in [19] for the optimal allocation of SVCs to maximize PV penetration while considering uncertainties related to photovoltaics (PVs) and loads.Quadratic power control for a central battery storage system was proposed in [20] to optimize the penetration of rooftop PVs.In [21], a linearized power flow model was formulated to determine the maximum HC of a 33-node distribution system while considering the load uncertainty.A probabilistic optimization approach was proposed in [22] to assess the effects of uncertainties on HC, while in [23], a strengthened second-order cone programming problem was formulated to maximize the HC of a 33-node distribution system using SOPs to replace tie-lines with fixed locations.An algorithm was proposed in [24] to assess the increase in the HC of a generic distribution system for SOP placement in the UK.HC has also been evaluated for various types of SOP placement, including two-terminal and multi-terminal SOPs and SOPs with energy storage [25].In these previous works, individual strategies were employed to maximize the HC of the distribution systems, but SOPs and DSR were not combined.In the present work, we put forward a novel approach for HC enhancement based on simultaneous SOP allocation and DSR to step on the benefits of the meshed networks in the presence of load uncertainties for two real distribution systems and also ensuring the radial structure while reconfiguring the non-SOP tie-lines which provide resiliency in allocating DGs.The main contributions of this work can be summarized as follows: (1) A multi-objective bilevel optimization problem is formulated to minimize the total active losses by introducing DSR to a lower level problem and then to maximize the HC and minimize the annual total cost for two real distribution systems.(2) A probabilistic HC maximization approach is proposed to illustrate the expected impact of load uncertainties on HC. (3) The proposed optimization approach ensures radiality among the studied distribution networks while reconfiguring the non-SOP tie-lines in a short time.(4) A combination of DSR and SOPs is successfully used to support the penetration of DGs in distribution systems while guaranteeing an economic planning framework.
The rest of this paper is organized as follows: Section 2 presents the problem statement, including the power flow equations, DSR methodology, DG, and SOP models.Section 3 introduces the problem formulation.Section 4 presents the obtained results and a discussion, while Section 5 contains the conclusions and identifies future work in this area.

Materials and Methods
In this section, the DSR algorithm, the power flow equations, DG, and SOP modeling are illustrated in detail.Figure 1 provides an overview of DG and SOP modeling in a distribution system.

Power Flow Equations
The power flow equations required to solve the distribution system under study are formulated as follows [26]: where  = {1,2, … ,  } is the set of nodes of the distribution system and  is the total number of nodes in the distribution system. is the impedance of the bth line joining nodes  and  + 1, and its real and imaginary components are  and  , respectively. is the apparent power injected at the jth node, where  and  are its real and imaginary components, respectively. is the jth node voltage. is the demanded active power at node  + 1.

Distribution System Reconfiguration
The graph theory-based DSR algorithm proposed by the current authors in [26] is employed in this work.This method is based on a graphical interpretation of the nodes (vertices) and lines (edges) of the distribution system.The primary advantage of this DSR method is that it gives the global/nearglobal solution within a short time for large distribution systems.The search procedure used to optimize the current best configuration of the distribution system ( ), hinges on fetching the various possibilities for exchanging the yth tie-line ( ) included in the set of the tie-lines () with its neighboring sectionalized lines, connected to its "From" ( ), and "To" ( ) nodes, as shown in Figure 2.Then, the candidate sectionalized lines required to be tie-lines are sorted in a descending order via an index called weighted voltage deviation index (), where the  for the Eth sectionalized line ( ) connected to the yth tie-line is calculated as follows: such that,  and  are the "From" and the "To" ends of the line  .Further, at this step, a loop starts here till finding a better configuration for the distribution network having the best fitness value, where a temporary test vector ( ) equals to  is initialized.Then, the loop continues by changing the open/close status of each tie-line  and its neighboring sectionalized line in  .

Power Flow Equations
The power flow equations required to solve the distribution system under study are formulated as follows [26]: where B = {1, 2, . . ., N n } is the set of nodes of the distribution system and N n is the total number of nodes in the distribution system.Z b is the impedance of the bth line joining nodes j and j + 1, and its real and imaginary components are r b and x b , respectively.S j is the apparent power injected at the jth node, where P j and Q j are its real and imaginary components, respectively.V j is the jth node voltage.P L j+1 is the demanded active power at node j + 1.

Distribution System Reconfiguration
The graph theory-based DSR algorithm proposed by the current authors in [26] is employed in this work.This method is based on a graphical interpretation of the nodes (vertices) and lines (edges) of the distribution system.The primary advantage of this DSR method is that it gives the global/near-global solution within a short time for large distribution systems.The search procedure used to optimize the current best configuration of the distribution system (R rec ), hinges on fetching the various possibilities for exchanging the yth tie-line (T y ) included in the set of the tie-lines (T) with its neighboring sectionalized lines, connected to its "From" (m T y ), and "To" (n T y ) nodes, as shown in Figure 2.Then, the candidate sectionalized lines required to be tie-lines are sorted in a descending order via an index called weighted voltage deviation index (WVD), where the WVD for the Eth sectionalized line (k E ) connected to the yth tie-line is calculated as follows: such that, m kE and n kE are the "From" and the "To" ends of the line k E .Further, at this step, a loop starts here till finding a better configuration for the distribution network having the best fitness value, where a temporary test vector (R done to ensure radiality while compiling this loop, including (1) if the power flow did not converge, thus, R temp rec is not radial because of the failure in one or more of the connected loads as shown in trail #2 in Figure 2, or (2) if a better objective function value is obtained, R rec is updated to R temp rec and the loop is terminated, or (3) the loop terminates, if no possibility for exchanging the tie-lines with their neighboring lines.This algorithm was tested on various real and large distribution networks up to the 4400-node distribution system, including 185 tie-lines, and has proven its ability to ensure radiality and fast convergence toward finding the best configuration in a short time, since it directly exchanges the existing tie-lines with its neighboring sectionalized lines, thus no randomness exists.In this paper, two case studies are conducted to maximize the HC of real distribution systems using DSR and SOPs placement, including the deterministic and the probabilistic case studies.In both cases, the SOP is placed instead of a certain tie-line, providing resiliency in delivering apparent powers between the SOP terminals.Regulating equations for SOPs are provided below for the deterministic and probabilistic cases.

Deterministic Case
An SOP is allocated if its allocation variable  is equal to one, and is not assigned if  is equal to zero.Each tie-line is connected to either two feeders or loop laterals, the first of which is denoted by  and the other by .Thus, the Ith feeder connected to the yth tie-line is denoted by  .

SOP Equality Constraint [29]:
SOPs are characterized by their ability in transferring the active powers between the adjacent feeders (i.e.,  and  feeders), where the sum of the injected SOP powers to the Ith and Jth feeders equals to zero in case of lossless SOP placement [29], whereas lossy SOP is represented in (7), where the internal active loss of the two VSCs is considered as follows [29,30]: SOP Capacity Limit Constraint [29,30]:

DG Modeling
In this study, the DG power factor is unity, where the power injected by the uth DG in the deterministic case study is formulated as follows: where as for the probabilistic case study at the sth scenario: and its allocation is constrained by the binary variables L DG u and L DG u,s in the deterministic and probabilistic cases, respectively, where L DG u and L DG u,s are equal to one in the case of DG allocation at the uth node.S DG is the maximum size of the installed DGs.S is the set of all scenarios studied.

SOP Modeling
SOPs were introduced in 2011 as a resilient alternative to normally open points (NOPs) and can provide the flexibility of meshed networks in sharing active and reactive power between adjacent feeders in active distribution systems [27].SOPs may have several integration topologies, including the back-to-back voltage source converter (VSC), static series synchronous compensator (SSSC), and unified power flow controller (UPFC) [28].In this work, the back-to-back VSC is used because of its ability to improve operational and power quality indices [29].
In this paper, two case studies are conducted to maximize the HC of real distribution systems using DSR and SOPs placement, including the deterministic and the probabilistic case studies.In both cases, the SOP is placed instead of a certain tie-line, providing resiliency in delivering apparent powers between the SOP terminals.Regulating equations for SOPs are provided below for the deterministic and probabilistic cases.

Deterministic Case
An SOP is allocated if its allocation variable L SOP y is equal to one, and is not assigned if L SOP y is equal to zero.Each tie-line is connected to either two feeders or loop laterals, the first of which is denoted by I and the other by J. Thus, the Ith feeder connected to the yth tie-line is denoted by I y .
SOP Equality Constraint [29]: SOPs are characterized by their ability in transferring the active powers between the adjacent feeders (i.e., I and J feeders), where the sum of the injected SOP powers to the Ith and Jth feeders equals to zero in case of lossless SOP placement [29], whereas lossy SOP is represented in (7), where the internal active loss of the two VSCs is considered as follows [29,30]: SOP Capacity Limit Constraint [29,30]: Each SOP is composed of two VSCs.These VSCs are connected back-to-back through a DC-link capable of transferring both active and reactive powers constantly, and their governing equations are formulated as follows [29,30]: where, P SOP are the active losses by each VSC [29,30].
To limit the allocated SOP size instead of the yth tie-line, the following capacity constraint is formulated as follows [29,30]: where, S SOP max is the maximum SOP size.SOP Internal Power Loss Equations [29]: Energies 2020, 13, 5446 6 of 20 , ∀y ∈ T (13) A SOP L is the loss coefficient of each VSC [29,30].P SOP−loss tot is the total SOP losses, and N SOP is the number of SOPs installed.

Probabilistic Case
In this case, a set of scenarios S = {1, 2, 3, . . ., N s } are applied to the distribution systems, where N s is the total number of scenarios.
SOP Equality Constraint [29]: P SOP I y ,s + P SOP J y ,s + P SOP−loss I y ,s SOP Capacity Limit Constraint [29]: In addition to the capacity constraints demonstrated in (10), and ( 11), the following constraints expressed in (16), and (17) are considered in this case study.

Scenario Reduction
In this work, the load uncertainties are represented by a preset number of scenarios, each of which has a certain corresponding probability.The 8760 hourly data are reduced to a relevant set of scenarios using the backward reduction technique developed by Growe-Kuska et al. in 2003 for stochastic programming [31].

Problem Formulation
In this work, a multi-objective bilevel optimization problem is developed to maximize HC and minimize the total annual cost (C A ) as the upper-level problem, whereas the lower-level problem is dedicated to minimization of the energy loss cost (C Loss ) using DSR.C A consists of the capital cost for SOP installations (C cap ), the annual operational cost of the SOPs (C op ), and the annual total loss cost (C loss ) [30].
where P L tot is the total active demand power for the distribution system under normal loading conditions, c is the electricity price, λ is the interest rate, q is the number of years, c SOP is the SOP capital cost per unit capacity, and α is the SOP annual operation cost coefficient.
The solution procedure of the proposed optimization problem starts by inserting DGs at specific DG nodes (indicated by L DG u ), and SOPs instead of specific candidate tie-lines (indicated by L SOP y ) at the upper level, then the rest of the tie-lines (non-SOP tie-lines) are changed at the lower level by the DSR algorithm to minimize C loss , as a result, the power loss minimization is improved and hence choosing the best tie-lines locations for the next planning stage at the upper level.At the end of the optimization process, the technique for order of preference by similarity to ideal solution (TOPSIS) [32] algorithm takes place to choose the best solution (alternative) among the Pareto solutions, including the HC and the total annual cost quantities by a ratio of 75% and 25%, respectively.

Deterministic HC
For the deterministic case, the total power loss (P Loss ) is formulated as follows: C Loss = 8760.c.P Loss (27) where, |I b | is the magnitude of the current passing through the bth line and N b is the total number of lines.

Upper Level
At this level, a multi-objective optimization problem is formulated to maximize HC and minimize C A .The objective functions f 1 , f 2 are formulated as follows: subject to (1)-( 3), ( 5), ( 7)-( 13), and the following operational limits: where HC is the hosting capacity, V L and V U are the upper and lower voltage limits of the jth node, respectively.I U b is the thermal current limit for the bth line.W is the set of lines.P slack is the active power delivered from the substation.

Lower Level
The objective function ( f 3 ) at this level aims to minimize C Loss by DSR to obtain the best operational point.The objective function is expressed as follows: subject to (1)-( 4), ( 26)- (27), and ( 29)-(31).

Probabilistic HC
In this case, the optimization problem is solved for several different scenarios S. The expected total annual energy loss (P exp Loss ) and the probabilistic HC (PHC) are expressed as follows: C Loss = 8760.c.P exp Loss (34) where I b,s is the magnitude of the current in the sth scenario.r b,s is the bth line resistance at the sth scenario.P SOP−loss tot,s is the total SOPs losses in the sth scenario.p s is the probability of the sth scenario.The bilevel multi-objective optimization for this case is formulated as follows:

Upper Level
At this level, the objective function f 4 = PHC is maximized and f 5 = C A is minimized subject to (1)-( 3), ( 6), ( 10)-( 11), ( 15)-( 19), ( 33)- (35), and the following constraints: where, V j,s is the magnitude of the jth node, at the sth scenario.P slack s is the injected power to the slack node, at the sth scenario.
The interaction between the upper and the lower optimization problems during the solution of the current case study is illustrated as follows: Step 1: Initialize the multi-objective optimization problem at the upper level, including the number of populations, number of iterations, also, the number of decision variables is given and their upper and lower limits.
Step 2: Input the optimization parameters, including R rec , L DG u,s , P DG u,s , L SOP y , S SOP I y , S SOP J y , P SOP I y ,s , Q SOP I y ,s , P SOP J y ,s , and Q SOP J y ,s .Step 3: Set SOPs locations according to L SOP y .
Step 4: For each scenario.

Results and Discussion
In this paper, a real 59-node distribution system in Cairo [33] and the 83-node distribution system [26] are used as test systems, shown in Figures 3 and 4, respectively.The input data configurations used for the upcoming case studies are provided in Table 1.In the subsequent case studies, three multi-objective optimization techniques are employed to solve the upper-level optimization problem: the non-dominated sorting genetic algorithm (NSGA-II) [34], multi-objective particle swarm optimization (MOPSO) [35], and multi-objective multi-verse optimization (MOMVO) [36].The authors in [34] proposed a new methodology to solve mixed-integer nonlinear programming problems via integrating the surrogate and the deterministic infeasibility sorting, where it has been tested against real-life building applications.A modified version of the MOPSO was proposed in [35] to handle multi-objective MINLP problems, where it has been tested against several benchmark functions and has proven its ability to find the best solution.In [36]      For the 59-node distribution system, the HC reaches >98%.Table 2 provides detailed results for these three optimization techniques.Tie-lines, SOPs size, and locations are given in Table 3.The installed DG nodes and their sizes are shown in Table 4 for MOMVO.For the 83-node distribution system, the HC reaches >98%.Table 5 gives detailed results for these three optimization techniques.It is evident that MOPSO outperformed against NSGA-II, and MOMVO by providing the best HC, and the lowest C A .Tie-lines, SOPs size, and locations are given in Table 6.The installed DG nodes and ratings are shown in Table 7 for MOPSO.Figure 5 shows the voltage profile improvement after allocating DGs and SOPs.

Probabilistic Case Study
For the probabilistic case study, the uncertainty in the load over the year is considered, as shown in Figure 6.The most crucial load scenarios are generated using the scenario reduction algorithm in [31].These scenarios are given in Table 8, including their loading levels and associated probabilities.

Probabilistic Case Study
For the probabilistic case study, the uncertainty in the load over the year is considered, as shown in Figure 6.The most crucial load scenarios are generated using the scenario reduction algorithm in [31].These scenarios are given in Table 8, including their loading levels and associated probabilities.For the 59-node distribution system, the HC reaches >62%, which is lower than the deterministic case study.Table 9 gives detailed results for these three optimization techniques.Tie-lines, SOPs size, and locations are given in Table 10.HCs and total power losses for the studied scenarios are shown in Figures 7 and 8, respectively, and the installed DG ratings for each node are shown in Figure 9.The voltage profiles before and after allocation of DGs and SOPs for the different scenarios are shown in Figure 10.

Real 59-Node Distribution System in Cairo
For the 59-node distribution system, the HC reaches >62%, which is lower than the deterministic case study.Table 9 gives detailed results for these three optimization techniques.Tie-lines, SOPs size, and locations are given in Table 10.HCs and total power losses for the studied scenarios are shown in Figures 7 and 8, respectively, and the installed DG ratings for each node are shown in Figure 9.The voltage profiles before and after allocation of DGs and SOPs for the different scenarios are shown in Figure 10.

83-Node Distribution System
For the 83-node distribution system, the HC reaches >57%.Table 11 gives detailed results for these three optimization techniques.Tie-lines, SOPs size, and locations are given in Table 12.HCs and total power losses for the studied scenarios are shown in Figures 11 and 12, respectively, and the installed DG ratings for each node are shown in Figure 13.The voltage profiles are shown in Figure 14 for the different scenarios before and after the allocation of DGs and SOPs.Also, the 83-node distribution network tie-lines, DGs and SOPs locations using MOPSO are shown in Figure 15.

83-Node Distribution System
For the 83-node distribution system, the HC reaches >57%.Table 11 gives detailed results for these three optimization techniques.Tie-lines, SOPs size, and locations are given in Table 12.HCs and total power losses for the studied scenarios are shown in Figures 11 and 12, respectively, and the installed DG ratings for each node are shown in Figure 13.The voltage profiles are shown in Figure 14 for the different scenarios before and after the allocation of DGs and SOPs.Also, the 83-node distribution network tie-lines, DGs and SOPs locations using MOPSO are shown in Figure 15.maximizing the expected HC for the 59-and the 83-node distribution systems, and minimizing the expected total annual cost better than NSGA-II and MOMVO for the 83-node distribution system.Thus, in most of the study cases, MOPSO has outperformed against NSGA-II and MOPSO.To sum up, SOPs insertion, along with the reconfiguration of the non-SOP tie-lines has provided a new perspective in promoting the concept of meshed networks in the distribution grids from the planning and the economic viewpoints.Also, the proposed planning strategy was applied to real distribution feeders to step on its validity before using it on real large distribution feeders like the 415-, 880-, 1760-, and 4400-node distribution systems, which will be considered in future works.

Conclusions
In this work, a bilevel multi-objective optimization approach is proposed to maximize the HCs of two real distribution systems using simultaneous DSR and SOP allocation.Two case studies are conducted using both systems, including both deterministic and probabilistic optimization approaches.From the obtained results, it is clear that HC was maximized efficiently while considering the load uncertainties, and the system is expected to accommodate greater than 62% and 57% DG penetration in the 59-and 83-node distribution systems.The use of DSR with SOPs As shown in this probabilistic study, it is clear that MOPSO outperformed against NSGA-II and MOPSO in maximizing the PHC. of the studied distribution networks up to >62%, and >60 for the 59-and the 83-node distribution systems, respectively.Thus, it is clear that the proposed DGs planning strategy has been successfully assisted via simultaneous DSR and SOPs allocation, where the 59-node distribution system accepted allocation of DGs among the distribution system nodes, and the 83-node distribution system accepted allocation of DGs among 68 nodes of the distribution network.Besides, the radiality of the non-SOP feeders is ensured, as shown in Figure 15.Also, MOPSO succeeded in maximizing the expected HC for the 59-and the 83-node distribution systems, and minimizing the expected total annual cost better than NSGA-II and MOMVO for the 83-node distribution system.Thus, in most of the study cases, MOPSO has outperformed against NSGA-II and MOPSO.To sum up, SOPs insertion, along with the reconfiguration of the non-SOP tie-lines has provided a new perspective in promoting the concept of meshed networks in the distribution grids from the planning and the economic viewpoints.Also, the proposed planning strategy was applied to real distribution feeders to step on its validity before using it on real large distribution feeders like the 415-, 880-, 1760-, and 4400-node distribution systems, which will be considered in future works.

Conclusions
In this work, a bilevel multi-objective optimization approach is proposed to maximize the HCs of two real distribution systems using simultaneous DSR and SOP allocation.Two case studies are conducted using both systems, including both deterministic and probabilistic optimization approaches.From the obtained results, it is clear that HC was maximized efficiently while considering the load uncertainties, and the system is expected to accommodate greater than 62% and 57% DG penetration in the 59-and 83-node distribution systems.The use of DSR with SOPs succeeded in improving the indices of the systems and also in minimizing the expected total annual costs by more than 75.5126% and 57.5967% compared to the initial annual cost for the 59-and 83-node distribution systems, respectively, at the maximum DGs penetration while considering the load uncertainties.It is therefore clear that allocating SOPs with the use of DSR succeeded in improving both the HC and the system
temp rec ) equals to R rec is initialized.Then, the loop continues by changing the open/close status of each tie-line T y and its neighboring sectionalized line in R temp rec .Various tests are Energies 2020, 13, 5446 4 of 20

Figure 2 .
Figure 2. Possible trails to exchange the tie-line  with its neighboring lines, including  ,  , and  .

Figure 2 .
Figure 2. Possible trails to exchange the tie-line T y with its neighboring lines, including K 1 , K 2 , and K 3 .

i
are the sending and receiving transferred complex apparent powers by the SOP installed instead of the tie-line y, respectively.P SOP−loss I y and P SOP−loss J y Set DGs locations according to L DG u,s .(c) If (1)-(

Figure 5 .
Figure 5. Voltage profile for the 83-node system in the deterministic case.

Figure 5 .
Figure 5. Voltage profile for the 83-node system in the deterministic case.

Figure 6 .
Figure 6.Load profile over one year.

Figure 10 .
Figure 10.Voltage profiles for the 59-node system in different scenarios.

Figure 10 .
Figure 10.Voltage profiles for the 59-node system in different scenarios.

Figure 14 .
Figure 14.Voltage profile for the 83-node system for different scenarios.

Figure 15 .
Figure 15.83-node distribution network after SOPs and DGs allocation.

Figure 15 .
Figure 15.83-node distribution network after SOPs and DGs allocation.
)-(38) violated, then return C A , and PHC equal to infinity, and zero, respectively.Then, the optimization process continues by setting new parameters' Loss , C cap , C op , C A , and PHC.After that, if the number of iterations is reached, jump to Step 6; otherwise, the upper-level optimization sets new parameters' values and continues from Step 2.

Table 2 .
Deterministic results for 59-node distribution system.

Table 5 .
Deterministic hosting capacity (HC) results for 83-node distribution system.

Table 9 .
Probabilistic results for the 59-node distribution system.

Table 11 .
Probabilistic results for the 83-node distribution system.

Table 11 .
Probabilistic results for the 83-node distribution system.