Biofilms are typically heterogeneous in morphology, structure, and composition, resulting in non-uniform mechanical properties. The distribution of mechanical properties, in turn, determines the biofilm mechanical behavior, such as deformation and detachment. Most past studies neglected heterogeneity of biofilms. In this study, an image-based modeling approach was developed to transform two-dimensional optical coherence tomography biofilm images to a pixel-scale non-Newtonian viscosity map of the biofilm. The spatial distribution of non-Newtonian viscosity was applied in an established Oldroyd-B constitutive model and implemented using the phase-field continuum approach for the deformation and stress analysis. The heterogeneous model was able to predict deformations and stresses more accurately than a homogenous one. This is the first time, to the best of our knowledge, that an image-based approach is used to map the mechanical heterogeneity of biofilms for computational studies. It provides an efficient method to characterize biofilm mechanical behavior.

The mechanical properties of biofilms can be used to predict biofilm deformation, for example under fluid flow. We used magnetic tweezers to spatially map the compliance of Pseudomonas aeruginosa biofilms at the micron scale, then used modeling to assess its effects on biofilm deformation. Biofilms were grown in capillary flow cells with Reynolds numbers (Re) ranging 0.28 to 13.9, bulk dissolved oxygen (DO) concentrations from 1 mg/L to 8 mg/L, and bulk calcium ion (Ca2+) concentrations of 0 and 100 mg CaCl2/L. Higher Re numbers resulted in more uniform biofilm morphologies. The biofilm was stiffer at the center of the flow cell than near the walls. Lower bulk DO led to more stratified biofilms. Higher Ca2+ led to increased stiffness and more uniform mechanical properties. Using the experimental mechanical properties, fluid-structure interaction models predicted up to 64% greater deformations for heterogeneous biofilms, compared to a homogeneous biofilms with the same average properties. However, the error depended on the biofilm morphology and flow regime. Our results show significant spatial mechanical variability exists at the micron scale, and that this variability can potentially affect biofilm deformation. The average mechanical properties, provided in many studies, should be used with caution when predicting biofilm deformation.

Biofilms commonly develop in flowing aqueous environments, where the flow causes the biofilm to deform. Because biofilm deformation affects the flow regime, and because biofilms behave as complex heterogeneous viscoelastic materials, few models are able to predict biofilm deformation. In this study, a phase field continuum model coupled with the Oldroyd-B constitutive equation was developed and used to simulate biofilm deformation. The accuracy of the model was evaluated using two types of biofilms: a synthetic biofilm, made from alginate mixed with bacterial cells, and a Pseudomonas aeruginosa biofilm. Shear rheometry was used to experimentally determine the mechanical parameters for each biofilm, as inputs for the model. Biofilm deformation under fluid flow was monitored experimentally using optical coherence tomography. The fit between the experimental and modeling geometries after fluid-driven deformation was very good, with relative errors of 12.8% for synthetic biofilm and 22.2% for homogenized P. aeruginosa biofilm. This is the first demonstration of the effectiveness of a viscoelastic phase field biofilm model. This model provides an important tool for predicting biofilm viscoelastic deformation. It also can benefit the design and control of biofilms in engineering systems.