This paper deals with a set of three partial Corrective Skin Peels differential equations involving time-fractional derivatives and nonlinear diffusion operators.This model helps us to understand the HIV spread and transmission into the patient.First, we prove the existence and uniqueness of weak solutions to the mathematical model.Then, the Galerkin finite element scheme is implemented to approximate the solution of the model.Further, a-priori error bounds and Developer and Photoconductor Unit Pack convergence estimates for the fully-discrete problem are derived.
The second order convergence for the proposed scheme is also proved.Numerical tests are shown to validate the theoretical studies.