Heart valve flow computation with the integrated Space-Time VMS, Slip Interface, Topology Change and Isogeometric Discretization methods

Kenji Takizawa, Tayfun E. Tezduyar, Takuya Terahara, Takafumi Sasaki

    Research output: Contribution to journalArticle

    33 Citations (Scopus)

    Abstract

    Heart valve flow computation requires accurate representation of boundary layers near moving solid surfaces, including the valve leaflet surfaces, even when the leaflets come into contact. It also requires dealing with a high level of geometric complexity. We address these computational challenges with a Space-Time (ST) method developed by integrating three special ST methods in the framework of the ST Variational Multiscale (ST-VMS) method. The special methods are the ST Slip Interface (ST-SI) and ST Topology Change (ST-TC) methods and ST Isogeometric Analysis (ST-IGA). The computations are for a realistic aortic-valve model with prescribed valve leaflet motion and actual contact between the leaflets. The ST-VMS method functions as a moving-mesh method, which maintains high-resolution boundary layer representation near the solid surfaces, including leaflet surfaces. The ST-TC method was introduced for moving-mesh computation of flow problems with TC, such as contact between the leaflets of a heart valve. It deals with the contact while maintaining high-resolution representation near the leaflet surfaces. The ST-SI method was originally introduced to have high-resolution representation of the boundary layers near spinning solid surfaces. The mesh covering a spinning solid surface spins with it, and the SI between the spinning mesh and the rest of the mesh accurately connects the two sides. In the context of heart valves, the SI connects the sectors of meshes containing the leaflets, enabling a more effective mesh moving. In that context, integration of the ST-SI and ST-TC methods enables high-resolution representation even when the contact is between leaflets that are covered by meshes with SI. It also enables dealing with contact location change or contact and sliding on the SI. By integrating the ST-IGA with the ST-SI and ST-TC methods, in addition to having a more accurate representation of the surfaces and increased accuracy in the flow solution, the element density in the narrow spaces near the contact areas is kept at a reasonable level. Furthermore, because the flow representation in the contact area has a wider support in IGA, the flow computation method becomes more robust. The computations we present for an aortic-valve model with two different modes of prescribed leaflet motion show the effectiveness of the ST-SI-TC-IGA method.

    Original languageEnglish
    JournalComputers and Fluids
    DOIs
    Publication statusAccepted/In press - 2016 Aug 19

    Fingerprint

    Topology
    Boundary layers

    Keywords

    • Contact
    • Heart valve flow
    • Isogeometric Analysis
    • Mesh generation
    • Space-Time Slip Interface method
    • Space-Time Topology Change method
    • Space-Time VMS method

    ASJC Scopus subject areas

    • Computer Science(all)
    • Engineering(all)

    Cite this

    @article{8f8a03569ccb4b6b9c2276ab737147ee,
    title = "Heart valve flow computation with the integrated Space-Time VMS, Slip Interface, Topology Change and Isogeometric Discretization methods",
    abstract = "Heart valve flow computation requires accurate representation of boundary layers near moving solid surfaces, including the valve leaflet surfaces, even when the leaflets come into contact. It also requires dealing with a high level of geometric complexity. We address these computational challenges with a Space-Time (ST) method developed by integrating three special ST methods in the framework of the ST Variational Multiscale (ST-VMS) method. The special methods are the ST Slip Interface (ST-SI) and ST Topology Change (ST-TC) methods and ST Isogeometric Analysis (ST-IGA). The computations are for a realistic aortic-valve model with prescribed valve leaflet motion and actual contact between the leaflets. The ST-VMS method functions as a moving-mesh method, which maintains high-resolution boundary layer representation near the solid surfaces, including leaflet surfaces. The ST-TC method was introduced for moving-mesh computation of flow problems with TC, such as contact between the leaflets of a heart valve. It deals with the contact while maintaining high-resolution representation near the leaflet surfaces. The ST-SI method was originally introduced to have high-resolution representation of the boundary layers near spinning solid surfaces. The mesh covering a spinning solid surface spins with it, and the SI between the spinning mesh and the rest of the mesh accurately connects the two sides. In the context of heart valves, the SI connects the sectors of meshes containing the leaflets, enabling a more effective mesh moving. In that context, integration of the ST-SI and ST-TC methods enables high-resolution representation even when the contact is between leaflets that are covered by meshes with SI. It also enables dealing with contact location change or contact and sliding on the SI. By integrating the ST-IGA with the ST-SI and ST-TC methods, in addition to having a more accurate representation of the surfaces and increased accuracy in the flow solution, the element density in the narrow spaces near the contact areas is kept at a reasonable level. Furthermore, because the flow representation in the contact area has a wider support in IGA, the flow computation method becomes more robust. The computations we present for an aortic-valve model with two different modes of prescribed leaflet motion show the effectiveness of the ST-SI-TC-IGA method.",
    keywords = "Contact, Heart valve flow, Isogeometric Analysis, Mesh generation, Space-Time Slip Interface method, Space-Time Topology Change method, Space-Time VMS method",
    author = "Kenji Takizawa and Tezduyar, {Tayfun E.} and Takuya Terahara and Takafumi Sasaki",
    year = "2016",
    month = "8",
    day = "19",
    doi = "10.1016/j.compfluid.2016.11.012",
    language = "English",
    journal = "Computers and Fluids",
    issn = "0045-7930",
    publisher = "Elsevier Limited",

    }

    TY - JOUR

    T1 - Heart valve flow computation with the integrated Space-Time VMS, Slip Interface, Topology Change and Isogeometric Discretization methods

    AU - Takizawa, Kenji

    AU - Tezduyar, Tayfun E.

    AU - Terahara, Takuya

    AU - Sasaki, Takafumi

    PY - 2016/8/19

    Y1 - 2016/8/19

    N2 - Heart valve flow computation requires accurate representation of boundary layers near moving solid surfaces, including the valve leaflet surfaces, even when the leaflets come into contact. It also requires dealing with a high level of geometric complexity. We address these computational challenges with a Space-Time (ST) method developed by integrating three special ST methods in the framework of the ST Variational Multiscale (ST-VMS) method. The special methods are the ST Slip Interface (ST-SI) and ST Topology Change (ST-TC) methods and ST Isogeometric Analysis (ST-IGA). The computations are for a realistic aortic-valve model with prescribed valve leaflet motion and actual contact between the leaflets. The ST-VMS method functions as a moving-mesh method, which maintains high-resolution boundary layer representation near the solid surfaces, including leaflet surfaces. The ST-TC method was introduced for moving-mesh computation of flow problems with TC, such as contact between the leaflets of a heart valve. It deals with the contact while maintaining high-resolution representation near the leaflet surfaces. The ST-SI method was originally introduced to have high-resolution representation of the boundary layers near spinning solid surfaces. The mesh covering a spinning solid surface spins with it, and the SI between the spinning mesh and the rest of the mesh accurately connects the two sides. In the context of heart valves, the SI connects the sectors of meshes containing the leaflets, enabling a more effective mesh moving. In that context, integration of the ST-SI and ST-TC methods enables high-resolution representation even when the contact is between leaflets that are covered by meshes with SI. It also enables dealing with contact location change or contact and sliding on the SI. By integrating the ST-IGA with the ST-SI and ST-TC methods, in addition to having a more accurate representation of the surfaces and increased accuracy in the flow solution, the element density in the narrow spaces near the contact areas is kept at a reasonable level. Furthermore, because the flow representation in the contact area has a wider support in IGA, the flow computation method becomes more robust. The computations we present for an aortic-valve model with two different modes of prescribed leaflet motion show the effectiveness of the ST-SI-TC-IGA method.

    AB - Heart valve flow computation requires accurate representation of boundary layers near moving solid surfaces, including the valve leaflet surfaces, even when the leaflets come into contact. It also requires dealing with a high level of geometric complexity. We address these computational challenges with a Space-Time (ST) method developed by integrating three special ST methods in the framework of the ST Variational Multiscale (ST-VMS) method. The special methods are the ST Slip Interface (ST-SI) and ST Topology Change (ST-TC) methods and ST Isogeometric Analysis (ST-IGA). The computations are for a realistic aortic-valve model with prescribed valve leaflet motion and actual contact between the leaflets. The ST-VMS method functions as a moving-mesh method, which maintains high-resolution boundary layer representation near the solid surfaces, including leaflet surfaces. The ST-TC method was introduced for moving-mesh computation of flow problems with TC, such as contact between the leaflets of a heart valve. It deals with the contact while maintaining high-resolution representation near the leaflet surfaces. The ST-SI method was originally introduced to have high-resolution representation of the boundary layers near spinning solid surfaces. The mesh covering a spinning solid surface spins with it, and the SI between the spinning mesh and the rest of the mesh accurately connects the two sides. In the context of heart valves, the SI connects the sectors of meshes containing the leaflets, enabling a more effective mesh moving. In that context, integration of the ST-SI and ST-TC methods enables high-resolution representation even when the contact is between leaflets that are covered by meshes with SI. It also enables dealing with contact location change or contact and sliding on the SI. By integrating the ST-IGA with the ST-SI and ST-TC methods, in addition to having a more accurate representation of the surfaces and increased accuracy in the flow solution, the element density in the narrow spaces near the contact areas is kept at a reasonable level. Furthermore, because the flow representation in the contact area has a wider support in IGA, the flow computation method becomes more robust. The computations we present for an aortic-valve model with two different modes of prescribed leaflet motion show the effectiveness of the ST-SI-TC-IGA method.

    KW - Contact

    KW - Heart valve flow

    KW - Isogeometric Analysis

    KW - Mesh generation

    KW - Space-Time Slip Interface method

    KW - Space-Time Topology Change method

    KW - Space-Time VMS method

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

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

    U2 - 10.1016/j.compfluid.2016.11.012

    DO - 10.1016/j.compfluid.2016.11.012

    M3 - Article

    JO - Computers and Fluids

    JF - Computers and Fluids

    SN - 0045-7930

    ER -