Penalty-Free SDDP: Feasibility Cuts for Robust Multi-Stage Stochastic Optimization in Energy Planning
Guilherme Freitas, Luiz Carlos da Costa Junior, Tiago Andrade, Alexandre Street
TL;DR
This paper introduces Penalty-Free SDDP, a general extension of stochastic dual dynamic programming that eliminates artificial penalties for infeasibilities by incorporating a Future Feasibility Function and a dual-stage recursion. Feasibility information is handled via a union of feasibility cuts, separating constraint satisfaction from economic optimization while preserving SDDP's decomposition and convergence properties. The approach is validated on a large-scale hydrothermal energy planning problem, demonstrating equivalent feasibility to a benchmark problem and improved interpretability by clearly identifying infeasible conditions without penalty calibration. The results suggest robust applicability to complex multi-stage energy planning and provide a foundation for extending to fully stochastic settings.
Abstract
Multi-stage decision problems under uncertainty can be efficiently solved with the Stochastic Dual Dynamic Programming (SDDP) algorithm. However, traditional implementations require all stage problems to be feasible. Feasibility is usually enforced by adding slack variables and penalizing them in the objective function, a process that depends on case-specific calibration and often distorts the economic interpretation of results. This paper proposes the Penalty-Free SDDP, an extension that introduces a Future Feasibility Function alongside the traditional Future Cost Function. The new recursion handles infeasibilities automatically, distinguishing between temporary and truly infeasible cases, and propagates feasibility information across stages through dedicated feasibility cuts. The approach was validated in a large-scale deterministic case inspired by the Brazilian hydrothermal system, achieving equivalent feasibility to the benchmark solution while eliminating miscalibrated artificial penalties. Results confirm its robustness and practicality as a foundation for future stochastic, multi-stage applications.
