摘要

光滑粒子流体动力学(SPH)作为一种拉格朗日型无网格粒子方法,已经成功地应用于多种类型耦合问题的求解。基于SPH方法探索求解二维求解域内含有圆形局部热源的水热耦合问题,其中揭示了程序前处理阶段已有粒子布置方案存在的不足。针对所要求解的具体问题,使用一种新的粒子布置方案以弥补数值计算过程中的精度损失,该粒子布置方案中径向粒子间距严格相等,环向间距大致等于径向粒子间距。通过对水热耦合方程数学特性的分析,表明了SPH算法可以用于此类耦合问题的数值模拟。最后,基于SPH算法自行编制的程序对一种非饱和土介质内含有局部热源的水热耦合问题进行数值模拟,对非饱和土介质中热能传输以及水分迁移的规律进行了分析。