[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
>
>
>