Derivation and relative performance of strings of line elements for modeling (un)confined and semi-confined flow

In the analytic element method, strings of line-sinks may be used to model streams and strings of line-doublets may be used to model impermeable walls or boundaries of inhomogeneities. The resulting solutions are analytic, but the boundary conditions are met approximately. Equations for line elements may be derived in two ways: through integration of point elements (the integral solution) and through application of separation of variables in elliptical coordinates (the elliptical solution). Using both approaches, two sets of line elements are presented for four flow problems: line-sinks and line-doublets in (un)confined flow, and line-sinks and line-doublets in semi-confined flow. Elliptical line elements have the advantage that they do not need a far-field expansion for accurate evaluation far away from the element. The derivation of elliptical line elements is new and applicable to both (un)confined flow and semi-confined flow; only the resulting expressions for elliptical line elements for semi-confined flow have not been found in the current groundwater literature. Existing solutions for elliptical line elements for (un)confined flow were intended for the modeling of isolated features. Four examples are presented, one for each flow problem, to demonstrate that strings of elliptical line elements may be used to obtain accurate solutions; elliptical line-doublets for semi-confined flow can only be strung together in combination with two integral line-doublets. For a zigzag canal in (un)confined flow, a string of elliptical line-sinks performed better than a string of integral line-sinks of the same order when the smallest angle between two adjacent segments is less than 130°. Elliptical line-doublets performed better than integral line-doublets for a square inhomogeneity in a uniform, confined flow field; the difference was smaller for an octagonal inhomogeneity. For semi-confined flow, the difference between the integral and elliptical line-sinks was small when modeling a zigzag canal.