<div dir="ltr"><div><div><div><div><div><div>Dear Frederic,<br><br></div>a big thank to you for helping me out of this. The model with your modifications works like a charm. I am amazed how could you figured out the problem with the initial conditions. How could you catch this? <br><br></div>I compared the results given by getdp with COMSOL and it's a perfect match! I am very pleased.<br><br></div>I continue developing this model. Because I still need to add a special constraint which should enforce not the external current density to the coil but the net current. That is int_Omega_coil (Jext + sigma*E) = I, where E is the electric field E= -dA/dt. It's a kind of iteration thing because the total net current depends on the result (A). Do you think such constraint is possible in getDP? <br><br></div>Thanks again,<br><br></div>Kind Regards<br><br><br></div><div>Peter<br></div><div><div><div>  <br></div></div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sun, Aug 9, 2015 at 1:52 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 Peter,<br>
<br>
Here it goes. I made some changes (transformation with an outer shell,<br>
killed the initial time value of A, you should not need it) and did some<br>
tests so the parameters of current density etc. have been changed. It<br>
should work fine now.<br>
<span class="im HOEnZb"><br>
Best,<br>
<br>
Frederic<br>
<br>
<br>
On Sun, 2015-08-09 at 00:27 +0200, Peter Kis wrote:<br>
</span><div class="HOEnZb"><div class="h5">> I am sending the model files. Thanks for your time.<br>
><br>
><br>
> Peter<br>
><br>
><br>
> On Sat, Aug 8, 2015 at 11:02 PM, Frederic Trillaud<br>
> <<a href="mailto:ftrillaudp@pumas.ii.unam.mx">ftrillaudp@pumas.ii.unam.mx</a>> wrote:<br>
>         Send me back your files, I will take a look at it over the<br>
>         weekend to<br>
>         see if I can pin the problem.<br>
><br>
>         Frederic<br>
><br>
><br>
>         -----Original Message-----<br>
>         From: Peter Kis <<a href="mailto:md2z34@gmail.com">md2z34@gmail.com</a>><br>
>         To: Frederic Trillaud <<a href="mailto:ftrillaudp@pumas.ii.unam.mx">ftrillaudp@pumas.ii.unam.mx</a>><br>
>         Cc: <a href="mailto:getdp@geuz.org">getdp@geuz.org</a><br>
>         Subject: Re: [Getdp] 2D axisym. coil<br>
>         Date: Sat, 8 Aug 2015 22:43:11 +0200<br>
><br>
>         Dear Frederic,<br>
><br>
><br>
>         thanks for your reply. There must be something else here,<br>
>         because<br>
>         VolAxiSqu hasn't helped.<br>
><br>
><br>
>         I think the time function in my code might be okay, because it<br>
>         seems to<br>
>         work fine for Jacobi Vol.<br>
><br>
><br>
>         Cheers,<br>
><br>
><br>
>         Peter<br>
><br>
><br>
>         On Sat, Aug 8, 2015 at 5:51 PM, Frederic Trillaud<br>
>         <<a href="mailto:ftrillaudp@pumas.ii.unam.mx">ftrillaudp@pumas.ii.unam.mx</a>> wrote:<br>
>                 Dear Peter,<br>
><br>
>                 I had a similar problem, you should try VolAxiSqu. It<br>
>         solves my<br>
>                 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<br>
>         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<br>
>         function that<br>
>                 will be<br>
>                 multiply to your current density. You have to<br>
>         implement it in<br>
>                 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<br>
>         function:<br>
><br>
>                 "timeFunction[] = If condition Then statement Else<br>
>         statement;"<br>
>                 is then<br>
>                 in conditional operators:<br>
>                 timeFunction[] = (condition) ? statement : statement;<br>
><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>
><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,<br>
>         which would<br>
>                 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<br>
>         get. I am<br>
>                 not<br>
>                 sure how to apply the time dependent excitation<br>
>         current density<br>
>                 to the<br>
>                 coil. Therefore I took something over from CoreMassive<br>
>         model<br>
>                 (found on<br>
>                 the Internet).<br>
><br>
><br>
>                 As you can see in the attached pictures the magnetic<br>
>         vector<br>
>                 potential<br>
>                 for the pure 2D model (with Jac Vol) looks good but<br>
>         for the<br>
>                 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>
><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" rel="noreferrer" target="_blank">http://www.geuz.org/mailman/listinfo/getdp</a><br>
><br>
><br>
><br>
><br>
><br>
><br>
><br>
</div></div></blockquote></div><br></div>