[论文解读] Continuation method for PDE-constrained global optimization: Analysis and application to the shallow water equations

Jorn Baayen, Teresa Piovesan+1 more|arXiv (Cornell University)|Jan 19, 2018

Advanced Control Systems Optimization参考文献 27被引用 5

一句话总结

本文提出一种同伦延续方法，通过路径稳定、无分岔的同伦路径，将凸松弛问题变形为非凸问题，以全局最优解求解PDE约束的全局优化问题。该方法通过嵌套同伦保持凸性特性，确保全局最优性，从而实现对河网等大规模系统（由浅水方程控制）的高效闭环模型预测控制。

ABSTRACT

This paper shows how a class of non-convex optimization problems constrained by discretized nonlinear partial differential equations may be solved to global optimality using an interior point continuation method. The solution procedure rests on a nested homotopy. The inner homotopy solves a barrier problem by driving the barrier parameter to zero. The outer homotopy deforms a convex relaxation to the original non-convex problem in a way that stays clear of bifurcations. A requirement for global optimality is that the objective is convex and that the search space remains path-connected. As a case study, a class of real-world optimization problems subject to the shallow water equations is analyzed. A benchmark as well as a practical implementation demonstrate that the approach is suitable for closed-loop non-convex model predictive control of large-scale cyber-physical systems.

研究动机与目标

• 解决大规模、非凸、PDE约束优化问题在实时控制场景中达到全局最优的挑战。
• 克服传统内点方法在非凸问题中无法保证全局最优性的局限。
• 开发一种计算高效的替代方法，以替代多项式层次方法，后者因复杂度高而难以应用于大规模闭环控制。
• 建立局部最优解即为全局最优解的充分条件，通过同伦延续方法实现。
• 通过荷兰水管理系统的浅水方程实际案例研究，展示方法的实际适用性。

提出的方法

• 构建嵌套同伦框架：内层同伦将障碍参数驱动至零以求解障碍问题，外层同伦将凸松弛问题变形为原始非凸问题。
• 外层同伦使用变形参数 θ ∈ [0,1]，其中 θ=0 对应凸松弛问题（保证全局最优解），θ=1 对应原始非凸问题。
• 该方法依赖两个新概念：零凸性与路径稳定性，以确保解路径连续且无分岔。
• 从 θ=0 处的凸解出发，连续追踪唯一解路径 f(θ)，直至 θ=1 处的最终解，同时保持全局最优性。
• 在每个 θ 步长上应用内点方法求解障碍问题，确保原始与对偶可行性，并保持变量的严格正性。
• 该方法避免增加变量数量（与Lasserre层次方法不同），使其可扩展至最多 500,000 个变量的问题。

实验结果

研究问题

• RQ1 在何种条件下，可通过基于同伦的方法对非凸、PDE约束优化问题实现全局最优解？
• RQ2 如何形式化定义路径稳定性和零凸性，并利用其确保解路径中无分岔？
• RQ3 同伦延续方法能否高效实现于大规模信息物理系统实时闭环模型预测控制？
• RQ4 与启发式或全局优化方法相比，该方法在计算效率和解质量方面对水力控制问题的性能提升程度如何？
• RQ5 如何在保持路径稳定性的同时，将非线性水力元件（如可调堰、水轮机）整合进优化框架？

主要发现

• 所提出的同伦方法通过从凸松弛到非凸问题的连续、无分岔路径，确保全局最优性，前提是目标函数为凸且问题具有路径稳定性。
• 该方法在不增加变量数量的前提下，实现对非凸PDE约束问题的全局优化，优于多项式层次方法。
• 基准测试显示，该方法可扩展至最多 500,000 个变量，证明其适用于大规模系统。
• 该方法成功应用于荷兰一个370公里运河系统的实时模型预测控制，通过四个泵站实现水位与水质的联合调控。
• 在一台MacBook Pro上，10个节点、72个时间步的示例问题在约0.4秒内求解，证明其具备实时可行性。
• 该方法支持在优化框架内对非线性水力结构（如可调堰、水轮机）进行建模，实现复杂水管理系统中的实际部署。