Finite-element analysis of the transport of water and solutes in title-drained soils

Abstract The finite-element method based on a Galerkin technique was used to formulate the problem of simulating the two-dimensional (cross-sectional) transient movement of water and solute in saturated or partially saturated nonuniform porous media. The numerical model utilizes linear triangular elements. Nonreactive, as well as reactive solutes whose behaviour can be described by a distribution coefficient or first-order reaction term were considered. The flow portion of the model was tested by comparison of the model results with experimental and finite-difference results for transient flow in an unsaturated sand column and the solute transport portion of the model was tested by comparison with analytical solution results. The model was applied to a hypothetical case involving movement of water and solutes in tile-drained soils. The simulation results showed the development of distinct solute leaching patterns in the soil as drainage proceeded. Although applied to a tile drainage problem in this study, the model should be equally useful in the study of a wide range of two-dimensional water and solute migration problems.

[1]  E. Bresler,et al.  Simultaneous transport of solutes and water under transient unsaturated flow conditions , 1973 .

[2]  G. Pinder,et al.  Functional coefficients in the analysis of groundwater flow , 1973 .

[3]  B. O'CONNOR,et al.  Rayleigh-Ritz and Galerkin finite elements for diffusion-convection problems , 1973 .

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

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

[6]  L. G. King,et al.  DESCRIPTION OF SOIL CHARACTERISTICS FOR PARTIALLY SATURATED FLOW , 1965 .

[7]  J. W. Biggar,et al.  Simultaneous Solute and Water Transfer for an Unsaturated Soil , 1971 .

[8]  J. Duffy,et al.  A Simulation Model of Biophysiochemical Transformations of Nitrogen in Tile-Drained Corn Belt Soil 1 , 1975 .

[9]  D. Russo,et al.  Two-Dimensional Solutes Transfer during Nonsteady Infiltration: Laboratory Test of Mathematical Model 1 , 1975 .

[10]  D. R. Baker,et al.  Simultaneous Transport of Water and Adsorbed Solutes Through Soil Under Transient Flow Conditions , 1975 .

[11]  R. Banks,et al.  A solution of the differential equation of longitudinal dispersion in porous media , 1961 .

[12]  J. Bear.  Dynamics of Fluids in Porous Media , 1975 .

[13]  Howard P. Johnson,et al.  Nitrate, Phosphorus, and Sulfate in Subsurface Drainage Water , 1975 .

[14]  O. Zienkiewicz.  The Finite Element Method In Engineering Science , 1971 .

[15]  W. D. Kemper,et al.  Diffusion of Salts in Clay-Water Systems1 , 1966 .

[16]  W. D. Kemper,et al.  Movement of Nutrients to Plant Roots , 1968 .

[17]  M. Mariño.  Distribution of Contaminants in Porous Media Flow , 1974 .

[18]  E. Bresler,et al.  Two-dimensional Transport of Solutes During Nonsteady Infiltration from a Trickle Source1 , 1975 .

[19]  A. Klute,et al.  Hydraulic Properties of a Porous Medium: Measurement and Empirical Representation 1 , 1976 .

[20]  George F. Pinder,et al.  A Galerkin‐finite element simulation of groundwater contamination on Long Island, New York , 1973 .