[Getdp] 2D axisym. coil

[Frederic Trillaud]

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