Konstantin Turitsyn

Thanh Long Vu, Konstantin Turitsyn

Analysis of transient stability of strongly nonlinear post-fault dynamics is one of the most computationally challenging parts of Dynamic Security Assessment. This paper proposes a novel approach for assessment of transient stability of the system. The approach generalizes the idea of energy methods, and extends the concept of energy function to a more general Lyapunov Functions Family (LFF) constructed via Semi-Definite-Programming techniques. Unlike the traditional energy function and its variations, the constructed Lyapunov functions are proven to be decreasing only in a finite neighborhood of the equilibrium point. However, we show that they can still certify stability of a broader set of initial conditions in comparison to the traditional energy function in the closest-UEP method. Moreover, the certificates of stability can be constructed via a sequence of convex optimization problems that are tractable even for large scale systems. We also propose specific algorithms for adaptation of the Lyapunov functions to specific initial conditions and demonstrate the effectiveness of the approach on a number of IEEE test cases.

Petr Vorobev, Po-Hsu Huang, Mohamed Al Hosani et al.

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.

Konstantin Turitsyn, Scott Backhaus, Maxim Ananyev et al.

We introduce and analyze Markov Decision Process (MDP) machines to model individual devices which are expected to participate in future demand-response markets on distribution grids. We differentiate devices into the following four types: (a) optional loads that can be shed, e.g. light dimming; (b) deferrable loads that can be delayed, e.g. dishwashers; (c) controllable loads with inertia, e.g. thermostatically-controlled loads, whose task is to maintain an auxiliary characteristic (temperature) within pre-defined margins; and (d) storage devices that can alternate between charging and generating. Our analysis of the devices seeks to find their optimal price-taking control strategy under a given stochastic model of the distribution market.

Thanh Long Vu, Hung Dinh Nguyen, Alexandre Megretski et al.

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.

René Pfitzner, Konstantin Turitsyn, Michael Chertkov

We study the evolution of fast blackout cascades in the model of the Polish (transmission) power grid (2700 nodes and 3504 transmission lines). The cascade is initiated by a sufficiently severe initial contingency tripping. It propagates via sequential trippings of many more overheated lines, islanding loads and generators and eventually arriving at a fixed point with the surviving part of the system being power-flow-balanced and the rest of the system being outaged. Utilizing an improved form of the quasi-static model for cascade propagation introduced in our earlier study (Statistical Classification of Cascading Failures in Power Grids, IEEE PES GM 2011), we analyze how the severity of the cascade depends on the order of tripping overheated lines. Our main observation is that the order of tripping has a tremendous effect on the size of the resulting outage. Finding the "best" tripping, defined as causing the least damage, constitutes a difficult dynamical optimization problem, whose solution is most likely computationally infeasible. Instead, here we study performance of a number of natural heuristics, resolving the next switching decision based on the current state of the grid. Overall, we conclude that controlled intentional tripping is advantageous in the situation of a fast developing extreme emergency, as it provides significant mitigation of the resulting damage.

Samuel Chevalier, Petr Vorobev, Konstantin Turitsyn

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.

Thanh Long Vu, Surour Al Araifi, Mohamed Elmoursi et al.

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.

Samuel Chevalier, Petr Vorobev, Konstantin Turitsyn

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.

Krishnamurthy Dvijotham, Hung Nguyen, Konstantin Turitsyn

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.

Thanh Long Vu, Spyros Chatzivasileiadis, Hsiao-Dong Chiang et al.

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.

Marzieh Parandehgheibi, Konstantin Turitsyn, Eytan Modiano

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.

Xiaozhe Wang, Konstantin Turitsyn

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.

Kathleen Cavanagh, Julia A. Belk, Konstantin Turitsyn

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.

Dmitry Podolsky, Konstantin Turitsyn

For a power system operating in the vicinity of the power transfer limit of its transmission system, effect of stochastic fluctuations of power loads can become critical as a sufficiently strong such fluctuation may activate voltage instability and lead to a large scale collapse of the system. Considering the effect of these stochastic fluctuations near a codimension 1 saddle-node bifurcation, we explicitly calculate the autocorrelation function of the state vector and show how its behavior explains the phenomenon of critical slowing-down often observed for power systems on the threshold of blackout. We also estimate the collapse probability/mean clearing time for the power system and construct a new indicator function signaling the proximity to a large scale collapse. The new indicator function is easy to estimate in real time using PMU data feeds as well as SCADA information about fluctuations of power load on the nodes of the power grid. We discuss control strategies leading to the minimization of the collapse probability.

Samuel Chevalier, Petr Vorobev, Konstantin Turitsyn et al.

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.

Thanh Long Vu, Konstantin Turitsyn

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.

Kathleen Cavanagh, Petr Vorobev, Konstantin Turitsyn

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.

Thanh Long Vu, Hsiao-Dong Chiang, Konstantin Turitsyn

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.

Thanh Long Vu, Konstantin Turitsyn

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).

Ashkan Haji Hosseinloo, Alexander Ryzhov, Aldo Bischi et al.

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.

Xiaozhe Wang, Janusz Bialek, Konstantin Turitsyn

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.

Julia A. Belk, Wardah Inam, David J. Perreault et al.

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.

Thanh Long Vu, Konstantin Turitsyn

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.

Xiaozhe Wang, Konstantin Turitsyn

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.

Thanh Long Vu, Spyros Chatzivasileiadis, Konstantin Turitsyn

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.

Xiaozhe Wang, Konstantin Turitsyn

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.

Dhagash Mehta, Daniel K Molzahn, Konstantin Turitsyn

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.

Ashkan Haji Hosseinloo, Thanh Long Vu, Konstantin Turitsyn

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.

Krishnamurthy Dvijotham, Konstantin Turitsyn

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.

Marzieh Parandehgheibi, Konstantin Turitsyn, Eytan Modiano

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.

Hung D. Nguyen, Konstantin Turitsyn

High levels of penetration of distributed generation and aggressive reactive power compensation may result in the reversal of power flows in future distribution grids. The voltage stability of these operating conditions may be very different from the more traditional power consumption regime. This paper focused on demonstration of multistability phenomenon in radial distribution systems with reversed power flow, where multiple stable equilibria co-exist at the given set of parameters. The system may experience transitions between different equilibria after being subjected to disturbances such as short-term losses of distributed generation or transient faults. Convergence to an undesirable equilibrium places the system in an emergency or \textit{in extremis} state. Traditional emergency control schemes are not capable of restoring the system if it gets entrapped in one of the low voltage equilibria. Moreover, undervoltage load shedding may have a reverse action on the system and can induce voltage collapse. We propose a novel pulse emergency control strategy that restores the system to the normal state without any interruption of power delivery. The results are validated with dynamic simulations of IEEE $13$-bus feeder performed with SystemModeler software. The dynamic models can be also used for characterization of the solution branches via a novel approach so-called the admittance homotopy power flow method.

Hung D. Nguyen, Konstantin Turitsyn

Dynamic response of loads has a significant effect on system stability and directly determines the stability margin of the operating point. Inherent uncertainty and natural variability of load models make the stability assessment especially difficult and may compromise the security of the system. We propose a novel mathematical "robust stability" criterion for the assessment of small-signal stability of operating points. Whenever the criterion is satisfied for a given operating point, it provides mathematical guarantees that the operating point will be stable with respect to small disturbances for any dynamic response of the loads. The criterion can be naturally used for identification of operating regions secure from the occurrence of Hopf bifurcation. Several possible applications of the criterion are discussed, most importantly the concept of Robust Stability Assessment (RSA) that could be integrated in dynamic security assessment packages and used in contingency screening and other planning and operational studies.