[Getdp] issue with axisymmetric magnetostatic problem

Frederic Trillaud Pighi ftrillaudp at pumas.iingen.unam.mx
Sun Dec 7 02:40:17 CET 2014


Thanks a lot using VolAxiSquSphShell and VolAxiSqu did the trick. Now it
works like a charm.

Best,

Frederic


On Fri, 2014-12-05 at 20:20 +0100, Jasper wrote:
> Hi Frederic,
> 
> 
> Have you tried using VolAxiSqu and VolAxiSquSphShell in the Jacobian?
> See:
> 
> 
> http://www.geuz.org/pipermail/getdp/2004/000581.html
> 
> 
> 
> Best regards,
> 
> 
> Jasper
> 
> 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
>         
> 
>