基于LTNE的二维饱和多孔介质平板热流固完全耦合数值分析

摘要

针对目前多孔介质热流固耦合数学模型多为不完全耦合模型,且物理模型多为一维或一维轴对称模型的现状,对二维饱和多孔介质平板热流固完全耦合问题进行数学建模和数值分析.采用强耦合方式实现热流固3个物理场的完全耦合.在平板左侧施加30℃流体和固体温度边界条件及固定位移边界条件,平板其他边界温度场和位移场假设为自然边界条件,平板四周加载0 Pa的孔隙压力边界条件.利用COMSOL Multiphysics有限元分析软件的偏微分方程(PDEs)模式实现上述完全耦合模型的求解,得到渗流场、应变场以及双温度场的数值解.数值结果表明,随时间增加,流体和固体温度、应变及孔隙压力沿x轴方向传递,同时发现x轴方向应变远大于y轴方向应变.此外,随时间增加,y轴方向应变和孔隙压力最大值逐渐减小.研究建立的数学模型和数值解有助于深入认识二维饱和多孔介质热流固完全耦合力学行为.