A new method for the numerical evaluation of domain integrals in a 3D boundary element method for transient heat conduction