Showing 1–10 of 10 results for author: Smith, K D

Contractivity of the Method of Successive Approximations for Optimal Control

Authors: Kevin D. Smith, Francesco Bullo

Abstract: Strongly contracting dynamical systems have numerous properties (e.g., incremental ISS), find widespread applications (e.g., in controls and learning), and their study is receiving increasing attention. This work starts with the simple observation that, given a strongly contracting system, its adjoint dynamical system is also strongly contracting, with the same rate, with respect to the dual norm,… ▽ More Strongly contracting dynamical systems have numerous properties (e.g., incremental ISS), find widespread applications (e.g., in controls and learning), and their study is receiving increasing attention. This work starts with the simple observation that, given a strongly contracting system, its adjoint dynamical system is also strongly contracting, with the same rate, with respect to the dual norm, under time reversal. As main implication of this dual contractivity, we show that the classic Method of Successive Approximations (MSA), an indirect method in optimal control, is a contraction mapping for short optimization intervals or large contraction rates. Consequently, we establish new convergence conditions for the MSA algorithm, which further imply uniqueness of the optimal control and sufficiency of Pontryagin's minimum principle under additional assumptions. △ Less

Submitted 16 December, 2022; v1 submitted 20 September, 2022; originally announced September 2022.

The "SPectrogram Analysis and Cataloguing Environment" (SPACE) Labelling Tool

Authors: C. K. Louis, C. M. Jackman, S. W. Mangham, K. D. Smith, E. P. O'Dwyer, A. Empey, B. Cecconi, P. Zarka, S. Maloney

Abstract: The SPectrogram Analysis and Cataloguing Environment (SPACE) tool is an interactive python tool designed to label radio emission features of interest in a time-frequency map (called 'dynamic spectrum'). The program uses Matplotlib's Polygon Selector widget to allow a user to select and edit an undefined number of vertices on top of the dynamic spectrum before closing the shape (polygon). Multiple… ▽ More The SPectrogram Analysis and Cataloguing Environment (SPACE) tool is an interactive python tool designed to label radio emission features of interest in a time-frequency map (called 'dynamic spectrum'). The program uses Matplotlib's Polygon Selector widget to allow a user to select and edit an undefined number of vertices on top of the dynamic spectrum before closing the shape (polygon). Multiple polygons may be drawn on any spectrum, and the feature name along with the coordinates for each polygon vertex are saved into a '.json' file as per the 'Time-Frequency Catalogue' (TFCat) format along with other data such as the feature id, observer name, and data units. This paper describes the first official stable release (version 2.0) of the tool. △ Less

Submitted 25 July, 2022; originally announced July 2022.

Comments: 9 pages, 1 Figure

Dual Seminorms, Ergodic Coefficients and Semicontraction Theory

Authors: Giulia De Pasquale, Kevin D. Smith, Francesco Bullo, Maria Elena Valcher

Abstract: Dynamical systems that are contracting on a subspace are said to be semicontracting. Semicontraction theory is a useful tool in the study of consensus algorithms and dynamical flow systems such as Markov chains. To develop a comprehensive theory of semicontracting systems, we investigate seminorms on vector spaces and define two canonical notions: projection and distance semi-norms. We show that t… ▽ More Dynamical systems that are contracting on a subspace are said to be semicontracting. Semicontraction theory is a useful tool in the study of consensus algorithms and dynamical flow systems such as Markov chains. To develop a comprehensive theory of semicontracting systems, we investigate seminorms on vector spaces and define two canonical notions: projection and distance semi-norms. We show that the well-known lp ergodic coefficients are induced matrix seminorms and play a central role in stability problems. In particular, we formulate a duality theorem that explains why the Markov-Dobrushin coefficient is the rate of contraction for both averaging and conservation flows in discrete time. Moreover, we obtain parallel results for induced matrix log seminorms. Finally, we propose comprehensive theorems for strong semicontractivity of linear and non-linear time-varying dynamical systems with invariance and conservation properties both in discrete and continuous time. △ Less

Submitted 20 December, 2022; v1 submitted 9 January, 2022; originally announced January 2022.

Convex Optimization of the Basic Reproduction Number

Authors: Kevin D. Smith, Francesco Bullo

Abstract: The basic reproduction number $R_0$ is a fundamental quantity in epidemiological modeling, reflecting the typical number of secondary infections that arise from a single infected individual. While $R_0$ is widely known to scientists, policymakers, and the general public, it has received comparatively little attention in the controls community. This note provides two novel characterizations of… ▽ More The basic reproduction number $R_0$ is a fundamental quantity in epidemiological modeling, reflecting the typical number of secondary infections that arise from a single infected individual. While $R_0$ is widely known to scientists, policymakers, and the general public, it has received comparatively little attention in the controls community. This note provides two novel characterizations of $R_0$: a stability characterization and a geometric program characterization. The geometric program characterization allows us to write $R_0$-constrained and budget-constrained optimal resource allocation problems as geometric programs, which are easily transformed into convex optimization problems. We apply these programs to allocating vaccines and antidotes in numerical examples, finding that targeting $R_0$ instead of the spectral abscissa of the Jacobian matrix (a common target in the controls literature) leads to qualitatively different solutions. △ Less

Submitted 2 September, 2022; v1 submitted 15 September, 2021; originally announced September 2021.

Complex Fault Geometry of the 2020 MWW6.5 Monte Cristo Range, Nevada Earthquake Sequence

Authors: Christine J. Ruhl, Emily A. Morton, Jayne M. Bormann, Rachel Hatch-Ibarra, Gene Ichinose, Kenneth D. Smith

Abstract: On 15 May 2020 an MWW 6.5 earthquake occurred beneath the Monte Cristo Range in the Mina Deflection region of western Nevada. Rapid deployment of eight temporary seismic stations enables detailed analysis of its productive and slowly decaying aftershock sequence (p=0.8) which included ~18,000 autodetected events in 3.5 months. Double-difference, waveform-based relative relocation of 16,714 earthqu… ▽ More On 15 May 2020 an MWW 6.5 earthquake occurred beneath the Monte Cristo Range in the Mina Deflection region of western Nevada. Rapid deployment of eight temporary seismic stations enables detailed analysis of its productive and slowly decaying aftershock sequence (p=0.8) which included ~18,000 autodetected events in 3.5 months. Double-difference, waveform-based relative relocation of 16,714 earthquakes reveals a complex network of faults, many of which cross the inferred 35-km long east-northeast-striking, left-lateral mainshock rupture. Seismicity aligns with left-lateral, right-lateral, and normal mechanism moment tensors of 128 of the largest earthquakes. The mainshock occurred near the middle of the aftershock zone at the intersection of two distinct zones of seismicity. In the western section, numerous subparallel, shallow, north-northeast-striking faults form a broad flower-structure-like fault mesh that coalesces at depth into a near-vertical, left-lateral fault. We infer the near-vertical fault to be a region of significant slip in the mainshock and an eastward extension of the left-lateral Candelaria fault. Near the mainshock hypocenter, seismicity occurs on a northeast-striking, west-dipping structure which extends north from the Eastern Columbus Salt Marsh normal fault. Together, these two intersecting structures bound the Columbus Salt Marsh tectonic basin. East of this intersection and the mainshock hypocenter, seismicity occurs in a narrow, near-vertical, east-northeast-striking fault zone through to its eastern terminus. At the eastern end, the aftershock zone broadens and extends northwest towards the southern extension of the northwest-striking, right-lateral Petrified Springs fault system. The eastern section hosts significantly fewer aftershocks than the western section, but has more moment release. △ Less

Submitted 28 February, 2021; v1 submitted 16 February, 2021; originally announced February 2021.

A Tutorial on Multivariate $k$-Statistics and their Computation

Authors: Kevin D. Smith

Abstract: This document aims to provide an accessible tutorial on the unbiased estimation of multivariate cumulants, using $k$-statistics. We offer an explicit and general formula for multivariate $k$-statistics of arbitrary order. We also prove that the $k$-statistics are unbiased, using Möbius inversion and rudimentary combinatorics. Many detailed examples are considered throughout the paper. We conclude… ▽ More This document aims to provide an accessible tutorial on the unbiased estimation of multivariate cumulants, using $k$-statistics. We offer an explicit and general formula for multivariate $k$-statistics of arbitrary order. We also prove that the $k$-statistics are unbiased, using Möbius inversion and rudimentary combinatorics. Many detailed examples are considered throughout the paper. We conclude with a discussion of $k$-statistics computation, including the challenge of time complexity, and we examine a couple of possible avenues to improve the efficiency of this computation. The purpose of this document is threefold: to provide a clear introduction to $k$-statistics without relying on specialized tools like the umbral calculus; to construct an explicit formula for $k$-statistics that might facilitate future approximations and faster algorithms; and to serve as a companion paper to our Python library PyMoments, which implements this formula. △ Less

Submitted 17 May, 2020; originally announced May 2020.

Topology Inference with Multivariate Cumulants: The Möbius Inference Algorithm

Authors: Kevin D. Smith, Saber Jafarpour, Ananthram Swami, Francesco Bullo

Abstract: Many tasks regarding the monitoring, management, and design of communication networks rely on knowledge of the routing topology. However, the standard approach to topology mapping--namely, active probing with traceroutes--relies on cooperation from increasingly non-cooperative routers, leading to missing information. Network tomography, which uses end-to-end measurements of additive link metrics (… ▽ More Many tasks regarding the monitoring, management, and design of communication networks rely on knowledge of the routing topology. However, the standard approach to topology mapping--namely, active probing with traceroutes--relies on cooperation from increasingly non-cooperative routers, leading to missing information. Network tomography, which uses end-to-end measurements of additive link metrics (like delays or log packet loss rates) across monitor paths, is a possible remedy. Network tomography does not require that routers cooperate with traceroute probes, and it has already been used to infer the structure of multicast trees. This paper goes a step further. We provide a tomographic method to infer the underlying routing topology of an arbitrary set of monitor paths using the joint distribution of end-to-end measurements, without making any assumptions on routing behavior. Our approach, called the Möbius Inference Algorithm (MIA), uses cumulants of this distribution to quantify high-order interactions among monitor paths, and it applies Möbius inversion to "disentangle" these interactions. In addition to MIA, we provide a more practical variant called Sparse Möbius Inference, which uses various sparsity heuristics to reduce the number and order of cumulants required to be estimated. We show the viability of our approach using synthetic case studies based on real-world ISP topologies. △ Less

Submitted 18 June, 2021; v1 submitted 16 May, 2020; originally announced May 2020.

Transient Stability of Droop-Controlled Inverter Networks with Operating Constraints

Authors: Kevin D. Smith, Saber Jafarpour, Francesco Bullo

Abstract: Due to the rise of distributed energy resources, the control of networks of grid-forming inverters is now a pressing issue for power system operation. Droop control is a popular control strategy in the literature for frequency control of these inverters. In this paper, we analyze transient stability in droop-controlled inverter networks that are subject to multiple operating constraints. Using a p… ▽ More Due to the rise of distributed energy resources, the control of networks of grid-forming inverters is now a pressing issue for power system operation. Droop control is a popular control strategy in the literature for frequency control of these inverters. In this paper, we analyze transient stability in droop-controlled inverter networks that are subject to multiple operating constraints. Using a physically-meaningful Lyapunov-like function, we provide two sets of criteria (one mathematical and one computational) to certify that a post-fault trajectory achieves frequency synchronization while respecting operating constraints. We show how to obtain less-conservative transient stability conditions by incorporating information from loop flows, i.e., net flows of active power around cycles in the network. Finally, we use these conditions to quantify the scale of parameter disturbances to which the network is robust. We illustrate our results with numerical case studies of the IEEE 24-bus system. △ Less

Submitted 6 September, 2020; v1 submitted 11 July, 2019; originally announced July 2019.

Flow and Elastic Networks on the $n$-torus: Geometry, Analysis, and Computation

Authors: Saber Jafarpour, Elizabeth Y. Huang, Kevin D. Smith, Francesco Bullo

Abstract: Networks with phase-valued nodal variables are central in modeling several important societal and physical systems, including power grids, biological systems, and coupled oscillator networks. One of the distinctive features of phase-valued networks is the existence of multiple operating conditions corresponding to critical points of an energy function or feasible flows of a balance equation. For n… ▽ More Networks with phase-valued nodal variables are central in modeling several important societal and physical systems, including power grids, biological systems, and coupled oscillator networks. One of the distinctive features of phase-valued networks is the existence of multiple operating conditions corresponding to critical points of an energy function or feasible flows of a balance equation. For networks with phase-valued states, it is not yet fully understood how many operating conditions exist, how to characterize them, and how to compute them efficiently. A deeper understanding of feasible operating conditions, including their dependence upon network structures, may lead to more reliable and efficient network systems. This paper introduces flow and elastic network problems on the $n$-torus and provides a rigorous and comprehensive framework for their study. Based on a monotonicity assumption, this framework localizes the solutions, bounds their number, and leads to an algorithm to compute them. Our analysis is based on a novel winding partition of the $n$-torus into winding cells, induced by Kirchhoff's angle law for undirected graphs. The winding partition has several useful properties, including notably that, each winding cell contains at most one solution. The proposed algorithm is based on a novel contraction mapping and is guaranteed to compute all solutions. Finally, we apply our results to numerically study the active power flow equations in several test cases and estimate power capacity and congestion of a power network. △ Less

Submitted 12 September, 2020; v1 submitted 30 January, 2019; originally announced January 2019.

Switchable collective pinning of flux quanta using magnetic vortex arrays

Authors: J. E. Villegas, K. D. Smith, Lei Huang, Yimei Zhu, R. Morales, Ivan K. Schuller

Abstract: We constructed a superconducting/ferromagnetic hybrid system in which the ordering of the pinning potential landscape for flux quanta can be manipulated. Flux pinning is induced by an array of magnetic nanodots in the magnetic vortex state, and controlled by the magnetic history. This allows switching on and off the collective pinning of the flux-lattice. In addition, we observed field-induced s… ▽ More We constructed a superconducting/ferromagnetic hybrid system in which the ordering of the pinning potential landscape for flux quanta can be manipulated. Flux pinning is induced by an array of magnetic nanodots in the magnetic vortex state, and controlled by the magnetic history. This allows switching on and off the collective pinning of the flux-lattice. In addition, we observed field-induced superconductivity that originates from the annihilation of flux quanta induced by the stray fields from the magnetic vortices. △ Less

Submitted 31 March, 2008; originally announced March 2008.

Comments: PDF file 18 pages including 5 figures, accepted for publication in Phys. Rev. B

Journal ref: Phys. Rev. B 77, 134510 (2008)