In power engineering, the powerflow study, or loadflow study, is a numerical analysis of the flow of electric power in an interconnected system. A powerflow study usually uses simplified notation such as a oneline diagram and perunit system, and focuses on various aspects of AC power parameters, such as voltages, voltage angles, real power and reactive power. It analyzes the power systems in normal steadystate operation.
Powerflow or loadflow studies are important for planning future expansion of power systems as well as in determining the best operation of existing systems. The principal information obtained from the powerflow study is the magnitude and phase angle of the voltage at each bus, and the real and reactive power flowing in each line.
Commercial power systems are usually too complex to allow for hand solution of the power flow. Special purpose network analyzers were built between 1929 and the early 1960s to provide laboratoryscale physical models of power systems. Largescale digital computers replaced the analog methods with numerical solutions.
In addition to a powerflow study, computer programs perform related calculations such as shortcircuit fault analysis, stability studies (transient & steadystate), unit commitment and economic dispatch.^{[1]} In particular, some programs use linear programming to find the optimal power flow, the conditions which give the lowest cost per kilowatt hour delivered.
A load flow study is especially valuable for a system with multiple load centers, such as a refinery complex. The power flow study is an analysis of the system’s capability to adequately supply the connected load. The total system losses, as well as individual line losses, also are tabulated. Transformer tap positions are selected to ensure the correct voltage at critical locations such as motor control centers. Performing a load flow study on an existing system provides insight and recommendations as to the system operation and optimization of control settings to obtain maximum capacity while minimizing the operating costs. The results of such an analysis are in terms of active power, reactive power, magnitude and phase angle.
The Open Energy Modelling Initiative promotes open source loadflow models, and other types of energy system models.
Contents

Model 1

Powerflow problem formulation 2

Newton–Raphson solution method 3

Other powerflow methods 4

References 5
Model
An alternating current powerflow model is a model used in electrical engineering to analyze power grids. It provides a nonlinear system which describes the energy flow through each transmission line. The problem is nonlinear because the power flow into load impedances is a function of the square of the applied voltages. Due to nonlinearity, in many cases the analysis of large network via AC powerflow model is not feasible, and a linear (but less accurate) DC powerflow model is used instead.
Usually analysis of a threephase system is simplified by assuming balanced loading of all three phases. Steadystate operation is assumed, with no transient changes in power flow or voltage due to load or generation changes. The system frequency is also assumed to be constant. A further simplification is to use the perunit system to represent all voltages, power flows, and impedances, scaling the actual target system values to some convenient base. A system oneline diagram is the basis to build a mathematical model of the generators, loads, buses, and transmission lines of the system, and their electrical impedances and ratings.
Powerflow problem formulation
The goal of a powerflow study is to obtain complete voltage angle and magnitude information for each bus in a power system for specified load and generator real power and voltage conditions.^{[2]} Once this information is known, real and reactive power flow on each branch as well as generator reactive power output can be analytically determined. Due to the nonlinear nature of this problem, numerical methods are employed to obtain a solution that is within an acceptable tolerance.
The solution to the powerflow problem begins with identifying the known and unknown variables in the system. The known and unknown variables are dependent on the type of bus. A bus without any generators connected to it is called a Load Bus. With one exception, a bus with at least one generator connected to it is called a Generator Bus. The exception is one arbitrarilyselected bus that has a generator. This bus is referred to as the slack bus.
In the powerflow problem, it is assumed that the real power P_{D} and reactive power Q_{D} at each Load Bus are known. For this reason, Load Buses are also known as PQ Buses. For Generator Buses, it is assumed that the real power generated P_{G} and the voltage magnitude V is known. For the Slack Bus, it is assumed that the voltage magnitude V and voltage phase Θ are known. Therefore, for each Load Bus, both the voltage magnitude and angle are unknown and must be solved for; for each Generator Bus, the voltage angle must be solved for; there are no variables that must be solved for the Slack Bus. In a system with N buses and R generators, there are then 2(N1)  (R1) unknowns.
In order to solve for the 2(N1)  (R1) unknowns, there must be 2(N1)  (R1) equations that do not introduce any new unknown variables. The possible equations to use are power balance equations, which can be written for real and reactive power for each bus. The real power balance equation is:

0 = P_{i} + \sum_{k=1}^N V_iV_k(G_{ik}\cos\theta_{ik}+B_{ik}\sin\theta_{ik})
where P_{i} is the net power injected at bus i, G_{ik} is the real part of the element in the bus admittance matrix Y_{BUS} corresponding to the ith row and kth column, B_{ik} is the imaginary part of the element in the Y_{BUS} corresponding to the ith row and kth column and \theta_{ik} is the difference in voltage angle between the ith and kth buses (\theta_{ik}=\delta_i\delta_k). The reactive power balance equation is:

0 = Q_{i} + \sum_{k=1}^N V_iV_k(G_{ik}\sin\theta_{ik}B_{ik}\cos\theta_{ik})
where Q_i is the net reactive power injected at bus i.
Equations included are the real and reactive power balance equations for each Load Bus and the real power balance equation for each Generator Bus. Only the real power balance equation is written for a Generator Bus because the net reactive power injected is not assumed to be known and therefore including the reactive power balance equation would result in an additional unknown variable. For similar reasons, there are no equations written for the Slack Bus.
In many transmission systems, the voltage angles \theta_{ik} are usually relatively small. There is thus a strong coupling between real power and voltage angle, and between reactive power and voltage magnitude, while the coupling between real power and voltage magnitude, as well as reactive power and voltage angle, is weak. As a result, real power is usually transmitted from the bus with higher voltage angle to the bus with lower voltage angle, and reactive power is usually transmitted from the bus with higher voltage magnitude to the bus with lower voltage magnitude. However, this approximation does not hold when the voltage angle is very large.^{[3]}
Newton–Raphson solution method
There are several different methods of solving the resulting nonlinear system of equations. The most popular is known as the Newton–Raphson method. This method begins with initial guesses of all unknown variables (voltage magnitude and angles at Load Buses and voltage angles at Generator Buses). Next, a Taylor Series is written, with the higher order terms ignored, for each of the power balance equations included in the system of equations . The result is a linear system of equations that can be expressed as:

\begin{bmatrix}\Delta \theta \\ \Delta V\end{bmatrix} = J^{1} \begin{bmatrix}\Delta P \\ \Delta Q \end{bmatrix}
where \Delta P and \Delta Q are called the mismatch equations:

\Delta P_i = P_i + \sum_{k=1}^N V_iV_k(G_{ik}\cos\theta_{ik}+B_{ik}\sin \theta_{ik})
\Delta Q_{i} = Q_{i} + \sum_{k=1}^N V_iV_k(G_{ik}\sin\theta_{ik}B_{ik}\cos\theta_{ik})
and J is a matrix of partial derivatives known as a Jacobian: J=\begin{bmatrix} \dfrac{\partial \Delta P}{\partial\theta} & \dfrac{\partial \Delta P}{\partial V} \\ \dfrac{\partial \Delta Q}{\partial \theta}& \dfrac{\partial \Delta Q}{\partial V}\end{bmatrix}.
The linearized system of equations is solved to determine the next guess (m + 1) of voltage magnitude and angles based on:

\theta^{m+1} = \theta^m + \Delta \theta\,

V^{m+1} = V^m + \Delta V\,
The process continues until a stopping condition is met. A common stopping condition is to terminate if the norm of the mismatch equations is below a specified tolerance.
A rough outline of solution of the powerflow problem is:

Make an initial guess of all unknown voltage magnitudes and angles. It is common to use a "flat start" in which all voltage angles are set to zero and all voltage magnitudes are set to 1.0 p.u.

Solve the power balance equations using the most recent voltage angle and magnitude values.

Linearize the system around the most recent voltage angle and magnitude values

Solve for the change in voltage angle and magnitude

Update the voltage magnitude and angles

Check the stopping conditions, if met then terminate, else go to step 2.
Other powerflow methods

Gauss–Seidel method: This is the earliest devised method. It shows slower rates of convergence compared to other iterative methods, but it uses very little memory and does not need to solve a matrix system.

Fastdecoupledloadflow method is a variation on NewtonRaphson that exploits the approximate decoupling of active and reactive flows in wellbehaved power networks, and additionally fixes the value of the Jacobian during the iteration in order to avoid costly matrix decompositions. Also referred to as "fixedslope, decoupled NR". Within the algorithm, the Jacobian matrix gets inverted only once, and there are three assumptions. Firstly, the conductance between the buses is zero. Secondly, the magnitude of the bus voltage is one per unit. Thirdly, the sine of phases between buses is zero. Fast decoupled load flow can return the answer within seconds whereas the Newton Raphson method takes much longer. This is useful for realtime management of power grids.

Holomorphic embedding load flow method: A recently developed method based on advanced techniques of complex analysis. It is direct and guarantees the calculation of the correct (operative) branch, out of the multiple solutions present in the power flow equations.
References

^ Low, S. H. (2013). "Convex relaxation of optimal power flow: A tutorial". 2013 IREP Symposium Bulk Power System Dynamics and Control  IX Optimization, Security and Control of the Emerging Power Grid. pp. 1–06.

^ Grainger, J.; Stevenson, W. (1994). Power System Analysis. New York: McGraw–Hill.

^ Andersson, G: Lectures on Modelling and Analysis of Electric Power Systems
This article was sourced from Creative Commons AttributionShareAlike License; additional terms may apply. World Heritage Encyclopedia content is assembled from numerous content providers, Open Access Publishing, and in compliance with The Fair Access to Science and Technology Research Act (FASTR), Wikimedia Foundation, Inc., Public Library of Science, The Encyclopedia of Life, Open Book Publishers (OBP), PubMed, U.S. National Library of Medicine, National Center for Biotechnology Information, U.S. National Library of Medicine, National Institutes of Health (NIH), U.S. Department of Health & Human Services, and USA.gov, which sources content from all federal, state, local, tribal, and territorial government publication portals (.gov, .mil, .edu). Funding for USA.gov and content contributors is made possible from the U.S. Congress, EGovernment Act of 2002.
Crowd sourced content that is contributed to World Heritage Encyclopedia is peer reviewed and edited by our editorial staff to ensure quality scholarly research articles.
By using this site, you agree to the Terms of Use and Privacy Policy. World Heritage Encyclopedia™ is a registered trademark of the World Public Library Association, a nonprofit organization.