Premium
Motion of the Seawater Interface in Coastal Aquifers: A Numerical Solution
Author(s) -
Shamir U.,
Dagan G.
Publication year - 1971
Publication title -
water resources research
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.863
H-Index - 217
eISSN - 1944-7973
pISSN - 0043-1397
DOI - 10.1029/wr007i003p00644
Subject(s) - phreatic , aquifer , geology , partial differential equation , linearization , method of characteristics , seawater intrusion , hydraulic head , mathematics , mechanics , mathematical analysis , geotechnical engineering , groundwater , nonlinear system , physics , quantum mechanics
The partial differential equations that describe the motion of the seawater interface and the free surface in a phreatic coastal aquifer (or the freshwater head replacing the latter, in the confined case) are presented. They are based on the Dupuit approximation and take into consideration the geometry of the vertical section through the aquifer, in whose plane the flow takes place, as well as the spatial variation of properties of the porous medium and the spatial and temporal distributions of accretion, recharge, and pumping. An implicit numerical scheme is presented to solve the set of simultaneous partial differential equations. The scheme is based on a linearization of the equations and employs a grid with one spacing over the intrusion length and a different spacing in the remainder of the field. Efficient solution of the resulting set of simultaneous linear equations for each time step is achieved by arranging them in a way that results in a 7 diagonal coefficient matrix. Examples are presented, for which the numerical solutions are compared with analytical solutions or laboratory experiments.