Hi Torquil, I have an element which is not an answer to your issue but it can help.I never used the RT0 elements but I used the Nedelec elements and the shape function associated to edge {ij} is :
1/length({ij})*(lambda_i grad lambda_j - lambda_j grad lambda_i).Maybe something similar is done for RT0 elements (for instance a dof associated to a face is not the flux through this face but the mean value of the normal component of the vector field on this face).
Best regards, Ronan Le 23/09/2010 20:01, Torquil Macdonald Sørensen a écrit :
Hi!
I'm having trouble with calculation of mass matrix elements on a 3d tetrahedron
mesh. I have two different cases. In the first case I want to use Whitney 1-form
edge elements, and in the second case Whitney 2-form face elements. I'm given to
understand that these are the same as the NEDELEC and RT0 elements offered by
GetFEM.
The Whitney elements are defined as follows: For a point i, lambda_i is the
usual affine function that is equal to 1 at i, and equal to zero at neighbouring
points.
For an edge {ij} going from point i to point j, the Whitney 1-form element is
lambda_ij := lambda_i grad lambda_j - lambda_j grad lambda_i
The "2-form" element for a face {ijk} is
lambda_{ijk} := 0.5*(lambda_i (grad lambda_j) x (grad lambda_k) + cyclic
permutations of {ijk}).
When varying the "fineness"/lattice spacing of the mesh for a given physical
domain [0,1]^3, these definitions lead to the following behaviour of the mass
matrix components:
m_ff' ~ 1/h
m_ee' ~ h
where m_ff' is a component of the mass matrix for the 2-form elements, m_ee' is
the same for the 1-form elements, and h is the lattice spacing.
However, in my case, the components m_ff' decreases when h decreases. Does this
indicate a difference between the definition of the RT0 element relative to the
Whitney 2-form element? If so, is it possible to get one from the other by a
simple multiplication by h?
Thanks
Torquil
_______________________________________________
Getfem-users mailing list
[email protected]
https://mail.gna.org/listinfo/getfem-users
<<attachment: ronan_perrussel.vcf>>
_______________________________________________ Getfem-users mailing list [email protected] https://mail.gna.org/listinfo/getfem-users
