A low diffusive flux difference splitting based kinetic scheme is developed on a discrete velocity Boltzmann equation, with novel three model. While two velocities are used for upwinding, the third utilized to introduce appropriate additional numerical diffusion only in expansion regions, identified using relative entropy (Kullback–Liebler divergence) at cell-interface, along estimation of phys...