The Finite-Difference Model of Fully Saturated Groundwater Contaminant Transport

  • Authors

    • NgakanPutu Purnaditya
    • Herr Soeryantono
    • Dwinanti Rika Marthanty
    • Jessica Sjah
  • Alternating Direction Implicit, Finite-Difference Method, Groundwater Contaminant Transport, Numerical Modelling
  • Groundwater quality is one of water resource problem. This problem is driven by contaminant transport phenomena and can be described as a mathematical model. Contaminant transports equation usually is composed by advection and dispersion flux. In the porous medium or aquifer, contaminant meet tortuosity effect, therefore hydrodynamic dispersion must be considered as the development of dispersion flux. This paper explains the mathematical model of groundwater contaminant on fully saturated condition. It starts from governing equation of contaminant transport. Advection flux is based on groundwater velocity. In steady state condition, groundwater velocity can be determined as certain value. In another hand, in the transient condition, groundwater velocity must be determined based on the solution of groundwater flow model. Dispersion flux is calculated through first Fick’s law and this component is distinguished into two parts follow mechanical dispersion and molecular diffusion. Mechanical dispersion affected by groundwater velocity and dispersivity. Contaminant transport equation is solved numerically using Finite Difference Method (FDM). This final model is validated theoretically and then this model is simulated into transient condition. The result of the simulation is described and explained graphically. Based on this research, the result of FDM model has similar physical behavior to FEM model from CTRAN example.

  • References

    1. [1] Bear, J. (1979). Hydraulics of Groundwater. New York: Dover Publication, inc.

      [2] Bear, J., & Cheng, A. (2010). Modeling Groundwater Flow adn Contaminant Transport. New York: Springer Science+Bussiness Media.

      [3] Chapra, S., & Canale, R. (2015). Numerical Methods for Engineers (7th ed.). New York: McGraw-Hill Education.

      [4] Corapcioglu, M. Y., & Wang, S. (1999). Dual-Porosity Groundwater Contaminant Transport in the Presence of Colloids. Water Resources Research, 35(11), 3261–3273.

      [5] Fischer, H., List, J., Koh, R., Imberger, J., & Brooks, N. (1979). Mixing in Inland and Coastal Waters. London: Academic Press, inc.

      [6] Freeze, A., & Cherry, J. (1979). Groundwater. New Jersey: Prentice-Hall, inc.

      [7] Geo-Slope. (n.d.). CTRAN Example File: Henry Saltwater Intrusion Problem. Alberta. Retrieved from

      [8] Guo, W., & Langevin, C. (2002). User’s Guide to Seawat: A Computer Program for Simulation of Three-Dimensional Variable-Density Ground-Water Flow. Florida: U.S Geological Survey.

      [9] Harbaugh, A. (2005). Modflow-2005, The U.S. Geological Survey Modular Ground-Water Model-the Ground-Water Flow Process. Virginia: U.S Geological Survey.

      [10] gboekwe, M., & Achi, N. . (2011). Finite Difference Method on Modelling Groundwater Flow. Journal of Water Resources and Protection, 3, 192–198.

      [11] Kumar, C. (2012). Numerical Model Parameters for Simulation of Contaminant Trasnport. International Journal of Scientific Research and Reviews, 1(3), 45–63.

      [12] Lalehzari, R., Tabatabaei, S. ., & Kholghi, M. (2013). Simulation of Nitrate Transport and Wastewater Seepage in Groundwater Flow System. International Journal of Environmental Science and Technology, 10, 1367–1376.

      [13] Notodarmojo, S. (2005). Pencemaran Tanah dan Air Tanah. Bandung: Penerbit ITB.

      [14] Okiro, J. ., Manyonge, A. ., Ongati, N. ., Shichikha, J. ., & Kimaiyo, J. . (2013). On the Solution of Confined Aquifer Flow Equations: Finite Difference Approximations. Applied Mathematical Science, 7(58), 2885–2896.

      [15] Patil, S. ., & Chore, H. . (2014). Contaminant Transport Through Porous Media: An Overview of Experimental and Numerical Studies. Advances in Environmental Research, 3(1), 45–69.

      [16] Rao, G. T., Rao, G., Surinaidu, L., Mahesh, J., & Padalu, G. (2011). Application of Numerical for Groundwater Flow and Contaminant Transport Analysis in the Basaltic Terrain, Bagalkot, India. Arabian Journal of Geosciences.

      [17] Rubin, H., & Atkinson, J. (2001). Environmental Fluid Mechanics. New York: Marcel Dekker, inc.

      [18] Sharma, P. ., Sawant, V. ., Shukla, S. ., & Khan, Z. (2014). Experimental and Numerical Simulation of Contaminant Transport Through Layered Soil. International Journal of Geotechnical Engineering, 8(4), 345–351.

      [19] Triatmodjo, B. (2002). Metode Numerik. Yogyakarta: Beta Offset.

      [20] Wang, H., & Anderson, M. (1982). Introduction to Groundwater Modeling: Finite Difference and Finite Element Methods. California: Academic Press, inc.

  • Downloads

  • How to Cite

    Purnaditya, N., Soeryantono, H., Marthanty, D. R., & Sjah, J. (2018). The Finite-Difference Model of Fully Saturated Groundwater Contaminant Transport. International Journal of Engineering & Technology, 7(4.35), 629-634.