Modeling of 2D Density-Dependent Flow and Transport in the Subsurface
Publication: Journal of Hydrologic Engineering
Volume 3, Issue 4
Abstract
A 2-dimensional finite-element model for density-dependent flow and transport through saturated-unsaturated porous media has been developed. The combined flow and transport model can handle a wide range of real-world problems, including the simulations of flow alone, contaminant transport alone, and combined flow and transport. The conventional finite-element methods and a hybrid Lagrangian-Eulerian finite-element method were incorporated in the transport module. Saltwater intrusion problems and instability caused by denser water on the top were investigated in this paper. Because the fundamental mechanism causing saltwater intrusion most likely is caused by density-induced convection and dispersion, the developed model was used to assess the interplay between density-driven flow and the subsurface media through which the saltwater intrusion occurs. The mathematical formulation of the model is comprised of fluid flow and solute transport equations, coupled by fluid density. In the specific case of saltwater intrusion and unstable brine transport problems, this set of governing equations is nonlinear and requires iterative methods to solve them simultaneously. Three case studies, which include a wide spectrum of physical conditions, show the verification and effectiveness of the model by comparing previously published solutions from other researchers with the simulation results of the present model. Two demonstrated problems examine the model capabilities to handle saltwater intrusion problems through unsaturated-saturated porous media and density-dependent flow and transport under unstable conditions.
Get full access to this article
View all available purchase options and get full access to this article.
References
1.
Andrews, R. W.(1981). “Salt-water intrusion in the Costa de Hermisillo, Mexico: A numerical analysis of water management proposals.”Ground Water, 19(6), 635–647.
2.
Bear, J. (1979). Hydraulics of groundwater. McGraw-Hill, Inc., New York.
3.
Contractor, D. N., and Srivastava, R.(1990). “Simulation of saltwater intrusion in the northern guam lens using a microcomputer.”J. Hydro., Amsterdam, The Netherlands, 118, 87–106.
4.
Duffy, C. J., and Al-Hassan, S.(1988). “Groundwater circulation in a closed desert basin: Topographic scaling and climatic forcing.”Water Resour. Res., 24(10), 1675–1688.
5.
Frind, E. O.(1982a). “Seawater intrusion in continuous coastal aquifer-aquitard systems.”Adv. Water Resour., 5(2), 89–97.
6.
Frind, E. O.(1982b). “Simulation of long-term transient density-dependent transport in groundwater.”Adv. Water Resour., 5(2), 73–88.
7.
Galeati, G., Gambolati, G., and Neuman, S. P.(1992). “Coupled and partially coupled Eulerian-Lagrangian model of freshwater-seawater mixing.”Water Resour. Res., 28(1), 149–165.
8.
Henry, H. R.(1959). “Salt intrusion into fresh water aquifers.”J. Geophys. Res., 64(11), 1911–1919.
9.
Henry, H. R. (1964). “Effects of dispersion on salt encroachment in coastal aquifers, seawater in coastal aquifers.”U.S. Geological Survey Water Supply Paper No. 1613-C, U.S. Geological Survey, Washington, D.C., 35–69.
10.
Huyakorn, P. S., Andersen, P. F., Mercer, J. W., and White Jr., H. O.(1987). “Saltwater intrusion in aquifers: Development and testing of a three-dimensional finite element model.”Water Resour. Res., 23(2), 293–312.
11.
Huyakorn, P. S., and Pinder, G. F. (1983). Computational methods in subsurface flow. Academic Press, New York.
12.
Huyakorn, P. S., Springer, E. P., Guvanasen, V., and Wadsworth, T. D.(1986). “A three-dimensional finite-element model for simulating water flow in variably saturated porous media.”Water Resour. Res., 22(13), 1790–1808.
13.
Lee, C. H., and Cheng, R. T.(1974). “On seawater encroachment in coastal aquifers.”Water Resour. Res., 10(5), 1039–1043.
14.
Liu, P. L.-F., Cheng, A. H.-D., Liggett, J. A., and Lee, J. H.(1981). “Boundary integral equation solutions to moving interface between two fluids in porous media.”Water Resour. Res., 17(5), 1445–1452.
15.
Pinder, G. F., and Cooper Jr., H. H.(1970). “A numerical technique for calculating the transient position of the saltwater front.”Water Resour. Res., 6(3), 875–880.
16.
Reilly, T. E., and Goodman, A. S.(1985). “Quantitative analysis of saltwater-freshwater relationships in groundwater systems—A historical perspective.”J. Hydro., Amsterdam, The Netherlands, 80, 125–160.
17.
Segol, G., Pinder, G. F., and Gray, W. G.(1975). “A Galerkin finite element technique for calculating the transient position of the salt water front.”Water Resour. Res., 11(2), 343–347.
18.
Shamir, V., and Dagan, G.(1971). “Motion of the seawater interface in coastal aquifers: A numerical solution.”Water Resour. Res., 7(3), 644–657.
19.
Strobl, R. O. (1993). “Modeling of two-dimensional density-dependent flow and transport through saturated-unsaturated porous media,” MS thesis, Department of Civil and Environmental Engineering, Pennsylvania State University, University Park, Pa.
20.
Wilson, J. L., and Sa da Costa, A. A. G.(1982). “Finite element simulation of a saltwater/freshwater interface with indirect toe tracking.”Water Resour. Res., 18(4), 1069–1080.
21.
Yeh, G. T.(1981). “On the computation of Darcian velocity and mass balance in the finite element modeling of groundwater flow.”Water Resour. Res., 17(5), 1529–1534.
22.
Yeh, G. T. (1987). “FEMWATER: A finite element model of water flow through saturated-unsaturated porous media.”First Revision, Rep. No. ORNL-5567/R1, Oak Ridge National Laboratory, Oak Ridge, Tenn.
23.
Yeh, G. T., and Chang, J. R. (1993). Users' manual of a hybrid Lagrangian-Eulerian finite element model of WASTE transport through saturated-unsaturated media: LEWASTE version 2.0. Pennsylvania State University, University Park, Pa.
24.
Yeh, G. T., Chang, J. R., Cheng, H. P., and Sung, C. H.(1995). “An adaptive local grid refinement based on the exact peak capture and oscillation free scheme to solve transport equations.”Comp. and Fluids, 24(3), 293–332.
25.
Yeh, G. T., Chang, J. R., Gwo, J. P., Lin, H. C., Richards, D. R., and Martin, W. D. (1994). “3DSALT: A three-dimensional finite element model of density-dependent flow and transport through saturated-unsaturated media.”Instruction Rep. No. HL-94-1, U.S. Army Corps of Engineers, Vicksburg, Miss.
26.
Yeh, G. T., Chang, J. R., and Short, T. E.(1992). “An exact peak capturing the oscillation free scheme to solve advection-dispersion transport equations.”Water Resour. Res., 28(11), 2937–2951.
Information & Authors
Information
Published In
Journal of Hydrologic Engineering
Volume 3 • Issue 4 • October 1998
Pages: 248 - 257
Copyright
Copyright © 1998 American Society of Civil Engineers.
History
Published online: Oct 1, 1998
Published in print: Oct 1998
Authors
Metrics & Citations
Metrics
Citations
Download citation
If you have the appropriate software installed, you can download article citation data to the citation manager of your choice. Simply select your manager software from the list below and click Download.