Compositional and Equilibrium-Free Stability Certification for Power Systems--Part II: Algorithms and Applications

arXiv:2506.11411 · eess.SY, cs.SY · Submitted 2025-06-13 · Read on arXiv

Listen

Transcript

Introduction to the show: ident: Robotics Radio. Generated commentary on the latest robotics and control papers.

Rosa: Today's paper: "Compositional and Equilibrium-Free Stability Certification for Power Systems--Part II".

Dev: The gist This two-part paper proposes a compositional and equilibrium-free approach to analyzing power system stability.

Rosa: First, who's behind it and why it matters.

Paper summary: Rosa: So we're looking at this paper today, "Compositional and Equilibrium-Free Stability Certification for Power Systems--Part II: Algorithms and Applications." Essentially, the authors are proposing a way to analyze power system stability that doesn't rely on finding a specific equilibrium first.

Dev: Right. They build on something they did in Part I, which established these stability conditions based on what they call delta dissipativity. The main claim here is that this approach helps us overcome some big limitations of the older, traditional methods.

Taro: Those limitations include scalability issues and privacy concerns when you're looking at huge, complex grids. It sounds like they're trying to create a framework that can handle those things better without getting bogged down in finding a single steady state.

Rosa: Exactly. In Part II, they focus on how to actually use this theory for real, complex power grids by proposing two main methods: one for checking the local condition of delta dissipativity and another for verifying the coupling condition using something called Alternating Direction Method of Multipliers, or ADMM.

Dev: That means they are moving beyond just the theory and giving us a concrete way to apply it to heterogeneous devices, which is when you have different types of equipment all interacting in the system. They also propose a distributed computational framework for checking that coupling condition.

Taro: So, what matters here for me is how this handles misbehavior. If the world misbehaves and the system shifts equilibria quickly, this method allows us to evaluate stability under those shifting conditions because it's equilibrium-free.

Rosa: That’s right. And they show off three key applications using modified IEEE benchmark systems—specifically the nine-bus, thirty-nine-bus, and one hundred eighteen-bus grids <ref:2506.11411#pg1,9-bus, 39-bus, and 118-bus>. These case studies really validate their theory and methods across different system sizes.

Dev: So, what we're seeing is a systematic process for verifying local delta dissipativity by first transforming device models into a standardized input-output form, then using a Krasovskii-type storage function to check the inequality.

Conclusion: Rosa: Looking at the whole paper, "Compositional and Equilibrium-Free Stability Certification for Power Systems--Part II: Algorithms and Applications," it really shows how you can build a stability analysis tool that is modular. The authors are using this compositional approach to tackle stability in massive, diverse power systems.

Dev: I agree. The implication is that we might be able to certify the stability of huge grids without having to solve for every possible equilibrium point beforehand, which saves a ton of computational effort and gives us more flexibility in testing different operating conditions.

Taro: For someone just listening, it means there's a way to check if a system is stable across its entire range of behavior dynamically, not just at one fixed point. That's what shifts the focus from finding static solutions to understanding the system's overall dynamic behavior.

Rosa: Right. And they show this works with multiple equilibria, meaning you can check stability for different possible steady states simultaneously using a theorem in Part I which is linked here in Part II.

Dev: The coupling condition verification using ADMM is pretty neat because it allows for a distributed computing framework. This means we can have subsystems check their local conditions independently without needing one big central computer to handle everything, which addresses those privacy and scalability concerns they mentioned upfront.

Taro: That's the practical part I care about. If you have thousands of devices, you don't want one bottleneck controlling the entire verification process, especially when trying to keep sensitive operational data private between subsystems.

Rosa: So, the overall message is that this framework provides a scalable and modular way to assess stability in modern power systems by separating the local device checking from the global coupling condition check. That's what they achieved with these case studies on those IEEE benchmarks.

eess.SY, cs.SY

Submitted: 2025-06-13

Updated: 2026-10-08

License: http://arxiv.org/licenses/nonexclusive-distrib/1.0/

Importance score: 90/100

The gist: The gist This two-part paper proposes a compositional and equilibrium-free approach to analyzing power system stability.

Key concepts

Delta Dissipativity: A mathematical property used to prove stability in power systems. It ensures that energy flows out of a system, which is crucial for proving asymptotic stability without needing to find a specific stable operating point.
Local Condition Verification: The process of checking if individual components (devices) satisfy the delta dissipativity requirement. This involves transforming device models into a standardized form and using a Krasovskii-type storage function to mathematically confirm the required dissipative inequality holds.
Alternating Direction Method of Multipliers (ADMM): A distributed optimization algorithm used to solve large-scale problems, like verifying the coupling condition for stability. It breaks down the complex global problem into smaller, manageable local problems solved by different subsystems, improving computational efficiency and scalability.
Coupling Condition: A necessary constraint in power system stability that ensures all interconnected parts behave correctly together. The paper uses ADMM to verify this condition efficiently in a distributed manner, allowing large systems to be analyzed without needing a single central coordinator.

Terminology

Verification — 13× in this paper · explained in 2 episodes
Benchmark — 12× in this paper · explained in 2 episodes
Alternating Direction Method of Multipliers — 4× in this paper
Core — 3× in this paper · explained in 2 episodes
Alternating Direction Method of Multipliers (ADMM) — 2× in this paper · explained in 2 episodes
Local Optimization — 1× in this paper · explained in 2 episodes
Steering — 1× in this paper · explained in 2 episodes
Trace — 1× in this paper · explained in 2 episodes

Summary

The gist This two-part paper proposes a compositional and equilibrium-free approach to analyzing power system stability.

How it works

The authors developed a compositional and equilibrium-free stability theory based on delta dissipativity in Part I, which provides a theoretical foundation for overcoming the limitations of traditional centralized and equilibrium-oriented methods such as poor scalability, inadequate privacy protection, and high computational demands. In Part II, they focus on methods for applying this theory to complex power grids by proposing a method to verify the local condition of delta dissipativity for heterogeneous devices and a method to verify the coupling condition based on Alternating Direction Method of Multipliers (ADMM).

Verification of Local Conditions

The process for verifying local delta dissipativity involves two primary steps: (i) transforming device models into a standardized input-output form, and (ii) validating dissipativity using a Krasovskii-type storage function. The authors introduce a model transformation step to reconcile existing device representations with the proposed framework, demonstrating this process using widely adopted dynamic models for synchronous generators and inverter-based resources.

Verification of dissipativity is achieved by finding a proper storage function that satisfies the dissipative inequality, specifically utilizing a Krasovskii’s type storage function, i.e., S(x, u) = f(x, u)TPf(x, u) for some positive definite matrix P. Proposition 1 provides a sufficient condition to verify the delta dissipativity with such a storage function. For linear input-state-output systems, this condition degenerates into a linear matrix equality independent of x and u.

Method for Verifying the Coupling Condition

The coupling condition can be formulated as a feasibility problem in optimization. To address scalability issues in large-scale systems, the authors propose a distributed computational framework based on the ADMM algorithm [12] to verify this condition.

The distributed algorithm involves several steps:

X-update:

Each dynamic bus i ∈ V1 solves a problem involving l d i(Xi,Pi) ⪯ -tiI, Pi ≻ 0, ti ≤ t¯.

p-check:

If for some error tolerance ε > 0, we have ti > ε for all i ∈ V, then the p-check step is performed by solving a problem involving l c(p1Xk+1,..., pN Xk+1) ⪯ tzI p1 + · · · + pN = N. If the optimal tz < −ϵ, the p-check passes and justifies the coupling condition and the iteration terminates.

Z-update:

If p-check fails or is not performed, the system coordinator solves a problem involving l c(Z1,..., ZN) ⪯ 0.

Applications

The paper demonstrates three key applications of the proposed framework using IEEE benchmark systems.

Application 1: Stability Assessment for Multiple Equilibria:

The method enables simultaneous stability evaluation across all equilibria within a specified region via Theorem 2 (Part I). For any equilibrium (x 0, u 0) of the system, each subsystem can independently verify whether its local state (x 0i, u 0i) lies within Di and hence certificate the system-wide asymptotic stability.

Application 2: Stability Assessment under Varying Operating Conditions:

Thanks to the equilibrium-free nature of delta dissipativity, our method enables rapid stability evaluation under shifting equilibria. The framework pre-computes the dissipativity matrix Xi and verifies the equilibrium-independent coupling condition.

Application 3: Distributed Stability Assessment:

Algorithm 1 enables a distributed cloud-edge-like framework to certificate stability of large-scale systems. This approach reduces the computational load on the central coordinator while preserving the privacy of individual subsystems, as each subsystem only needs to share its Xi matrix rather than its detailed model. The distributed nature of our algorithm also reduces the computational burden on central coordinators and enhances privacy.

The case studies on modified IEEE 9-bus, 39-bus, and 118-bus benchmarks validated the methods and applications. The two-part study offers a comprehensive framework for scalable and modular stability analysis in modern power systems.

--- Page 9 ---

The flow chart of Algorithm 1 shows the iterative process involving X-update, p-check, Z-update, Y-update, Residuals update, and ρ-update. The p-check process is activated when min ti > epsilon. The iteration successfully terminates at k = 259 as the p-check process is passed.

--- Page 8 ---

The system equilibrium stably shifts from the nominal equilibrium to another after a disturbance occurred at t = 10s. This shows the asymptotic stability of these two equilibria.

--- Page 7 ---

The red dots represent the nominal inputs u 0i of each dynamic subsystem in Figure 3. The cross sections of Di, i = 33, 16 on the 2-dimensional input plane with xi = x 0i are shown in Figure 6.

--- Page 5 ---

The constant impedance load is delta-D(X, R 2) for any X with Q >= -Zlp00-Zlp. The proof for the constant impedance load involves showing that I∂h(u)T Q S S T R I∂h(u) >= Q + Zlp00 Zlp ⪯ 0.

--- Page 4 ---

The intermediate node in Example 6 and the constant voltage source in Example 7 are delta-D(X, R 2) for any X with Q >= 0. The proof for these static subsystems involves invoking Proposition 2.

--- Page 3 ---

The system (28) is transformed into a matrix inequality involving local constraints l d i(Xi,Pi) and l s i(Xi). The augmented Lagrangian of (28) is Lrho(X, Z, Y) = X i∈V L iρ(Xi, Zi, Yi), where with a little abuse of notation we let X, Z, and Y denote the collection of Xi, Zi, and Yi.

--- Page 6 ---

The optimal solution yields Xk+1. The augmented Lagrangian of (28) is Lrho(X, Z, Y) = X i∈V L iρ(Xi, Zi, Yi), where with a little abuse of notation we let X, Z, and Y denote the collection of Xi, Zi, and Yi.