非饱和冻土水汽迁移与相变过程的光滑粒子法模拟

为了研究多年冻土表层的水热分布情况,在非饱和冻土的能量守恒方程和水分迁移的质量控制方程的基础上考虑冰水相变和水汽相变过程,并考虑水汽运移传热及温度势对水汽迁移的影响,建立了非饱和冻土的水-热-汽耦合模型。采用光滑粒子流体动力学(smoothed particle hydrodynamics,简称SPH)方法可方便地计算它们的演化过程。为此,在计算中先求解能量守恒方程的含冰量及气态水含量,再对未冻水含量和温度场进行求解,从而实现了温度场与水汽场的耦合。在此基础上,模拟计算了第1类热边界条件下半无限空间介质内非稳态温度场、体积含水率及水汽通量的分布情况,并将计算结果与未考虑耦合的解析解进行比较,结果显示水汽耦合的作用不容忽略。最后,针对处于季节性周期温度边界下路基的水热场的分布情况进行计算。研究表明,相比于水-热耦合模型,所建立的水-热-汽耦合模型得到的计算结果更为接近实际监测结果,可很好地揭示非饱和冻土中的水热汽迁移特征及其相变过程。