ABSTRACT
Due to advances in medical imaging, computational fluid dynamics algorithms and high performance computing, computer simulation is developing into an important tool for understanding the relationship between cardiovascular diseases and intraventricular blood flow. The field of cardiac flow simulation is challenging and highly interdisciplinary. We apply a computational framework for automated solutions of partial differential equations using Finite Element Methods where any mathematical description directly can be translated to code. This allows us to develop a cardiac model where specific properties of the heart such as fluid-structure interaction of the aortic valve can be added in a modular way without extensive efforts. In previous work, we simulated the blood flow in the left ventricle of the heart. In this paper, we extend this model by placing prototypes of both a native and a mechanical aortic valve in the outflow region of the left ventricle. Numerical simulation of the blood flow in the vicinity of the valve offers the possibility to improve the treatment of aortic valve diseases as aortic stenosis (narrowing of the valve opening) or regurgitation (leaking) and to optimize the design of prosthetic heart valves in a controlled and specific way. The fluid-structure interaction and contact problem are formulated in a unified continuum model using the conservation laws for mass and momentum and a phase function. The discretization is based on an Arbitrary Lagrangian-Eulerian space-time finite element method with streamline diffusion stabilization, and it is implemented in the open source software Unicorn which shows near optimal scaling up to thousands of cores. Computational results are presented to demonstrate the capability of our framework.
ABSTRACT
The combination of medical imaging with computational fluid dynamics (CFD) has enabled the study of 3-D blood flow on a patient-specific level. However, with models based on gated high-resolution data, the study of transient flows, and any model implementation into routine cardiac care, is challenging. This paper presents a novel pathway for patient-specific CFD modelling of the left ventricle (LV), using 4-D transthoracic echocardiography (TTE) as input modality. To evaluate the clinical usability, two sub-studies were performed. First, a robustness evaluation was performed, where repeated models with alternating input variables were generated for six subjects and changes in simulated output quantified. Second, a validation study was carried out, where the pathway accuracy was evaluated against pulsed-wave Doppler (100 subjects), and 2-D through-plane phase-contrast magnetic resonance imaging measurements over seven intraventricular planes (6 subjects). The robustness evaluation indicated a model deviation of <12%, with highest regional and temporal deviations at apical segments and at peak systole, respectively. The validation study showed an error of <11% (velocities <10 cm/s) for all subjects, with no significant regional or temporal differences observed. With the patient-specific pathway shown to provide robust output with high accuracy, and with the pathway dependent only on 4-D TTE, the method has a high potential to be used within future clinical studies on 3-D intraventricular flow patterns. To this, future model developments in the form of e.g., anatomically accurate LV valves may further enhance the clinical value of the simulations.
