Strong DG likely cannot handle unstructured grids which may have flipped, rotated, or oriented faces.

dougshidong / PHiLiP

Parallel High-Order Library for PDEs through hp-adaptive Discontinuous Galerkin methods

Other

45 stars 36 forks source link

Strong DG likely cannot handle unstructured grids which may have flipped, rotated, or oriented faces. #251

Open dougshidong opened 9 months ago

dougshidong commented 9 months ago

Discovery while investigating: https://github.com/dougshidong/PHiLiP/issues/249

Weak DG uses deal.II to obtain correct face projections: https://github.com/dougshidong/PHiLiP/blob/a29ecce53f33a3f0bafc14ef8db767fe967e29b8/src/dg/weak_dg.cpp#L1468-L1501

However strong DG does not and will likely fail for unstructured grids. Attaching a Gmsh grid (had to change suffix from .msh to .log, please change it back before reading into PHiLiP) that could be used to test.

naca0012.log

AlexanderCicchino commented 9 months ago

https://www.dealii.org/developer/doxygen/deal.II/structGeometryInfo.html#a41eabdf96722189263e5ede8ab8d566a Attached this here for reference on the the line/face numering would get changed for each condition.

AlexanderCicchino commented 7 months ago

naca0012.log As an update: there is an error within gmsh reader when using this mesh.