Radial basis function-generated finite difference scheme for simulating the brain cancer growth model under radiotherapy in various types of computational domains

被引:9
|
作者
Dehghan, Mehdi [1 ]
Narimani, Niusha [1 ]
机构
[1] Amirkabir Univ Technol, Fac Math & Comp Sci, Dept Appl Math, 424 Hafez Ave, Tehran 15914, Iran
关键词
Numerical simulations; Radial basis function-generated finite difference scheme; Multiquadric radial functionbrain tumor; Reaction-diffusion models of radiation therapy; LINEAR-QUADRATIC MODEL; SOLID TUMOR-GROWTH; PATIENTS IN-VIVO; MATHEMATICAL-MODEL; GLIOMA GROWTH; NUMERICAL-SIMULATION; MESHLESS METHOD; CELL INVASION; DIFFUSION; INTERPOLATION;
D O I
10.1016/j.cmpb.2020.105641
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
Background and Objectives: We extend the original mathematical model, i.e., Swanson's reaction-diffusion equation to the surfaces with no boundary, and we find a new numerical method based on a meshless approach for solving numerically Swanson's reaction-diffusion model in the square and on the sphere. Methods: To solve numerically the Swanson's reaction-diffusion model and its extension version, a collocation meshless technique, namely radial basis function-generated finite difference (RBF-FD) scheme is employed for approximating the spatial variables in the square domain and on the sphere, respectively. Also, to approximate the time variable of the studied models, a first-order semi-implicit backward Euler scheme is used. The resulting fully discrete scheme is a linear system of algebraic equations per time step that is solved via the biconjugate gradient stabilized (BiCGSTAB) iterative algorithm with a zero-fill incomplete lower-upper (ILU) preconditioner. Results: The numerical simulations show the growth of untreated and treated brain tumors with radiotherapy using estimated and clinical data (given from magnetic resonance imaging (MRI) scans of patients). Moreover, the results reported here can be used for improving the treatment strategies of the invasive brain tumor. Conclusions: Using the developed numerical scheme in this paper, we can simulate the behavior of the invasive form of brain tumor response to radiotherapy. Also, we can see the effects of radiation response on the brain tumor cell concentration of individual patients. The proposed meshless technique, which is applied for solving numerically the studied model, does not depend on any background mesh or triangulation for approximation in comparison with mesh-dependent methods. Moreover, we apply this technique to the sphere via any set of distributed points easily. (C) 2020 Elsevier B.V. All rights reserved.
引用
收藏
页数:15
相关论文
共 9 条