Open AccessThe Finite-Difference Model of Fully Saturated Groundwater Contaminant Transport
Author(s) -
N. P. Purnaditya,
Herr Soeryantono,
Dwinanti Rika Marthanty,
Jessica Sjah
Publication year - 2018
Publication title -
international journal of engineering and technology
Language(s) - English
Resource type - Journals
ISSN - 2227-524X
DOI - 10.14419/ijet.v7i4.35.23073
Subject(s) - groundwater , aquifer , dispersion (optics) , mechanics , groundwater flow equation , advection , flux (metallurgy) , groundwater flow , groundwater model , diffusion , finite element method , diffusion equation , finite difference , environmental science , geotechnical engineering , geology , materials science , mathematics , thermodynamics , physics , engineering , mathematical analysis , metric (unit) , operations management , metallurgy , optics
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.