Abstract: Image data is able to provide increasingly detailed information of anatomy and flow, but imaging alone is not a predictive tool. Patient-specific blood flow simulations have the potential to provide quantitative predictive tools for virtual surgery, treatment planning, and risk stratification. To accurately resolve the blood flows based on the patient-specific geometry and parameters is still a big challenge because of the complex geometry and the turbulence, and it is also important to obtain the results in a short amount of computing time so that the simulation can be used in surgery planning. In this talk, we will discuss some scalable parallel methods for the simulation of blood flow in compliant arteries on large scale supercomputers. The blood flow is modeled by 3D unsteady incompressible Navier-Stokes equations with a lumped parameter boundary condition, which are discretized with a stabilized finite element based on unstructured meshes in space and a fully implicit method in time. The large scale discretized nonlinear systems are solved by a parallel Newton-Krylov method preconditioned by linear and nonlinear domain decomposition preconditioners.