skip to main content


Title: Computational Methods for Fluid-Structure Interaction Simulation of Heart Valves in Patient-Specific Left Heart Anatomies
Given the complexity of human left heart anatomy and valvular structures, the fluid–structure interaction (FSI) simulation of native and prosthetic valves poses a significant challenge for numerical methods. In this review, recent numerical advancements for both fluid and structural solvers for heart valves in patient-specific left hearts are systematically considered, emphasizing the numerical treatments of blood flow and valve surfaces, which are the most critical aspects for accurate simulations. Numerical methods for hemodynamics are considered under both the continuum and discrete (particle) approaches. The numerical treatments for the structural dynamics of aortic/mitral valves and FSI coupling methods between the solid Ωs and fluid domain Ωf are also reviewed. Future work toward more advanced patient-specific simulations is also discussed, including the fusion of high-fidelity simulation within vivo measurements and physics-based digital twining based on data analytics and machine learning techniques.  more » « less
Award ID(s):
1946202 2019077
NSF-PAR ID:
10341868
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Fluids
Volume:
7
Issue:
3
ISSN:
2311-5521
Page Range / eLocation ID:
94
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Venous valves are bicuspidal valves that ensure that blood in veins only flows back to the heart. To prevent retrograde blood flow, the two intraluminal leaflets meet in the center of the vein and occlude the vessel. In fluid‐structure interaction (FSI) simulations of venous valves, the large structural displacements may lead to mesh deteriorations and entanglements, causing instabilities of the solver and, consequently, the numerical solution to diverge. In this paper, we propose an arbitrary Lagrangian‐Eulerian (ALE) scheme for FSI simulations designed to solve these instabilities. A monolithic formulation for the FSI problem is considered, and due to the complexity of the operators, the exact Jacobian matrix is evaluated using automatic differentiation. The method relies on the introduction of a staggered in time velocity to improve stability, and on fictitious springs to model the contact force of the valve leaflets. Because the large structural displacements may compromise the quality of the fluid mesh as well, a smoother fluid displacement, obtained with the introduction of a scaling factor that measures the distance of a fluid element from the valve leaflet tip, guarantees that there are no mesh entanglements in the fluid domain. To further improve stability, a streamline upwind Petrov‐Galerkin (SUPG) method is employed. The proposed ALE scheme is applied to a two‐dimensional (2D) model of a venous valve. The presented simulations show that the proposed method deals well with the large structural displacements of the problem, allowing a reconstruction of the valve behavior in both the opening and closing phase.

     
    more » « less
  2. Abstract

    A new general contact model is proposed for preventing inter‐leaflet penetration of bio‐prosthetic heart valves (BHV) at the end of the systole, which has the advantage of applying kinematic constraints directly and creating smooth free edges. At the end of each time step, the impenetrability constraints and momentum exchange between the impacting bodies are applied separately based on the coefficient of restitution. The contact method is implemented in a rotation‐free, large deformation, and thin shell finite‐element (FE) framework based on loop's subdivision surfaces. A nonlinear, anisotropic material model for a BHV is employed which uses Fung‐elastic constitutive laws for in‐plane and bending responses, respectively. The contact model is verified and validated against several benchmark problems. For a BHV‐specific validation, the computed strains on different regions of a BHV under constant pressure are compared with experimentally measured data. Finally, dynamic simulations of BHV under physiological pressure waveform are performed for symmetrical and asymmetrical fiber orientations incorporating the new contact model and compared with the penalty contact method. The proposed contact model provides the coaptation area of a functioning BHV during the closing phase for both of the fiber orientations. Our results show that fiber orientation affects the dynamic of leaflets during the opening and closing phases. A swirling motion for the BHV with asymmetrical fiber orientation is observed, similar to experimental data. To include the fluid effects, fluid–structure interaction (FSI) simulation of the BHV is performed and compared to the dynamic results.

     
    more » « less
  3. Abstract Computational fluid dynamics (CFD) modeling of left ventricle (LV) flow combined with patient medical imaging data has shown great potential in obtaining patient-specific hemodynamics information for functional assessment of the heart. A typical model construction pipeline usually starts with segmentation of the LV by manual delineation followed by mesh generation and registration techniques using separate software tools. However, such approaches usually require significant time and human efforts in the model generation process, limiting large-scale analysis. In this study, we propose an approach toward fully automating the model generation process for CFD simulation of LV flow to significantly reduce LV CFD model generation time. Our modeling framework leverages a novel combination of techniques including deep-learning based segmentation, geometry processing, and image registration to reliably reconstruct CFD-suitable LV models with little-to-no user intervention.1 We utilized an ensemble of two-dimensional (2D) convolutional neural networks (CNNs) for automatic segmentation of cardiac structures from three-dimensional (3D) patient images and our segmentation approach outperformed recent state-of-the-art segmentation techniques when evaluated on benchmark data containing both magnetic resonance (MR) and computed tomography(CT) cardiac scans. We demonstrate that through a combination of segmentation and geometry processing, we were able to robustly create CFD-suitable LV meshes from segmentations for 78 out of 80 test cases. Although the focus on this study is on image-to-mesh generation, we demonstrate the feasibility of this framework in supporting LV hemodynamics modeling by performing CFD simulations from two representative time-resolved patient-specific image datasets. 
    more » « less
  4. Faithful, accurate, and successful cardiac biomechanics and electrophysiological simulations require patient-specific geometric models of the heart. Since the cardiac geometry consists of highly-curved boundaries, the use of high-order meshes with curved elements would ensure that the various curves and features present in the cardiac geometry are well-captured and preserved in the corresponding mesh. Most other existing mesh generation techniques require computer-aided design files to represent the geometric boundary, which are often not available for biomedical applications. Unlike such methods, our technique takes a high-order surface mesh, generated from patient medical images, as input and generates a high-order volume mesh directly from the curved surface mesh. In this paper, we use our direct high-order curvilinear tetrahedral mesh generation method [1] to generate several second-order cardiac meshes. Our meshes include the left ventricle myocardia of a healthy heart and hearts with dilated and hypertrophic cardiomyopathy. We show that our high-order cardiac meshes do not contain inverted elements and are of sufficiently high quality for use in cardiac finite element simulations. 
    more » « less
  5. Abstract Transcatheter aortic valve replacement (TAVR) first received FDA approval for high-risk surgical patients in 2011 and has been approved for low-risk surgical patients since 2019. It is now the most common type of aortic valve replacement, and its use continues to accelerate. Computer modeling and simulation (CM&S) is a tool to aid in TAVR device design, regulatory approval, and indication in patient-specific care. This study introduces a computational fluid-structure interaction (FSI) model of TAVR with Medtronic’s CoreValve Evolut R device using the immersed finite element-difference (IFED) method. We perform dynamic simulations of crimping and deployment of the Evolut R, as well as device behavior across the cardiac cycle in a patient-specific aortic root anatomy reconstructed from computed tomography (CT) image data. These IFED simulations, which incorporate biomechanics models fit to experimental tensile test data, automatically capture the contact within the device and between the self-expanding stent and native anatomy. Further, we apply realistic driving and loading conditions based on clinical measurements of human ventricular and aortic pressures and flow rates to demonstrate that our Evolut R model supports a physiological diastolic pressure load and provides informative clinical performance predictions. 
    more » « less