Numerical simulation of red blood cell distributions in three-dimensional microvascular bifurcations

Toru Hyakutake, Shinya Nagai

    • 9 Citations

    Abstract

    We constructed three-dimensional microvascular bifurcation models using a parent vessel of diameter 10. μm and investigated the flow behavior of the red blood cells (RBCs) through bifurcations. We considered symmetric and asymmetric model types. Two cases of equal daughter vessel diameter were employed for the asymmetric models, where the first was 10. μm, which is the same as the parent vessel and the second was 7.94. μm, which satisfies Murray's law. Simulated blood flow was computed using the lattice Boltzmann method in conjunction with the immersed boundary method for incorporating fluid-membrane interactions between the flow field and deformable RBCs. First, we investigated the flow behavior of a single RBC through microvascular bifurcations. In the case of the symmetric bifurcation, the turning point of the fractional plasma flow wherein the RBC flow changed from one daughter vessel to the other was 0.50. This turning point was however different for asymmetric bifurcations. Additionally, we varied the initial offset of RBCs from the centerline of the parent vessel. The simulation results indicated that the RBCs preferentially flow through the branch of a larger flow ratio. Next, we investigated the distribution characteristics of multiple RBCs. Simulations indicated that the results of the symmetric model were similar to those predicted by a previously published empirical model. On the other hand, results of asymmetric models deviated from those of the symmetric and empirical models. These results suggest that the distribution of RBCs varies according to the bifurcation angle and daughter vessel diameter in a microvascular bifurcation of the size considered.

    Original languageEnglish
    Pages (from-to)115-123
    Number of pages9
    JournalMicrovascular Research
    Volume97
    DOIs
    StatePublished - 2015 Jan 1

    Fingerprint

    Erythrocytes
    Membranes

    Keywords

    • Immersed boundary method
    • Lattice Boltzmann method
    • Microvascular bifurcation
    • Numerical simulation
    • Red blood cell

    ASJC Scopus subject areas

    • Biochemistry
    • Cardiology and Cardiovascular Medicine
    • Cell Biology

    Cite this

    Numerical simulation of red blood cell distributions in three-dimensional microvascular bifurcations. / Hyakutake, Toru; Nagai, Shinya.

    In: Microvascular Research, Vol. 97, 01.01.2015, p. 115-123.

    Research output: Contribution to journalArticle

    Hyakutake, Toru; Nagai, Shinya / Numerical simulation of red blood cell distributions in three-dimensional microvascular bifurcations.

    In: Microvascular Research, Vol. 97, 01.01.2015, p. 115-123.

    Research output: Contribution to journalArticle

    @article{a6f54a87905e4473bee9c6f97355c727,
    title = "Numerical simulation of red blood cell distributions in three-dimensional microvascular bifurcations",
    keywords = "Immersed boundary method, Lattice Boltzmann method, Microvascular bifurcation, Numerical simulation, Red blood cell",
    author = "Toru Hyakutake and Shinya Nagai",
    year = "2015",
    month = "1",
    doi = "10.1016/j.mvr.2014.10.001",
    volume = "97",
    pages = "115--123",
    journal = "Microvascular Research",
    issn = "0026-2862",
    publisher = "Academic Press Inc.",

    }

    TY - JOUR

    T1 - Numerical simulation of red blood cell distributions in three-dimensional microvascular bifurcations

    AU - Hyakutake,Toru

    AU - Nagai,Shinya

    PY - 2015/1/1

    Y1 - 2015/1/1

    N2 - We constructed three-dimensional microvascular bifurcation models using a parent vessel of diameter 10. μm and investigated the flow behavior of the red blood cells (RBCs) through bifurcations. We considered symmetric and asymmetric model types. Two cases of equal daughter vessel diameter were employed for the asymmetric models, where the first was 10. μm, which is the same as the parent vessel and the second was 7.94. μm, which satisfies Murray's law. Simulated blood flow was computed using the lattice Boltzmann method in conjunction with the immersed boundary method for incorporating fluid-membrane interactions between the flow field and deformable RBCs. First, we investigated the flow behavior of a single RBC through microvascular bifurcations. In the case of the symmetric bifurcation, the turning point of the fractional plasma flow wherein the RBC flow changed from one daughter vessel to the other was 0.50. This turning point was however different for asymmetric bifurcations. Additionally, we varied the initial offset of RBCs from the centerline of the parent vessel. The simulation results indicated that the RBCs preferentially flow through the branch of a larger flow ratio. Next, we investigated the distribution characteristics of multiple RBCs. Simulations indicated that the results of the symmetric model were similar to those predicted by a previously published empirical model. On the other hand, results of asymmetric models deviated from those of the symmetric and empirical models. These results suggest that the distribution of RBCs varies according to the bifurcation angle and daughter vessel diameter in a microvascular bifurcation of the size considered.

    AB - We constructed three-dimensional microvascular bifurcation models using a parent vessel of diameter 10. μm and investigated the flow behavior of the red blood cells (RBCs) through bifurcations. We considered symmetric and asymmetric model types. Two cases of equal daughter vessel diameter were employed for the asymmetric models, where the first was 10. μm, which is the same as the parent vessel and the second was 7.94. μm, which satisfies Murray's law. Simulated blood flow was computed using the lattice Boltzmann method in conjunction with the immersed boundary method for incorporating fluid-membrane interactions between the flow field and deformable RBCs. First, we investigated the flow behavior of a single RBC through microvascular bifurcations. In the case of the symmetric bifurcation, the turning point of the fractional plasma flow wherein the RBC flow changed from one daughter vessel to the other was 0.50. This turning point was however different for asymmetric bifurcations. Additionally, we varied the initial offset of RBCs from the centerline of the parent vessel. The simulation results indicated that the RBCs preferentially flow through the branch of a larger flow ratio. Next, we investigated the distribution characteristics of multiple RBCs. Simulations indicated that the results of the symmetric model were similar to those predicted by a previously published empirical model. On the other hand, results of asymmetric models deviated from those of the symmetric and empirical models. These results suggest that the distribution of RBCs varies according to the bifurcation angle and daughter vessel diameter in a microvascular bifurcation of the size considered.

    KW - Immersed boundary method

    KW - Lattice Boltzmann method

    KW - Microvascular bifurcation

    KW - Numerical simulation

    KW - Red blood cell

    UR - http://www.scopus.com/inward/record.url?scp=84918513220&partnerID=8YFLogxK

    UR - http://www.scopus.com/inward/citedby.url?scp=84918513220&partnerID=8YFLogxK

    U2 - 10.1016/j.mvr.2014.10.001

    DO - 10.1016/j.mvr.2014.10.001

    M3 - Article

    VL - 97

    SP - 115

    EP - 123

    JO - Microvascular Research

    T2 - Microvascular Research

    JF - Microvascular Research

    SN - 0026-2862

    ER -