Cardiac modeling: coupling of elastomechanics and fluid dynamics
In silico modeling of the human heart and the simulation based on this are gaining an increased meaning in medicine and science. Therefore, computer models can be consulted to support a doctor’s diagnosis, to enhance clinical therapies or to plan a medical intervention.
For this reason, at the Institute of Biomedical Engineering an in house simulation tool has been created and is continuously improved since several years.
The aim of this research project is the further development of the existing simulation tool regarding hemodynamics with the desired outcome to provide recognition of diseases concerning the cardiovascular system. This is possible because fluid dynamics in the human heart are immediately affected by changes of the cardiac function.
Therefore, a coupling of hemodynamics, elastomechanics and circulatory system should accomplish a contribution to the improvement of medical diagnosis.
All Publications, sorted by years
Journal Articles (2)
Influence of pressure boundary condition definition on flow patterns in cardiac simulations.
In Modeling the Cardiac Function, 2022
Sequential Coupling Shows Minor Effects of Fluid Dynamics on Myocardial Deformation in a Realistic Whole-Heart Model..
In Frontiers in cardiovascular medicine, vol. 8, pp. 768548, 2021
The human heart is a masterpiece of the highest complexity coordinating multi-physics aspects on a multi-scale range. Thus, modeling the cardiac function to reproduce physiological characteristics and diseases remains challenging. Especially the complex simulation of the blood's hemodynamics and its interaction with the myocardial tissue requires a high accuracy of the underlying computational models and solvers. These demanding aspects make whole-heart fully-coupled simulations computationally highly expensive and call for simpler but still accurate models. While the mechanical deformation during the heart cycle drives the blood flow, less is known about the feedback of the blood flow onto the myocardial tissue. To solve the fluid-structure interaction problem, we suggest a cycle-to-cycle coupling of the structural deformation and the fluid dynamics. In a first step, the displacement of the endocardial wall in the mechanical simulation serves as a unidirectional boundary condition for the fluid simulation. After a complete heart cycle of fluid simulation, a spatially resolved pressure factor (PF) is extracted and returned to the next iteration of the solid mechanical simulation, closing the loop of the iterative coupling procedure. All simulations were performed on an individualized whole heart geometry. The effect of the sequential coupling was assessed by global measures such as the change in deformation and-as an example of diagnostically relevant information-the particle residence time. The mechanical displacement was up to 2 mm after the first iteration. In the second iteration, the deviation was in the sub-millimeter range, implying that already one iteration of the proposed cycle-to-cycle coupling is sufficient to converge to a coupled limit cycle. Cycle-to-cycle coupling between cardiac mechanics and fluid dynamics can be a promising approach to account for fluid-structure interaction with low computational effort. In an individualized healthy whole-heart model, one iteration sufficed to obtain converged and physiologically plausible results.
Conference Contributions (5)
Cardiac fluid dynamics based on immersed boundary method for application in hypertrophic cardiomyopathy.
In 7th International Conference on Computational & Mathematical Biomedical Engineering - CMBE2021, pp. 439 - 442, 2022
Computational models of the fluid dynamics in the human heart are a powerful tool to investigate disease mechanisms and their impact on the blood flow patterns. These models can for example be used to assess alterations occurring in hypertrophic cardiomyopathy, which is a genetic disease that increases the risk of sudden cardiac death. To overcome the challenges of a moving mesh approach, we modeled the movement of the endocardial surface based on an immersed boundary method. The verification on a simple moving 2D geometry proved plausible results. The application to the dis- eased, hypertrophic heart geometry confirmed that the computation of the mesh movement is made possible with this approach.
Influence of Geometrical Properties for the Calculation of a Pressure-Free Whole Heart Geometry.
In 14th WCCM-ECCOMAS Congress 2020, pp. 1-9, 2021
Individualized computer models of the geometry of the human heart are often based on mag- netic resonance images (MRI) or computed tomography (CT) scans. The stress distribution in the imaged state cannot be measured but needs to be estimated from the segmented geometry, e.g. by an iterative algorithm. As the convergence of this algorithm depends on different geometrical conditions, we system- atically studied their influence. Beside various shape alterations, we investigated the chamber volume, as well as the effect of material parameters. We found a marked influence of passive material parameters: increasing the model stiffness by a factor of ten halved the residual norm in the first iteration. Flat and concave areas led to a reduced robustness and convergence rate of the unloading algorithm. With this study, the geometric effects and modeling aspects governing the unloading algorithm’s convergence are identified and can be used as a basis for further improvement.
Fluid dynamics in the human heart: Altered vortex formation and wash-out in mitral regurgitation simulations.
In Current Directions in Biomedical Engineering, vol. 7(2) , pp. 199-202, 2021
Mitral regurgitation alters the flow conditions in the left ventricle. To account for quantitative changes and to investigate the behavior of different flow components, a realistic computational model of the whole human heart was employed in this study. While performing fluid dynamics simulations, a scalar transport equation was solved to analyze vortex formation and ventricular wash-out for different regurgitation severities. Additionally, a particle tracking algorithm was implemented to visualize single components of the blood flow. We confirmed a significantly lowered volume of the direct flow component as well as a higher vorticity in the diseased case.
A sequential coupling approach for fluid-structure interaction in a patient specific whole heart geometry.
In iHEART Congress – Modelling the Cardiac Function, 2021
In order to be used in a clinical context, numerical simulation tools have to strike a balance between accuracy and low computational effort. For re- producing the pumping function of the human heart numerically, the physical domains of cardiac continuum mechanics and fluid dynamics have a significant relevance. In this context, fluid-structure interaction between the heart muscle and the blood flow is particularly important: Myocardial tension development and wall deformation drive the blood flow. However, the degree to which the blood flow has a retrograde effect on the cardiac mechanics in this multi-physics problem remains unclear up to now. To address this question, we implemented a cycle-to-cycle coupling based on a finite element model of a patient-specific whole heart geometry. The deforma- tion of the cardiac wall over one heart cycle was computed using our mechanical simulation framework. A closed loop circulatory system model as part of the simulation delivered the chamber pressures. The displacement of the endo- cardial surfaces and the pressure courses of one cycle were used as boundary conditions for the fluid solver. After solving the Navier-Stokes equations, the relative pressure was extracted for all endocardial wall elements from the three dimensional pressure field. These local pressure deviations were subsequently returned to the next iteration of the continuum mechanical simulation, thus closing the loop of the iterative coupling procedure. Following this sequential coupling approach, we simulated three iterations of mechanic and fluid simulations. To characterize the convergence, we evaluated the time course of the normalized pressure field as well as the euclidean distance between nodes of the mechanic simulation in subsequent iterations. For the left ventricle (LV), the maximal euclidean distance of all endocardial wall nodes was smaller than 2mm between the first and second iteration. The maximal distance between the second and third iteration was 70μm, thus the limit of necessary cycles was already reached after two iterations. In future work, this iterative coupling approach will have to prove its abil- ity to deliver physiologically accurate results also for diseased heart models. Altogether, the sequential coupling approach with its low computational effort delivered promising results for modeling fluid-structure interaction in cardiac simulations.
Characterization of the Fluid Dynamic Pressure Field in the Human Heart as a Basis for Coupled Fluid-Structure Simulations.
In BMT 2020, vol. Poster Session, pp. 259, 2020