This second part of the work describes the numerical aspects of the developed contact domain method for large deformation frictional contact problems. The theoretical basis of this contact method is detailed in the first part of this work. Starting from this, the present contribution focuses on describing important algorithmic details that go along with the finite element implementation for two-dimensional problems. Important aspects are the construction of the contact domain mesh, via a constraint Delaunay triangulation, the linearization of the discretized contact contributions and some important technical aspects about the extrapolation procedure used for the predictive active set strategy. Finally a set of numerical examples is presented to demonstrate the performance of the developed contact strategy. Demanding static and dynamic contact problems in the context of large deformations, including frictional effects as well as self contact, show the wide applicability and the robustness of the proposed method.