besides the missing weighting factor of 2*Pi for cylindrically symmetric cases (see your post: viewtopic.php?f=3&t=6959&start=10#p22098
), the resulting component magnetic forces seems to be slightly different to the single process result but acceptable (Good results are
archived, if the distance between the loops <=1, for larger distances the mesh resultion must be better, I guess). I would say this is due to the fact, that there is virtually only one line with small errors within a wire.
The situation is different when looking at the vacuum component: As the picture indicates the resulting force in y-direction is about two orders of magnitudes different compared to the single process result (which is virtually close to zero). And this is -I would say- an absolute unphysical result. (Additional to this, there exist a resulting x-force at all components in both, parallel and single processed cases. This I do really not understand, but this is another issue.)
I compared the results for a larger test case including two thick coils. The results in y-direction are different but comparable. Here, the missing weighting factor is blowing up the correct solution, I guess.
But here comes the real problem:
If the nodal forces are coupled to the linear elasticity solver via Displacement Loads, the solver might not reach convergence if these sharing artefacts are accepted. In the small test cases the solution converges, but it takes much more iterations. In a huge case, it does not converge.