Numerical simulation of colloid-facilitated impurity transport in porous media

A numerical model is developed to calculate characteristics of impurity transport (including non-classical regimes) in double porous media in the presence of colloids. A governing equation set in the space-time representation is formulated basing on earlier proposed by one of the authors (L. M.) model. Analytical expressions for transport characteristics in limiting time intervals in the one-dimensional case are derived. The features of the effect of accelerating the impurity transport due to the presence of colloids are revealed. An interpolation form is proposed for the integral kernel entering in one of the transport equations. On a basis of this expression, an algorithm is developed for numerical modeling the impurity transport in double porous media in the presence of colloids. Trial numerical 1D calculations are performed basing on this algorithm. Good agreement was found between the numerical simulation results and the limiting analytical expressions.