[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>