裂隙岩体温度场数值流形方法初步研究

基于三角形有限覆盖数值流形理论,提出了裂隙岩体二维不稳定温度场求解的数值流形方法(NMM),并推导了第一类温度边界条件处理的罚方法,给出了计算格式。数值流形方法采用两套覆盖,其中数学覆盖独立于求解区域,生成数学覆盖时不用考虑裂隙数量、位置和方向,避免了常规方法处理裂隙处网格划分的不便。对于裂隙岩体,裂隙两侧对应不同流形单元,可以实现温度的不连续模拟。裂隙作为内部流形单元或外部边界条件叠加进入到温度场的总体求解矩阵中,实现了裂隙岩体温度场的数值流形求解。该方法具有较高精度,且能求解任意多条裂隙的岩体温度场问题。最后对部分文献算例进行计算,结果具有较好的一致性。