add support of optimal high-order discrete velocity set

[zhulianhua]

Used optimized discrete velocity for 2D and 3D flows, instead of using tensor product of 1D Gauss-Hermite discrete velocity to construct for 2D&3D cases. Add option to use discrete equilibrium distribution function to save computing time and improve conservation property of collision term.