Safe Trajectory Tracking of the Stefan Problem with Second-Order Moving Boundary Dynamics

技术摘要：具有二阶移动边界动力学的 Stefan 问题的安全轨迹跟踪

问题陈述
本文解决了一相 Stefan 问题的安全轨迹跟踪控制，该问题是描述熔化和凝固等热相变现象的经典模型。所针对的具体挑战是控制一个移动边界（相界面）由二阶动力学（代表热惯性）而非更常见的一阶动力学支配的系统。该系统由定义在时变空间域上的温度分布的抛物型偏微分方程（PDE）建模，并与描述边界位置的二阶常微分方程（ODE）耦合。

控制目标是设计一个边界热通量输入 $q_c(t)$，驱动系统的移动边界和温度分布跟踪期望的参考轨迹 $s_r(t)$ 和 $T^r(x,t)$。关键在于，控制器必须确保安全性，即模型的物理约束必须始终得到满足：温度必须保持在熔点以上（$T(x,t) \ge T_m$），以防止热入口处的液相冻结，且边界位置必须保持在物理域内（$0 < s(t) < L$）。

方法论
所提出的控制策略结合了前馈分量与通过能量整形和 PDE 反步法设计的反馈控制器。

1. 前馈控制（逆 Stefan 问题）：
   为了生成参考轨迹，作者针对给定的参考边界运动 $s_r(t)$ 求解逆 Stefan 问题。他们采用级数展开方法，将参考温度分布表示为相对于移动边界的空间坐标的幂级数。该级数的系数通过递归确定。本文证明了在假设参考轨迹具有 Gevrey 阶的条件下，该无穷级数的收敛性，从而确保了定义良好的参考温度分布以及相应的前馈热通量 $q_c^r(t)$ 的存在。

2. 反馈控制（能量整形与安全性）：
   利用能量整形方法推导了轨迹跟踪控制器。系统的总能量被定义为包含热能以及移动边界的动能/势能。控制律设计为：
   $$q_c(t) = q_c^r(t) - c(E(t) - E_r(t))$$
   其中 $E(t)$ 和 $E_r(t)$ 分别为实际能量和参考能量，$c$ 为正增益。
   为了确保安全性（具体而言，即热通量保持非负且温度保持在熔点以上），作者利用了最大值原理。他们对初始能量和参考轨迹施加了特定假设，以确保闭环系统满足物理约束。这种方法有效地充当了安全过滤器，确保即使参考轨迹要求冷却（负热通量），实际控制输入仍保持非负，从而防止违反相变条件。

3. 稳定性分析：
   利用结合PDE 反步变换的Lyapunov 分析，建立了跟踪误差的全局指数稳定性。跟踪误差系统被转换为具有理想稳定性特性的目标系统。通过构建一个考虑误差状态能量和移动边界动力学的 Lyapunov 函数，作者证明了跟踪误差以指数形式收敛至零。该分析仔细处理了时变域以及 PDE 与二阶 ODE 之间的耦合。

主要贡献

• 二阶动力学： 本文将 Stefan 问题的控制扩展至包含二阶移动边界动力学，从而更好地捕捉相变过程中的热惯性。
• 严格的安全性证明： 与以往通常假设安全性或静态参考的研究不同，本文提供了严格的证明，表明所提出的控制律能保证针对动态参考轨迹满足物理约束（不冻结和域边界）。
• 级数解的收敛性： 作者证明了在二阶边界动力学下逆 Stefan 问题级数展开解的收敛性，这是生成前馈控制的必要步骤。
• 全局指数稳定性： 本文利用能量整形、反步法和 Lyapunov 分析的新颖组合，建立了闭环跟踪误差的全局指数稳定性，无需依赖线性化或关于误差动力学的简化假设。

结果
使用锌熔化的参数进行了数值模拟。参考轨迹被设计为一个随时间衰减的非平凡正弦函数，确保边界单调移动但速度变化。

• 跟踪性能： 闭环系统成功跟踪了参考边界位置和温度分布，跟踪误差在模拟时间内（约 60 分钟）收敛至零。
• 安全性验证： 模拟证实，虽然参考热通量偶尔变为负值（表明理论上需要冷却），但实际控制输入保持非负，满足了安全约束。在整个模拟过程中，温度分布严格保持在熔点以上。
• 瞬态行为： 结果显示，控制器最初施加了较大的热通量，以补偿系统初始能量低于参考值的情况，导致在进入精确跟踪之前出现短暂的超调。

意义
作者声称，这项工作首次为具有轨迹跟踪控制的 Stefan 问题提供了安全性和跟踪稳定性的严格证明。通过解决几何非线性以及 PDE 与二阶 ODE 之间的耦合，本文提供了一种系统的控制设计，在实现精确跟踪的同时确保了解的物理有效性（安全性）。该方法架起了 PDE 理论控制设计与相变应用中满足硬性物理约束的实际需求之间的桥梁。所提出方法的代码已公开，以促进进一步的研究和验证。