Thromboembolic complications in Bileaflet mechanical heart valves (BMHVs) are believed to be due to the combination of high shear stresses and large recirculation regions. Relating blood damage to design geometry is therefore essential to ultimately optimize the design of BMHVs. The aim of this research is to quantitatively study the effect of 3D channel geometry on shear-induced platelet activation and aggregation, and to choose an appropriate blood damage index (BDI) model for future numerical simulations. The simulations in this study use a recently developed lattice-Boltzmann with external boundary force (LBM-EBF) method [Wu, J., and C. K. Aidun. Int. J. Numer. Method Fluids 62(7):765-783, 2010; Wu, J., and C. K. Aidun. Int. J. Multiphase flow 36:202-209, 2010]. The channel geometries and flow conditions are re-constructed from recent experiments by Fallon [The Development of a Novel in vitro Flow System to Evaluate Platelet Activation and Procoagulant Potential Induced by Bileaflet Mechanical Heart Valve Leakage Jets in School of Chemical and Biomolecular Engineering. Atlanta: Georgia Institute of Technology] and Fallon et al. [Ann. Biomed. Eng. 36(1):1]. The fluid flow is computed on a fixed regular 'lattice' using the LBM, and each platelet is mapped onto a Lagrangian frame moving continuously throughout the fluid domain. The two-way fluid-solid interactions are determined by the EBF method by enforcing a no-slip condition on the platelet surface. The motion and orientation of the platelet are obtained from Newtonian dynamics equations. The numerical results show that sharp corners or sudden shape transitions will increase blood damage. Fallon's experimental results were used as a basis for choosing the appropriate BDI model for use in future computational simulations of flow through BMHVs.