A two‐dimensional continuum model of biofilm growth incorporating fluid flow and shear stress based detachment

We present a two‐dimensional biofilm growth model in a continuum framework using an Eulerian description. A computational technique based on the eXtended Finite Element Method (XFEM) and the level set method is used to simulate the growth of the biofilm. The model considers fluid flow around the biofilm surface, the advection–diffusion and reaction of substrate, variable biomass volume fraction and erosion due to the interfacial shear stress at the biofilm–fluid interface. The key assumptions of the model and the governing equations of transport, biofilm kinetics and biofilm mechanics are presented. Our 2D biofilm growth results are in good agreement with those obtained by Picioreanu et al. (Biotechnol Bioeng 69(5):504–515, 2000). Detachment due to erosion is modeled using two continuous speed functions based on: (a) interfacial shear stress and (b) biofilm height. A relation between the two detachment models in the case of a 1D biofilm is established and simulated biofilm results with detachment in 2D are presented. The stress in the biofilm due to fluid flow is evaluated and higher stresses are observed close to the substratum where the biofilm is attached. Biotechnol. Bioeng. 2009;103: 92–104. © 2008 Wiley Periodicals, Inc.

[1]  A new mathematical model for chemotactic bacterial colony growth. , 2004, Water science and technology : a journal of the International Association on Water Pollution Research.

[2]  J. Wimpenny,et al.  A unifying hypothesis for the structure of microbial biofilms based on cellular automaton models , 1997 .

[3]  Paul Stoodley,et al.  A model of fluid–biofilm interaction using a Burger material law , 2007, Biotechnology and bioengineering.

[4]  I. Ozturk,et al.  The effect of anaerobic pre-treatment on the inert soluble COD of fermentation industry effluents , 1995 .

[5]  J J Heijnen,et al.  Effect of diffusive and convective substrate transport on biofilm structure formation: a two-dimensional modeling study. , 2000, Biotechnology and bioengineering.

[6]  D. Chopp,et al.  A mathematical model of quorum sensing in a growing bacterial biofilm , 2002, Journal of Industrial Microbiology and Biotechnology.

[7]  Hideki Harada,et al.  Adhesion strength of biofilm developed in an attached-growth reactor , 1994 .

[8]  J. J. Heijnen,et al.  Discrete-differential modelling of biofilm structure , 1999 .

[9]  B. Rittmann,et al.  A unified theory for extracellular polymeric substances, soluble microbial products, and active and inert biomass. , 2002, Water research.

[10]  Pedro Díez,et al.  Numerical modelling of tectonic plates subduction using X-FEM , 2007 .

[11]  P. Stewart,et al.  Biofilm accumulation model that predicts antibiotic resistance of Pseudomonas aeruginosa biofilms , 1994, Antimicrobial Agents and Chemotherapy.

[12]  Anders Heyden,et al.  Three‐dimensional biofilm model with individual cells and continuum EPS matrix , 2006, Biotechnology and bioengineering.

[13]  B. Peyton,et al.  A statistical analysis of the effect of substrate utilization and shear stress on the kinetics of biofilm detachment , 1993, Biotechnology and bioengineering.

[14]  S. Karamanos,et al.  Finite element modeling to expand the UMCCA model to describe biofilm mechanical behavior , 2005 .

[15]  David L. Chopp,et al.  Influence of the Hydrodynamic Environment on Quorum Sensing in Pseudomonas aeruginosa Biofilms , 2007, Journal of bacteriology.

[16]  P. M. Gresho,et al.  Incompressible flow and the finite element method. Volume 2: Incompressible flow and finite element , 1998 .

[17]  D. Chopp,et al.  The impact of quorum sensing and swarming motility on Pseudomonas aeruginosa biofilm formation is nutritionally conditional , 2006, Molecular microbiology.

[18]  J J Heijnen,et al.  Two-dimensional model of biofilm detachment caused by internal stress from liquid flow. , 2001, Biotechnology and bioengineering.

[19]  J. J. Heijnen,et al.  A three-dimensional numerical study on the correlation of spatial structure, hydrodynamic conditions, and mass transfer and conversion in biofilms , 2000 .

[20]  J. Keener,et al.  A mathematical model for quorum sensing in Pseudomonas aeruginosa , 2001, Bulletin of mathematical biology.

[21]  D. Chopp,et al.  A combined extended finite element and level set method for biofilm growth , 2008 .

[22]  Cory J. Rupp,et al.  Viscoelastic fluid description of bacterial biofilm material properties. , 2002, Biotechnology and bioengineering.

[23]  Benjamin L. Vaughan.  Application of the eXtended finite element method in mathematical biology , 2007 .

[24]  W Gujer,et al.  A multispecies biofilm model , 1986, Biotechnology and bioengineering.

[25]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[26]  I. Klapper,et al.  A Multidimensional Multispecies Continuum Model for Heterogeneous Biofilm Development , 2007 .

[27]  B. Rittman,et al.  The effect of shear stress on biofilm loss rate. , 1982, Biotechnology and bioengineering.

[28]  J J Heijnen,et al.  A theoretical study on the effect of surface roughness on mass transport and transformation in biofilms. , 2000, Biotechnology and bioengineering.

[29]  Wing Kam Liu,et al.  Nonlinear Finite Elements for Continua and Structures , 2000 .

[30]  H. Eberl,et al.  Exposure of biofilms to slow flow fields: the convective contribution to growth and disinfection. , 2008, Journal of theoretical biology.

[31]  D. Chopp,et al.  The dependence of quorum sensing on the depth of a growing biofilm , 2003, Bulletin of mathematical biology.

[32]  J. Kreft,et al.  Microbial motility involvement in biofilm structure formation--a 3D modelling study. , 2007, Water science and technology : a journal of the International Association on Water Pollution Research.

[33]  Bruce E. Rittmann,et al.  Simulation of multispecies biofilm development in three dimensions , 1999 .

[34]  P. Stewart,et al.  A model of biofilm detachment , 1993, Biotechnology and bioengineering.

[35]  Isaac Klapper,et al.  Finger Formation in Biofilm Layers , 2002, SIAM J. Appl. Math..

[36]  I. Klapper,et al.  Description of mechanical response including detachment using a novel particle model of biofilm/flow interaction. , 2007, Water science and technology : a journal of the International Association on Water Pollution Research.

[37]  Ted Belytschko,et al.  Arbitrary discontinuities in finite elements , 2001 .

[38]  A B Cunningham,et al.  Modeling biofilm accumulation and mass transport in a porous medium under high substrate loading , 1995, Biotechnology and bioengineering.

[39]  Bruce E Rittmann,et al.  Non-steady state modeling of extracellular polymeric substances, soluble microbial products, and active and inert biomass. , 2002, Water research.

[40]  David L. Chopp,et al.  The extended finite element method for boundary layer problems in biofilm growth , 2007 .

[41]  D. F. Parker,et al.  A New Deterministic Spatio-Temporal Continuum Model for Biofilm Development , 2001 .

[42]  G. Booth,et al.  BacSim, a simulator for individual-based modelling of bacterial colony growth. , 1998, Microbiology.

[43]  B. Hawkins.  A Framework for the , 2004 .

[44]  Mark C M van Loosdrecht,et al.  A general description of detachment for multidimensional modelling of biofilms. , 2005, Biotechnology and bioengineering.

[45]  David L. Chopp,et al.  A comparison of the extended finite element method with the immersed interface method for elliptic equations with discontinuous coefficients and singular sources , 2006 .

[46]  David L. Chopp,et al.  Some Improvements of the Fast Marching Method , 2001, SIAM J. Sci. Comput..

[47]  Qi Wang,et al.  PHASE-FIELD MODELS FOR BIOFILMS. I. THEORY AND 1-D SIMULATIONS ∗ , 2008 .

[48]  Tianyu Zhang,et al.  Phase-Field Models for Biofilms II. 2-D Numerical Simulations of Biofilm-Flow Interaction , 2008 .

[49]  Cristian Picioreanu,et al.  Constrained discounted Markov decision processes and Hamiltonian cycles , 1997, Proceedings of the 36th IEEE Conference on Decision and Control.

[50]  Mark C M van Loosdrecht,et al.  A framework for multidimensional modelling of activity and structure of multispecies biofilms. , 2005, Environmental microbiology.

[51]  B. Rittmann,et al.  Evaluating trends in biofilm density using the UMCCA model. , 2004, Water research.

[52]  Mark C.M. van Loosdrecht,et al.  A more unifying hypothesis for biofilm structures , 1997 .

[53]  S W Hermanowicz.  A simple 2D biofilm model yields a variety of morphological features. , 2001, Mathematical biosciences.