[Getdp] Dual Formulation in the wiki
David Colignon
David.Colignon at ulg.ac.be
Wed Jan 31 14:11:21 CET 2007
Hi,
from http://geuz.org/pipermail/getdp/2006/000869.html in the mailing list:
...
the condition on the coefficients associated to the second order basis functions must be homegeneous ( 0 value ) otherwise the imposed condition considers a boundary condition other than constant ( with second degree variation ). This is due to the fact that the second order elements are built with hierarchical basis functions.
...
Cheers,
Dave
--
David Colignon, Ph.D.
ELAP - Service d'Electricité Appliquée
Institut Montefiore B28
Université de Liège
4000 Liège - BELGIQUE
Tél: +32 (0)4 366 37 32
Fax: +32 (0)4 366 29 10
http://elap.montefiore.ulg.ac.be
Amit.Itagi at seagate.com wrote:
> Hi,
>
> I have a question regarding the dual formulation example (
> http://www.geuz.org/getdp/wiki/DualFormulations ). I find that the first
> order and the higher order elements satisfy different constraints (See
> excerpt below). Could someone kindly explain the motivation for defining
> separate constraints ?
>
> The higher order elements are being used as a result of P1 adaptation. I
> would think that since the same problem is being solved, the constraints
> would remain the same.
>
> Thanks
>
> Rgds,
> Amit
>
>
> //----------------------------------------------------------------------------------------------------------------------------
> FunctionSpace {
> { Name Hgrad_phi_MagSta; Type Form0;
> BasisFunction {
> { Name sn1n ; NameOfCoef wn1n ; Function BF_Node ; Support Tot_Mag;
> Entity NodesOf[All]; }
> { Name sn2e ; NameOfCoef wn2e ; Function BF_Node_2E ; Support Tot_Mag
> ; Entity EdgesOf[ All ] ; }
> .
> .
> .
> }
> Constraint {
> { NameOfCoef wn1n; EntityType NodesOf ; NameOfConstraint Sta_phi; }
> { NameOfCoef wn2e ; EntityType EdgesOf ; NameOfConstraint
> Sta_phi_Zero ; }
> .
> .
> .
> }
> }
>
> Constraint {
> { Name Sta_phi ;
> Case {
> { Region left ; Value -100. ; }
> { Region right ; Value 100. ; }
> }
> }
> { Name Sta_phi_Zero ;
> Case {
> { Region left ; Value 0. ; }
> { Region right ; Value 0. ; }
> }
> }
> //---------------------------------------------------------------------------------
>
>
>
> _______________________________________________
> getdp mailing list
> getdp at geuz.org
> http://www.geuz.org/mailman/listinfo/getdp