A Multiple-Relaxation-Time Lattice-Boltzmann Model for Bacterial Chemotaxis: Effects of Initial Concentration, Diffusion, and Hydrodynamic Dispersion on Traveling Bacterial Bands

Bacterial chemotaxis can enhance the bioremediation of contaminants in aqueous and subsurface environments if the contaminant is a chemoattractant that the bacteria degrade. The process can be promoted by traveling bands of chemotactic bacteria that form due to metabolism-generated gradients in chemoattractant concentration. We developed a multiple-relaxation-time (MRT) lattice-Boltzmann method (LBM) to model chemotaxis, because LBMs are well suited to model reactive transport in the complex geometries that are typical for subsurface porous media. This MRT-LBM can attain a better numerical stability than its corresponding single-relaxation-time LBM. We performed simulations to investigate the effects of substrate diffusion, initial bacterial concentration, and hydrodynamic dispersion on the formation, shape, and propagation of bacterial bands. Band formation requires a sufficiently high initial number of bacteria and a small substrate diffusion coefficient. Uniform flow does not affect the bands while shear flow does. Bacterial bands can move both upstream and downstream when the flow velocity is small. However, the bands disappear once the velocity becomes too large due to hydrodynamic dispersion. Generally bands can only be observed if the dimensionless ratio between the chemotactic sensitivity coefficient and the effective diffusion coefficient of the bacteria exceeds a critical value, that is, when the biased movement due to chemotaxis overcomes the diffusion-like movement due to the random motility and hydrodynamic dispersion.