我正在尝试在单位圆上使用均质 Neumann BC 求解 PDE 的径向对称极坐标形式：

$$u_t = \Delta u - \nabla \cdot (u \nabla h)$$ $$h_t = \Delta h$$

1. 我不确定如何处理原点以及在这种情况下双线性形式如​​何变化。

2. 另外，如何在deal.II中使用网格和三角剖分。

有人可以解释或建议类似的文献吗？