A new approximation scheme is presented for the mathematical model of convection-diffusion and adsorption. The method is based on the relaxation method and the method of characteristics. We prove the convergence of the method and present some numerical experiments in 1D. The results can be applied to the model of contaminant transport in porous media with multi-site, equilibrium and non-equilibrium type of adsorption.