LBM is a recent technique that has been shown to be as accurate as traditional CFD methods. Furthermore LBM is able to integrate arbitrarily complex geometries at a reduced computational cost. LBM simplifies Boltzmann’s original formulation by discretizing both time and space. In this discretization, fluid particle positions are confined to the node of a lattice. Instead of considering the entire continua of velocity directions and magnitudes, and varying particle mass, variations of fluid particle momenta are discretized into several finite directions. In addition to a fluid solver such as LBM, IBM has been widely used in the simulation of blood flow. In this method, the effect of a moving boundary is accounted for by a force density that is distributed to the Cartesian mesh in the vicinity of the moving boundary. The force density is calculated from the boundarys constitutive law. IBM is suitable for modeling fluid-structure interaction in blood flow simulations where plasma and RBC mechanics are solved in a partitioned approach.