Testifying the divergence free property of the pillar free torsion boundary element method integrand, warping functions can be converted into the numerical evaluation of torsion potential functions at the boundary points and the singular numerical integrals are not needed at all. Numerical results are present to show the excellent accuracy of the new method.