Stability analysis of a continuum-based constrained mixture model for vascular growth and remodeling

A stabilizing criterion is derived for equations governing vascular growth and remodeling. We start from the integral state equations of the continuum-based constrained mixture theory of vascular growth and remodeling and obtain a system of time-delayed differential equations describing vascular growth. By employing an exponential form of the constituent survival function, the delayed differential equations can be reduced to a nonlinear ODE system. We demonstrate the degeneracy of the linearized system about the homeostatic state, which is a fundamental cause of the neutral stability observations reported in prior studies. Due to this degeneracy, stability conclusions for the original nonlinear system cannot be directly inferred. To resolve this problem, a sub-system is constructed by recognizing a linear relation between two states. Subsequently, Lyapunov's indirect method is used to connect stability properties between the linearized system and the original nonlinear system, to rigorously establish the neutral stability properties of the original system. In particular, this analysis leads to a stability criterion for vascular expansion in terms of growth and remodeling kinetic parameters, geometric quantities and material properties. Numerical simulations were conducted to evaluate the theoretical stability criterion under broader conditions, as well as study the influence of key parameters and physical factors on growth properties. The theoretical results are also compared with prior numerical and experimental findings in the literature.