<div dir="ltr"><div><div><div><div>Dear Frederic,<br><br></div>thanks for your reply. There must be something else here, because VolAxiSqu hasn't helped.<br><br></div>I think the time function in my code might be okay, because it seems to work fine for Jacobi Vol.<br><br></div>Cheers,<br><br></div>Peter<br></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Aug 8, 2015 at 5:51 PM, Frederic Trillaud <span dir="ltr"><<a href="mailto:ftrillaudp@pumas.ii.unam.mx" target="_blank">ftrillaudp@pumas.ii.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 Peter,<br>
<br>
I had a similar problem, you should try VolAxiSqu. It solves my issue on<br>
the axis of symmetry.<br>
<br>
Jacobian<br>
{<br>
 {// Dimension ND regions in plane<br>
  Name AxiJacobian;<br>
  Case<br>
  {<br>
   {<br>
    Region infiniteBoundaryRegion;<br>
    Jacobian VolAxiSquSphShell{0.95*infiniteBoundaryInnerRadius,<br>
1.05*infiniteBoundaryOuterRadius};<br>
   }<br>
   {<br>
    Region All;<br>
    Jacobian VolAxiSqu;<br>
   }<br>
  }<br>
 }<br>
}<br>
<br>
To apply the current density, you need to use time function that will be<br>
multiply to your current density. You have to implement it in the<br>
Constraint.<br>
<br>
Constraint<br>
{<br>
{<br>
  Name currentDensityCondition;<br>
  Case<br>
  {<br>
   {<br>
    Region inductorRegion;<br>
    Value currentDensity[];<br>
    TimeFunction timeFunction[];<br>
   }<br>
  }<br>
 }<br>
}<br>
<br>
I use the conditional operators to create the time function:<br>
<br>
"timeFunction[] = If condition Then statement Else statement;" is then<br>
in conditional operators:<br>
timeFunction[] = (condition) ? statement : statement;<br>
<br>
<a href="http://www.c4learn.com/c-programming/c-conditional-operators/" rel="noreferrer" target="_blank">http://www.c4learn.com/c-programming/c-conditional-operators/</a><br>
<br>
Good luck,<br>
<br>
Frederic<br>
<div><div class="h5"><br>
-----Original Message-----<br>
From: Peter Kis <<a href="mailto:md2z34@gmail.com">md2z34@gmail.com</a>><br>
To: <a href="mailto:getdp@geuz.org">getdp@geuz.org</a><br>
Subject: [Getdp] 2D axisym. coil<br>
Date: Sat, 8 Aug 2015 16:57:49 +0200<br>
<br>
Dear getDP list members,<br>
<br>
<br>
I am still struggling with my first model in getDP, which would be a<br>
single turn coil on 2D axisymmetric plane.<br>
<br>
<br>
Now, at least I get something, but I don't like what I get. I am not<br>
sure how to apply the time dependent excitation current density to the<br>
coil. Therefore I took something over from CoreMassive model (found on<br>
the Internet).<br>
<br>
<br>
As you can see in the attached pictures the magnetic vector potential<br>
for the pure 2D model (with Jac Vol) looks good but for the axisym model<br>
(with Jac VolAxi) it looks terrible.<br>
<br>
<br>
Definitely there is still something here what I miss.<br>
<br>
<br>
Any comment is highly appreciated.<br>
<br>
<br>
Thanks in advance,<br>
<br>
<br>
Peter<br>
<br>
</div></div>_______________________________________________<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" rel="noreferrer" target="_blank">http://www.geuz.org/mailman/listinfo/getdp</a><br>
<br>
</blockquote></div><br></div>