Optimal Power Flow Formulations for Coordinating Controllable Loads in Low-Voltage Grids: An Overview of Constraint Handling and Hyper Parameter Tuning When Using Metaheuristic Solvers

1. Introduction

Higher penetrations of electrical loads – like for example charging stations for charging electric vehicles or electrical heat pumps for heating houses – will lead to higher electrical loads in future electrical low-voltage distribution grids. Besides this, more fluctuating renewable generators, like for example photovoltaic power systems, will also need to be integrated in these distribution grids. Existing grids might at some point be over-loaded by these additional electrical loads and fluctuating renewable generators. One possible solution to prevent this overload, might be to reinforce existing distribution grids or build new ones. However, this might take a long time and also come along with high costs. Another solution to prevent this overload might be the coordination of those additional electrical loads – given those loads are controllable [1]. Examples for these kinds of loads might be charging stations to charge electric vehicles (EVs) or electrical heat pumps together with thermal energy storages.

For this purpose of coordinating multiple controllable electrical loads in electrical low voltage distribution grids, an optimization problem can be formulated. The goal of this optimization problem formulation is to capture physical laws of grid operation (like power balance, keeping permissable nodal voltage limits or loading limits of power lines and transformers) and also the controllable components (for example minimum run/shut down time, maximum gradients or minimum/maximum power). The outcome of solving the formulated optimization problem would be an optimized time series telling what controllable load might be charged with what power at what timestep. This way an optimal coordination of multiple controllable electrical loads can be realized; keeping requirements (like given heat demands that need to be fulfilled by heat pumps or desired charging levels of EVs) and preventing grid overload, also attributing for fluctuating generation.

These kind of optimization problems are known as optimal power flow (OPF) problems. These OPF problems are non-linear, non-convex (this will be discussed in subsection 3.1) and generally time intensive to solve – especially if integer or binary variables are involved in the OPF formulation (which might be the case for realistic modelling of certain components, like for example minimum/maximum output power of heat pumps) [2]. Analytical solvers – typically based on gradient descent methods – cannot guarantee to find the global optimum of such non-convex optimization problems and might take very long to converge [3]. Alternative solutions for OPF problems are therefore being researched. One very prominent approach is the usage of metaheuristics (this will be discussed in Section 3.2) for solving such OPF problems. These metaheuristics cannot guarantee to find the global optimimum – however, classical analitycal solvers also cannot – and can also find solutions quicker than those analytical solvers [4]. Often, these metaheuristics are based on swarm intelligence – taking inspiration from natural swarms like birds or ants – or also from evolution.

However, metaheuristics only solve unconstrained optimization problems [5] – so can only be used to maximize or minimize a certain objective function without consideration of constraints. But OPF problems consist of many constraints (to ensure proper operation of the considered grid without overload or also proper operation of connected loads and keeping other requirements). So in order to utilize metaheuristics for solving OPF problems, one first has to ensure that constraints are being taken into account. Therefore, there are constraint handling techniques [6]. Besides this, metaheuristics have certain hyper parameters to be tuned [7]. These hyper parameters determine the “performance” of the metaheuristic – so how fast the metaheuristic converges to a result and also the optimality of that result. So besides constraint handling techniques one also needs to tune hyper parameters in order to utilize metaheuristics for solving OPF problems (or generally optimization problems) and to achieve results as optimal as possible.

The goal of this review paper is not to show how the mentioned optimization problem formulation for controlling multiple controllable electrical loads in electrical low voltage distribution grids might look like. And also not to give an overview of metaheuristics typically used for solving these kind of optimal power flow problems; [8] already gives here an up to date comprehensive overview. However, in [8] it is not mentioned which constraint handling techniques were used in the reviewed articles, or how hyper parameters of the metaheuristics were tuned. In contrast, an overview of available constraint handling techniques is given in [5] – however, without connection to OPF problems. So the goal of this review paper is to show what constraint handling techniques are typically used when solving OPF problems using metaheuristics. And how hyper parameters of those metaheuristics are tuned. Besides this, the inclusion of controllable loads and fluctuating renewable generation, like wind and PV is examined. Therefore, a comprehensive literature review is carried out.

The remainder of this review is structured as follows: In Section 2 the methodology for the literature review is described, showing how research was conducted and which criteria are important for the review. Section 3 describes some background regarding optimal power flow, metaheuristics and constraint handling techniques – but just as much as needed to be self-contained. In Section 4 the results of the literature review are shown. Finally, in Section 5 the results are discussed and a conclusion is drawn.

2. Materials and Methods

Relevant literature is found using scientific literature data bases. The data bases considered for this review are:

• Scopus and
• IEEE Xplore.

In these literature data bases, the results are further refined using advanced queries. To keep relevant literature up to date, the results are constrained to years from 2016 to 2024 (both included). Furthermore, only open accessible (gold open access) articles are considered. Relevant keywords are:

• “optimal power flow”,
• “meta heuristic” and
• “constraint handling technique”

The so found articles afterwards are scanned for certain criteria:

• The used optimal power flow formulation (what is the objective of optimization?)
• the used metaheuristic
• the used constrained handling technique
• the used technique for tuning hyper parameters
• are violations of constraints (node voltages, transformer load) monitored?
• are controllable loads or fluctuating renewables considered?
• is a multiple timestep or only a single timestep problem considered?
• which grid topology, at which voltage level is considered?
• how is the load flow formulated?

3. Background

It is not the scope of this article to give a comprehensive overview of existing metaheuristics and constraint handling techniques or the concept of optimal power flow. However, for the sake of being self-contained, these points are introduced here.

3.1. Optimal Power Flow

Optimal power flow problems are a family of optimisation problems that all try to optimize power flows in an electrical grid. There might be many different goals for this optimisation – each associated to a corresponding objective function. However, all formulations have one thing in common: a safe operation of the electrical grid is to be ensured; for example power balance and keeping equipment in the grid in tolerable operational limits. This leads to corresponding constraints. The decision variables are made up of controllable variables and state variables. The control variables can be directly acted upon by grid operators: for example the output power of controllable, thermal power plants or the tap setting of transformers. The state variables are coupled to the control variables: via constraints that describe the operation of the grid. For example node voltages or voltage angles depend on the power drawn from or fed into the grid. Parameters are given as demands that need to be fulfilled and information about the grid itself – which is gathered in form of an admittance matrix, reflecting the physical properties of the grid.

By the nature of physical laws describing the grid operation, these optimal power flow optimisation problems are [9]

non-linear:

either the objective function or constraints cannot be formulated exclusively as linear combination of the decision variables. For example $f\left(\mathbf{x}\right)=x_1+x_2$ is linear, but $f\left(\mathbf{x}\right)=x_1{x}_2$ is not linear. Non-linearity might lead to non-convexnes.

non-convex:

there is not one global optimum, but many local optima (in which solvers might “get stuck”), an example is given in Figure 1.

constrained:

there are constraints involved in the formulation of the optimization problem, for example $x_1<2x_2$.

Figure 1. Example graph of a non-convex function, here $f\left(x\right)=x^4-4x^2+2x$, with local and global optima for minimization. Solvers might “get stuck” in local optima. For ease of comprehension and visualization, this is just a one dimensional example. However, the same also holds for n dimensional examples.

Figure 1. Example graph of a non-convex function, here $f\left(x\right)=x^4-4x^2+2x$, with local and global optima for minimization. Solvers might “get stuck” in local optima. For ease of comprehension and visualization, this is just a one dimensional example. However, the same also holds for n dimensional examples.

A generic non-linear, constrained OPF optimisation problem might look like shown in equation (1):

$$\begin{array}{c}\hfill \begin{array}{cc}\hfill \min .& f\left(\mathbf{x}\right)\hfill \\ \hfill \mathrm{s}.\mathrm{t}.& g_i\left(\mathbf{x}\right)\le 0\forall i\in I\hfill \\ \hfill & h_j\left(\mathbf{x}\right)=0\forall j\in J\hfill \end{array}\end{array}$$

(1)

Here, $f:{\mathbb{R}}^{m+n}\to \mathbb{R}$ is the objective function – in the context of OPF this might typically be some kind of cost metric to be minimised. g and h are inequality and equality constraints, respectively, and I and J are the corresponding sets of all inequality and equality constraints, respectively. $\mathbf{x}\in {\mathbb{R}}^{m+n}$ is a vector of decision variables that may be sub divided in a vector of control variables $\mathbf{u}\in {\mathbb{R}}^m$ and a vector of state variables $\mathbf{v}\in {\mathbb{R}}^n$, such that: $\mathbf{x}={(\mathbf{u},\mathbf{v})}^{⊺}$ [9]. Control variables are those variables that can be directly acted upon – for example active power output of thermal generators, tap settings of tap changing transformers or reactive power fed by shunt var compensators. State variables, on the other hand, are those variables that are affected by the settings of control variables – for example node voltage magnitudes and angles or apparent powers flowing through the lines.

3.1.1. Ways to Calculate Load Flow

To determine optimal power flow, it is important to calculate node voltages and currents on the lines in dependence of load/generation in the grid. This is usually referred to as load flow. From the view point of the optimization problem, the load flow links the state variables to the control variables. In general, there are two ways to calculate this load flow, which differ in their applicability: (1) The “complete” power flow is valid for arbitrary grid topologies and (2) the forward-backward-sweep, which is only valid for radial, non-meshed grids.

Complete Power Flow

The formulation for the complete power flow is valid for arbitrary grid topologies. It is based on the nodel voltage analysis [10]. Provided information on grid topology, given in the form of an admittance matrix $\underline{\mathbf{Y}}\in {\mathbb{C}}^{n\times n}$ for a grid consisting of n buses and currents $\underline{\mathbf{I}}\in {\mathbb{C}}^n$ drawn from/fed into the grid at the n nodes, the voltages $\underline{\mathbf{V}}\in {\mathbb{C}}^n$ at those n nodes, can be calculated according to equation (2).

$$\underline{\mathbf{I}}=\underline{\mathbf{Y}}\underline{\mathbf{V}}$$

(2)

Typically, admittances will be given in cartesian form such that ${\underline{Y}}_{ij}=G_{ij}+\mathrm{j}{B}_{ij}$, where $G_{ij}$ is the conductance and $B_{ij}$ is the susceptance of the $ij$-th element of the admittance matrix. Node votlages are typically given in polar form such that ${\underline{V}}_i=V_{i}exp\left(\mathrm{j}(\omega t+{\delta }_i)\right)$, where $V_i$ is the magnitude and ${\delta }_i$ is the angle of voltage at node i. Formulating active and reactive power balances for each node in the grid, this will lead to formulations as given in equations (3) and (4):

$$P_{G,i}-P_{L,i}-V_i\sum _{j\in N\setminus i}{V}_j\left(G_{ij}cos({\delta }_i-{\delta }_j)+B_{ij}sin({\delta }_i-{\delta }_j)\right)=0\forall i\in N,$$

(3)

where N is the set of all nodes in the grid. $P_{G,i}$ and $P_{L,i}$ denote active power fed into and taken from the grid at node i, respectively. The remaining term calculates the active power flowing through the lines attached to node i.

$$Q_{G,i}-Q_{L,i}-V_i\sum _{j\in N\setminus i}{V}_j\left(G_{ij}cos({\delta }_i-{\delta }_j)-B_{ij}sin({\delta }_i-{\delta }_j)\right)=0\forall i\in N,$$

(4)

where $Q_{G,i}$ and $Q_{L,i}$ denote reactive power fed into and taken from the grid at node i, respectively. The remaining term calculates the reactive power flowing through the lines attached to node i.

Forward-Backward-Sweep

The forward-backward-sweep is only valid for radial grids with no meshes. It is an iterative calculation [11]. Each iteration consists of a forward sweep and a backward sweep. In the first iteration, it is assumed that all node voltage magnitudes are at the nominal level. Given active and reactive power demands/generations ${\underline{S}}_i$ at the nodes as well as assumed voltages ${\underline{V}}_i$ at the nodes, active and reactive currents ${\underline{I}}_i$ flowing out/in the grid at those nodes are computed according to equation (5).

$${\underline{I}}_i={\displaystyle \frac{{\underline{S}}_i}{{\underline{V}}_i}}$$

(5)

Now in the forward sweep, starting at the last node moving towards the first node, the currents ${\underline{I}}_{ik}$ flowing on the lines are calculated according to equation (6):

$${\underline{I}}_{ik}=\sum _{\begin{array}{c}i\in N\\ i\ge k\end{array}}{\underline{I}}_i.$$

(6)

In the following backward sweep, startig at the fist node moving towards the last node, the node voltages are updated using the previously calculated line currents and line impedances ${\underline{Z}}_{ik}=R_{ik}+\mathrm{j}{X}_{ik}$ according to equation (7):

$${\underline{V}}_k={\underline{V}}_i-(R_{ik}+\mathrm{j}{X}_{ik}){\underline{I}}_{ik}.$$

(7)

This will result in node voltages lower than the ones initially assumed. In the following iterations this will lead to accordingly updated line currents. The end of iteration is reached once the difference in node voltages between two consecutive iterations is smaller than a certain delta. A visualization of a radial grid line and corresponding quantities is given in Figure 2

3.1.2. Typical Objectives for Optimal Power Flow Problems

There is not the one and only OPF formulation; there may be different objectives that lead to different formulations of the objective function. Some objectives commonly found in literature are summarized here. At the end of this subsection will be an outlook what a formulation for an objective for coordinating controllable loads might look like.

Minimizing Fuel Costs

The goal here is to minimize the fuel costs for thermal power plants to supply all the required load. This is formulated according to equation (8):

$$\min .costs\left({\mathbf{P}}_G\right)=\sum _{i\in N_{\mathrm{G}}}\left(a_i+b_i{P}_{G,i}+c_i{P}_{G,i}^2\right),$$

(8)

where $N_{\mathrm{G}}$ is the set of all controllable generators attached to the grid. $a_i$, $b_i$ and $c_i$ are coefficients describing the fuel cost for the i-th generator with output power $P_{G,i}$.

Minimizing Active Power Losses

Here, the objective is to minimize active power losses in all transmission lines of the grid. This is formulated according to equation (9)

$$\min .P_{\mathrm{loss}}\left(\underline{\mathbf{V}}\right)=\sum _{i\in N}\sum _{j\in N\setminus i}{G}_{ij}\left(V_i^2+V_j^2-2V_i{V}_{j}cos({\delta }_i-{\delta }_j)\right),$$

(9)

where $V_i$ and $V_j$ are voltage magnitudes as nodes i and j, respectively. Similarly ${\delta }_i$ and ${\delta }_j$ are voltage angles and $G_{ij}$ is the conductance of the $ij$-th element of the system admittance matrix $\underline{\mathbf{Y}}$. N is the set of nodes in the grid.

Minimizing Voltage Deviations

The objective is to minimize voltage deviations across all nodes in the considered grid. This is formulated according to equation (10)

$$\min .\mathrm{Voltage}\mathrm{Deviation}\left(\mathbf{V}\right)=\sum _{i\in N}|1-V_i|$$

(10)

where $V_i$ is the voltage magnitude (in p.u.) at node i and N is the set of nodes in the grid.

Outlook: Possible Objective for Coordinating Controllable Loads

An objective for coordinating multiple controllable loads, as mentioned in the introduction (e. g. charging EVs or electrical heat pumps) might for example look like equation (11)

$$\max .\sum _{i\in N_{\mathrm{EV}}}\left(P_{\mathrm{EV},i}\right)-p\sum _{i\in N_{\mathrm{HP}}}{P}_{\mathrm{slack},i}.$$

(11)

Here, $N_{\mathrm{EV}}$ and $N_{\mathrm{HP}}$ are the sets of nodes with a wallbox and a heat pump connected to it, respectively. $P_{\mathrm{EV},i}$ is the charging power of the EV charging at node i and $P_{\mathrm{slack},i}$ is the amount of thermal energy missing to fulfil the thermal demand of the houshold at node i. p is an according penalty, which will help minimize $P_{\mathrm{slack}}$. Of course, this will also require corresponding equality constraints for balancing thermal demands of households. However, this is just supposed to be a quick outlook, so these additional constraints will not be described here.

3.1.3. Equality Constraints

Equality constraints ensure power balance in the grid, thus enforcing all the required demand is fulfilled by appropriate generation. There are two separate sets of equality constraints: one set for active power balance and another set for reactive power balance. Those constraints were already described in equations (3) and (4) in Section 3.1.1.

3.1.4. Inequality Constraints

Inequality constraints ensure safe operation of the grid, keeping all decision variables in tolerable operating limits. This way it can be ensured that no grid infrastructure is damaged and a safe grid operation is achieved. These constraints are formulated according to equations (12) to (17).

$$\begin{array}{cc}\hfill P_{\mathrm{G},i}^{\min }\le & P_{\mathrm{G},i}\le P_{\mathrm{G},i}^{\max }\forall i\in N_{\mathrm{G}}\hfill \end{array}$$

(12)

$$\begin{array}{cc}\hfill Q_{\mathrm{G},i}^{\min }\le & Q_{\mathrm{G},i}\le Q_{\mathrm{G},i}^{\max }\forall i\in N_{\mathrm{G}}\hfill \end{array}$$

(13)

$$\begin{array}{cc}\hfill V_{\mathrm{G},i}^{\min }\le & V_{\mathrm{G},i}\le V_{\mathrm{G},i}^{\max }\forall i\in N_{\mathrm{G}}\hfill \end{array}$$

(14)

$$\begin{array}{cc}\hfill a_{\mathrm{T},i}^{\min }\le & a_{\mathrm{T},i}\le a_{\mathrm{T},i}^{\max }\forall i\in N_{\mathrm{T}}\hfill \end{array}$$

(15)

$$\begin{array}{cc}\hfill V_i^{\min }\le & V_i\le V_i^{\max }\forall i\in N\hfill \end{array}$$

(16)

$$\begin{array}{cc}\hfill S_i^{\min }\le & S_i\le S_i^{\max }\forall i\in L\hfill \end{array}$$

(17)

Here, $P_{\mathrm{G},i}$, $Q_{\mathrm{G},i}$ and $V_{\mathrm{G},i}$ are active power, reactive power and voltage of the i-th Generator, $N_{\mathrm{G}}$ is the set of all generators. $a_{\mathrm{T},i}$ is the tap setting of the i-th transformer and $N_{\mathrm{T}}$ is the set of all generators. $V_i$ and $S_i$ are voltage at i-th node and apparent power flowing through i-th line, respectively. N and L are the sets of all nodes and all lines, respectively.

3.2. Population-Based Metaheuristics and Why to Use Them

Population-based metaheuristics can be used as an alternative to a classical solver. So a metaheuristic instead of a classical solver is used to solve an optimization problem. To optimally control multiple controllable electrical loads in a low-voltage distribution grid, a corresponding OPF optimization problem can be formulated and solved. As described in subsection 3.1, these OPF problems are constrained, non-linear and non-convex, and thus generally hard to solve. If integer or binary variables are involved (which might be the case for realisitic modeling of components like heat pumps) the time for classical solvers to converge rises exponentially with the number of binary variables [2]. Besides this, classical solvers cannot guarantee to find the global optimum of non-convex problems [12]. Metaheuristic solvers also cannot guarantee to find the global optimum of such problems. However, metaheuristics might converge faster than classical solvers. So metaheuristics might be a viable alternative to classical solvers when it comes to large, non-convex optimization problems including integer or binary variables. Figure 3 shows the connection between the physical grid (including devices to be optimally coordinated), the corresponding OPF formulation and the metaheuristic used to solve the OPF problem.

In contrast to trajectory-based metaheuristics which only consider one solution per Iteration (just like classical solvers) population-based metaheuristics consider multiple solutions per iteration, where each solution equals one member of the population. So, for a population-based metaheuristic a population is initialized. Each member of this population represents one possible solution to the considered OPF problem. So for an OPF problem with n decision variables, one member $\mathbf{x}$ represents one point in an n dimensional search space. So, each member can be thought of as a vector consisting of n elements, where each element represents a concrete value for one of the n decision variables: $\mathbf{x}={\left(x_1,x_2,\cdots ,x_n\right)}^{⊺}$. To draw a connection to the grid, one might think of active/reactive output of controllable generators ($P_{\mathrm{G},i}$ or $Q_{\mathrm{G},i}$), transformer tap settings ($a_{\mathrm{T},i}$), node voltages ($V_i$) or line apparent powers ($S_{ij}$) for example, as depicted for a simple grid in Figure 4.

Typically, control variables are used to make up the members of the population, such that each member represents one vector of concrete values for those control variables [13,14]. So, for the sample grid in Figure 4 each member of the population consists of concrete numerical values for each of the control variables, for example:

$$\mathbf{x}=\left(P_{\mathrm{G},1},P_{\mathrm{G},2},P_{\mathrm{G},3},Q_{\mathrm{G},1},Q_{\mathrm{G},3},Q_{\mathrm{G},3},a_{\mathrm{T}}\right).$$

Each of these members can be assigned a value of the objective function, evaluated at the position of this member: $f\left(\mathbf{x}\right)$.

As an example serves the particle swarm optimization (PSO), first described by [15] in 1995. Based on the swarm behaviour of for example birds, the members of the population “navigate through the topology of the objective function”. Each member updates its current position ${\mathbf{x}}_{i,t}$ based on its own best position ${\mathbf{x}}_{i,t}^{\mathrm{p}}$ found so far and the best position ${\mathbf{x}}_t^{\mathrm{g}}$ found by the entire population so far. According to equation (18) the “velocity” ${\mathbf{v}}_{i,t}$ of member i at iteration t is updated:

$${\mathbf{v}}_{i,t+1}=\xi {\mathbf{v}}_{i,t}+rand{c}_1({\mathbf{x}}_{i,t}^{\mathrm{p}}-{\mathbf{x}}_{i,t})+rand{c}_2({\mathbf{x}}_t^{\mathrm{g}}-{\mathbf{x}}_{i,t}),$$

(18)

where $\xi$, $c_1$ and $c_2$ are referred to as inertia, cognitive parameter and social parameter respectively. These factors determine how much the members orient on their own best solutions or the whole populations best solutions found so far. Thus, these parameters have influence on convergence of the PSO and need to be tuned for proper performance. These parameters are also referred to as hyper parameters. A stochastic element is introduced with rand, which stands for a random number. Then the new position ${\mathbf{x}}_{i,t+1}$ of each member is calculated according to equation (19):

$${\mathbf{x}}_{i,t+1}={\mathbf{x}}_{i,t}+{\mathbf{v}}_{i,t+1}$$

(19)

This way, the members of the population explore the landscape of the objective function. Depending on the hyper parameter settings the members of the population may take different paths, which might lead to more exploration of their environment or exploitation of discovered optima. With help of the stochastic element, members may also be able to leave local optima in case they get stuck. However, constraints are not yet considered. So the optimization problem (which represents the physical grid, see Figure 3) cannot be solved satisfying all the constraints yet. Therefore, there are constraint handling techniques.

3.3. Constraint Handling Techniques

As already stated in Section 3.2, metaheuristics can only optimize unconstrained optimization problems. However, many problems – including OPF – need constraints to be modelled realistically. Therefore, constraint handling techniques (CHTs) are used to account for constraints that get violated by proposed solutions.

3.3.1. Static Penalty Function

With help of a static penalty function the violation of constraints – introduced by solutions proposed by the metaheuristic – can be integrated into the objective function. This way, such solutions get penalized and less attractive [6]. For a constrained optimization problem as shown in equation (1), the penalized objective $F\left(\mathbf{x}\right)$ can be written as equation (20):

$$F\left(\mathbf{x}\right)=f\left(\mathbf{x}\right)+\sum _{i\in I}\left({\lambda }_i^{\mathrm{in}}max(0,g_i\left(\mathbf{x}\right))\right)+\sum _{j\in J}\left({\lambda }_j^{\mathrm{eq}}\left|h_j\left(\mathbf{x}\right)\right|\right),$$

(20)

where ${\lambda }_i^{\mathrm{in}}$ and ${\lambda }_j^{\mathrm{eq}}$ are penalty factors for inequality and equality constraints, respectively. These penalty factors determine how much constraint violations are penalized. The penalty factors also need to be tuned to ensure good performance of the metaheuristic solver. In general, constraints that are deemed more important are associated with a higher penalty factor. The penalized objective $F\left(\mathbf{x}\right)$ is usually referred to as “fitness”.

3.3.2. Superiority of Feasible Solution

With the help of feasibility rules, as described in [16], infeasible solutions can be excluded from the population. Here, a tournament selection operator is used, where two solutions are compared at a time. The “winning” solution to be kept in the population is determined as follows:

• feasible solutions are always preferred to infeasible ones,
• given two feasible solutions, the one with better objective value is preferred and
• given two infeasible solutions, the one with less constraint violation is preferred.

4. Results of the Literature Review

Including all the relevant keywords in the query only lead to two articles in Scopus and zero articles in IEEE Xplore. However, skipping the keyword “constraint handling technique” already 100 articles were found in Scopus and 30 articles in IEEE Xplore. However, some articles were listed in both Scopus and IEEE Xplore. Among those articles the ones relevant have been scanned for the relevant criteria, as outlined in Section 2. The results are listed in Table 1, Table 2, Table 3, Table 4 and Table 5. Each of the criteria mentioned in Section 2 is devoted a separate column in Table 1, Table 2, Table 3, Table 4 and Table 5. Inside those tables, only abbreviations are used. Here is an overview of the used abbreviations:

MH:

The metaheuristic used to solve the OPF problem (as it is not the scope of this review to give a comprehensive overview of metaheuristics, they are just mentioned for the sake of completeness. A prependen (I) means an “improved” version of the base algorithm – according to the authors of the respective article). The following values might appear in tbis column: GWO: Gray Wolf Optimizer, HHO: Harris Hawk Opimizer, MSA: Moth Swarm Algorithm, SSA: Salp Swarm Algorithm, MRF: Manta Ray Foraging Algorithm, SGA: Search Group Algorithm, JAYA: Jaya Algorithm, GA: Genetic Algorithm, SGO: Social Group Optimization, TSO: Transient Search Optimization, GMO: Geometric Mean Optimization, FFO: Firefly Optimization, TFW: Turbulent Flow of Water Optimization, ACO: Ant Colony Optimization, DE: Differential Evolution, CO: Coati Optimization, WSO: War Strategy Optimization, FHO: Fire Hawk Optimization, FPA: Flower Pollination Algorithm, SOA: Skill Optimization Algorithm, PSO: Particle Swarm optimiation, ABC: Artificial Bee Colony Optimization, MGO: Mountain Gazelle Optimizer, GBE: Gradient Bald Eagle Search, BSA: Bird Swarm Algorithm, HSA: Harmony Search Algorithm, GSA: Gravitational Search Algorithm, WHA: Wild Horse Optimization, SFS: Stochastic Fractal Search, MFA: Moth Flame Algorithm, MVO: Multi-Verse Optimization, WOA: Whale Optimization Algorithm, SBB: Satin Bowerbird Optimization, ALO: Ant-lion Optimizer, KHA: Krill Herd Algorithm, AO: Aquila Optimizer, SMO: Slime Mould Optimizer, CO: Coyote Optimization, GHO: Grasshopper Optimization, POA Peafowl Optimization, HGS: Hunger Games Search, AHB: Artificial Hummingbird Optimization, ISA: Interior Search Algorithm, EO: Equilibrium Optimizer, VND: Variable Neighbourhood Descent Algorithm, SFL: Shuffled Frog leaping Optimization, TSA: Tree Seed Optimization, SOS: Symbiotic Organisms Search Algorithm, GNDO: Generalized Normal Distribution Optimizer, COO: Coot Optimizer.

CHT:

The used constraint handling technique to ensure constraints are taken into account when solving the OPF problem with a metaheuristic solver. The following values might appear in this column: SPF: static penalty function, SFS: superiority of feasible solution, LPIM: linear penalty incremental method, ACC: archive-based constraint correction, SAP: self adaptive penalty, ROP: robust oracle penalty, N!: not even mentioned.

HPT:

Hyper parameter tuning to ensure good convergence of the meat heuristic when solving the OPF problem. The following values might appear in this column: SV: at least static values are given, TE: trial and error, N!: not even mentioned.

ObjOPF:

The objective of the OPF problem. The following values might appear in this column: MFC: minimize fuel costs, MFC*: minimize fuel costs (considering valve point effect), MIO: minimize invest and operational costs, MOC: minimize operational costs, MIC: minimize investment costs, MVD: minimize voltage deviation, MPL: minimize active power losses, MPL*: minimize reactive power losses, VSI: maximize voltage stability index, ME: minimize emissions, MES: minimize energy not served. MCC: minimize congestion costs, MPP: maximize PV penetration, TLL: maximize total loadability limit, MCP: minimize costs for changing power output, MSR: minimize system risk, MPI: minimize power import.

CVM:

Whether constraint violations are monitored; for example node voltages or transformer load. The following values might appear in this column: NV: node voltages, N!: not even considered.

MTC:

Whether multiple time steps are considered for the formulation of the OPF problem. The following values might appear in this column: DHR: one day in hourly resolution, YHR: one year in hourly resolution, N!: not considered.

CLC:

Whether controllable loads are considered in the OPF problem formulation (as opposed to just controllable thermal generators). The following values might appear in this column: EV: electric vehicles, Alu: Aluminium plant, N!: not considered.

FRC:

Whether fluctuating renewables are considered in the OPF problem formulation. The following values might appear in this column: PV: photovoltaic, WE: wind energy, HE: hydro energy, BG: Bio gas, N!: not considered.

CGr:

The grid considered in the study. The following values might appear in this column: ERG: existing real-world grid, TGr: test grid (like for example IEEE xxx-bus grids).

GVL:

The voltage level of the considered grid. The following values might appear in this column: HMV: high to medium-voltage (e. g. everything above low-voltage), LV: low-voltage.

CRC:

Whether the results of solving the OPF problem were compared when using different constraint handling techniques (yes or no).

FLF:

How the load flow is formulated. The following values might appear in this column: FB: forward-backward-sweep, PF: complete power flow, N!: not even mentioned.

Table 1. Results of reviewing literature for relevant criteria (MH: metaheuristic, CHT: constraint handling technique, HPT: hyper parameter tuning, ObjOPF: objective of optimal power flow, CVM: constraint violations monitored, MTC: multiple time steps considered, CLC: controllable loads considered, FRC: fluctuating renewables considered, CGr: considered grid, GVL: grid voltage level, CRC: comparison of results for different CHTs, FLF: formulation of load flow)

Ref.	MH	CHT	HPT	ObjOPF	CVM	MTC	CLC	FRC	CGr	GVL	CRC	FLF
[17]	GWO	SPF	N!	MPL, MFC, VSI, MVD	N!	N!	N!	N!	TGr	HMV	No	PF
[18]	MA, AO	N!	N!	MFC, MPL, ME, VSI, MOC	N!	N!	N!	WE	TGr	HMV	No	PF
[19]	GWO	SPF	N!	MFC*, MOC, ME	NV	N!	N!	PV, WE	TGr	HMV	No	PF
[20]	(I)HHO	SPF	N!	MFC*	NV	N!	N!	N!	TGr	HMV	No	PF
[21]	(I)GWO	LPIM, ACC	SV	MFC, VSI, MPL, MVD, ME	N!	N!	N!	N!	TGr	HMV	No	N!
[22]	(I)MSA	SPF	N!	MFC, MVD, VSI	N!	N!	N!	N!	TGr	HMV	No	PF
[23]	(I)SSA	SPF	N!	MPL, MVD, VSI	N!	N!	N!	N!	TGr	HMV	No	PF
[24]	(I)MRF	N!	SV	MCC	NV	N!	N!	N!	TGr	HMV	No	PF
[25]	(I)SGA	SPF	N!	MPL, VSI, MVD	NV	N!	N!	N!	TGr	HMV	No	N!
[26]	JAYA	SPF	N!	MPL	N!	N!	N!	WE	TGr	HMV	No	PF
[27]	GA	SPF	N!	MIO	N!	DHR	N!	PV, WE	ERG	LV	No	PF
[28]	(I)SGO	SPF	SV	MFC*, VSI, MPL, MVD	N!	DHR	EV	N!	TGr	HMV	No	PF
[29]	TSO	N!	SV	MOC	N!	DHR	EV	PV, WE	TGr	HMV	No	PF
[30]	GMO	N!	N!	VSI, MPL, MVD	NV	N!	N!	N!	TGr	HMV	No	N!
[31]	FFO	N!	SV	MPL, MES	NV	YHR	N!	WE	TGr	HMV	No	PF
[32]	TFW	N!	SV	MFC*, ME	N!	N!	N!	N!	N!	N!	No	N!
[33]	(I)ACO	ROP	N!	MPI	N!	DMIN	EV	PV	ERG	LV	No	PF
[34]	DE	SFS, SAP	SV	MFC, MFC*, VSI, MPL, ME	NV	N!	N!	N!	TGr	HMV	Yes	PF

Table 1. Results of reviewing literature for relevant criteria (MH: metaheuristic, CHT: constraint handling technique, HPT: hyper parameter tuning, ObjOPF: objective of optimal power flow, CVM: constraint violations monitored, MTC: multiple time steps considered, CLC: controllable loads considered, FRC: fluctuating renewables considered, CGr: considered grid, GVL: grid voltage level, CRC: comparison of results for different CHTs, FLF: formulation of load flow)

Ref.	MH	CHT	HPT	ObjOPF	CVM	MTC	CLC	FRC	CGr	GVL	CRC	FLF
[17]	GWO	SPF	N!	MPL, MFC, VSI, MVD	N!	N!	N!	N!	TGr	HMV	No	PF
[18]	MA, AO	N!	N!	MFC, MPL, ME, VSI, MOC	N!	N!	N!	WE	TGr	HMV	No	PF
[19]	GWO	SPF	N!	MFC*, MOC, ME	NV	N!	N!	PV, WE	TGr	HMV	No	PF
[20]	(I)HHO	SPF	N!	MFC*	NV	N!	N!	N!	TGr	HMV	No	PF
[21]	(I)GWO	LPIM, ACC	SV	MFC, VSI, MPL, MVD, ME	N!	N!	N!	N!	TGr	HMV	No	N!
[22]	(I)MSA	SPF	N!	MFC, MVD, VSI	N!	N!	N!	N!	TGr	HMV	No	PF
[23]	(I)SSA	SPF	N!	MPL, MVD, VSI	N!	N!	N!	N!	TGr	HMV	No	PF
[24]	(I)MRF	N!	SV	MCC	NV	N!	N!	N!	TGr	HMV	No	PF
[25]	(I)SGA	SPF	N!	MPL, VSI, MVD	NV	N!	N!	N!	TGr	HMV	No	N!
[26]	JAYA	SPF	N!	MPL	N!	N!	N!	WE	TGr	HMV	No	PF
[27]	GA	SPF	N!	MIO	N!	DHR	N!	PV, WE	ERG	LV	No	PF
[28]	(I)SGO	SPF	SV	MFC*, VSI, MPL, MVD	N!	DHR	EV	N!	TGr	HMV	No	PF
[29]	TSO	N!	SV	MOC	N!	DHR	EV	PV, WE	TGr	HMV	No	PF
[30]	GMO	N!	N!	VSI, MPL, MVD	NV	N!	N!	N!	TGr	HMV	No	N!
[31]	FFO	N!	SV	MPL, MES	NV	YHR	N!	WE	TGr	HMV	No	PF
[32]	TFW	N!	SV	MFC*, ME	N!	N!	N!	N!	N!	N!	No	N!
[33]	(I)ACO	ROP	N!	MPI	N!	DMIN	EV	PV	ERG	LV	No	PF
[34]	DE	SFS, SAP	SV	MFC, MFC*, VSI, MPL, ME	NV	N!	N!	N!	TGr	HMV	Yes	PF

Table 2. Results of reviewing literature for relevant criteria; continuation of Table 1 (MH: metaheuristic, CHT: constraint handling technique, HPT: hyper parameter tuning, ObjOPF: objective of optimal power flow, CVM: constraint violations monitored, MTC: multiple time steps considered, CLC: controllable loads considered, FRC: fluctuating renewables considered, CGr: considered grid, GVL: grid voltage level, CRC: comparison of results for different CHTs, FLF: formulation of load flow)

Ref.	MH	CHT	HPT	ObjOPF	CVM	MTC	CLC	FRC	CGr	GVL	CRC	FLF
[35]	MSA	SPF	SV	MFC, MFC*, MPL, VSI, MVD	NV	N!	N!	N!	TGr	HMV	No	PF
[36]	CO, WSO	N!	N!	MFC, MFC*, MVD	NV	N!	N!	WE	TGr	HMV	No	N!
[37]	FHO	N!	N!	MFC, MVD, MPL	NV	N!	N!	PV, WE	TGr	HMV	No	N!
[38]	(I)FPA	SFS	N!	MFC*, MPL, ME, MVD	NV	N!	N!	PV, WE	TGr	HMV	No	PF
[39]	GWO, HHO	N!	N!	MFC, ME, MPL, MVD	NV	N!	N!	N!	TGr	HMV	No	PF
[40]	SOA	N!	N!	MFC	NV	N!	N!	N!	TGr	HMV	No	N!
[41]	PSO, ABC, DE	SPF	N!	MPP, MVD	NV	N!	N!	PV	TGr	HMV	No	PF
[42]	(I)ACO	SPF	N!	MFC	N!	N!	N!	N!	TGr	HMV	No	N!
[43]	GWO, FPA	N!	SV	TLL	NV	N!	N!	N!	TGr	HMV	No	PF
[44]	MGO	N!	N!	MFC, MPL, MVD	N!	N!	N!	N!	TGr	HMV	No	PF
[45]	(I)HHO	SPF	SV	MFC, ME, MPL, MVD	N!	N!	N!	N!	TGr	HMV	No	PF
[46]	ACO	N!	SV	MPL, MVD	NV	N!	N!	N!	TGr	HMV	No	N!
[47]	GBE	N!	SV	MFC, MFC*, MRC, MIO	NV	N!	EV	WE, PV	TGr	HMV	No	PF
[48]	BSA, JAYA	N!	N!	MFC, MFC*, ME, MPL, MVD	NV	N!	N!	N!	TGr	HMV	No	PF
[49]	(I)GA, HSA	N!	N!	MPL, TLL	NV	DHR	N!	N!	TGr	HMV	No	PF

Table 2. Results of reviewing literature for relevant criteria; continuation of Table 1 (MH: metaheuristic, CHT: constraint handling technique, HPT: hyper parameter tuning, ObjOPF: objective of optimal power flow, CVM: constraint violations monitored, MTC: multiple time steps considered, CLC: controllable loads considered, FRC: fluctuating renewables considered, CGr: considered grid, GVL: grid voltage level, CRC: comparison of results for different CHTs, FLF: formulation of load flow)

Ref.	MH	CHT	HPT	ObjOPF	CVM	MTC	CLC	FRC	CGr	GVL	CRC	FLF
[35]	MSA	SPF	SV	MFC, MFC*, MPL, VSI, MVD	NV	N!	N!	N!	TGr	HMV	No	PF
[36]	CO, WSO	N!	N!	MFC, MFC*, MVD	NV	N!	N!	WE	TGr	HMV	No	N!
[37]	FHO	N!	N!	MFC, MVD, MPL	NV	N!	N!	PV, WE	TGr	HMV	No	N!
[38]	(I)FPA	SFS	N!	MFC*, MPL, ME, MVD	NV	N!	N!	PV, WE	TGr	HMV	No	PF
[39]	GWO, HHO	N!	N!	MFC, ME, MPL, MVD	NV	N!	N!	N!	TGr	HMV	No	PF
[40]	SOA	N!	N!	MFC	NV	N!	N!	N!	TGr	HMV	No	N!
[41]	PSO, ABC, DE	SPF	N!	MPP, MVD	NV	N!	N!	PV	TGr	HMV	No	PF
[42]	(I)ACO	SPF	N!	MFC	N!	N!	N!	N!	TGr	HMV	No	N!
[43]	GWO, FPA	N!	SV	TLL	NV	N!	N!	N!	TGr	HMV	No	PF
[44]	MGO	N!	N!	MFC, MPL, MVD	N!	N!	N!	N!	TGr	HMV	No	PF
[45]	(I)HHO	SPF	SV	MFC, ME, MPL, MVD	N!	N!	N!	N!	TGr	HMV	No	PF
[46]	ACO	N!	SV	MPL, MVD	NV	N!	N!	N!	TGr	HMV	No	N!
[47]	GBE	N!	SV	MFC, MFC*, MRC, MIO	NV	N!	EV	WE, PV	TGr	HMV	No	PF
[48]	BSA, JAYA	N!	N!	MFC, MFC*, ME, MPL, MVD	NV	N!	N!	N!	TGr	HMV	No	PF
[49]	(I)GA, HSA	N!	N!	MPL, TLL	NV	DHR	N!	N!	TGr	HMV	No	PF

Table 3. Results of reviewing literature for relevant criteria; continuation of Table 2 (MH: metaheuristic, CHT: constraint handling technique, HPT: hyper parameter tuning, ObjOPF: objective of optimal power flow, CVM: constraint violations monitored, MTC: multiple time steps considered, CLC: controllable loads considered, FRC: fluctuating renewables considered, CGr: considered grid, GVL: grid voltage level, CRC: comparison of results for different CHTs, FLF: formulation of load flow)

Ref.	MH	CHT	HPT	ObjOPF	CVM	MTC	CLC	FRC	CGr	GVL	CRC	FLF
[50]	JAYA	SPF	N!	MFC, MPL, VSI	N!	N!	N!	N!	TGr	HMV	No	PF
[51]	FPA	N!	N!	MPL, MIC	NV	N!	N!	N!	TGr	HMV	No	FB
[52]	(I)GSA	SPF	N!	MFC, MPL, MVD	NV	N!	N!	WE, PV	TGr	HMV	No	PF
[53]	(I)WHA	N!	SV	MPL	NV	N!	N!	N!	TGr	HMV	No	FB, PF (?!)
[14]	FFA	SPF	N!	MCP	NV	N!	N!	N!	TGr	HMV	No	PF
[54]	(I)HHO	SPF	SV	MFC, ME, MPL	N!	N!	N!	N!	TGr	HMV	No	PF
[55]	SFS	SPF	N!	MPL, MVD, VSI	N!	N!	N!	N!	TGr	HMV	No	PF
[56]	(I)GWO, MFA, SSA, MVO	N!	N!	MIO, MPL, ME	NV	N!	N!	WE, PV, HE	TGr	HMV	No	N!
[57]	WOA	N!	N!	MIC, MPL	N!	N!	N!	N!	TGr	HMV	No	PF
[58]	(I)GWO	SPF	N!	MFC, MFC*	NV	N!	N!	N!	TGr	HMV	No	PF
[13]	(I)GWO	N!	N!	MPL, MIO	NV	N!	N!	N!	TGr	HMV	No	PF
[59]	SBB	SPF	SV	MCC	NV	N!	N!	N!	TGr	HMV	No	PF
[60]	FPA	N!	N!	MPL, MVD	NV	N!	N!	N!	TGr	HMV	No	FB
[61]	(I)FFO	SPF	N!	MFC, MVD, VSI, MPL, MPL*	NV	N!	N!	N!	TGr	HMV	No	PF
[62]	ALO	SPF	SV	MFC, MVD, VSI, MPL, MPL*	NV	N!	N!	N!	TGr	HMV	No	PF
[63]	MRF	N!	N!	MPL, MVD, VSI	NV	N!	N!	N!	TGr	HMV	No	FB
[64]	KHA	N!	N!	MFC*, MPL, ME	N!	N!	N!	WE	TGr	HMV	No	PF

Table 3. Results of reviewing literature for relevant criteria; continuation of Table 2 (MH: metaheuristic, CHT: constraint handling technique, HPT: hyper parameter tuning, ObjOPF: objective of optimal power flow, CVM: constraint violations monitored, MTC: multiple time steps considered, CLC: controllable loads considered, FRC: fluctuating renewables considered, CGr: considered grid, GVL: grid voltage level, CRC: comparison of results for different CHTs, FLF: formulation of load flow)

Ref.	MH	CHT	HPT	ObjOPF	CVM	MTC	CLC	FRC	CGr	GVL	CRC	FLF
[50]	JAYA	SPF	N!	MFC, MPL, VSI	N!	N!	N!	N!	TGr	HMV	No	PF
[51]	FPA	N!	N!	MPL, MIC	NV	N!	N!	N!	TGr	HMV	No	FB
[52]	(I)GSA	SPF	N!	MFC, MPL, MVD	NV	N!	N!	WE, PV	TGr	HMV	No	PF
[53]	(I)WHA	N!	SV	MPL	NV	N!	N!	N!	TGr	HMV	No	FB, PF (?!)
[14]	FFA	SPF	N!	MCP	NV	N!	N!	N!	TGr	HMV	No	PF
[54]	(I)HHO	SPF	SV	MFC, ME, MPL	N!	N!	N!	N!	TGr	HMV	No	PF
[55]	SFS	SPF	N!	MPL, MVD, VSI	N!	N!	N!	N!	TGr	HMV	No	PF
[56]	(I)GWO, MFA, SSA, MVO	N!	N!	MIO, MPL, ME	NV	N!	N!	WE, PV, HE	TGr	HMV	No	N!
[57]	WOA	N!	N!	MIC, MPL	N!	N!	N!	N!	TGr	HMV	No	PF
[58]	(I)GWO	SPF	N!	MFC, MFC*	NV	N!	N!	N!	TGr	HMV	No	PF
[13]	(I)GWO	N!	N!	MPL, MIO	NV	N!	N!	N!	TGr	HMV	No	PF
[59]	SBB	SPF	SV	MCC	NV	N!	N!	N!	TGr	HMV	No	PF
[60]	FPA	N!	N!	MPL, MVD	NV	N!	N!	N!	TGr	HMV	No	FB
[61]	(I)FFO	SPF	N!	MFC, MVD, VSI, MPL, MPL*	NV	N!	N!	N!	TGr	HMV	No	PF
[62]	ALO	SPF	SV	MFC, MVD, VSI, MPL, MPL*	NV	N!	N!	N!	TGr	HMV	No	PF
[63]	MRF	N!	N!	MPL, MVD, VSI	NV	N!	N!	N!	TGr	HMV	No	FB
[64]	KHA	N!	N!	MFC*, MPL, ME	N!	N!	N!	WE	TGr	HMV	No	PF

Table 4. Results of reviewing literature for relevant criteria; continuation of Table 3 (MH: metaheuristic, CHT: constraint handling technique, HPT: hyper parameter tuning, ObjOPF: objective of optimal power flow, CVM: constraint violations monitored, MTC: multiple time steps considered, CLC: controllable loads considered, FRC: fluctuating renewables considered, CGr: considered grid, GVL: grid voltage level, CRC: comparison of results for different CHTs, FLF: formulation of load flow)

Ref.	MH	CHT	HPT	ObjOPF	CVM	MTC	CLC	FRC	CGr	GVL	CRC	FLF
[65]	ABC, MFO	N!	N!	MSR, MIO	NV	DHR	N!	WE, HE	TGr	HMV	No	PF
[66]	AO	N!	SV	MOC	N!	N!	N!	WE	TGr	HMV	No	PF
[67]	(I)DE	N!	SV, TE	MFC, MVD, VSI, MPL	NV	N!	N!	WE, PV	TGr	HMV	No	PF
[68]	SMO	SFS	N!	MOC, ME	NV	N!	N!	PV, WE	TGr	HMV	No	PF
[69]	(I)CO	SPF	N!	MFC*, MPL	NV	N!	N!	PV	TGr	HMV	No	PF
[70]	MVO, GHO, HHO	SPF	N!	MFC, MPL, MVD	NV	N!	N!	N!	TGr	HMV	No	PF
[71]	POA	SPF	N!	MFC, MPL, MVD, ME	NV	N!	N!	N!	TGr	HMV	No	PF
[72]	HGS	N!	N!	MFC, MPL, ME, MVD, VSI	NV	N!	N!	N!	TGr	HMV	No	PF
[73]	ALO, MFO, SSO	N!	SV	MIC, TLL	N!	N!	N!	N!	TGr	HMV	No	N!
[74]	(I)AHB	SPF	SV	MVD, MPL, ME, MFC	NV	N!	N!	N!	TGr	HMV	No	PF
[75]	MRF	SPF	SV, TE	MPL, ME, MFC, MVD	NV	N!	N!	WE, PV	TGr	HMV	No	PF
[76]	WOA, GA	N!	N!	MFC	NV	N!	N!	N!	TGr	HMV	No	PF
[77]	ISA	SPF	SV	MFC, MFC*, MVD	NV	N!	N!	N!	TGr	HMV	No	PF
[78]	EO	SPF	SV	MPP	N!	DHR	N!	PV	TGr	HMV	No	N!

Table 4. Results of reviewing literature for relevant criteria; continuation of Table 3 (MH: metaheuristic, CHT: constraint handling technique, HPT: hyper parameter tuning, ObjOPF: objective of optimal power flow, CVM: constraint violations monitored, MTC: multiple time steps considered, CLC: controllable loads considered, FRC: fluctuating renewables considered, CGr: considered grid, GVL: grid voltage level, CRC: comparison of results for different CHTs, FLF: formulation of load flow)

Ref.	MH	CHT	HPT	ObjOPF	CVM	MTC	CLC	FRC	CGr	GVL	CRC	FLF
[65]	ABC, MFO	N!	N!	MSR, MIO	NV	DHR	N!	WE, HE	TGr	HMV	No	PF
[66]	AO	N!	SV	MOC	N!	N!	N!	WE	TGr	HMV	No	PF
[67]	(I)DE	N!	SV, TE	MFC, MVD, VSI, MPL	NV	N!	N!	WE, PV	TGr	HMV	No	PF
[68]	SMO	SFS	N!	MOC, ME	NV	N!	N!	PV, WE	TGr	HMV	No	PF
[69]	(I)CO	SPF	N!	MFC*, MPL	NV	N!	N!	PV	TGr	HMV	No	PF
[70]	MVO, GHO, HHO	SPF	N!	MFC, MPL, MVD	NV	N!	N!	N!	TGr	HMV	No	PF
[71]	POA	SPF	N!	MFC, MPL, MVD, ME	NV	N!	N!	N!	TGr	HMV	No	PF
[72]	HGS	N!	N!	MFC, MPL, ME, MVD, VSI	NV	N!	N!	N!	TGr	HMV	No	PF
[73]	ALO, MFO, SSO	N!	SV	MIC, TLL	N!	N!	N!	N!	TGr	HMV	No	N!
[74]	(I)AHB	SPF	SV	MVD, MPL, ME, MFC	NV	N!	N!	N!	TGr	HMV	No	PF
[75]	MRF	SPF	SV, TE	MPL, ME, MFC, MVD	NV	N!	N!	WE, PV	TGr	HMV	No	PF
[76]	WOA, GA	N!	N!	MFC	NV	N!	N!	N!	TGr	HMV	No	PF
[77]	ISA	SPF	SV	MFC, MFC*, MVD	NV	N!	N!	N!	TGr	HMV	No	PF
[78]	EO	SPF	SV	MPP	N!	DHR	N!	PV	TGr	HMV	No	N!

Table 5. Results of reviewing literature for relevant criteria; continuation of Table 4 (MH: metaheuristic, CHT: constraint handling technique, HPT: hyper parameter tuning, ObjOPF: objective of optimal power flow, CVM: constraint violations monitored, MTC: multiple time steps considered, CLC: controllable loads considered, FRC: fluctuating renewables considered, CGr: considered grid, GVL: grid voltage level, CRC: comparison of results for different CHTs, FLF: formulation of load flow)

Ref.	MH	CHT	HPT	ObjOPF	CVM	MTC	CLC	FRC	CGr	GVL	CRC	FLF
[79]	(I)JAYA	SPF	N!	MFC, ME, MPL, MVD	NV	N!	N!	N!	TGr	HMV	No	N!
[80]	VND	SPF	N!	MFC	NV	N!	N!	N!	TGr	HMV	No	PF
[81]	(I)PSO	N!	SV	MPL, MOC, MVD	NV	N!	N!	N!	TGr	HMV	No	PF
[82]	SFL, TSA	N!	N!	MPL, MVD, VSI	NV	N!	N!	N!	TGr	HMV	No	PF
[83]	(I)ACO	N!	N!	MFC	NV	N!	N!	N!	TGr	HMV	No	PF
[84]	SOS	N!	N!	MFC, MPL, MVD, VSI	NV	N!	N!	N!	TGr	HMV	No	PF
[85]	(I)GNDO	SPF	N!	MOC, MVD, VSI, ME, MPL	NV	N!	N!	WE	TGr	HMV	No	PF
[86]	(I)FPA	N!	N!	MPL	N!	N!	N!	PV, WE, BG	TGr	HMV	No	FB
[87]	COO	SPF	SV, TE	MPL, ME, MVD	N!	N!	Alu	PV, HE	TGr	HMV	No	PF

Table 5. Results of reviewing literature for relevant criteria; continuation of Table 4 (MH: metaheuristic, CHT: constraint handling technique, HPT: hyper parameter tuning, ObjOPF: objective of optimal power flow, CVM: constraint violations monitored, MTC: multiple time steps considered, CLC: controllable loads considered, FRC: fluctuating renewables considered, CGr: considered grid, GVL: grid voltage level, CRC: comparison of results for different CHTs, FLF: formulation of load flow)

Ref.	MH	CHT	HPT	ObjOPF	CVM	MTC	CLC	FRC	CGr	GVL	CRC	FLF
[79]	(I)JAYA	SPF	N!	MFC, ME, MPL, MVD	NV	N!	N!	N!	TGr	HMV	No	N!
[80]	VND	SPF	N!	MFC	NV	N!	N!	N!	TGr	HMV	No	PF
[81]	(I)PSO	N!	SV	MPL, MOC, MVD	NV	N!	N!	N!	TGr	HMV	No	PF
[82]	SFL, TSA	N!	N!	MPL, MVD, VSI	NV	N!	N!	N!	TGr	HMV	No	PF
[83]	(I)ACO	N!	N!	MFC	NV	N!	N!	N!	TGr	HMV	No	PF
[84]	SOS	N!	N!	MFC, MPL, MVD, VSI	NV	N!	N!	N!	TGr	HMV	No	PF
[85]	(I)GNDO	SPF	N!	MOC, MVD, VSI, ME, MPL	NV	N!	N!	WE	TGr	HMV	No	PF
[86]	(I)FPA	N!	N!	MPL	N!	N!	N!	PV, WE, BG	TGr	HMV	No	FB
[87]	COO	SPF	SV, TE	MPL, ME, MVD	N!	N!	Alu	PV, HE	TGr	HMV	No	PF

4.1. Used Constraint Handling Techniques

Most of the reviewed articles use a static penalty function as shown in equation (20), to integrate constraint violation into the objective function, then the resulting fitness function is optimized. Only very few articles use different constraint handling techniques. Also, many articles don’t even mention constraint handling techniques – even though they have to be used in order to achieve feasible solutions to OPF problems. The results are shown in Figure 5.

A typical example for an article using static penalty function for handling constraints is [22]. Bounds of decision variables are formulated as inequality constraints. Violation of those bounds are penalized via penalty factors and added on top of the original objective function, as described by equation (20).

Only one article compares the results of optimization for different constraint handling techniques: In [34] authors compare “superiority of feasible solution” and “self-adadptive penalty” for handling constraints as well as an ensemble of both methods. Optimizations are carried out on IEEE 30-, 57- and 118-bus systems with either single or also multiple objectives like for example minimizing fuel costs, power losses or voltage deviations. A differential evolution was used for solving the OPF problem, combined with different CHTs. Statistical analyses of the objective function values are carried out using Wilcoxon signed rank test. However, the results show that no single CHT is able to produce best results for all considered cases.

Other CHTs like linear penalty incremental method (LPIM), archive-based constraint correction (ACC), robust oracle penalty (ROP) were only used very seldomly. In [21] the authors used LPIM and ACC as CHTs. Optimizations were carried out on IEEE 30- and 57-bus systems using a grey wolf optimizer based on symbiotic learning. Objectives of the OPF formulation were, for example, minimization of fuel costs, emissions, power losses or voltage deviation – either considered as single objectives or multiple objectives via weighted sum. For constraint handling the authors use a combination of LPIM and ACC, with LPIM as primary CHT. After a certain number of iterations ACC might be used as secondary CHT, based on probability. In the LPIM, the penalty factors – as described in equation (20) – are linearly increased with the number of iterations. Thus, in the first iterations penalties are smaller, leading to more exploration, while at higher iterations penalties are bigger, leading to more exploitation of found optima. In ACC an archive stores solutions with the least constraint violations so far. In order for new solutions to enter this archive, their constraint violation must be less than of those solutions stored in the archive.

In [33] the authors use a robust oracle penalty for handling constraints. A micro grid consisting of five nodes is considered for optimization. The objective of the OPF formulation is to minimize power import from the main grid. The authors don’t expand on the explanation of ROP, however, further information can be found in [88]: ROP is based on only one parameter, the oracle. Ideally, this oracle is to be set slightly greater than the optimal feasible solution – which of course is not known beforehand. Starting with an initial guess, the oracle gets updated each iteration, moving closer to an optimum, thus resulting in a kind of self-tuning effect. For further details, the reader is referred to [88].

4.2. Used Techniques for Hyper Parameter Tuning

Many of the reviewed articles don’t even mention how the hyper parameters for the metaheuristic solver were set. The other articles at least showed what values they have set for the hyper parameters of their used metaheuristic. Only three articles stated that “optimal” values for hyper parameters were determined based on trial and error. The results are shown in Figure 6.

In [67] an IEEE 30-bus grid is considered for optimization. The objectives of the OPF are to simultaneously minimize fuel costs and power losses in transmission lines. For solving the OPF problem the authors use a hybrid of differential evolution and symbiotic organisms search. The authors state that the optimal hyper parameter settings were obtained “[a]fter several trial runs of the algorithm”.

In [75] IEEE 30- and 118-bus systems are considered for optimization. Here, the objectives of the OPf formulation include minimizing fule costs, emissions, voltage deviation and power losses in transmission lines. A manta ray foraging optimization is used to solve the OPF problem. According to the authors “[t]he best solution (optimal values) for each parameter was chosen”. The authors provide a table with a testing range for each hyper parameter of different metaheuristics. However, it is not explained, how the best solution was determined – probably also by trial and error in the given ranges.

4.3. Objectives of the Optimal Power Flow Formulation

Most of the reviewed articles consider minimization of power losses in transmission lines, voltage deviation or fuel costs as objective. The results are shown in Figure 7. In most of the cases either single objectives or multiple objectives are considered. That is the reason why the number of articles in Figure 7 sum up to more than the total number of reviewed articles.

In the case of multiple objectives, two options are possible. Either a weighted sum of the single objectives, which leads to a scalar multi-objective function, as is done for example in [39]. Or calculating a pareto front based on a vector-valued objective function, as for example described in [18]. The mostly used objectives are to minimize power losses in transmission lines, voltage deviation or fuel costs for thermal generators. But also sometimes objectives like minimizing congestion costs or increasing total loadability limit are included in OPF formulations. Most of the times these objectives are achieved by corresponding coordination of controllable thermal units. But also sometimes for example placement or sizing of flexible AC transmission systems (FACTS) devices is used to achieve those objectives or distributed generation is considered.

In [35] different objectives like minimizing fule costs, power losses in transmission lines, voltage deviations and emissions or maximizing voltage stability index are considered. The authors consider different cases, either single objectives or multiple objectives as a weighted sum. This way up to four objectives are considered simultaneously. Optimizations are carried out on IEEE 30-, 57- and 18-bus systems using moth swarm algorithm as solver. The objectives are achieved using controllable thermal generators.

In [61] the authors consider different objectives like minimizing fuel costs, voltage deviation, active and reactive power losses in transmission lines or maximizing voltage stability index. However, here only one objective at a time is considered. Optimizations are carried out for an IEEE 30-bus grid. A hybrid firefly particle swarm optimization is used to solve the OPF problem. Here, also controllable thermal generators are used to achieve the objectives.

In [49] the objectives of the OPF formulation are to minimize power losses in transmission lines and to maximize total loadability limit of the grid. Both objectives are considered simultaneously using a pareto front. The Optimizations are carried out on an IEEE 30-bus grid. The authors use a genetic and a harmony search algorithm. The objectives are achieved by optimal placement and sizing of FACTS devices.

Similarly, in [57] the optimal placement and sizing of FACTS devices is used to achieve certain objectives. Here, this objective is to minimize investment and operating costs for those FACTs devices. In a first step, the optimal placement of those FACTS devices is manually determined. In a next step, the optimal sizing of those FACTS devices is determined using whale optimization. Here, IEEE 14- and 30-bus systems are considered for the optimization.

In [50] the objectives of the OPF formulation are to minimize fuel costs, power losses in transmission lines and maximizing voltage stability index. The objectives are considered one at a time. Here, IEEE 30- and 118-bus grids are considered for the optimization. To solve the OPF problem, the authors use the Jaya algorithm. Using statistical analyses, also the effect of distributed generation on the performance of the solver is considered.

It is already noticeable here that none of the articles examined pursues the goal described in the introduction (to enable the coordination of as many controllable loads as possible, such as electric heat pumps and charging stations in existing low-voltage grids)

4.4. Methods for Formulating the Load Flow

Besides the different objectives for the OPF problem, as mentioned in Section 4.3, there are also differences in the way how the load flow is incorporated into the OPF formulation. However, there are only two different approaches in the reviewed literature: 1.) for radial non-meshed grids a forward-backward sweep according to equations (5)–(7) can be employed or 2.) for arbitrary (also meshed) grids a “complete” power flow formulation according to equations (3) and (4) can be employed. Finally, both formulations yield the state variables in dependence of the control variables. The results of the literature review are shown in Figure 8

Articles like [63] or [60] use the forward-backward-sweep because they only consider radial, non-meshed grids. In those grids, solving the forward-backward-sweep iterations converges faster than solving Gauss-Seidel or Newton-Raphson for the corresponding full power flow formulation. Most other articles that consider meshed grids, have to use the full power flow formulation, because meshed grids cannot be formulated using forward-backward-sweep.

However, there are also some articles that don’t even mention how load flow is formulated, like for example [30] or [42]. Such papers just formulate balances for active and reactive powers, but don’t describe how this effects the load flow (there is no connection between control and state variables). Also one article [53] describes both approaches and it is not obvious which one is used (probably forward-backward-sweep because of the radial grid).

4.5. Consideration of Controllable Loads

Most of the reviewed articles consider the optimization of controllable, thermal generators. Only in four articles electric vehicles were considered – however, not as controllable load. But just as additional load that has to be supplied. Only one article deviates from this pattern. All those articles have in common that demand for EV charging is determined by probability distribution functions. In [28] additional EV charging in IEEE 30 bus and IEEE 57 bus networks are considered. A certain penetration of EVs (and also PHEVs) with battery capacities and initial SOCs are assumed according to normal distributions. This results in an additional power demand added on top of residential load profiles. The OPF problem with objectives like minimization of fuel costs, voltage deviations, power losses and maximization of voltage stabiliy index is solved with an improved social group optimization. Since charging profiles are fixed in advance according to probability density functions, EV charging cannot be optimized anymore. Therefore, the work does not allow a coordination of EVs in the sense of controllable loads.

In [29] additional EV charging as well as fluctuating production of renewables is considered in an islanded AC-DC hybrid microgrid based on an IEEE 33-bus test grid. Production of renewables as well as demand for charging EVs are based on according probability distribution functions. Using monte-carlo simulations and a fast forward scenario reduction technique additional loads and charging demands are projected over 24 hours for multiple scenarios. The optimization with the objective of minimizing operational costs is carried out using a transient search optimization. Since charging profiles for EVs are calculated beforehand using probability density functions, for the optimization those charging profiles are just fixed parameters and cannot be optimized. Thus, the work in this article does not allow a coordination of EVs in the sense of controllable loads.

In [47] EV charging and also vehicle to grid (V2G) together with fluctuating wind and PV production are considered in an IEEE 30-bus grid. Again, renewables production and EV load profiles are determined beforehand the optimization according to given probability density functions obtained after multiple monte-carlo simulations. The objectives of the OPF formulation are to minimize power losses on transmission lines, voltage deviations, and total operation costs. A gradient bald eagle optimization is used for solving the OPF problem. Since the charging pattern for the EV is already predefined according to a given probability distribution, the EV charging cannot be optimized. Thus, also this work does not allow a coordination of EVs as controllable loads.

In [33] EV charging in a 5-bus microgrid is considered. There are also fluctuating PV production and other components like a storage in this micro grid. The goal of optimization is to minimize the energy import from the main grid. For this purpose a storage and charging EVs are utilized. Time slots when cars are available as well as starting and desired finish SOCs are randomly generated based on normal distributions. As a results, EVs and the storage are charged in such a way that the energy imported from the main grid is minimized. In this sense this article is unique in the fact that it tries to utilize EVs as controllable load – instead of just considering predefined load profiles for EVs.

In [87] the operation of aluminium plants is considered as controllable load. The grid under consideration is an IEEE 57-bus grid. There are also renewable sources like PV and hydro energy feeding into the grid. The goal of optimization is to minimize carbon emissions, power losses and voltage deviation. The aluminium plants can adapt their power factor in order to consume more or less reactive power, thus affecting voltage levels. However, it is not apparent how exactly this is integrated into the OPF formulation. Still, this is considered as controllable loads. So this paper in unique in the fact that it tries to optimize power factor of controllable loads in order to achieve the optimization goals.

4.6. Consideration of Fluctuating Generation of Renewables

Most of the reviewed articles only consider generation of thermal generators. However, some articles also include fluctuating generation of renewables like photovoltaic, wind energy or hydro energy. The results are shown in Figure 9. Typically, wind turbine and photovoltaic generations are modelled according to specific probability distribution functions.

In [18] stochastic wind energy generation in a IEEE 30-bus grid with two wind parks is considered. The wind energy generation is modelled according to a Weibull distribution. The parameters of the Weibull distribution are calculated are calculated using mayfly algorithm and aquila optimizer. Based on the resulting probability distribution function an OPF problem with multiple objectives like minimizing emissions of thermal generators or power losses in transmission lines and total operation costs is formulated. The costs for wind energy are divided in for over- and underestimating the actual wind energy generation. The OPF problem is solved by the mayfly algorithm. As a result, fluctuating wind energy generation is incorporated by means of according control of thermal generators.

In [41] the authors consider optimal placement and sizing of PV in IEEE 13- and 37-bus systems. For the PV generation a given load profile is used. The objective of the OPF formulation is to minimize voltage deviation and maximize the PV penetration. The authors use different metaheuristics like artificial bee colony, particle swarm optimization or differential evolution to solve the OPF problem. As a result, renewable generation of PV is incorporated by means of optimal placement of PV in the grid.

In [37] stochastic wind and PV generation in an IEEE 69-bus grid are considered. Furthermore, a diesel generator and a micro turbine are included in the grid. Wind is modelled according to a Weibull distribution and solar generation is modelled using a normal probability density function. The objective of the OPF problem is to minimize generation costs, power losses in transmission lines and voltage deviation. The authors don’t describe closer how the costs for PV and wind production are calculated. The locations of PV and wind turbine in the grid are are fixed. A fire hawk optimization is used to solve the OPF problem. Wind and PV production are thus incorporated by adapting other controllable loads.

4.7. Considered Grids and Voltage Levels

Most of the reviewed articles don’t consider real, existing grids, but rather test grids, like for example IEEE xxx-bus (where xxx denotes the number of buses, there are many different test grids with different numbers of buses). Many articles also consider multiple such test grids with varying numbers of nodes, to show scalability of their considered metaheuristic. Also, most articles don’t explicitly mention the voltage level. However, those test grids are mostly used to simulate high-voltage grids. Only very few papers consider for example real micro grids.

One example for an article with a micro grid is [27]: the authors consider a micro grid consisting of 8 nodes. However, there is no further information on the voltage level. Here, it is assumed this is a low-voltage grid. There is also no information whether this is a real existing grid, the authors call it a “generalizable [micro grid]”.

In [29] the authors consider a hybrid AC-DC radial distribution grid based on an IEEE 30-bus test grid.

[33] considers a real existing microgrid at the Wroclaw University in Poland. The microgrid consists of 5 nodes, however, the authors don’t explicitly mention the voltage level, but it is assumed that it is a low-voltage grid.

In [46] the authors consider an IEEE 30 bus test grid working at 12.6 kV

5. Discussion/Conclusion

Most of the reviewed articles focus on an optimal coordination of controllable, thermal generators, given certain inflexible loads. Typical objectives of the OPF problem are minimization of voltage deviation, power losses in transmission lines or fuel costs. The results of such optimizations are settings for control variables, like generator active power or transformer tap settings as well as state variables like node voltages or angles. However, almost all OPF formulations in the reviewed literature only consider a single time step problem.

Most of those articles compare the performance of their used metaheuristic solvers to other benchmarks in the literature, thus developing efficient solvers. However, none of the reviewed articles research the effect of hyper parameter tuning on the performance of the metaheuristic solver. Often there are only fixed values given for hyper parameters, but not how those values were determined. In fact, only three papers stated that “optimal” values for hyper parameters were determined by trial and error. And only one article examines the effects of different constraint handling techniques on the performance of the metaheuristic solver.

A few of the reviewed articles also included fluctuating production of renewables like PV, wind energy or hydro energy in their OPF formulations. The stochastic character of renewable production is incorporated by using according probability density functions. However, all articles consider either optimal placement of newly to be install renewables, subject to typical objectives like minimizing operational costs or effects on grid operation or integration of renewables by means of shifting thermal generation accordingly. None of the reviewed articles consider how to include existing plants into operation of existing low-voltage distribution grids.

Regarding the considered grid for optimization, almost all of the reviewed articles exclusively focus on test grids like IEEE xxx-bus (where xxx stands for a concrete number of buses). These are no real existing grids, but are mostly used to evaluate the performance of algorithms used for solving OPF problems. Only two of the reviewed articles examine a real existing grid. Regarding the voltage level of the considered grids, there is no clear indication to be found. This is because it is customary to specify the voltages in p.u. However, due to the structure (multiple branches, partly meshed) of the test grids, it is assumed that those grids are operated in high to middle-voltage. It is therefore also not clear whether the optimization techniques used in the articles reviewed can be directly transferred to low-voltage grids.

Also, almost none of the reviewed articles consider controllable loads in low-voltage distribution grids. Few articles consider electric vehicles, however, almost exclusively as static loads added on top of residential load profiles. Only one article considers electric vehicles as real controllable load – by including charging power as decision variables in the OPF formulation. Also one article considers aluminium plants as controllable loads, by means of adapting their power factor. Only a few other articles consider flexible AC transmission system (FACTS) devices which might be considered as “controllable loads”. However, loads like electric vehicles or electric heat pumps – with associated requirements like departure time or SOC or heat demands to be covered – were never considered. There exists literature like for example [89,90,91], that examine optimal coordination of controllable loads, however such articles don’t consider the grid load.

For an optimal scheduling of controllable electric loads in low-voltage distribution grids, accordingly extended OPF formulations must be investigated. These formulations should include requirements associated with controllable loads: for example desired departure time and SOC for charging electric vehicles or heat demands that need to be covered by electric heat pumps. But also fluctuating production of local renewable generation like PV should be incorporated in such formulations. To allow a meaningful coordination of controllable loads, the optimization problem also has to cover multiple time steps. Thus, additional decision variables and constraints have to be added to the optimization problem. For example constraints to ensure power balance for all time steps, calculate SOCs of EVs or operating limits of other state variables and grid infrastructure.

However, such expanded OPF formulations also require more time for potential solvers to converge to a solution. This is especially true for formulations that involve integer or binary variables. A possible solution can be the usage of metaheuristic solvers – as is already the current state of research. However, only few articles research the effect of different constraint handling techniques or hyper parameter settings on the performance of the metaheuristic solver. Also, there are no efforts in finding an automated way to determine optimal hyper parameter settings.

Put all together, directions for future research could include:

• whether the reviewed optimization techniques (which were almost always applied in high to middle-voltage grids) can also be applied in low-voltage grids as they are,
• extended OPF formulations (with according constraints) that also account for controllable loads and consider multiple time steps,
• statistical evaluation of the performance of metaheuristic solvers for different constraint handling techniques and settings of hyper parameters and
• research on automated – e. g. based on machine learning – methods to determine optimal hyper parameters that maximize the performance of the metaheuristic solver.

Covering those research gaps might help to achieve an optimal coordination of multiple, controllable loads in low-voltage distribution grids – in a manner so that all requirements of those loads are considered and also a safe operation of the grid is ensured.