土壤溶质迁移方程的迎风稳定有限点法

针对土壤溶质迁移方程具有对流项的特性,构造了迎风稳定有限点方法.该方法采用自适应迎风格式使其支持域偏向于迎风侧,以获取到上游信息,在对流占优的情况下,避免数值震荡现象.通过对一维和二维土壤溶质迁移方程的数值计算,详细分析了在不同布点、不同时间步长、不同影响因子的作用下,新算法数值解的收敛性、收敛阶与稳定性.数值结果表明本文方法在边界及梯度变化大的区域内,可以有效地提高计算精度,达到消除数值震荡的目的.