pith. sign in

arxiv: 2511.04777 · v2 · submitted 2025-11-06 · 📡 eess.SP

OPF-Based Optimal Power System Network Restoration Considering Frequency Dynamics

Dawn Virginillo , Asja Dervi\v{s}kadi\'c , Mario Paolone This is my paper

Pith reviewed 2026-05-18 00:24 UTC · model grok-4.3

classification 📡 eess.SP
keywords power system restorationfrequency dynamicsoptimal power flowIEEE 9-bus modeldynamic stabilityrestoration planningswitching sequencelow inertia
0
0 comments X

The pith

Static optimal restoration sequences for power grids can violate frequency stability constraints during switching.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper formulates the optimal power system restoration problem using optimal power flow while adding frequency dynamics constraints to ensure stability during switching. It first validates the static version of the problem for global optimality via brute-force tree search on switching sequences. When the dynamic formulation is applied to the IEEE 9-bus model, the sequence that is optimal under static constraints produces frequency violations, demonstrating that low-inertia restoration scenarios require dynamic modeling. The work addresses growing interest in re-evaluating restoration plans amid blackouts and the energy transition toward lower system inertia.

Core claim

The authors pose and solve an OPF-based formulation of the optimal power system restoration problem that includes frequency dynamics. They validate the switching constraints for global optimality in a static version using brute-force tree search. Application of the dynamic formulation to the IEEE 9-bus model shows that the switching sequence optimal under the static formulation violates the dynamic frequency constraints, establishing the importance of dynamic considerations in PSR planning.

What carries the argument

OPF-based formulation of optimal power system restoration that incorporates frequency dynamics constraints on switching sequences.

If this is right

  • Restoration planning must integrate dynamic frequency constraints to avoid unstable sequences in low-inertia islands.
  • Optimal switching orders can change once frequency dynamics are enforced rather than checked after the fact.
  • Static-only optimization is no longer sufficient for reliable PSR in systems with reduced inertia.
  • Brute-force validation methods can confirm global optimality for the static subproblem before adding dynamics.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The approach could be scaled to larger test systems or real utility networks to quantify how often static plans fail dynamic checks.
  • Similar dynamic extensions might apply to other restoration objectives such as voltage stability or protection coordination.
  • Real-time implementation would require fast enough solvers and accurate inertia estimates from online measurements.

Load-bearing premise

The frequency dynamics model accurately represents real system behavior during the restoration switching sequence.

What would settle it

A high-fidelity dynamic simulation of the static optimal switching sequence on the IEEE 9-bus model that shows no frequency limit violations would falsify the claim that static plans are insufficient.

Figures

Figures reproduced from arXiv: 2511.04777 by Asja Dervi\v{s}kadi\'c, Dawn Virginillo, Mario Paolone.

Figure 1
Figure 1. Figure 1: IEEE 9-bus model, adapted from [27]. • θ is the matrix of bus voltage angles. • Sb the energization of buses (integer variable), defined by the number of live connections to each bus. Of these, it is noted that branch flows Pl , Pw and bus energization Sb are dependent variables and deterministically constrained. The constraints can be divided into four parts: power balance (OPF) constraints, initializatio… view at source ↗
Figure 2
Figure 2. Figure 2: Block diagram of the implemented governor, turbine, and generator dynamics. [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Block diagram of implemented governor/turbine model IEEEG3 [35]. [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Results of the inference experiments: (a) Experiment 1 (b) Experiment 2 (c) Experiment 3. [PITH_FULL_IMAGE:figures/full_fig_p012_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Results of DynOPF-R. the optimized switching sequence can be practically infeasible. The DVLP can be further used to validate the optimality of the dynamic setpoints, given the opti￾mized energization order as parameters. Indeed, the optimal frequency setpoints computed using DVLP with a frequency tolerance ∆fmax = 1.5Hz are ∆ωref = [5.330, 1.565, 0] rad/s. Finally, it is noted that as the problem and form… view at source ↗
Figure 6
Figure 6. Figure 6: Results of DVLP with load energizations at Bus 5. [PITH_FULL_IMAGE:figures/full_fig_p014_6.png] view at source ↗
read the original abstract

Due to recent blackout and system split incidents in power grids worldwide, as well as increased system complexity in view of the energy transition, there has been increasing interest in re-evaluating existing Power System Restoration (PSR) plans. In restoration scenarios, due to low island inertia, it is necessary to ensure not only the static, but also the dynamic stability of the system. In this paper, we pose and solve a formulation of the optimal PSR problem including frequency dynamics. We validate the switching constraints for global optimality within a static version of the formulation using a brute-force tree search method. We apply the dynamic problem formulation to the IEEE 9-Bus model, and show that the optimal switching sequence using the static formulation would violate dynamic constraints, illustrating the importance of dynamic considerations in PSR planning.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The paper formulates an OPF-based optimization problem for power system restoration (PSR) that incorporates frequency dynamics to ensure both static and dynamic stability under low-inertia conditions. It validates global optimality of switching constraints in the static formulation via brute-force tree search and applies the dynamic version to the IEEE 9-bus test system, demonstrating that the static optimum violates frequency constraints and thereby underscoring the need for dynamic considerations in PSR planning.

Significance. If the central demonstration holds, the work is significant for highlighting how static-only PSR plans can fail under realistic frequency transients, especially relevant to grids with high renewable penetration. The brute-force validation of the static case supplies an independent optimality check that strengthens the static-dynamic comparison.

major comments (2)
  1. [§4 and §5] §4 (Dynamic Formulation) and §5 (IEEE 9-bus Results): The claim that the static optimal sequence violates dynamic constraints rests on the fidelity of the frequency model (swing equations, inertia constants, governor response, and chosen nadir/ROCOF limits). The manuscript applies a standard simplified generator model to the IEEE 9-bus case but provides no cross-validation against EMT-level simulation or measured restoration transients; any mismatch in effective inertia or load-frequency dependence during islanding could alter the reported violations.
  2. [§3] §3 (Static Validation): While brute-force tree search is used to confirm global optimality of the static switching constraints, the manuscript does not report the size of the search space explored, the number of feasible sequences evaluated, or quantitative error metrics between the OPF solution and the brute-force optimum.
minor comments (2)
  1. Notation for frequency deviation and ROCOF limits should be defined explicitly at first use and kept consistent across equations and figures.
  2. Figure captions for the IEEE 9-bus results should include the exact frequency trajectory values at the reported violation points rather than qualitative statements.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for the constructive comments on our manuscript. We address each major comment below and indicate the revisions made.

read point-by-point responses

  1. Referee: [§4 and §5] §4 (Dynamic Formulation) and §5 (IEEE 9-bus Results): The claim that the static optimal sequence violates dynamic constraints rests on the fidelity of the frequency model (swing equations, inertia constants, governor response, and chosen nadir/ROCOF limits). The manuscript applies a standard simplified generator model to the IEEE 9-bus case but provides no cross-validation against EMT-level simulation or measured restoration transients; any mismatch in effective inertia or load-frequency dependence during islanding could alter the reported violations.

    Authors: We agree that the simplified swing-equation model with governor dynamics is an approximation whose quantitative predictions depend on the chosen parameters. This model is standard in frequency-stability studies of low-inertia restoration (we cite representative references in the revision). In the revised Section 4 we have expanded the discussion of modeling assumptions, including the treatment of load-frequency dependence and the selection of nadir/ROCOF thresholds, and we explicitly note the absence of EMT-level cross-validation. While a full EMT comparison would strengthen the numerical results, the central qualitative finding—that a statically optimal sequence can violate frequency limits—remains robust because the model captures the dominant inertial and primary-frequency-response dynamics that drive the reported violations. revision: partial

  2. Referee: [§3] §3 (Static Validation): While brute-force tree search is used to confirm global optimality of the static switching constraints, the manuscript does not report the size of the search space explored, the number of feasible sequences evaluated, or quantitative error metrics between the OPF solution and the brute-force optimum.

    Authors: We thank the referee for this observation. The revised Section 3 now reports that the brute-force enumeration considered all 256 possible binary switching sequences for the eight controllable elements in the IEEE 9-bus restoration problem. Of these, 42 sequences satisfied the static feasibility constraints. The OPF solution coincides exactly with the best brute-force solution: the objective values differ by less than 1e-6 (numerical tolerance) and the switching decisions are identical, confirming global optimality of the reported static solution. revision: yes

standing simulated objections not resolved
  • Full EMT-level or measurement-based cross-validation of the frequency transients for the specific restoration sequence, which would require simulation tools and data sets not available within the scope of the present study.

Circularity Check

0 steps flagged

No circularity; formulation, brute-force validation, and dynamic application are independent

full rationale

The paper defines a static OPF-based PSR formulation, validates its switching constraints for global optimality via an independent brute-force tree search on the IEEE 9-bus case, then applies a separate dynamic formulation (incorporating frequency dynamics) to the same model and demonstrates that the static optimum violates dynamic limits. No quoted step reduces a claimed prediction or result to a fitted parameter, self-citation, or definitional equivalence by construction. The central illustration relies on external model application and brute-force enumeration rather than any load-bearing self-reference or renaming of inputs.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Abstract-only review provides no explicit free parameters, axioms, or invented entities; the frequency dynamics model itself is an unstated modeling choice whose accuracy is the load-bearing assumption.

pith-pipeline@v0.9.0 · 5668 in / 1175 out tokens · 36316 ms · 2026-05-18T00:24:50.612807+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

39 extracted references · 39 canonical work pages

  1. [1]

    (2025, 05) ENTSO-E expert panel initiates the investigation into the causes of Iberian blackout

    ENTSO-E. (2025, 05) ENTSO-E expert panel initiates the investigation into the causes of Iberian blackout. [Online]. Available: https://www.entsoe.eu/news/2025/05/09/entso-e-expert- panel-initiates-the-investigation-into-the-causes-of-iberian-blackout/

    work page 2025

  2. [2]

    A distributed black-start optimization method for global transmission and distribution network,

    J. Zhao, Q. Zhang, Z. Liu, and X. Wu, “A distributed black-start optimization method for global transmission and distribution network,”IEEE Transactions on Power Systems, vol. 36, no. 5, 2021

    work page 2021

  3. [3]

    Participation of DERs to the bottom-up power system frequency restoration processes,

    F. Sanniti, R. Benato, and F. Milano, “Participation of DERs to the bottom-up power system frequency restoration processes,”IEEE Transactions on Power Systems, vol. 38, no. 3, 2023

    work page 2023

  4. [4]

    Borghetti, C

    A. Borghetti, C. Nucci, and M. Paolone,Handbook of Electrical Power System Dynamics. John Wiley & Sons, Ltd., 2013, ch. 14, pp. 864–899

    work page 2013

  5. [5]

    Power system restoration - the second task force report,

    M. M. Adibi, J. N. Borkoski, and R. J. Kafka, “Power system restoration - the second task force report,”IEEE Transactions on Power Systems, vol. 2, no. 4, pp. 927–932, 1987

    work page 1987

  6. [6]

    System Restoration Procedure and Practices,

    B. Li et.al., “System Restoration Procedure and Practices,” CIGRE WG C2.23, Tech. Rep. 712, December 2017

    work page 2017

  7. [7]

    Dynamic security con- strained AC optimal power flow for microgrids,

    W. S. Alvarez, J. C. L´ opez, F. W. Liederer, D. Dotta, and M. J. Rider, “Dynamic security con- strained AC optimal power flow for microgrids,”Electric Power Systems Research, vol. 236, 2024

    work page 2024

  8. [8]

    A unified framework for frequency control and con- gestion management,

    C. Zhao, E. Mallada, S. Low, and J. Bialek, “A unified framework for frequency control and con- gestion management,” in2016 Power Systems Computation Conference (PSCC), 2016

    work page 2016

  9. [9]

    Development of a black start decision supporting system for isolated power systems,

    Y. T. Chou, C. W. Liu, Y. J. Wang, C. C. Wu, and C. C. Lin, “Development of a black start decision supporting system for isolated power systems,”IEEE Transactions on Power Systems, vol. 28, no. 3, 2013

    work page 2013

  10. [10]

    A post-event generator start-up strategy for renewable penetrated transmission system considering dynamic frequency regulation,

    L. Jiang, Z. Bie, T. Long, X. Wang, H. Xie, and G. Li, “A post-event generator start-up strategy for renewable penetrated transmission system considering dynamic frequency regulation,”IEEE Transactions on Sustainable Energy, vol. 14, no. 2, 2023

    work page 2023

  11. [11]

    Kundur,Power System Stability and Control

    P. Kundur,Power System Stability and Control. McGraw-Hill, 1994

    work page 1994

  12. [12]

    P. Bonami. Handling non-linearities with Gurobi. Gurobi Optimization. [Online]. Available: https://www.gurobi.com/wp-content/uploads/gurobiLiveBarcelona2023-thurs-5-Gurobi-pierre.pdf

    work page

  13. [13]

    Optimal black start allocation for power system restoration,

    G. Patsakis, D. Rajan, I. Aravena, J. Rios, and S. Oren, “Optimal black start allocation for power system restoration,”IEEE Transactions on Power Systems, vol. 33, no. 6, 2018. 14

    work page 2018

  14. [14]

    Optimal generator start-up sequence for bulk system restoration with active distribution networks,

    J. Zhao, H. Wang, Q. Wu, N. D. Hatziargyriou, and F. Shen, “Optimal generator start-up sequence for bulk system restoration with active distribution networks,”IEEE Transactions on Power Sys- tems, vol. 36, no. 3, 2021

    work page 2021

  15. [15]

    Availability assessment based case-sensitive power system restoration strategy,

    W. Liu, J. Zhan, C. Y. Chung, and L. Sun, “Availability assessment based case-sensitive power system restoration strategy,”IEEE Transactions on Power Systems, vol. 35, no. 2, 2020

    work page 2020

  16. [16]

    An optimal power flow function to aid restoration studies of long transmission segments,

    E. M. Viana, E. J. de Oliveira, N. Martins, J. L. R. Pereira, and L. W. de Oliveira, “An optimal power flow function to aid restoration studies of long transmission segments,”IEEE Transactions on Power Systems, vol. 28, no. 1, 2013

    work page 2013

  17. [17]

    A hierarchical response-based approach to the load restoration prob- lem,

    A. Gholami and F. Aminifar, “A hierarchical response-based approach to the load restoration prob- lem,”IEEE Transactions on Smart Grid, vol. 8, no. 4, 2017

    work page 2017

  18. [18]

    Incorporating wind energy in power system restoration planning,

    A. Golshani, W. Sun, Q. Zhou, Q. P. Zheng, and Y. Hou, “Incorporating wind energy in power system restoration planning,”IEEE Transactions on Smart Grid, vol. 10, no. 1, 2019

    work page 2019

  19. [19]

    Rolling optimization of transmission network recovery and load restoration considering hybrid wind-storage system and cold load pickup,

    L. Sun, W. Liu, C. Chung, M. Ding, and J. Ding, “Rolling optimization of transmission network recovery and load restoration considering hybrid wind-storage system and cold load pickup,”Inter- national Journal of Electrical Power & Energy Systems, vol. 141, 2022

    work page 2022

  20. [20]

    An improved two-stage optimization for network and load recovery during power system restoration,

    S. Liao, W. Yao, X. Han, J. Fang, X. Ai, J. Wen, and H. He, “An improved two-stage optimization for network and load recovery during power system restoration,”Applied Energy, vol. 249, 2019

    work page 2019

  21. [21]

    System restoration for low-inertia power systems incorporating fast frequency response via distributionally robust optimization,

    Z. Qin, Y. Li, X. Chen, and H. Liu, “System restoration for low-inertia power systems incorporating fast frequency response via distributionally robust optimization,”IEEE Transactions on Power Systems, vol. 40, no. 3, 2025

    work page 2025

  22. [22]

    Power system restoration with transient stability,

    H. Hijazi, T. W. Mak, and P. V. Hentenryck, “Power system restoration with transient stability,” inProceedings of the Twenty-Ninth AAAI Conference on Artificial Intelligence, 2015

    work page 2015

  23. [23]

    A network reconfiguration approach for power system restoration based on preference-based multiobjective optimization,

    R. Sun, Y. Liu, H. Zhu, R. Azizipanah-Abarghooee, and V. Terzija, “A network reconfiguration approach for power system restoration based on preference-based multiobjective optimization,”Ap- plied Soft Computing, vol. 83, 2019

    work page 2019

  24. [24]

    Multiobjective evolution programming method for feeder reconfiguration,

    Y. T. Hsiao, “Multiobjective evolution programming method for feeder reconfiguration,”IEEE Transactions on Power Systems, vol. 19, 2004

    work page 2004

  25. [25]

    Dynamic frequency-constrained load restoration considering multi-phase cold load pickup behaviors,

    D. Xie, Y. Xu, S. Nadarajan, V. Viswanathan, and A. K. Gupta, “Dynamic frequency-constrained load restoration considering multi-phase cold load pickup behaviors,”IEEE Transactions on Power Systems, vol. 39, no. 1, pp. 107–118, 2024

    work page 2024

  26. [26]

    Distribution system blackstart and restoration using DERs and dynamically formed microgrids,

    S. Maharjan, C. Bai, H. Wang, Y. Yao, F. Ding, and Z. Wang, “Distribution system blackstart and restoration using DERs and dynamically formed microgrids,”IEEE Transactions on Smart Grid, vol. 16, no. 3, pp. 2100–2114, 2025

    work page 2025

  27. [27]

    Anderson and A

    P. Anderson and A. Fouad,Power System Control and Stability. Wiley-IEEE Press, 2003

    work page 2003

  28. [28]

    Stevenson and J

    W. Stevenson and J. Grainger,Power System Analysis. New York: McGraw-Hill, 1994

    work page 1994

  29. [29]

    Yalmip: A toolbox for modeling and optimization in MATLAB,

    J. L¨ ofberg, “Yalmip: A toolbox for modeling and optimization in MATLAB,” inProceedings of the CACSD Conference, Taipei, 2004

    work page 2004

  30. [30]

    Natick, Massachusetts: The MathWorks Inc., 2023

    MATLAB,version R2023b. Natick, Massachusetts: The MathWorks Inc., 2023

    work page 2023

  31. [31]

    Gurobi Optimizer Reference Manual,

    Gurobi Optimization, LLC, “Gurobi Optimizer Reference Manual,” 2024. [Online]. Available: https://www.gurobi.com

    work page 2024

  32. [32]

    Boyd and L

    S. Boyd and L. Vandenberghe,Convex Optimization. Cambridge University Press, 2004

    work page 2004

  33. [33]

    Traversing binary trees simply and cheaply,

    J. M. Morris, “Traversing binary trees simply and cheaply,”Information Processing Letters, vol. 9, no. 5, 1979

    work page 1979

  34. [34]

    Transmission Code 2019,

    Swissgrid, “Transmission Code 2019,” Tech. Rep., May 2020. 15

    work page 2019

  35. [35]

    Review of Existing Hydroelectric Turbine-Governor Simulation Models,

    V. Koritarov et.al., “Review of Existing Hydroelectric Turbine-Governor Simulation Models,” Ar- gonne National Laboratory, Tech. Rep. ANL/DIS-13/05, August 2013

    work page 2013

  36. [36]

    IEEE 9-bus dynamic model,

    D. Virginillo, A. Derviˇ skadi´ c, and M. Paolone, “IEEE 9-bus dynamic model,” https://github.com/DESL-EPFL/ieee-9bus-dynamic-model, 2025

    work page 2025

  37. [37]

    Synchronous machine,

    S. Denneti` ere, J. Mahseredjian, and U. Karaagac, “Synchronous machine,”EMTP User Manual, 2012

    work page 2012

  38. [38]

    Technische Anforderungen an die automatische Frequenzentlastung unter Ber¨ ucksichtigung ver¨ anderter Vorgaben,

    Verband Schweizerischer Elektrizit¨ atsunternehmen, “Technische Anforderungen an die automatische Frequenzentlastung unter Ber¨ ucksichtigung ver¨ anderter Vorgaben,” VSE Netztechnik & Betrieb, Tech. Rep., September 2016

    work page 2016

  39. [39]

    ELIA’s restoration plan - non-confidential version,

    ELIA, “ELIA’s restoration plan - non-confidential version,” Elia Group, Tech. Rep., January 2024. 16

    work page 2024