Calculating Nodal Voltages Using the Admittance Matrix Spectrum of an Electrical Network

Calculating nodal voltages and branch current flows in a meshed network is fundamental to electrical engineering. This work demonstrates how such calculations can be performed using the eigenvalues and eigenvectors of the Laplacian matrix which describes the connectivity of the electrical network. These insights should permit the functioning of electrical networks to be understood in the context of spectral analysis.


Introduction
Electrical power system calculations rely heavily on the bus admittance matrix, Y bus , which is a Laplacian matrix weighted by the complex-valued admittance of each branch in the network.It is well established that the eigenvalues and eigenvectors (deemed the spectrum) of a Laplacian matrix encode meaningful information about a network's structure [1].Recent work in [2,3] indicates that, in electrical networks, this spectrum can be directly related to nodal voltages and branch current flows.The purpose of the present paper is to clarify the derivations provided in [3].The scope of the present work is narrowly theoretical: Linear algebra is used to articulate the correct relationship between the variables treated in [3].
Notwithstanding these modest ambitions, a key motivation for the present work is to begin to link power flow analysis with the mature literature [4] on spectral graph theory.Extant efforts to apply spectral graph theory to electrical networks are scarce, but include [5,6].The use of graph theory more generally in this role is reviewed in [7,8].Notably, simplistic topological approaches do not properly account for the physical realities of electrical power flow, and can thereby fail to identify the critical components in an electrical network [9][10][11].The present work seeks to articulate one particular linkage between spectral graph theory and circuit theory, which may offer new ways to understand how power flows in meshed electrical networks.
The rest of this paper is organized as follows: In Section 2 we establish the necessary preliminaries, including electrical flow basics and notation.The main results are presented in Section 3.

Electrical Flow Basics and Notation
Ohm's law linearly relates the current flowing through an edge in a circuit with the voltage difference between the nodes that the edge connects.Specifically, I kj = ∆V kj Z kj , and ∑ N j=1 I kj = F k , k, j = 1, 2, . . ., N, where I kj is the current passing from the k-th node to the j-th in a (typically sparsely connected) network of N nodes, ∆V kj = V k − V j is the voltage difference between the k-th node and the j-th, Z kj = Z jk is branch impedance and F k are complex-valued net current injections or withdrawals.
From the above notation we arrive easily at ∑ N j=1 In the article we will denote with δ ij for the Kronecker delta, i.e., δ ii = 1 and δ ij = 0 for i = j.With ū we will denote the complex conjugate of u, and with T , * the conjugate transpose, and conjugate transpose tensor respectively.

An Exemplary Electrical Network
To provide some context, a nation-spanning electrical power system is shown in Figure 1.This diagram of the nesta_case2224_edin test system [12] was created using the techniques described in [13], which uses electrical distances measures, rather than physical geography, to positions nodes.Note the relative spareseness of its connective struture, and how lower nominal voltage levels (< 143 kV) correspond to more tree-like structures.This network of 2224 nodes supplies a total load of up to 60 GW, supplied from 378 different generating sites.

Derivations
In this section, first we rewrite (1) in matrix form and define the relevant Laplacian matrix.Then we provide a formula which explicitly relates the voltage differences to the eigenvalues and eigenvectors of the Laplacian matrix for meshed electrical networks.We can state now the following theorem.
Theorem 1.Consider an electrical network with branch currents I kj , ∀ k, j = 1, 2, . . ., N passing from node k to node j, a complex impedance describing each branch Z kj = Z jk , and F k being the complex-valued net current flow at each bus with ∑ N k=1 F k = 0. Then the voltage difference ∆V mn between two arbitrary nodes m and n is given by: where λ k , k = 2, 3, . . ., N are the non-zero eigenvalues of the G matrix (equivalent to the Y bus matrix in the power systems context) which describes the connectivity of the electrical network: and u 1k u 2k . . .u Nk T is an eigenvector of the eigenvalue λ k .
, or, equivalently, by setting V = V i i=1,2,...,N , F = F i i=1,2,...,N , and , We arrive at GV = F.We observe that if G kj , k, j = 1, 2, . . ., N is an element of G, then for k = j, G kk = ∑ N j=1 1 . Hence G is given by (3).We will refer to G as the Laplacian matrix.Note that the rows of G sum to zero, i.e., the matrix has the zero eigenvalue (see [3,14]).The algebraic multiplicity of the zero eigenvalue in the Laplacian is the number of connected components in the network.In the power systems case we deal with only one network which means the algebraic multiplicity of the zero eigenvalue is one.Since the matrix G is symmetric it can be written in the following form: G = PDP * , where P = [u kj ] j=1,2,...,N k=1,2,...,N , P * is the conjugate transpose of P such that PP * is the N × N identity matrix and D is the diagonal matrix D = diag{0, λ 2 , λ 3 , . . ., λ N }.By applying the above expression into the system we get: and since P * is the inverse of P we have: .
Let 1 N be a column vector that contains exactly N 1's.From the fact that every row of G sums to zero we have the eigenspace of the zero eigenvalue.Indeed G • 1 N = 0 • 1 N which means that < 1 N > is the eigenspace of the zero eigenvalue.Hence there exist c ∈ C such that From ( 4) u k1 = c, ∀k = 1, 2, . . ., N, or, equivalently, ūk1 = c, ∀k = 1, 2, . . ., N. In addition, By ignoring the first row of each column of the above expression we get: Which can be rewritten in the following form: , we have Let V m , V n be two arbitrary nodal voltages, i.e., Then, the difference between them is given by

Conclusions
This work has clarified the relationship between the admittance matrix spectrum, the current inflows & withdrawals prevailing in an electrical network and the resulting nodal voltage profile.Applying these spectral relationships to practical electrical engineering problems is left to future work.

Figure 1 .
Figure 1.This diagram shows the nesta_case2224_edin test power system