Nondimensionalization, Variational Inequalities and Uncertainty Quantification in a Morphoelastic Model for Burn Injury Healing
Abstract
Mathematical modelling of burn injury healing involves coupling biological processes with mechanical deformation, resulting in computationally stiff, nonlinear partial differential equations that severely limit the scalability of existing numerical solutions. We present a computational framework addressing these challenges. By nondimensionalizing the model and performing scaling analysis, we show that chemotaxis and inertia are negligible, permitting a quasi-static momentum approach. To enforce non-negativity constraints, we implement a Variational Inequality formulation using the Finite Element Method in FEniCSx, and verify it via numerical convergence studies and mass balances. We identify high-order parameter interactions, specifically equilibrium collagen density and signalling dynamics, as the primary drivers of solver stability. A sensitivity analysis of nearly four million simulations reveals that isolated parameter variations predominantly govern the maximum wound contraction. In contrast, higher-order interactions determine the timing of maximum contraction and the permanent contracture. However, no significant pairwise parameter interactions were found, suggesting that the model's behaviour is governed by higher-order interactions, highlighting the inherent complexity of the wound-healing process.
Related
- Enhanced Uncertainty Quantification Variational Autoencoders for the Solution of Bayesian Inverse Problems
- Sharp logarithmic Sobolev inequalities on free groups by variational comparison
- Lp-Integral Affine Surface Area-Variational Capacity Inequalities
- Some Properties of Mordukhovich Derivatives with Applications to Locally Variational Inequalities in Banach Spaces
- Solving Quasi-Variational Inequalities Using the Progressive Decoupling of Linkages
BibTeX
@misc{hertog2026nondimensionalization,
title = {{Nondimensionalization, Variational Inequalities and Uncertainty Quantification in a Morphoelastic Model for Burn Injury Healing}},
author = {Wilbert den Hertog and Fred Vermolen},
howpublished = {arXiv},
year = {2026},
month = oct,
doi = {10.48550/arxiv.2610.23177},
eprint = {2610.23177},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2610.23177},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Nondimensionalization, Variational Inequalities and Uncertainty Quantification in a Morphoelastic Model for Burn Injury Healing
AU - Wilbert den Hertog
AU - Fred Vermolen
PY - 2026
DA - 2026/10/16
PB - arXiv
DO - 10.48550/arxiv.2610.23177
UR - https://arxiv.org/abs/2610.23177
AB - Mathematical modelling of burn injury healing involves coupling biological processes with mechanical deformation, resulting in computationally stiff, nonlinear partial differential equations that severely limit the scalability of existing numerical solutions. We present a computational framework addressing these challenges. By nondimensionalizing the model and performing scaling analysis, we show that chemotaxis and inertia are negligible, permitting a quasi-static momentum approach. To enforce non-negativity constraints, we implement a Variational Inequality formulation using the Finite Element Method in FEniCSx, and verify it via numerical convergence studies and mass balances. We identify high-order parameter interactions, specifically equilibrium collagen density and signalling dynamics, as the primary drivers of solver stability. A sensitivity analysis of nearly four million simulations reveals that isolated parameter variations predominantly govern the maximum wound contraction. In contrast, higher-order interactions determine the timing of maximum contraction and the permanent contracture. However, no significant pairwise parameter interactions were found, suggesting that the model's behaviour is governed by higher-order interactions, highlighting the inherent complexity of the wound-healing process.
ER -