Hi colleagues,
When testing my codes with FENedelec and FENedelecSZ elements, I saw that
the shape functions change versus the mesh size.
I think that for each edge of a cell, the shape functions for the degree of
freedom related to that edge are scaled with the inverse edge length.
For example, assuming I have a cell that is a cube, whose edge length is 1,
the shape function at the quadrature point located on the edge is {0, 0, 1}.
For the smaller cell whose edge length is 0.5, the shape function values is
{0, 0, 2}.
In my limited opinion, the shape function is defined in unit cell over
[0,1]^3, and the shape function should not change when the real cell change
its size.
This problem leads to the change of my solution when using the finer mesh
compared to that of coarser one.
Because both my cell matrix and cell rhs change versus the edge length.
For FeNedelecSZ, I tested with a cell limited in the region [0,1]^3, which
means that the edge length is 1 everywhere.
When printing out the integral over the cell for one degree of freedom by
\int_Omega \varphi_i(x) Js(x) dx
For simplicity, I set Omega =1 and Js(x) = {1,1,1} for testing the integral
over the cell at degree of freedom i_th when source magnitude is 1.
The result for this is 0.25 rather than 1.
Is this means that for a source of magnitude 1., I actually can model it as
a source edge at degree of freedom i_th with the magnitude 0.25?
I think there should be some idea that we change the shape function values
like this.
Could you please give me some advice?
Thank you very much for your great support.
Best regards,
--
The deal.II project is located at http://www.dealii.org/
For mailing list/forum options, see
https://groups.google.com/d/forum/dealii?hl=en
---
You received this message because you are subscribed to the Google Groups
"deal.II User Group" group.
To unsubscribe from this group and stop receiving emails from it, send an email
to [email protected].
For more options, visit https://groups.google.com/d/optout.
