Numerical approximation of a BGK-type relaxation model for reactive mixtures