<div dir="ltr">Hi Frederic,<div><br></div><div>Have you tried using VolAxiSqu and VolAxiSquSphShell in the Jacobian? See:</div><div><br></div><div><a href="http://www.geuz.org/pipermail/getdp/2004/000581.html">http://www.geuz.org/pipermail/getdp/2004/000581.html</a><br></div><div><br></div><div>Best regards,</div><div><br></div><div>Jasper</div></div><div class="gmail_extra"><br><div class="gmail_quote">On Thu, Dec 4, 2014 at 3:23 AM, Frederic Trillaud Pighi <span dir="ltr"><<a href="mailto:ftrillaudp@pumas.iingen.unam.mx" target="_blank">ftrillaudp@pumas.iingen.unam.mx</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear all,<br>
<br>
I am scratching my head over the computation of the magnetic flux<br>
density B using magnetic vector potential formulation. The code computes<br>
the vector potential in 2D axisymmetry but when I tried to derive the<br>
magnetic flux density from B = Curl A, the result shows diverging values<br>
of B close to the axis of symmetry. I cannot get it work as expected.<br>
Here are the geometry and the *pro file attached if someone could<br>
pinpoint the issue!<br>
<br>
Best,<br>
<br>
Frederic<br>
<br>_______________________________________________<br>
getdp mailing list<br>
<a href="mailto:getdp@geuz.org">getdp@geuz.org</a><br>
<a href="http://www.geuz.org/mailman/listinfo/getdp" target="_blank">http://www.geuz.org/mailman/listinfo/getdp</a><br>
<br></blockquote></div><br></div>