A discontinuous Galerkin finite element method is proposed for solving batch crystallization models. The suggested method has the capability of capturing sharp discontinuities and narrow peaks of the crystal size distribution (CSD). The accuracy of the method can be improved by introducing additional nodes in the same solution element and, hence, avoids the expansion of mesh stencils which is normally observed in the high order finite volume schemes. For that reason, the method can be easily applied up to boundary cells without losing accuracy. The method is robust and well suited for large-scale time-dependent computations in which a high degree of accuracy is demanded. Several test cases are carried out in this paper. The numerical results verify the efficiency and accuracy of the proposed method.