non delaunay mesh elements

[tcaduser]

devsim assumes all elements are delaunay. A bad mesh element can result in poor conservation of current, because of the large current cross section introduced. This appears to be handled in Genius using a truncation scheme in 2d:
https://github.com/cogenda/Genius-TCAD-Open/blob/master/src/geom/face_tri3_fvm.cc

For 3d, a similar approach may be used. It would be much more complicated since it would require handling obtuse triangles on the tetrahedron faces. It would also have to handle the center of the circumsphere being outside of the tetrahedral element.

[tcaduser]

a nice description of the problem in 2d is here:
http://www.iue.tuwien.ac.at/phd/klima/node21.html

[tcaduser]

http://www.iue.tuwien.ac.at/phd/fleischmann/

[tcaduser]

also consider an alternative like some kind of weighted center as opposed to truncation