Modeling electromagnetics on cylindrical meshes with applications to steel-cased wells

j=I(δ(rrs+)δ(rrs))e=ϕ\begin{split} \nabla \cdot \vec{j} &= I\left(\delta(\vec{r} - \vec{r}_{s^{+}}) - \delta(\vec{r} - \vec{r}_{s^{-}})\right) \\ \vec{e} &= - \nabla \phi \end{split}