# American Institute of Mathematical Sciences

June  2019, 12(3): 591-613. doi: 10.3934/dcdss.2019038

## Efficient numerical method for a model arising in biological stoichiometry of tumour dynamics

 Department of Mathematics and Applied Mathematics, University of the Western Cape, Private Bag X17, Bellville, 7535, South Africa

* Corresponding author: mkowolax@yahoo.com (K. M. Owolabi)

Received  May 2017 Revised  September 2017 Published  September 2018

Fund Project: The research contained in this report is supported by South African National Research Foundation.

In this paper, we extend a system of coupled first order non-linear system of delay differential equations (DDEs) arising in modeling of stoichiometry of tumour dynamics, to a system of diffusion-reaction system of partial delay differential equations (PDDEs). Since tumor cells are further modified by blood supply through the vascularization process, we determine the local uniform steady states of the homogeneous tumour growth model with respect to the vascularization process. We show that the steady states are globally stable, determine the existence of Hopf bifurcation of the homogeneous tumour growth model with respect to the vascularization process. We derive, analyse and implement a fitted operator finite difference method (FOFDM) to solve the extended model. This FOFDM is analyzed for convergence and we observe seen that it has second-order accuracy. Some numerical results confirming theoretical observations are also presented. These results are comparable with those obtained in the literature.

Citation: Kolade M. Owolabi, Kailash C. Patidar, Albert Shikongo. Efficient numerical method for a model arising in biological stoichiometry of tumour dynamics. Discrete & Continuous Dynamical Systems - S, 2019, 12 (3) : 591-613. doi: 10.3934/dcdss.2019038
##### References:

show all references

##### References:
Numerical solution for the dynamics of homogeneous tumour growth model, when $a = 3,d_x = 2,b_1 = 6,d_1 = 0.5$
Numerical solution for the dynamics of homogeneous tumour growth model, when $a = 6,d_x = 1,b_1 = 6,d_1 = 1$
Numerical solution for the dynamics of homogeneous tumour growth model, when $a = 6,d_x = 0.5,b_1 = 3,d_1 = 2$
Parameter values[7]
 $m=20.00$ $n=10.00$ $k_h=10.00$ $k_t=3.00$ $f=0.6667$ $P=150$ $m_1=20.00$ $\beta_1=1.00$ $c=0.005$ $dz=0.20$ $g=100.00$ $\alpha=0.05$
 $m=20.00$ $n=10.00$ $k_h=10.00$ $k_t=3.00$ $f=0.6667$ $P=150$ $m_1=20.00$ $\beta_1=1.00$ $c=0.005$ $dz=0.20$ $g=100.00$ $\alpha=0.05$
 [1] Yang Kuang, John D. Nagy, James J. Elser. Biological stoichiometry of tumor dynamics: Mathematical models and analysis. Discrete & Continuous Dynamical Systems - B, 2004, 4 (1) : 221-240. doi: 10.3934/dcdsb.2004.4.221 [2] Joseph A. Connolly, Neville J. Ford. Comparison of numerical methods for fractional differential equations. Communications on Pure & Applied Analysis, 2006, 5 (2) : 289-307. doi: 10.3934/cpaa.2006.5.289 [3] Nicola Guglielmi, Christian Lubich. Numerical periodic orbits of neutral delay differential equations. Discrete & Continuous Dynamical Systems - A, 2005, 13 (4) : 1057-1067. doi: 10.3934/dcds.2005.13.1057 [4] Iasson Karafyllis, Lars Grüne. Feedback stabilization methods for the numerical solution of ordinary differential equations. Discrete & Continuous Dynamical Systems - B, 2011, 16 (1) : 283-317. doi: 10.3934/dcdsb.2011.16.283 [5] Alfonso Ruiz-Herrera. Chaos in delay differential equations with applications in population dynamics. Discrete & Continuous Dynamical Systems - A, 2013, 33 (4) : 1633-1644. doi: 10.3934/dcds.2013.33.1633 [6] Xiaoying Han, Jinglai Li, Dongbin Xiu. Error analysis for numerical formulation of particle filter. Discrete & Continuous Dynamical Systems - B, 2015, 20 (5) : 1337-1354. doi: 10.3934/dcdsb.2015.20.1337 [7] Tomás Caraballo, P.E. Kloeden, Pedro Marín-Rubio. Numerical and finite delay approximations of attractors for logistic differential-integral equations with infinite delay. Discrete & Continuous Dynamical Systems - A, 2007, 19 (1) : 177-196. doi: 10.3934/dcds.2007.19.177 [8] Eugen Stumpf. Local stability analysis of differential equations with state-dependent delay. Discrete & Continuous Dynamical Systems - A, 2016, 36 (6) : 3445-3461. doi: 10.3934/dcds.2016.36.3445 [9] Roman Srzednicki. A theorem on chaotic dynamics and its application to differential delay equations. Conference Publications, 2001, 2001 (Special) : 362-365. doi: 10.3934/proc.2001.2001.362 [10] Gang Huang, Yasuhiro Takeuchi, Rinko Miyazaki. Stability conditions for a class of delay differential equations in single species population dynamics. Discrete & Continuous Dynamical Systems - B, 2012, 17 (7) : 2451-2464. doi: 10.3934/dcdsb.2012.17.2451 [11] Qiumei Huang, Xiuxiu Xu, Hermann Brunner. Continuous Galerkin methods on quasi-geometric meshes for delay differential equations of pantograph type. Discrete & Continuous Dynamical Systems - A, 2016, 36 (10) : 5423-5443. doi: 10.3934/dcds.2016039 [12] Evelyn Sander, E. Barreto, S.J. Schiff, P. So. Dynamics of noninvertibility in delay equations. Conference Publications, 2005, 2005 (Special) : 768-777. doi: 10.3934/proc.2005.2005.768 [13] Ilknur Koca. Numerical analysis of coupled fractional differential equations with Atangana-Baleanu fractional derivative. Discrete & Continuous Dynamical Systems - S, 2019, 12 (3) : 475-486. doi: 10.3934/dcdss.2019031 [14] Marek Bodnar, Monika Joanna Piotrowska, Urszula Foryś, Ewa Nizińska. Model of tumour angiogenesis -- analysis of stability with respect to delays. Mathematical Biosciences & Engineering, 2013, 10 (1) : 19-35. doi: 10.3934/mbe.2013.10.19 [15] Philip Gerlee, Alexander R. A. Anderson. Diffusion-limited tumour growth: Simulations and analysis. Mathematical Biosciences & Engineering, 2010, 7 (2) : 385-400. doi: 10.3934/mbe.2010.7.385 [16] H.M. Byrne, S.M. Cox, C.E. Kelly. Macrophage-tumour interactions: In vivo dynamics. Discrete & Continuous Dynamical Systems - B, 2004, 4 (1) : 81-98. doi: 10.3934/dcdsb.2004.4.81 [17] Sun Yi, Patrick W. Nelson, A. Galip Ulsoy. Delay differential equations via the matrix lambert w function and bifurcation analysis: application to machine tool chatter. Mathematical Biosciences & Engineering, 2007, 4 (2) : 355-368. doi: 10.3934/mbe.2007.4.355 [18] Evelyn Buckwar, Girolama Notarangelo. A note on the analysis of asymptotic mean-square stability properties for systems of linear stochastic delay differential equations. Discrete & Continuous Dynamical Systems - B, 2013, 18 (6) : 1521-1531. doi: 10.3934/dcdsb.2013.18.1521 [19] Abdelhai Elazzouzi, Aziz Ouhinou. Optimal regularity and stability analysis in the $\alpha-$Norm for a class of partial functional differential equations with infinite delay. Discrete & Continuous Dynamical Systems - A, 2011, 30 (1) : 115-135. doi: 10.3934/dcds.2011.30.115 [20] Iacopo P. Longo, Sylvia Novo, Rafael Obaya. Topologies of continuity for Carathéodory delay differential equations with applications in non-autonomous dynamics. Discrete & Continuous Dynamical Systems - A, 2019, 39 (9) : 5491-5520. doi: 10.3934/dcds.2019224

2019 Impact Factor: 1.233