跳到主要内容
Sandbox Physics
动态界面

L3 · 第 7 章

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

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

先检查不该变化的量

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

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

Vfilm=∫h dx,Vthread=π∫r2 dz,eV=V(t)−V(0)V(0)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)}
薄膜体积按单位横向宽度计;液丝体积来自轴对称截面。

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

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

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

dFdt=∫μchemht dx=−∫M(h)[(μchem)x]2 dx≤0\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,256N=64,128,256 个格点得到增长率 0.240,0.2475,0.2493750.240,0.2475,0.249375。相邻差为 0.00750.0075 和 0.0018750.001875,相差四倍。这组教学数字符合二阶误差随网格间距平方缩小。

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

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

pobs=log⁡2∣QΔx−QΔx/2QΔx/2−QΔx/4∣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

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

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

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

先预测,再动手解释

观察停止之前的颈缩

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

练习与解释性答案

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

    展开解释

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

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

    展开解释

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

来源与计算范围

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

  1. U. Thiele, Thin film evolution equations from (evaporating) dewetting liquid layers to epitaxial growth, J. Phys.: Condens. Matter 22, 084019 (2010)

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

  2. Bowen, Sur, Bertozzi & Behringer, Nonlinear dynamics of two-dimensional undercompressive shocks, Physica D 209, 36–48 (2005)

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

  3. Eggers & Dupont, Drop formation in a one-dimensional approximation of the Navier–Stokes equation, J. Fluid Mech. 262, 205–221 (1994)

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

反馈这个页面有错误、看不明白,或者缺了什么?告诉我们
反馈类型

提交时会附带本页地址、Lab ID、网站语言和设备类型。请不要填写个人信息。反馈最多保存 12 个月。隐私政策