A Higher-order accurate Adaptive Mesh Refinement method is developed for blood flow simulation. The method based on the Interpolated Differential Operator scheme has third-order accuracy for Poisson equation. We integrate. Navier-Stokes equation by second-stage implicit Runge-Kutta method for numerical stability. We demonstrate the blood flow simulation with many branches, and show the effectiveness of the method.