从增长率读到纳米薄膜的涨落谱

先复算 Thiele 2010 的图 1

打开作者提供的已发表稿，找到第二节式 (5) 和图 1。横轴是波数，纵轴是平膜小扰动的指数增长率。先找零波数、过零点和峰值，分别对应质量守恒、稳定性边界和最快模式。稳定支线在非零波数处位于零线下方；零波数仍由质量守恒固定为零。

对失稳支线定义 $K=k/k_c$，再把增长率除以 $M_0\gamma k_c^4$。所有采用该线性模型的参数组都收成下面同一条无量纲曲线。取 $K=0,1/2,1/\sqrt2,1,5/4$，得到 $0,3/16,1/4,0,-225/256$。这些点来自公式复算，原图没有数字化采样。

来源：U. Thiele

再看液丝公式为什么多出一个时间尺度

打开 Eggers–Dupont 1994 论文的存档版本 arXiv:physics/0110081v1。第二节式 (27) 是增长率的二次方程；式 (28) 给出增长支，式 (30)、(31) 检查低黏与高黏极限。这里保留惯性，速度需要单独演化，不能把薄膜的过阻尼式直接套过来。

以初始半径和惯性毛细时间无量纲化，用 $q=kR_0$、$s=\omega t_c$。在 $\mathrm{Oh}=0$ 的细丝模型中，最强模式为 $q=1/\sqrt2$，波长为 $8.89R_0$；原文也给出完整无黏柱体的参考值约 $9.01R_0$。这个小差别能提醒我们：约化方程也有自己的误差。

原文图 1 则是有喷嘴的水射流形状比较，图注写出 $\lambda/R_0=14.57$。它需要入口条件、射流速度与实验照片对齐。本课程的周期液丝验证式 (27)，图 1 留作理解真实实验几何的阅读对象。

来源：Eggers & Dupont

Zhang 2024 的图 4 画的是谱，不是增长率

继续打开 Zhang 的 arXiv:2406.08179v1，并对照已发表版本 JFM 996, A19。第二节先定义强滑移随机润滑方程；第四节说明分子模拟怎样控制壁面作用。图 3 用受压流动的速度分布标定滑移长度，这是进入去润湿比较前的重要一步。

图 4 固定膜厚 $1.2\,\mathrm{nm}$，比较强滑移 $400\,\mathrm{nm}$ 与弱滑移 $0.2\,\mathrm{nm}$。符号来自 40 次独立分子模拟的均方根表面振幅谱，线条是随机连续理论的预测。它们随时间改变，不能当成第三章的单条确定性增长率曲线。

图 4(a) 还比较保留惯性和删去惯性的理论。强滑移改变流动的主导平衡，原来对弱滑移适用的无惯性近似在这些参数下不够。研究问题由“哪条波长长得快”走向“哪一种速度场与噪声模型能给出观测到的谱”。

来源：Yixin Zhang

把下一步研究写成能回答的问题

第一个扩展可以只比较确定性模型：保持润湿能、平均厚度和初始轮廓相同，改变滑移条件，分别推导相应通量。先比较早期谱和能量，再考虑引入热噪声。进入强滑移范围时，还要重新检查惯性与长波标度。

要计算持续热涨落，就需要方程中的随机通量、噪声与耗散的关系，以及多次独立运行。随机初态只在开始时给一次扰动，不能代替整个过程中持续输入的热噪声。跨越破裂还需要适合拓扑改变的方法，分子模拟则需要材料相互作用、体系大小和统计误差的独立核查。

读下一篇论文时，先写下几何、守恒量、壁面条件、观测量与平均次数，再找理论线条用了哪些参数。最后才比较曲线。这样可以把一个吸引人的图变成下一次可执行的计算。

来源：Yixin Zhang · Eggers & Dupont