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