跳到主要内容
Sandbox Physics

M099 · 流体力学

表面张力与毛细实验室

液体在细管里升多高,能分清表面张力和润湿吗? 分不清。细管里液面升多高,只取决于表面张力乘以接触角的余弦,也就是润湿张力。张力小但很润湿的液体,和张力大但不太润湿的液体,可能升得一样高,光看高度拆不开这两样。

比较五根粗细不同的毛细管,看每根管里真实的球面弯月面,并在施加的前进角和后退角之间切换。采集体积平均高度,推出润湿张力。

交互模型表面张力与毛细实验室
体积平均高度就绪\text{就绪}
弯月面中心高度就绪\text{就绪}
接触线高度就绪\text{就绪}
施加的分支接触角就绪\text{就绪}
液体减气体压力跃变就绪\text{就绪}
管径 Bond 数就绪\text{就绪}
有符号竖直张力就绪\text{就绪}
力平衡残差就绪\text{就绪}
高度推断润湿张力就绪\text{就绪}
独立已知接触角时的张力就绪\text{就绪}
采集测点数就绪\text{就绪}
原始坐标拟合残差就绪\text{就绪}
推断状态就绪\text{就绪}
模型适用域就绪\text{就绪}

物理教程

弯月面也是力平衡

背景两台测量装置彼此独立,采集数据的方式相同,但物理定律、几何形状和近似成立的判据各管各的。

为什么值得理解解读估计值之前,先把它和实际几何、参考剖面、原始传感器读数对应起来。

先抓重点

核心问题
力平衡里该用哪个液面高度?
一句话直觉
先确认近似成立,再用彼此独立的原始读数推断参数。

核心数学模型

平均液柱的力平衡

ρgπr2h‾=2πrγcos⁡θ\rho g\pi r^2\overline h=2\pi r\gamma\cos\theta

高度是液柱的体积平均高度,从液位不变的液池算起,低于液池时为负。

球面弯月面

h‾=2γcos⁡θρgrΔp=pl−pg=−2γcos⁡θr\begin{aligned}\overline h&=\frac{2\gamma\cos\theta}{\rho gr}\\\Delta p&=p_l-p_g=-\frac{2\gamma\cos\theta}{r}\end{aligned}

压力跃变取小 Bond 数下球面弯月面的参考值。液柱上升和下降时,它的符号相反。

重力与毛细力比较

Bo=ρgr2/γ\mathrm{Bo}=\rho gr^2/\gamma

把弯月面当成球面,就忽略了重力在弯月面上各处的影响差别。这个数不超过 0.1 的管才采集;这条界线不代表通用的误差容限。

高度数据能定出什么

h‾i=b+S/riγcos⁡θ^=ρgS/2\begin{aligned}\overline h_i&=b+S/r_i\\\widehat{\gamma\cos\theta}&=\rho gS/2\end{aligned}

拟合时让高度截距自由变化,推出润湿张力。再除以另外已知的接触角余弦,才得到张力;中性润湿附近不行。

施加的接触角

θ=clip⁡ ⁣(θ0+bΔθ/2,0∘,180∘)\theta=\operatorname{clip}\!\left(\theta_0+b\Delta\theta/2,0^\circ,180^\circ\right)

接触角取后退、中心或前进三种值之一。这是人为设定的静态比较,不模拟钉扎或润湿动力学。

常见难点

参考曲线不是估计值

容易误解不记录数据,也能读出估计值。

正确理解清空记录,估计值就没了。传感器扰动可能拟合出不合物理的结果;估计值不做截断,也没有独立校准。

动手实验

  1. 01

    比较管径

    选择润湿上升,查看五个弯月面和测量位置视角。

    你会看到: 管越细,液面升得越高。中心处的高度和体积平均高度不一样。
  2. 02

    反转润湿

    比较中性润湿与不润湿下降。

    你会看到: 曲率、带符号的高度和压力跃变一起改变。
  3. 03

    比较前进角与后退角

    选择角度滞后,切换前进与后退。

    你会看到: 液体和管径都相同,施加的接触角不同,静止高度也不同。
  4. 04

    测量整排管子

    采集合格管,增加高度扰动并导出 CSV。

    你会看到: 高度对半径倒数做直线拟合,可以估出润湿张力。中性润湿时,单靠高度定不出张力。

相关实验与文章

动态界面 · 阅读路径

曲率、薄膜、温差流与液丝

这条阅读路径里用到本实验的文章,后面跟着它们推荐的其他实验。

  1. 文章表面张力为什么能把水拉成珠
反馈这个页面有错误、看不明白,或者缺了什么?告诉我们
反馈类型

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