Farzaneh Nasresfahani,M.R. Eslahchi
Abstract:The main target of this paper is to present an efficient method to solve a nonlinear free boundary mathematical model of prostate tumor. This model consists of two parabolics, one elliptic and one ordinary differential equations that are coupled together and describe the growth of a prostate tumor. We start our discussion by using the front fixing method to fix the free domain. Then, after employing a nonclassical finite difference and the collocation methods on this model, their stability and convergence are proved analytically. Finally, some numerical results are considered to show the efficiency of the mentioned methods.
What problem does this paper attempt to address?
### What problems does this paper attempt to solve?
This paper aims to solve the problems of numerical convergence and stability in the mathematical model of prostate cancer. Specifically, the paper proposes an effective method to solve the nonlinear free - boundary mathematical model of prostate tumors. This model consists of two parabolic equations, one elliptic equation and one ordinary differential equation. These equations are coupled with each other and describe the growth process of prostate tumors.
#### Main objectives:
1. **Fix the free boundary**: Use the front - fixing method to transform the free - boundary problem into a fixed - domain problem for more convenient numerical analysis.
2. **Construct a numerical method**: Apply the non - classical finite - difference method and the collocation method to solve the coupled partial differential equations and prove the convergence of the constructed sequence.
3. **Prove stability and convergence**: Prove the stability and convergence of the proposed method through theoretical analysis.
4. **Show numerical results**: Verify the effectiveness of the proposed method through numerical experiments.
#### Model background:
- **Prostate cancer**: Prostate cancer is the second most common type of cancer in men. Its cells are sensitive to androgens and depend on androgens for growth and survival. Treatment methods such as androgen - deprivation therapy (ADT) and total androgen blockade (TAB) are effective for a certain period of time, but often recur because tumor cells may evolve into androgen - independent (AI) cells.
- **Mathematical model**: The model used in this paper considers the mutation of androgen - dependent (AD) cells to androgen - independent (AI) cells and the influence of inhibitors. The model assumes that the tumor is a radially symmetric sphere with a radius \(R(t)\) that changes over time.
#### Equations of the mathematical model:
The model includes the following four equations:
1. **Parabolic equation (AD cells)**:
\[
\frac{\partial P}{\partial t}+\frac{1}{r^{2}}\frac{\partial}{\partial r}(r^{2}u(r,t)P(r,t))=\frac{D_{p}}{r^{2}}\frac{\partial}{\partial r}\left(r^{2}\frac{\partial P}{\partial r}\right)+\alpha_{p}(a(t))P(r,t)-\delta_{p}(a(t))P(r,t)-(1 - I)\beta(a(t))P(r,t)
\]
Initial and boundary conditions:
\[
P(r,0)=P_{0}(r),\quad\frac{\partial P}{\partial r}(0,t)=\frac{\partial P}{\partial r}(R(t),t)=0
\]
2. **Parabolic equation (AI cells)**:
\[
\frac{\partial q}{\partial t}+\frac{1}{r^{2}}\frac{\partial}{\partial r}(r^{2}u(r,t)q(r,t))=\frac{D_{q}}{r^{2}}\frac{\partial}{\partial r}\left(r^{2}\frac{\partial q}{\partial r}\right)+\alpha_{q}(a(t))q(r,t)-\delta_{q}(a(t))q(r,t)+(1 - I)\beta(a(t))p(r,t)
\]
Initial and boundary conditions:
\[
q(r,0)=q_{0}(r),\quad\frac{\partial q}{\partial r}(0,t)=\frac{\partial q}{\partial r}(R(t),t)=0
\]
3. **Elliptic equation**: