[Getdp] function space for 3D a-v-formulation, and higher order approximation ?

Christophe Geuzaine cag32 at case.edu
Sun May 14 01:14:57 CEST 2006


Thomas Jung wrote:
> Hi,
> 
> I would be grateful if someone could give me a little help on this:
> 
> First, would the following be a correct declaration of a function space for 
> the 3D magnetic vector potential ?
> 
> FunctionSpace {
> 
>   // Magnetic vector potential a (b = curl a)
>   { Name Hcurl_a_Mag_3D; Type Form1;
>     BasisFunction {
> 	{ Name se; NameOfCoef ae; Function BF_Edge;
>         Support Domain_Mag; Entity EdgesOf[ All ]; }
>     }
>     Constraint {
> 	{ NameOfCoef ae; EntityType EdgesOf;NameOfConstraint a; }
> 	{ NameOfCoef ae; EntityType EdgesOfTreeIn; EntitySubType 
> StartingOn;NameOfConstraint Gauge; }
>     }
>   }
> 
> }
> 
> with the Gauge constraint like this:
> 
>   { Name Gauge; Type Assign;
>   Case {
> 	{ Region Domain_Mag ; SubRegion OuterBoundary; Value 0.; }
>     }
>   }


Looks good to me.

> 
> 
> If that would be o.k., how would one do a higher order approximation?
> 

You need to add the higher order basis functions in the function space: 
BF_Edge_2E, BF_Edge_2F, etc.


> 
> Thank you very much !
> 
> 
> Thomas Jung
> _______________________________________________
> getdp mailing list
> getdp at geuz.org
> http://www.geuz.org/mailman/listinfo/getdp
> 


-- 
Christophe Geuzaine
Assistant Professor, Case Western Reserve University, Mathematics
http://www.case.edu/artsci/math/geuzaine