海洋试点国家实验室海洋矿产资源评价与探测技术功能实验室科研人员在水合物多场耦合数值模拟方法研究方面取得新进展，提出了一种基于控制体有限元法和有限元混合的流-固耦合数值模拟方法，为水合物开采的产能和力学稳定性评价提供了一种新的数值模拟研究手段。
天然气水合物的开采是一个涉及渗流、传热、水合物相变及沉积物力学变形的THMC多物理场耦合过程。从研究对象来看，其中渗流、传热和水合物相变可视作流体力学问题，而沉积物变形则是固体力学问题。由于研究对象的不同，流体和固体系统应当使用不同的数值方法来处理。目前，在水合物THMC耦合数值模拟方法有两种研究思路：第一种是使用有限体积法计算流体系统，而使用有限元法处理固体系统；第二种思路则是全部使用有限元方法同时处理流体系统和固体系统。然而，第一种方法在流固耦合过程中存在需要插值以及流体系统的网格不适用于固体系统的问题；第二种方法则存在着流体系统局部不守恒从而造成不收敛和非物理意义解的问题。
本研究提出了一种控制体有限元和有限元混合的新思路，在全非结构网格的基础上构建一套对偶的控制体网格，在控制体网格上用有限体积法结合有限元的插值函数进行流体系统的求解，而固体系统则在原网格上使用有限元法求解。这样一套新的研究思路有几个方面的优势：①构建的控制体使得流体系统和固体系统的求解变量位于相同的节点上，在保证流体系统求解局部守恒性的同时避免了流-固耦合的插值和网格不适应的问题；②引入有限元插值函数进行多点流量近似(Multiple-point flux approximation, MPFA)，避免了传统有限体积法对网格正交性的要求，可以使用完全非结构网格构建复杂的几何模型。基于上述算法，该团队开发了一套具有完全自主知识产权的数值模拟器QIMGHyd-THMC，利用室内模拟实验的结果和国际主流数值模拟软件对该模拟器进行了验证，为水合物多场耦合数值模拟方法提供了一种新的思路和框架。
控制体有限元和有限元混合的流固耦合模拟方法及在复杂结构井中的应用
开发的QIMGHyd-THMC模拟器与国际主流模拟器计算结果的对比
成果以“Coupled thermal-hydrodynamic-mechanical–chemical numerical simulation for gas production from hydrate-bearing sediments based on hybrid finite volume and finite element method（基于控制体有限元(CVFEM)和有限元(FEM)混合的天然气水合物开采渗流、传热、相变和变形的多场耦合数值模拟）”为题在线发表于岩土工程领域顶级期刊Computers and Geotechnics(计算机和岩土工程)，海洋矿产资源评价与探测技术功能实验室万义钊副研究员为第一作者，吴能友研究员为通讯作者。该研究得到了青岛海洋科学与技术试点国家实验室海洋矿产资源评价与探测技术功能实验室自主课题(MMRZZ201807)、国家自然科学基金(41906187、41976205)、山东省自然科学基金(ZR2019BD058)和泰山学者特聘专家计划(ts201712079)等共同资助。
文章链接：https://www.sciencedirect.com/science/article/pii/S0266352X2200057X