In this presentation, I will present a highly parallel method for the incompressible Navier-Stokes and Darcy equations for the simulation of the blood flows in the full three-dimensional patient-specific human liver, which include hepatic artery, portal vein, hepatic vein and hepatic tissue. To compute the blood flows, a scalable parallel method is used to implicitly solve the unsteady incompressible Navier-Stokes and Darcy equations discretised with a stabilized finite element method on fully unstructured meshes. The parallel algebraic solver includes a Newton method, a Krylov subspace method (GMRES) and an overlapping Schwarz preconditioner. As applications, I also simulate the flow in a patient with hepatectomy and calculate the Portal Pressure Gradient (PPG)， where PPG is a gold standard value to assess the portal hypertension. Moreover, the robustness and scalability of the algorithm are also investigated. At last, some extensions about nonlocal fluid model will also be introduced.