Coupled and partially coupled Eulerian‐Lagrangian Model of freshwater‐seawater mixing

The problem of density-dependent transport of salt in unconfined coastal aquifers is solved numerically by means of an implicit Eulerian-Lagrangian finite element formulation. Such a formulation leads to symmetric positive definite finite element matrices which are ideally suited for efficient solution by preconditioned conjugate gradient methods. Additional known advantages of the formulation are unconditional stability, reduced numerical dispersion and suitability for parallel computation. The method has been used to study the effect of dewatering on seawater intrusion within a vertical cross section through an aquifer in southern Italy, related to the construction of a thermoelectric power plant. To investigate the extent to which the dependence of fluid density on salt concentration affects the numerical solution, the flow and advection-dispersion equations were solved in coupled (iterative), partially coupled (noniterative) and completely decoupled modes. Partial coupling was found to yield results very close to those obtained by full coupling but at great savings in computer time; the less rigorous decoupled approach led to results substantially different from those obtained through coupling and partial coupling. Effects of aquifer heterogeneity and the construction of a cutoff wall on seawater intrusion are discussed.

[1]  H. H. Cooper,et al.  Sea water in coastal aquifers , 1964 .

[2]  D. W. Peaceman,et al.  Numerical Calculation of Multidimensional Miscible Displacement by the Method of Characteristics , 1964 .

[3]  S. Schmork,et al.  Upconing of Fresh Water—Sea Water Interface Below Pumping Wells, Field Study , 1969 .

[4]  S. P. Neuman,et al.  Finite Element Method of Analyzing Steady Seepage with a Free Surface , 1970 .

[5]  G. Pinder,et al.  A Numerical Technique for Calculating the Transient Position of the Saltwater Front , 1970 .

[6]  M. A. Collins,et al.  Seawater Intrusion in Layered Aquifers , 1971 .

[7]  R. T. Cheng,et al.  On seawater encroachment in coastal aquifers , 1974 .

[8]  J. Bear,et al.  The shape of the interface in steady flow in a stratified aquifer , 1974 .

[9]  W. G. Gray,et al.  A Galerkin‐finite element technique for calculating the transient position of the saltwater front , 1975 .

[10]  G. Gambolati.  Transient free surface flow to a well: An analysis of theoretical solutions , 1976 .

[11]  George F. Pinder,et al.  Transient Simulation of Saltwater Intrusion in Southeastern Florida , 1976 .

[12]  Ad McCowan.  Numerical Simulation of Shallow Water Waves , 1978 .

[13]  J. Bredehoeft,et al.  Computer model of two-dimensional solute transport and dispersion in ground water , 1978 .

[14]  John L. Wilson,et al.  A numerical model of seawater intrusion in aquifers , 1979 .

[15]  Charles R. Faust,et al.  Simulation of Salt‐Water Interface Motion , 1980 .

[16]  Gour-Tsyh Yeh,et al.  On the computation of Darcian velocity and mass balance in the finite element modeling of groundwater flow , 1981 .

[17]  James A. Liggett,et al.  Boundary integral equation solutions to moving interface between two fluids in porous media , 1981 .

[18]  S. P. Neuman,et al.  A Eulerian-Lagrangian numerical scheme for the dispersion-convection equation using conjugate space-time grids , 1981 .

[19]  E. Frind.  Simulation of long-term transient density-dependent transport in groundwater , 1982 .

[20]  E. Frind.  Seawater intrusion in continuous coastal aquifer-aquitard systems , 1982 .

[21]  K. Rushton,et al.  An assessment of the importance of some parameters for seawater intrusion in aquifers and a comparison of dispersive and sharp-interface modelling approaches , 1982 .

[22]  T. F. Russell,et al.  NUMERICAL METHODS FOR CONVECTION-DOMINATED DIFFUSION PROBLEMS BASED ON COMBINING THE METHOD OF CHARACTERISTICS WITH FINITE ELEMENT OR FINITE DIFFERENCE PROCEDURES* , 1982 .

[23]  The Conjugate Gradients in Subsurface Flow and Land Subsidence Modelling , 1984 .

[24]  James A. Liggett,et al.  Boundary Integral Solution to Seawater Intrusion Into Coastal Aquifers , 1984 .

[25]  S. P. Neuman.  Adaptive Eulerian–Lagrangian finite element method for advection–dispersion , 1984 .

[26]  C. Voss,et al.  SUTRA (Saturated-Unsaturated Transport). A Finite-Element Simulation Model for Saturated-Unsaturated, Fluid-Density-Dependent Ground-Water Flow with Energy Transport or Chemically-Reactive Single-Species Solute Transport. , 1984 .

[27]  K. Stolzenbach,et al.  The 2-D, Unsteady, Transport Equation Solved by the Combined Use of the Finite Element Method and the Method of Characteristics , 1984 .

[28]  Randall J. Charbeneau,et al.  Saltwater upconing in unconfined aquifers , 1985 .

[29]  W. Sanford,et al.  A two-constituent solute-transport model for ground water having variable density , 1985 .

[30]  B. Thomas,et al.  Runge-Kutta methods , 2007, Scholarpedia.

[31]  W. R. Souza,et al.  Variable density flow and solute transport simulation of regional aquifers containing a narrow freshwater‐saltwater transition zone , 1987 .

[32]  Clifford I. Voss,et al.  Analysis of an anisotropic coastal aquifer system using variable-density flow and solute transport simulation , 1987 .

[33]  James R. Fair,et al.  Applied numerical methods with personal computers , 1987 .

[34]  P. F. Andersen,et al.  Saltwater intrusion in aquifers: Development and testing of a three‐dimensional finite element model , 1987 .

[35]  M. Hill.  A Comparison of Coupled Freshwater-Saltwater Sharp-Interface and Convective-Dispersive Models of Saltwater Intrusion in a Layered Aquifer System , 1988 .

[36]  The behavior of groundwater with dispersion in coastal aquifers , 1988 .

[37]  S. P. Neuman,et al.  Three-Dimensional Adaptive Eulerian-Lagrangian Finite Element Method for Advection-Dispersion , 1988 .

[38]  H.-J. Diersch,et al.  Finite element modelling of recirculating density-driven saltwater intrusion processes in groundwater , 1988 .

[39]  M. Sherif,et al.  A two-dimensional finite element model for dispersion (2D-FED) in coastal aquifers , 1988 .

[40]  George F. Pinder,et al.  Can the Sharp Interface Salt-Water Model Capture Transient Behavior? , 1988 .

[41]  S. Majid Hassanizadeh,et al.  On the modeling of brine transport in porous media , 1988 .

[42]  A Comparative Analysis of Mathematical Mass Transport Codes for Groundwater Pollution Studies , 1988 .

[43]  Giuseppe Gambolati,et al.  On boundary conditions and point sources in the finite element integration of the transport equation , 1989 .

[44]  H. Essaid.  A multilayered sharp interface model of coupled freshwater and saltwater flow in coastal systems: Model development and application , 1990 .

[45]  M. A. Celia,et al.  A Eulerian-Lagrangian localized adjoint method for reactive transport in groundwater. , 1990 .