Numerical schemes for multi-species BGK equations based on a variational procedure

applied to multi-species BGK equations with velocity-dependent collision frequency and to quantum multi-species BGK equations