COMSOL Multiphysics® 软件提供了四种可用于模拟自由液面的方法:水平集、相场、动网格和稳态自由表面。作为系列的第一节,我们将讨论水平集和相场法,这两种基于场的方法几乎可以描述任何类型的自由液面。在第二节中,我们计划将本文的求解结果与通过动网格接口获得的结果进行比较。
何为水平集和相场法?
水平集和相场法都是基于场的方法,这类方法将自由液面表征为水平集或相场函数的等值面。自由液面对应固定网格框架下的液体和气体之间的相分界面。
下图为管道内两颗液滴的表面,摘自附加产品“微流体模块”的“案例库”所提供的液滴破碎模型。从这张图中可以看出,尽管液滴的表面非常明显,但液滴周围的单元并没有贴合到液滴表面上。
水平集和相场函数都是由纳维-斯托克斯方程计算的速度矢量进行对流传输的。在水平集和相场法中,对应公式为:
需要注意的是,水平集和相场函数都使用了 Φ。二者的不同在于方程右侧的 F。在初始水平集方法中 F = 0,因此得到纯对流传输方程。然而当 F = 0 时,数值解不仅不稳定,而且大部分情况下实用性很小。所以为了保持相界面清晰,我们在水平集方法的 F 中添加了高阶导数项 Φ 。
在相场法中,F 代表设法将系统的自由能最小化的项。此项也引入了高阶导数 Φ。实际上,相场方程中的源项中包含了四阶项。这意味着,出于实用性考虑,方程经常被分解为两个方程,与此同时,辅助因变量被定义为 Φ 的二阶导数函数形式。COMSOL Multiphysics 中也采取了这种做法。
两种方法均将自由液面的表面张力引入到纳维-斯托克斯方程的源项中。水平集方法利用表征自由边界的水平集等值面的曲率来描述表面张力。相场法根据化学势计算出表面张力对纳维-斯托克斯方程的源项贡献,因为正是化学势导致了界面附近产生的表面张力和相场函数梯度的产生。
在给定的水平集或相场函数值范围内,即自由表面的数值范围内,流体特性从液体平滑地过渡到气体。
穿过自由表面,水平集函数 Φ 在 0 和 1 之间变化,在两种流体内部分别为常数 0 或 1。比如在液相中为 0,气相中为 1。至于自由液体面处,也就是液体和气体之间的分界面上,对应的水平集函数值为 Φ = 0.5。根