Stability of a one-dimensional morphoelastic model for post-burn contraction

Journal Article (2021)
Author(s)

Ginger Egberts (Universiteit Hasselt, TU Delft - Electrical Engineering, Mathematics and Computer Science)

Fred Vermolen (Universiteit Hasselt)

Paul van Zuijlen (Red Cross Hospital, Amsterdam UMC)

Research Group
Numerical Analysis
DOI related publication
https://doi.org/10.1007/s00285-021-01648-5 Final published version
More Info
expand_more
Publication Year
2021
Language
English
Related content
Research Group
Numerical Analysis
Journal title
Journal of Mathematical Biology
Issue number
3
Volume number
83
Article number
24
Pages (from-to)
1-24
Downloads counter
206
Collections
Institutional Repository
Reuse Rights

Other than for strictly personal use, it is not permitted to download, forward or distribute the text or part of it, without the consent of the author(s) and/or copyright holder(s), unless the work is under an open content license such as Creative Commons.

Abstract

To deal with permanent deformations and residual stresses, we consider a morphoelastic model for the scar formation as the result of wound healing after a skin trauma. Next to the mechanical components such as strain and displacements, the model accounts for biological constituents such as the concentration of signaling molecules, the cellular densities of fibroblasts and myofibroblasts, and the density of collagen. Here we present stability constraints for the one-dimensional counterpart of this morphoelastic model, for both the continuous and (semi-) discrete problem. We show that the truncation error between these eigenvalues associated with the continuous and semi-discrete problem is of order O(h2). Next we perform numerical validation to these constraints and provide a biological interpretation of the (in)stability. For the mechanical part of the model, the results show the components reach equilibria in a (non) monotonic way, depending on the value of the viscosity. The results show that the parameters of the chemical part of the model need to meet the stability constraint, depending on the decay rate of the signaling molecules, to avoid unrealistic results.