An expanded computational model of surface induced thrombin generation was developed that includes hemodynamic effects, 22 biochemical reactions and 44 distinct chemical species. Surface binding of factors V, VIII, IX, and X was included in order to more accurately simulate the formation of the surface complexes tenase and prothrombinase. In order to model these reactions, the non-activated, activated and inactivated forms were all considered. This model was used to investigate the impact of surface bound heparin on thrombin generation with and without the additive effects of thrombomodulin (TM). In total, 104 heparin/TM pairings were evaluated (52 under venous conditions, 52 under arterial conditions), the results demonstrating the synergistic ability of heparin and TM to reduce thrombin generation. Additionally, the role of circulating tissue factor (TFp) was investigated and compared to that of surface-bound tissue factor (TFs). The numerical results suggest that circulating TF has the power to amplify thrombin generation once the coagulation cascade is already initiated by surface-bound TF. TFp concentrations as low as 0.01 nM were found to have a significant impact on total thrombin generation.