Creation of ordered cuts (cohomology-homology) in 2D
Dear all,
I have created the same issue in Models for Gmsh and GetDP. I thought that it may be more relevant here after all.
I have a 2D axisymmetric case of a coil made of 10 conductors in series. The coil is part of a small electric circuit (power supply, resistance in series with the coil). The same current +I is impressed to the series-connected conductors through cuts associated with branches in the circuit, from node 2 to 11.
{ Region cut_0; Branch{2, 3}; }
{ Region cut_1; Branch{3, 4}; }
{ Region cut_2; Branch{4, 5}; }
{ Region cut_3; Branch{5, 6}; }
{ Region cut_4; Branch{6, 7}; }
{ Region cut_5; Branch{7, 8}; }
{ Region cut_6; Branch{8, 9}; }
{ Region cut_7; Branch{9, 10}; }
{ Region cut_8; Branch{10, 11}; }
{ Region cut_9; Branch{11, 0}; }
The issue is that the cuts are not in the same direction for all the branches and it is somehow random leading to a wrong direction for the current in each branch.
I know that a colleague of mine made a manual fixed by reordering the edges in the cuts saving the mesh in ascii format. I think that it would be more appropriate to have it built in the source code to avoid having to reorder the cuts or to give manually a particular sign to the current or even the direction of the branch for the matter after the creation of the mesh and the cuts.
I am joining a picture of some of the cuts and the current density resulting from the basic code without manual ordering to point out the issue.
Best,
Frederic