I am solving Laplace equation on circular domain from which i get solution for potential.
The code I used:
w = ImplicitRegion[x^2 + y^2 <= 4, {{x, -2, 2}, {y, -2, 2}}]
op = -Laplacian[u[x, y], {x, y}]
h = {DirichletCondition[u[x, y] == 50, x^2 + y^2 == 4 && x > 1.99],
DirichletCondition[u[x, y] == -50, x^2 + y^2 == 4 && x < -1.99]};
uif = NDSolveValue[{op ==
NeumannValue[0, x^2 + y^2 == 4 && -1.9 < x < 1.9], h},
u, {x, y} \[Element] w]
ContourPlot[uif[x, y], {x, y} \[Element] w,
ColorFunction -> "Temperature", AspectRatio -> Automatic,
Contours -> 30]
Now i would like to calculate current density (gradient of potential) from solution for potential and plot it as vector field on the same graph. Can this be done with some built in function or how would i go about this?