Structured Mesh Element Orientation #15
Description
For structured mesh portions, the element orientation theta_11 does not appear to always be calculated correctly for triangular elements. This issue comes up for the first/last element of each layer with the rounded edges. Also, potentially a problem for sharp edge refinement of quad elements to triangles (e.g., sharp trailing edges).
Notably, normal quad elements placed around the perimeter appear to be handled correctly and before the elements get split to triangles.
For debugging can add the line:
print(len(self.nodes), np.degrees(np.arctan2(v1.Y(), v1.X())))
Running the example for the 15 MW at the 0.7 station appears to show the behavior by the triangular elements having notably different angles than neighboring triangles/quads. Specifically for meshing segment 3, layer 2. Also segment 0, layer 2 appears to shows this.