This work presents a computational study of coronary artery disease using finite element simulations of incompressible blood flow. The Navier–Stokes equations are solved with a custom scheme implemented in FEniCSx to analyze key hemodynamic variables such as velocity, pressure, wall shear stress and vorticity. Blood flow is simulated through healthy and stenosed coronary artery geometries to assess the influence of morphological changes on local hemodynamics.