[论文解读] Spatially Heterogeneous Biofilm Simulations using an Immersed Boundary Method with Lagrangian Nodes Defined by Bacterial Locations

Jason F. Hammond, Elizabeth Stewart+3 more|arXiv (Cornell University)|Feb 15, 2013

Lattice Boltzmann Simulation Studies参考文献 3被引用 3

一句话总结

本文提出了一种新颖的浸入边界法，通过在三维实验测得的细菌细胞位置直接放置拉格朗日节点，对空间异质性生物膜进行建模，实现了流体-结构相互作用、黏弹性变形及流动下脱落的精确模拟。关键结果表明，生物膜密度和黏度显著影响脱落动力学，其中黏度增加可延迟脱落，而密度差异则极大改变结果——相比以往模型，该方法提供了更具生物学真实性的框架。

ABSTRACT

In this work we consider how surface-adherent bacterial biofilm communities respond in flowing systems. We simulate the fluid-structure interaction and separation process using the immersed boundary method. In these simulations we model and simulate different density and viscosity values of the biofilm than that of the surrounding fluid. The simulation also includes breakable springs connecting the bacteria in the biofilm. This allows the inclusion of erosion and detachment into the simulation. We use the incompressible Navier-Stokes (N-S) equations to describe the motion of the flowing fluid. We discretize the fluid equations using finite differences and use a geometric multigrid method to solve the resulting equations at each time step. The use of multigrid is necessary because of the dramatically different densities and viscosities between the biofilm and the surrounding fluid. We investigate and simulate the model in both two and three dimensions. Our method differs from previous attempts of using IBM for modeling biofilm/flow interactions in the following ways: the density and viscosity of the biofilm can differ from the surrounding fluid, and the Lagrangian node locations correspond to experimentally measured bacterial cell locations from 3D images taken of Staphylococcus epidermidis in a biofilm.

研究动机与目标

• 通过将实验测得的细菌细胞位置作为拉格朗日节点，开发一种更符合生物实际的生物膜-流体相互作用模型。
• 研究可变生物膜密度和黏度对流动系统中流体诱导变形与脱落动力学的影响。
• 利用几何多重网格法与有限差分法，提高高黏度、高密度生物膜模拟的数值精度与稳定性。
• 实现对不同流动条件下生物膜碎裂与脱落的预测性模拟，以供未来实验验证。

提出的方法

• 将浸入边界法扩展至使用三维成像数据获取的*Staphylococcus epidermidis*生物膜细菌细胞位置作为拉格朗日节点，而非网格细化节点。
• 在结构笛卡尔欧拉网格上采用有限差分法与几何多重网格求解器求解纳维-斯托克斯方程，以处理大黏度与密度对比。
• 通过细菌细胞间的线性弹簧模拟黏弹性行为，并修改应力本构方程以包含生物膜特异性黏度。
• 将狄拉克δ函数近似方法调整为基于细菌半径而非网格宽度进行缩放，提升力传递精度。
• 当键长超过其原始长度两倍时应用弹簧断裂准则，模拟脱落与侵蚀过程。
• 在二维与三维模拟中均进行了方法验证，包括复杂形状生物膜，重点考察收敛性与稳定性。

实验结果

研究问题

• RQ1 细菌分布的空间异质性如何影响生物膜对流体流动的响应？
• RQ2 改变生物膜密度与黏度对变形与脱落动力学有何影响？
• RQ3 不同数值策略——尤其是多重网格与网格自适应——如何提升高对比度流体-结构系统模拟的精度与稳定性？
• RQ4 能否通过实验获取的细菌位置提升生物膜模拟的真实感与预测能力？
• RQ5 生物膜黏弹性特性的变化如何影响脱落时间与碎片化模式？

主要发现

• 生物膜密度的微小变化显著改变脱落结果，表明生物膜与流体之间密度差异是以往模型常忽略的关键因素。
• 提高生物膜黏度可延迟脱落，证明黏弹性特性在生物膜抗流能力中起关键作用。
• 将实验测得的细菌位置作为拉格朗日节点，可实现更真实的生物膜形态与力学响应模拟。
• 几何多重网格法能有效求解高对比度流体-结构系统中的纳维-斯托克斯方程，但高生物膜黏度下收敛速度减慢。
• 对复杂三维形状生物膜的模拟显示出与实验观察一致的现实变形与碎片化模式，符合生物膜侵蚀现象。
• 该模型揭示，当前假设的均匀密度与线性弹性弹簧不足以描述真实生物膜，因为真实生物膜在断裂前表现出塑性变形与屈服行为。