[schradc]

Hello,

I am trying to calculate the magnetic field induced by a current density field given by



Where B is the magnetic field, J is the current density and mu is the magnetic permiability.

I have implemented this equation as

solve(fvm::laplacian(B) == -magPerm * fvc::curl(Je));

I set up a simple case with uniform current density traveling down a circular conductor. The anticipated results are a magnetic field which approaches 0 on the centerline and increases in strength to the outer surface. Additionally, the magnetic field is expected to have a clockwise rotation.

This expected behavior is observed when run in serial as shown in the first two images.

However, when run in parallel, I get unphysical results with a pattern consistent with the borders between processors. The subdomains for the 4 processors can be seen in the third attachment and the magnetic field results in the last two figures.

Does anyone have any ideas as to why this behavior would occur in parallel? Is there anything wrong with the line I implemented?

Many thanks.