Objective: The aim of this study is to numerically evaluate thrombosis risk within occluded coronary arterial fistulas (CAF) with terminal aneurysms, and provide guidance in choosing occlusion positions, with clinical observations as reference.Method: Four patients with CAF were studied, with different occlusion positions in actual treatments. Hemodynamics simulations were conducted, with blood residue predicted using the blood stasis model. Three types of models (untreated model, aneurysm-reserved model and aneurysm-removed model) were studeid for each patient. Four metrics, i.e., proportion of high oscillatory shear index (OSI), area of high OSI, old blood volume fraction (OBVF)) and old blood volume (OBV) was obtained to distinguish the thrombosis risk of different treatments (proximal or distal occlusion), comparing with the follow-up CTA.Results: For all the postopertive models, the high OBVF, high OSI(>0.3) and low time-averaged wall shear stress (TAWSS) regions were mainly at the distal fistula, indicating these regions were prone to thrombosis. The regions where blood residue remains are roughly regions of high OSI, corresponding well with clinical observations. In contrast, TAWSS failed to distinguish the difference in thrombosis risk. Absolute values (area of high OSI, OBV) can better reflect the degree of thrombosis risk between treatment types compared with percentage values (proportion of high OSI, OBVF). By comparing with the actual clinical treatments and observations, the OBV is superior to the area of high OSI in determining treatment type.Conclusion: The OBV, a volumetric parameter for blood stasis, can better account for the CAF thrombosis and reflect the degree of blood stasis compared with OSI or TAWSS, is a more appropriate metric for thrombosis in the fistula. Together with morphological parameters, the OBV could guide clinicians to formulate more appropriate surgical plans, which is of great significance for the preoperative evaluation and treatment prognosis of CAF patients.