Numerical simulation of a prostate tumor growth model by the RBF-FD scheme and a semi-implicit time discretization

• Published in: Journal of computational and applied mathematics Vol. 388; p. 113314
• Main Authors: Mohammadi, Vahid, Dehghan, Mehdi, De Marchi, Stefano
• Format: Journal Article
• Language: English
• Published: Elsevier B.V 01.05.2021
• Subjects: A prostate tumor growth model based on time-dependent partial differential equations; Biconjugate gradient stabilized method; Mathematical oncology; Moving interface problem; Non-dimensionalization technique; Radial basis function-generated finite difference scheme; Non-dimensionalization technique; Moving interface problem; 35Q35; Mathematical oncology; Biconjugate gradient stabilized method; 92C17; A prostate tumor growth model based on time-dependent partial differential equations; 35Q92; Radial basis function-generated finite difference scheme
• ISSN: 0377-0427; 1879-1778
• DOI: 10.1016/j.cam.2020.113314

The aim of this work consists of finding a suitable numerical method for the solution of the mathematical model describing the prostate tumor growth, formulated as a system of time-dependent partial differential equations (PDEs), which plays a key role in the field of mathematical oncology. In the literature on the subject, there are a few numerical methods for solving the proposed mathematical model. Localized prostate cancer growth is known as a moving interface problem, which must be solved in a suitable stable way. The mathematical model considered in this paper is a system of time-dependent nonlinear PDEs that describes the interaction between cancer cells, nutrients, and prostate-specific antigen (PSA). Here, we first derive a non-dimensional form of the studied mathematical model using the well-known non-dimensionalization technique, which makes it easier to implement different numerical techniques. Afterward, the analysis of the numerical method describing the two-dimensional prostate tumor growth problem, based on radial basis function-generated finite difference (RBF-FD) scheme, in combination with a first-order time discretization has been done. The numerical technique we use, does not need the use of any adaptivity techniques to capture the features in the interface. The discretization leads to solving a linear system of algebraic equations solved via the biconjugate gradient stabilized (BiCGSTAB) algorithm with zero-fill incomplete lower–upper (ILU) preconditioner. Comparing the results obtained in this investigation with those reported in the recent literature, the proposed approach confirms the ability of the developed numerical scheme. Besides, the effect of choosing constant parameters in the mathematical model is verified by many simulations on rectangular and circular domains.