A pressure model of immune response to mycobacterium tuberculosis infection in several space dimensions

Mycobacterium tuberculosis (Mtb) is a widely diffused infection. However, in general, the human immune system is able to contain it. In this work, we propose a mathematical model which describes the early immune response to the Mtb infection in the lungs, also including the possible evolution of the infection in the formation of a granuloma. The model is based on coupled reaction-diffusion-transport equations with chemotaxis, which take into account the interactions among bacteria, macrophages and chemoattractant. The novelty of this approach is in the modeling of the velocity field, proportional to the gradient of the pressure developed between the cells, which makes possible to deal with a full multidimensional description and efficient numerical simulations. We perform a linear stability analysis of the model and propose a robust implicit-explicit scheme to deal with long time simulations. Both in one and two-dimensions, we find that there are threshold values in the parameters space, between a contained infection and the uncontrolled bacteria growth, and the generation of granuloma-like patterns can be observed numerically.