Multiscale Modeling of Uranium Bioreduction in Porous Media by One-Dimensional Biofilms

We formulate a multiscale mathematical model that describes the bioreduction of uranium in porous media. On the mesoscale we describe the bioreduction of uranium [VI] to uranium [IV] using a multispecies one-dimensional biofilm model with suspended bacteria and thermodynamic growth inhibition. We upscale the mesoscopic (colony scale) model to the macroscale (reactor scale) and investigate the behavior of substrate utilization and production, attachment and detachment processes, and thermodynamic effects not usually considered in biofilm growth models. Simulation results of the reactor model indicate that thermodynamic inhibition quantitatively alters the dynamics of the model and neglecting thermodynamic effects may over- or underestimate chemical concentrations in the system. Furthermore, we numerically investigate uncertainties related to the specific choice of attachment and detachment rate coefficients and find that while increasing the attachment rate coefficient or decreasing the detachment rate coefficient leads to thicker biofilms, performance of the reactor remains largely unaffected.