NUMERICAL SOLUTION OF ADVECTION-DIFFUSION EQUATION

Abstract

Accurate solution of the Advection-Dispersion equation by numerical methods is difficult due to the presence of artificial diffusion induced while approximating the advective term. Many numerical schemes have been developed to solve Advection-Dispersion equation. In the present study, the performance of MacCormack scheme has been investigated while solving Advection-Dispersion equation for wide range of Courant and Peclet numbers. MacCormack scheme is chosen due to its simplicity. The advective part is solved by MacCormack scheme and. the dispersion part by an implicit. finite difference method. Three forms of MacCormack scheme have been used. Simulations are carried out for the constant concentration source at the inlet. It has been found that treating the change in concentration due to advection on a sink/ source term in the dispersive part gives accurate results.