Numerical Analysis of the Immersed Boundary Method for Cell-Based Simulation

Abstract

Mathematical modelling provides a useful framework within which to investigate the organization of biological tissues. With advances in experimental biology leading to increasingly detailed descriptions of cellular behaviour, models that consider cells as individual objects are becoming a common tool to study how processes at the single-cell level affect collective dynamics and determine tissue size, shape and function. However, there often remains no comprehensive account of these models, their method of solution, computational implementation or analysis of parameter scaling, hindering our ability to utilise and accurately compare different models. Here we present an effcient, open-source implementation of the immersed boundary method (IBM), tailored to simulate the dynamics of cell populations. This approach considers the dynamics of elastic membranes, representing cell boundaries, immersed in a viscous Newtonian fluid. The IBM enables complex and emergent cell shape dynamics, spatially heterogeneous cell properties and precise control of growth mechanisms. We solve the model numerically using an established algorithm, based on the fast Fourier transform, providing full details of all technical aspects of our implementation. The implementation is undertaken within Chaste, an open-source C++ library that allows one to easily change constitutive assumptions. Our implementation scales linearly with time step, and subquadratically with mesh spacing and immersed boundary node spacing. We identify the relationship between the immersed boundary node spacing and fluid mesh spacing required to ensure fluid volume conservation within immersed boundaries, and the scaling of cell membrane stiffness and cell-cell interaction strength required when refining the immersed boundary discretization. This study provides a recipe allowing consistent parametrisation of IBM models.