Mathematical modeling of tumor surface growth with necrotic kernels

Abstract

A two‐dimensional tumor‐immune model with the time delay of the adaptive immune response is considered in this paper. The model is designed to account for the interaction between cytotoxic T lymphocytes (CTLs) and cancer cells on the surface of a solid tumor. The model considers the surface growth as a major growth pattern of solid tumors in order to describe the existence of necrotic kernels. The qualitative analysis shows that the immune‐free equilibrium is unstable, and the behavior of positive equilibrium is closely related to the ratio of the immune killing rate to tumor volume growth rate. The positive equilibrium is locally asymptotically stable when the ratio is smaller than a critical value. Otherwise, the occurrence of the delay‐driven Hopf bifurcation at the positive equilibrium is proved. Applying the center manifold reduction and normal form method, we obtain explicit formulas to determine the properties of the Hopf bifurcation. The global continuation of a local Hopf bifurcation is investigated based on the coincidence degree theory. The results reveal that the time of the adaptive immune system taken to response to tumors can lead to oscillation dynamics. We also carry out detailed numerical analysis for parameters and numerical simulations to illustrate our qualitative analysis. Numerically, we find that shorter immune response time can lead to longer patient survival time, and the period and amplitude of a stable periodic solution increase with the increasing immune response time. When CTLs recruitment rate and death rate vary, we show how the ratio of the immune killing rate to tumor volume growth rate and the first bifurcation value change numerically, which yields further insights to the tumor‐immune dynamics.