---
title: "从增长率读到纳米薄膜的涨落谱"
description: "复算 Thiele 的图 1，核读 Eggers–Dupont 的色散，再比较强滑移与热涨落论文中的观测量。"
canonical_url: "https://sandboxphysics.com/zh/stories/dynamic-interfaces/paper-reading/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/dynamic-interfaces/paper-reading/"
last_modified: "2026-10-28"
---

[动态界面](https://sandboxphysics.com/zh/stories/dynamic-interfaces/)

L4 · 第 8 章

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

Sandbox Physics · 2026-10-28

读界面论文，先问图里的纵轴究竟是什么。同样是一条随波数起伏的曲线，它可能是确定性增长率，也可能是多次分子模拟平均得到的振幅谱。把这两者分清，才能判断一个简化模型应该与哪张图比较。

## 开始前需要什么

第三章的增长率提取、第六章的细丝色散与第七章的误差分析。需要把原文与课程并排打开；本章按具体公式和图号阅读。

- 独立复算两种连续模型的解析基准。
- 识别强滑移和热涨落增加了哪些未知量与观测要求。

## 先复算 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$。这些点来自公式复算，原图没有数字化采样。

$$
\widehat\sigma=K^2(1-K^2)
$$

课程图与可复算基准使用这条解析曲线；原文非线性分支和稳定性延拓不包含在这里。

解析增长率 $\widehat\sigma=K^2(1-K^2)$

$K_*=1/\sqrt2$ $K_c=1$

- 失稳薄膜
- 稳定薄膜（另取增长率尺度）

按 Thiele（2010）图 1 的解析式重新计算的示意图。横坐标为波数除以失稳截止波数；图中没有实验数据点。

来源：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 留作理解真实实验几何的阅读对象。

$$
s^2+3\mathrm{Oh}\,q^2s-\frac{q^2(1-q^2)}2=0
$$

这是轴对称一维细丝模型的连续线性基准，稳定边界在单位无量纲波数。

来源：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

## 先预测，再动手解释

### [把一条保存记录变成可比较的增长率](https://sandboxphysics.com/zh/labs/thin-film-instability/)

1. 导出小振幅薄膜记录，按第三章计算每个模式的斜率，保留拟合时间区间。
2. 与归一化解析曲线比较，注明有限周期域只允许离散波数。
3. 回看 Zhang 图 4，列出若要比较那张图，现有记录还缺哪些物理和统计信息。

## 练习与解释性答案

1. 验证 $K=5/4$ 的图 1 基准，并解释其符号。

   展开解释

   代入得到 $(25/16)(1-25/16)=-225/256$。负值表示该短波受曲率稳定，会指数衰减。

2. 一条有固定随机种子的确定性薄膜记录，能否直接代表热涨落谱？

   展开解释

   不能。它只演化一次初始扰动；热涨落在演化过程中持续输入扰动，论文还把多次独立模拟的振幅平方取平均，再开平方。必须匹配随机动力学和统计定义。

3. 本课程已经复算的结果和仍属阅读的结果分别是什么？

   展开解释

   已复算平膜解析色散与周期细丝线性关系，并用模拟检验它们。Thiele 非线性延拓分支、喷嘴实验照片，以及 Zhang 强滑移分子模拟和随机谱属于明确定位的阅读材料。

[上一章：细颈越来越小，计算什么时候该停](https://sandboxphysics.com/zh/stories/dynamic-interfaces/resolution-and-singularities/)

## 来源与计算范围

例题与示意图由课程编写。论文比较复算指定解析关系；阅读指南另行解释实验照片与分子模拟各自回答的问题。

1. [Yixin Zhang, Linear stability theory and molecular simulations of nanofilm dewetting with disjoining pressure, strong liquid–solid slip, and thermal fluctuations, J. Fluid Mech. 996, A19 (2024)](https://arxiv.org/html/2406.08179v1)

   对照已发表版 DOI 10.1017/jfm.2024.701，核读图 3、4 与第二、四、五节，分清滑移标定与系综谱。这部分为阅读指南；没有复现分子数据或强滑移随机曲线。

2. [U. Thiele, Thin film evolution equations from (evaporating) dewetting liquid layers to epitaxial growth, J. Phys.: Condens. Matter 22, 084019 (2010)](https://www.uni-muenster.de/Physik.TP/~thiele/Paper/Thie2010jpm.pdf)

   作者提供的已发表论文：式 (1)–(3)、第二节势函数、式 (5) 与图 1。复算限于解析线性色散；图 2、3 的非线性延拓分支作为阅读内容。

3. [Eggers & Dupont, Drop formation in a one-dimensional approximation of the Navier–Stokes equation, J. Fluid Mech. 262, 205–221 (1994)](https://arxiv.org/pdf/physics/0110081v1)

   arXiv:physics/0110081v1 是 1994 年论文后来存入的版本。核读第二节式 (17)–(18)、曲率式 (15)、色散式 (27)–(32) 及图 1 的水射流条件。课程计算周期细丝基准；喷嘴实验和卫星滴留作论文阅读。