In this paper novel finite element algorithms for robustly solving flow equations of asystemic circulation are presented and compared. A novel, Locally Conservative Galerkin (LCG)method is developed by adopting semi- and fully- implicit time discretisations to address the behaviourof blood flow in the human circulatory system. These techniques are efficient for nonlinearsystem of equations used in blood flow models and achieve rapid convergence. Severalcomparisons are made between the methods to demonstrate the validity and stability of the currentnumerical models to solve the blood flow characteristics in a human body
