Numerical Methods to Compute Stresses and Displacements from Cellular Forces: Application to the Contraction of Tissue

by   Qiyao Peng, et al.

We consider a mathematical model for wound contraction, which is based on solving a momentum balance under the assumptions of isotropy, homogeneity, Hooke's Law, infinitesimal strain theory and point forces exerted by cells. However, point forces, described by Dirac Delta distributions lead to a singular solution, which in many cases may cause trouble to finite element methods due to a low degree of regularity. Hence, we consider several alternatives to address point forces, that is, whether to treat the region covered by the cells that exert forces as part of the computational domain or as 'holes' in the computational domain. The formalisms develop into the immersed boundary approach and the 'hole approach', respectively. Consistency between these approaches is verified in a theoretical setting, but also confirmed computationally. However, the 'hole approach' is much more expensive and complicated for its need of mesh adaptation in the case of migrating cells while it increases the numerical accuracy, which makes it hard to adapt to the multi-cell model. Therefore, for multiple cells, we consider the polygon that is used to approximate the boundary of cells that exert contractile forces. It is found that a low degree of polygons, in particular triangular or square shaped cell boundaries, already give acceptable results in engineering precision, so that it is suitable for the situation with a large amount of cells in the computational domain.


page 1

page 2

page 3

page 4


Point Forces in Elasticity Equation and Their Alternatives in Multi Dimensions

We consider several mathematical issues regarding models that simulate f...

Distributed model order reduction of a model for microtubule-based cell polarization using HAPOD

In this contribution we investigate in mathematical modeling and efficie...

Zyxin is all you need: machine learning adherent cell mechanics

Cellular form and function emerge from complex mechanochemical systems w...

Quality of approximating a mass-emitting object by a point source in a diffusion model

For the sake of computational efficiency and for theoretical purposes, i...

Quadrature for Implicitly-defined Finite Element Functions on Curvilinear Polygons

H^1-conforming Galerkin methods on polygonal meshes such as VEM, BEM-FEM...

What do cells regulate in soft tissues on short time scales?

Cells within living soft biological tissues seem to promote the maintena...

Please sign up or login with your details

Forgot password? Click here to reset