[Getdp] 2D axisym. coil
Peter Kis
md2z34 at gmail.com
Sun Aug 9 23:28:38 CEST 2015
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> 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> 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 --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/getdp/attachments/20150809/fd479b13/attachment.html>