Bloodstain pattern analysis is most important part of forensic science to investigate the homicide crime scene. Bloodstain pattern is determined by impact conditions of blood drop such as drop diameter, impact angle and velocity. Source of blood drop which represents the position of victim can be detected with these impact conditions of blood drop. So, it is significant to find out the correlation between impact conditions of blood drop and bloodstain pattern. In the present study, impacting motion of blood drop was numerically simulated to predict the bloodstain pattern. 3-D transient simulations of two phase interfacial flow were conducted by using the VOF method. Non-Newtonian viscosity and dynamic contact angle were used to get the accurate prediction of bloodstain pattern. Simulation results showed a good agreement with the experimental data of bloodstain on various impact angles. With the validated numerical methods, bloodstain patterns were predicted with various Reynolds numbers and Weber numbers. Characteristics of bloodstain were analyzed with bloodstain area, oval shape and fingering phenomena around the rim of bloodstain.