REVIEW 2 major objections 5 minor 32 references
Optimizing Flexibility in Power Systems by Maximizing the Region of Manageable Uncertainties
T0 review · 2 major / 5 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read The paper establishes that a power grid's flexibility—the range of uncertain loads and generations it can safely absorb—can be maximized by solving an existence-constrained semi-infinite optimization problem that returns rigorous bounds…
desk verdict A credible extension of semi-infinite optimization to power-system flexibility maximization, with two caveats: an unverified existence assumption and 5% solver gaps that weaken the claim of rigorous global bounds. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the flexibility index $\delta$, which scales a user-chosen parametric uncertainty region $T(\delta,x)$; an existence-constrained semi-infinite program is an optimization problem with a constraint that must hold for every uncertainty value and, for each value, asserts the existence of a control. The paper's Proposition 1 is the identity that carries the argument: the existential constraint over the variable region $T(\delta,x)$ is replaced by a min-of-two-objectives constraint over the fixed host set $\bar Y$, making the problem solvable by standard discretization. At each iteration the algorithm solves a relaxed mixed-integer linear master problem and a maxmin worst-case problem that adds the most violating uncertainty point to the discretization; two specializations drop redundant discretization points and transform retained points to the boundary of the shrinking hyperbox. The numerical engine is the uniqueness and continuity of the implicit grid-state function $s(x,y,z)$ together with the load-distribution model, which fixes the generator control offsets uniquely once the uncertainty and set-points are given.
What would settle it
Take a returned solution $(\delta,x)$ and evaluate the worst-case generation problem over the full host set; if any point $y$ inside $T(\delta,x)$ has no feasible control $z$ satisfying $g(x,y,z)\le0$, the lower-bound certificate is wrong. A simpler direct check on the motivating example is to enumerate the polyhedron of manageable injections shown in the paper and verify that the returned hyperbox is truly contained in it.
Extended reading notes
Core claim
The central discovery is that the maximal manageable uncertainty region can be optimized, not merely evaluated, by turning the problem into an existence-constrained semi-infinite program and solving it with an adaptive discretization scheme. The key step is a reformulation: when the uncertainty region has the form $T(x,\delta)=\{y\in\bar Y:h(x,y)\le\delta\}$ with continuous $h$, the constraint that every $y$ in $T(x,\delta)$ admits a feasible control is equivalent, up to an arbitrarily small objective error, to requiring for all $y$ in the fixed host set $\bar Y$ that $\min_{z\in Z(x,y)}\min\{\alpha(\delta-h(x,y)),g(x,y,z)\}\le0$. This removes the dependence of the uncertainty set on the preventive actions and lets the algorithm iterate between a relaxed master problem and a worst-case generation problem. For the power-system model, which uses the DC flow approximation, piecewise-linear phase shifters, bus merging, and a unique implicit load-distribution response, each subproblem is a mixed-integer linear program. The result is a rigorous interval for the flexibility index and guaranteed preventive actions for the returned uncertainty region.
Load-bearing premise
Everything rests on the assumption that, for every preventive action and every uncertainty in the host set, the grid state is uniquely determined and at least one feasible control action exists; the paper assumes this rather than proving it for all iterates.
Editorial extensions
If this is right
- With the scaled-hyperbox parameterization, the returned $\delta$ is a guaranteed inner approximation: every uncertainty inside the box $[y_0-\Delta^-\delta,\,y_0+\Delta^+\delta]$ can be managed, and no larger scaled box can be certified.
- With the transfer parameterization, the result guarantees that every additional power transfer between two regions up to the returned capacity is safe, not merely that the maximal transfer is safe.
- For fixed preventive actions, the auxiliary flexibility-index problem can be solved on its own, giving a flexibility measure for an already chosen operating point.
- The method can in principle be adapted to any continuous parameterization $h$ of the uncertainty region, which the medium-scale experiments exploit by comparing the hyperbox and transfer formulations.
- On the medium-scale French-network instance, the transfer formulation returns approximately $[4530,4760]$ MW of guaranteed France-to-Spain additional transfer capacity in about 300 seconds; the hyperbox formulation returns $\delta\in[0.546,0.573]$ in about 83 minutes.
Reading between the lines
- The strong dependence of run time on the scaling parameter $\alpha$, which the paper documents but does not resolve, suggests an adaptive $\alpha$-update rule as a natural extension; the authors leave the choice of a systematic value open.
- Because Proposition 1 only needs $h$ to be continuous, the same discretization could be tested on rotated ellipsoids, zonotopes, or other convex bodies, giving less conservative approximations of the manageable region than an axis-aligned hyperbox.
- The paper's stated assumption that grid states are uniquely determined by $(x,y,z)$ could be checked numerically on returned iterates by computing the implicit function residual; if residuals are not small, the certified bounds would not be valid for the actual model.
- A direct practical extension is to treat $\alpha$ and the integrality tolerance as tunable hyperparameters on a benchmark set, since the observed time-outs indicate that solver settings materially change which instances can be solved.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes a method for maximizing a scalar flexibility index δ that measures the size of a parametrized uncertainty region T(δ, x) that a DC-power-flow grid can manage in the worst case through preventive actions x and corrective controls z, including mixed-integer controls such as bus merging. The problem is formulated as an existence-constrained semi-infinite optimization problem, and the authors prove an exact reformulation (Proposition 1) that replaces the x-dependent uncertainty set with a fixed host set, then specialize the discretization algorithm of Djelassi and Mitsos (2021). The method is tested on a two-generator motivating example, a modified 30-bus IEEE system, and a 6,800-node RTE-based instance, for two parameterizations: a scaled hyperbox of uncertain injections and an additional net power-transfer capacity between two regions. The reported results are bracketed intervals, e.g., δ ∈ [0.925, 0.954] and δ ∈ [0.590, 0.614] on the small instance, and δ ∈ [4530, 4760] MW and δ ∈ [0.546, 0.573] on the medium instance.
Significance. If the algorithmic guarantees are valid, the paper represents a useful methodological contribution: it connects the process-systems flexibility index to existence-constrained semi-infinite optimization in a power-system setting, accommodates mixed-integer controls, and demonstrates computational tractability on a medium-scale DC-flow instance. The reformulation in Proposition 1 is a helpful device, and the transformation-based specialization of Section 6.3.2 is a practical enhancement. The numerical section is transparent about wall-clock variability and timeouts, and the medium-scale instance data are made available. The central caveat is that the rigorous validity of the reported bounds depends on a nonempty-control feasibility condition that is not verified for all reported instances; this is the main issue that must be resolved before the paper's claims are fully supported.
major comments (2)
- [Section 4 and Section 5.3.2, Eq. (9)] The global-validity guarantee stated in the abstract and Section 6.2 is conditional on the assumption, made in Section 4, that Z(x, y) is nonempty on X × Ybar. In the concrete DC model, this reduces to the load-distribution feasibility condition (9). The paper neither includes (9) as a constraint in the optimization problems (13) and (14) nor verifies it over the full host set for the transfer-capacity parametrization of Section 7.2. For the hyperbox parametrization, Section 8.1 chooses δ_UB to guarantee (9), and Section 8.2 states that generator bounds were increased 'to avoid a violation of (9)', but no such verification is reported for the medium-scale instance of Section 8.3. If (9) fails for some y ∈ Ybar, the expression min_{z ∈ Z(x,y)} ... is not well defined, the worst-case generation subproblem (12f) can encounter infeasible inner problems, and the intervals δ ∈ [4530, 4760] MW and δ ∈ [0.546, 0.573] are not rigorously justified. Please either impose (9) explicitly, verify it for all (x, y) encountered, or restrict Ybar and state the resulting limitation.
- [Section 6.1, Eq. (12g)] The displayed equivalence '0 ≥ sup_y min_z min{α(δ−h), g} ⇔ ∀y∈Ybar [∃z∈Z(x,y): g(x,y,z) ≤ 0]' is not correct. The right-hand side should be '∀y∈Ybar [∃z∈Z(x,y): min{α(δ−h(x,y)), g(x,y,z)} ≤ 0]' (equivalently, with the min expression inside the quantifier). As printed, the proof asserts a strictly stronger statement, since it would require g ≤ 0 for every y in the host set, which is not equivalent to the original constraint. The proposition statement is correct, so this is a fixable proof error, but it must be corrected before the central reformulation can be considered proved.
minor comments (5)
- [Section 6.1, Proposition 1 statement] The displayed reformulation in the proposition statement is missing the '≤ 0' after the min expression; compare with equation (12e).
- [Section 6.3.2, Eq. (14)] The transformation divides by h(y_d); the case h(y_d) = 0 (for example y_d = y0) is not discussed, and a non-zero assumption should be stated explicitly.
- [Section 8.1 and Fig. 5] The reported δ ∈ [-1.857, -1.857] is the algorithm's bracket; it would be clearer to state explicitly that this is a computed interval up to the chosen tolerances, and to explain the negative sign convention (the objective is -δ).
- [Section 7.2, Eq. (19)] The role of the αh(x, y) term inside the constraint g is not immediately clear; a sentence explaining that this encodes the min-form of the Proposition 1 reformulation would improve readability.
- [Tables 2 and 3] The paper reports ranges over three repetitions but does not state whether the repetitions differ in random seeds or solver settings; adding this information would help readers interpret the variability.
Circularity Check
No significant circularity: the flexibility bounds are produced by solving the stated optimization, and the cited prior algorithms are independent support rather than inputs that force the results.
full rationale
The paper's central derivation is self-contained. The flexibility index δ is the decision variable of the optimization, not a fitted parameter; the reported intervals [δ_LB, δ_UB] are computed by lower- and upper-bounding MILP subproblems. Proposition 1 is proved in Eqs. (12a)-(12g), including the exactness argument (δ−ε), so the reformulation to a fixed uncertainty set does not assume the conclusion. The discretization algorithm and grid model are taken from the authors' prior work (Djelassi and Mitsos 2021; Djelassi et al. 2018), but those citations supply an independently published algorithmic framework and model equations; they are not used to forbid alternatives or to define the present result into existence. No quantity is fitted and then renamed as a prediction, and no ansatz is smuggled in through citation: the PST piecewise-linear model and load-distribution equation (8) are explicit modeling assumptions. The main caveat is validity, not circularity: Section 4 assumes Z(x,y) nonempty and states unique on X×Ybar, and Section 5 notes this requires (9); the paper verifies (9) for the hyperbox parameterization via δ_UB but not for the transfer-capacity case. If (9) fails for some y, the inner control offsets do not exist and the claimed bounds would not be valid for those points. That is an unverified assumption affecting correctness, not a case where the derivation reduces to its own inputs.
Assumptions & free parameters
free parameters (3)
- Alpha scaling constant in the exact reformulation =
alpha' = 0.5 (small instance), alpha' = 10 (medium instance), alpha = alpha' / delta_norm
- Initial restriction epsilon_R^0 =
0.05 for the main algorithm, 0.005 for the auxiliary problem
- Relative optimality tolerance =
0.05 overall, 0.025 auxiliary
assumptions (5)
- domain assumption DC power flow approximation is an adequate model of grid operation for this analysis (affine line flows, Eq. (2))
- domain assumption Grid state is uniquely given by a continuous implicit function s(x,y,z) and Z(x,y) is nonempty on X times Ybar
- domain assumption The load-distribution control via the mid-function (8) has a unique solution whenever condition (9) holds
- domain assumption MILP subproblems are solved to global optimality within stated relative gaps by Gurobi 10.0.2
- standard math Convergence and finite termination of the adaptive discretization algorithm (Blankenship-Falk 1976; Mitsos 2011)
Cite this review
Pith. "Pith review of Optimizing Flexibility in Power Systems by Maximizing the Region of Manageable Uncertainties." pith.science (2026). https://pith.science/paper/VE2LTRGX
@misc{pith2026241118178,
author = {Pith},
title = {Pith review of: Optimizing Flexibility in Power Systems by Maximizing the Region of Manageable Uncertainties},
year = {2026},
howpublished = {\url{https://pith.science/paper/VE2LTRGX}},
note = {Machine review of arXiv:2411.18178}
}
read the original abstract
Motivated by the increasing need to hedge against load and generation uncertainty in the operation of power grids, we propose flexibility maximization during operation. We consider flexibility explicitly as the amount of uncertainty that can be handled while still ensuring nominal grid operation in the worst-case. We apply the proposed flexibility optimization in the context of a DC flow approximation. By using a corresponding parameterization, we can find the maximal range of uncertainty and a range for the manageable power transfer between two parts of a network subject to uncertainty. We formulate the corresponding optimization problem as an (existence-constrained) semi-infinite optimization problem and specialize an existing algorithm for its solution.
Figures
Figures from the paper (3 more)
Reference graph
Works this paper leans on
-
[1]
Akrami A, Doostizadeh M, Aminifar F (2019) Power system flexibility: an overview of emergence to evolution. Journal of Modern Power Systems and Clean Energy 7(5):987--1007, doi:https://doi.org/10.1007/s40565-019-0527-4
-
[2]
Energy Reports 6:101--106, doi:https://doi.org/10.1016/j.egyr.2019.11.048
Babatunde OM, Munda JL, Hamam Y (2020) Power system flexibility: A review. Energy Reports 6:101--106, doi:https://doi.org/10.1016/j.egyr.2019.11.048
-
[3]
Blankenship JW, Falk JE (1976) Infinitely constrained optimization problems. Journal of Optimization Theory and Applications 19(2):261--281, doi:https://doi.org/10.1007/bf00934096
-
[4]
In: 2015 IEEE Eindhoven PowerTech, IEEE, IEEE, pp 1--6, doi:https://doi.org/10.1109/ptc.2015.7232514
Bucher MA, Delikaraoglou S, Heussen K, Pinson P, Andersson G (2015) On quantification of flexibility in power systems. In: 2015 IEEE Eindhoven PowerTech, IEEE, IEEE, pp 1--6, doi:https://doi.org/10.1109/ptc.2015.7232514
-
[5]
Capitanescu F (2021) Power system flexibility region under uncertainty with respect to congestion and voltage constraints. In: 2021 IEEE Madrid PowerTech, IEEE, pp 1--6, doi:https://doi.org/10.1109/powertech46648.2021.9494935
arXiv 2021
-
[6]
Djelassi H (2020) D iscretization-based algorithms for the global solution of hierarchical programs. Dissertation, Rheinisch-Westfälische Technische Hochschule Aachen, doi:https://doi.org/10.18154/RWTH-2020-09163
-
[7]
Journal of Global Optimization 68(2):227--253, doi:https://doi.org/10.1007/s10898-016-0476-7
Djelassi H, Mitsos A (2017) A hybrid discretization algorithm with guaranteed feasibility for the global solution of semi-infinite programs. Journal of Global Optimization 68(2):227--253, doi:https://doi.org/10.1007/s10898-016-0476-7
-
[8]
Djelassi H, Mitsos A (2021) Global solution of semi-infinite programs with existence constraints. Journal of Optimization Theory and Applications 188(3):863--881, doi:https://doi.org/10.1007/s10957-021-01813-2
Show all 32 references
-
[9]
In: 2018 Power Systems Computation Conference (PSCC), IEEE, pp 1--7, doi:https://doi.org/10.23919/pscc.2018.8444136
Djelassi H, Fliscounakis S, Mitsos A, Panciatici P (2018) Hierarchical programming for worst-case analysis of power grids. In: 2018 Power Systems Computation Conference (PSCC), IEEE, pp 1--7, doi:https://doi.org/10.23919/pscc.2018.8444136
2018
-
[10]
Journal of Global Optimization 75(2):341--392, doi:https://doi.org/10.1007/s10898-019-00764-3
Djelassi H, Glass M, Mitsos A (2019) Discretization-based algorithms for generalized semi-infinite and bilevel programs with coupling equality constraints. Journal of Global Optimization 75(2):341--392, doi:https://doi.org/10.1007/s10898-019-00764-3
2019 doi
-
[11]
Naval Research Logistics Quarterly 24(3):441--450, doi:https://doi.org/10.1002/nav.3800240307
Falk JE, Hoffman K (1977) A nonconvex max‐min problem. Naval Research Logistics Quarterly 24(3):441--450, doi:https://doi.org/10.1002/nav.3800240307
1977 doi
-
[12]
In: 2019 IEEE Power & Energy Society General Meeting (PESGM), IEEE, pp 1--5, doi:https://doi.org/10.1109/pesgm40551.2019.8974032
G \'o mez JD, Borda CE (2019) Operational flexibility in power systems—a geometrical view. In: 2019 IEEE Power & Energy Society General Meeting (PESGM), IEEE, pp 1--5, doi:https://doi.org/10.1109/pesgm40551.2019.8974032
2019
-
[13]
Computers & Chemical Engineering 11(6):675--693, doi:https://doi.org/10.1016/0098-1354(87)87011-4
Grossmann IE, Floudas CA (1987) Active constraint strategy for flexibility analysis in chemical processes. Computers & Chemical Engineering 11(6):675--693, doi:https://doi.org/10.1016/0098-1354(87)87011-4
1987 doi
-
[14]
Computers & Chemical Engineering 7(4):439--462, doi:https://doi.org/10.1016/0098-1354(83)80022-2
Grossmann IE, Halemane KP, Swaney RE (1983) Optimization strategies for flexible chemical processes. Computers & Chemical Engineering 7(4):439--462, doi:https://doi.org/10.1016/0098-1354(83)80022-2
1983 doi
-
[15]
OECD, doi:https://doi.org/10.1787/9789264111394-en, ://www.oecd-ilibrary.org/content/publication/9789264111394-en
International Energy Agency (2011) Harnessing Variable Renewables: A Guide to the Balancing Challenge. OECD, doi:https://doi.org/10.1787/9789264111394-en, ://www.oecd-ilibrary.org/content/publication/9789264111394-en
2011
-
[16]
IEEE Transactions on Power Systems 36(6):4953--4966, doi:https://doi.org/10.1109/tpwrs.2021.3075925
Lee D, Turitsyn K, Molzahn DK, Roald LA (2021) Robust ac optimal power flow with robust convex restriction. IEEE Transactions on Power Systems 36(6):4953--4966, doi:https://doi.org/10.1109/tpwrs.2021.3075925
2021
-
[17]
In: 2013 IEEE power & energy society general meeting, IEEE, pp 1--11, doi:https://doi.org/10.1109/tste.2012.2212471
Ma J, Silva V, Belhomme R, Kirschen DS, Ochoa LF (2013) Evaluating and planning flexibility in sustainable power systems. In: 2013 IEEE power & energy society general meeting, IEEE, pp 1--11, doi:https://doi.org/10.1109/tste.2012.2212471
2013
-
[18]
In: 2011 IEEE power and energy society general meeting, IEEE, pp 1--8, doi:https://doi.org/10.1109/pes.2011.6039617
Menemenlis N, Huneault M, Robitaille A (2011) Thoughts on power system flexibility quantification for the short-term horizon. In: 2011 IEEE power and energy society general meeting, IEEE, pp 1--8, doi:https://doi.org/10.1109/pes.2011.6039617
2011
-
[19]
Optimization 60(10-11):1291--1308, doi:https://doi.org/10.1080/02331934.2010.527970
Mitsos A (2011) Global optimization of semi-infinite programs via restriction of the right-hand side. Optimization 60(10-11):1291--1308, doi:https://doi.org/10.1080/02331934.2010.527970
2011
-
[20]
Journal of Global Optimization 61(1):1--17, doi:https://doi.org/10.1007/s10898-014-0146-6
Mitsos A, Tsoukalas A (2015) Global optimization of generalized semi-infinite programs via restriction of the right hand side. Journal of Global Optimization 61(1):1--17, doi:https://doi.org/10.1007/s10898-014-0146-6
2015 doi
-
[21]
IEEE Transactions on Sustainable Energy 6(3):800--809, doi:https://doi.org/10.1109/tste.2015.2410760
Nosair H, Bouffard F (2015) Flexibility envelopes for power system operational planning. IEEE Transactions on Sustainable Energy 6(3):800--809, doi:https://doi.org/10.1109/tste.2015.2410760
2015
-
[22]
European Journal of Operational Research 119(2):301--313, doi:https://doi.org/10.1016/s0377-2217(99)00132-0
Still G (1999) Generalized semi-infinite programming: Theory and methods. European Journal of Operational Research 119(2):301--313, doi:https://doi.org/10.1016/s0377-2217(99)00132-0
1999 doi
-
[23]
European Transactions on Electrical Power 22(6):846--865, doi:https://doi.org/10.1002/etep.617
Wei-qing S, Cheng-min W, Yan Z (2012) Flexibility evaluation and flexible comprehensive optimization in power systems. European Transactions on Electrical Power 22(6):846--865, doi:https://doi.org/10.1002/etep.617
2012 doi
-
[24]
Computers & Chemical Engineering 166:107969, doi:https://doi.org/10.1016/j.compchemeng.2022.107969
Zhao F, Ochoa MP, Grossmann IE, Garc \' a-Mu \ n oz S, Stamatis SD (2022) Novel formulations of flexibility index and design centering for design space definition. Computers & Chemical Engineering 166:107969, doi:https://doi.org/10.1016/j.compchemeng.2022.107969
2022
-
[25]
IEEE Transactions on Power Systems 31(1):339--347, doi:https://doi.org/10.1109/TPWRS.2015.2390038
Zhao J, Zheng T, Litvinov E (2015) A unified framework for defining and measuring flexibility in power system. IEEE Transactions on Power Systems 31(1):339--347, doi:https://doi.org/10.1109/TPWRS.2015.2390038
2015
-
[26]
://publications.rwth-aachen.de/record/828516
Zingler A, Fliscounakis S (2021) M edium-scale power grid instance with 6800 nodes for DC -flow approximation. ://publications.rwth-aachen.de/record/828516
2021
-
[27]
, " * write output.state after.block = add.period write newline
ENTRY address archive author booktitle chapter doi edition editor eid eprint howpublished institution journal key month note number organization pages publisher school series title type url volume year label extra.label sort.label short.list INTEGERS output.state before.all mi...
-
[28]
write newline
" write newline "" before.all 'output.state := FUNCTION add.period duplicate empty 'skip "." * add.blank if FUNCTION if.digit duplicate "0" = swap duplicate "1" = swap duplicate "2" = swap duplicate "3" = swap duplicate "4" = swap duplicate "5" = swap duplicate "6" = swap dupl...
-
[29]
, " * write output.state after.block = add.period write newline
ENTRY address author booktitle chapter doi edition editor eid howpublished institution journal key month note number organization pages publisher school series title type url volume year label INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION in...
-
[30]
write newline
" write newline "" before.all 'output.state := FUNCTION if.digit duplicate "0" = swap duplicate "1" = swap duplicate "2" = swap duplicate "3" = swap duplicate "4" = swap duplicate "5" = swap duplicate "6" = swap duplicate "7" = swap duplicate "8" = swap "9" = or or or or or or...
-
[31]
, " * write output.state after.block = add.period write newline
ENTRY address author booktitle chapter doi edition editor eid howpublished institution journal key month note number organization pages publisher school series title type url volume year label INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION in...
-
[32]
write newline
" write newline "" before.all 'output.state := FUNCTION if.digit duplicate "0" = swap duplicate "1" = swap duplicate "2" = swap duplicate "3" = swap duplicate "4" = swap duplicate "5" = swap duplicate "6" = swap duplicate "7" = swap duplicate "8" = swap "9" = or or or or or or...
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.