Tumour dynamics and necrosis: surface tension and stability

A model is developed for the motion of cells within a multicell spherical tumour. The model allows either for the intercellular forces to be in compression and cells to be compacted to a fixed number density, or for the cell number density to fall and cells to become isolated from each other. The model develops necrotic regions naturally due to force balances rather than being directly attributable to a critical oxygen concentration. These necrotic regions may result in a gradual reduction in local cell density rather than jump to a completely dead region. Numerical and analytical analysis of the spherically symmetric model shows that the long time behaviour of the spheroid depends on any surface tension effects created by cells on the outer surface. For small surface tension the spheroid grows linearly in time developing a large necrotic region, while for larger surface tension the growth can be halted. The linear stability to spherically symmetric perturbations of all the possible resulting steady states is revealed.