Motion of leukocytes in a pulsatile blood flow through a stenotic artery is numerically simulated by the Eulerian-Lagrangian approach. Blood is modeled by a dilute particle-suspending liquid, the particle phase representing the leukocytes and the liquid phase representing the blood plasma. The particle phase is approximated as monodispersed solid spheres and the liquid phase as one-phase Newtonian fluid. The pulsatile liquid flow is computed by PISO algorithm based on the FEM(Finite Element Method) which has second-order time accuracy. The particle trajectories are obtained by solving the equation of motion in which the force term consists of Stokes’ drag force acting on a sphere in a creeping flow. Particle motions are simulated by the SPT(Simultaneous Particle Tracking) method. The leukocyte number density(LND) on the endothelium of the blood vessel is regarded important in the atherogenesis. The LND has contradicting aspects depending on what leukocyte-endothelium interaction model we take, namely, depending on whether the particle adheres to or reflected from the wall upon its impact. In this thesis, both the two leukocyte-endothelium interaction models are taken into consideration. Time-dependent shear stress and leukocyte distribution in the pulsatile blood flow field are presented. The time-averaged values of wall shear stress and LND are also evaluated.