热-孔隙弹性系统的混合化间断伽辽金方法

《Applied Numerical Mathematics》：Hybridizable discontinuous galerkin methods for thermo-poroelastic systems

　　

本文提出了一种用于全动态线性热-孔隙弹性问题的高阶混合化间断伽辽金（HDG）格式。控制方程被表述为一阶双曲系统，其中包含固体速度、流体速度、热通量、有效应力、孔隙压力和温度作为状态变量。研究人员利用半群理论建立了连续问题的适定性，并发展了一种能量一致的HDG离散格式。该方法利用了HDG的计算优势，包括局部性和静态凝聚，同时保持了耦合系统的能量守恒。研究人员给出了hp收敛性分析，并通过全面的数值实验予以验证，确认了理论收敛速率，并展示了该方法在非均匀介质中热-孔隙弹性波传播的有效性。

**论文解读：热-孔隙弹性系统的混合化间断伽辽金方法**

**研究背景与问题**

热-孔隙弹性耦合了固体力学变形、孔隙流体流动和热传输三个物理过程，在地热能开采、核废料存储、水力压裂和强化采油等众多工程与地球物理应用中至关重要。然而，该系统的控制方程涉及多个时间尺度，且材料非均匀性或局部源会导致解出现尖锐梯度，精确数值模拟面临重大挑战。经典计算方法（如有限差分、有限元、边界元、有限体积和混合方法）通常需要极细的网格才能充分分辨这些特征，计算开销巨大。高阶方法（如谱方法和高阶时空伽辽金格式）虽能缓解部分问题，但先进的多面体间断伽辽金（DG）方案虽然具有局部守恒、几何灵活性和自然的hp自适应能力，其较大的计算模板导致每个自由度成本较高。混合化间断伽辽金（HDG）方法通过将单元未知量静态凝聚到面上，在保持标准DG方法精度的同时，产生更小的全局系统。此前HDG已成功应用于弹性和孔隙弹性问题，但尚未有HDG策略被发展用于热-孔隙弹性系统。因此，本文旨在填补这一空白，将HDG框架扩展至更复杂的热-孔隙弹性情形。

**研究内容与结论**

研究人员以Bonetti等人和Carcione等人提出的线性全动态热-孔隙弹性模型为出发点，该模型耦合了固体力学、达西流和有限速度热传输，其中热传导遵循麦克斯韦-卡塔尼奥（Maxwell-Cattaneo）定律，消除了经典傅里叶扩散中非物理的无限传播速度。与既有文献中将问题表述为时间二阶偏微分方程不同，研究人员将控制方程重新表述为一阶双曲系统，这便于利用半群理论建立连续问题的适定性。基于此，研究人员发展了一种能量一致的HDG离散格式，适用于一般单纯形网格，并给出了hp收敛性分析。数值实验验证了理论收敛速率，并展示了该方法在非均匀介质中波传播的有效性，包括与纯孔隙弹性模型的对比，突出了热耦合效应的重要性。

**主要技术方法**

本文采用的核心技术方法包括：1）将热-孔隙弹性控制方程转化为一阶双曲系统，引入状态变量（固体速度、流体速度、热通量、有效应力、孔隙压力、温度）；2）利用强连续半群理论（Lumer-Phillips定理）证明连续问题的适定性；3）构造HDG空间离散，采用分片多项式空间（速度与热通量用k+1次，应力与标量用k次），并引入数值通量进行静态凝聚；4）发展能量一致的半离散格式，通过稳定性分析和先验误差估计获得hp收敛速率；5）数值实验采用制造解验证收敛阶，并模拟均匀和非均匀介质中的波传播问题。样本来源为文献[14]和[9]中的典型地球物理参数，未涉及实际实验样本。

**研究结果**

**6.1 收敛速率验证**：通过制造解（精确解为三角函数组合）在方形域上测试，采用Crank-Nicolson时间离散，结果显示应力、孔隙压力和温度误差达到O(h^k+1)的最优阶，速度误差达到O(h^k+2)，与理论预测一致。在近不可压缩和小存储系数（c₀=10^-6）的极端参数下，方法仍保持鲁棒性。固定网格和时间步长，改变多项式阶数k从1到7，观察到指数收敛。时间步长细化实验确认了O(Δt²)的二阶时间精度。

**6.2 热弹性介质中的波传播**：在1500m×1500m的均匀域中模拟剪切源激发，源用高斯正则化狄拉克函数近似。结果显示速度场呈现特征的四极辐射模式，清晰分辨出快压缩波（P波）和剪切波（S波），未见慢P波（Biot波）符合粘性流体理论预期。温度场展现了传播速度显著慢于弹性波的热波（T波），与参考PolyDG方法结果定性一致。

**6.3 孔隙弹性与热-孔隙弹性波传播对比**：在三个接收点（两个主方向、一个对角方向）比较固体速度和过滤速度的时程。主方向接收点两种模型差异较小，但对角方向接收点热-孔隙弹性模型出现明显幅值和相位偏移，表明热效应对剪切波主导区域影响显著，纯孔隙弹性模型不足以准确预测完整波场。

**6.4 非均匀热-孔隙弹性介质中的波传播**：将域分为左右两半，右侧采用更刚性的材料参数。结果显示材料界面破坏了对称波型，产生复杂的反射和透射现象，以及沿边界传播的头波。热通量和温度场在界面处出现尖锐不连续和增强的边界效应，热耦合系数差异导致沿界面的优先热通道，形成持续存在的边界层效应。

**讨论与结论**

本文建立的HDG格式为全动态热-孔隙弹性问题提供了稳健且计算高效的工具。通过半群理论确立了连续问题的适定性，并证明了方法的最优hp误差估计。数值实验验证了理论收敛速率，并在非均匀介质中有效捕捉了P波、S波和热波的传播。该框架可推广至地质力学、生物力学和材料科学中的多物理场模拟。未来工作将扩展至非线性全耦合热-孔隙弹性问题，特别是包含能量方程中对流输运项的系统。