Numerical Simulation of Salt-Water Intrusion in Coastal Aquifers. Part 1.

摘要

The set of non-linear partial differential equations which describe the movement of the saltwater front in a coastal aquifer is solved by the Galerkin-finite element method. Pressure and velocities are obtained simultaneously in order to guarantee continuity of velocities between elements. A layered aquifer is modeled either with a functional representation of permeability or by a constant value of permeability over each element. The Galerkin-finite element method is subsequently used to simulate the two-dimensional movement of the saltwater front in the Cutler area of the Biscayne aquifer near Miami,Florida. The concentration distributions for steady-state and transient conditions are obtained using pressure head data measured in the field in formulating boundary conditions.