[Getdp] Stranded inductors + circuit relation (3D a-v formulation)

[Zacharie De Grève]

Dear GetDP users and developpers,

I have a short question about the GetDP implementation of the stranded 
inductors model in a 3D a-v formulation, according to reference [1].

More precisely, I am interested in the way the circuit relation (9) and the 
equivalent current density js (section IV) are implemented in GetDP. In the 
case of an explicit current density distribution (which seems simpler to begin 
with), we have : 

js=(Ns/Sc)*t,

where t is a unit vector tangent to the coil direction. My question is the 
following. How could I define the function spaces related to such a quantity ?
I was thinking to something like :

/////////// Function Space ////////////////////

// a
...

// Stranded inductors
  { Name Hregion_i_Mag_3D ; Type Vector ;
    BasisFunction {
      { Name ss ; NameOfCoef is ; Function ?????? ;
        Support DomainS_Mag ; Entity DomainS_Mag ; }
    }
    GlobalQuantity {
      { Name Is ; Type AliasOf        ; NameOfCoef is ; }
      { Name Us ; Type AssociatedWith ; NameOfCoef is ; }
    }
    Constraint {
      { NameOfCoef Us ; EntityType Region; NameOfConstraint Voltage_3D ;}
      { NameOfCoef Is ; EntityType Region; NameOfConstraint Current_3D ;}
    }
  }
  
  // Uz et Iz for circuit relations 
  { Name Hregion_Cir ; Type Scalar ;
    BasisFunction {
      { Name sr ; NameOfCoef ir ; Function BF_Region ;
        Support DomainZtot_Cir ; Entity DomainZtot_Cir ; }
    }
    GlobalQuantity {
      { Name Iz ; Type AliasOf        ; NameOfCoef ir ; }
      { Name Uz ; Type AssociatedWith ; NameOfCoef ir ; }
    }
    Constraint {
      { NameOfCoef Uz ; EntityType Region ; NameOfConstraint VoltageInit_Cir ; 
}
      { NameOfCoef Uz ; EntityType Region ; NameOfConstraint Voltage_Cir ; }
      { NameOfCoef Iz ; EntityType Region ; NameOfConstraint Current_Cir ; }
    }
  }


Then, it would be easy to write the terms related to the circuit relation in 
the formulation object.

Thank you very much in advance for your help.

Best regards,


Zacharie De Grève
UMons


[1] P. Dular, N. Sadowski, J.P.A. Bastos and W. Legros, Dual Complete 
Procedures to Take Stranded Inductors into Account in Magnetic Vector 
Potential Formulations, IEEE Trans Mag, vol 36(4), July 2000


-- 
---------------------------------------------------
Ir Zacharie De Grève
PhD student - FNRS Research Fellow
Dept. of Power Electrical Engineering
University of Mons - Faculty of Engineering

Bd Dolez, 31
BE - 7000 Mons
Belgique
Tél : +32 65 374125
---------------------------------------------------
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/getdp/attachments/20110321/c55d7399/attachment.html>