Axisymmetric solutions for a chemotaxis model of Multiple Sclerosis

In this paper we study radially symmetric solutions for our recently proposed reaction–diffusion–chemotaxis model of Multiple Sclerosis. Through a weakly nonlinear expansion we classify the bifurcation at the onset and derive the amplitude equations ruling the formation of concentric demyelinating patterns which reproduce the concentric layers observed in Balò sclerosis and in the early phase of Multiple Sclerosis. We present numerical simulations which illustrate and fit the analytical results.