如何计算搅拌器和旋转机械中的自由液面?
搅拌器或搅拌反应器中的叶轮平稳搅动时，液体表面会产生波动。如果此波动的高度远小于容器中流体的高度，那么可以在一个专门的求解步骤中根据速度场显式地计算出自由液面的形状和高度。COMSOL Multiphysics 的“搅拌器模块”中包含了一个“稳态自由表面”特征，专门用于这类计算量不大的运算。
计算自由表面的形状：各种不同算法
可以有多种方法来计算与气体接触的自由液体的表面形状。举例来说，自由表面可发生程度很大的变形，并引起拓扑变化（即波浪破碎或者液滴脱离表面），那么便可以通过水平集方法和相场方法有效地预测自由表面的形状。这两种方法适用于瞬态仿真和稳态仿真。
采用水平集方法仿真的填充封闭场的过程。仿真过程产生了新的边，图中清晰地显示了所引起的拓扑变化。
现在我们来看另一种情况：自由表面变形程度很大，但没有发生拓扑变化。任意拉格朗日–欧拉（arbitrary Lagrangian-Eulerian，简称 ALE）移动网格方法能够精确地预测此类自由表面的形状，这种方法同样适用于瞬态仿真和稳态仿真，是“搅拌器模块”中用于计算自由表面形状的默认方法。
用 ALE 方法模拟气泡上升。气泡在外半径处即将分裂时停止仿真，因为分裂会生成新的边，引起拓扑变化。
如果计算时选择了合适的网格解析度，那么水平集、相场和 ALE 移动网格这几种方法都可以非常准确地预测自由表面的形状。不过，要获得这种精确度，计算量会很大。而且，这只是预测了自由表面的形状，当这些方法还要用于之后的其他运算时，计算量之大令人难以想象。
在 COMSOL Multiphysics 5.2 版本中，“搅拌器模块”使用了一种计算量相对较少的方法来计算稳态流动中自由流体的表面形状。这一方法适用于自由表面变形程度较小的情形，例如表面只有波纹。这时会用到一个方程，其中使用速度场和压力场作为输入，通过求解表面平均高程来预测自由表面的形状，见方程 (1)。
左图：近似稳态冻结转子中自由表面模型的仿真结果。右图：全瞬态仿真中自由表面模型的结果。
搅拌器 App 展示了这一新特征。请注意，对于涉及旋转机械的稳态流体流动问题，通常都可以利用这一方法计算自由表面形状。下图展示了一个相关应用。
战中使用的鱼雷仿真模型。鱼雷行进时接近水面，使水面产生波动。
接下来，我们将进一步探讨稳态自由表面计算背后的物理原理及稳态自由表面这个特征本身。
稳态自由表面计算背后的原理
要计算自由表面的变形情况，根据 Navier-Stokes 方程求出流场的解，得到边界上的压力分布
其中
利用稳态自由表面特征
稳态自由表面特征位于旋转机械，流体流动接口，可轻松地预测由于体积力和流体流动力共同作用而产生的自由表面变形。
如下方的截图所示，使用“搅拌器模块”时，模型树中加入了“稳态自由表面”特征。选中该特征后，可以进一步启用“稳态自由表面”研究步骤。这个研究步骤的目的是，利用上一个稳态研究步骤求得的速度和压力场，求解上文提到的表面高程方程。
模型树中的“稳态自由表面”特征及其对应的研究步骤。
正如其名所示，“稳态自由表面”特征只适用于平稳的流动。（我们在上文提到，水平集、相场和 ALE 自由表面计算功能都可用于计算流动的瞬态仿真。）扩展“稳态自由表面”特征的一个很自然的做法是对瞬态仿真添加类似功能。
延伸阅读
从“案例下载”中下载相关教程
具有自由表面和底部叶轮的搅拌器
搅拌器
想了解 COMSOL Multiphysics 5.2 版本中还有哪些新增功能吗？赶紧前往 5.2 版本发布亮点页面
本文内容来自 COMSOL 博客，点击“阅读原文”，下载相关教程。
