The Gaussian-BGK model of Boltzmann equation with small Prandtl number