[Getdp] issue with axisymmetric magnetostatic problem

Frederic Trillaud ftrillaudp at pumas.ii.unam.mx
Fri Dec 5 23:11:48 CET 2014


I have not since my geometry is spherical. I will try out to see, if it
comes from there.

Best,

Frederic


On vie, 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
>         
> 
>