A Multi-Dimensional Mathematical Model for Surface Exerted Point Forces in Elastic Media

Sabia Asghar,Qiyao Peng,Fred J. Vermolen

2024-10-14

Abstract:Some phenotypes of biological cells exert mechanical forces on their direct environment during their development and progression. In this paper the impact of cellular forces on the surrounding tissue is considered. Assuming the size of the cell to be much smaller than that of the computational domain, and assuming small displacements, linear elasticity (Hooke's Law) with point forces described by Dirac delta distributions is used in momentum balance equation. Due to the singular nature of the Dirac delta distribution, the solution does not lie in the classical $H^1$ finite element space for multi-dimensional domains. We analyze the $L^2$-convergence of forces in a superposition of line segments across the cell boundary to an integral representation of the forces on the cell boundary. It is proved that the $L^2$-convergence of the displacement field away from the cell boundary matches the quadratic order of convergence of the midpoint rule on the forces that are exerted on the curve or surface that describes the cell boundary.

Numerical Analysis

What problem does this paper attempt to address?

The paper attempts to address the problem of modeling and numerical analysis of the mechanical effects caused by point forces applied on the cell surface in an elastic medium. Specifically, the paper focuses on the mechanical interactions between cells and their surrounding environment (such as the extracellular matrix) and describes the impact of these forces on the surrounding tissue through a mathematical model. The main objectives of the study include:

1. **Modeling**: Establish a multi-dimensional mathematical model to describe the propagation and impact of point forces applied on the cell surface in an elastic medium. To simplify the problem, it is assumed that the cell size is much smaller than the size of the computational domain, and the displacements are small, thus allowing the use of linear elasticity theory (Hooke's law).

2. **Numerical Analysis**: Analyze the application of point forces described by the Dirac δ distribution in the momentum balance equation and explore the singularity of the solution caused by this distribution. Due to the singular nature of the Dirac δ distribution, classical finite element methods cannot directly handle such problems. Therefore, the paper proposes techniques to remove singularities to analyze the convergence of the solution.

3. **Convergence Analysis**: Prove the L2 convergence between the discrete representation (integral form approximation) and the continuous representation (integral form) of the point forces applied on the cell boundary. Specifically, it is shown that the convergence order of the displacement field away from the cell boundary matches the quadratic convergence order of the midpoint rule applied to the forces described on the curve or surface of the cell boundary.

Through these studies, the paper aims to better understand and quantify the mechanical interactions between cells and their surrounding environment, particularly in biological phenomena such as wound healing and cancer. This is of great significance for developing new therapeutic methods and predicting tumor growth and metastasis.