[Getdp] issue with axisymmetric magnetostatic problem

[Frederic Trillaud Pighi]

Dear all,

I am scratching my head over the computation of the magnetic flux
density B using magnetic vector potential formulation. The code computes
the vector potential in 2D axisymmetry but when I tried to derive the
magnetic flux density from B = Curl A, the result shows diverging values
of B close to the axis of symmetry. I cannot get it work as expected.
Here are the geometry and the *pro file attached if someone could
pinpoint the issue!

Best,

Frederic
-------------- next part --------------
A non-text attachment was scrubbed...
Name: axiGeometry2D.geo
Type: text/x-csrc
Size: 6043 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20141203/d0ad6572/attachment.geo>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: axiGeometry2D.pro
Type: text/x-csrc
Size: 7033 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20141203/d0ad6572/attachment.pro>