在 COMSOL 中简化地下水流建模

COMSOL Multiphysics® 软件中的地下水流模块包含了一些有用的功能,能够更高效地设置复杂的模拟任务。例如,在对井进行建模时,使用井功能进行设置,网格划分明显变得简单,也更加直观。在这篇文章中,我们将介绍井功能,并讨论如何使用这项功能以及它如何增强建模过程。

在 COMSOL Multiphysics® 中对井进行建模

对地下水流问题进行建模,通常需要处理暴露在相对较小的源或汇中的大型建模域。
以前,COMSOL Multiphysics 中将井作为大型分层地下域中的小型三维圆柱体引入,例如关于 地热回灌中传热和多孔介质流的耦合仿真 的文章中所介绍的。这种方法需要使用合适的边界条件,并且会涉及对小对象进行网格划分。

使用 COMSOL Multiphysics 井边界条件,也就是用边替换圆柱体对于网格划分算法更好,并且需要的网格划分也更少。与精确求解细节相比,这个功能可以提供了准确的解。接下来,让我们来更详细地讨论井边界条件。

井边界条件的设置

首先要熟悉新的边界条件及其设置。井边界条件可用作二维中的点特征和三维中的边特征,并可与达西定律、理查兹方程和两相达西定律接口一起使用。使用这个边界条件,可以选择井是注入井还是生产井,并指定压力或质量流量。下图显示了一些不同的可用选项。
在 COMSOL 中简化地下水流建模的图1 注入井建模的设置达西定律、理查兹方程接口(左)和两相达西定律接口(右),其中还必须指定饱和度。

比较模拟井的两种方法

现在,让我们看看井边界条件与其他用于模拟井的选项相比如何。为了便于说明,我们使用了一个基本模型,如下图所示。

在 COMSOL 中简化地下水流建模的图2
半径为 20m、高度为 3m 的水库中,半径为 0.5m 的井的几何模型,其周围是一个无限的单元域。

使用无限元是为了使我们可以在离井很远的地方施加压力而不增加建模域。这里显示的几何图形将井解析为一个圆柱形的表面。为了能够应用边界条件,必须将井的圆柱体从储层中切割出来。另外,也可以使用质量通量的边条件,但前提是我们要应用质量通量而不是压力。我们可以使用井边界条件,它适用于压力和质量通量条件。
我们用完全相同的网格设置来比较这两种情况下的网格。在这个案例中,我们划分了 65,674 个域单元,而使用井边界条件,仅划分了 28,728 个域单元。这还不到网格单元数量的一半。

在 COMSOL 中简化地下水流建模的图3
使用相同设置的在完全解析井时和使用井边界条件时的网格比较。

这个优势只有在我们得到一个准确的解时才有用。继续使用这个测试案例,我们在井口施加一个 1 kg/s 的质量流速,M0。这相当于在面积为 A 的边界处的质量通量为  在 COMSOL 中简化地下水流建模的图4 。在长度为 l 的边处的质量通量为  在 COMSOL 中简化地下水流建模的图5 。压力在外部无限元的边界是固定的。
一维绘图显示沿中心线的压力与井外的方法几乎完全一致。

在 COMSOL 中简化地下水流建模的图6
沿截线的压力比较。

与指定边的质量通量相反,井功能考虑了井半径,即使没有明确解析也考虑了。质量通量 边特征不合适计算井中的压力,因为它不考虑径向的膨胀。井边界条件提供了一个变量 dl.well1.p,该变量给出了井压。

传热与多孔介质流的耦合应用:地热回灌

井边界条件可用于前面提到的地热回灌示例。我们在之前的文章中介绍的模型的略微修改版本。在这种情况下,地热地下水以 150 l/s 的速度通过生产井生产。在用于产生热量后,水再以相同的速率重新注入,此时水温为 5°C。在外部边界施加 2 mm/m 的水平水力梯度。

在 COMSOL 中简化地下水流建模的图7
地热回灌模型的模型设置(左)和网格(右)。

在 COMSOL Multiphysics 5.3 之前,生产井和注入井被绘制为嵌入地质构造中的圆柱体,并使用达西定律的入口和出口 边界条件以及传热 的温度和流出边界条件在圆柱体表面定义质量和热通量。
现在,井由单边定义,新的井 特征和线热源特征用于定义质量和热通量。这两个功能的设置如下图所示。

在 COMSOL 中简化地下水流建模的图8
井特征(左)和 线热源特征(右)的设置。

在井特征中,通过设置 M0= 150 l/s ρwater。在材料节点中指定水的密度,该节点可通过表达式 mat5.def.rho 访问。对于线热源,我们根据以下公式定义每单位长度的源项 Ql = MlCpΔT。这里 Ml 是每单位长度的质量流量,由井 特征 (dl.well1.Ml)计算; Cp 是水的比热容(mat5.def.Cp); ΔT = Tinj – T 是注入温度与实际温度之间的温差。
使用井特征,与使用圆柱体作为钻孔的模型相比,网格单元(118,000 vs. 126,000)减少了约 8%;模拟运行速度提高了大约 10%(31 min vs. 26 min)。下面的动画显示了 5 年内的温度演变。

在 COMSOL 中简化地下水流建模的图9

五年来地热储层温度的演变。

为了证明井特征给出的结果与圆柱表面上的相应边界条件相同,我们比较了生产温度的结果。可以看到,两个结果具有很好的一致性。

在 COMSOL 中简化地下水流建模的图10
两种建模选项的生产温度比较。

结束语

在这篇文章中,我们看到新的井边界条件可以提高模型性能,并使井建模变得更加容易。我们还了解了相关的仿真背景以及如何设置与传热的耦合。

本文来自: COMSOL 博客

COMSOL地下水流模块传热

在 COMSOL 中简化地下水流建模的评论0条

    暂无评论

    在 COMSOL 中简化地下水流建模的相关视频课程

    在 COMSOL 中简化地下水流建模的相关案例教程

    在 COMSOL Multiphysics 中可以使用 AC/DC 模块中的非线性磁性材料数据库中的非线性磁饱和曲线进行频域仿真。您也可以使用有效非线性磁曲线计算器仿真 App 将关联的 B-H 或 H-B 曲线(以前仅支持稳态和瞬态研究)转换为有效的 B-H 或 H-B 曲线。这篇文章我们将讨论如何在频域仿真中使用这个仿真 App。 频域中的非线性磁性材料 一个常见的建模假设是在本构关系中指定线
    业界领先的多物理场仿真软件COMSOL Multiphysics® 最新版本推出用于仿真数据管理的平台通用功能,以及用于不确定性量化分析的新模块,并为已有产品带来全面的功能更新和性能提升。 美国马萨诸塞州,伯灵顿(2021年12月14日)—— 业界领先的多物理场仿真解决方案提供商 COMSOL 公司发布了全新的COMSOL Multiphysics® 6.0版本。 新版本推出的“模型管理器”为 C
    3 总结和展望 在锂离子电池的研究中,仍存在许多科学问题尚未解决,这些问题严重影响着锂离子电池的安全性能和使用寿命。例如,锂枝晶的生长演化、SEI膜的形成和破裂演化、正极颗粒在循环中的破裂、电池寿命预测、热失控、以及电池组的电池状态实施监测和管理等问题。这些问题涉及到电场、浓度场、力场和温度场等多个物理场之间的耦合,很难通过单一的实验表征手段对各个驱动力进行分别观测,更难以给出多场耦合的综合结果。
    如果你抬头看看夜空,特别是在远离城市灯光的地方,就会发现很多闪烁的星星。古往今来,美丽的星空是诗人和浪漫主义者的灵感来源,但也一直是研究夜空和宇宙的天文学家们面临的一个挑战。它使自适应光学得以发展。如今,自适应光学被用来提高光学系统的功率,以消除光学介质所带来的不良影响。 自适应光学:消除光学介质带来的不必要影响 当一束来自天文物体(例如恒星)的光,穿过地球大气层时,大气层的湍流会导致这些恒星得图
    摘 要:作为一种具有前景的能量存储系统,锂离子电池需要进一步提高能量密度、功率密度、可靠性和循环稳定性,以满足不断增长的大型能源存储、电动汽车和便携式电子设备需求。当前对锂离子电池的实验研究仍然面临多个挑战,这些挑战包括电解液的导电性和安全性、高能量负极的沉积-剥离机制的优化、高能量正极的循环电压和容量维持、高电流条件下的界面极化和容量释放,以及在极端电流-温度-针刺条件下的热失控管理等问题。这些
    影响力
    粉丝
    内容
    获赞
    收藏
      项目客服
      培训客服
      0 6