Development of non-boundary-fitted Cartesian grid method for numerical simulation of mechanical heart valve and the potential for blood clotting

Computational fluid dynamics (CFD) simulations are becoming a reliable tool in understanding disease progression, investigating blood flow patterns and evaluating medical device performance such as mechanical heart valves (MHV). Previous studies indicated that the non-physiological flow patter...

Full description

Saved in:
Bibliographic Details
Main Author: Zakaria, Mohamad Shukri
Format: Thesis
Language:English
Published: 2018
Subjects:
Online Access:http://psasir.upm.edu.my/id/eprint/71429/1/FK%202018%2095%20-%20IR.pdf
Tags: Add Tag
No Tags, Be the first to tag this record!
id my-upm-ir.71429
record_format uketd_dc
spelling my-upm-ir.714292019-11-13T04:35:22Z Development of non-boundary-fitted Cartesian grid method for numerical simulation of mechanical heart valve and the potential for blood clotting 2018-04 Zakaria, Mohamad Shukri Computational fluid dynamics (CFD) simulations are becoming a reliable tool in understanding disease progression, investigating blood flow patterns and evaluating medical device performance such as mechanical heart valves (MHV). Previous studies indicated that the non-physiological flow pattern (i.e. recirculation, stagnation, and vortex) might cause a trapped platelet and be responsible for the formation of blood clots in MHV. Accurate simulation of this flow requires a high order accuracy numerical scheme together with a scale resolving turbulence model such as large eddy simulation (LES). This requires the use of uniform orthogonal grids for the descretisation process, which is not able to handle complex branching arterial domains that contain MHV, where the generation are usually boundary-fitted (BF) grid with non-orthogonality and distortions. Therefore, nonboundary fitted (NBF) Cartesian grid method is an alternative solution. The objective of this study is to develop a new NBF method based on the volume of fluid (VOF), containing the colour function, namely NBF-VOF Cartesian grid method. A single set of governing equation is used for both solid and fluids identified by unity colour function and zero colour function respectively. The solid was treated as a fluid with very high viscosity to theoretically reduce its deformability, and subsequently satisfy a no-slip condition at the boundary. In the first attempt, we found that in prior, the treatment was not satisfied. To suppress the fluid velocities in the solid, we introduced the artificial term derived from the colour function into an algebraic system of momentum equations, which had a significant impact on the originality of this study. The developed solver, NBF-VOF, is then thoroughly validated using a variety of numerical and experimental results available in the literature which is Hagen-Poiseuille flow, lid-driven cavity, flow over a cylinder, 90o tube flow, and pulsatile flow through the real anatomic aorta. Opensource CFD software was used as our simulation platform. Although the second order method degenerates the spatial accuracy of convergence rate as function of the grid size from 2 to 1.5, an agreement was found for all cases qualitative and quantitatively. The grid uncertainty obtained was less than 5%, which was within the acceptable range. The computational time was lower when the viscosity of solid was higher. However, higher solid viscosity gives lagging in the result for transient cases. Despite this, using higher time step, until the maximum Courant number of 4.0, can speed up the simulation time and preserved the stability. Finally, another breakthrough in this study was the application of the solver to simulate pulsatile blood flow of MHV placed in an axisymmetric and real patient anatomic aorta with the sinus, which reveals complex blood flow patterns, shear stress loading, and history of particles age in the local domain, that consequently can identified the potential of blood clotting. Heart valve prosthesis Aortic valve - Transplantation Tissue engineering 2018-04 Thesis http://psasir.upm.edu.my/id/eprint/71429/ http://psasir.upm.edu.my/id/eprint/71429/1/FK%202018%2095%20-%20IR.pdf text en public doctoral Universiti Putra Malaysia Heart valve prosthesis Aortic valve - Transplantation Tissue engineering
institution Universiti Putra Malaysia
collection PSAS Institutional Repository
language English
topic Heart valve prosthesis
Aortic valve - Transplantation
Tissue engineering
spellingShingle Heart valve prosthesis
Aortic valve - Transplantation
Tissue engineering
Zakaria, Mohamad Shukri
Development of non-boundary-fitted Cartesian grid method for numerical simulation of mechanical heart valve and the potential for blood clotting
description Computational fluid dynamics (CFD) simulations are becoming a reliable tool in understanding disease progression, investigating blood flow patterns and evaluating medical device performance such as mechanical heart valves (MHV). Previous studies indicated that the non-physiological flow pattern (i.e. recirculation, stagnation, and vortex) might cause a trapped platelet and be responsible for the formation of blood clots in MHV. Accurate simulation of this flow requires a high order accuracy numerical scheme together with a scale resolving turbulence model such as large eddy simulation (LES). This requires the use of uniform orthogonal grids for the descretisation process, which is not able to handle complex branching arterial domains that contain MHV, where the generation are usually boundary-fitted (BF) grid with non-orthogonality and distortions. Therefore, nonboundary fitted (NBF) Cartesian grid method is an alternative solution. The objective of this study is to develop a new NBF method based on the volume of fluid (VOF), containing the colour function, namely NBF-VOF Cartesian grid method. A single set of governing equation is used for both solid and fluids identified by unity colour function and zero colour function respectively. The solid was treated as a fluid with very high viscosity to theoretically reduce its deformability, and subsequently satisfy a no-slip condition at the boundary. In the first attempt, we found that in prior, the treatment was not satisfied. To suppress the fluid velocities in the solid, we introduced the artificial term derived from the colour function into an algebraic system of momentum equations, which had a significant impact on the originality of this study. The developed solver, NBF-VOF, is then thoroughly validated using a variety of numerical and experimental results available in the literature which is Hagen-Poiseuille flow, lid-driven cavity, flow over a cylinder, 90o tube flow, and pulsatile flow through the real anatomic aorta. Opensource CFD software was used as our simulation platform. Although the second order method degenerates the spatial accuracy of convergence rate as function of the grid size from 2 to 1.5, an agreement was found for all cases qualitative and quantitatively. The grid uncertainty obtained was less than 5%, which was within the acceptable range. The computational time was lower when the viscosity of solid was higher. However, higher solid viscosity gives lagging in the result for transient cases. Despite this, using higher time step, until the maximum Courant number of 4.0, can speed up the simulation time and preserved the stability. Finally, another breakthrough in this study was the application of the solver to simulate pulsatile blood flow of MHV placed in an axisymmetric and real patient anatomic aorta with the sinus, which reveals complex blood flow patterns, shear stress loading, and history of particles age in the local domain, that consequently can identified the potential of blood clotting.
format Thesis
qualification_level Doctorate
author Zakaria, Mohamad Shukri
author_facet Zakaria, Mohamad Shukri
author_sort Zakaria, Mohamad Shukri
title Development of non-boundary-fitted Cartesian grid method for numerical simulation of mechanical heart valve and the potential for blood clotting
title_short Development of non-boundary-fitted Cartesian grid method for numerical simulation of mechanical heart valve and the potential for blood clotting
title_full Development of non-boundary-fitted Cartesian grid method for numerical simulation of mechanical heart valve and the potential for blood clotting
title_fullStr Development of non-boundary-fitted Cartesian grid method for numerical simulation of mechanical heart valve and the potential for blood clotting
title_full_unstemmed Development of non-boundary-fitted Cartesian grid method for numerical simulation of mechanical heart valve and the potential for blood clotting
title_sort development of non-boundary-fitted cartesian grid method for numerical simulation of mechanical heart valve and the potential for blood clotting
granting_institution Universiti Putra Malaysia
publishDate 2018
url http://psasir.upm.edu.my/id/eprint/71429/1/FK%202018%2095%20-%20IR.pdf
_version_ 1747813000861974528