土体冻结过程中相变导热的光滑粒子法分析

相变导热在矿山工程、城市地铁、河底隧道等领域有着广泛的应用.由于该问题模型方程独特的数学性质及网格背景数值方法的局限性,引入光滑粒子算法求解相变导热问题.给出基于Fortran语言开发的数值分析程序,并通过推导得到的半无限域凝固问题的解析解验证了其精度和可靠性.最后针对低温冻结管作用下土体冻结过程进行了分析,计算结果表明:土体冻结过程中冻结锋面的移动速度随时间逐渐减小,温度场曲线上的"折线"变化是相变时潜热释放引起的;冻结温度对土层冻结的发展影响显著,冻结管的温度越低冻结锋面的移动速度越快;冻土导热系数k_s的取值大小对冻结锋面的移动速度没有显著影响;然而未冻土导热系数k_l取值越大,土层中冻结锋面的移动速度越快.