1. 声表面波微流体控制的核心原理
声表面波（Surface Acoustic Wave, SAW）在微流体控制中的应用，本质上是通过压电效应将电能转化为机械振动，进而驱动流体运动。当我们在压电基底上施加交变电压时，会产生特定频率的机械波沿着表面传播。这种波遇到液滴时，能量会通过"泄漏波"机制耦合到液体内部。
关键点在于声流（Acoustic Streaming）现象的产生机制。声波在液体中传播时，由于黏性耗散会产生一个稳定的二阶体积力，这个力可以用著名的Nyborg公式表示：
F = -ρ₀< (v·∇)v + v(∇·v) >
其中ρ₀是流体密度，v是声速场的一阶振动速度，尖括号表示时间平均。这个体积力就是驱动微流体运动的"引擎"。
2. COMSOL中的多物理场耦合建模
2.1 基础物理场接口配置
在COMSOL中实现这一过程，需要建立四个核心物理场的耦合：
- 静电接口：模拟IDT电极上的电势分布
- 固体力学接口：计算压电基底中的应力应变
- 压力声学接口：模拟液体中的声波传播
- 蠕动流接口：计算稳态的声流场
特别要注意的是，在压力声学接口中必须选择"黏性"流体模型，否则无法正确模拟声能耗散导致的体积力。
2.2 体积力公式的实现技巧
体积力的植入是模型成功的关键。在COMSOL中有两种实现方式：
- 直接公式输入法：
matlab复制// 在体积力节点中输入Nyborg公式的COMSOL表达式
-(spf.Ux*spf.Ux_x + spf.Uy*spf.Uy_x + spf.Uz*spf.Uz_x)*rho
- 声流多物理场耦合：
使用内置的"声流场耦合"多物理场接口，它会自动计算：
- 声辐射力
- Eckart流
- 边界层流
注意：使用直接公式法时需要手动添加所有速度分量梯度项，而耦合接口会自动处理这些非线性项。
3. 模型构建的实战细节
3.1 几何处理技巧
对于液滴建模，建议采用以下步骤：
- 先用参数化曲线定义液滴轮廓
- 使用"旋转"操作生成3D液滴
- 对液滴-基底接触区域进行局部网格加密
matlab复制// 示例：定义半球形液滴轮廓
radius = 1[mm];
contact_angle = 78[deg];
h = radius*(1-cos(contact_angle));
3.2 材料参数设置
压电材料需要特别注意晶体取向。对于常用的128° YX切型LiNbO₃：
- 创建材料坐标系
- 设置欧拉角(0, -38°, 0)
- 在压电材料属性中选择对应坐标系
流体材料建议使用甘油-水混合物的实测参数：
- 密度：1120 kg/m³
- 声速：1920 m/s
- 动力黏度：0.0012 Pa·s
4. 求解器配置优化
4.1 内存管理策略
三维模型极易出现内存不足问题，推荐采用：
-
几何简化：
- 利用对称性减少计算域
- 将IDT简化为二维边界条件
-
网格策略：
matlab复制// 非共形网格配对设置
physics('pair').feature('pair1').set('master', 'bnd1');
physics('pair').feature('pair1').set('slave', 'bnd2');
- 求解器配置：
- 频域研究使用GMRES迭代求解器
- 开启移位拉普拉斯预处理器
- 对固体力学部分使用直接求解器
4.2 分步求解技巧
- 先运行频域研究获取声场
- 在稳态研究中使用"冻结"功能固定声场
- 仅求解流场方程
实测数据：在64GB内存工作站上，3D液滴模型约需2小时完成计算。
5. 常见问题排查指南
5.1 典型报错处理
| 错误类型 | 可能原因 | 解决方案 |
|---|---|---|
| 不收敛 | 网格太粗 | 在声波波长λ/6处加密网格 |
| 内存不足 | 自由度太多 | 使用迭代求解器+多重网格 |
| 无流动 | 未启用黏性 | 检查压力声学中的流体模型 |
5.2 结果验证方法
-
能量守恒检查：
比较输入电功率与声流动能matlab复制electric_power = abs(0.5*real(V*conj(I))); flow_energy = 0.5*rho*spf.U^2; -
特征速度验证：
理论预测的声流速度应满足：
$$ U_{theory} = \frac{αv^2}{ω} $$
其中α是衰减系数，v是声速幅值
6. 进阶应用：颗粒操控
在基础模型上添加颗粒追踪：
- 新增"流体流动颗粒追踪"接口
- 添加声泳力节点：
matlab复制F_acoustic = f_scattering*grad(p^2)/(2*rho*c^2) - 设置双向耦合考虑颗粒对流场影响
实测发现直径5μm的PS颗粒在20MHz声场中会产生明显的聚集效应，这与文献报道的声辐射力理论相符。