<html><head><meta http-equiv="Content-Type" content="text/html charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class="">Dear Peter,<div class=""><br class=""></div><div class="">It reminds me of the integral constraint in COMSOL. I do not know in GetDP. As far as I know, you can provide constant values for the constraint but I am not sure how you can use out of the box a function to implement the integral constraint. It may require some extra-implementation. The other option is to solve it as an extra equation: Int (Jext+sigma*E-J) = 0 where J is an homogeneous current density J = I/Omega_coil. Not sure if it makes sense for your problem.</div><div class=""><br class=""></div><div class="">Hopefully, someone has an appropriate answer to your problem. I would be interested to know how it could be implemented. I may need it down the road.</div><div class=""><br class=""></div><div class="">Best,</div><div class=""><br class=""></div><div class="">Frederic</div><div class=""><br class=""></div><div class=""><br class=""></div><div class=""><br class=""><div><blockquote type="cite" class=""><div class="">On Aug 9, 2015, at 4:28 PM, Peter Kis <<a href="mailto:md2z34@gmail.com" class="">md2z34@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div class=""><div class=""><div class=""><div class=""><div class=""><div class="">Dear Frederic,<br class=""><br class=""></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 class=""><br class=""></div>I compared the results given by getdp with COMSOL and it's a perfect match! I am very pleased.<br class=""><br class=""></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 class=""><br class=""></div>Thanks again,<br class=""><br class=""></div>Kind Regards<br class=""><br class=""><br class=""></div><div class="">Peter<br class=""></div><div class=""><div class=""><div class=""> <br class=""></div></div></div></div><div class="gmail_extra"><br class=""><div class="gmail_quote">On Sun, Aug 9, 2015 at 1:52 AM, Frederic Trillaud Pighi <span dir="ltr" class=""><<a href="mailto:ftrillaudp@pumas.iingen.unam.mx" target="_blank" class="">ftrillaudp@pumas.iingen.unam.mx</a>></span> wrote:<br class=""><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Dear Peter,<br class="">
<br class="">
Here it goes. I made some changes (transformation with an outer shell,<br class="">
killed the initial time value of A, you should not need it) and did some<br class="">
tests so the parameters of current density etc. have been changed. It<br class="">
should work fine now.<br class="">
<span class="im HOEnZb"><br class="">
Best,<br class="">
<br class="">
Frederic<br class="">
<br class="">
<br class="">
On Sun, 2015-08-09 at 00:27 +0200, Peter Kis wrote:<br class="">
</span><div class="HOEnZb"><div class="h5">> I am sending the model files. Thanks for your time.<br class="">
><br class="">
><br class="">
> Peter<br class="">
><br class="">
><br class="">
> On Sat, Aug 8, 2015 at 11:02 PM, Frederic Trillaud<br class="">
> <<a href="mailto:ftrillaudp@pumas.ii.unam.mx" class="">ftrillaudp@pumas.ii.unam.mx</a>> wrote:<br class="">
> Send me back your files, I will take a look at it over the<br class="">
> weekend to<br class="">
> see if I can pin the problem.<br class="">
><br class="">
> Frederic<br class="">
><br class="">
><br class="">
> -----Original Message-----<br class="">
> From: Peter Kis <<a href="mailto:md2z34@gmail.com" class="">md2z34@gmail.com</a>><br class="">
> To: Frederic Trillaud <<a href="mailto:ftrillaudp@pumas.ii.unam.mx" class="">ftrillaudp@pumas.ii.unam.mx</a>><br class="">
> Cc: <a href="mailto:getdp@geuz.org" class="">getdp@geuz.org</a><br class="">
> Subject: Re: [Getdp] 2D axisym. coil<br class="">
> Date: Sat, 8 Aug 2015 22:43:11 +0200<br class="">
><br class="">
> Dear Frederic,<br class="">
><br class="">
><br class="">
> thanks for your reply. There must be something else here,<br class="">
> because<br class="">
> VolAxiSqu hasn't helped.<br class="">
><br class="">
><br class="">
> I think the time function in my code might be okay, because it<br class="">
> seems to<br class="">
> work fine for Jacobi Vol.<br class="">
><br class="">
><br class="">
> Cheers,<br class="">
><br class="">
><br class="">
> Peter<br class="">
><br class="">
><br class="">
> On Sat, Aug 8, 2015 at 5:51 PM, Frederic Trillaud<br class="">
> <<a href="mailto:ftrillaudp@pumas.ii.unam.mx" class="">ftrillaudp@pumas.ii.unam.mx</a>> wrote:<br class="">
> Dear Peter,<br class="">
><br class="">
> I had a similar problem, you should try VolAxiSqu. It<br class="">
> solves my<br class="">
> issue on<br class="">
> the axis of symmetry.<br class="">
><br class="">
> Jacobian<br class="">
> {<br class="">
> {// Dimension ND regions in plane<br class="">
> Name AxiJacobian;<br class="">
> Case<br class="">
> {<br class="">
> {<br class="">
> Region infiniteBoundaryRegion;<br class="">
> Jacobian<br class="">
> VolAxiSquSphShell{0.95*infiniteBoundaryInnerRadius,<br class="">
> 1.05*infiniteBoundaryOuterRadius};<br class="">
> }<br class="">
> {<br class="">
> Region All;<br class="">
> Jacobian VolAxiSqu;<br class="">
> }<br class="">
> }<br class="">
> }<br class="">
> }<br class="">
><br class="">
> To apply the current density, you need to use time<br class="">
> function that<br class="">
> will be<br class="">
> multiply to your current density. You have to<br class="">
> implement it in<br class="">
> the<br class="">
> Constraint.<br class="">
><br class="">
> Constraint<br class="">
> {<br class="">
> {<br class="">
> Name currentDensityCondition;<br class="">
> Case<br class="">
> {<br class="">
> {<br class="">
> Region inductorRegion;<br class="">
> Value currentDensity[];<br class="">
> TimeFunction timeFunction[];<br class="">
> }<br class="">
> }<br class="">
> }<br class="">
> }<br class="">
><br class="">
> I use the conditional operators to create the time<br class="">
> function:<br class="">
><br class="">
> "timeFunction[] = If condition Then statement Else<br class="">
> statement;"<br class="">
> is then<br class="">
> in conditional operators:<br class="">
> timeFunction[] = (condition) ? statement : statement;<br class="">
><br class="">
><br class="">
> <a href="http://www.c4learn.com/c-programming/c-conditional-operators/" rel="noreferrer" target="_blank" class="">http://www.c4learn.com/c-programming/c-conditional-operators/</a><br class="">
><br class="">
> Good luck,<br class="">
><br class="">
> Frederic<br class="">
><br class="">
> -----Original Message-----<br class="">
> From: Peter Kis <<a href="mailto:md2z34@gmail.com" class="">md2z34@gmail.com</a>><br class="">
> To: <a href="mailto:getdp@geuz.org" class="">getdp@geuz.org</a><br class="">
> Subject: [Getdp] 2D axisym. coil<br class="">
> Date: Sat, 8 Aug 2015 16:57:49 +0200<br class="">
><br class="">
> Dear getDP list members,<br class="">
><br class="">
><br class="">
> I am still struggling with my first model in getDP,<br class="">
> which would<br class="">
> be a<br class="">
> single turn coil on 2D axisymmetric plane.<br class="">
><br class="">
><br class="">
> Now, at least I get something, but I don't like what I<br class="">
> get. I am<br class="">
> not<br class="">
> sure how to apply the time dependent excitation<br class="">
> current density<br class="">
> to the<br class="">
> coil. Therefore I took something over from CoreMassive<br class="">
> model<br class="">
> (found on<br class="">
> the Internet).<br class="">
><br class="">
><br class="">
> As you can see in the attached pictures the magnetic<br class="">
> vector<br class="">
> potential<br class="">
> for the pure 2D model (with Jac Vol) looks good but<br class="">
> for the<br class="">
> axisym model<br class="">
> (with Jac VolAxi) it looks terrible.<br class="">
><br class="">
><br class="">
> Definitely there is still something here what I miss.<br class="">
><br class="">
><br class="">
> Any comment is highly appreciated.<br class="">
><br class="">
><br class="">
> Thanks in advance,<br class="">
><br class="">
><br class="">
> Peter<br class="">
><br class="">
><br class="">
> _______________________________________________<br class="">
> getdp mailing list<br class="">
> <a href="mailto:getdp@geuz.org" class="">getdp@geuz.org</a><br class="">
> <a href="http://www.geuz.org/mailman/listinfo/getdp" rel="noreferrer" target="_blank" class="">http://www.geuz.org/mailman/listinfo/getdp</a><br class="">
><br class="">
><br class="">
><br class="">
><br class="">
><br class="">
><br class="">
><br class="">
</div></div></blockquote></div><br class=""></div>
</div></blockquote></div><br class=""></div></body></html>