A Semidiscrete Model for Water and Solute Movement in Tile‐Drained Soils: 1. Governing Equations and Solution

A finite element model has been developed to simulate solute transport in tile-drained soil-aquifer systems. Water flow in the unsaturated zone and to drains in the saturated zone was assumed to be at steady state. The model considers the transport of nonreactive solutes, as well as of reactive solutes whose behavior can be described by a distribution coefficient. The exact-in-time numerical solution yields explicit expressions for the concentration field at any future point in time without having to compute concentrations at intermediate times. The semidiscrete method involves the determination of an eigensystem of eigenvalues and eigenvectors of the coefficient matrix. The eigensystem may be complex (i.e., it may have imaginary components) due to asymmetry created by the convection term in the governing convection-dispersion equation. The proposed approach facilitates long-term predictions of concentrations in drainage effluents and of salt distributions in soil and groundwater. The accuracy of the model was verified by comparing model results with those based on an analytical solution for two-dimensional solute transport in groundwater.

[1]  R. Genuchten CALCULATING TEE UNSATURATED HYDRAULIC CONDUCTIVITY WITH A NEW CLOSED-FORM ANALYTICAL MODEL , 1978 .

[2]  The Problem of Complex Eigensystems in the Semianalytical Solution For Advancement of Time in Solute Transport Simulations: A New Method Using Real Arithmetic , 1986 .

[3]  Analytic solution of spatially discretized groundwater flow equations , 1973 .

[4]  W. G. Gray,et al.  Finite Element Simulation in Surface and Subsurface Hydrology , 1977 .

[5]  A General Numerical Solution of the Two‐Dimensional Diffusion‐Convection Equation by the Finite Element Method , 1970 .

[6]  Martin Braun Differential equations and their applications , 1976 .

[7]  D. K. Sunada,et al.  Finite element method for the hydrodynamic dispersion equation with mixed partial derivatives , 1972 .

[8]  G. Yeh,et al.  An eigenvalue solution continuous in time to the spatially discretized solute transport equation in steady groundwater flow. , 1984 .

[9]  D. Suarez Impact of agricultural practices on groundwater salinity , 1989 .

[10]  R. Allan Freeze,et al.  Three-Dimensional, Transient, Saturated-Unsaturated Flow in a Groundwater Basin , 1971 .

[11]  J. O. Duguid,et al.  Material transport through porous media: a finite-element Galerkin model , 1976 .

[12]  Ian P. King,et al.  Salinity Management Model: I. Development , 1987 .

[14]  Richard Bellman,et al.  Introduction to Matrix Analysis , 1972 .

[15]  Two-Dimensional Dispersion , 1967 .

[16]  Robert W. Gillham,et al.  Finite-element analysis of the transport of water and solutes in title-drained soils , 1979 .

[17]  A. Scheidegger General Theory of Dispersion in Porous Media , 1961 .

[18]  D. Kirkham,et al.  STEADY DRAINAGE OF LAYERED SOILS: I, THEORY , 1971 .

[19]  G. L. Guymon,et al.  A finite element solution of the one-dimensional diffusion-convection equation , 1970 .

[20]  Gideon Kruseman,et al.  Water table control, reuse and disposal of drainage water in Haryana , 1988 .

[21]  R. Skaggs,et al.  Experimental evaluation of theoretical solutions for subsurface drainage and irrigation , 1977 .

[22]  D. Kirkham,et al.  Seepage of steady rainfall through soil into drains , 1958 .

[23]  H. Saunders,et al.  Elementary Finite Element Method , 1978 .

[24]  S. P. Neuman,et al.  SATURATED-UNSATURATED SEEPAGE BY FINITE ELEMENTS , 1973 .

[25]  M. Todsen Numerical studies of two-dimensional saturated/ unsaturated drainage models , 1973 .

[26]  D. S. Ward,et al.  FEMWASTE: a Finite-Element Model of Waste transport through porous saturated-unsaturated media , 1981 .

[27]  P. J. Wierenga,et al.  Solute Distribution Profiles Computed with Steady‐State and Transient Water Movement Models , 1977 .

[28]  Andrés Sahuquillo,et al.  An eigenvalue numerical technique for solving unsteady linear groundwater models continuously in time , 1983 .