A computer-implemented method of reducing 4D Digital Subtracted Angiography (DSA) reconstruction artifacts using a computational fluid dynamics (CFD) simulation includes a computer receiving first DSA time sequence data comprising a representation of a plurality of vessels and segmenting a vessel of interest from the first DSA time sequence data. The computer uses the CFD simulation to simulate fluid dynamics across the vessel of interest to yield a flow field and determines a plurality of simulated time activity curve parameters for each voxel inside the vessel of interest using the flow field. Then, the computer applies a reconstruction process to second DSA time sequence data to yield a DSA volume. This reconstruction process is constrained by the plurality of simulated time activity curve parameters for each voxel inside the vessel of interest.