Sandra Vaz,Delfim F. M. Torres
Abstract:In this work, we derive a nonstandard finite difference scheme for the SICA (Susceptible-Infected-Chronic-AIDS) model and analyze the dynamical properties of the discretized system. We prove that the discretized model is dynamically consistent with the continuous, maintaining the essential properties of the standard SICA model, namely, the positivity and boundedness of the solutions, equilibrium points, and their local and global stability.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop a non - standard finite - difference (NSFD) scheme that is dynamically consistent with the continuous model to discretize the SICA (Susceptible - Infected - Chronic - AIDS) model and analyze the dynamic characteristics of the discretized system. Specifically, the authors aim to:
1. **Maintain the key properties of the continuous model**: Ensure that the discretized model can maintain the basic properties of the continuous SICA model, including the positivity and boundedness of the solution, equilibrium points and their local and global stabilities.
2. **Eliminate numerical instabilities**: By using the NSFD scheme constructed by the Mickens method, reduce or eliminate the oscillation and instability phenomena that may occur when other numerical methods (such as the Euler method and the Runge - Kutta method) are used to solve nonlinear systems.
### Model background
The SICA model describes the transmission dynamics of HIV/AIDS in the population and divides the population into four mutually exclusive categories:
- Susceptible individuals ($S$)
- HIV - infected individuals (without AIDS clinical symptoms, denoted as $I$)
- HIV - infected individuals receiving ART treatment (with low viral load, denoted as $C$)
- HIV - infected individuals with AIDS clinical symptoms (denoted as $A$)
### Main research content
To achieve the above goals, the authors carried out the following work:
1. **Propose a discretization scheme**: Based on the Mickens method, a dynamically consistent NSFD scheme is proposed to discretize the SICA model.
2. **Prove key properties**:
- **Positivity and boundedness**: It is proved that the discretized solution always maintains positivity and boundedness.
- **Existence of equilibrium points**: The equilibrium points of the discrete model are found and the existence of these equilibrium points is proved.
- **Stability analysis**: Using tools such as the Schur - Cohn criterion, the local and global stabilities of the discrete model are analyzed.
3. **Numerical simulation**: Through numerical simulation, the effectiveness of the proposed NSFD discretization scheme is verified, especially its superiority in maintaining the dynamic characteristics of the continuous model.
### Formula display
The following are some important formulas involved in the paper:
1. **Differential equations of the SICA model**:
\[
\begin{cases}
\dot{S}(t)=\Lambda-\lambda(t)S(t)-\mu S(t),\\
\dot{I}(t)=\lambda(t)S(t)-(\rho + \phi+\mu)I(t)+\alpha A(t)+\omega C(t),\\
\dot{C}(t)=\phi I(t)-(\omega + \mu)C(t),\\
\dot{A}(t)=\rho I(t)-(\alpha + \mu + d)A(t),
\end{cases}
\]
where,
\[
\lambda(t)=\beta\frac{I(t)+\eta_C C(t)+\eta_A A(t)}{N(t)}, \quad N(t)=S(t)+I(t)+C(t)+A(t).
\]
2. **Discretization scheme**:
\[
\begin{cases}
S_{n + 1}=S_n+\psi(\Lambda-\tilde{\lambda}_n S_{n + 1}-\mu S_{n + 1}),\\
I_{n + 1}=I_n+\psi(\tilde{\lambda}_n S_{n + 1}-(C_1)I_{n + 1}+\alpha A_{n + 1}+\omega C_{n + 1}),\\
C_{n + 1}=C_n+\psi(\phi I_{n + 1}-C_3 C_{n + 1})
\end{cases}
\]