No abstract

We study multistability in phase locked states in networks of phase oscillators under both Kuramoto dynamics and swing equation dynamics-a popular model for studying coarse-scale dynamics of an electrical AC power grid. We first establish the existence of geometrically frustrated states in such systems-where although a steady state flow pattern exists, no fixed point exists in the dynamical variables of phases due to geometrical constraints. We then describe the stable fixed points of the system with phase differences along each edge not exceeding π/2 in terms of cycle flows-constant flows along each simple cycle-as opposed to phase angles or flows. The cycle flow formalism allows us to compute tight upper and lower bounds to the number of fixed points in ring networks. We show that long elementary cycles, strong edge weights, and spatially homogeneous distribution of natural frequencies (for the Kuramoto model) or power injections (for the oscillator model for power grids) cause such networks to have more fixed points. We generalize some of these bounds to arbitrary planar topologies and derive scaling relations in the limit of large capacity and large cycle lengths, which we show to be quite accurate by numerical computation. Finally, we present an algorithm to compute all phase locked states-both stable and unstable-for planar networks.

We introduce the concept of network susceptibilities quantifying the response of the collective dynamics of a network to small parameter changes. We distinguish two types of susceptibilities: vertex susceptibilities and edge susceptibilities, measuring the responses due to changes in the properties of units and their interactions, respectively. We derive explicit forms of network susceptibilities for oscillator networks close to steady states and offer example applications for Kuramoto-type phase-oscillator models, power grid models, and generic flow models. Focusing on the role of the network topology implies that these ideas can be easily generalized to other types of networks, in particular those characterizing flow, transport, or spreading phenomena. The concept of network susceptibilities is broadly applicable and may straightforwardly be transferred to all settings where networks responses of the collective dynamics to topological changes are essential.

A new graph dual formalism is presented for the analysis of line outages in electricity networks. The dual formalism is based on a consideration of the flows around closed cycles in the network. After some exposition of the theory is presented, a new formula for the computation of Line Outage Distribution Factors (LODFs) is derived, which is not only computationally faster than existing methods, but also generalizes easily for multiple line outages and arbitrary changes to line series reactance. In addition, the dual formalism provides new physical insight for how the effects of line outages propagate through the network. For example, in a planar network a single line outage can be shown to induce monotonically decreasing flow changes, which are mathematically equivalent to an electrostatic dipole field.

Networks of phase oscillators are studied in various contexts, in particular in the modeling of the electric power grid. A functional grid corresponds to a stable steady state, such that any bifurcation can have catastrophic consequences up to a blackout. But also the existence of multiple steady states is undesirable, as it can lead to sudden transitions or circulatory flows. Despite the enormous practical importance there is still no general theory of the existence and uniqueness of steady states in such systems. Analytic results are mostly limited to grids without Ohmic losses. In this article, we introduce a method to systematically construct the solutions of the real power load-flow equations in the presence of Ohmic losses and explicitly compute them for tree and ring networks. We investigate different mechanisms leading to multistability and discuss the impact of Ohmic losses on the existence of solutions.The stable operation of the electric power grid relies on a precisely synchronized state of all generators and machines. All machines rotate at exactly the same frequency with fixed phase differences, leading to steady power flows throughout the grid. Whether such a steady state exists for a given network is of eminent practical importance. The loss of a steady state typically leads to power outages up to a complete blackout. But also the existence of multiple steady states is undesirable, as it can lead to sudden transitions, circulating flows and eventually also to power outages. Steady states are typically calculated numerically, but this approach gives only limited insight into the existence and (non-)uniqueness of steady states. Analytic results are available only for special network configuration, in particular for grids with negligible Ohmic losses or radial networks without any loops. In this article, we introduce a method to systematically construct the solutions of the real power loadflow equations in the presence of Ohmic losses. We calculate the steady states explicitly for elementary networks demonstrating different mechanisms leading to multistability. Our results also apply to models of coupled oscillators which are widely used in theoretical physics and mathematical biology.