Numerical schemes for a multi-species BGK model with velocity-dependent collision frequency