A Single-Loop Bilevel Deep Learning Method for Optimal Control of Obstacle Problems
Yongcun Song, Shangzhi Zeng, Jin Zhang, Lvgang Zhang
TL;DR
This work tackles the challenging problem of optimal control for obstacle problems, which are nonsmooth and have a bilevel structure. It introduces constraint-embedding neural networks to preserve feasibility for state and control, and a single-loop stochastic bilevel training algorithm (S2-FOBA) built on a Moreau envelope reformulation to avoid nested optimization. A two-stage lower-level refinement ensures numerical feasibility of the state after training. The approach is mesh-free, scalable to high dimensions and complex domains, and yields competitive accuracy with reduced computational cost compared with classical discretization-based methods. The framework is extended to complex geometries and elliptic variational inequalities, highlighting its versatility and potential for broad applicability in PDE-constrained optimization.
Abstract
Optimal control of obstacle problems arises in a wide range of applications and is computationally challenging due to its nonsmoothness, nonlinearity, and bilevel structure. Classical numerical approaches rely on mesh-based discretization and typically require solving a sequence of costly subproblems. In this work, we propose a single-loop bilevel deep learning method, which is mesh-free, scalable to high-dimensional and complex domains, and avoids repeated solution of discretized subproblems. The method employs constraint-embedding neural networks to approximate the state and control and preserves the bilevel structure. To train the neural networks efficiently, we propose a Single-Loop Stochastic First-Order Bilevel Algorithm (S2-FOBA), which eliminates nested optimization and does not rely on restrictive lower-level uniqueness assumptions. We analyze the convergence behavior of S2-FOBA under mild assumptions. Numerical experiments on benchmark examples, including distributed and obstacle control problems with regular and irregular obstacles on complex domains, demonstrate that the proposed method achieves satisfactory accuracy while reducing computational cost compared to classical numerical methods.
