Numerical solutions for the advection-diffusion-reaction equation of Bacteroidales with an arbitrary, solenoidal, velocity field, with A = 1. Contours indicate concentration of Bacteroidales; streamlines of the steady-state velocity field have been also indicated. (a) Power-law scheme. (b) Upwind scheme.