Search

DC microgrids are prone to small-signal instabilities due to the presence of tightly regulated loads. This paper develops a decentralized stability certificate which is capable of certifying the small-signal stability of an islanded DC network containing such loads. Utilizing a novel homotopy approach, the proposed standards ensure that no system eigenmodes are able to cross into the unstable right half plane for a continuous range of controller gain levels. The resulting "standards" can be applied to variety of grid components which meet the specified, but non-unique, criteria. These standards thus take a step towards offering plug-and-play operability of DC microgrids. The proposed theorems are explicitly illustrated and numerically validated on multiple DC microgrid test-cases containing both buck and boost converter dynamics.

Electric power grids regularly experience uncertain fluctuations from load demands and renewables, which poses a risk of violating operational limits designed to safeguard the system. In this paper, we consider the robust AC OPF problem that minimizes the generation cost while requiring a certain level of system security in the presence of uncertainty. The robust AC OPF problem requires that the system satisfy operational limits for all uncertainty realizations within a specified uncertainty set. Guaranteeing robustness is particularly challenging due to the non-convex, nonlinear AC power flow equations, which may not always have a solution. In this work, we extend a previously developed convex restriction to a robust convex restriction, which is a convex inner approximation of the non-convex feasible region of the AC OPF problem that accounts for uncertainty in the power injections. We then use the robust convex restriction in an algorithm that obtains robust solutions to AC OPF problems by solving a sequence of convex optimization problems. We demonstrate our algorithm and its ability to control robustness versus operating cost trade-offs using PGLib test cases.

We present an algorithm for robust model predictive control with consideration of uncertainty and safety constraints. Our framework considers a nonlinear dynamical system subject to disturbances from an unknown but bounded uncertainty set. By viewing the system as a fixed point of an operator acting over trajectories, we propose a convex condition on control actions that guarantee safety against the uncertainty set. The proposed condition guarantees that all realizations of the state trajectories satisfy safety constraints. Our algorithm solves a sequence of convex quadratic constrained optimization problems of size n*N, where n is the number of states, and N is the prediction horizon in the model predictive control problem. Compared to existing methods, our approach solves convex problems while guaranteeing that all realizations of uncertainty set do not violate safety constraints. Moreover, we consider the implicit time-discretization of system dynamics to increase the prediction horizon and enhance computational accuracy. Numerical simulations for vehicle navigation demonstrate the effectiveness of our approach.

Smart buildings have great potential for shaping an energy-efficient, sustainable, and more economic future for our planet as buildings account for approximately 40% of the global energy consumption. Future of the smart buildings lies in using sensory data for adaptive decision making and control that is currently gloomed by the key challenge of learning a good control policy in a short period of time in an online and continuing fashion. To tackle this challenge, an event-triggered -- as opposed to classic time-triggered -- paradigm, is proposed in which learning and control decisions are made when events occur and enough information is collected. Events are characterized by certain design conditions and they occur when the conditions are met, for instance, when a certain state threshold is reached. By systematically adjusting the time of learning and control decisions, the proposed framework can potentially reduce the variance in learning, and consequently, improve the control process. We formulate the micro-climate control problem based on semi-Markov decision processes that allow for variable-time state transitions and decision making. Using extended policy gradient theorems and temporal difference methods in a reinforcement learning set-up, we propose two learning algorithms for event-triggered control of micro-climate in buildings. We show the efficacy of our proposed approach via designing a smart learning thermostat that simultaneously optimizes energy consumption and occupants' comfort in a test building.

This paper presents a convex sufficient condition for solving a system of nonlinear equations under parametric changes and proposes a sequential convex optimization method for solving robust optimization problems with nonlinear equality constraints. By bounding the nonlinearity with concave envelopes and using Brouwer's fixed point theorem, the sufficient condition is expressed in terms of closed-form convex inequality constraints. We extend the result to provide a convex sufficient condition for feasibility under bounded uncertainty. Using these conditions, a non-convex optimization problem can be solved as a sequence of convex optimization problems, with feasibility and robustness guarantees. We present a detailed analysis of the performance and complexity of the proposed condition. The examples in polynomial optimization and nonlinear network are provided to illustrate the proposed method.

Nonconvexity induced by the nonlinear AC power flow equations challenges solution algorithms for AC optimal power flow (OPF) problems. While significant research efforts have focused on reliably computing high-quality OPF solutions, it is not always clear that there exists a feasible path to reach the desired operating point. Transitioning between operating points while avoiding constraint violations can be challenging since the feasible space of the OPF problem is nonconvex and potentially disconnected. To address this problem, we propose an algorithm that computes a provably feasible path from an initial operating point to a desired operating point. Given an initial feasible point, the algorithm solves a sequence of convex quadratically constrained optimization problems over conservative convex inner approximations of the OPF feasible space. In each iteration, we obtain a new, improved operating point and a feasible transition from the operating point in the previous iteration. In addition to computing a feasible path to a known desired operating point, this algorithm can also be used to improve the operating point locally. Extensive numerical studies on a variety of test cases demonstrate the algorithm and the ability to arrive at a high-quality solution in few iterations.

This paper develops a systematic framework for analyzing how low frequency forced oscillations propagate in electric power systems. Using this framework, the paper shows how to mathematically justify the so-called Dissipating Energy Flow (DEF) forced oscillation source location technique. The DEF's specific deficiencies are pinpointed, and its underlying energy function is analyzed via incremental passivity theory. This analysis is then used to prove that there exists no passivity transformation (i.e. quadratic energy function) which can simultaneously render all components of a lossy classical power system passive. The paper goes on to develop a simulation-free algorithm for predicting the performance of the DEF method in a generalized power system, and it analyzes the passivity of three non-classical load and generation components. The proposed propagation framework and performance algorithm are both tested and illustrated on the IEEE 39-bus New England system and the WECC 179-bus system., Comment: 14 pages, re-submitted to IEEE TPWRS

Energy storage has been proven to yield positive effects on planning, operation and control of electric grids. It has become a crucial task to properly model the energy storage systems (ESS) under the framework of grid optimization on transmission and distribution networks including microgrids. This paper presents a review on mathematical models and test cases of ESSs used for grid optimization studies, where the network constraints of power systems are included. The existing ESS models are mainly classified into two categories; linear and nonlinear models. The two main categories are further divided into several subcategories respectively; such as mixed integer linear and convex nonlinear subcategories. Based on the review and discussions, this paper aims at providing suggestions for choosing proper ESS models for specific grid optimization studies considering the chosen power network model., Comment: Accepted for publication ahead of the PES General Meeting 2019

The small-signal stability is an integral part of the power system security analysis. The introduction of renewable source related uncertainties is making the stability assessment difficult as the equilibrium point is varying rapidly. This paper focuses on the Differential Algebraic Equation (DAE) formulation of power systems and bridges the gap between the conventional reduced system and the original one using logarithmic norm. We propose a sufficient condition for stability using Bilinear Matrix Inequality and its inner approximation as Linear Matrix Inequality. Another contribution is the construction of robust stability regions in state-space in contrast to most existing approaches trying same in the parameter space. Performance evaluation of the sufficiency condition and its inner approximation has been given with the necessary and sufficient condition for a small-scale test case. The paper provides a necessary base to develop tractable construction techniques for the robust stability region of power systems.

Despite wide-scale deployment of phasor measurement unit technology, locating the sources of low frequency forced oscillations in power systems is still an open research topic. The dissipating energy flow method is one source location technique which has performed remarkably well in both simulation and real time application at ISO New England. The method has several deficiencies, though, which are still poorly understood. This paper borrows the concepts of passivity and positive realness from the controls literature in order to interpret the dissipating energy flow method, pinpoint the reasons for its deficiencies, and set up a framework for improving the method. The theorems presented in this paper are then tested via simulation on a simple infinite bus power system model., Comment: Submitted to PESGM19

This paper presents a novel scalable framework to solve the optimization of a nonlinear system with differential algebraic equation (DAE) constraints that enforce the asymptotic stability of the underlying dynamic model with respect to certain disturbances. Existing solution approaches to analogous DAE-constrained problems are based on discretization of DAE system into a large set of nonlinear algebraic equations representing the time-marching schemes. These approaches are not scalable to large size models. The proposed framework, based on LaSalle's invariance principle, uses convex Lyapunov functions to develop a novel stability certificate which consists of a limited number of algebraic constraints. We develop specific algorithms for two major types of nonlinearities, namely Lur'e, and quasi-polynomial systems. Quadratic and convex-sum-of-square Lyapunov functions are constructed for the Lur'e-type and quasi-polynomial systems respectively. A numerical experiment is performed on a 3-generator power network to obtain a solution for transient-stability-constrained optimal power flow., Comment: 10 pages, 7 figures

Since forced oscillations are exogenous to dynamic power system models, the models by themselves cannot predict when or where a forced oscillation will occur. Locating the sources of these oscillations, therefore, is a challenging problem which requires analytical methods capable of using real time power system data to trace an observed oscillation back to its source. The difficulty of this problem is exacerbated by the fact that the parameters associated with a given power system model can range from slightly uncertain to entirely unknown. In this paper, a Bayesian framework, via a two-stage Maximum A Posteriori optimization routine, is employed in order to locate the most probable source of a forced oscillation given an uncertain prior model. The approach leverages an equivalent circuit representation of the system in the frequency domain and employs a numerical procedure which makes the problem suitable for real time application. The derived framework lends itself to successful performance in the presence of PMU measurement noise, high generator parameter uncertainty, and multiple forced oscillations occurring simultaneously. The approach is tested on a 4-bus system with a single forced oscillation source and on the WECC 179-bus system with multiple oscillation sources., Comment: 9 pages; submitted to IEEE Transactions on Power Systems

This paper investigates the water network's potential ability to provide demand response services to the power grid under the framework of a distribution-level water-energy nexus (micro-WEN). In particular, the hidden controllability of water loads, such as irrigation systems, was closely studied to improve the flexibility of electrical grids. A optimization model is developed for the demand-side management (DSM) of micro-WEN, and the simulation results assert that grid flexibility indeed benefits from controllable water loads. Although the proposed optimal DSM model is an intractable mixed-integer nonlinear programming (MINLP) problem, quasi-convex hull techniques were developed to relax the MINLP into a mixed-integer convex programming (MICP) problem. The numerical study shows that the quasi-convex hull relaxation is tight and that the resulting MICP problem is computationally efficient.

The solvability set of a power network - the set of all power injection vectors for which the corresponding Power Flow equations admit a solution - is central to power systems stability and security, as well as to the tightness of Optimal Power Flow relaxations. Whenever the solvability set is convex, this allows for substantial simplifications of various optimization and risk assessment algorithms. In this paper we focus on the solvability set of power distribution networks and prove convexity of the full solvability set (real and reactive powers) for tree homogeneous networks with the same r/x ratio for all elements. We also show this result can not be improved: once the network is not homogeneous, the convexity is immediately lost. It is nevertheless the case that if the network is almost homogeneous, a substantial practically-important part of the solvability set is still convex. Finally, we prove convexity of real solvability set (only real powers) for any tree network as well as for purely resistive networks with arbitrary topology., Comment: 6 pages

The Primal-Dual (PD) algorithm is widely used in convex optimization to determine saddle points. While the stability of the PD algorithm can be easily guaranteed, strict contraction is nontrivial to establish in most cases. This work focuses on continuous, possibly non-autonomous PD dynamics arising in a network context, in distributed optimization, or in systems with multiple time-scales. We show that the PD algorithm is indeed strictly contracting in specific metrics and analyze its robustness establishing stability and performance guarantees for different approximate PD systems. We derive estimates for the performance of multiple time-scale multi-layer optimization systems, and illustrate our results on a primal-dual representation of the Automatic Generation Control of power systems., Comment: 6 pages, 1 figures, published on LCSS and CDC 2018

The aggressive integration of distributed renewable sources is changing the dynamics of the electric power grid in an unexpected manner. As a result, maintaining conventional performance specifications, such as transient stability, may not be sufficient to ensure its reliable operation in stressed conditions. In this paper, we introduce a novel criteria in transient stability with consideration of operational constraints over frequency deviation and angular separation. In addition, we provide a robustness measure of the region of attraction, which can quantify the ability of the post-fault system to remain synchronized even under disturbances. To assess this new stability specification, we adopt the notion of Input-to-State Stability (ISS) to the context of power systems and introduce a new class of convex Lyapunov functions, which will result in tractable convex-optimization-based stability certificates. As a result, we are able to quantify the level of disturbance a power system can withstand while maintaining its safe operation. We illustrate the introduced stability specification and certificate on the IEEE 9 bus system., Comment: Accepted for publication in the IEEE Control Systems Letters (L-CSS)

DC grids are prone to small-signal instabilities due to the presence of tightly controlled loads trying to keep the power consumption constant over range of terminal voltage variations. Th, so-called, constant power load (CPL) represents a classical example of this destabilizing behavior acting as an incremental negative resistance. Real-life DC loads represented by controlled power converters exhibit the CPL behavior over a finite frequency range. There exist a number of methods for stability certification of DC grids which are primarily concerned with the source-load interaction and do not explicitly account for the influence of network. In the present manuscript, we develop a method for stability assessment of arbitrary DC grids by introducing the Augmented Power Dissipation and showing that it's positive definiteness is a sufficient condition for stability. We present an explicit expression for this quantity through load and network impedances and show how it could be directly used for stability certification of networks with arbitrary configuration., Comment: 10 pages

The wide deployment of renewable generation and the gradual decrease in the overall system inertia make modern power grids more vulnerable to transient instabilities and unacceptable frequency fluctuations. Time-domain simulation-based assessment of the system robustness against uncertain and stochastic disturbances is extremely time-consuming. In this paper, we develop an alternative approach, which has its roots in the input-output stability analysis for Lur'e systems. Our approach consists of a mathematically rigorous characterization of the external disturbances that the power system is transiently stable and the frequency constraints are not violated. The derived certificate is efficiently constructed via convex optimization and is shown to be non-conservative for different IEEE test cases.

The convex restriction of the power flow feasible sets identifies the convex subset of power injections where the solution for power flow is guaranteed to exist and satisfy the operational constraints. In contrast to convex relaxations, the convex restriction provides a sufficient condition for power flow feasibility and is particularly useful for problems involving uncertainty in the power generation and demand. In this paper, we present a general framework of constructing convex restriction of an algebraic set defined by equality and inequality constraints and apply the framework to power flow feasibility problem. The procedure results in convex quadratic constraints that provide a sufficiently large region for practical operation of the grid.

This paper introduces a micro-nexus of water and energy which can be considered as one of the physical infrastructures of the future building/city/village systems. For the electricity side, an alternating current (AC) power flow model integrated with battery energy storage and renewable generation is adopted. The nonlinear hydraulic characteristics in pipe networks is also considered in the proposed micro water-energy nexus (WEN) model. Integer variables are involved to represent the on/off state of pumps. Base on the proposed nexus model, a co-optimization framework of water and energy networks is developed. The overall co-optimization model is a mixed-integer nonlinear programming problem which is tested on a waterenergy nexus which consists of the IEEE 13-bus distribution system and a 8-node water distribution network. The simulation results demonstrate that the cost-efficiency of the co-optimization framework is higher than optimizing two systems independently., Comment: Submitted to IEEE Power & Energy Society 2018 General Meeting

For many nonlinear control systems, the chosen equilibrium determines both the steady-state efficiency and the dynamic performance. This paper addresses the issue of obtaining an optimal equilibrium in terms of some steady-state operation criteria for a Lur'e-type system and such an equilibrium can also guarantee a sufficiently large stability region in the dynamic domain such that the system can tolerate some given disturbance. For this purpose, a set of computationally tractable algebraic constraints, which can properly represent the stability certificate under the optimization framework, are proposed. The existing methods formulate the dynamic performance under the optimization framework by discretizing the differential-algebraic equations, which are computationally intractable for large-scale Lur'e systems like power grids. Dissimilarly, the introduced approach first constructs the stability region based on quadratic Lyapunov functions. Then, a novel method is proposed to project the stability region onto the feasible domain of the optimization problem such that the stability certificate can be incorporated into the optimization framework easily. In the transient stability-constrained optimal power flow (TSCOPF) problem of power systems, researchers look for a steady-state operating point with the minimum generation costs that can maintain system stability under some given transient disturbances. The proposed approach is applied to develop a scalable TSCOPF framework for power systems. The TSCOPF model is tested on the IEEE 118-Bus power system., Comment: Submitted to IEEE Transactions on Automatic Control

We propose a scalable optimization framework for estimating convex inner approximations of the steady-state security sets. The framework is based on Brouwer fixed point theorem applied to a fixed-point form of the power flow equations. It establishes a certificate for the self-mapping of a polytope region constructed around a given feasible operating point. This certificate is based on the explicit bounds on the nonlinear terms that hold within the self-mapped polytope. The shape of the polytope is adapted to find the largest approximation of the steady-state security region. While the corresponding optimization problem is nonlinear and non-convex, every feasible solution found by local search defines a valid inner approximation. The number of variables scales linearly with the system size, and the general framework can naturally be applied to other nonlinear equations with affine dependence on inputs. Test cases, with the system sizes up to $1354$ buses, are used to illustrate the scalability of the approach. The results show that the approximated regions are not unreasonably conservative and that they cover substantial fractions of the true steady-state security regions for most medium-sized test cases., Comment: 11 pages, 5 figures

Locating the sources of forced low-frequency oscillations in power systems is an important problem. A number of proposed methods demonstrate their practical usefulness, but many of them rely on strong modeling assumptions and provide poor performance in certain cases for reasons still not well understood. This paper proposes a systematic method for locating the source of a forced oscillation by considering a generator's response to fluctuations of its terminal voltages and currents. It is shown that a generator can be represented as an effective admittance matrix with respect to low-frequency oscillations, and an explicit form for this matrix, for various generator models, is derived. Furthermore, it is shown that a source generator, in addition to its effective admittance, is characterized by the presence of an effective current source thus giving a natural qualitative distinction between source and nonsource generators. Detailed descriptions are given of a source detection procedure based on this developed representation, and the method's effectiveness is confirmed by simulations on the recommended testbeds (eg. WECC 179-bus system). This method is free of strong modeling assumptions and is also shown to be robust in the presence of measurement noise and generator parameter uncertainty., Comment: 13 pages

Ad hoc electrical networks are formed by connecting power sources and loads without planning the interconnection structure (topology) in advance. They are designed to be installed and operated by individual communities---without central oversight---and as a result are well-suited to addressing the lack of electricity access in rural and developing areas. However, ad hoc networks are not widely used, and a major technical challenge impeding their development (and deployment) is the difficulty of certifying network stability without a priori knowledge of the topology. We develop conditions on individual power sources and loads such that a microgrid comprised of many units will be stable. We use Brayton-Moser potential theory to develop design constraints on individual microgrid components that certify transient stability---guaranteeing that the system will return to a suitable equilibrium after load switching events. Our central result is that stability can be ensured by installing a parallel capacitor at each constant power load, and we derive an expression for the required capacitance., Comment: 6 pages

In this paper, a hybrid measurement- and model-based method is proposed which can estimate the dynamic state Jacobian matrix and the dynamic system state matrix in near real-time utilizing statistical properties extracted from PMU measurements. The proposed method can be used to detect and identify network topology changes that have not been reflected in an assumed network model. Additionally, an application of the estimated system state matrix in online dynamic stability monitoring is presented., Comment: To appear in IEEE Transactions on Power Systems

Power injection uncertainties in distribution power grids, which are mostly induced by aggressive introduction of intermittent renewable sources, may drive the system away from normal operating regimes and potentially lead to the loss of long-term voltage stability (LTVS). Naturally, there is an ever increasing need for a tool for assessing the LTVS of a distribution system. This paper presents a fast and reliable tool for constructing \emph{inner approximations} of LTVS regions in multidimensional injection space such that every point in our constructed region is guaranteed to be solvable. Numerical simulations demonstrate that our approach outperforms all existing inner approximation methods in most cases. Furthermore, the constructed regions are shown to cover substantial fractions of the true voltage stability region. The paper will later discuss a number of important applications of the proposed technique, including fast screening for viable injection changes, constructing an effective solvability index and rigorously certified loadability limits., Comment: 11 pages, 9 figs, published on Transactions on Smart Grid

Transient stability assessment of power systems needs to account for increased risk from uncertainties due to the integration of renewables and distributed generators. The uncertain operating condition of the power grid hinders reliable assessment of transient stability. Conventional approaches such as time-domain simulations and direct energy methods are computationally expensive to take account of uncertainties. This paper proposes a reachability analysis approach that computes bounds of the possible trajectories from uncertain initial conditions. The eigenvalue decomposition is used to construct a polytopic template with a scalable number of hyperplanes that is guaranteed to converge near the equilibrium. The proposed algorithm bounds the possible states at a given time with a polytopic template and solves the evolution of the polytope over time. The problem is solved with linear programming relaxation based on outer-approximations of nonlinear functions, which is scalable for large scale systems. We demonstrate our method on IEEE test cases to certify the stability and bound the state trajectories.

Quadratic systems of equations appear in several applications. The results in this paper are motivated by quadratic systems of equations that describe equilibrium behavior of physical infrastructure networks like the power and gas grids. The quadratic systems in infrastructure networks are parameterized- the parameters can represent uncertainty (estimation error in resistance/inductance of a power transmission line, for example)or controllable decision variables (power outputs of generators,for example). It is then of interest to understand conditions on the parameters under which the quadratic system is guaranteed to have a solution within a specified set (for example, bounds on voltages and flows in a power grid). Given nominal values of the parameters at which the quadratic system has a solution and the Jacobian of the quadratic system at the solution i snon-singular, we develop a general framework to construct convex regions around the nominal value such that the system is guaranteed to have a solution within a given distance of the nominal solution. We show that several results from recen tliterature can be recovered as special cases of our framework,and demonstrate our approach on several benchmark power systems.

In modern power systems, the operating point, at which the demand and supply are balanced, may take different values due to changes in loads and renewable generation levels. Understanding the dynamics of stressed power systems with a range of operating points would be essential to assuring their reliable operation, and possibly allow higher integration of renewable resources. This letter introduces a non-traditional way to think about the stability assessment problem of power systems. Instead of estimating the set of initial states leading to a given operating condition, we characterize the set of operating conditions that a power grid converges to from a given initial state under changes in power injections and lines. We term this problem as "inverse stability", a problem which is rarely addressed in the control and systems literature, and hence, poorly understood. Exploiting quadratic approximations of the system's energy function, we introduce an estimate of the inverse stability region. Also, we briefly describe three important applications of the inverse stability notion: (i) robust stability assessment of power systems w.r.t. different renewable generation levels, (ii) stability-constrained optimal power flow (sOPF), and (iii) stability-guaranteed corrective action design.

In this paper, a hybrid measurement- and model-based method is proposed which can estimate the dynamic state Jacobian matrix in near real-time. The proposed method is computationally efficient and robust to the variation of network topology. A numerical example is given to show that the proposed method is able to provide good estimation for the dynamic state Jacobian matrix and is superior to the model-based method under undetectable network topology change. The proposed method may also help identify big discrepancy in the assumed network model., Comment: IEEE International Conference on Circuits and Systems (ISCAS) 2017

This paper studies the contraction properties of nonlinear differential-algebraic equation (DAE) systems. Specifically we develop scalable techniques for constructing the attraction regions associated with a particular stable equilibrium, by establishing the relation between the contraction rates of the original systems and the corresponding virtual extended systems. We show that for a contracting DAE system, the reduced system always contracts faster than the extended ones; furthermore, there always exists an extension with contraction rate arbitrarily close to that of the original system. The proposed construction technique is illustrated with a power system example in the context of transient stability assessment., Comment: 9 pages, 3 figures, submitted to TAC

Proper modeling of inverter-based microgrids is crucial for accurate assessment of stability boundaries. It has been recently realized that the stability conditions for such microgrids are significantly different from those known for large- scale power systems. While detailed models are available, they are both computationally expensive and can not provide the insight into the instability mechanisms and factors. In this paper, a computationally efficient and accurate reduced-order model is proposed for modeling the inverter-based microgrids. The main factors affecting microgrid stability are analyzed using the developed reduced-order model and are shown to be unique for the microgrid-based network, which has no direct analogy to large-scale power systems. Particularly, it has been discovered that the stability limits for the conventional droop-based system (omega - P/V - Q) are determined by the ratio of inverter rating to network capacity, leading to a smaller stability region for microgrids with shorter lines. The theoretical derivation has been provided to verify the above investigation based on both the simplified and generalized network configurations. More impor- tantly, the proposed reduced-order model not only maintains the modeling accuracy but also enhances the computation efficiency. Finally, the results are verified with the detailed model via both frequency and time domain analyses.

Power grids normally operate at some stable operating condition where power supply and demand are balanced. In response to emergency situations, load shedding is a prevailing approach where local protective devices are activated to cut a suitable amount of load to quickly rebalance the supply demand and hopefully stabilize the system. This traditional emergency control results in interrupted service with severe economic damage to customers. Also, such control is usually less effective due to the lack of coordination among protective devices. In this paper, we propose a novel structural emergency control to render post-fault dynamics from the critical/emergency fault-cleared state to the stable equilibrium point. This is a new control paradigm that does not rely on any continuous measurement or load shedding, as in the classical setup. Instead, the grid is made stable by discretely relocating the equilibrium point and its stability region such that the system is consecutively attracted from the fault-cleared state back to the original equilibrium point. The proposed control is designed by solving linear and convex optimization problems, making it possibly scalable to large-scale power grids. Finally, this emergency control scheme can be implemented by exploiting transmission facilities available on the existing grids.

In this paper, we analyze the impact of communication failures on the performance of optimal distributed frequency control. We consider a consensus-based control scheme, and show that it does not converge to the optimal solution when the communication network is disconnected. We propose a new control scheme that uses the dynamics of power grid to replicate the information not received from the communication network, and prove that it achieves the optimal solution under any single communication link failure. In addition, we show that this control improves cost under multiple communication link failures. Next, we analyze the impact of discrete-time communication on the performance of distributed frequency control. In particular, we will show that the convergence time increases as the time interval between two messages increases. We propose a new algorithm that uses the dynamics of the power grid, and show through simulation that it improves the convergence time of the control scheme significantly., Comment: 8 pages, 7 figures

Interconnected networks describe the dynamics of important systems in a wide range such as biological systems and electrical power grids. Some important features of these systems were successfully studied and understood through simplified model of linear interconnection of linear subsystems, where provably global properties, e.g. global convergence to a specific state, usually hold true. However, in severely disturbed conditions many of those systems exhibit strongly nonlinear behaviour. Particularly, multiple equilibrium points may coexist and make the dynamical behavior of the system difficult to predict. Aiming at understanding the fragility of interconnected systems, we will provide a hierarchical framework to assess the metastability and resilience of such systems. This framework is based on independently characterizing stability of individual subsystems when they are uncoupled from the network, and then enforcing the diagonal dominance property on a structure matrix capturing the subsystems stability and the input-to-output gains of interconnection network. Since the subsystems are usually of low order and the structure matrix has size equal to the number of subsystems, this framework is easy to implement and thus scalable to large scale interconnected systems. Possible application of this framework in assessing stability of microgrids will be discussed at the end of this paper.

Ad hoc electrical networks are formed by connecting power sources and loads without pre-determining the network topology. These systems are well-suited to addressing the lack of electricity in rural areas because they can be assembled and modified by non-expert users without central oversight. There are two core aspects to ad hoc system design: 1) designing source and load units such that the microgrid formed from the arbitrary interconnection of many units is always stable and 2) developing control strategies to autonomously manage the microgrid (i.e., perform power dispatch and voltage regulation) in a decentralized manner and under large uncertainty. To address these challenges we apply a number of nonlinear control techniques---including Brayton-Moser potential theory and primal-dual dynamics---to obtain conditions under which an ad hoc dc microgrid will have a suitable and asymptotically stable equilibrium point. Further, we propose a new decentralized control scheme that coordinates many sources to achieve a specified power dispatch from each. A simulated comparison to previous research is included., Comment: To be presented at the 2016 Conference on Decision and Control. Changes from v1: refined proposed control strategy, updated Lyapunov function to generalize stability analysis to broader class of networks, fixed errors and typos

Power systems normally operate at their stable operating conditions where the power supply and demand are balanced. In emergency situations, the operators proceed to cut a suitable amount of loads to rebalance the supply-demand and hopefully stabilize the system. This traditional emergency control scheme results in interrupted service with severely economic damages to customers. In order to provide seamless electricity service to customers, this paper proposes a viable alternative for traditional remedial controls of power grids by exploiting the plentiful transmission facilities. In particular, we consider two emergency control schemes involving adjustment of the susceptance of a number of selected transmission lines to drive either fault-on dynamics or post-fault dynamics, and thereby stabilize the system under emergency situations. The corresponding emergency control problems will be formulated and partly solved in some specific cases. Simple numerical simulation will be used to illustrate the concept of this paper., Comment: arXiv admin note: substantial text overlap with arXiv:1510.07325

In this paper, a hybrid measurement and model-based method is proposed which can estimate the dynamic state Jacobian matrix in near real-time. The proposed method is computationally efficient and robust to the variation of network topology. Since the estimated Jacobian matrix carries significant information on system dynamics and states, it can be utilized in various applications. In particular, two application of the estimated Jacobian matrix in online oscillation analysis, stability monitoring and control are illustrated with numerical examples. In addition, a side-product of the proposed method can facilitate model validation by approximating the damping of generators., Comment: submitted to IEEE Transactions on Power Systems

Traditional emergency control schemes in power systems usually accompany with power interruption yielding severely economic damages to customers. This paper sketches the ideas of a viable alternative for traditional remedial controls for power grids with high penetration of renewables, in which the renewables are integrated with synchronverters to mimic the dynamics of conventional generators. In this novel emergency control scheme, the power electronics resources are exploited to control the inertia and damping of the imitated generators in order to quickly compensate for the deviations caused by fault and thereby bound the fault-on dynamics and stabilize the power system under emergency situations. This emergency control not only saves investments and operating costs for modern and future power systems, but also helps to offer seamless electricity service to customers. Simple numerical simulation will be used to illustrate the concept of this paper., Comment: arXiv admin note: text overlap with arXiv:1504.04684

The power flow equations are at the core of most of the computations for designing and operating electric power systems. The power flow equations are a system of multivariate nonlinear equations which relate the power injections and voltages in a power system. A plethora of methods have been devised to solve these equations, starting from Newton-based methods to homotopy continuation and other optimization-based methods. While many of these methods often efficiently find a high-voltage, stable solution due to its large basin of attraction, most of the methods struggle to find low-voltage solutions which play significant role in certain stability-related computations. While we do not claim to have exhausted the existing literature on all related methods, this tutorial paper introduces some of the recent advances in methods for solving power flow equations to the wider power systems community as well as bringing attention from the computational mathematics and optimization communities to the power systems problems. After briefly reviewing some of the traditional computational methods used to solve the power flow equations, we focus on three emerging methods: the numerical polynomial homotopy continuation method, Groebner basis techniques, and moment/sum-of-squares relaxations using semidefinite programming. In passing, we also emphasize the importance of an upper bound on the number of solutions of the power flow equations and review the current status of research in this direction., Comment: 13 pages, 2 figures. Submitted to the Tutorial Session at IEEE 2016 American Control Conference

Ease of miniaturization and minimal maintenance are among the advantages for replacing conventional batteries with vibratory energy harvesters in a wide of range of disciplines and applications, from wireless communication sensors to medical implants. However, the current harvesters do not extract energy from the ambient vibrations in a very efficient and robust fashion, and hence, there need to be more optimal harvesting approaches. In this paper, we introduce a generic architecture for vibration energy harvesting and delineate the key challenges in the field. Then, we formulate an optimal control problem to maximize the harvested energy. Though possessing similar structure to that of the standard LQG problem, this optimal control problem is inherently different from the LQG problem and poses theoretical challenges to control community. As the first step, we simplify it to a tractable problem of optimizing control gains for a linear system subjected to Gaussian white noise excitation, and show that this optimal problem has non-trivial optimal solutions in both time and frequency domains.

This paper presents a tensor-recovery method to solve probabilistic power flow problems. Our approach generates a high-dimensional and sparse generalized polynomial-chaos expansion that provides useful statistical information. The result can also speed up other essential routines in power systems (e.g., stochastic planning, operations and controls). Instead of simulating a power flow equation at all quadrature points, our approach only simulates an extremely small subset of samples. We suggest a model to exploit the underlying low-rank and sparse structure of high-dimensional simulation data arrays, making our technique applicable to power systems with many random parameters. We also present a numerical method to solve the resulting nonlinear optimization problem. Our algorithm is implemented in MATLAB and is verified by several benchmarks in MATPOWER $5.1$. Accurate results are obtained for power systems with up to $50$ independent random parameters, with a speedup factor up to $9\times 10^{20}$., Comment: 8 pages, 10 figures, submitted to IEEE Trans. Power Systems

Narrow bandwidth and easy detuning, inefficiency in broadband and non-stationary excitations, and difficulties in matching linear harvester's resonance frequency to low-frequency excitations at small scales, have convinced the researchers to investigate the nonlinear, in particular the bistable energy harvesters in the recent years. However, the bistable harvesters suffer from co-existing low and high energy orbits, and sensitivity to initial conditions, and have been recently proven inefficient when subjected to many real-world random and non-stationary excitations. Here, we propose a novel buy-low-sell-high strategy that can significantly improve the harvester's efficiency at low-frequencies in a much more robust fashion. This strategy could be realized by a passive adaptive bistable system. Simulation results confirm high efficiency of the adaptive bistable system following a buy-low-sell-high logic when subjected to harmonic and random non-stationary walking excitations compared to its conventional bistable and linear counterparts.

Contingency screening for transient stability of large-scale, strongly nonlinear, interconnected power systems is one of the most computationally challenging parts of Dynamic Security Assessment and requires huge resources to perform time-domain simulations-based assessment. To reduce computational cost of time-domain simulations, direct energy methods have been extensively developed. However, these methods, as well as other existing methods, still rely on time-consuming numerical integration of the fault-on dynamics. This task is computationally hard, since possibly thousands of contingencies need to be scanned and thousands of accompanied fault-on dynamics simulations need to be performed and stored on a regular basis. In this paper, we introduce a novel framework to eliminate the need for fault-on dynamics simulations in contingency screening. This simulation-free framework is based on bounding the fault-on dynamics and extending the recently introduced Lyapunov Function Family approach for transient stability analysis of structure-preserving model. In turn, a lower bound of the critical clearing time (CCT) is obtained by solving convex optimization problems without relying on any time-domain simulations. A comprehensive analysis is carried out to validate this novel technique on a number of IEEE test cases.

We develop a new approach for construction of convex analytically simple regions where the AC power flow equations are guaranteed to have a feasible solutions. Construction of these regions is based on efficient semidefinite programming techniques accelerated via sparsity exploiting algorithms. Resulting regions have a simple geometric shape in the space of power injections (polytope or ellipsoid) and can be efficiently used for assessment of system security in the presence of uncertainty. Efficiency and tightness of the approach is validated on a number of test networks.

Sustained oscillations observed in power systems can damage equipment, degrade the power quality and increase the risks of cascading blackouts. There are several mechanisms that can give rise to oscillations, each requiring different countermeasure to suppress or eliminate the oscillation. This work develops mathematical framework for analysis of sustained oscillations and identifies statistical signatures of each mechanism, based on which a novel oscillation diagnosis method is developed via real-time processing of phasor measurement units (PMUs) data. Case studies show that the proposed method can accurately identify the exact mechanism for sustained oscillation, and meanwhile provide insightful information to locate the oscillation sources., Comment: The paper has been accepted by IEEE Transactions on Power Systems

We study the interaction between the power grid and the communication network used for its control. We design a centralized emergency control scheme under both full and partial communication support, to improve the performance of the power grid. We use our emergency control scheme to model the impact of communication loss on the grid. We show that unlike previous models used in the literature, the loss of communication does not necessarily lead to the failure of the correspondent power nodes; i.e. the "point-wise" failure model is not appropriate. In addition, we show that the impact of communication loss is a function of several parameters such as the size and structure of the power and communication failure, as well as the operating mode of power nodes disconnected from the communication network. Our model can be used to design the dependency between the power grid and the communication network used for its control, so as to maximize the benefit in terms of intelligent control, while minimizing the risks due to loss of communication., Comment: 6 pages, 8 figures, Submitted to SmartGridComm 2015

Security assessment of large-scale, strongly nonlinear power grids containing thousands to millions of interacting components is a computationally expensive task. Targeting at reducing the computational cost, this paper introduces a framework for constructing a robust assessment toolbox that can provide mathematically rigorous certificates for the grids' stability in the presence of variations in power injections, and for the grids' ability to withstand a bunch sources of faults. By this toolbox we can "off-line" screen a wide range of contingencies or power injection profiles, without reassessing the system stability on a regular basis. In particular, we formulate and solve two novel robust stability and resiliency assessment problems of power grids subject to the uncertainty in equilibrium points and uncertainty in fault-on dynamics. Furthermore, we bring in the quadratic Lyapunov functions approach to transient stability assessment, offering real-time construction of stability/resiliency certificates and real-time stability assessment. The effectiveness of the proposed techniques is numerically illustrated on a number of IEEE test cases.

Power flow solvable boundary plays an important role in contingency analysis, security assessment, and planning processes. However, to construct the real solvable boundary in multidimensional parameter space is burdensome and time consuming. In this paper, we develop a new technique to approximate the solvable boundary of distribution systems based on Banach fixed point theorem. Not only the new technique is fast and non-iterative, but also the approximated boundary is more valuable to system operators in the sense that it is closer to the feasible region. Moreover, a simple solvable criterion is also introduced that can serve as a security constraint in various planning and operational problems., Comment: 5 pages, 6 figures, accepted to PESGM 2015

The increasing development of the electric power grid, the largest engineered system ever, to an even more complicated and larger system requires a new generation of stability assessment methods that are computationally tractable and feasible in real-time. In this paper we first extend the recently introduced Lyapunov Functions Family (LFF) transient stability assessment approach, that has potential to reduce the computational cost on large scale power grids, to structure-preserving power grids. Then, we introduce a new geometry-based method to construct the stability region estimate of power systems. Our conceptual demonstration shows that this new method can certify stability of a broader set of initial conditions compared to the minimization-based LFF method and the energy methods (closest UEP and controlling UEP methods).