饱和多孔弹性杆热传导的广义多辛方法及其数值实现

首先根据多孔介质理论,利用饱和多孔介质的能量方程和本构关系,推导出饱和多孔弹性杆局部热平衡的热传导方程;继而引入正交变量,将热传导方程导入Hamilton系统,得到饱和多孔弹性杆热传导方程的广义多辛形式和多种局部守恒律形式;接着采用中点离散方法对热传导方程的广义多辛形式进行数值离散;最后利用计算机数值实现了饱和多孔弹性杆的热传导过程,并且讨论了参数取值的不同对热传导过程的影响,同时在数值模拟过程中记录了广义多辛格式的局部动量误差。研究结果表明,构造的广义多辛方法能够很好地模拟系统的热传导过程和耗散效应,同时也可长时间保持系统的固有几何性质。