The Use of Black-Box Optimization Method for Determination of the Bus Connection Capacity in Electric Power Grid

One of the main tasks of the Power System Operators is to ensure a proper, safe, and trouble-free operation and development of power grids. Growth of power system is inseparably linked to a connection of both renewable and classical new energy sources. For network operators and potential investors, it is essential to know the place and volume of generation capacity that can be connected to the grid. A publication of such data is currently a legal obligation for many operators. This paper proposes a method of determining the bus connection capacity in power grid of any type with the use of black-box optimization. Calculations and analyses were performed with a full, nonlinear model of the analyzed network. Obtained results show the effectiveness of this method for both single and multiple nodes in any configuration. All the constraints relevant for the proper and safe system operation, such as bus voltages, line loads, and short-circuit currents, both in a steady-state and (n-1) contingency states, are taken into consideration. Calculations confirmed the good convergence and repeatability of the method for all three tested computational algorithms. This has also confirmed the possibility of use of open source software to extend the functionality of Siemens PSS®E commercial power system calculation software.


Introduction
One of the most important factors taken into account in the network development forecast is knowledge of the network bus connection capacity. This information is essential to any Power System Operator (PSO) and a potential investor in order to identify the development and investment opportunities in new power sources.
The bus connection capacity can be defined as the maximum active power of generating units which can be connected to a single node or to a selected group of nodes without deteriorating the existing status of the power grid. Knowledge of these data allows PSO to determine the development opportunities in the network and to identify weak elements requiring an upgrading.
The literature lacks descriptions of methods allowing PSOs to determine, in a simple and quick way, the bus connection capacity for the selected nodes of the operated network, without the necessity of laborious creating of auxiliary computational models. The problem of determining the bus connection capacity is resolved indirectly in many publications but is not the main objective of the articles. In most cases authors use optimization methods with calculations based on deterministic models. Network models used in calculations have to be, however, each time created anew, mostly by transformation or reduction of existing power network models used in power flow calculations.
Most of the papers that have been published so far, however, has been focused on the aspects that are not very useful for PSOs in the day-to-day operation of the power grid. Some authors focus on economic aspects of capacities to be connected, useful for short-and long-term planning.
Afful-Dadzie et al. [1] present a multi-period stochastic optimization model for studying long-term power generation capacity planning in developing countries. The problem is modeled using a stochastic mixed-integer linear programming technique under demand uncertainty. Wogrin et al. [2] apply multi-year bi-level equilibrium model to analyse the generation capacity expansion problem in a liberalized framework. The model is based on a market equilibrium equation with constraints, which are transformed into a mixed-integer linear program and solved using diagonalization in order to maintain the market equilibrium condition. Hesamzadeh and Yazdani [3] propose a mathematical model for transmission planning in an environment where there is imperfect competition in the electricity supply industry. The model is developed based on the concept of the leader-followers game in applied mathematics. Tohidi et al. [4] apply mathematical models for progressive transmission expansion planning including new generation resources. The paper study both the proactive coordination modeled as a mixed-integer bi-level linear program (MIBLP) and the reactive coordination, which is modeled as a mixed-integer linear program (MILP).
The authors use simplified linear system models in the above articles. These models, although useful for forecasting and making strategic decisions, do not fully reflect all important phenomena in actual power systems.
Other papers analyze bus connection capacity referring to generation technology. These methods provide an opportunity to identify and select technologies that maximize profits of potential investors.
Centeno et al. [5] deal with the case of generation company, which need to to choose between a number of alternative technologies. The main outcome of the paper consist in the two stage market simulation model using a conjectured-price-variation approach. In the first step the net present value is maximized assuming that the company is in an incumbent position and knows the new capacity that will be built by the others. Pisciella et al. [6] present a multi-stage stochastic optimization model for the generation capacity expansion problem of a price-taker power producer. The objective function represents a trade-off between expected profit and risk in the presented model. Chen et al. [7] present a copula-based fuzzy chance-constrained programming model and applied it to electric power generation systems planning under multiple uncertainties.
Few publications focus on economic issues related to determining the bus connection capacity in the power grid while providing an optimal mix of connected sources. Mokgonyana et al. [8] propose a planning approach, which maximizes profits of distribution companies. The model is used to define both optimal location and capacity for renewable energy sources including both self-producers and grid exporters.
This approach is characteristic for planning and defining of the energy mix, it is not very useful, however, from a practical point of view of PSOs, who are obliged to show the connection capacities regardless of the power source technology.
Determining the bus connection capacity in relation to accommodation of renewable and distributed energy sources in existing transmission and distribution networks is gaining weight in recent years, due to the high degree of penetration of such sources.
Theo et al. [9] present a review of distributed generation system planning and optimization techniques. Numerical and mathematical modeling optimization methods are reviewed, analyzed and criticized with recommendations for their improvements. Mohtashami et al. [10] focus on a multi-year planning of distribution networks with large share of distributed, renewable sources. The optimization model considers both investment in traditional network technologies and smart grid solutions including dynamic line rating, active network management, quadrature-booster transformers. This optimization considers also adjustment of the set points for various network devices, the curtailment of distributed generation output taking into account take into account different network access options (firm or non-firm contracts).
Due to the specifics of the location of renewable sources in distributed generation networks, methods used in these cases are limited to radial networks.
Koutroumpezis and Safigianni [11] propose a method to determine the optimum allocation of the maximum distributed generation penetration in medium voltage power distribution networks. Technical constraints, such as thermal rating, transformer capacity, voltage profile, and short-circuit level, are considered. Papaioannou and Purvins [12] present an easily applicable methodology to calculate maximum distributed generation capacity in radial low-voltage feeders. Kolenc et al. [13] present a probabilistic approach to network planning, which has many advantages compared to the traditional approaches using estimated peak values and empirically defined simultaneity factors. The method takes into account the stochastic natures of future distributed generation location and loads consumption.
Xu and Zhuan [14] are searching for the optimal wind generation capacity considering grid coupling assets. The authors maximize yearly profits of the entire power system accounting for operational constraints and the reliability requirement. The problem is solved using combined Monte Carlo simulation, neural networks and genetic algorithm optimization technique.
Only a few papers deal with the determination of bus connection capacity taking into account only technical limitations, separate from economic issues. Sobierajski et al. [15] bring down the task of determination of bus connection capacity to the minimization of a linear goal function subjected to linear equality (generation and demand balance) and inequality constraints (permissible branch flows and thermal generation minimum). A linear programming method has been used, however, to solve so formulated problem. Bajor [16] describes the use of coherent node method for bus capacity determination, but this method can be applied only for a group of nodes selected in a very specific way.
Harrison and Wallace [17] present a method for determining the installed capacity of new generation sources in the distribution network. In the method, named by the authors as "reverse load-ability", new generation sources are modeled as loads of negative value. Using the optimal power flow method, the authors search for the minimum value of negative load. In the calculation process, technical constraints are considered, such as: maximum capacity of distribution lines, voltage tolerance boundary at network nodes. Cornélusse et al. [18] also propose a method for determining the connection capacity of new generation sources to the distribution network. Additionally, in this method the authors use typical tools that are usually used by Network Operators in the planning and operation of distribution systems: long-term load forecasting, long-term network expansion/reconstruction planning, network reconfiguration, and power flow calculation tools. In both papers, the authors analyze the problem of, "network sterilization". This problem occurs in the case of improper distribution of the connected generation between network nodes, e.g., too high generation capacity connected to one of the network nodes. Consequently, it may block the possibility of connecting additional power to other nodes of this network. Both methods use a power flow model of the grid offering the exact mathematical description of the power grid operation parameters. In both methods, the calculation model must be created from scratch (based on network data).
Georgilakis and Hatziargyriou [19] in their article attempted to summarize the methods used to determine the connection capacity of new generation sources to distribution networks. It presents an overview of the most advanced methods and models used to solve the problem. A key contribution of this paper consists in the analysis and clustering of research trends in this field.

Proposed Method and Its Novelty
One of the main tasks of the PSOs is to ensure a proper, safe and trouble-free operation of power grids and their successive and economically justified development. One of the key factors in the successful accomplishment of these tasks is profound knowledge of the bus connection capacity. For this purpose optimization methods are used, which allow PSOs to more and more efficiently and rationally utilize technical capabilities of the operated networks. This problem can be generally solved in two ways:

•
Developing an own, fully functional power system calculation software with access to all variables necessary in the optimization process. It is a labor-intensive method requiring a long time and multiple verifications.

•
Using available open source or freeware power flow calculation software, which allows for performing power flow calculations with simultaneous access to all variables necessary in the optimization process. Unfortunately, such software most often does not meet all requirements for professional power flow calculations. For example, an inability to adjust a transformer ratio or lack of reactive power regulation significantly limits the application of such programs in actual power networks [20].
Optimization of operation or structure of a power system is a known and well described problem in literature [21][22][23][24]. Various methods and optimization algorithms are used to solve this problem. The main purpose aims at minimizing the costs of system operation or at the improvement of technical parameters of system operation, such as active or passive power losses or a system stability reserve. Currently, primarily commercial software enabling high-speed computing, efficiency, and reliability of the results is used for calculations and analyses of complex engineering systems. Many programs also have optimization modules. The basic disadvantage of these modules is that the user' possibilities are limited to the use of one of the optimization models implemented in the program, without an opportunity to formulate their own optimization problem. The main obstacle is the lack of direct and free access to all system variables and to the full mathematical model of the analyzed power system. As an output from such software, the user receives a set of values describing the state of the system, such as line and transformer loads and bus voltages. For this reason, it is not possible to directly use any power flow software known to the authors to determine the bus connection capacity in the analyzed power grid. But these results can be used to solve many optimization problems related to power systems by use of Simulation Optimization methods [25] or a black-box optimization [26,27].
This paper proposes a new method, which uses black-box optimization to determine the bus connection capacity in electric power network of any type. One of the main elements determining the uniqueness of the proposed method is the use, in calculations of a power system flow, of the full nonlinear power system model, without any simplification or linearization. There is no need to create a new network model or to convert an existing one, PSO can apply a model which is used in daily routine power flow calculations.
The full, nonlinear power system model is applied in all calculations used in the method presented in this article. The necessary condition is the possibility of obtaining numerical results from the flow calculation program. All the constraints relevant for proper and safe system operation, such as bus voltages, line loads, and short-circuit currents, both in a steady-state and in (n-1) contingency states, can be taken into consideration. The type and scope of included constraints can be arbitrary determined by the user. This essential quality ensures that the obtained solution meets the constraints of the actual power system.
In addition, the proposed algorithm is an example of the so-called technical optimization, which determines only the values of the generation power that can be connected in the analyzed nodes regardless of the generation cost (does not determine an optimal "mix" of generating sources). This approach is preferred by all PSOs, who primarily wish to determine how large a source can be connected to the node(s), without worsening the existing network state or define a clear indication of network investments necessary to connect the source of a certain size in the node(s). The problems of cost and choice of generation technologies are usually on the site of the investor, who in turn should not be burdened with technical problems regarding the network.
Obtained results show the effectiveness of this method both for single nodes and for their groups in any configuration. Calculations confirmed good convergence and repeatability of the method for all three tested computational algorithms. The method creates the possibility to use open source software to extend the functionality of Siemens PSS R E commercial power system calculation software.
Authors of the article do not know any publications regarding methods that directly determine bus connection capacity with the use of power flow calculation software employing full, nonlinear network model. This is especially important for PSOs who, on the basis of such models, simulate and analyze power system operation on a daily basis. This article fills this gap, giving simultaneously the opportunity to simplify and expand the functionality of already used software.

Black-Box Optimization
An optimization algorithm used in the black-box method is usually based on a simple search algorithm [26,28]. This is mainly due to the difficulty or inability to calculate gradients or second derivatives matrices used by more efficient gradient or Newton methods. The problem of determination of the bus connection capacity in the power grid is nonlinear and, in its general form, can be presented as below: Adaptation of the above equations for the black-box optimization method is as follows: for each iteration t = 1, . . . , t max , the optimization algorithm calculates the current value of searched variables vector X t and passes it to the black-box object. The black-box in the presented problem constitutes a model of the power system under consideration, implemented in the power system calculation software. The software calculates the power flow for given vector X t and passes numerical values of the objective function ( f (X t )), all inequality constraints g j (X t ) and/or all equality constraints h k (X t ) to the optimization algorithm. Basing on these values, the optimization algorithm calculates the new value of the searched variables vector X t+1 , which is subsequently used in the next iteration of the optimization process. The knowledge of functional dependencies between the state vector X and the constraints is not necessary and in practice is often lacking. The result of this method is the optimal variables vector X opt , satisfying defined constraints.
The searched variables of the optimization problem are capacities (P Add i ) to be connected at selected power system nodes and the objective function is the maximum sum of these capacities. The most important values determining the proper and safe operation of the power system (such as: admissible transformer and line loads (I r j ), admissible bus voltages (U min j , U max j ), allowable short-circuit current values (I SC j ) in network nodes, etc.) are the problem constraints. The user can also create his own set of constraints, taking into account the specific requirements of the PSO. A characteristic feature of the presented method is that the above-presented constraints are not explicitly described by mathematical equations. The constraints values are obtained from the power flow calculation software (a black-box in the method) in a numeric form. The optimal solution is found by the optimization algorithm to which constraints values are passed (from a black-box) in each step of the calculations. A detailed model of the optimization problem is presented in Section 4.
The problem of determining the bus connection capacity in power grid is a nonlinear multi-dimensional and in general case a non-convex optimization problem. Therefore and due to applied optimization algorithms, the solution X opt = P Add opt obtained using the proposed method is only a local optimum, not a global one. In practice, this may result in finding many solutions (local optimums), depending on, e.g., the selection of the starting point or chosen optimization algorithms.
The main fields of black-box optimization application are research areas based on experiments, simulations or real values measurements. Researchers in many cases have a problem with the creation of appropriately accurate models that reflect the phenomena they are investigating or even fail to build these models at all, because of the complexity of the studied phenomena. Black-box optimization methods allow for efficient and accurate calculations in such cases, despite lack or incomplete knowledge of the exact form of the objective function or constraints. Many examples of the effectiveness of black-box optimization in solving complex engineering problems can be found in the literature.
Schaul [29] presents and compares the latest methods in the field of continuous black-box optimization. He also shows possibilities and implementations of these methods in an open source Python environment. Bajaj and Hasan [30] propose an effective strategy for the sampling, modeling, and optimization of black-box problems with constraints. Santarelli and Pellegrino [31] use black-box method to optimize a stand-alone energy system based only on renewable sources integrated with a system for the production of hydrogen, with the aim to supply the electricity needs of a residential user. Brus and Zambrano [32] adapt a black-box method for nonlinear systems for identification of the solar collectors at an air conditioning plant. Brabec et al. [33] present black-box models for forecasting hourly average solar irradiance. The results are compared with classical statistical methods. Paulescu et al. [34] compare two advanced models for forecasting the output power of photovoltaic plants: a black-box fuzzy model and a physically inspired, semiparametric statistical model (gray box). Ruschenburg et al. [35] describe the applicability and the ways to validate a black-box heat pump simulation model. The results are backed by research for five real objects. Afram and Janabi-Sharifi [36] show the possibilities of the use of black-box models for effective control of the residential HVAC systems. Tang et al. [37] solve the problem of distribution system feeder reconfiguration using black box Mesh Adaptive Direct Search (MADS). Black and white box method is used by Ma et al. [38] for diagnostic of hub permanent-magnet synchronous motors based on abnormal noise analysis. Bent et al. [39] apply black-box approach based on a discrepancy-bounded local search (DBLS) to optimize expansion of transmission assets and considering power flow related constraints. Valdivia et al. [40] apply a black-box models for the control of DC-DC converters used for fuel-cell power conditioning. A black-box approach for the identification of parameters of proton exchange membranes in fuel cell stack predictive control has been proposed by Lopes et al. [41]. It uses neural networks with nonlinear autoregressive, exogenous input and nonlinear output error.

Mathematical Model
The basic aim of the method presented in this paper is to find the total, maximum active power of generating units that can be attached to any selected group of power grid nodes while maintaining safe system operation. The searched variables are active powers P Add i , which can be connected to an arbitrarily chosen group of N network nodes: ( The objective function for the described problem is: In practice, each newly connected energy source changes the volume of total active power losses in the network, to which it is connected. For the power grid operation, the real gain is the net power, which is a difference between the sum of the active power of the added sources (3) and the change of total active power losses caused by their connection. This can be presented as: Equation (4) is an objective function of the considered optimization task. This may lead to a situation in which the net power is greater than the sum of capacities of connected generating units. This is due to such an allocation of newly connected generating units that results in lower total active power losses in the network. For the complete definition of the optimization task, it is necessary to formulate problem constraints-mathematical relationships describing the technical limitations existing in the actual network. They are: • permissible line and transformer loads: • range of permissible voltage levels at power network nodes: • range of permissible short-circuit current values at network nodes: • limit of the total added generation: All the above constraints should be met both in a steady state and in contingency states ((n − 1) criteria). Usually, a power network fails to meet one or more (n − 1) criteria even without any additional power sources. Then some of the constraints (5) and/or (6) are not fulfilled at the beginning of calculation. In such a case, the constraints (5) and/or (6) must be modified. They take a more general form: and Equations (9) and (10) mean that the connection of an additional generation cannot worsen the existing overloads of lines or transformers and/or voltage transgressions. Results of calculations carried out so far, indicate, that connection of additional generation in many cases lowers or even eliminates the previously described exceedances in contingency states.
Voltage levels specified in the relevant regulations, regarding the correct operation of the power system are used as maximum and minimum bus voltage levels U max j and U min j in constraints (6) and (10). Constraint (7) results from the necessity of maintaining the short-circuit strength of existing in the analyzed network switching, measuring and protection equipment, in face of changing short-circuit conditions, resulting from the connection of new generation sources.
Constraint (8) arises from the requirement to keep the power balance in analyzed power system unchanged; the total output of the designated existing centrally controlled power units has to be reduced by the total power of new connected units. Centrally controlled units are selected in a manner determined by the PSO. In each calculation step their output is reduced in proportion to currently generated power. In practice, due to technical constraints or market mechanisms, centrally controlled power units should remain operating. It means that their output should be within the limits P Gmin j , P G j , giving adjustable power value P Greg j = P G j − P Gmin j . This will ensure, among other things, the minimum requirements for voltage regulation and stability of the entire power system.
The range of reactive power generation of the newly connected power units can be determined arbitrarily as Q min , Q max , and in practice should be related to the unit characteristics Q = f (P) or be determined by PSO's requirements.
Verification of constraints (5) to (10) and power balance corrections are carried out in each step of the proposed method, both for the steady state and for all considered contingency states.
Calculations are performed successively for the steady state n0 and for all k considered (n − 1) contingency states: (n − 1) 1 , . . . , (n − 1) k . From obtained k + 1 solutions P Add n0 , P Add (n−1)1 , . . . , P Add (n−1)k , the one with the lowest objective value is the final solution of the problem. Because of the nonlinearity of the problem, the final solution is additionally verified for compliance with all the constraints in the remaining states of the system.

Validation and Testing
The proposed method is presented on the IEEE 118-bus test network [42,43], shown in Appendix B in Figure A1. The full, nonlinear model of the power system was used.
The largest generating units with output exceeding 300 MW were selected as the centrally controlled units (used to balance the power in the network). Their details are presented in Appendix B in Table A1. A group of nodes was arbitrarily chosen within the considered network, for which the total maximum power of the connected generating sources was to be determined. The group consisted of 15 nodes, and its data are given in Appendix B in Table A2. Calculations were made on a PC characterized by the following parameters: INTEL R Core i7-4790K 4 GHz processor, 32 GB RAM, running on Windows 10. The following software was used in the calculations:

1.
Three optimization algorithms:  Figure 1 shows the block diagram of the presented calculation method. The left part of the algorithm is implemented in the MATLAB R environment and is responsible for optimization. The current calculation point X t = P Add t is computed in this block after each iteration, the optimization algorithms are called from this block and the criterion for calculation termination is tested. The right part of the algorithm is implemented in the Python environment and is responsible for data transfer from the tested black-box object, which in the presented case is the network flow calculation software (Siemens PSS R E) with implemented full nonlinear model of the analyzed power network. Ultimately, MATLAB R can be replaced by Python what will allow eliminating the cost of MATLAB R license. All three optimization algorithms were run with default settings and two stopping criteria: • obtained calculation accuracy ε = 10 −6 , or • a maximum number of iterations (depending on algorithm). The following values have been here adopted: COBYLA-2500; NOMAD-5000; Pattern Search (MATLAB R )-5000. In practical calculations, further increase of the maximum iteration number did not change the accuracy of the obtained results. Calculations were carried out for following network configurations: • steady state, • (n − 1) contingency states, including the single outage of all lines (179 cases) and transformers (9 cases). This together makes 188 different configurations of the tested power grid.
From the solutions obtained for all such defined network configurations, the one with the lowest sum of additionally connected powers described by (4), was taken as the final solution.
Calculations in contingency states were made for two levels of acceptable overload of power system components (lines and transformers): • no overload accepted (max load = 100% rated value), • 20% overload in contingency states accepted (max load = 120% rated value).
The second case results from a common network operation practice, in which small, short-term overloads are allowed. All calculations were made: • separately for each node belonging to a chosen group, • for the group of nodes as a whole.
The results of calculations made with the above assumptions are presented in Tables 2-4. Calculation time was below 1 min for the single nodes and from several min up to a few h for the group of nodes. For example, the following algorithms performance indicators in one of the calculation variants are presented in Table 1. Larger values (both: function evaluation and CPU time) for NOMAD and Pattern Search algorithms result from the specifics of these algorithms that require a larger number of function calculations. Table 1. Example of algorithm performance indicators-variant "Group as a whole", the max load of the lines, and transformers in contingency states = 100%. COBYLA = Constrained Optimization BY Linear Approximations. General comments on calculation results:

Number of Function Evaluation
• All three tested optimization algorithms have confirmed their suitability for solving of the presented problem.

•
Calculations show that the net power resulting from the additional connection of power units sometimes is greater than the sum of power values of the connected generators. This is due to such an allocation of newly connected units that results in lower total active power losses in the network.

•
The results of total bus connection capacity are very close for all three tested algorithms, but the power distribution between the nodes may be different. This is because several different power distributions can give practically the same results. A simple example of such a situation is presented in Figure 2. Both bus capacity distributions in Figure 2 are equivalent. In each of buses 1 and 2, the connection capacity is limited by the rated load of identical lines 1 and 2 (P r 1 = P r 2 =150 MW), while the total connection capacity is bounded by the rated load of line 3 The efficiency of the algorithm depends on the starting point selection. The standard starting point for calculations is point P Add 0 = 0. Such a choice allows avoiding any violation of the constraints (e.g., overloading of network elements) just after the beginning of the calculations. Starting point P Add 0 = 0 also meets the condition of starting optimization calculations from the permitted area, which is required for most of the optimization algorithms incorporating constraints. When an approximate upper power range for considered network nodes is known to a PSO (what often happens in practice), sometimes better results can be acquired by starting calculations from the point P Add

Summary and Conclusions
Conducted calculations and analyses have demonstrated the possibility of using the presented black-box optimization method to determine the bus connection capability in the power grid of any type and for any configuration of nodes. The presented results allow for drawing the following conclusions: • Analysis and calculation of the bus connection capacity can be made using a full, nonlinear grid model, without any simplification or linearization. A very important factor is that the following network flow analyses are performed on the same model.

•
The proposed method is convergent and relatively quick. The number of iterations and calculation time depends on the algorithm used, network size and a number of selected nodes. Sample graphs showing the convergence process for the IEEE 118-bus test network (steady state) are presented in Figure 3.  • The most effective in the conducted analysis was the COBYLA algorithm, which is also the simplest one. The optimum was reached with the smallest number of iterations and in the shortest CPU time. Larger values (both iterations and CPU time) for NOMAD and Pattern Search algorithms result from the specifics of these algorithms that require a larger number of function calculations. Additionally, the effectiveness of the Pattern Search algorithm depends in a significant way on the proper selection of the starting point and the control parameters of the algorithm. In the case of bad selection of control parameters or/and starting point, calculations result was often much worse than for other tested algorithms.

•
Convergence of the proposed method results from convergence of used optimization algorithm, provided the convergence of the power network model.

•
The method can be used for any configuration of buses. PSO can define arbitrarily the composition of such node groups, according to its guidelines or needs.

•
The proposed method makes it possible to take into account all the constraints relevant for safe operation of the system, such as bus voltages, line loads, and short-circuits currents, both in steady state and in contingency states.

•
If the PSO has flow calculation software (which is very common), the proposed method can be implemented at relatively low cost. It is based on open source optimization algorithms, and requires only access to results from flow calculation software. MATLAB R software used by the authors as a computational interface can be fully replaced with procedures written in other programming language, e.g., Python.
The presented method also allows to solve the problem of "network sterilization" described in [17,18]. Correct analysis of the results presented in Tables 2-4 allows to determine such parameters as: • "strongest" node of the analyzed area of the network to which the highest power can be connected (LINKOLN node); • maximum power that can be connected to the entire considered network area (Group as a whole = 67.531 MW); • maximum powers in individual nodes of the analyzed grid area Tables 2 and 3. However, it should be taken into consideration that the installed capacity connected to individual nodes have to respect the following relationship: Group as a whole max ; • the method allows distributing the installed capacity of newly connected generating sources among grid nodes in the considered area by formulating appropriate constraints linked to the investment strategy of the grid operator. This issue was not presented and analyzed in the article.
The proposed method can be successfully applied and used by all PSOs, who own power flow calculation software programs on condition that this software can export calculation results to a file or pass them to other software as a memory variable. Additional benefits for PSOs, software, or designing companies are: • the method can be used for identification of network bottlenecks that should be upgraded to increase the grid connection capacity, • the proposed method can be also implemented by developers of flow calculation software to increase the functionality of the offered tool, • due to the use of the open source software, costs of implementing the method can be relatively low.
For the more practical testing of the proposed method, calculations were also made for full, nonlinear model of the whole Polish Power System, the model of which consists of over 4500 nodes.
The results obtained for the selected case study presented in [52] fully confirmed the suitability of the proposed method and its superiority over algorithm used by PSOs in Poland.

Future Research
The proposed method, thanks to its open structure, can be further developed and refined, without interfering in complicated flow calculation software. The most important proposed further development cover: • testing new optimization algorithms that could increase the speed and accuracy of the method, • identification of additional constraints that will contribute to a more unambiguous power distribution between the nodes, • replacement of the commercial software (MATLAB R ) by routines created in Python programming language.
All future work will focus on developing a fully independent, free computing tool that could be integrated with any power flow calculation software.   Figure A1. IEEE 118-bus test system (PSS R E diagram).