Hi all, I'm simulating ferrohydrodynamics and I'm trying to impose a Volume Force in the Laminar Flow physics. This volume force has a part like this:
((1/r)(d(r(d(A,r))^2,r))+d(d(A,r)d(A,z),z))mu0_const
where "A" is the magnetic scalar potential. When I plot this term, the closest layer to the axis has not been taken into account actually (take a look at the attached image please, the black vertical line is the axis.) How can I prevent this problem? I have instabilities on the axis because of this issue. Thanks