[Getdp] issue with axisymmetric magnetostatic problem

Jasper funkybob at gmail.com
Fri Dec 5 20:20:54 CET 2014

Hi Frederic,

Have you tried using VolAxiSqu and VolAxiSquSphShell in the Jacobian? See:


Best regards,


On Thu, Dec 4, 2014 at 3:23 AM, Frederic Trillaud Pighi <
ftrillaudp at pumas.iingen.unam.mx> wrote:

> 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
> _______________________________________________
> getdp mailing list
> getdp at geuz.org
> http://www.geuz.org/mailman/listinfo/getdp
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/getdp/attachments/20141205/b43ec513/attachment.html>