论文部分内容阅读
本文在对现有求解潜水含水层剖面二维流的多种方法及其特点进行分析的基础上，探讨新的求解方法。提出的方法只需进行一次网格剖分，在网格节点编号保持不变的情况下，把潜水面与单元相交的形状概化为五边形、任意四边形和三角形，按五边形、任意四边形和三角形子程序准确地求出潜水面穿越的单元的渗透矩阵。同时不用计算完全在潜水面之上的单元的渗透矩阵。本文给出具体的思路、算法步骤，计算潜水面的位置，通过比较数值解与解析解，两者虽然存在一定的误差，但是该方法总体结果是比较令人满意的。