I am trying to solve the radially symmetric polar form of the PDE with homogeneous Neumann BC in deal.II on a unit circle: $$ u_t = \Delta u - \nabla \cdot (u \nabla h) $$ $$ h_t = \Delta h $$
I am not sure how to handle the origin & how the bilinear form change in this case.
Also, how to use the grid and triangulation in deal.II.
Could someone please explain or suggest similar literature?