焦　甜,王军霞,唐仲华，刘　耿

(1.中国地质大学(武汉)环境学院，湖北 武汉 430074；2.中国地质大学(武汉)数学与物理学院，湖北 武汉 430074；3.北京得力合环境治理有限公司，北京 100025)



地下水溶质运移方程有限差分格式的实证研究

焦甜1,王军霞2,唐仲华1，刘耿3

(1.中国地质大学(武汉)环境学院，湖北 武汉 430074；2.中国地质大学(武汉)数学与物理学院，湖北 武汉 430074；3.北京得力合环境治理有限公司，北京 100025)

摘要：在求解对流占优的地下水溶质运移方程时，为克服有限差分法和有限单元法产生的过量和值弥散等虚假数值现象，出现了许多数值方法，其中包括Crank-Nicholson格式、上游加权法和人工扩散量法等。以一维地下水溶质运移方程为例，分析三种差分格式的截断误差、稳定性和收敛性，并通过数值试验，在不同的对流强度下，将不同算法的数值解与解析解进行对比，研究不同算法的优缺点和适用条件。结果表明：在对流占优较弱时，三种方法都适用，但Crank-Nicholson差分格式误差最小；当对流强度增大时，适用时间步长较小的人工扩散量法；当对流强度较大时，适用时间步长较大的人工扩散量法，或者空间步长较小的上游加权法。该研究可为解决实际问题时算法的选择提供依据。

关键词：地下水溶质运移方程；对流占优强度；有限差分格式；解析解；实证对比

随着经济的快速发展，地下水污染问题日益严重，由于地下水污染具有滞后性、隐蔽性以及一旦污染便难以清除的特点[1]，因此地下水污染问题得到了广泛关注。……