解冻结壁温度场的一个数值方法

相变热传导问题是一个弱间断的非线性问题。借助于一种数学变换将原控制微分方程变为空间变量为无量纲的微分方程后,采用三时间层隐式有限差分格式(Bonacina,1973)求数值解是非常方便的。本文用这种方法求解了冻结壁温度场。但这种方法只能用于相变温度等于零且水分迁移可以被忽略的粗颗粒土,对于不同相变温度区间的各类粘土则不能使用。