[Getdp] 2D axisym. coil

Frederic Trillaud Pighi ftrillaudp at pumas.iingen.unam.mx
Mon Aug 10 19:04:29 CEST 2015


Thanks to point out this option.

Best,

Frederic

> On Aug 10, 2015, at 11:55 AM, David Colignon <david.colignon at ulg.ac.be> wrote:
> 
> 
> On 10/08/15 17:14, Frederic Trillaud Pighi wrote:
> 
>> Dear Peter,
>> 
>> 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.
>> 
>> 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.
> 
> 
> Hi,
> 
> maybe should you have a look at "global quantities" in the manual ?
> 
> Regards,
> 
> Dave
> 
> -- 
> David Colignon, Ph.D.
> 1er logisticien de recherche
> Université de Liège
> ACE - Applied & Computational Electromagnetics
> Sart-Tilman B28
> 10, Grande Traverse
> 4000 Liège - BELGIQUE
> Tél: +32 (0)4 366 37 32
> http://www.ulg.ac.be/nic4
>> 
>> Best,
>> 
>> Frederic
>> 
>> 
>> 
>>> On Aug 9, 2015, at 4:28 PM, Peter Kis <md2z34 at gmail.com <mailto:md2z34 at gmail.com>> wrote:
>>> 
>>> Dear Frederic,
>>> 
>>> 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?
>>> 
>>> I compared the results given by getdp with COMSOL and it's a perfect match! I am very pleased.
>>> 
>>> 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?
>>> 
>>> Thanks again,
>>> 
>>> Kind Regards
>>> 
>>> 
>>> Peter
>>> 
>>> 
>>> On Sun, Aug 9, 2015 at 1:52 AM, Frederic Trillaud Pighi <ftrillaudp at pumas.iingen.unam.mx
>>> <mailto:ftrillaudp at pumas.iingen.unam.mx>> wrote:
>>> 
>>>    Dear Peter,
>>> 
>>>    Here it goes. I made some changes (transformation with an outer shell,
>>>    killed the initial time value of A, you should not need it) and did some
>>>    tests so the parameters of current density etc. have been changed. It
>>>    should work fine now.
>>> 
>>>    Best,
>>> 
>>>    Frederic
>>> 
>>> 
>>>    On Sun, 2015-08-09 at 00:27 +0200, Peter Kis wrote:
>>>    > I am sending the model files. Thanks for your time.
>>>    >
>>>    >
>>>    > Peter
>>>    >
>>>    >
>>>    > On Sat, Aug 8, 2015 at 11:02 PM, Frederic Trillaud
>>>    > <ftrillaudp at pumas.ii.unam.mx <mailto:ftrillaudp at pumas.ii.unam.mx>> wrote:
>>>    >         Send me back your files, I will take a look at it over the
>>>    >         weekend to
>>>    >         see if I can pin the problem.
>>>    >
>>>    >         Frederic
>>>    >
>>>    >
>>>    >         -----Original Message-----
>>>    >         From: Peter Kis <md2z34 at gmail.com <mailto:md2z34 at gmail.com>>
>>>    >         To: Frederic Trillaud <ftrillaudp at pumas.ii.unam.mx <mailto:ftrillaudp at pumas.ii.unam.mx>>
>>>    >         Cc: getdp at geuz.org <mailto:getdp at geuz.org>
>>>    >         Subject: Re: [Getdp] 2D axisym. coil
>>>    >         Date: Sat, 8 Aug 2015 22:43:11 +0200
>>>    >
>>>    >         Dear Frederic,
>>>    >
>>>    >
>>>    >         thanks for your reply. There must be something else here,
>>>    >         because
>>>    >         VolAxiSqu hasn't helped.
>>>    >
>>>    >
>>>    >         I think the time function in my code might be okay, because it
>>>    >         seems to
>>>    >         work fine for Jacobi Vol.
>>>    >
>>>    >
>>>    >         Cheers,
>>>    >
>>>    >
>>>    >         Peter
>>>    >
>>>    >
>>>    >         On Sat, Aug 8, 2015 at 5:51 PM, Frederic Trillaud
>>>    >         <ftrillaudp at pumas.ii.unam.mx <mailto:ftrillaudp at pumas.ii.unam.mx>> wrote:
>>>    >                 Dear Peter,
>>>    >
>>>    >                 I had a similar problem, you should try VolAxiSqu. It
>>>    >         solves my
>>>    >                 issue on
>>>    >                 the axis of symmetry.
>>>    >
>>>    >                 Jacobian
>>>    >                 {
>>>    >                  {// Dimension ND regions in plane
>>>    >                   Name AxiJacobian;
>>>    >                   Case
>>>    >                   {
>>>    >                    {
>>>    >                     Region infiniteBoundaryRegion;
>>>    >                     Jacobian
>>>    >         VolAxiSquSphShell{0.95*infiniteBoundaryInnerRadius,
>>>    >                 1.05*infiniteBoundaryOuterRadius};
>>>    >                    }
>>>    >                    {
>>>    >                     Region All;
>>>    >                     Jacobian VolAxiSqu;
>>>    >                    }
>>>    >                   }
>>>    >                  }
>>>    >                 }
>>>    >
>>>    >                 To apply the current density, you need to use time
>>>    >         function that
>>>    >                 will be
>>>    >                 multiply to your current density. You have to
>>>    >         implement it in
>>>    >                 the
>>>    >                 Constraint.
>>>    >
>>>    >                 Constraint
>>>    >                 {
>>>    >                 {
>>>    >                   Name currentDensityCondition;
>>>    >                   Case
>>>    >                   {
>>>    >                    {
>>>    >                     Region inductorRegion;
>>>    >                     Value currentDensity[];
>>>    >                     TimeFunction timeFunction[];
>>>    >                    }
>>>    >                   }
>>>    >                  }
>>>    >                 }
>>>    >
>>>    >                 I use the conditional operators to create the time
>>>    >         function:
>>>    >
>>>    >                 "timeFunction[] = If condition Then statement Else
>>>    >         statement;"
>>>    >                 is then
>>>    >                 in conditional operators:
>>>    >                 timeFunction[] = (condition) ? statement : statement;
>>>    >
>>>    >
>>>    > http://www.c4learn.com/c-programming/c-conditional-operators/
>>>    >
>>>    >                 Good luck,
>>>    >
>>>    >                 Frederic
>>>    >
>>>    >                 -----Original Message-----
>>>    >                 From: Peter Kis <md2z34 at gmail.com <mailto:md2z34 at gmail.com>>
>>>    >                 To: getdp at geuz.org <mailto:getdp at geuz.org>
>>>    >                 Subject: [Getdp] 2D axisym. coil
>>>    >                 Date: Sat, 8 Aug 2015 16:57:49 +0200
>>>    >
>>>    >                 Dear getDP list members,
>>>    >
>>>    >
>>>    >                 I am still struggling with my first model in getDP,
>>>    >         which would
>>>    >                 be a
>>>    >                 single turn coil on 2D axisymmetric plane.
>>>    >
>>>    >
>>>    >                 Now, at least I get something, but I don't like what I
>>>    >         get. I am
>>>    >                 not
>>>    >                 sure how to apply the time dependent excitation
>>>    >         current density
>>>    >                 to the
>>>    >                 coil. Therefore I took something over from CoreMassive
>>>    >         model
>>>    >                 (found on
>>>    >                 the Internet).
>>>    >
>>>    >
>>>    >                 As you can see in the attached pictures the magnetic
>>>    >         vector
>>>    >                 potential
>>>    >                 for the pure 2D model (with Jac Vol) looks good but
>>>    >         for the
>>>    >                 axisym model
>>>    >                 (with Jac VolAxi) it looks terrible.
>>>    >
>>>    >
>>>    >                 Definitely there is still something here what I miss.
>>>    >
>>>    >
>>>    >                 Any comment is highly appreciated.
>>>    >
>>>    >
>>>    >                 Thanks in advance,
>>>    >
>>>    >
>>>    >                 Peter
>>>    >
>>>    >
>>>    >                 _______________________________________________
>>>    >                 getdp mailing list
>>>    > getdp at geuz.org <mailto:getdp at geuz.org>
>>>    > http://www.geuz.org/mailman/listinfo/getdp
>>>    >
>>>    >
>>>    >
>>>    >
>>>    >
>>>    >
>>>    >
>>> 
>>> 
>> 
>> 
>> 
>> _______________________________________________
>> getdp mailing list
>> getdp at geuz.org
>> http://www.geuz.org/mailman/listinfo/getdp
>> 
> 
> _______________________________________________
> getdp mailing list
> getdp at geuz.org
> http://www.geuz.org/mailman/listinfo/getdp