AbstractLower limb peripheral artery disease is a prevalent chronic non-communicable disease without obvious symptoms. However, the effect of ischemic lower limb peripheral arteries on hemodynamics remains unclear. In this study, we investigated the variation of the hemodynamics caused by patient-specific structural artery characteristics. Computational fluid dynamic simulations were performed on seven lower limb (including superficial femoral, deep femoral and popliteal) artery models that were reconstructed from magnetic resonance imaging. We found that increased wall shear stress (WSS) was mainly caused by the increasing severity of stenosis, bending, and branching. Our results showed that the increase in the WSS value at a stenosis at the bifurcation was 2.7 Pa. In contrast, the isolated stenosis and branch caused a WSS increase of 0.7 Pa and 0.5 Pa, respectively. The WSS in the narrow popliteal artery was more sensitive to a reduction in radius. Our results also demonstrate that the distribution of the velocity and pressure gradient are highly structurally related. At last, Ultrasound Doppler velocimeter measured result was presented as a validation. In conclusion, the distribution of hemodynamics may serve as a supplement for clinical decision-making to prevent the occurrence of a morbid or mortal ischemic event.