Modeling of impurity effect on drift instabilities in plasmas with many ion species