[论文解读] Finite-time singularities in the dynamical evolution of contact lines

Dmitry E. Pelinovsky, Ayrat Giniyatullin|arXiv (Cornell University)|Feb 5, 2013

Fluid Dynamics and Thin Films参考文献 6被引用 4

一句话总结

本文研究了在180°接触角条件下，接触线动力学的线性对流-扩散方程中的有限时间奇点。通过先验能量估计和自相似解，证明了若变速度 $ V(t) $ 不振荡，则解将在有限时间内爆破，并推导出爆破速率为 $ V(t) \sim \sqrt{a_4 / (t_0 - t)} $，与早期数值结果所声称的对数发散相矛盾。

ABSTRACT

We study finite-time singularities in the linear advection-diffusion equation with a variable speed on a semi-infinite line. The variable speed is determined by an additional condition at the boundary, which models the dynamics of a contact line of a hydrodynamic flow at a 180 contact angle. Using apriori energy estimates, we derive conditions on variable speed that guarantee that a sufficiently smooth solution of the linear advection--diffusion equation blows up in a finite time. Using the class of self-similar solutions to the linear advection-diffusion equation, we find the blow-up rate of singularity formation. This blow-up rate does not agree with previous numerical simulations of the model problem.

研究动机与目标

• 分析变速度线性对流-扩散方程解在何种条件下会出现有限时间奇点。
• 评估先前数值结果的有效性，即在爆破时间附近 $ V(t) \sim C_1 \log(t_0 - t) + C_2 $ 的假设。
• 利用能量估计和自相似解，确定 $ V(t) $ 和 $ \beta(t) $ 的正确爆破速率。
• 研究第三类边界条件 $ h_{xxx}|_{x=0} = -1/2 $ 在奇点形成和解存在性中的作用。
• 调和简化薄膜流模型中数值模拟与解析预测之间的差异。

提出的方法

• 应用先验能量估计，证明若对所有时间都有 $ V(t) > 0 $，则解将在有限时间内爆破。
• 通过尺度变换 $ x \sim (1-H)^{-1/3} $，$ t \sim (1-H)^{-4/3} $ 推导自相似解，得到关于 $ G(z) $ 的简化常微分方程。
• 使用Heun方法对解 $ G(z) $ 进行数值逼近，以分析其在 $ z_0 $ 附近的行为，其中 $ z_0 $ 是 $ G(z) $ 的最大负实根。
• 从偏微分方程推导点态方程：$ u_{xxxx}|_{x=0} = V(t)\beta(t) $ 和 $ \frac{d\beta}{dt} + u_{xxxxx}|_{x=0} = -\frac{1}{2}V(t) $。
• 在 $ t \uparrow t_0 $ 附近对系统进行渐近分析，假设 $ \beta(t) \to 0 $，$ u_{xxxx}|_{x=0} \to a_4 \neq 0 $，且 $ u_{xxxxx}|_{x=0} \to a_5 < \infty $。
• 将推导出的爆破速率与自相似解速率及数值上声称的对数速率进行比较。

实验结果

研究问题

• RQ1 在何种变速度 $ V(t) $ 条件下，线性对流-扩散方程的解会在有限时间内爆破？
• RQ2 自相似解类是否满足原始第三类边界条件 $ h_{xxx}|_{x=0} = -1/2 $？若不满足，这对物理相关性有何含义？
• RQ3 在奇点附近，$ V(t) $ 的正确渐近爆破速率是什么？其与数值报告的对数发散速率相比如何？
• RQ4 是否存在全局光滑解？若存在，$ V(t) $ 和 $ \beta(t) $ 必须表现出何种行为？
• RQ5 为何推导出的爆破速率与自相似解速率及数值模拟结果均不一致？

主要发现

• 若对所有 $ t $ 都有 $ V(t) > 0 $，或 $ \beta(t) $ 与 $ V(t)\beta^2(t) $ 均保持为负，则解将在有限时间内爆破，这意味着全局光滑解的存在要求 $ V(t) $ 必须振荡。
• 自相似解类给出 $ \beta(t) > 0 $ 对所有 $ t \in [0, t_0) $ 成立，与初始时物理要求 $ \beta(t) < 0 $ 相矛盾，因此其与原始问题无物理相关性。
• 点态渐近分析得到爆破速率 $ V(t) \sim \sqrt{a_4 / (t_0 - t)} $，其中 $ a_4 > 0 $，与数值上声称的对数速率 $ V(t) \sim C_1 \log(t_0 - t) + C_2 $ 不一致。
• 与数值结果的差异表明渐近假设中可能存在 $ a_4 = 0 $ 或 $ a_5 \to \infty $ 的情况，提示数值分辨率或建模中可能存在缺陷。
• 推导出的爆破速率 $ \sim \sqrt{1/(t_0 - t)} $ 与自相似解速率及数值预测均不一致，凸显该边值问题动力学的复杂性。
• 结果表明，需开展更精确的数值模拟与解析研究，以澄清此接触线模型中奇点形成的本质。