Simple 2D thermal solver on regular mesh