Search

The goal of this study is to develop a computational fluid dynamics (CFD) modeling approach to better estimate the blood flow dynamics in the bundles of the hollow fiber membrane based medical devices (i.e., blood oxygenators, artificial lungs, and hemodialyzers). Three representative types of arrays, square, diagonal, and random with the porosity value of 0.55, were studied. In addition, a 3D array with the same porosity was studied. The flow fields between the individual fibers in these arrays at selected Reynolds numbers (Re) were simulated with CFD modeling. Hemolysis is not significant in the fiber bundles but the platelet activation may be essential. For each type of array, the average wall shear stress is linearly proportional to the Re. For the same Re but different arrays, the average wall shear stress also exhibits a linear dependency on the pressure difference across arrays, while Darcy's law prescribes a power-law relationship, therefore, underestimating the shear stress level. For the same Re, the average wall shear stress of the diagonal array is approximately 3.1, 1.8, and 2.0 times larger than that of the square, random, and 3D arrays, respectively. A coefficient C is suggested to correlate the CFD predicted data with the analytical solution, and C is 1.16, 1.51, and 2.05 for the square, random, and diagonal arrays in this paper, respectively. It is worth noting that C is strongly dependent on the array geometrical properties, whereas it is weakly dependent on the flow field. Additionally, the 3D fiber bundle simulation results show that the three-dimensional effect is not negligible. Specifically, velocity and shear stress distribution can vary significantly along the fiber axial direction.

Ventricular assist devices (VADs) have already helped many patients with heart failure but have the potential to assist more patients if current problems with blood damage (hemolysis, platelet activation, thrombosis and emboli, and destruction of the von Willebrand factor (vWf)) can be eliminated. A step towards this goal is better understanding of the relationships between shear stress, exposure time, and blood damage and, from there, the development of numerical models for the different types of blood damage to enable the design of improved VADs. In this study, computational fluid dynamics (CFD) was used to calculate the hemodynamics in three clinical VADs and two investigational VADs and the shear stress, residence time, and hemolysis were investigated. A new scalar transport model for hemolysis was developed. The results were compared with in vitro measurements of the pressure head in each VAD and the hemolysis index in two VADs. A comparative analysis of the blood damage related fluid dynamic parameters and hemolysis index was performed among the VADs. Compared to the centrifugal VADs, the axial VADs had: higher mean scalar shear stress (sss); a wider range of sss, with larger maxima and larger percentage volumes at both low and high sss; and longer residence times at very high sss. The hemolysis predictions were in agreement with the experiments and showed that the axial VADs had a higher hemolysis index. The increased hemolysis in axial VADs compared to centrifugal VADs is a direct result of their higher shear stresses and longer residence times. Since platelet activation and destruction of the vWf also require high shear stresses, the flow conditions inside axial VADs are likely to result in more of these types of blood damage compared with centrifugal VADs.