Spatiotemporal dynamics of RNA viruses in the presence of immunity and treatment: case of SARS-CoV-2

In this paper, we develop a mathematical model using partial differential equations to investigate the behavior of RNA viruses in the presence of antiviral treatment.  The developed model includes both cell-to-cell and virus-to-cell modes of transmission.  Initially, we establish the well-posedness of the model by demonstrating the existence and uniqueness of solutions, as well as their positivity and boundedness.  Additionally, we identify and analyze the stable equilibrium states, their global stability depending on specific threshold parameters, using Lyapunov functions.  To corroborate our theoretical findings, we provide illustrations through numerical simulations.