Skip to content

Draft: Resolve "Magnetic-mechanic coupling with non-conforming interfaces in 3D"

Vladimir Filkin requested to merge nonconforming_magnetic_issue into master

cfs!57

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

Merge request reports