Patient-Specific Finite Element Modeling of the Blood Flow in the Left Ventricle of a Human Heart

• Heart disease is the leading cause of death in the world. Therefore, numerous studies are undertaken to identify indicators which can be applied to discover cardiac dysfunctions at an early age. Among others, the fluid dynamics of the blood flow (hemodymanics) is considered to contain relevant information related to abnormal performance of the heart.This thesis presents a robust framework for numerical simulation of the fluid dynamics of the blood flow in the left ventricle of a human heart and the fluid-structure interaction of the blood and the aortic leaflets.We first describe a patient-specific model for simulating the intraventricular blood flow. The motion of the endocardial wall is extracted from data acquired with medical imaging and we use the incompressible Navier-Stokes equations to model the hemodynamics within the chamber. We set boundary conditions to model the opening and closing of the mitral and aortic valves respectively, and we apply a stabilized Arbitrary Lagrangian-Eulerian (ALE) space-time finite element method to simulate the blood flow. Even though it is difficult to collect in-vivo data for validation, the available data and results from other simulation models indicate that our approach possesses the potential and capability to provide relevant information about the intraventricular blood flow.To further demonstrate the robustness and clinical feasibility of our model, a semi-automatic pathway from 4D cardiac ultrasound imaging to patient-specific simulation of the blood flow in the left ventricle is developed. The outcome is promising and further simulations and analysis of large data sets are planned.In order to enhance our solver by introducing additional features, the fluid solver is extended by embedding different geometrical prototypes of both a native and a mechanical aortic valve in the outflow area of the left ventricle.Both, the contact as well as the fluid-structure interaction, are modeled as a unified continuum problem using conservation laws for mass and momentum. To use this ansatz for simulating the valvular dynamics is unique and has the expedient properties that the whole problem can be described with partial different equations and the same numerical methods for discretization are applicable.All algorithms are implemented in the high performance computing branch of Unicorn, which is part of the open source software framework FEniCS-HPC. The strong advantage of implementing the solvers in an open source software is the accessibility and reproducibility of the results which enhance the prospects of developing a method with clinical relevance.

Ämnesord

NATURVETENSKAP  -- Matematik -- Beräkningsmatematik (hsv//swe)

NATURAL SCIENCES  -- Mathematics -- Computational Mathematics (hsv//eng)

Nyckelord

Finite element method

Arbitrary Lagrangian-Eulerian method

Fluid-Structure interaction

Contact model

parallel algorithm

blood flow

left ventricle

aortic valves

patient-specific heart model

Publikations- och innehållstyp

vet (ämneskategori)

dok (ämneskategori)