Draft: Resolve "Magnetic-mechanic coupling with non-conforming interfaces in 3D"
The problem
Non-conforming interfaces significantly increase RAM consumption and CPU time for some models.
We know
- Non-conforming interface introduces zeros as non-zero entities into a system matrix.
- Gmesh- and Cubit-based models seem to have different behaviour (NOT SURE, not easy to build exactly the same non-hex meshes)
- Each DoF of master surface is connected with each dof of slave surface (and vice versa). If there is no real connection, we have a zero non-zero element in the system matrix.
- For some cases (significant increase of matrix density) we have a huge increase in computational time (solver) and memory consumption.
Open questions
-
Check how a non-conforming interface affects the number of non-zero entities. For Cubit and gmsh. -
Only for Cubit go back to only hex elements (3D edge) - DONE!
-
Investigate different intersections and how many nodes and non-zero entities we have. What is the pattern? DONE! (see jupyter file fore comments and patterns (separate folder))
-
Check the same for gmsh. Does it work in the same way? NOT DONE
-
Check if zero non-zero entities affect the performance of the solvers significantly. (tested with scipy and pypardiso) -
Implement changes?
Edited by Vladimir Filkin