A numerical algorithm for computational modelling of coupled advection-diffusion-reaction systems