Abstract:
Solution of transport equation by conventional numerical methods such as finite difference and finite element analysis exhibit either excessive diffusion and/or oscillations near the concentration front. This is due to the presence of advection term in the solute transport equation. In the present work, a computationally simple, numerical algorithm is developed to solve the solute transport equation in groundwater. The governing equation is solved using finite differences employing the modified Picard iteration scheme to determine the temporal derivative of the solute concentration. The total solute concentration is expanded in a Taylor series with respect to the solution concentration to linearize the transport equation, which then solved with conventional finite difference method. The algorithm avoids mass-balance errors and is numerically stable. The numerical solution is compared with analytical solution. The model is further being used for virus transport in ground water.
