In this article we propose a finite-volume discretization of a one-dimensional nonlinear reaction kinetic model proposed in Neumann & Schmeiser (2016), which describes a two-species recombination-generation process. Specifically, we establish the long-time convergence of approximate solutions towards equilibrium, at exponential rate. The study is based on an adaptation for a discretization of the linearized problem of the $L^{2}$ hypocoercivity method introduced in Dolbeault et al. (2015). From this we can deduce a local result for the discrete nonlinear problem, in the sense that small initial perturbations from the steady state are considered. As in the continuous framework this result requires the establishment of a maximum principle, which necessitates the use of monotone numerical fluxes.