Intrinsic remeshing preserves the original geometry
Numerical quadrature estimates the contact area
Meshless 3D restricted voronoi diagram discretizes the volume
Heat visualized on each intrinsic triangle and flux visualized on the surface