---
title: "细颈越来越小，计算什么时候该停"
description: "用体积、能量、网格收敛和物理尺度判断界面计算，理解数值停止与真实断裂之间的距离。"
canonical_url: "https://sandboxphysics.com/zh/stories/dynamic-interfaces/resolution-and-singularities/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/dynamic-interfaces/resolution-and-singularities/"
last_modified: "2026-10-28"
---

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

L3 · 第 7 章

# 细颈越来越小，计算什么时候该停

Sandbox Physics · 2026-10-28

细颈缩小以后，决定下一步运动的曲率也集中在更窄的地方。屏幕还能画出一条光滑曲线，网格却可能早已看不清它。可信的计算需要说明最后能分辨到哪里，并保留到那里为止的结果。

## 开始前需要什么

薄膜的守恒与增长率，以及液丝半径方程。需要理解相对误差；下面用两次网格加密介绍收敛阶，不要求预先学过数值分析。

- 把离散误差、参数误差和模型误差分开检验。
- 设计一个能重复的网格比较，并解释停止时间的含义。

## 先检查不该变化的量

周期薄膜没有流入、流出或蒸发，平均厚度必须保持不变。周期液丝守恒的是截面积积分，因此要平方半径再积分。两种几何看起来相似，直接相加半径却会给出错误的质量检验。

把当前体积与初始体积相减，再除以初始体积，得到相对体积误差。误差很小只能说明这项守恒做对了；一种把液体搬错方向的方法也可能严格守恒，还要看动力学。

$$
V_{\rm film}=\int h\,\mathrm dx,\qquad V_{\rm thread}=\pi\int r^2\,\mathrm dz,\qquad e_V=\frac{V(t)-V(0)}{V(0)}
$$

薄膜体积按单位横向宽度计；液丝体积来自轴对称截面。

## 没有外加驱动时，能量给出第二条线索

无驱动薄膜沿广义压强梯度流动。把能量对高度的变分记为 $\mu_{\rm chem}=-\gamma h_{xx}+f\prime(h)$。这里的下标区别于温差流一章的黏度符号。把高度方程乘以 $\mu_{\rm chem}$ 再分部积分，周期边界项消失，余下的是一个负的平方积分。能量下降时，局部凸起仍可能增长：降低的润湿能可以支付增加的曲率能。

Marangoni 实验持续维持温差，外界可以向液体输入机械功；只检查毛细能是否单调下降会误判。应把表面牵引做功与黏性耗散放在一起比较。液丝还含有动能，不能只看表面积。

$$
\frac{\mathrm dF}{\mathrm dt}=\int\mu_{\rm chem} h_t\,\mathrm dx=-\int M(h)[(\mu_{\rm chem})_x]^2\,\mathrm dx\le0
$$

正迁移率、周期或无通量边界、无外驱动共同保证这个耗散式。

来源：U. Thiele · Eggers & Dupont · Bowen

## 完整例题：结果接近了多少

沿用第三章最大增长率为四分之一的无量纲例子。假设同一物理初态、同一观察时刻，在 $N=64,128,256$ 个格点得到增长率 $0.240,0.2475,0.249375$。相邻差为 $0.0075$ 和 $0.001875$，相差四倍。这组教学数字符合二阶误差随网格间距平方缩小。

若已确认时间误差更小，Richardson 外推给出 $0.249375+(0.249375-0.2475)/3=0.250$。它与第三章独立线性式的最大增长率相符。只看两张外观差不多的图，会漏掉这个可以量化的差别。

真正比较时要保留同一连续初态。若每次换一套随机网格噪声，细网格上多出新的短波，差异就混入初态变化。还要在同一物理时间比较；用各自“断裂前一帧”比较，观察时间已经不同。

$$
p_{\rm obs}=\log_2\left|\frac{Q_{\Delta x}-Q_{\Delta x/2}}{Q_{\Delta x/2}-Q_{\Delta x/4}}\right|
$$

只有处于渐近收敛区、且分母没有被舍入误差淹没时，观测阶才有解释力。

## 数值极限和物理极限要分别记录

网格只能表示比格距更宽的结构。细颈、前驱薄层附近的弯曲区和热点宽度都需要足够的格点；半径小于格距的某个阈值就停止，是分辨率规则。坡度过大、原先忽略的惯性变强、厚度接近分子尺度时，约化方程也可能失效。继续增加格点也不能修复缺少的物理。

停止时保存最后接受的有限、正厚度轮廓，以及触发停止的原因。不要把负高度截成很小的正数继续演化；那会改变质量、曲率和后续断裂位置。液丝模型用单值半径描述一个连通表面，跨过断裂要另行定义拓扑改变后的问题。

来源：Eggers & Dupont · U. Thiele

## 正则化参数本身也是实验条件

前驱膜厚度变了，既可能改变短程力，又可能改变表观接触角。应先说明保持哪些宏观量相同，再比较液滴移动和薄膜变细。把它只当成“让计算不崩”的旋钮，会把物理敏感性藏起来。

一个可重复记录至少包含模型版本、单位、边界、格数、时间步规则、初态种子、全部参数与未修改的高度或半径序列。另存解析参考、拟合区间与误差指标，读者就可以从同一记录重算增长率。

## 先预测，再动手解释

### [观察停止之前的颈缩](https://sandboxphysics.com/zh/labs/filament-necking/)

1. 选择“细颈越变越细”，记录细颈半径的下降，再查看停止提示。指出它描述的是几何极限还是分辨率极限。
2. 增大黏度，重新计算。比较同一毛细时间的半径与线性增长率，解释黏度改变了哪一项。

## 练习与解释性答案

1. 一个模拟体积误差接近舍入精度，就能说明物理正确吗？

   展开解释

   不能。守恒检验只约束总量；还要比较平衡、耗散、独立线性色散与加密后的运动。守恒但反向输运仍然可能给出错误解。

2. 加密网格后停止时间推迟，能否直接宣布断裂时间收敛？

   展开解释

   不能。先比较共同可分辨时间内的半径和轮廓，再研究明确停止阈值的变化。停止时间记录的是这一网格的分辨能力；拓扑断裂时刻需要额外的物理与计算方法。

[上一章：细水流为什么会长成一串水珠](https://sandboxphysics.com/zh/stories/dynamic-interfaces/liquid-thread/) [下一章：从增长率读到纳米薄膜的涨落谱](https://sandboxphysics.com/zh/stories/dynamic-interfaces/paper-reading/)

## 来源与计算范围

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

1. [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 的非线性延拓分支作为阅读内容。

2. [Bowen, Sur, Bertozzi & Behringer, Nonlinear dynamics of two-dimensional undercompressive shocks, Physica D 209, 36–48 (2005)](https://www.math.ucla.edu/~bertozzi/papers/BowenPhysicaD05.pdf)

   第二节式 (2.1) 给出热应力、重力与毛细通量。浏览器模型取水平基底、外加温度场，保留热应力与毛细项。论文中的斜面前沿与实验照片不在复现范围。

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 的水射流条件。课程计算周期细丝基准；喷嘴实验和卫星滴留作论文阅读。