Numerical solution of a 3D mathematical model for the progression of tumor angiogenic factor in a tissue

Melike Keleş Duman; Serdal Pamuk

doi:10.59400/adecp2489

Melike Keleş Duman Department of Mathematics, University of Kocaeli (Umuttepe Campus), Kocaeli 41001, Turkey
Serdal Pamuk Department of Mathematics, University of Kocaeli (Umuttepe Campus), Kocaeli 41001, Turkey

Article ID: 2489

DOI: https://doi.org/10.59400/adecp2489

Keywords: method of lines; porosity; tumor angiogenic factor; numerical solution; porous medium

Abstract

In this work, the movement of tumor angiogenic factor in a three-dimensional tissue is obtained by the Method of Lines. This method transforms a partial differential equation into a system of ordinary differential equations together with the initial and boundary conditions. The more the number of lines is increased, the more the accuracy of the method increases. This method results in very accurate numerical solutions for linear and non-linear problems in contrast with other existing methods. We present Matlab-generated figures, which are the movement of tumor angiogenic factor in porous medium and explain the biological importance of this progression. The computer codes are also provided.

References

[1]Emamjome M, Azarnavid B, Ghehsareh HR. A reproducing kernel Hilbert space pseudospectral method for numerical investigation of a two-dimensional capillary formation model in tumor angiogenesis problem. Neural Computing and Applications. 2049; 31: 2233–2241.

[2]Pamuk S, Erdem A. The method of lines for the numerical solution of a mathematical model for capillary formation: The role of endothelial cells in the capillary. Applied Mathematics and Computation. 2007; 186(1): 831–835.

[3]Shakeri F, Dehghan M. The method of lines for solution of the one-dimensional wave equation subject to an integral conservation condition. Computers and Mathematics with Applications. 2008; 56(9): 2175–2188.

[4]Sharaf AA, Bakodah HO. A good spatial discretisation in the method of lines. Applied Mathematics and Computation. 2005; 171(2): 1253–1263.

[5]Verwer JG, Sanz-Serna JM. Convergence of method of lines approximations to partial differential equations. Computing. 1984; 33: 297–313. doi: 10.1007/BF02242274.

[6]Mousa MM. Efficient numerical scheme based on the method of lines for the shallow water equations. Journal of Ocean Engineering and Science. 2018; 3(4): 303–309.

[7]Ali MR. The Method of Lines Analysis of Heat Transfer of Ostwald-de Waele Fluid Generated by a Non-uniform Rotating Disk with a Variable Thickness. Journal of Applied and Computational Mechanics. 2021;7(2): 432–441. doi: 10.22055/JACM.2020.30890.1787

[8]Agud Albesa L, Boix García M, Pla Ferrando ML, Navarrete SCC. A study about the solution of convection–diffusion–reaction equation with Danckwerts boundary conditions by analytical, method of lines and Crank–Nicholson techniques. Mathematical Methods in the Applied Sciences. 2023; 46(2): 2133–2164. doi: 10.1002/mma.8633

[9]Shi A, Persson PO, Zahr MJ. Implicit shock tracking for unsteady flows by the method of lines. Journal of Computational Physics. 2022; 454: 110906. doi: 10.1016/j.jcp.2021.110906

[10]Spiller W. An enhancement of instruments for solution of general transmission line equations with method of lines, impedance-/admittance and field transformation in combination with finite differences. Optical and Quantum Electronics. 2022; 54: 591. doi: 10.1007/s11082-022-03894-3

[11]Kazem S, Dehghan M. Semi-analytical solution for time-fractional diffusion equation based on finite difference method of lines (MOL). Engineering with Computers. 2019; 35: 229–241.

[12]Salehi Y, Darvishi MT, Schiesser WE. Numerical solution of space fractional diffusion equation by the method of lines and splines. Applied Mathematics and Computation. 2018; 336: 465–480.

[13]Soradi-Zeid S, Mesrizadeh M. The method of lines for parabolic integro-differential equations. Journal of Mathematical Modeling. 2020; 8(3): 291–308.

[14]Bratsos AG. The solution of the two-dimensional sine-Gordon equation using the method of lines. Journal of Computational and Applied Mathematics. 2007; 206(1): 251–277.

[15]Subramanian VR, White RE. Semianalytical method of lines for solving elliptic partial differential equations. Chemical Engineering Science. 2004; 59(4): 781–788.

[16]Manshoor B, Salleh H, Khalid A, Abdelaal MAS. Method of Lines and Runge-Kutta Method in Solving Partial Differential Equation for Heat Equation. Journal of Complex Flow. 2021; 3(1): 21–25.

[17]Ahmad I, Berzins M. MOL solvers for hyperbolic PDEs with source terms. Mathematics and Computers in Simulation. 2001; 56(2): 115–125.

[18]Kazem S, Dehghan M. Application of finite difference method of lines on the heat equation. Numerical Methods Partial Differential Eq. 2018; 34(2): 626–660.

[19]Erdem A, Pamuk S. The method of lines for the numerical solution of a mathematical model for capillary formation: The role of tumor angiogenic factor in the extra-cellular matrix. Applied Mathematics and Computation. 2007; 186(1): 891–897.

[20]Eikenberry SE, Sankar T, Preul MC, et al. Virtual Glioblastoma: Growth, Migration and Treatment in a Three-Dimensional Mathematical Model. Cell Cell Proliferation. 2009; 42(4): 511–528.

[21]Perfahl H, Byrne HM, Chen T, et al. Multiscale Modelling of Vascular Tumour Growth in 3D: The Roles of Domain Size and Boundary Conditions. PLoS One. 2011; 6(4): e14790.

[22]Sakai K, Hayashi T, Sakai Y, et al. A three-dimensional model with two-body interactions for endothelial cells in angiogenesis. Scientific Reports. 2023; 13: 20549. doi: 10.1038/s41598-023-47911-1

[23]Shirinifard A, Gens JS, Zaitlen BL, et al. 3D Multi-Cell Simulation of Tumor Growth and Angiogenesis. PLoS One. 2009; 4(10): e7190.

[24]Stéphanou A, McDougall SR, Anderson ARA, Chaplain MAJ. Mathematical Modelling of Flow in 2D and 3D Vascular Networks: Applications to Anti-Angiogenic and Chemotherapeutic Drug Strategies. Mathematical and Computer Modelling. 2005; 41(10): 1137–1156.

[25]Tang L, van de Ven AL, Guo D, et al. Computational Modeling of 3D Tumor Growth and Angiogenesis for Chemotherapy Evaluation. PLoS One. 2014; 9(1): e83962.

[26]Xie H, Jiao Y, Fan Q, et al. Modeling Three-Dimensional Invasive Solid Tumor Growth in Heterogeneous Microenvironment under Chemotherapy. PLoS One. 2018; 13(10): e0206292.

[27]Xu J, Vilanova G, Gomez H. Phase-Field Model of Vascular Tumor Growth: Three-Dimensional Geometry of the Vascular Network and Integration with Imaging Data. Computer Methods in Applied Mechanics and Engineering. 2020; 359: 112648.

[28]Zhao G, Chen E, Yu X, et al. Three-Dimensional Model of Metastatic Tumor Angiogenesis in Response to Anti-Angiogenic Factor Angiostatin. Journal of Mechanics in Medicine and Biology. 2017; 17(10): 1750094.

[29]Alias N, bin Masseri MIS, Islam MR, Khalid SN. The Visualization of Three Dimensional Brain Tumors’ Growth on Distributed Parallel Computer Systems. Journal of Applied Sciences. 2009; 9(3): 505–512.

[30]Karaa S, Zhang J, Yang F. A numerical study of a 3D bioheat transfer problem with different spatial heating. Mathematics and Computers in Simulation. 2005; 68(4): 375–388.

[31]Al Qubeissi M. Proposing a Numerical Solution for the 3D Heat Conduction Equation. In: Proceedings of the 2012 Sixth Asia Modelling Symposium; 29–31 May 2012; Bali, Indonesia.

[32]Tsega EG. Numerical Solution of Three-Dimensional Transient Heat Conduction Equation in Cylindrical Coordinates. Journal of Applied Mathematics. 2022; 1993151.

[33]Tsega EG. A Numerical Solution of Three-Dimensional Unsteady State Heat Equation. International Journal of Mathematical Modelling & Computations. 2021; 11(01): 49–60.

[34]Wang TY, Lee YM, Chen CP. 3D thermal-ADI—An efficient chip-level transient thermal simulator. In: Proceedings of the ISPD’ 03: International Symposium on Physical Design; 6–9 April 2003; Monterey, CA, USA.

[35]He CH, Liu HW, Liu C. A Fractal-Based Approach to the Mechanical Properties of Recycled Aggregate Concretes. Series: Mechanical Engineering. 2024; 22(2): 329–342. doi: 10.22190/FUME240605035H