COMSOL液滴与基板碰撞变形，包含2个模型：单液滴碰撞铺展，双液滴碰撞融合铺展。
当咖啡滴落在办公桌的瞬间，我们肉眼可见的液滴变形背后藏着复杂的流体力学机制。今天咱们用COMSOL复现两个典型场景：单液滴撞击铺展和双液滴碰撞融合，看看数值模拟如何捕捉这些微尺度流动的细节。
先说单液滴撞击。初始设置需要定义直径0.5mm的水滴以3m/s速度垂直撞击疏水基板。建模时建议采用轴对称简化，这样计算量能减少80%以上。关键是要处理好移动接触线的处理，这里采用动态接触角模型：
contact_angle = 80[deg]*(1 - exp(-t/1e-5[s]));这个指数衰减公式模拟了接触角从初始值到平衡值的过渡过程，时间常数1e-5秒的设置需要根据实际材料特性调整。计算时记得开启Level Set方法追踪气液界面，相场参数的平滑宽度建议取网格尺寸的1.2倍左右。
网格划分有讲究：撞击中心区域采用边界层网格，基板附近设置5层渐缩网格，最小单元尺寸0.5μm。运行后会发现最大铺展直径出现在撞击后0.8ms左右，此时液膜边缘会出现明显的冠状飞溅现象——这验证了Ohnesorge数大于0.1时惯性主导的特征。
COMSOL液滴与基板碰撞变形，包含2个模型：单液滴碰撞铺展，双液滴碰撞融合铺展。
当两个液滴以相位差撞击基板时，情况变得有趣起来。双液滴模型需要构建三维几何，但可以通过对称面减少计算量。碰撞速度差控制在0.2m/s以内时，融合过程会出现明显的"桥接"现象。这时候要特别注意表面张力项的设置：
sigma = 0.072[N/m]*(1 + 0.01*(phi1^3 + phi2^3));这种非线性表面张力系数能更好反映融合瞬间的界面突变。计算结果显示，当液滴间距小于初始直径的1.5倍时，融合后的液膜会呈现明显的非对称波动，这在实际喷涂工艺中会导致涂层不均匀问题。
在求解器设置上，建议先用瞬态求解器跑前50μs，再切换到相场稳态求解。遇到计算发散时，试试把人工黏度系数调高到3e-4 kg/(m·s)，这招能有效稳定界面振荡。后处理阶段，用表面参数提取器抓取接触线移动速度，和经典的Renardy公式对比误差能控制在5%以内。
模拟中最反直觉的是基板温度的影响。当把基板设为60℃时，看似应该加快铺展，但实际模拟显示最大铺展半径反而减小了8%——这是因为温度升高降低了表面张力但增大了黏性耗散，这种竞争关系需要结合无量纲数分析。
最后给个实用技巧：在参数化扫描时，把撞击速度设为对数分布（比如从1到10 m/s分5个量级），这样能一次性捕捉到铺展、飞溅、破碎多种流态。记得输出视频时打开粒子追踪功能，动态显示涡旋结构的演变过程，这样的可视化结果发文章绝对加分。
