[论文解读] Heteroclinic orbits, mobility parameters and stability for thin film type equations

Richard S. Laugesen, M. C. Pugh|ArXiv.org|Mar 29, 2000

Fluid Dynamics and Thin Films参考文献 11被引用 3

一句话总结

本文通过数值方法研究了具有非线性四阶与二阶扩散项的薄膜方程中异宿轨、稳定性以及迁移率参数的影响。结合能量分析与数值模拟，结果表明：不稳定周期稳态的扰动会演化为常数或液滴状稳态，长期行为由能量最小化决定；尽管时间尺度发生变化，异宿连接在迁移率参数改变时仍能保持存在。

ABSTRACT

We study numerically the phase space of the evolution equation h_t = -(h^n h_{xxx})_x - B (h^m h_x)_x . Here h(x,t) is nonnegative, n>0 and m is real, and the Bond number B is positive. We pursue three goals: to investigate the nonlinear stability of the positive periodic and constant steady states; to locate heteroclinic connecting orbits between these steady states and the compactly supported 'droplet' steady states; and to determine how these orbits change when the 'mobility' exponents n and m are changed. For example, we change the mobility coefficients in such a way that the steady states are unchanged and find evidence that heteroclinic orbits between steady states are perturbed but not broken. We also find that when there appear to be touch-down singularities, the exponent n affects whether they occur in finite or infinite time. It also can affect whether there is one touch-down or two touch-downs per period.

研究动机与目标

• 研究具有幂律迁移率的薄膜方程中常数与周期稳态的非线性稳定性。
• 定位并分析连接周期稳态、常数稳态与紧支集液滴稳态的异宿轨。
• 确定迁移率指数 $n$ 与 $m$ 的变化如何影响解的动力学行为，包括有限时间奇点与稳定性。
• 评估在迁移率参数变化时，异宿连接是否仍能保持存在，同时保持稳态与能量水平不变。
• 探索控制有限时间接触塌陷奇点出现的临界迁移率指数。

提出的方法

• 在周期域上对薄膜方程 $h_t = -(h^n h_{xxx})_x - \mathcal{B}(h^m h_x)_x$ 进行数值延拓与模拟，采用自适应时间步长。
• 利用方程关于时变加权 $H^{-1}$ 内积的梯度流结构进行基于能量的分析。
• 通过分岔图与弱非线性分析，根据 $m-n$ 对稳态的稳定性进行分类。
• 系统性地改变迁移率指数 $n$ 与 $m$，同时通过固定 $q = m - n + 1$ 保持能量泛函与稳态不变。
• 采用时间推进格式，在奇点附近自适应提高分辨率，以检测有限时间接触塌陷事件。
• 通过比较扰动下的解轨迹，判断其是否收敛至常数、周期或液滴稳态。

实验结果

研究问题

• RQ1 周期稳态的线性稳定性是否能准确预测其在数值模拟中的非线性稳定性或不稳定性？
• RQ2 能否在数值上观测到连接不同稳态（如周期到常数或液滴）的异宿轨？其行为如何依赖于 $n$ 与 $m$？
• RQ3 迁移率指数 $n$ 与 $m$ 的变化如何影响有限时间接触塌陷奇点的发生时机与频率？
• RQ4 当迁移率参数改变时，即使能量景观保持不变，异宿连接是否仍能保持存在？
• RQ5 能量最小化在决定扰动解的长期行为中起何种作用？

主要发现

• 常数稳态在扰动下是渐近稳定的，解会随时间恢复至该状态。
• 当正周期稳态在某一方向受扰动时，解会演化至常数稳态；但在相反方向扰动时，会导致有限时间接触塌陷。
• 数值上观测到了连接周期、常数与液滴稳态的异宿轨，支持了基于能量的山口情景预测。
• 当迁移率参数 $n$ 与 $m$ 变化但保持 $q = m - n + 1$ 不变时，异宿连接的形状在定性上保持相似，尽管时间尺度发生显著变化。
• 指数 $n$ 决定接触塌陷发生在有限时间还是无限时间，且可影响每周期内接触塌陷事件的数量（一个或两个）。
• 未发现虚假奇点的证据，自适应时间步长确认：仅在真实奇点附近使用小时间步长，而非由数值不稳定性引起。