Numerical simulation for solute transport in fractal porous media

Abstract

A modified Fokker-Planck equation with continuous source for solute transport in fractal porous media is considered. The dispersion term of the governing equation uses a fractional-order derivative and the diffusion coefficient can be time and scale dependent. In this paper, numerical solution of the modified Fokker-Planck equation is proposed. The effects of different fractional orders and fractional power functions of time and distance are numerically investigated. The results show that motions with a heavy tailed marginal distribution can be modelled by equations that use fractional-order derivatives and/or time and scale dependent dispersivity.