[Getdp] 2D axisym. coil

[Frederic Trillaud Pighi]

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> 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>
>         To: Frederic Trillaud <ftrillaudp at pumas.ii.unam.mx>
>         Cc: 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> 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>
>                 To: 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
>                 http://www.geuz.org/mailman/listinfo/getdp
>         
>         
>         
>         
>         
> 
> 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: coil.geo
Type: text/x-csrc
Size: 1987 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20150808/723670bb/attachment.geo>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: coil.pro
Type: text/x-csrc
Size: 6026 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20150808/723670bb/attachment.pro>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: param.geo
Type: text/x-csrc
Size: 259 bytes
Desc: not available
URL: <http://www.geuz.org/pipermail/getdp/attachments/20150808/723670bb/attachment-0001.geo>