This work presents some numerical aspects of isogeometric boundary element methods (IGABEM). The behavior of hyper-singular and nearly-singular integration is first explored on the distorted NURBS surface. Several numerical treatments are proposed to enhance the quadrature in the framework of isogeometric analysis. Then a numerical implementation of IGABEM on the trimmed NURBS is detailed. Based on this idea, the surface crack problem is modeled incorporation with the phantom element method. The proposed method allows the crack to intersect with the boundary of the body while preserving the original parametrization of the NURBS-based CAD geometry.