comsol两相流之水平集法

两相流是流体力学中一类包含物质两相状态的流动问题，可分为固液、液气、气固、液液（两种不均混液体）。两相流问题非常复杂，通常情况下，多相流体系总是由若干种连续介质或一种连续介质和若干种不连续介质组成，在物质的两相界面处各项参数呈断崖式变化，相比于单相流，未知参数成倍增加。多相流的数值模拟难点还在于不同相之间的传热传质交换，不同相分布的确定等等。

在较小的尺度上，可以对相边界的形状进行详细建模；例如，气泡与液体之间的气液界面的形状。在软件中，这种模型称为分离多相流模型，通常使用表面追踪法来描述此类模型。

在较大尺度上，如果必须详细描述相边界，则模型方程无法求解。相反，可以使用场（例如体积分数）描述不同的相。分散多相流模型方程中，相间效应（例如表面张力、浮力和跨越相边界的传递）被视为源和汇。

上图显示了分离和分散多相流模型的主要区别。在上述两种示例中，均使用函数 Φ 来描述气相和液相。但是，在分离多相流模型中，不同相之间相互排斥，并存在一个清晰的相边界，在此边界上相场函数 Φ 发生突变。除了追踪相边界的位置以外，相场函数没有任何物理意义。

在分散多相流模型中，函数 Φ 描述了气相（分散相）和液相（连续相）的局部平均体积分数。通过平均体积分数可以在该区域的任一点顺利地找到介于 0 和 1 之间的值，这预示着在其他均质域中是存在少量还是大量气泡。也就是说，在分散多相流模型中，可以在同一时间和空间点上定义气相和液相；而在分离多相流模型中，在给定的时间和空间点上，只能定义气相或液相。

分离多相流模型

comsol提供了三种分离多相流模型，水平集法，相场法，移动网格法。

水平集和相场都是基于场的方法，其中相之间的界面代表水平集或相场函数的等值面。移动网格法与上述两种方法完全不同，它将相界面模拟为分隔两个域的几何表面，每个域对应不同的相。

基于场的问题通常是在固定的网格上解决，而使用移动的网格可以解决移动网格问题。

对于不希望发生拓扑变化的微流体系统，通常首选移动网格法；如果需要拓扑变化，则必须使用相场法：当表面张力的影响较大时，首选相场法如果可以忽略表面张力，首选水平集法。

但当你的问题涉及湍流模型时，这时一般都会涉及拓扑变化，要选择水平集法。

分散多相流模型

气泡流模型
气泡流模型显然适用于液体中的气泡。由于忽略了分散相的动量贡献，因此该模型仅在分散相的密度比连续相小几个数量级时才有效。

混合物模型
混合物模型与气泡流模型相似，但考虑了分散相的动量贡献。它通常用于模拟分散在液相中的气泡或固体颗粒。混合物模型还可以处理任意数量的分散相。混合物模型和气泡流模型均假设分散相与连续相处于平衡状态，即分散相不能相对于连续相加速。因此，混合物模型无法处理分散在气体中的大固体颗粒。
欧拉–欧拉模型
欧拉-欧拉模型是最精确的分散多相流模型，也是用途最多的模型。它可以处理任何类型的分散多相流。它允许分散相加速，并且对不同相的体积分数没有限制。但是，它为每个相定义了一组 Navier-Stokes 方程。

在实践中，欧拉-欧拉模型仅适用于两相流，并且其计算成本（CPU 时间和内存）较高。因此，使用它也相对困难，并且需要良好的初始条件才能在数值解中收敛。

欧拉–拉格朗日模型
当连续流体中悬浮有一些（成千上万但不是十亿）非常小的气泡、液滴或颗粒时，我们也许可以使用欧拉–拉格朗日模型模拟多相流系统。该方法的优点是计算成本相对较低。从数值的角度来看，这些模型通常也“不错”。因此，当连续流体中分散相的颗粒数量相对较少时，优选欧拉–拉格朗日模型。

此外，还有一些方法可以使用欧拉-拉格朗日模型来模拟大量粒子，它们使用的相互作用项和体积分数可以模拟具有数十亿个粒子的系统。这些方法可以在 COMSOL Multiphysics 中实现，但在预定义的物理接口中无法实现。

案例-上升气泡

这是一个油泡在水中向上移动，最终与顶部的油相融，导致拓扑变化。使用的是水平集方法（非常适合求解几何拓扑随时间变化的移动边界问题）。


1.选择物理场二维轴对称，多相流-两相流，水平集

2.包含初始化的瞬态
3.几何模型构建—→构建第一个矩形，构建一个多边形（水油初始界面），构建一个圆（油滴），构建联合体

3.添加材料：water 和 transformer oil

4.设置水平集接口——添加初始界面

5.设置多物理场接口

6.设置水平集初始化参数
7.设置水平集初始值2

8.设置层流接口

9.设置压力点约束

10.多物理场——润湿壁

11.创建网格之前，添加一个变量来用于计算模型域中油的质量。稍后您将使用该变量
来测试质量守恒。

12.网格划分-自由三角形网格

12.求解器设置，步骤2步长

13.显示默认求解器设置，解1计算

14.后处理——测试油质量守恒——计算

15.表格1中绘制表图


不知道为什么跟表格中绘制的数据对不上。。。存疑1.

16.流体1的体积分数

最终呈现结果——流体1的体积分数

17.向旋转到三维的轴对称模型添加速度大小的切面图。

、
再加个等值面
18.流体1的体积分数二维三维绘图设置

本文参考了多相流建模与仿真，侵删。