[Getdp] Magnetodynamic 3D - Electromagnetic circuit element

Florin CIUPRINA florin at elmat.pub.ro
Mon Apr 14 10:59:40 CEST 2003


Hello Christophe,

Thank you for your help!

>A correct definition of the discrete vector potential involves the edges
>of the _complement_ of the spanning tree. The easiest way to do this
>is to define the vector potential everywhere, and use the constraint to
>fix the dofs to zero on the edges of the tree (so that only the edges of
>the complement lead to actual unknowns).

If I understand well, this means:

FunctionSpace {
   { Name Hcurl_a_Gauge; Type Form1;
     BasisFunction {
      // a = a  s
       //     e  e
       { Name se; NameOfCoef ae; Function BF_Edge;
         Support DomainC_Mag ; Entity EdgesOf[All]; }
     }
     Constraint {
      //gauge condition ae = 0 on tree edges and on some boundaries
       { NameOfCoef ae;
         EntityType EdgesOfTreeIn; EntitySubType StartingOn;
         NameOfConstraint GaugeCondition_a_Mag_3D; }
     }
   }

Constraint {
{ Name GaugeCondition_a_Mag_3D; Type Assign;
     Case {
       { Region DomainC_Mag ; SubRegion Surf; Value 0.; }
     }
}
Is this correct? If yes, probably something else is wrong because
unfortunately with this the problem does not converge:

GetDP :  998  1.0825175e+06  9.4973969e-01
GetDP :  999  1.0559604e+06  9.2643995e-01
GetDP : Warning   : Iterative solver has iterated too many times
GetDP : 1000  1.0548821e+06  9.2549394e-01
GetDP : Solver    : 1000 Iterations / Residual: 1.21005e+06

What can I do?

Florin



Florin Ciuprina, Ph.D.
Associate Professor
----------------------------------------------------------------------------
POLITEHNICA University of Bucharest
Electrical Engineering Faculty
Electrical Materials Laboratory
Spl.Independentei 313, 77206, Bucharest, Romania
Tel:  +40.21.402 92 91
Fax: +40.21.410 43 55
e-mail: florin at elmat.pub.ro
---------------------------------------------------------------------------- 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://www.geuz.org/pipermail/getdp/attachments/20030414/edef1f1e/attachment.html>