Finite volume simulation of a semi-linear Neumann problem (Keller-Segel model) on rectangular domains

In this study, the finite volume method is implemented for solving the problem of the semilinear equation: -d δ u+ u=uq (d, q>0)with a homogeneous Neumann boundary condition. This problem is equivalent to the known stationary Keller-Segel model, which arises in chemotaxis.After discretization, a nonlinear algebraic system is obtained and solved on the platform Matlab. As a result, many single peaked and multi-peaked shapes in 3D and contour plots can be drawn depending on the parameters d and q.