In this work a coupled (LDG-CFEM) method for singularly perturbed Volterra integro-differential equations with a smooth kernel is implemented. The existence and uniqueness of the coupled solution is given, provided that the source function and the kernel function are sufficiently smooth. Furthermore, the coupled solution achieves the optimal convergence rate p + 1 in the L 2 norm and a superconvergence rate 2p at nodes for the numerical solutionÛ N with the one-sided flux inside the boundary layer region under layer-adapted meshes uniformly with respect to the singular perturbation parameter .