In this thesis, the validity of the adopted MRT lattice Boltzmann model [14] is examined by computing two-dimensional Poiseuille flow, lid-driven cavity flow, and three-dimensional Poiseuille flow in a square duct. We also use the density fluctuation and assume the mean density to reduce effects due to compressibility [48]. The present simulation indicates that SRT method fails to simulate higher Reynolds number flows due to the limitation of tau value. In contrast to the SRT method, the tau value for MRT method can be close to 0.5 for 2-D and 3-D simulations owing to the different relaxation rates. It should be noted that 3D MRT computations are very sensitive to boundaries, especially the treatment of the corners and edges. As adopting smaller value of tau combined with MRT method, higher Reynolds number flow can be realized. Finally, parallel-MRT and parallel-SRT models are implemented for large grid density simulations, and simulated results are presented.