A spectral method is developed to study the steady and unsteady flow of fluid into a line sink from a horizontally confined aquifer, and the results are compared to solutions obtained implementing the finite element package COMSOLTM. The aquifer or drain is considered to be confined below so that the solutions are fundamentally unsteady. Comparison is made between the two methods in determining the drawdown of the surface.

The withdrawal of water with a free surface through a line sink from a two-dimensional, vertical sand column is considered using the hodograph method and a novel spectral method. Hodograph solutions are presented for slow flow and for critical, limiting steady flows, and these are compared with spectral solutions to the steady problem. The spectral method is then extended to obtain unsteady solutions and hence the evolution of the phreatic surface to the steady solutions when they exist. It is found that for each height of the interface there is a unique critical coning value of flow rate, but also that the value obtained is dependent on the flow history.

The steady response of a fluid with two layers of different density in a porous medium is considered during extraction through a point sink. Supercritical withdrawal in which both layers are being withdrawn is investigated using a spectral method. We show that for each withdrawal rate, there is a single entry angle of the interface into the point sink. As the flow rate decreases the angle of entry steepens until it becomes almost vertical, at which point the method fails. This limit is shown to correspond to the upper bound on sub-critical (single-layer) flow.