Modeling of three-dimensional blood flow in microchannels using a two-fluid method

Abstract

This work presents a novel two-fluid method based on our recently proposed viscosity model for red blood cells (RBCs)—for simulating three-dimensional (3D) blood flow in a microchannel of dimension comparable to the diameter of red blood cells and larger. Toward this, whole blood is assumed as a suspension of red blood cells in blood plasma, with each phase considered as interpenetrating continua having its separate mass and momentum conservation equations. The proposed approach-based performance study is presented after comprehensively validating it with experimental data for blood flow in a uniform, sudden expansion-constriction, and Y-shaped bifurcated rectangular microchannels over—an extensive range of size (25–330 μm), flow rates (11.8 μl/h–30 ml/h), and inlet hematocrit (0%–45%). The proposed approach effectively captures significant biophysical and biomechanical insights into blood flow. It highlights a migration of red blood cells toward the center of the microchannel and the formation of a cell-free layer near the wall. Notably, with the introduction of constriction and expansion in the microchannel, it predicts a fivefold enhancement of the cell-free layer. The Fahraeus and Fahraeus–Lindquist effects are also demonstrated in microchannels, with less than 300 μm characteristic dimensions. These findings are consistent with experimental evidence. In addition to experimentally evident phenomena, our simulations unveil several additional flow phenomena and features of blood flow in the microchannel. It is observed that the presence of confluence (merging flow) is more disturbing to the blood flow than the presence of diverging bifurcations (splitting flow). Furthermore, after the confluence, velocity profiles exhibit a local peak that persists up to the microchannel outlet. Primary contribution of this work lies in the proposal of a two-fluid method for simulating 3D blood flow in complex geometries. This approach provides a comprehensive understanding of blood flow dynamics in microchannels and can be applied to optimize dimensions and geometries during the initial phases of plasma separation microdevices development.