abstract: The power system is the largest man-made physical network in the world. Performing analysis of a large bulk system is computationally complex, especially when the study involves engineering, economic and environmental considerations. For instance, running a unit-commitment (UC) over a large system involves a huge number of constraints and integer variables. One way to reduce the computational expense is to perform the analysis on a small equivalent (reduced) model instead on the original (full) model. The research reported here focuses on improving the network reduction methods so that the calculated results obtained from the reduced model better approximate the performance of the original model. An optimization-based Ward reduction (OP-Ward) and two new generator placement methods in network reduction are introduced and numerical test results on large systems provide proof of concept. In addition to dc-type reductions (ignoring reactive power, resistance elements in the network, etc.), the new methods applicable to ac domain are introduced. For conventional reduction methods (Ward-type methods, REI-type methods), eliminating external generator buses (PV buses) is a tough problem, because it is difficult to accurately approximate the external reactive support in the reduced model. Recently, the holomorphic embedding (HE) based load-flow method (HELM) was proposed, which theoretically guarantees convergence given that the power flow equations are structure in accordance with Stahl’s theory requirements. In this work, a holomorphic embedding based network reduction (HE reduction) method is proposed which takes advantage of the HELM technique. Test results shows that the HE reduction method can approximate the original system performance very accurately even when the operating condition changes. / Dissertation/Thesis / Doctoral Dissertation Electrical Engineering 2017

This thesis considers three separate but connected problems regarding energy networks: the load flow problem, the optimal power flow problem, and the islanding problem. All three problems are non-convex non linear problems, and so have the potential of returning local solutions. The goal of this thesis is to find solution methods to each of these problems that will minimize the chances of returning a local solution. The thesis first considers the load ow problem and looks into a novel approach to solving load flows, the Holomorphic Embedding Load Flow Method (HELM). The current literature does not provide any HELM models that can accurately handle general power networks containing PV and PQ buses of realistic sizes. This thesis expands upon previous work to present models of HELM capable of solving general networks efficiently, with computational results for the standard IEEE test cases provided for comparison. The thesis next considers the optimal power flow problem, and creates a framework for a load flow-based OPF solver. The OPF solver is designed with incorporating HELM as the load flow solver in mind, and is tested on IEEE test cases to compare it with other available OPF solvers. The OPF solvers are also tested with modified test cases known to have local solutions to show how a LF-OPF solver using HELM is more likely to find the global optimal solution than the other available OPF solvers. The thesis finally investigates solving a full AC-islanding problem, which can be considered as an extension of the transmission switching problem, using a standard MINLP solver and comparing the results to solutions obtained from approximations to the AC problem. Analysing in detail the results of the AC-islanding problem, alterations are made to the standard MINLP solver to allow better results to be obtained, all the while considering the trade-off between results and elapsed time.

abstract: For a (N+1)-bus power system, possibly 2N solutions exists. One of these solutions is known as the high-voltage (HV) solution or operable solution. The rest of the solutions are the low-voltage (LV), or large-angle, solutions. In this report, a recently developed non-iterative algorithm for solving the power- flow (PF) problem using the holomorphic embedding (HE) method is shown as being capable of finding the HV solution, while avoiding converging to LV solutions nearby which is a drawback to all other iterative solutions. The HE method provides a novel non-iterative procedure to solve the PF problems by eliminating the non-convergence and initial-estimate dependency issues appeared in the traditional iterative methods. The detailed implementation of the HE method is discussed in the report. While published work focuses mainly on finding the HV PF solution, modified holomorphically embedded formulations are proposed in this report to find the LV/large-angle solutions of the PF problem. It is theoretically proven that the proposed method is guaranteed to find a total number of 2N solutions to the PF problem and if no solution exists, the algorithm is guaranteed to indicate such by the oscillations in the maximal analytic continuation of the coefficients of the voltage power series obtained. After presenting the derivation of the LV/large-angle formulations for both PQ and PV buses, numerical tests on the five-, seven- and 14-bus systems are conducted to find all the solutions of the system of nonlinear PF equations for those systems using the proposed HE method. After completing the derivation to find all the PF solutions using the HE method, it is shown that the proposed HE method can be used to find only the of interest PF solutions (i.e. type-1 PF solutions with one positive real-part eigenvalue in the Jacobian matrix), with a proper algorithm developed. The closet unstable equilibrium point (UEP), one of the type-1 UEP’s, can be obtained by the proposed HE method with limited dynamic models included. The numerical performance as well as the robustness of the proposed HE method is investigated and presented by implementing the algorithm on the problematic cases and large-scale power system. / Dissertation/Thesis / Doctoral Dissertation Electrical Engineering 2015

abstract: Power flow calculation plays a significant role in power system studies and operation. To ensure the reliable prediction of system states during planning studies and in the operating environment, a reliable power flow algorithm is desired. However, the traditional power flow methods (such as the Gauss Seidel method and the Newton-Raphson method) are not guaranteed to obtain a converged solution when the system is heavily loaded. This thesis describes a novel non-iterative holomorphic embedding (HE) method to solve the power flow problem that eliminates the convergence issues and the uncertainty of the existence of the solution. It is guaranteed to find a converged solution if the solution exists, and will signal by an oscillation of the result if there is no solution exists. Furthermore, it does not require a guess of the initial voltage solution. By embedding the complex-valued parameter α into the voltage function, the power balance equations become holomorphic functions. Then the embedded voltage functions are expanded as a Maclaurin power series, V(α). The diagonal Padé approximant calculated from V(α) gives the maximal analytic continuation of V(α), and produces a reliable solution of voltages. The connection between mathematical theory and its application to power flow calculation is described in detail. With the existing bus-type-switching routine, the models of phase shifters and three-winding transformers are proposed to enable the HE algorithm to solve practical large-scale systems. Additionally, sparsity techniques are used to store the sparse bus admittance matrix. The modified HE algorithm is programmed in MATLAB. A study parameter β is introduced in the embedding formula βα + (1- β)α^2. By varying the value of β, numerical tests of different embedding formulae are conducted on the three-bus, IEEE 14-bus, 118-bus, 300-bus, and the ERCOT systems, and the numerical performance as a function of β is analyzed to determine the “best” embedding formula. The obtained power-flow solutions are validated using MATPOWER. / Dissertation/Thesis / Flow chart of the HE algorithm / Presentation for mater's thesis defense / Masters Thesis Electrical Engineering 2015

abstract: Recently, a novel non-iterative power flow (PF) method known as the Holomorphic Embedding Method (HEM) was applied to the power-flow problem. Its superiority over other traditional iterative methods such as Gauss-Seidel (GS), Newton-Raphson (NR), Fast Decoupled Load Flow (FDLF) and their variants is that it is theoretically guaranteed to find the operable solution, if one exists, and will unequivocally signal if no solution exists. However, while theoretical convergence is guaranteed by Stahl’s theorem, numerical convergence is not. Numerically, the HEM may require extended precision to converge, especially for heavily-loaded and ill-conditioned power system models. In light of the advantages and disadvantages of the HEM, this report focuses on three topics: 1. Exploring the effect of double and extended precision on the performance of HEM, 2. Investigating the performance of different embedding formulations of HEM, and 3. Estimating the saddle-node bifurcation point (SNBP) from HEM-based Thévenin-like networks using pseudo-measurements. The HEM algorithm consists of three distinct procedures that might accumulate roundoff error and cause precision loss during the calculations: the matrix equation solution calculation, the power series inversion calculation and the Padé approximant calculation. Numerical experiments have been performed to investigate which aspect of the HEM algorithm causes the most precision loss and needs extended precision. It is shown that extended precision must be used for the entire algorithm to improve numerical performance. A comparison of two common embedding formulations, a scalable formulation and a non-scalable formulation, is conducted and it is shown that these two formulations could have extremely different numerical properties on some power systems. The application of HEM to the SNBP estimation using local-measurements is explored. The maximum power transfer theorem (MPTT) obtained for nonlinear Thévenin-like networks is validated with high precision. Different numerical methods based on MPTT are investigated. Numerical results show that the MPTT method works reasonably well for weak buses in the system. The roots method, as an alternative, is also studied. It is shown to be less effective than the MPTT method but the roots of the Padé approximant can be used as a research tool for determining the effects of noisy measurements on the accuracy of SNBP prediction. / Dissertation/Thesis / Masters Thesis Electrical Engineering 2018

abstract: With the power system being increasingly operated near its limits, there is an increasing need for a power-flow (PF) solution devoid of convergence issues. Traditional iterative methods are extremely initial-estimate dependent and not guaranteed to converge to the required solution. Holomorphic Embedding (HE) is a novel non-iterative procedure for solving the PF problem. While the theory behind a restricted version of the method is well rooted in complex analysis, holomorphic functions and algebraic curves, the practical implementation of the method requires going beyond the published details and involves numerical issues related to Taylor's series expansion, Padé approximants, convolution and solving linear matrix equations. The HE power flow was developed by a non-electrical engineer with language that is foreign to most engineers. One purpose of this document to describe the approach using electric-power engineering parlance and provide an understanding rooted in electric power concepts. This understanding of the methodology is gained by applying the approach to a two-bus dc PF problem and then gradually from moving from this simple two-bus dc PF problem to the general ac PF case. Software to implement the HE method was developed using MATLAB and numerical tests were carried out on small and medium sized systems to validate the approach. Implementation of different analytic continuation techniques is included and their relevance in applications such as evaluating the voltage solution and estimating the bifurcation point (BP) is discussed. The ability of the HE method to trace the PV curve of the system is identified. / Dissertation/Thesis / Masters Thesis Electrical Engineering 2014