Smoothed particle hydrodynamics method applied to pulsatile flow inside a rigid two-dimensional model of left heart cavity

S. Shahriari, L. Kadem, B. D. Rogers, Ibrahim Hassan

Research output: Contribution to journalArticle

17 Citations (Scopus)

Abstract

This paper aims to extend the application of smoothed particle hydrodynamics (SPH), a meshfree particle method, to simulate flow inside a model of the heart's left ventricle (LV). This work is considered the first attempt to simulate flow inside a heart cavity using a meshfree particle method. Simulating this kind of flow, characterized by high pulsatility and moderate Reynolds number using SPH is challenging. As a consequence, validation of the computational code using benchmark cases is required prior to simulating the flow inside a model of the LV. In this work, this is accomplished by simulating an unsteady oscillating flow (pressure amplitude: A=2500N/m3 and Womersley number: Wo=16) and the steady lid-driven cavity flow (Re=3200,5000). The results are compared against analytical solutions and reference data to assess convergence. Then, both benchmark cases are combined and a pulsatile jet in a cavity is simulated and the results are compared with the finite volume method. Here, an approach to deal with inflow and outflow boundary conditions is introduced. Finally, pulsatile inlet flow in a rigid model of the LV is simulated. The results demonstrate the ability of SPH to model complex cardiovascular flows and to track the history of fluid properties. Some interesting features of SPH are also demonstrated in this study, including the relation between particle resolution and sound speed to control compressibility effects and also order of convergence in SPH simulations, which is consistently demonstrated to be between first-order and second-order at the moderate Reynolds numbers investigated.

Original languageEnglish
Pages (from-to)1121-1143
Number of pages23
JournalInternational Journal for Numerical Methods in Biomedical Engineering
Volume28
Issue number11
DOIs
Publication statusPublished - Nov 2012
Externally publishedYes

Fingerprint

Pulsatile Flow
Pulsatile flow
Hydrodynamics
Cavity
Heart Ventricles
Left Ventricle
Benchmarking
Meshfree Method
Particle Method
Reynolds number
Cardiovascular Models
Oscillating flow
Inlet flow
Model
Benchmark
Driven Cavity Flow
Finite volume method
Lid-driven Cavity
Compressibility
Order of Convergence

Keywords

  • Left heart ventricle
  • Moderate Reynolds number
  • Oscillating and pulsatile flows
  • Particle resolution
  • Smoothed particle hydrodynamics

ASJC Scopus subject areas

  • Software
  • Modelling and Simulation
  • Biomedical Engineering
  • Molecular Biology
  • Computational Theory and Mathematics
  • Applied Mathematics

Cite this

Smoothed particle hydrodynamics method applied to pulsatile flow inside a rigid two-dimensional model of left heart cavity. / Shahriari, S.; Kadem, L.; Rogers, B. D.; Hassan, Ibrahim.

In: International Journal for Numerical Methods in Biomedical Engineering, Vol. 28, No. 11, 11.2012, p. 1121-1143.

Research output: Contribution to journalArticle

@article{cef9cef6f1f24a44a142f79630894c10,
title = "Smoothed particle hydrodynamics method applied to pulsatile flow inside a rigid two-dimensional model of left heart cavity",
abstract = "This paper aims to extend the application of smoothed particle hydrodynamics (SPH), a meshfree particle method, to simulate flow inside a model of the heart's left ventricle (LV). This work is considered the first attempt to simulate flow inside a heart cavity using a meshfree particle method. Simulating this kind of flow, characterized by high pulsatility and moderate Reynolds number using SPH is challenging. As a consequence, validation of the computational code using benchmark cases is required prior to simulating the flow inside a model of the LV. In this work, this is accomplished by simulating an unsteady oscillating flow (pressure amplitude: A=2500N/m3 and Womersley number: Wo=16) and the steady lid-driven cavity flow (Re=3200,5000). The results are compared against analytical solutions and reference data to assess convergence. Then, both benchmark cases are combined and a pulsatile jet in a cavity is simulated and the results are compared with the finite volume method. Here, an approach to deal with inflow and outflow boundary conditions is introduced. Finally, pulsatile inlet flow in a rigid model of the LV is simulated. The results demonstrate the ability of SPH to model complex cardiovascular flows and to track the history of fluid properties. Some interesting features of SPH are also demonstrated in this study, including the relation between particle resolution and sound speed to control compressibility effects and also order of convergence in SPH simulations, which is consistently demonstrated to be between first-order and second-order at the moderate Reynolds numbers investigated.",
keywords = "Left heart ventricle, Moderate Reynolds number, Oscillating and pulsatile flows, Particle resolution, Smoothed particle hydrodynamics",
author = "S. Shahriari and L. Kadem and Rogers, {B. D.} and Ibrahim Hassan",
year = "2012",
month = "11",
doi = "10.1002/cnm.2482",
language = "English",
volume = "28",
pages = "1121--1143",
journal = "International Journal for Numerical Methods in Biomedical Engineering",
issn = "2040-7939",
publisher = "Wiley-Blackwell",
number = "11",

}

TY - JOUR

T1 - Smoothed particle hydrodynamics method applied to pulsatile flow inside a rigid two-dimensional model of left heart cavity

AU - Shahriari, S.

AU - Kadem, L.

AU - Rogers, B. D.

AU - Hassan, Ibrahim

PY - 2012/11

Y1 - 2012/11

N2 - This paper aims to extend the application of smoothed particle hydrodynamics (SPH), a meshfree particle method, to simulate flow inside a model of the heart's left ventricle (LV). This work is considered the first attempt to simulate flow inside a heart cavity using a meshfree particle method. Simulating this kind of flow, characterized by high pulsatility and moderate Reynolds number using SPH is challenging. As a consequence, validation of the computational code using benchmark cases is required prior to simulating the flow inside a model of the LV. In this work, this is accomplished by simulating an unsteady oscillating flow (pressure amplitude: A=2500N/m3 and Womersley number: Wo=16) and the steady lid-driven cavity flow (Re=3200,5000). The results are compared against analytical solutions and reference data to assess convergence. Then, both benchmark cases are combined and a pulsatile jet in a cavity is simulated and the results are compared with the finite volume method. Here, an approach to deal with inflow and outflow boundary conditions is introduced. Finally, pulsatile inlet flow in a rigid model of the LV is simulated. The results demonstrate the ability of SPH to model complex cardiovascular flows and to track the history of fluid properties. Some interesting features of SPH are also demonstrated in this study, including the relation between particle resolution and sound speed to control compressibility effects and also order of convergence in SPH simulations, which is consistently demonstrated to be between first-order and second-order at the moderate Reynolds numbers investigated.

AB - This paper aims to extend the application of smoothed particle hydrodynamics (SPH), a meshfree particle method, to simulate flow inside a model of the heart's left ventricle (LV). This work is considered the first attempt to simulate flow inside a heart cavity using a meshfree particle method. Simulating this kind of flow, characterized by high pulsatility and moderate Reynolds number using SPH is challenging. As a consequence, validation of the computational code using benchmark cases is required prior to simulating the flow inside a model of the LV. In this work, this is accomplished by simulating an unsteady oscillating flow (pressure amplitude: A=2500N/m3 and Womersley number: Wo=16) and the steady lid-driven cavity flow (Re=3200,5000). The results are compared against analytical solutions and reference data to assess convergence. Then, both benchmark cases are combined and a pulsatile jet in a cavity is simulated and the results are compared with the finite volume method. Here, an approach to deal with inflow and outflow boundary conditions is introduced. Finally, pulsatile inlet flow in a rigid model of the LV is simulated. The results demonstrate the ability of SPH to model complex cardiovascular flows and to track the history of fluid properties. Some interesting features of SPH are also demonstrated in this study, including the relation between particle resolution and sound speed to control compressibility effects and also order of convergence in SPH simulations, which is consistently demonstrated to be between first-order and second-order at the moderate Reynolds numbers investigated.

KW - Left heart ventricle

KW - Moderate Reynolds number

KW - Oscillating and pulsatile flows

KW - Particle resolution

KW - Smoothed particle hydrodynamics

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

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

U2 - 10.1002/cnm.2482

DO - 10.1002/cnm.2482

M3 - Article

VL - 28

SP - 1121

EP - 1143

JO - International Journal for Numerical Methods in Biomedical Engineering

JF - International Journal for Numerical Methods in Biomedical Engineering

SN - 2040-7939

IS - 11

ER -