地下水是灌溉水和饮用水的重要来源，因此保持它的纯净性至关重要。然而，一旦进入污染物，将它们排出去就不是那么容易了。要建立从地下水中去除污染物的系统，需要了解非混相流体在土壤中迁移的复杂行为。借助 COMSOL® 软件中的相传递功能，工程师可以对这种行为进行深入探究。

预测地下流体的迁移

流体，特别是非混相流体在地下的迁移行为对许多行业都很重要。以前，有研究人员对这种传递方式进行过研究，用于改进石油的开采和回收方法。如今，研究人员更加关注非混相溶流体对环境产生的负面影响。例如，废弃物管理和环境保护等行业的工程师正在设计能够更好地预防、检测和清理土壤中污染物的系统。如果不加以控制，污染物会进入地下水，并影响农作物的质量。

地下水插图

然而，由于这个过程所涉及的复杂性和非水相液体的性质，理解流体流动是一件困难的事情。油和其他非水相液体污染物的密度通常与水不同，即使它们混合在一起，也会分开流动（多相流）。非水相液体还可能被困在如淤泥和黏土等渗透性低的土壤区域，使得这些污染物更难被检测和清除。此外，非混相流体如何流过土壤（以及任何低渗透区）取决于土壤的饱和度。      为了深入理解流体在土壤中的流动，工程师可以借助数值仿真技术进行模拟与分析。本文，我们来看一个模拟多相流通过含有低渗透性晶体的多孔介质的示例模型。

模拟多孔介质中的多相流

我们的示例模型包含一个圆形域，代表中间含低渗透性晶体的多孔介质。利用轴对称性，我们只需要建立一个二维横截面模型，从而降低仿真的计算成本。首先，我们用一种比非水相流体轻的流体（比方说水）充满整个域。模拟开始后，顶部的小孔允许较重的非水相液体以恒定且均匀的质量通量进入该域，持续 100 分钟。

在这里我们不会详细介绍所有的建模细节，您可以在 COMSOL 案例库中的低渗透性晶体上方的两相流教程模型中找到流体和土壤的属性，也可以下载相关的 MPH 文件。

我们可以使用 COMSOL Multiphysics® 软件中的 CFD 模块和地下水流模块中的相传递功能轻松模拟两相流。使用多孔介质相传递接口，可以分析非混相流体如何在土壤中流动，从而确定饱和度，并考虑相界面发生的微观现象，以及重力的影响。我们还可以通过多孔介质多相流多物理场耦合节点将这个接口的功能与达西定律接口相结合，这样就可以更轻松地模拟多孔介质中多相流体的流动和传递。

虽然这个示例中没有使用，但我们也可以使用相传递接口模拟自由流动的多相流动，无需追踪相界面，从而使计算要求低于其他方法。

仿真结果显示了低渗透性晶体周围的两相流。

刚开始，预期重相会在低渗透性区域周围流动，但不会渗透进去，因为晶体的入口毛细管压力高于周围土壤的压力，这意味着只有当其饱和度达到临界值时，密度更大的非润湿流体才能进入晶体。要模拟这个过程，我们可以使用多孔介质不连续性边界条件，它允许土壤-晶体界面处不连续，并帮助确定临界饱和度。

可视化结果

运行模拟后，我们可以看到重相的流动与预期的结果一致。流体不会立即渗入晶体，而是仅在其周围的土壤达到临界饱和度时才进入，这发生在大约第 12 分钟。一个小时后，重相到达晶体的底部。

正如文  中的示例所展示的，工程师可以使用 COMSOL Multiphysics 中的相传递功能分析单相和多相流体在多孔介质中的传输。通过更深入地理解所涉及的物理过程，他们可以设计更加高效地去除地下水污染，以及开采石油等系统。