[论文解读] Backstepping Control of the One-Phase Stefan Problem

Shumon Koga, Mamadou Diagne+2 more|arXiv (Cornell University)|Jul 14, 2016

Stability and Controllability of Differential Equations参考文献 14被引用 5

一句话总结

本文提出了一种针对单相Stefan问题的非线性反推边界控制器，该问题为具有移动边界的1D PDE-ODE系统。通过引入一种新颖的非线性反推变换，避免了域重标度，该方法实现了移动界面和温度误差的H¹-范数的指数稳定性，同时相对于初始条件保持了物理约束。

ABSTRACT

In this paper, a backstepping control of the one-phase Stefan Problem, which is a 1-D diffusion Partial Differential Equation (PDE) defined on a time varying spatial domain described by an ordinary differential equation (ODE), is studied. A new nonlinear backstepping transformation for moving boundary problem is utilized to transform the original coupled PDE-ODE system into a target system whose exponential stability is proved. The full-state boundary feedback controller ensures the exponential stability of the moving interface to a reference setpoint and the ${\cal H}_1$-norm of the distributed temperature by a choice of the setpint satisfying given explicit inequality between initial states that guarantees the physical constraints imposed by the melting process.

研究动机与目标

• 解决时间变域中具有移动边界的1D扩散PDE的控制方法缺乏的问题。
• 开发一种基于反推的控制器，确保指数稳定性，且无需域重标度或固定域假设。
• 通过确保控制输入保持为正，保证物理有效性，维持模型完整性。
• 将反推方法扩展至移动边界问题，克服了先前有限维近似方法的局限性。
• 建立参考设定点相对于初始数据的条件，以在熔化/凝固过程中保持物理约束。

提出的方法

• 为移动边界PDE-ODE系统专门设计了一种新型非线性反推变换，避免了将域重标度至固定空间区间的需要。
• 该变换将原始耦合PDE-ODE系统映射到一个具有已知指数稳定性特性的目标系统。
• 控制器由该变换的逆推导得出，确保对移动界面和温度分布的全状态反馈控制。
• 构造了一个李雅普诺夫泛函以证明闭环系统的指数稳定性，其衰减速率依赖于控制器增益和系统参数。
• 通过仔细估计包含温度H¹-范数和界面偏差的能量类范数，将时变域的影响纳入稳定性分析。
• 控制器增益的选择使得控制输入保持为正，从而确保与熔化/凝固动力学的物理一致性。

实验结果

研究问题

• RQ1 能否在不进行域重标度的情况下，将反推控制框架扩展至具有移动边界的PDE-ODE系统？
• RQ2 在单相Stefan问题中，参考设定点相对于初始数据需满足何种条件，才能确保物理可行性和稳定性？
• RQ3 如何实现移动界面和温度误差H¹-范数的指数稳定性？
• RQ4 控制输入的正性在保持相变过程中的物理有效性和稳定性方面起什么作用？
• RQ5 与先前方法相比，该方法在收敛速度和约束处理方面表现如何？

主要发现

• 所提出的反推变换成功实现了在不进行域重标度的情况下，单相Stefan问题中移动界面和温度误差H¹-范数的稳定。
• 通过李雅普诺夫泛函证明了闭环系统的指数稳定性，衰减速率依赖于控制器增益和系统参数。
• 控制器确保了控制输入为正，这对物理一致性至关重要，可避免出现非物理解的熔化或凝固行为。
• 设定点必须满足相对于初始条件的显式不等式（如公式(10)所示），以保持物理约束。
• 数值仿真验证了理论结果：当 $ s_0 = 0.01 $m 时，界面收敛至 $ s_r = 0.35 $m 且 $ ilde{q}(t) > 0 $；当 $ s_0 = 0.25 $m 时，控制器变为负值，与理论条件一致。
• 该方法实现了界面误差与温度误差H¹-范数之和的指数稳定性，优于先前仅实现渐近稳定性或L²-范数稳定性的方法。