The application of the lattice Boltzmann method to the large vessel bifurcation blood flow is investigated in awide range of Reynolds numbers. The velocity, shear stress and pressure distributions at the bifurcation arepresented in detail. The flow separation zones revealed with increase of Reynolds number are located in theareas of the daughter branches distal to the outer coers of the bifurcation where some deposition of particularblood components might occur to form arteriosclerosis. The results also demonstrate that the lattice Boltzmannmethod is adaptive to simulating the flow in larger vessels under a high Reynolds number.