Hi Sebastian FE_NedelecSZ does not support hanging nodes, i.e. the mesh must be conforming.
Implementing handing-node constraints for this element is possible, but would require a special approach since constraints are not static. If you feel like implementing this, I would be glad to discuss -- I happen to miss this functionality too, but I've never come to a point where I can devote enough time to this. Best Alexander On Wednesday, September 22, 2021 at 10:37:08 AM UTC+2 [email protected] wrote: > Thank you for the first answer! > > When I implement the computation of the embedding matrix , as done in > FE_Nedelec, in FE_NedelecSZ the program compiles and runs without any > warnings or errors (even in debug mode). > But from the results it appears that there still is a problem with the > hanging nodes. > > Maybe it is the best when I attach some pictures of the solution. For the > test case I computed a reference solution to which I know the analytical > solution. > All faces in that test case are oriented (therefore there should be no > sign conflict) moreover, it works when I use FE_Nedelec (see appendix). > With just global refinement also FE_NedelecSZ works fine. > > So my best guess is, that there is still a problem with the hanging nodes. > (Are there maybe more functions that need to be called for the hanging > nodes to work correctly? At the moment I am only define > this->interface_constraints in FE_NedelecSZ.) > > Best, > Sebastian > > [image: FE_NedelecSZ.png][image: FE_Nedelec.png] > > Wolfgang Bangerth schrieb am Dienstag, 21. September 2021 um 18:55:57 > UTC+2: > >> On 9/21/21 8:35 AM, Sebastian Kinnewig wrote: >> > >> > To make the switch to FE_NedelecSZ I need to compute the face >> embedding. >> > I tried to compute the face embedding via the function >> > FETools::compute_face_embedding_matricies() and pass the resulting face >> > embedding matrix to the interface_constraints (similar to what is done >> > in the FE_Nedelec class). But that doesn't seem to work, does anyone >> > have any ideas or suggestions as to what the problem could be? >> >> Can you describe what happens (i.e., what "that doesn't seem to work" >> actually means)? >> >> Best >> W. >> >> >> -- >> ------------------------------------------------------------------------ >> Wolfgang Bangerth email: [email protected] >> www: http://www.math.colostate.edu/~bangerth/ >> > -- 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]. To view this discussion on the web visit https://groups.google.com/d/msgid/dealii/9ec3439d-bcec-4514-a7bc-cb9c148402e0n%40googlegroups.com.
