Hi Shahab, My best guess is that the attached manifold is not consistent to what you need after the grid transformations. You have to be aware that the attached manifold is used for generating new vertices (as it is needed during AMR/load-balancing in regions that is newly owned by a rank). Maybe you could attach a new `CylindricalManifold` via `tria.set_manifold(0, CylindricalManifold<3>(/*TODO: appropriate arguments*/));` (with arguments consistent to the transformation). Hopefully, this fixes your problem!
PM On Monday, 2 August 2021 at 21:33:59 UTC+2 [email protected] wrote: > Shahab. > > You can loop over the cells, get the vertices (see > https://dealii.org/current/doxygen/deal.II/classTriaAccessor.html#a3dd6518eb0cf5fccc5926470128415d9) > > and print the Point associated to each vertex. Then you can compare the > vertices before and after load balancing. > > Best, > > Bruno > > On Monday, August 2, 2021 at 1:12:09 PM UTC-4 [email protected] wrote: > >> Dear Bruno, >> >> Thank you for your reply. >> I am 100% sure that this is not a rendering issue, because in my solver >> the particles that are located on top of this boundary are influenced by >> the deformation of the cells. >> How can I exactly quantify the deformation? >> >> Best, >> Shahab >> On Thursday, July 29, 2021 at 2:50:09 PM UTC-4 [email protected] >> wrote: >> >>> Shahab, >>> >>> Can you quantify the deformation? It would be useful if you could output >>> the coordinates of the vertices of the deformed cells before and after load >>> balancing. It is possible that's a rendering issue but it's hard to tell >>> without having numbers. >>> >>> Best, >>> >>> Bruno >>> >>> On Wednesday, July 28, 2021 at 6:48:50 PM UTC-4 [email protected] >>> wrote: >>> >>>> Dear All, >>>> >>>> I have noticed that some moving grids get deformed after >>>> load-balancing. I should mention that this only happens when I move the >>>> grid using GridTools::rotate or GridTools::shift, and only for curved >>>> geometries. For instance, it doesn't happen in a box, but I've observed >>>> this problem in cylinders. Furthermore, only the cells which are exchanged >>>> between the processors during load-balancing get deformed. I've attached >>>> an >>>> image of this problem. It shows a rotating cylinder. On the left, it shows >>>> the triangulation output before load-balancing. On the right, you can see >>>> the triangulation after load-balancing and how the exchanged cells are >>>> deformed (inside the red circle). >>>> >>>> I would appreciate it if you give me any suggestions about the reason >>>> for this problem. Thank you in advance. >>>> >>>> Best regards, >>>> Shahab >>>> >>> -- 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/86cbabe2-57f2-4474-aa9c-966aa8cc191cn%40googlegroups.com.
