求解潜水含水层剖面二维稳定流问题探讨

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

Unconfined aquifer sectional two-dimensional steady flow problem discussion

Based on the analyses of the methods and their characters that solve sectional two--dimensional steady flow in unconfined aquifer, this paper discusses new method for the problem. The new method proposed is only needed to make once mesh grid, retains the node serial number fixed, and generalizes the element shapes crossed by the shallow surface to pentagon, and quadrangle, triangle. A FORTRAN program is developed for each shapes, and the programs can accurately calculate the element transmit matrix. The matrix of the element fully above the water table does not need to be calculated. This paper presents detailed process and arithmetic steps, calculates the position of the shallow surface, compares the numerical value with the analytic value,although there is error between the numerical value and analytic value, the result is satisfied.