Error analysis of finite difference/collocation method for the nonlinear coupled parabolic free boundary problem modeling plaque growth in the artery
Pith reviewed 2026-05-25 00:18 UTC · model grok-4.3
The pith
Front-fixing transformation followed by nonclassical finite difference and collocation methods produces stable convergent solutions for the nonlinear coupled plaque growth model.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
After the front-fixing transformation and change of variables, the nonclassical finite difference and collocation methods applied to the resulting system of equations are stable and converge to the solution of the original free-boundary model.
What carries the argument
Front-fixing transformation that immobilizes the free boundary and converts the mixed boundary condition into a Neumann condition, followed by nonclassical finite difference and collocation discretization of the fixed-domain system.
If this is right
- The discrete solutions can be used to compute the time evolution of plaque thickness and the concentrations of the biological species inside the artery wall.
- Error bounds between the numerical approximation and the true solution follow directly from the convergence proof.
- The same discretization framework applies to the coupled system without further modification once the transformation has been performed.
Where Pith is reading between the lines
- The same front-fixing plus collocation strategy could be tested on other free-boundary models that couple reaction-diffusion equations to an ordinary differential equation for interface motion.
- If the convergence rates remain high under realistic parameter ranges, the scheme could support inverse problems that estimate growth coefficients from medical imaging sequences.
Load-bearing premise
The front-fixing transformation and the subsequent change of variables that converts the mixed boundary condition into a Neumann condition remain valid for the full nonlinear coupled system.
What would settle it
A sequence of numerical solutions on successively refined meshes that fails to approach a known exact solution, or that exhibits instability for parameter values inside the model's physical range, would show that the stability or convergence claim does not hold.
Figures
read the original abstract
The main target of this paper is to present a new and efficient method to solve a nonlinear free boundary mathematical model of atherosclerosis. This model consists of three parabolics, one elliptic and one ordinary differential equations that are coupled together and describe the growth of a plaque in the artery. We start our discussion by using the front fixing method to fix the free domain and simplify the model by changing the mix boundary condition to a Neumann one by applying suitable changes of variables. Then, after employing a nonclassical finite difference and the collocation method on this model, we prove the stability and convergence of methods. Finally, some numerical results are considered to show the efficiency of the method.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops a numerical method for a nonlinear free boundary problem modeling plaque growth, consisting of three coupled parabolic PDEs, one elliptic PDE, and one ODE. It applies a front-fixing transformation to fix the domain and a change of variables to convert mixed boundary conditions to Neumann type, then employs nonclassical finite difference and collocation methods, proves stability and convergence of the methods, and presents numerical results to demonstrate efficiency.
Significance. If the transformation is shown to preserve the structure of the full nonlinear coupled system and the stability/convergence proofs are rigorous, the work provides a useful contribution to numerical analysis of free-boundary problems in mathematical biology, with potential relevance to biomedical modeling of atherosclerosis.
major comments (1)
- [Transformation and change of variables (preceding the method application)] The front-fixing transformation and change of variables (described prior to the numerical discretization) are presented at a high level. It is not shown explicitly that the elliptic equation and the ODE remain consistent with the original nonlinear couplings, or that no additional singularities or coupling terms are introduced that would affect the subsequent finite-difference/collocation analysis. This verification is load-bearing for the stability and convergence claims.
minor comments (1)
- [Proof sections] The abstract and introduction assert proofs of stability and convergence, but the manuscript would benefit from clearer cross-references between the transformed system equations and the specific estimates used in the proofs.
Simulated Author's Rebuttal
We thank the referee for their careful reading of the manuscript and for highlighting this important point regarding the transformation step. We address the concern below and are prepared to revise the manuscript accordingly.
read point-by-point responses
-
Referee: [Transformation and change of variables (preceding the method application)] The front-fixing transformation and change of variables (described prior to the numerical discretization) are presented at a high level. It is not shown explicitly that the elliptic equation and the ODE remain consistent with the original nonlinear couplings, or that no additional singularities or coupling terms are introduced that would affect the subsequent finite-difference/collocation analysis. This verification is load-bearing for the stability and convergence claims.
Authors: We agree that the presentation of the transformed system can be strengthened by greater explicitness. In the revised manuscript we will add a dedicated subsection that derives the fully transformed equations for all five components of the model (the three parabolic PDEs, the elliptic PDE, and the ODE). This derivation will verify term-by-term that the original nonlinear couplings are preserved, that the change of variables maps the mixed boundary conditions to Neumann conditions without introducing singularities, and that the resulting system remains suitable for the subsequent stability and convergence analysis. The proofs will then reference these explicit transformed equations. revision: yes
Circularity Check
No significant circularity; derivation self-contained on transformed system
full rationale
The paper describes applying the front-fixing transformation and variable change to convert the free-boundary model into a fixed-domain system with Neumann conditions, followed by nonclassical finite differences and collocation, then proving stability and convergence on that transformed system. No quoted equations, fitted parameters, or self-citations reduce the stability/convergence claims to tautologies or inputs by construction. The transformation is presented as a modeling step whose validity is presupposed for the subsequent analysis; this is an assumption, not a circular reduction within the derivation chain itself. The numerical analysis is therefore independent of the original free-boundary formulation once the transformed equations are accepted.
Axiom & Free-Parameter Ledger
Reference graph
Works this paper leans on
-
[1]
Inflammation in atherosclerotic plaque formation [video file] youtube; [cited 2013 jan 19]. retrieved from:. https://www.youtube.com/watch?v=Na6-kP9VYCU
work page 2013
-
[2]
Math- ematical modelling of tumour invasion and metastasis
Alexander RA Anderson, Mark AJ Chaplain, E Luke Newman, Robert JC Steele, and Alastair M Thompson. Math- ematical modelling of tumour invasion and metastasis. Computational and mathematical methods in medicine , 2(2):129–154, 2000
work page 2000
-
[3]
Role of free radicals in cancer and atherosclerosis
DD Bankson, M Kestin, and N Rifai. Role of free radicals in cancer and atherosclerosis. Clinics in laboratory medicine, 13(2):463–480, 1993
work page 1993
-
[4]
Mechanisms of plaque formation and rupture
Jacob Fog Bentzon, Fumiyuki Otsuka, Renu Virmani, and Erling Falk. Mechanisms of plaque formation and rupture. Circulation research, 114(12):1852–1866, 2014. Title Suppressed Due to Excessive Length 33
work page 2014
-
[5]
The early history of partial differential equations and of partial differentiation and integration
Florian Cajori. The early history of partial differential equations and of partial differentiation and integration. The American Mathematical Monthly, 35(9):459–467, 1928
work page 1928
-
[6]
Mathematical modelling of the atheroscle- rotic plaque formation
Vincent Calvez, Abderrhaman Ebde, Nicolas Meunier, and Annie Raoult. Mathematical modelling of the atheroscle- rotic plaque formation. In ESAIM: Proceedings, volume 28, pages 1–12. EDP Sciences, 2009
work page 2009
-
[7]
Mathematical and numerical modeling of early atherosclerotic lesions
Vincent Calvez, Jean Gabriel Houot, Nicolas Meunier, Annie Raoult, and Gabriela Rusnakova. Mathematical and numerical modeling of early atherosclerotic lesions. In ESAIM: Proceedings, volume 30, pages 1–14. EDP Sciences, 2010
work page 2010
-
[8]
Claudio Canuto, M Youssuff Hussaini, Alfio Quarteroni, and Thomas A Zang. Spectral methods. Springer, 2006
work page 2006
-
[9]
A mathematical analysis of a model for tumour angiogenesis
Mark AJ Chaplain, Susan M Giles, BD Sleeman, and Richard J Jarvis. A mathematical analysis of a model for tumour angiogenesis. Journal of mathematical biology , 33(7):744–770, 1995
work page 1995
-
[10]
A mathematical model of cell-to-cell spread of hiv-1 that includes a time delay
Rebecca V Culshaw, Shigui Ruan, and Glenn Webb. A mathematical model of cell-to-cell spread of hiv-1 that includes a time delay. Journal of mathematical biology , 46(5):425–444, 2003
work page 2003
-
[11]
Zhong-Shan Deng and Jing Liu. Mathematical modeling of temperature mapping over skin surface and its imple- mentation in thermal disease diagnostics. Computers in biology and medicine , 34(6):495–521, 2004
work page 2004
-
[12]
EH Doha, WM Abd-Elhameed, and YH Youssri. Fully legendre spectral galerkin algorithm for solving linear one- dimensional telegraph type equation. International Journal of Computational Methods , page 1850118, 2018
work page 2018
-
[13]
Numerical simulation of growth of an atherosclerotic lesion with a moving boundary
Jens P Eberhard and Peter Frolkovic. Numerical simulation of growth of an atherosclerotic lesion with a moving boundary. 2006
work page 2006
-
[14]
Sakine Esmaili and Mohammad Reza Eslahchi. Optimal control for a parabolic–hyperbolic free boundary problem modeling the growth of tumor with drug application. Journal of Optimization Theory and Applications , 173(3):1013– 1041, 2017
work page 2017
-
[15]
Sakine Esmaili and MR Eslahchi. Application of fixed point-collocation method for solving an optimal control problem of a parabolic–hyperbolic free boundary problem modeling the growth of tumor with drug application. Computers & Mathematics with Applications , 75(7):2193–2216, 2018
work page 2018
-
[16]
Wound healing angiogenesis: the clinical implications of a simple mathematical model
Jennifer A Flegg, Helen M Byrne, Mark B Flegg, and DL Sean McElwain. Wound healing angiogenesis: the clinical implications of a simple mathematical model. Journal of theoretical biology, 300:309–316, 2012
work page 2012
-
[17]
Modelling the early growth of ductal carcinoma in situ of the breast
SJ Franks, HM Byrne, JR King, JCE Underwood, and CE Lewis. Modelling the early growth of ductal carcinoma in situ of the breast. Journal of mathematical biology , 47(5):424–452, 2003
work page 2003
-
[18]
Avner Friedman and Wenrui Hao. A mathematical model of atherosclerosis with reverse cholesterol transport and associated risk factors. Bulletin of mathematical biology , 77(5):758–781, 2015
work page 2015
-
[19]
A free boundary problem for steady small plaques in the artery and their stability
Avner Friedman, Wenrui Hao, and Bei Hu. A free boundary problem for steady small plaques in the artery and their stability. Journal of Differential Equations , 259(4):1227–1255, 2015
work page 2015
-
[20]
Walter Gautschi. Numerical analysis. Springer Science & Business Media, 1997
work page 1997
-
[21]
The role of the monocyte in atherogenesis: I
Ross G Gerrity. The role of the monocyte in atherogenesis: I. transition of blood-borne monocytes into foam cells in fatty lesions. The American journal of pathology , 103(2):181, 1981. 34 F. Nasresfahani, M. R. Eslahchi
work page 1981
-
[22]
The ldl-hdl profile determines the risk of atherosclerosis: a mathematical model
Wenrui Hao and Avner Friedman. The ldl-hdl profile determines the risk of atherosclerosis: a mathematical model. PloS one, 9(3):e90497, 2014
work page 2014
-
[23]
Mendelian proportions in a mixed population
Godfrey Harold Hardy et al. Mendelian proportions in a mixed population. Classic papers in genetics. Prentice-Hall, Inc.: Englewood Cliffs, NJ , pages 60–62, 1908
work page 1908
-
[24]
Exercitatio anatomica de motu cordis et sanguinis in animalibus , volume 2
William Harvey. Exercitatio anatomica de motu cordis et sanguinis in animalibus , volume 2. 1737
-
[25]
A mathematical model of tumour angiogenesis incorporating cellular traction and viscoelastic effects
MJ Holmes and BD Sleeman. A mathematical model of tumour angiogenesis incorporating cellular traction and viscoelastic effects. Journal of theoretical biology, 202(2):95–112, 2000
work page 2000
-
[26]
Spectral methods for substantial fractional differential equations
Can Huang, Zhimin Zhang, and Qingshuo Song. Spectral methods for substantial fractional differential equations. Journal of Scientific Computing , 74(3):1554–1574, 2018
work page 2018
-
[27]
A mathematical model for atherosclerotic plaque formation and arterial wall remodelling
Md Hamidul Islam and Peter Johnston. A mathematical model for atherosclerotic plaque formation and arterial wall remodelling. ANZIAM Journal, 57:320–345, 2016
work page 2016
-
[28]
A mathematical analysis of physiological and morphological aspects of wound closure
E Javierre, FJ Vermolen, C Vuik, and S Van der Zwaag. A mathematical analysis of physiological and morphological aspects of wound closure. Journal of mathematical biology , 59(5):605–630, 2009
work page 2009
-
[29]
Generalized bessel functions: Theory and their appli- cations
Hassan Khosravian-Arab, Mehdi Dehghan, and MR Eslahchi. Generalized bessel functions: Theory and their appli- cations. Mathematical Methods in the Applied Sciences , 40(18):6389–6410, 2017
work page 2017
-
[30]
Mark P Little, Anna Gola, and Ioanna Tzoulaki. A model of cardiovascular disease giving a plausible mechanism for the effect of fractionated low-dose ionizing radiation exposure. PLoS computational biology, 5(10):e1000539, 2009
work page 2009
-
[31]
Mutations of bacteria from virus sensitivity to virus resistance
Salvador E Luria and Max Delbr¨ uck. Mutations of bacteria from virus sensitivity to virus resistance. Genetics, 28(6):491, 1943
work page 1943
-
[32]
A model of partial differential equations for hiv propagation in lymph nodes
EBS Marinho, Flora S Bacelar, and RFS Andrade. A model of partial differential equations for hiv propagation in lymph nodes. Physica A: Statistical Mechanics and its Applications , 391(1-2):132–141, 2012
work page 2012
-
[33]
Mark M Meerschaert. Mathematical modeling. Academic press, 2013
work page 2013
-
[34]
Methodus fluxionum et seriarum infinitarum
Isaac Newton. Methodus fluxionum et seriarum infinitarum. Opuscula mathematica, philosophica et philologica , 1, 1774
-
[35]
International statistical classification of diseases and related health problems , volume 1
World Health Organization. International statistical classification of diseases and related health problems , volume 1. World Health Organization, 2004
work page 2004
-
[36]
An ode model of early stages of atherosclerosis: mechanisms of the inflammatory response
Anna Ougrinovskaia, Rosemary S Thompson, and Mary R Myerscough. An ode model of early stages of atherosclerosis: mechanisms of the inflammatory response. Bulletin of mathematical biology , 72(6):1534–1561, 2010
work page 2010
-
[37]
Mathematical modelling of macrophage dynamics in tumours
R Owen and Jonathan A Sherratt. Mathematical modelling of macrophage dynamics in tumours. Mathematical Models and Methods in Applied Sciences , 9(04):513–539, 1999
work page 1999
-
[38]
A mathematical system for human implantable wound model studies
Salomonsky Paul-Michael. A mathematical system for human implantable wound model studies. 2013
work page 2013
-
[39]
Mehdi Ramezani, Mehdi Dehghan, and Mohsen Razzaghi. Combined finite difference and spectral methods for the numerical solution of hyperbolic equation with an integral condition. Numerical Methods for Partial Differential Equations: An International Journal , 24(1):1–8, 2008
work page 2008
-
[40]
A mathematical model for brain tumor response to radiation therapy
R Rockne, EC Alvord, JK Rockhill, and KR Swanson. A mathematical model for brain tumor response to radiation therapy. Journal of mathematical biology , 58(4-5):561, 2009. Title Suppressed Due to Excessive Length 35
work page 2009
-
[41]
Models and the history of modeling
Hermann Schichl. Models and the history of modeling. In Modeling languages in mathematical optimization , pages 25–36. Springer, 2004
work page 2004
-
[42]
Of mice and men: Sparse statistical modeling in cardiovascular genomics
David M Seo, Pascal J Goldschmidt-Clermont, Mike West, et al. Of mice and men: Sparse statistical modeling in cardiovascular genomics. The Annals of Applied Statistics , 1(1):152–178, 2007
work page 2007
-
[43]
Spectral methods: algorithms, analysis and applications, volume 41
Jie Shen, Tao Tang, and Li-Lian Wang. Spectral methods: algorithms, analysis and applications, volume 41. Springer Science & Business Media, 2011
work page 2011
-
[44]
Role of free radical in atherosclerosis, diabetes and dyslipidaemia: larger-than-life
Randhir Singh, Sushma Devi, and Rakesh Gollen. Role of free radical in atherosclerosis, diabetes and dyslipidaemia: larger-than-life. Diabetes/metabolism research and reviews, 31(2):113–126, 2015
work page 2015
-
[45]
Ang´ elique St´ ephanou, Steven R McDougall, Alexander RA Anderson, and Mark AJ Chaplain. Mathematical modelling of the influence of blood rheological properties upon adaptative tumour-induced angiogenesis. Mathematical and Computer Modelling, 44(1-2):96–123, 2006
work page 2006
-
[46]
Emmanuelle Terry, Jacqueline Marvel, Christophe Arpin, Olivier Gandrillon, and Fabien Crauste. Mathematical model of the primary cd8 t cell immune response: stability analysis of a nonlinear age-structured system. Journal of mathematical biology, 65(2):263–291, 2012
work page 2012
-
[47]
Mathematical modeling and simulation of the evolution of plaques in blood vessels
Yifan Yang, Willi J¨ ager, Maria Neuss-Radu, and Thomas Richter. Mathematical modeling and simulation of the evolution of plaques in blood vessels. Journal of mathematical biology , 72(4):973–996, 2016
work page 2016
-
[48]
The growth of population and the factors which control it
Ge Yule. The growth of population and the factors which control it . Harrison & Sons, 1925
work page 1925
-
[49]
Yanmin Zhao, Yadong Zhang, Fawang Liu, I Turner, Yifa Tang, and V Anh. Convergence and superconvergence of a fully-discrete scheme for multi-term time fractional diffusion equations. Computers & Mathematics with Applications, 73(6):1087–1099, 2017
work page 2017
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.