多相流通常包括气-液、液-液、液-固、气-固、气-液-液、气-液-固或气-液-液-固混合物的流动。本系列博客主要讨论气-液和液-液混合物，并简要讨论固-气和固-液混合物。此外，我们还将介绍 COMSOL 软件中 CFD 模块和微流体模块中的模型和模拟策略。
不同尺度的多相流建模
使用数学建模可以对不同尺度的多相流进行研究。最小尺度约几分之一微米，而最大尺度可达几米或几十米。由于尺度甚至可以相差大约八个数量级，最大尺度可能比最小尺度大一亿倍,在整个尺度范围内使用相同的力学模型，在数值上无法解析最小尺度到最大尺度的多相流。因此，多相流的建模通常分为不同的尺度。
在较小的尺度上，可以对相边界的形状进行详细建模；例如，气泡与液体之间的气液界面的形状。在软件中，这种模型称为分离多相流模型，通常使用表面追踪法来描述此类模型。
在较大尺度上，如果必须详细描述相边界，则模型方程无法求解。相反，可以使用场（例如体积分数）描述不同的相。分散多相流模型方程中，相间效应（例如表面张力、浮力和跨越相边界的传递）被视为源和汇。
分离多相流模型详细描述了相边界，分散多相流模型则只考虑分散在连续相中的一个相的体积分数。
上图显示了分离和分散多相流模型的主要区别。在上述两种示例中，均使用函数 Φ 来描述气相和液相。但是，在分离多相流模型中，不同相之间相互排斥，并存在一个清晰的相边界，在此边界上相场函数 Φ 发生突变。除了追踪相边界的位置以外，相场函数没有任何物理意义。
在分散多相流模型中，函数 Φ 描述了气相（分散相）和液相（连续相）的局部平均体积分数。通过平均体积分数可以在该区域的任一点顺利地找到介于 0 和 1 之间的值，这预示着在其他均质域中是存在少量还是大量气泡。也就是说，在分散多相流模型中，可以在同一时间和空间点上定义气相和液相；而在分离多相流模型中，在给定的时间和空间点上，只能定义气相或液相。
分离多相流模型
对于分离多相流的模拟，COMSOL Multiphysics® 软件提供了3种不同的界面追踪方法：
水平集法
相场法
移动网格法
水平集和相场都是基于场的方法，其中相之间的界面代表水平集或相场函数的等值面。移动网格法与上述两种方法完全不同，它将相界面模拟为分隔两个域的几何表面，每个域对应不同的相。
基于场的问题通常是在固定的网格上解决，而使用移动的网格可以解决移动网格问题。
下面的动画为一个T型微通道中生产乳液的模拟结果，该模型使用了相场法进行求解。在动画中，我们可以看到相边界与网格的平面和边缘不一致，相边界由相场函数的等值面表示。
在相场法和水平集法中，有限元网格不必与两个相的边界一致。
相反，下图显示了带有移动网格的上升气泡的验证模型。网格与相边界的形状保持一致，并且网格边缘与相边界重合。但是，移动网格模型也有缺点，即气泡的变形使两个次级气泡从母气泡分离。此时，必须将原始相边界划分为几个边界。该方法太复杂，并且尚未在 COMSOL® 软件中实现。因此，COMSOL® 软件中的移动网格法无法处理拓扑变化。而相场法不存在这个缺点，可以处理相边界形状的任何变化。
上升气泡的验证。当两个次级气泡脱离母气泡时，发生拓扑变化。
什么时候使用相场法和移动网格法？
对于给定的网格，移动网格法具有更高的精度。基于这一优势，我们可以直接在相边界上施加力和通量。基于相场的方法需要围绕相边界表面建立密集网格，以解析该表面的等值面。由于很难定义一个精确贴合等值面的自适应网格，因此通常必须在等值面周围建立大量密集网格。在具有相同精度的情况下，与移动网格相比，这样做会降低基于场的方法的表现。那么，什么时候使用这些不同的方法呢？
对于不希望发生拓扑变化的微流体系统，通常首选移动网格法；
如果需要拓扑变化，则必须使用相场法：当表面张力的影响较大时，首选相场法
如果可以忽略表面张力，首选水平集法
分离多相流模型和湍流模型
在湍流模型中，由于仅解析平均速度和压力，流体的细节会丢失。从这一点来看，表面张力效应在流体的宏观描述中也变得不那么重要。由于湍流表面的流动也比较剧烈，因此几乎不可能避免拓扑变化。所以对于湍流模型和分离多相流模型的组合，最好使用水平集法。水平集法和相场法都可以与 COMSOL Multiphysics 中的所有湍流模型结合使用，如下图和动画所示。
在COMSOL Multiphysics中，所有湍流模型都可以与相场法和水平集法相结合来模拟两相流。
将水平集法与 k-e 湍流模型相结合模拟反应堆中水和空气的两相流。
分散多相流模型
万一相边界过于复杂而无法解析，则必须使用分散多相流模型。
CFD 模块提供了 4 种不同的模型（原理上）：
气泡流模型适合高密度相中包含较小体积分数低密度相
混合模型适合连续相中包含较小体积分数的分散相（或几个分散相），其密度与一个或多个分散相相近
欧拉–欧拉模型适用于任何类型的多相流
可以处理任何类型的多相流，且气体中有密集颗粒，例如流化床
欧拉–拉格朗日模型适合包含相对较少（成千上万，而不是数十亿）的气泡、液滴或悬浮颗粒流体
适合气泡、颗粒、液滴或使用方程模拟的颗粒，该方程假定流体中每个颗粒的力平衡
什么时候使用不同的分散多相流模型？
气泡流模型
气泡流模型显然适用于液体中的气泡。由于忽略了分散相的动量贡献，因此该模型仅在分散相的密度比连续相小几个数量级时才有效。
混合物模型
混合物模型与气泡流模型相似，但考虑了分散相的动量贡献。它通常用于模拟分散在液相中的气泡或固体颗粒。混合物模型还可以处理任意数量的分散相。混合物模型和气泡流模型均假设分散相与连续相处于平衡状态，即分散相不能相对于连续相加速。因此，混合物模型无法处理分散在气体中的大固体颗粒。
当多相流混合物被迫通过孔口时，用混合物模型模拟了5种不同大小的气泡。流动中的剪切力导致较大的气泡破裂成较小的气泡。
欧拉–欧拉模型
欧拉-欧拉模型是最精确的分散多相流模型，也是用途最多的模型。它可以处理任何类型的分散多相流。它允许分散相加速，并且对不同相的体积分数没有限制。但是，它为每个相定义了一组 Navier-Stokes 方程。
在实践中，欧拉-欧拉模型仅适用于两相流，并且其计算成本（CPU 时间和内存）较高。因此，使用它也相对困难，并且需要良好的初始条件才能在数值解中收敛。
使用欧拉-欧拉多相流模型模拟流化床中固体颗粒的体积分数。
欧拉–拉格朗日模型
当连续流体中悬浮有一些（成千上万但不是十亿）非常小的气泡、液滴或颗粒时，我们也许可以使用欧拉–拉格朗日模型模拟多相流系统。该方法的优点是计算成本相对较低。从数值的角度来看，这些模型通常也“不错”。因此，当连续流体中分散相的颗粒数量相对较少时，优选欧拉–拉格朗日模型。
此外，还有一些方法可以使用欧拉-拉格朗日模型来模拟大量粒子，它们使用的相互作用项和体积分数可以模拟具有数十亿个粒子的系统。这些方法可以在 COMSOL Multiphysics 中实现，但在预定义的物理接口中无法实现。
附加的 CFD 模块和粒子追踪模块可在 COMSOL Multiphysics 中模拟欧拉-拉格朗日多相流模型。
混合物模型能够处理任何相的组合，并且计算成本较低。在大多数情况下，我们可以使用此模型模拟。对于流化床（具有高密度和高体积分数的大颗粒分散相）之类的系统，只能使用 Euler-Euler 模型模拟。
分散多相流模型和湍流模型
分散多相流模型本质上是近似的，并且也与近似的湍流模型非常吻合。可以在分散相和连续相之间以及在分散相中的气泡、液滴和颗粒之间引入相互作用。这些相互作用的起源可以是用湍流模型模拟的湍流。气泡流、混合物流和欧拉-拉格朗日多相流模型可以与 COMSOL Multiphysics 中的所有湍流模型结合使用。仅针对具有可实现性约束的标准 k-e 湍流模型预定义了 Euler-Euler 多相流模型。
混合物模型可以与 COMSOL Multiphysics 中的任何湍流模型结合使用。
结语
即使可以使用超级计算机求解多相流数值模型方程，这也可能是一项非常艰巨的任务。如果没有计算能力的限制，则表面追踪方法将用于所有类型的混合。实际上，这些模型仅限于微流体以及用于黏性液体自由表面的研究。
分散多相流法允许研究包含数以亿计的气泡、液滴或颗粒的系统。但是，即使是最简单的分散多相流模型，也可以生成非常复杂且要求很高的模型方程。上述几种不同的模型非常适合描述特定的混合物，并能满足工程师和科学家以相对较高的精度和合理的计算成本研究多相流的要求。
我们将继续讨论有关多相流的模拟，敬请关注 COMSOL 博客更新！
