Droplets sitting on thin elastic sheets: A study with the boundary element method

技术摘要：沉积在薄弹性片上的液滴：基于边界元法的研究

问题陈述
本文研究了液体液滴与薄弹性片之间的弹塑性毛细相互作用（elasto-capillary interaction）。与刚性基底不同，柔性固体会在毛细力的影响下发生变形，从而导致复杂的平衡形状，即液滴会凹陷进入片层中，形成透镜状结构。虽然之前的研究已经探索了皱纹图案、趋硬性（durotaxis，向较软区域移动）以及应力下的各向异性形状，但仍需要一个统一的数值框架，用以模拟在各种边界条件和张力协议下，沉积在弹性片上的液滴的平衡形状。具体而言，作者旨在量化片层厚度、各向同性张力和单轴拉伸如何影响液滴形态、应力分布以及光学特性。

方法论
作者扩展了其先前开发的边界元法（Boundary Element Method, Bм），以模拟液滴与弹性片的耦合动力学。该方法的核心包括：

1. 动力学方程： 系统被建模为过阻尼动力学，其中液滴表面（$\dot{R}_0$）与弹性片（$\dot{R}_{sh}$）的速度场是耦合的。控制方程为一个线性系统 $G\dot{R} = -\nabla_R F$，其中 $G$ 是摩擦矩阵，$F$ 是总自由能。
2. 自由能公式：
   • 液滴 ($F_0$)： 包括液-气表面张力和克服基底-液体界面所做的功，受杨氏接触角支配。模型中包含了一项线张力项，但在考虑的弹性片中，该项被设定为零，因为其影响可忽略不计。
   • 弹性片 ($F_{sh}$)： 使用 Skalak 模型进行面内剪切和面积扩张（超弹性）建模，并使用 Helfrich 模型进行面外弯曲建模。弹性模量（剪切模量 $\kappa_S$、面积模量 $\kappa_A$ 和弯曲模量 $\kappa_b$）由体材料属性（杨氏模量 $E$ 和泊松比 $\nu$）以及片层厚度 $d$ 推导得出。
   • 约束 ($F_{cons}$)： 通过拉格朗日乘子强制执行，以保持恒定的液滴体积，并确保液滴网格与弹性片网格在界面处的几何连通性。
3. 摩擦矩阵： 矩阵 $G$ 整合了粘性剪切应力（源自 Oseen 张量）、接触线摩擦（基于 Moffatt 的结果）以及受 Navier 滑移条件支配的液-固摩擦。
4. 数值实现： 系统使用由 gmsh 生成的三角形网格进行离散化。作者采用加权插值方案（使用 Surrogates.jl）在液滴和弹性片的不同网格之间映射力和速度。模拟过程使用接触线半径 $R_0$ 以及基于 Cox-Voinov 定律推导出的特征时间尺度进行无量纲化处理。

关键结果
本研究将扩展后的 BEM 应用于三种特定的协议：

1. 夹紧弹性片：

   • 液滴的平衡形状是一个透镜，随着片层厚度的增加，其相对于接触线平面的不对称性日益增强。
   • 厚度依赖性： 随着厚度 ($d$) 的增加，液滴下半部分的高度 ($h_l$) 减小，而上半部分的高度 ($h_u$) 增加，尽管总高度保持相对恒定。
   • 接触角： 接触线处的微观接触角与杨氏角一致，验证了模型的有效性。然而，宏观接触角（由整体透镜形状决定）与杨氏角显著不同，并取决于片层厚度。
   • 解耦现象： 液滴下半部分的形状主要由弹性片的弹性决定，且几乎独立于平衡接触角 ($\theta_{eq}$)，而 $\theta_{eq}$ 则主要影响液滴上半部分。

2. 各向同性拉伸弹性片：

   • 施加各向同性张力会抬升液滴，减少片层的变形，直到在高延伸率（$\epsilon \approx 0.1$）时片层恢复为平面。
   • 应力分布： 作者量化了方位角应力 ($T_{\phi\phi}$) 和径向应力 ($T_{rr}$) 的径向分布。应力随延伸而增加，且接触线位于应力曲线的拐点附近。
   • 光学调控： 透镜状液滴充当液体光学透镜。通过改变施加的张力，可以实现约 15% 的可逆调节，从而调节焦距（通过透镜制造家方程计算）。这一调节范围与现代智能手机中的透镜调节范围相当。

3. 单轴拉伸弹性片：

   • 沿单一方向拉伸片层会导致液滴发生拉长。拉长的程度取决于侧边边界条件（夹紧、受控向内运动或自由）。
   • 形态特征： 自由侧边会导致最显著的拉长以及显著的垂直位移，从而在片层上产生向上褶皱和凹陷。接触线变得非平面化，在侧向方向抬升，而在拉伸方向下沉。
   • 各向异性： 研究量化了接触线足迹的纵横比以及宏观接触角的各向异性。这种拉长现象归因于液滴在拉伸方向和侧向方向感受到了不同的有效片层刚度。

意义与主张
本文声称扩展后的 BEM 提供了一个鲁棒的数值工具，用于模拟弹塑性毛细系统，成功捕捉了流体动力学与弹性变形之间的相互作用。其主要贡献在于：

• 方法论层面： 建立了一个形式化的 BEM，将液滴表面动力学与弹性片变形（包括面内和弯曲能量）耦合在一起。
• 物理洞察： 证明了片层厚度是控制液滴不对称性的关键参数，并且尽管存在宏观变形，微观接触角仍与杨氏定律保持一致。
• 应用层面： 研究明确展示了通过对弹性片施加张力，可以实现液体透镜焦距的动态且可逆的调节。这为控制微型光流控器件中的光学特性提供了一种潜在机制。
• 未来方向： 作者指出，该方法为研究皱纹图案、张力梯度中的趋硬性以及多液滴间的弹塑性毛细相互作用开辟了道路；此外，单轴张力可能会导致各向异性的相互作用势和复杂的动力学过程。