Thromboembolic complications (TECs) of bileaflet mechanical heart valves (BMHVs) are believed to be due to the nonphysiologic mechanical stresses imposed on blood elements by the hinge flows. Relating hinge flow features to design features is, therefore, essential to ultimately design BMHVs with lower TEC rates. This study aims at simulating the pulsatile three-dimensional hinge flows of three BMHVs and estimating the TEC potential associated with each hinge design. Hinge geometries are constructed from micro-computed tomography scans of BMHVs. Simulations are conducted using a Cartesian sharp-interface immersed-boundary methodology combined with a second-order accurate fractional-step method. Leaflet motion and flow boundary conditions are extracted from fluid–structure-interaction simulations of BMHV bulk flow. The numerical results are analyzed using a particle-tracking approach coupled with existing blood damage models. The gap width and, more importantly, the shape of the recess and leaflet are found to impact the flow distribution and TEC potential. Smooth, streamlined surfaces appear to be more favorable than sharp corners or sudden shape transitions. The developed framework will enable pragmatic and cost-efficient preclinical evaluation of BMHV prototypes prior to valve manufacturing. Application to a wide range of hinges with varying design parameters will eventually help in determining the optimal hinge design.

Thromboembolic complications of bileaflet mechanical heart valves (BMHV) are believed to be due to detrimental stresses imposed on blood elements by the hinge flows. Characterization of these flows is thus crucial to identify the underlying causes for complications. In this study, we conduct three-dimensional pulsatile flow simulations through the hinge of a BMHV under aortic conditions. Hinge and leaflet geometries are reconstructed from the Micro-Computed Tomography scans of a BMHV. Simulations are conducted using a Cartesian sharp-interface immersed-boundary methodology combined with a second-order accurate fractional-step method. Physiologic flow boundary conditions and leaflet motion are extracted from the Fluid–Structure Interaction simulations of the bulk of the flow through a BMHV. Calculations reveal the presence, throughout the cardiac cycle, of flow patterns known to be detrimental to blood elements. Flow fields are characterized by: (1) complex systolic flows, with rotating structures and slow reverse flow pattern, and (2) two strong diastolic leakage jets accompanied by fast reverse flow at the hinge bottom. Elevated shear stresses, up to 1920 dyn/cm2 during systole and 6115 dyn/cm2 during diastole, are reported. This study underscores the need to conduct three-dimensional simulations throughout the cardiac cycle to fully characterize the complexity and thromboembolic potential of the hinge flows.