Numerical schemes for the BGK kinetic model