基于COMSOL注气驱替瓦斯的多物理问题耦合分析
一、热流固thm耦合下注co2/n2驱替瓦斯多场耦合
一般来说，注co2/n2驱替瓦斯多场耦合涉及流固耦合、热流固耦合及热流固化耦合等多种耦合形式，本视频主要介绍热流固(thm)多场耦合方法，以期为相关研究方向的老师学生提供参考。此方法具有以下几方面的技术难点：（1）煤层中注入co2，涉及到双组分（co2/n2和ch4），气体运移更复杂；（2）煤层的渗透率、孔隙率方程增添由co2吸附扩散引起的变化项；（3）涉及到的物理场增多，方程更复杂，数值求解中模型收敛性存在一定困难。
本文构建的物理场方程来自于公开发表的文献，对于具体的数值求解方法，限于篇幅，会做出一部分解释，主要从co2-ecbm的机理角度出发。
首先构建模型的物理场方程，如图1。该物理场方程主要分为气体扩散对流方程、温度场方程、煤体变形控制方程，其中还有一些辅助方程，如渗透率方程、孔隙率方程等。煤体的有效应力方程考虑了基质、裂隙中的孔压，基质变形引起的应力、煤层温度变化引起的热应力。同时在煤体变形控制方程中，考虑有效应力变化的煤体变形方程。煤体的对流扩散方程分为扩散项、对流项。
此过程，将裂隙和基质假设为一个整体，在这个整体上获得统一的co2与甲烷的对流扩散方程，其中该系统的源项为0。温度场需要考虑煤层本身的传热以及内部对流换热与基质、煤体变形引起的温度变化。
将三个物理场方程耦合解算，是该数值模拟的一个难点。
本案列选择多物理场求解工具comsol，其在多场求解方面广泛应用。
图1 co2-ecbm多场耦合物理场方程
二、comsol设置简要介绍
comsol中求解步骤主要为参数、变量设置，几何模型设置，物理场设置，网格划分，求解器设置，后处理。参数变量设置中，需要把co2-ecbm耦合方程中，相关的参数、变量设置到全局参数中。同时把一些物理场方程用到的变量设置到局部变量中。几何模型中，采用1/4煤层进行构建，在中心位置设置注气孔，在右上角设置抽采孔。 物理场选择2个系数型偏微分方程和1个一般形式偏微分方程，其中co2与甲烷的对流扩散方程、温度场方程采用系数型偏微分方程，煤层变形控制方程采用一般形式偏微分方程。
对流扩散方程的边界条件设置中，在注气孔边界只设置co2注气边界，在抽采边界只设置甲烷抽采边界。煤层变形控制方程中，需要设置对称边界，即零通量。在上边界设置应力载荷，在右边界设置位移边界。煤层变形控制方程在固体力学中设置，也可以使用pde模块，编写相应表达式。温度场可在多孔介质传热方程设置，也可以用pde方程变形相应表达式。 本模型全部选用pde方程编写表达式求解，其好处在于可在同一求解器中求解，方便方程收敛。求解器采用全耦合隐式算法，采用自动(牛顿)非线性方法终止。
图2 参数、几何模型设置
后处理设置：后处理主要展示煤层渗透率、co2，甲烷的压力、煤层应力、位移变化等。
MAX7219驱动共阴极LED原理图及程序
智能手机的电路结构图分享
诺基亚8上市时间和外形照片曝光，情怀+实力或有惊喜！
智能电销机器人vs电销员工，5项能力PK
导热塑料散热能力超越金属铜、铝之10倍以上
基于COMSOL注气驱替瓦斯的多物理问题耦合分析
低噪声高精度磁头放大器(AD797)
移动电视新模式发展方向
小米7效果图曝光：前方高能，喜欢存粮啦！
明天正式首售！OPPO A93首批搭载高通骁龙480处理器
苹果推出的Airtag究竟能否带火UWB？
移动设备性能与健康状况：安卓设备故障率高于IOS，三星是最容易坏的
开发者圆桌会议：HTML5误解
三星GalaxyFold折叠屏手机高清图赏
SDSoC就绪，它如何在实际的设计开发工作中支持你的“表演”
通过弹片微针模组来完美解决FPC测试的难题
如何计算方管的重量，它的计算公式是怎样的
苹果A11芯片将采用10nm制程 A12/A13采用台积电7nm
什么是嵌入式软件开发
OPPO Reno5 Pro究竟值不值得买
一般来说，注co2/n2驱替瓦斯多场耦合涉及流固耦合、热流固耦合及热流固化耦合等多种耦合形式，本视频主要介绍热流固(thm)多场耦合方法，以期为相关研究方向的老师学生提供参考。此方法具有以下几方面的技术难点：（1）煤层中注入co2，涉及到双组分（co2/n2和ch4），气体运移更复杂；（2）煤层的渗透率、孔隙率方程增添由co2吸附扩散引起的变化项；（3）涉及到的物理场增多，方程更复杂，数值求解中模型收敛性存在一定困难。
本文构建的物理场方程来自于公开发表的文献，对于具体的数值求解方法，限于篇幅，会做出一部分解释，主要从co2-ecbm的机理角度出发。
首先构建模型的物理场方程，如图1。该物理场方程主要分为气体扩散对流方程、温度场方程、煤体变形控制方程，其中还有一些辅助方程，如渗透率方程、孔隙率方程等。煤体的有效应力方程考虑了基质、裂隙中的孔压，基质变形引起的应力、煤层温度变化引起的热应力。同时在煤体变形控制方程中，考虑有效应力变化的煤体变形方程。煤体的对流扩散方程分为扩散项、对流项。
此过程，将裂隙和基质假设为一个整体，在这个整体上获得统一的co2与甲烷的对流扩散方程，其中该系统的源项为0。温度场需要考虑煤层本身的传热以及内部对流换热与基质、煤体变形引起的温度变化。
将三个物理场方程耦合解算，是该数值模拟的一个难点。
本案列选择多物理场求解工具comsol，其在多场求解方面广泛应用。
图1 co2-ecbm多场耦合物理场方程
二、comsol设置简要介绍
comsol中求解步骤主要为参数、变量设置，几何模型设置，物理场设置，网格划分，求解器设置，后处理。参数变量设置中，需要把co2-ecbm耦合方程中，相关的参数、变量设置到全局参数中。同时把一些物理场方程用到的变量设置到局部变量中。几何模型中，采用1/4煤层进行构建，在中心位置设置注气孔，在右上角设置抽采孔。 物理场选择2个系数型偏微分方程和1个一般形式偏微分方程，其中co2与甲烷的对流扩散方程、温度场方程采用系数型偏微分方程，煤层变形控制方程采用一般形式偏微分方程。
对流扩散方程的边界条件设置中，在注气孔边界只设置co2注气边界，在抽采边界只设置甲烷抽采边界。煤层变形控制方程中，需要设置对称边界，即零通量。在上边界设置应力载荷，在右边界设置位移边界。煤层变形控制方程在固体力学中设置，也可以使用pde模块，编写相应表达式。温度场可在多孔介质传热方程设置，也可以用pde方程变形相应表达式。 本模型全部选用pde方程编写表达式求解，其好处在于可在同一求解器中求解，方便方程收敛。求解器采用全耦合隐式算法，采用自动(牛顿)非线性方法终止。
图2 参数、几何模型设置
后处理设置：后处理主要展示煤层渗透率、co2，甲烷的压力、煤层应力、位移变化等。
MAX7219驱动共阴极LED原理图及程序
智能手机的电路结构图分享
诺基亚8上市时间和外形照片曝光，情怀+实力或有惊喜！
智能电销机器人vs电销员工，5项能力PK
导热塑料散热能力超越金属铜、铝之10倍以上
基于COMSOL注气驱替瓦斯的多物理问题耦合分析
低噪声高精度磁头放大器(AD797)
移动电视新模式发展方向
小米7效果图曝光：前方高能，喜欢存粮啦！
明天正式首售！OPPO A93首批搭载高通骁龙480处理器
苹果推出的Airtag究竟能否带火UWB？
移动设备性能与健康状况：安卓设备故障率高于IOS，三星是最容易坏的
开发者圆桌会议：HTML5误解
三星GalaxyFold折叠屏手机高清图赏
SDSoC就绪，它如何在实际的设计开发工作中支持你的“表演”
通过弹片微针模组来完美解决FPC测试的难题
如何计算方管的重量，它的计算公式是怎样的
苹果A11芯片将采用10nm制程 A12/A13采用台积电7nm
什么是嵌入式软件开发
OPPO Reno5 Pro究竟值不值得买