Implicit and Semi-implicit Numerical Schemes for the Gradient Flow of the Formation of Biological Transport Networks
The SMAI journal of computational mathematics, Volume 5 (2019) , pp. 229-249.

Implicit and semi-implicit time discretizations are developed for the Cai–Hu model describing the formation of biological transport networks. The model couples a nonlinear elliptic equation for the pressure with a nonlinear reaction-diffusion equation for the network conductance vector. Numerical challenges include the nonlinearity and the stiffness, thus an explicit discretization puts severe constraints on the time step. We propose an implicit and a semi-implicit discretizations, which decays the energy unconditionally or under a condition independent of the mesh size respectively, as will be proven in 1D and verified numerically in 2D.

Published online:
DOI: https://doi.org/10.5802/smai-jcm.59
Classification: 65M06,  92B99
Keywords: biological transport networks, gradient flow, numerical schemes
@article{SMAI-JCM_2019__5__229_0,
author = {Di Fang and Shi Jin and Peter Markowich and Beno{\^\i}t Perthame},
title = {Implicit and {Semi-implicit} {Numerical} {Schemes} for the {Gradient} {Flow} of the {Formation} of {Biological} {Transport} {Networks}},
journal = {The SMAI journal of computational mathematics},
pages = {229--249},
publisher = {Soci\'et\'e de Math\'ematiques Appliqu\'ees et Industrielles},
volume = {5},
year = {2019},
doi = {10.5802/smai-jcm.59},
language = {en},
url = {https://smai-jcm.centre-mersenne.org/articles/10.5802/smai-jcm.59/}
}
Di Fang; Shi Jin; Peter Markowich; Benoît Perthame. Implicit and Semi-implicit Numerical Schemes for the Gradient Flow of the Formation of Biological Transport Networks. The SMAI journal of computational mathematics, Volume 5 (2019) , pp. 229-249. doi : 10.5802/smai-jcm.59. https://smai-jcm.centre-mersenne.org/articles/10.5802/smai-jcm.59/

[1] G. Albi; M. Artina; M. Foransier; P. Markowich Biological transportation networks: modeling and simulation, Anal. Appl. (Singap.), Volume 14 (2016) no. 1, pp. 185-206 | Article | MR 3438650 | Zbl 1329.35313

[2] G. Albi; M. Burger; J. Haskovec; P. Markowich; M. Schlottbom Continuum modeling of biological network formation, Model. Simul. Sci. Eng. Technol. (2017), pp. 1-48 | MR 3644587

[3] S. Bohn; M. O. Magnasco Structure, Scaling, and Phase Transition in the Optimal Transport Network, Phys. Rev. Lett., Volume 98 (2007), 088702 pages | Article

[4] M. Burger; J. Haskovec; P. Markowich; H. Ranetbauer A mesoscopic model of biological transportation networks (2018) (https://arxiv.org/abs/1806.00120) | Zbl 07148389

[5] S. D. Conte Elementary numerical analysis: An algorithmic approach, McGraw-Hill Book Co., 1965, x+278 pages | MR 0202267 | Zbl 0213.41501

[6] F. Corson Fluctuations and Redundancy in Optimal Transport Networks, Phys. Rev. Lett., Volume 104 (2010), 048703 pages | Article

[7] P. Dimitrov; S. W. Zucker A constant production hypothesis guides leaf venation patterning, Proceedings of the National Academy of Sciences, Volume 103 (2006) no. 24, pp. 9363-9368 http://www.pnas.org/content/103/24/9363 | arXiv:http://www.pnas.org/content/103/24/9363.full.pdf | Article

[8] Q. Du; W. Zhu Stability analysis and application of the exponential time differencing schemes, J. Comput. Math., Volume 22 (2004) no. 2, pp. 200-209 (Special issue dedicated to the 70th birthday of Professor Zhong-Ci Shi) | MR 2058932 | Zbl 1052.65081

[9] Q. Du; W. Zhu Analysis and applications of the exponential time differencing schemes and their contour integration modifications, BIT, Volume 45 (2005) no. 2, pp. 307-328 | Article | MR 2176196 | Zbl 1080.65074

[10] X. Feng; H. Song; T. Tang; J. Yang Nonlinear stability of the implicit-explicit methods for the Allen-Cahn equation, Inverse Probl. Imaging, Volume 7 (2013) no. 3, pp. 679-695 | Article | MR 3105349

[11] J. Haskovec; L. M. Kreusser; P. Markowich ODE and PDE based modeling of biological transportation networks (2018) (https://arxiv.org/abs/1805.08526) | Zbl 07148390

[12] J. Haskovec; L. M. Kreusser; P. Markowich Rigorous Continuum Limit for the Discrete Network Formation Problem (2018) (https://arxiv.org/abs/1808.01526)

[13] J. Haskovec; P. Markowich; B. Perthame Mathematical analysis of a PDE system for biological network formation, Commun. Partial Differ. Equations, Volume 40 (2015) no. 5, pp. 918-956 | Article | MR 3306822 | Zbl 1345.35120

[14] J. Haskovec; P. Markowich; B. Perthame; M. Schlottbom Notes on a PDE system for biological network formation, Nonlinear Anal., Theory Methods Appl., Volume 138 (2016), pp. 127-155 | Article | MR 3485142 | Zbl 1334.35101

[15] D. Hu; D. Cai Adaptation and Optimization of Biological Transport Networks, Phys. Rev. Lett., Volume 111 (2013), 138701 pages | Article

[16] E. Katifori; Gergely J. Szöllősi; M. O. Magnasco Damage and Fluctuations Induce Loops in Optimal Transport Networks, Phys. Rev. Lett., Volume 104 (2010), 048704 pages | Article

[17] A. Koller; G. Kaley Endothelial regulation of wall shear stress and blood flow in skeletal muscle microcirculation, Am. J. Physiol. Heart Circ. Physiol., Volume 260 (1991) no. 3, p. H862-H868 (PMID: 2000980) | Article

[18] M. F. Laguna; S. Bohn; E. A. Jagla The Role of Elastic Stresses on Leaf Venation Morphogenesis, PLoS Comput. Biol., Volume 4 (2008) no. 4, pp. 1-10 | Article | MR 2399338

[19] D. Lu; G. S. Kassab Role of shear stress and stretch in vascular mechanobiology, J. R. Soc. Interface, Volume 8 (2011) no. 63, pp. 1379-1385 | Article

[20] C. D. Murray The Physiological Principle of Minimum Work: I. The Vascular System and the Cost of Blood Volume, Proc. Natl. Acad. Sci. USA, Volume 12 (1926) no. 3, pp. 207-214 | Article

[21] T. Nelson; N. Dengler Leaf Vascular Pattern Formation., The Plant Cell, Volume 9 (1997) no. 7, pp. 1121-1135 | Article

[22] U. Pohl; J. Holtz; R. Busse; E. Bassenge Crucial role of endothelium in the vasodilator response to increased flow in vivo., Hypertension, Volume 8 (1986) no. 1, pp. 37-44 | Article