Cycle flows and multistabilty in oscillatory networks: an overview

The functions of many networked systems in physics, biology or engineering rely on a coordinated or synchronized dynamics of its constituents. In power grids for example, all generators must synchronize and run at the same frequency and their phases need to appoximately lock to guarantee a steady power flow. Here, we analyze the existence and multitude of such phase-locked states. Focusing on edge and cycle flows instead of the nodal phases we derive rigorous results on the existence and number of such states. Generally, multiple phase-locked states coexist in networks with strong edges, long elementary cycles and a homogeneous distribution of natural frequencies or power injections, respectively. We offer an algorithm to systematically compute multiple phase- locked states and demonstrate some surprising dynamical consequences of multistability

Network recovery after massive failures

This paper addresses the problem of efficiently restoring sufficient resources in a communications network to support the demand of mission critical services after a large scale disruption. We give a formulation of the problem as an MILP and show that it is NP-hard. We propose a polynomial time heuristic, called Iterative Split and Prune (ISP) that decomposes the original problem recursively into smaller problems, until it determines the set of network components to be restored. We performed extensive simulations by varying the topologies, the demand intensity, the number of critical services, and the disruption model. Compared to several greedy approaches ISP performs better in terms of number of repaired components, and does not result in any demand loss. It performs very close to the optimal when the demand is low with respect to the supply network capacities, thanks to the ability of the algorithm to maximize sharing of repaired resources

Network Survivability Analysis: Coarse-Graining And Graph-Theoretic Strategies

In this dissertation, the interplay between geographic information about the network and the principal properties and structure of the underlying graph are used to quantify the structural and functional survivability of the network. This work focuses on the local aspect of survivability by studying the propagation of loss in the network as a function of the distance of the fault from a given origin-destination node pair

Network Flows

Not Availabl

Faster Algorithms for All-Pairs Bounded Min-Cuts

The All-Pairs Min-Cut problem (aka All-Pairs Max-Flow) asks to compute a minimum s-t cut (or just its value) for all pairs of vertices s,t. We study this problem in directed graphs with unit edge/vertex capacities (corresponding to edge/vertex connectivity). Our focus is on the k-bounded case, where the algorithm has to find all pairs with min-cut value less than k, and report only those. The most basic case k=1 is the Transitive Closure (TC) problem, which can be solved in graphs with n vertices and m edges in time O(mn) combinatorially, and in time O(n^{omega}) where omega<2.38 is the matrix-multiplication exponent. These time bounds are conjectured to be optimal. We present new algorithms and conditional lower bounds that advance the frontier for larger k, as follows: - A randomized algorithm for vertex capacities that runs in time {O}((nk)^{omega}). This is only a factor k^omega away from the TC bound, and nearly matches it for all k=n^{o(1)}. - Two deterministic algorithms for edge capacities (which is more general) that work in DAGs and further reports a minimum cut for each pair. The first algorithm is combinatorial (does not involve matrix multiplication) and runs in time {O}(2^{{O}(k^2)}* mn). The second algorithm can be faster on dense DAGs and runs in time {O}((k log n)^{4^{k+o(k)}}* n^{omega}). Previously, Georgiadis et al. [ICALP 2017], could match the TC bound (up to n^{o(1)} factors) only when k=2, and now our two algorithms match it for all k=o(sqrt{log n}) and k=o(log log n). - The first super-cubic lower bound of n^{omega-1-o(1)} k^2 time under the 4-Clique conjecture, which holds even in the simplest case of DAGs with unit vertex capacities. It improves on the previous (SETH-based) lower bounds even in the unbounded setting k=n. For combinatorial algorithms, our reduction implies an n^{2-o(1)} k^2 conditional lower bound. Thus, we identify new settings where the complexity of the problem is (conditionally) higher than that of TC. Our three sets of results are obtained via different techniques. The first one adapts the network coding method of Cheung, Lau, and Leung [SICOMP 2013] to vertex-capacitated digraphs. The second set exploits new insights on the structure of latest cuts together with suitable algebraic tools. The lower bounds arise from a novel reduction of a different structure than the SETH-based constructions

Designing and Expanding Electrical Networks – Complexity and Combinatorial Algorithms

The transition from conventional to renewable power generation has a large impact on when and where electricity is generated. To deal with this change the electric transmission network needs to be adapted and expanded. Expanding the network has two benefits. Electricity can be generated at locations with high renewable energy potentials and then transmitted to the consumers via the transmission network. Without the expansion the existing transmission network may be unable to cope with the transmission needs, thus requiring power generation at locations closer to the energy demand, but at less well-suited locations. Second, renewable energy generation (e.g., from wind or solar irradiation) is typically volatile. Having strong interconnections between regions within a large geographical area allows to the smooth the generation and demand over that area. This smoothing makes them more predictable and the volatility of the generation easier to handle. In this thesis we consider problems that arise when designing and expanding electric transmission networks. As the first step we formalize them such that we have a precise mathematical problem formulation. Afterwards, we pursue two goals: first, improve the theoretical understanding of these problems by determining their computational complexity under various restrictions, and second, develop algorithms that can solve these problems. A basic formulation of the expansion planning problem models the network as a graph and potential new transmission lines as edges that may be added to the graph. We formalize this formulation as the problems Flow Expansion and Electrical Flow Expansion, which differ in the flow model (graph-theoretical vs. electrical flow). We prove that in general the decision variants of these problems are $\mathcal{NP}$-complete, even if the network structure is already very simple, e.g., a star. For certain restrictions, we give polynomial-time algorithms as well. Our results delineate the boundary between the $\mathcal{NP}$-complete cases and the cases that can be solved in polynomial time. The basic expansion planning problems mentioned above ignore that real transmission networks should still be able to operate if a small part of the transmission equipment fails. We employ a criticality measure from the literature, which measures the dynamic effects of the failure of a single transmission line on the whole transmission network. In a first step, we compare this criticality measure to the well-used $N-1$ criterion. Moreover, we formulate this criticality measure as a set of linear inequalities, which may be added to any formulation of a network design problem as a mathematical program. To exemplify this usage, we introduce the criticality criterion in two transmission network expansion planning problems, which can be formulated as mixed-integer linear programs (MILPs). We then evaluate the performance of solving the MILPs. Finally, we develop a greedy heuristic for one of the two problems, and compare its performance to solving the MILP. Microgrids play an important role in the electrification of rural areas. We formalize the design of the cable layout of a microgrid as a geometric optimization problem, which we call Microgrid Cable Layout. A key difference to the network design problems above is that there is no graph with candidate edges given. Instead, edges and new vertices may be placed anywhere in the plane. We present a hybrid genetic algorithm for Microgrid Cable Layout and evaluate it on a set of benchmark instances, which include a real microgrid in the Democratic Republic of the Congo. Finally, instead of expanding electrical networks one may place electric equipment such as FACTS (flexible AC transmission system). These influence the properties of the transmission lines such that the network can be used more efficiently. We apply a model of FACTS from the literature and study the problem whether a given network with given positions and properties of the FACTS admits an electrical flow provided that FACTS are set appropriately. We call such a flow a FACTS flow. In this thesis we prove that in general it is $\mathcal{NP}$-complete to determine whether a network admits a FACTS flow, and we present polynomial-time algorithms for two restricted cases

All-pairs min-cut in sparse networks

Algorithms are presented for the all-pairs min-cut problem in bounded tree-width, planar and sparse networks. The approach used is to preprocess the input $n$-vertex network so that, afterwards, the value of a min-cut between any two vertices can be efficiently computed. A tradeoff is shown between the preprocessing time and the time taken to compute min-cuts subsequently. In particular, after an $O(n\log n)$ preprocessing of a bounded tree-width network, it is possible to find the value of a min-cut between any two vertices in constant time. This implies that for such networks the all-pairs min-cut problem can be solved in time $O(n^2)$. This algorithm is used in conjunction with a graph decomposition technique of Frederickson to obtain algorithms for sparse and planar networks. The running times depend upon a topological property, $\gamma$, of the input network. The parameter $\gamma$ varies between 1 and $\Theta(n)$; the algorithms perform well when $\gamma = o(n)$. The value of a min-cut can be found in time $O(n + \gamma^2 \log \gamma)$ and all-pairs min-cut can be solved in time $O(n^2 + \gamma^4 \log \gamma)$ for sparse networks. The corresponding running times4 for planar networks are $O(n+\gamma \log \gamma)$ and $O(n^2 + \gamma^3 \log \gamma)$, respectively. The latter bounds depend on a result of independent interest: outerplanar networks have small ``mimicking'' networks which are also outerplanar