Algorithms for Optimal Power Flow Extended to Controllable Renewable Systems and Loads

In an effort to quantify and manage uncertainties inside power systems with penetration of renewable energy, uncertainty costs have been defined and different uncertainty cost functions have been calculated for different types of generators and electric vehicles. This article seeks to use the uncertainty cost formulation to propose algorithms and solve the problem of optimal power flow extended to controllable renewable systems and controllable loads. In a previous study, the first and second derivatives of the uncertainty cost functions were calculated and now an analytical and heuristic algorithm of optimal power flow are used. To corroborate the analytical solution, the optimal power flow was solved by means of metaheuristic algorithms. Finally, it was found that analytical algorithms have a much higher performance than metaheuristic methods, especially as the number of decision variables in an optimization problem grows.


Introduction
In recent years, the advances in technology, the reduction in manufacturing costs and the growing concern for the environment, unconventional renewable energy sources have had a great penetration in electrical power systems. For example, in Colombia, it is expected that 1050 MW of wind energy will come into operation by 2023, according to information published by the Mining Energy Planning Unit (UPME) [1]. The start-up of these new non-conventional renewable energy sources poses a challenge to power systems, with regard to the dispatch of generation resources. One of the biggest drawbacks of energy sources such as solar or wind is their high variation over time and their stochastic nature. In addition, advances in technology such as large-scale energy storage, the massification of electric vehicles and the increase in distributed generation must be taken into account when making economic dispatches.
According to the review of the state of the art, it is possible to address the problem of uncertainties of power systems with different approaches, whether probabilistic, possibilistic or probabilistic-possibilistic hybrids [2]. Economic dispatch problems are approached mainly by probabilistic methods, while possibilistic and probabilistic-possibilistic methods are especially useful for generating generation forecasts [3]. In an effort to quantify the impact of uncertainties in power systems, from the calculation of expected values it has been possible to determine the cost of uncertainty in non-conventional renewable energy sources, starting from the probability distributions associated with the random variables that affect the behavior of the system. Arévalo et al. [4] presents the calculation of the cost of uncertainty associated with solar and wind plants and the behavior of electric vehicles; for this case it was assumed that the irradiance for solar plants was associated with a lognormal function, the wind speed with a Rayleigh distribution and the power delivered by electric vehicles had a Gaussian behavior. Complementing the work on the cost of uncertainty, in the article written by Molina et al. [5], the cost of uncertainty was calculated for a run-of-river hydroelectric power station, whose generation depends on the flow of the river from which the station feeds; it is exposed that the behavior of river flows is associated with a Gumbell distribution.
In many cases, you may not have historical measurements of the resource you want to take advantage of (for example, in remote, non-interconnected areas). For these cases, a uniform distribution of the resource to be exploited could be assumed. Bernal et al. [6] propose the calculation of the cost of uncertainty for uniform distributions. The uncertainty in the loads can be modeled by a normal distribution [2]. Vargas et al. [7] describe the calculation of the uncertainty cost for controllable loads. The uncertainty cost calculation for controllable loads is also valid for any other load, since controllable loads have normal behavior. It is necessary to clarify that the results of the uncertainty cost functions in the previously mentioned works [4][5][6][7] were validated using the Monte Carlo method. Now, to achieve a reliable operation, and so that the beginning of the operation about non-conventional renewable energy sources can be reflected in an economic benefit, it is necessary to perform Optimal Power Flows (OPFs) which consider costs of uncertainty. In the literature, there are models of optimal power flow that consider different costs of uncertainty and are solved analytically [3,8] or through heuristic techniques [9][10][11][12][13][14].
The OPFs have different applications and are generally differentiated by objective functions, restrictions, the type of network to be optimized, and so forth. In relation to OPF in Hetzer et al. [3] an optimal power dispatch is presented considering two wind generators and their associated uncertainty cost in a system with two conventional generators. The work from Hetzer et al. was supplemented by Zhao et al. [11] when solving an OPF using particle swarm optimization (PSO), considering the uncertainty costs of electric vehicles and wind generators in a 118 node IEEE power system. From the previously mentioned works [3,4,11], Arévalo et al. [9] show the calculation of an OPF that takes into account the cost of uncertainty in wind and solar generators and that of the large-scale entry of electric vehicles. In this work, the 118-node IEEE system was used and the PSO algorithm was used to solve the optimization problem.
Torres and Rivera [15] propose the calculation of an OPF in various operating scenarios of the system. In their work, a simplified model of the Colombian network was considered, as well as the future agglomeration of solar and wind resources on the Caribbean coast [16]. These same authors, in a different work [12], showed the calculation of an OPF using the DEEPSO method (combination of particle swarm and differential evolution) for the IEEE system of 118 nodes with costs of uncertainty from wind generators, solar plants and electric vehicles. The optimal power flows mentioned so far have been formulated for large systems and their sources of uncertainty have been solar radiation, the speed of time and/or the delivery of energy to the grid by electric vehicles. However, with the penetration of distributed generation, the possibility arises that some clients of the system satisfy part or all of their demand for a certain time, becoming a controllable load. Guzmán et al. [13] describe the controllable loads and their types of contract with the energy marketer. In addition to this, they calculate an OPF solved by the DEEPSO algorithm where they take into account controllable loads, position of the transformer taps and compensation in the IEEE system of 118 nodes.
The uncertainty cost functions, the basis for calculating the uncertainty costs, are usually complex since they appear to be non-elementary integrals [4,5,7]. One possible way to deal with this difficulty is the one proposed by Martinez and Rivera [8], where a quadratic approximation of the uncertainty cost of electric vehicles, solar, wind and run-ofriver plants is made to carry out an analytical calculation of the OPF using MATPOWER. Due to the reduction in the costs of equipment for the generation of energy from nonconventional renewable sources, it is now possible that small consumers of energy at medium and low voltages can also generate energy on a small scale, giving rise to the appearance of microgrids. Peña et al. [14] performed the calculation of an OPF for a microgrid with different sources of uncertainty where there was also energy storage. Using a genetic algorithm (NSGA-II) they managed to minimize the cost of power generation and maximize the useful life of the batteries. Similarly, Li et al. [10] propose an OPF where the operating cost of a microgrid with batteries is minimized. Unlike the other works mentioned so far, Li et al. did not use the Monte Carlo method for the validation of their results, but instead used the method of the estimated point. In this way, the novelty of this paper is in including the marginal uncertainty cost functions (MUCFs) in the algorithms for optimal power flow extended to controllable renewable systems and controllable loads, since MUCFs have been tested with Monte Carlo simulations.

Optimization Problem Statement
Optimal power flow is, as its name implies, an optimization problem defined as follows [17,18]: where u is a set of decision variables and x is a set of dependent variables. F is the scalar objective function, h represents a set of equality constraints and g a set of inequality constraints.

Target Function
The objective function evaluated in this proposal is shown below, in the expression (4): The variables presented in the objective function are described below: N gen,C Number of conventional generators, N gen,B Number of batteries, N gen,S Number of solar generators, N gen,E Number of wind generators, N gen,H Number of run-of-river hydraulic generators, N gen,CL Number of controllable loads or nodes with electric vehicle connection, P j Active power delivered by the conventional generator j, P k Active power delivered by battery k, P m Active power delivered by the solar generator m, P n Active power delivered by wind generator n, P r Active power delivered by the hydraulic generator r, P t Active power dispatched by controllable load or electric vehicle station t.
The functions f j and f k refer to the cost functions for conventional generators and batteries, respectively. On the other hand, the functions f m , f n , f r and f t refer to uncertainty cost functions for solar, wind, run-of-the-river hydro generators and electric vehicles/controllable loads, respectively.
The cost functions for conventional generators correspond to quadratic polynomial functions of the form [19]: On the other hand, the cost functions for batteries correspond to linear functions. In this work, only batteries are taken into account in their discharge cycle, which means, working as generators [20,21]: In the expression (6), cs m is the battery operation and maintenance cost and PD is the power delivered by the battery.
The functions f m , f n , f r and f t refer to cost functions of uncertainty due to the variable behavior of the resources from which energy is extracted, such as solar radiation, wind speed, water flow or the number of connected electric vehicles.
The aforementioned functions are the sum of the costs of uncertainty due to underestimation and the costs of uncertainty due to overestimation, as shown below: In the expressions from (7)-(10), the subscript u refers to the cost function of uncertainty due to underestimation while the subscript o refers to the cost function of uncertainty due to overestimation.

Concept of Uncertainty Cost Functions from Previous Studies
In order to calculate Uncertainty Costs Functions (UCFs), it is necessary to define underestimation and overestimation costs developed in previous studies [4,11] : Uncertainty Cost Due to Underestimation Costs due to underestimation refer to the power that a renewable generation unit cannot deliver to the grid when the scheduled power value of the plant is smaller than the available generation power: where P Sch and P Av are the scheduled power and the available power, respectively. In this case, penalty cost due to underestimation is given by: where c u is the penalty cost coefficient due to underestimation and P max is the generator maximum output power. Now, because of the variability of the renewable power sources, the power generated by these sources has an associated Probability Density Function (PDF) f n (P). The uncertainty cost due to underestimation is defined as the expected value of C sub (developed from Expression (12)): Uncertainty Cost Due to Overestimation Costs due to overestimation refer to the power that cannot be supplied by a renewable generator because the available power is smaller than the previously scheduled power: In this case, penalty cost due to overestimation is given by: where c o is the penalty cost coefficient due to overestimation and P min is the generator minimum output power.
In the same way as the underestimation condition, and based on the stochastic nature of renewable sources, the uncertainty cost due to overestimation is given by the expected value of C so (developed from Expression 15): Finally, the Uncertainty Cost Function (UCF) for a given renewable source is equal to the sum of underestimation and overestimation costs (developed from Expressions (13) and (16) respectively): In [11], the authors of this paper presented the development of the Uncertainty Cost functions of PVG, WEG, PEV and RHG and the Formulation and application of marginal cost functions of PVG, WEG, PEV, and RHG. They are useful as target functions in this paper, and the marginal cost functions are used in the analytical algorithm applied in this paper.

Constraints
Below are the constraints that are taken into account for optimal power flow [18]. The active power limits of each generator must be respected: The reactive power limits of each generator must be respected: The voltage limits in each system bar must be respected: The chargeability of any line or transformer should not be exceeded. For this restriction, the maximum apparent power of each element, the input power and the output power are taken into account: S p,out ≤ S max,p p = 1...n lines (22) S t,in ≤ S max,t p = 1...n trans f ormers (23) S t,out ≤ S max,t p = 1...n trans f ormers .
Finally, for each node k of the system, the balance of active and reactive power must be maintained, that is, the load flow [19]:

DEEPSO Algorithm
In this paper, two metaheuristic algorithms were used for the validation of the results obtained by the analytical method. The first one is the DEEPSO (Differential Evolutionary Particle Swarm Optimization) algorithm and will be studied in this subsection. The other algorithm used in the validation of results is the genetic algorithm, which will be described in the next subsection.
The DEEPSO algorithm is an improvement of the PSO (Particle Swarm Optimization) algorithm, to which evolutionary and differential computing techniques are added [22]; it is for this reason that, for a better understanding of the algorithm implemented in this thesis, we will begin by describing the algorithm PSO.

PSO Algorithm
The Particle Swarm Optimization (PSO) algorithm is an algorithm based on the behavior of fish or bird populations, which either seek to avoid a predator or search for food [23]. There is an initial population of particles of size (number of elements) N and dimension (number of variables) D, which is denoted as follows: The Particle X moves to a new position X (2) according to the following rule [22][23][24]: where V is known as the velocity of the particle and is defined as: The expressions (29) and (30) can be rewritten in a vector manner as follows for each element of the population [24]: In the expression (32), w represents a scalar value that descends from w max = 0.9 to w min = 0.4 as the number of iterations of the algorithm progresses [24], c 1 and c 2 are constants inherent to the method and r 1 and r 2 are random values uniformly distributed among 0 and 1.
The operation of the PSO algorithm consists of making a search for new optimal values by each particle using the expressions (31) and (32), where a new position of particle i , taking into account the best historical position of each particle P (k) best i . That is, the value of X i that has given the best value of the objective function and the best historical position of the total set of particles . This is illustrated in Figure 1.

EPSO Algorithm
The EPSO (Evolutionary Particle Swarm Optimization) algorithm is a hybrid of the PSO algorithm and evolutionary techniques. According to what was proposed by Vladimiro Miranda, the algorithm consists of the following steps [22,25]: Each particle is replicated r times, usually r is equal to 1; 2.
Parameters A, B, C of expression (30) for the r replicas are mutated; 3.
Each of the r + 1 particles creates a new generation from the particle motion rules (29) and (30); 4.
The value of the objective function is evaluated for each particle of the new generation; 5.
Through some selection procedure, the best descendants of each ancestor are selected to form a new generation.

DEEPSO Algorithm
The DEEPSO algorithm is a variation of the EPSO algorithm in which a modified form of the speed equation (30) shown below is used [22]: In the expression (33), b * G is given by: That is to say, b * G is affected by normal noise with average 0 and standard deviation w G . P is a random matrix of ones and zeros, where each element has a 75% probability of being 1 and a 25% probability of being 0 [22]. X r1 and X r2 are any two points of the population that define four types of implementation of the DEEPSO algorithm that are defined in [22]. For the case of this study, the implementation of the algorithm used was called DEEPSO Pb-rnd, in which X and each row of the matrix b (k) r1 are chosen from P b , the particle vector with the best historical values; on the other hand, X (k) r2 = X (k) and for the case of minimization, the following must be fulfilled [22]: The flow chart for the DEEPSO algorithm implemented to solve the OPF is shown in Figure 2.

Genetic Algorithm
The genetic algorithm is a metaheuristic technique initially conceived for solving discrete variable problems. A modification of the genetic algorithm designed to solve continuous variable problems was used for this work [26]. In the following subsections, the implemented genetic algorithm will be described.

Genetic Algorithm Pseudocode
The pseudocode of the genetic algorithm implemented in order to solve the optimal power flow is presented below: 1.
Generate an initial random population; 2.
Evaluate the objective function for each element of the initial population; 3.
Carry out a selection scheme to choose the individuals with the greatest probability of reproduction; 4.
Evaluation of the objective function for each element of the new population; 7.
Selection of the individuals with the best objective function values.

Generation of the Initial Population
To generate the initial population, we proceeded in a similar way to the case of the DEEPSO algorithm, that is, values of variables that met the proposed restrictions were randomly selected.

Evaluation of the Objective Function of the Initial Population
For the genetic algorithm, the implemented objective function, which from now on will be known as the fitness function, is shown in the expression (36): Thus, if a restriction is violated, the fitness function will assign an infinite value to the cost function, which is known as a penalty.

Selection
For the selection, the tournament scheme was used in references [27,28], in which k groups of m elements were established, where each element m is the value of the objective function associated with a given individual of the population. Since this work seeks to solve a minimization problem, in each group k the m elements are ordered from highest to lowest and a probability of recombination and reproduction is assigned to each element depending on the position they occupy in the array according to the expression (37). P(pos) = 2pos n(n + 1) ; pos = 1, 2, 3..., m.

Reproduction and Recombination
For the reproduction phase, two different parents are randomly selected from the entire population according to the reproductive probabilities given in the selection phase. The actual recombination strategy used is shown below, which was taken from [26]. There is a vector α, shown in the following expression.
where n is the number of decision variables of the optimization problem. Each value α i is randomly and uniformly distributed in the interval [−γ, 1 + γ] with γ = 0.1. Now, we have the parent vectors x1 and x2 that are shown in the expressions (39) and (40).
The recombination procedure for each element of y 1 and y 2 is presented in expressions (43) and (44): Finally, it is clarified that as many parents as elements in the initial population are selected and each pair of parents has a pair of children, so that the new generation (children) has the same size as the initial population and the total population would be twice the size of the initial population size.

Mutation
The mutation takes place for a percentage of individuals of the new generation, that is, not all individuals of the new generation are mutated. For mutated individuals, each variable is modified by adding a normally distributed value with mean 0 and standard deviation σ.

Evaluation and Selection
For the new generation, the fitness function is assessed. Finally, all the existing particles are ordered, that is, first and second generation from highest to lowest and 50% of the population with the lowest objective function value is selected. This reduces the size of the population to its initial value.

Genetic Algorithm Flow Diagram
The flow chart of the genetic algorithm is shown in Figure 3.

Analytical Algorithm: Interior Point Method
In this subsection, a description of the interior point method will be presented, since this method was the analytical technique used to solve the optimal power flow in view of its ability to solve non-linear optimization problems. The description presented here is illustrative of the Matpower MIPS solver method [18]. It should be clarified that the solvers used to solve the optimal power flow can use different variations of the method [18,29,30].
The choice of the interior point method is due to a limitation on the part of the software that was used to solve the optimal power flow, since solvers that use Matlab code must be used to implement the non-linear cost functions, its gradients (marginal costs) and its Hessians (derivatives of marginal costs) [18].

Karush-Kuhn-Tucker Conditions
Karush-Kuhn-Tucker conditions are necessary conditions for a point to be a restricted local optimum [31].
If x * is a local maximum value for the following problem [31]: The gradients of the constraints at the optimum ∇g i (x * ) are independent, then there is a vector λ * = [λ 1 , λ 2 , ..., λ m ] T such that [31]:

Lagrangian Function
The Lagrangian function is defined as shown below [31]: Thus, if we only have equality restrictions, it is possible to rewrite the Karush-Kuhn-Tucker conditions in terms of the Lagrangian function (53) as shown in the expressions (54) and (55) [31]. This fact is the basis of the interior point method.

MIPS based Optimization
The optimization problem under study is redefined below [18]: Subject to: G(X) = 0 (57) This optimization problem then becomes a problem where there are only equality restrictions through the inclusion of slack variables Z and the sum of a barrier function to the objective function: Subject to: G(X) = 0 (60) As the perturbation parameter γ approaches zero, the solution of the previous problem approaches the solution of the original problem. The Lagrangian for a given value of γ from the previous optimization problem is then: From expression (63), the gradient to each of the variables is obtained: It should be noted that e is a vector of ones and [A] represents a diagonal matrix with vector A on the diagonal.
To satisfy the Karush-Kuhn-Tucker conditions, the values of X, Z, λ and µ must be found which make expressions (64) to (67) equal to zero as γ approaches zero. To determine these values, the Newton-Raphson method is used, which produces the following matrix equation: Since it is possible that the variations ∆X, ∆Z ∆λ and ∆µ take the variables to a new infeasible point, they are usually truncated by α coefficients as shown below: There are different ways to determine the values of α such as functions of merit, confidence regions, and filtering methods [31]. Finally, the value of α must be updated which, as mentioned, should decrease as the number of iterations increases.

Results: Optimal Power Flow Solution Considering Heuristic and Analytical Algorithms
For the solution of the optimal power flow consider three test systems: The IEEE system of nine nodes and the IEEE system of 57 nodes. This section shows the results of the optimal power flow solution in each of the aforementioned systems and a validation and comparison of the results obtained by analytical and heuristic methods.

Scenarios Analyzed Scenario 1
Scenario 1 consists of keeping the generators connected to nodes 1 and 2 as they are established in the original case, that is, with quadratic cost functions; however, the generator connected to node 3 is assumed to be a hydraulic generator with the cost functions defined in [32].

Scenario 2
Scenario 2 now contemplates the connection of a solar plant to node 3, replacing the existing conventional generator. The other generators are not modified.

Scenario 3
In the same way as in scenarios 1 and 2, the conventional generator connected to node 3 is replaced by a wind generator.

Scenario 4
In this scenario, a new controllable load is added to node 5 with a maximum value equal to −30 MW with an average of −19.54 MW and a standard deviation of 0.54 MW.
Tables 1-4 show the parameters of the renewable generators and the controllable load.    The results obtained for all the scenarios in the nine-node system are observed in Tables 5-8 where the active and reactive power values, the voltage at the slack node, the execution time and the value of the objective function are shown for each optimization technique or solver with which the optimal power flow was solved.  It is noteworthy that, for scenario 1, as seen in Table 5, the DEEPSO algorithm yielded a slightly lower value than the analytical solver FMINCON. It is also observed that the analytical solver FMINCON outperforms heuristic algorithms in terms of computation times. On the other hand, it is also observed that the IPOPT and MIPS solvers did not converge for this scenario.
In scenarios 2 and 3 shown in Tables 6 and 7, the three solvers analyzed reached the local optimum value. The MIPS solver had the best performance in terms of execution time. For the evaluation of the objective function the three solvers reached very similar answers.
In scenario 4, shown in Table 8, the IPOPT solver did not converge, while FMINCON had an excellent performance in terms of computational time. From the analyzed scenarios, the superiority in terms of computation times of the analytical methods with respect to the metaheuristic algorithms is clearly observed.
Although the genetic algorithm converged much faster than the DEEPSO algorithm, the results of the DEEPSO algorithm are much closer to the results obtained by the analytical methods; this is mainly observed in the reactive power values obtained for each generator. It is also observed in the analytical methods that, although the active power values in generators produced by each method are too close to each other, there are small differences in the reactive power values.

Fifty Seven-Node IEEE System Description
The IEEE 57-bus test system is used to evaluate the optimization problem of this study. Based on details given in [4] for system buses and branches, the data of the system have been structured in the MATPOWER [6] data format. Branch thermal limits were defined based on reference values given in [5]. A summary of the characteristics of the test system are: seven generators, 42 loads, 63 line/cables, 15 transformers step/wise, two transformers fixed taps and three shunt compensation binary On/Off. Figure 5 below shows the 57-node system on which solar, wind, run-of-the-river generators, controllable loads and batteries were mounted.

Scenarios Analyzed
For the 57-node system, only one scenario was analyzed that took into account all types of generation. It begins by clarifying that for this system the slack node is node 1. The nodes whose generators were changed by generators with uncertain costs are listed below. The nodes to which some type of generation was added and that previously did not have it are also listed.

•
The generator connected to node 3 is replaced by a wind generator with uncertain costs; • The generator connected to node 6 is replaced by a solar generator with uncertain costs; • The generator connected to node 9 is replaced by a run-of-river hydraulic generator with uncertain costs; • A battery bank is connected to node 57; • A controllable charging or electric vehicle charging station is connected to node 47.
The parameters of the cost functions of the generators with cost of uncertainty and of the batteries are shown in Tables 9-12 (data from previous research: [32]):   The results obtained are presented in Table 13, where it is observed that the optimal power flow solution could be found using all optimization techniques. It is also observed that, for the heuristic methods, 17 decision variables are being taken, namely, active and reactive powers in all generators, except for the slack and the voltage of the slack node. Again the interior point method has superiority with respect to the speed of solving the problem, since the metaheuristic techniques took more than 10 min to solve the problem. From the results shown in Table 13, it is possible to see that, as the number of variables in a problem grows, analytical methods have a greater advantage over the heuristics in terms of the speed with which they solve the optimal flow of power and due to the fact that they are not subject to random variation computation of the first iteration of the variable decisions and the computationally high costs for large systems.
It is also observed that the analytical optimum presents a value of the objective function that is lower than that presented by heuristic methods. Furthermore, in this case it is observed that the DEEPSO algorithm does not have such good behavior as to determine the highest value of the objective function, as is observed in the 57 buses system. Additionally, it is also observed that the genetic algorithm has a lower value of the objective function and takes less time than the DEEPSO algorithm. In this way, the appropriate optimization technique for solving the OPF was determined. The inclusion of the cost functions of renewable systems was also considered.

Conclusions
In this paper, an optimal power flow was developed that takes into account the uncertainties produced by the penetration of renewable energy sources and controllable loads. It was observed that including the Jacobian (First derivatives) and the Hessian (Second derivatives) of the uncertainty cost functions in the optimal power flow model brings great benefits, since it allows the use of analytical techniques for the determination of active and reactive power dispatches. The main benefit of analytical techniques over the metaheuristic techniques used in the state-of-the-art for the solution of the optimal power flow is the reduction of the computation times that was evidenced in all the cases analyzed, as shown in Table 13. Additionally, it was observed that the IPOPT solver did not converge to a solution in several of the systems and scenarios proposed. On the other hand, the Matlab FMINCON solver converged in all cases to a solution. The MIPS solver did not converge in scenario 1 of the nine-node system.
On the other hand, heuristic methods were used in order to validate the results obtained through analytical methods and very close responses were achieved for the active power dispatches in the IEEE systems of nine and 57 nodes, as shown in Tables 5-8 and 13. However, for the 118-node system, the heuristic methods yielded significantly different dispatch values with respect to the analytical methods; in addition, the dispatches produced by the analytical methods produced a lower value in the objective function. Finally, the presented formulation and its respective solution allow network operators to have a tool to manage the controllable renewable resources found in the network, satisfying the different physical restrictions of the system.
In this way, in this paper the marginal costs of uncertainty are used to determine the active power values that a particular generator with uncertainty must deliver to minimize the uncertainty costs. The field of optimization in electrical power systems is in continuous development. In the case of controllable renewable sources and controllable loads, it is possible to extend this analysis to optimal power flows in DC. It is also possible to use the flow extension optimal power for controllable renewable systems and controllable loads and problems such as Unit Commitment or Security Constraints Optimal Power Flow.
To solve the proposed OPF, when controllable renewable systems and loads are considered, it is also necessary to take into account all the constraints that can affect the power system. The solution would allow the system to work between specified constraints considering uncertainty cost functions, making this problem very complex. Including all constraints in the traditional OPF interior point solution method could prevent convergence and tackling the problem using metaheuristic strategies makes the problem hard to solve within short periods of time (critical clearing times of the contingency in the case of the security constraint optimal power flow) without high computational power.
In this study, the developed research relies on two aspects: the first is analytical methods, using marginal uncertainty cost functions which allow us to use the interior point method for the optimization problem. The second is the solution of the problem using metaheuristic methods, in cases where it is not possible to determine the marginal cost functions. The use of marginal uncertainty cost functions (MUCF) strategies can be employed in small power systems, using strong assumptions of the analytical modeling of cost functions, which makes some of the applications infeasible in real operations or very costly in terms of hardware implementation for real systems. As a solution for those limitations, future research will address the problem using the potential of Parallel and Heterogeneous Computing (PHC).