Integrating the normal current density over an insulating surface should give zero as you expect. However, if you integrate the norm (magnitude) of the current density it should not be zero. That’s because there can be a current parallel to the insulating surface that will affect the integral.
Now reacf() is a tricky operator but very powerful as Ivar mentioned. I wish it was better documented. It provided the nodal “reaction forces” corresponding to fixed degrees of freedom. In solid mechanics, reacf(displacement-x) give the nodal reaction in the x direction, in fluid flow reacf(velocity-x) also gives the reaction in the x-direction and in your case reacf(Voltage) should give the nodal current normal to the boundary. That should be zero on an electrically insulated boundary.
Nagi Elabbasi
Veryst Engineering
Now reacf() is a tricky operator but very powerful as Ivar mentioned. I wish it was better documented. It provided the nodal “reaction forces” corresponding to fixed degrees of freedom. In solid mechanics, reacf(displacement-x) give the nodal reaction in the x direction, in fluid flow reacf(velocity-x) also gives the reaction in the x-direction and in your case reacf(Voltage) should give the nodal current normal to the boundary. That should be zero on an electrically insulated boundary.
Nagi Elabbasi
Veryst Engineering