In this thesis, implementation of the Lattice Boltzmann method for simulation of non-Newtonian blood fluid flows is considered. Several non-Newtonian models are introduced and analytic solutions, where possible, are derived. Among these non-Newtonian models, those which have been used for blood flow are considered in particular. Furthermore, some abnormal characteristics of blood flow and effects of blood cells on blood flow properties are explained. In the present work, numerical results of lattice Boltzmann simulation of Newtonian fluid flow between two parallel plates are presented. Several boundary conditions namely periodic, entrance velocity and constant pressure are used and mid grid bounce back boundary condition is implemented on the solid walls. In addition, momentum equation is derived from the lattice Boltzmann equation for non-Newtonian fluids. Lattice Boltzmann simulations of five non-Newtonian models are presented and in cases which analytic solution exists; analytic solutions are compared with numerical results. Predicted velocity profiles by the non-Newtonian LBM for a channel flow have an excellent agreement with those of exact solutions, presenting the capability of the LB method for the non-Newtonian fluid flow simulation. Finally, the procedure of parallel programming using LBM is explained and Parallel Distribution toolbox of Matlab software is introduced.