Effect of bending rigidity in a dynamic model of a polyurethane prosthetic mitral valve

On the Lagrangian-Eulerian Coupling in the Immersed Finite Element/Difference Method

The immersed boundary (IB) method is a non-body conforming approach to fluid-structure interaction (FSI) that uses an Eulerian description of the momentum, viscosity, and incompressibility of a coupled fluid-structure system and a Lagrangian description of the deformations, stresses, and resultant forces of the immersed structure. Integral transforms with Dirac delta function kernels couple the Eulerian and Lagrangian variables, and in practice, discretizations of these integral transforms use regularized delta function kernels. Many different kernel functions have been proposed, but prior numerical work investigating the impact of the choice of kernel function on the accuracy of the methodology has often been limited to simplified test cases or Stokes flow conditions that may not reflect the method's performance in applications, particularly at intermediate-to-high Reynolds numbers, or under different loading conditions. This work systematically studies the effect of the choice of regularized delta function in several fluid-structure interaction benchmark tests using the immersed finite element/difference (IFED) method, which is an extension of the IB method that uses a finite element structural discretization combined with a Cartesian grid finite difference method for the incompressible Navier-Stokes equations. Whereas the conventional IB method spreads forces from the nodes of the structural mesh and interpolates velocities to those nodes, the IFED formulation evaluates the regularized delta function on a collection of interaction points that can be chosen to be denser than the nodes of the Lagrangian mesh. This opens the possibility of using structural discretizations with wide node spacings that would produce gaps in the Eulerian force in nodally coupled schemes (e.g., if the node spacing is comparable to or broader than the support of the regularized delta functions). Earlier work with this methodology suggested that such coarse structural meshes can yield improved accuracy for shear-dominated cases and, further, found that accuracy improves when the structural mesh spacing is increased. However, these results were limited to simple test cases that did not include substantial pressure loading on the structure. This study investigates the effect of varying the relative mesh widths of the Lagrangian and Eulerian discretizations in a broader range of tests. Our results indicate that kernels satisfying a commonly imposed even-odd condition require higher resolution to achieve similar accuracy as kernels that do not satisfy this condition. We also find that narrower kernels are more robust, in the sense that they yield results that are less sensitive to relative changes in the Eulerian and Lagrangian mesh spacings, and that structural meshes that are substantially coarser than the Cartesian grid can yield high accuracy for shear-dominated cases but not for cases with large normal forces. We verify our results in a large-scale FSI model of a bovine pericardial bioprosthetic heart valve in a pulse duplicator.

An Immersed Interface Method for Discrete Surfaces

Fluid-structure systems occur in a range of scientific and engineering applications. The immersed boundary (IB) method is a widely recognized and effective modeling paradigm for simulating fluid-structure interaction (FSI) in such systems, but a difficulty of the IB formulation of these problems is that the pressure and viscous stress are generally discontinuous at fluid-solid interfaces. The conventional IB method regularizes these discontinuities, which typically yields low-order accuracy at these interfaces. The immersed interface method (IIM) is an IB-like approach to FSI that sharply imposes stress jump conditions, enabling higher-order accuracy, but prior applications of the IIM have been largely restricted to numerical methods that rely on smooth representations of the interface geometry. This paper introduces an immersed interface formulation that uses only a C 0 representation of the immersed interface, such as those provided by standard nodal Lagrangian finite element methods. Verification examples for models with prescribed interface motion demonstrate that the method sharply resolves stress discontinuities along immersed boundaries while avoiding the need for analytic information about the interface geometry. Our results also demonstrate that only the lowest-order jump conditions for the pressure and velocity gradient are required to realize global second-order accuracy. Specifically, we demonstrate second-order global convergence rates along with nearly second-order local convergence in the Eulerian velocity field, and between first- and second-order global convergence rates along with approximately first-order local convergence for the Eulerian pressure field. We also demonstrate approximately second-order local convergence in the interfacial displacement and velocity along with first-order local convergence in the fluid traction along the interface. As a demonstration of the method's ability to tackle more complex geometries, the present approach is also used to simulate flow in a patient-averaged anatomical model of the inferior vena cava, which is the large vein that carries deoxygenated blood from the lower extremities back to the heart. Comparisons of the general hemodynamics and wall shear stress obtained by the present IIM and a body-fitted discretization approach show that the present method yields results that are in good agreement with those obtained by the body-fitted approach.

Some Effects of Different Constitutive Laws on FSI Simulation for the Mitral Valve

In this paper, three different constitutive laws for mitral leaflets and two laws for chordae tendineae are selected to study their effects on mitral valve dynamics with fluid-structure interaction. We first fit these three mitral leaflet constitutive laws and two chordae tendineae laws with experimental data. The fluid-structure interaction is implemented in an immersed boundary framework with finite element extension for solid, that is the hybrid immersed boundary/finite element(IB/FE) method. We specifically compare the fluid-structure results of different constitutive laws since fluid-structure interaction is the physiological loading environment. This allows us to look at the peak jet velocity, the closure regurgitation volume, and the orifice area. Our numerical results show that different constitutive laws can affect mitral valve dynamics, such as the transvalvular flow rate, closure regurgitation and the orifice area, while the differences in fiber strain and stress are insignificant because all leaflet constitutive laws are fitted to the same set of experimental data. In addition, when an exponential constitutive law of chordae tendineae is used, a lower closure regurgitation flow is observed compared to that of a linear material model. In conclusion, combining numerical dynamic simulations and static experimental tests, we are able to identify suitable constitutive laws for dynamic behaviour of mitral leaflets and chordae under physiological conditions.

On the chordae structure and dynamic behaviour of the mitral valve

We develop a fluid-structure interaction (FSI) model of the mitral valve (MV) that uses an anatomically and physiologically realistic description of the MV leaflets and chordae tendineae. Three different chordae models-complex, 'pseudo-fibre' and simplified chordae-are compared to determine how different chordae representations affect the dynamics of the MV. The leaflets and chordae are modelled as fibre-reinforced hyperelastic materials, and FSI is modelled using an immersed boundary-finite element method. The MV model is first verified under static boundary conditions against the commercial finite element software ABAQUS and then used to simulate MV dynamics under physiological pressure conditions. Interesting flow patterns and vortex formulation are observed in all three cases. To quantify the highly complex system behaviour resulting from FSI, an energy budget analysis of the coupled MV FSI model is performed. Results show that the complex and pseudo-fibre chordae models yield good valve closure during systole but that the simplified chordae model leads to poorer leaflet coaptation and an unrealistic bulge in the anterior leaflet belly. An energy budget analysis shows that the MV models with complex and pseudo-fibre chordae have similar energy distribution patterns but the MV model with the simplified chordae consumes more energy, especially during valve closing and opening. We find that the complex chordae and pseudo-fibre chordae have similar impact on the overall MV function but that the simplified chordae representation is less accurate. Because a pseudo-fibre chordal structure is easier to construct and less computationally intensive, it may be a good candidate for modelling MV dynamics or interaction between the MV and heart in patient-specific applications.

Hybrid finite difference/finite element immersed boundary method

The immersed boundary method is an approach to fluid-structure interaction that uses a Lagrangian description of the structural deformations, stresses, and forces along with an Eulerian description of the momentum, viscosity, and incompressibility of the fluid-structure system. The original immersed boundary methods described immersed elastic structures using systems of flexible fibers, and even now, most immersed boundary methods still require Lagrangian meshes that are finer than the Eulerian grid. This work introduces a coupling scheme for the immersed boundary method to link the Lagrangian and Eulerian variables that facilitates independent spatial discretizations for the structure and background grid. This approach uses a finite element discretization of the structure while retaining a finite difference scheme for the Eulerian variables. We apply this method to benchmark problems involving elastic, rigid, and actively contracting structures, including an idealized model of the left ventricle of the heart. Our tests include cases in which, for a fixed Eulerian grid spacing, coarser Lagrangian structural meshes yield discretization errors that are as much as several orders of magnitude smaller than errors obtained using finer structural meshes. The Lagrangian-Eulerian coupling approach developed in this work enables the effective use of these coarse structural meshes with the immersed boundary method. This work also contrasts two different weak forms of the equations, one of which is demonstrated to be more effective for the coarse structural discretizations facilitated by our coupling approach.

Modelling mitral valvular dynamics-current trend and future directions

Dysfunction of mitral valve causes morbidity and premature mortality and remains a leading medical problem worldwide. Computational modelling aims to understand the biomechanics of human mitral valve and could lead to the development of new treatment, prevention and diagnosis of mitral valve diseases. Compared with the aortic valve, the mitral valve has been much less studied owing to its highly complex structure and strong interaction with the blood flow and the ventricles. However, the interest in mitral valve modelling is growing, and the sophistication level is increasing with the advanced development of computational technology and imaging tools. This review summarises the state-of-the-art modelling of the mitral valve, including static and dynamics models, models with fluid-structure interaction, and models with the left ventricle interaction. Challenges and future directions are also discussed.

Modeling Patient-Specific Deformable Mitral Valves

Medical imaging has advanced enormously over the last few decades, revolutionizing patient diagnostics and care. At the same time, additive manufacturing has emerged as a means of reproducing physical shapes and models previously not possible. In combination, they have given rise to 3-dimensional (3D) modeling, an entirely new technology for physicians. In an era in which 3D imaging has become a standard for aiding in the diagnosis and treatment of cardiac disease, this visualization now can be taken further by bringing the patient's anatomy into physical reality as a model. The authors describe the generalized process of creating a model of cardiac anatomy from patient images and their experience creating patient-specific dynamic mitral valve models. This involves a combination of image processing software and 3D printing technology. In this article, the complexity of 3D modeling is described and the decision-making process for cardiac anesthesiologists is summarized. The management of cardiac disease has been altered with the emergence of 3D echocardiography, and 3D modeling represents the next paradigm shift.

Immersed boundary-finite element model of fluid-structure interaction in the aortic root

It has long been recognized that aortic root elasticity helps to ensure efficient aortic valve closure, but our understanding of the functional importance of the elasticity and geometry of the aortic root continues to e-volve as increasingly detailed in vivo imaging data become available. Herein, we describe a fluid-structure interaction model of the aortic root, including the aortic valve leaflets, the sinsuses of Valsalva, the aortic annulus, and the sinotubular junction, that employs a version of Peskin's immersed boundary (IB) method with a finite element (FE) description of the structural elasticity. As in earlier work, we use a fiber-based model of the valve leaflets, but this study extends earlier IB models of the aortic root by employing an incompressible hyperelastic model of the mechanics of the sinuses and ascending aorta using a constitutive law fit to experimental data from human aortic root tissue. In vivo pressure loading is accounted for by a backward displacement method that determines the unloaded configurations of the root model. Our model yields realistic cardiac output at physiological pressures, with low transvalvular pressure differences during forward flow, minimal regurgitation during valve closure, and realistic pressure loads when the valve is closed during diastole. Further, results from high-resolution computations indicate that although the detailed leaflet and root kinematics show some grid sensitivity, our IB model of the aortic root nonetheless produces essentially grid-converged flow rates and pressures at practical grid spacings for the high-Reynolds number flows of the aortic root. These results thereby clarify minimum grid resolutions required by such models when used as stand-alone models of the aortic valve as well as when used to provide models of the outflow valves in models of left ventricular fluid dynamics.

Emerging Trends in Heart Valve Engineering: Part IV. Computational Modeling and Experimental Studies

In this final portion of an extensive review of heart valve engineering, we focus on the computational methods and experimental studies related to heart valves. The discussion begins with a thorough review of computational modeling and the governing equations of fluid and structural interaction. We then move onto multiscale and disease specific modeling. Finally, advanced methods related to in vitro testing of the heart valves are reviewed. This section of the review series is intended to illustrate application of computational methods and experimental studies and their interrelation for studying heart valves.

A finite strain nonlinear human mitral valve model with fluid-structure interaction

A computational human mitral valve (MV) model under physiological pressure loading is developed using a hybrid finite element immersed boundary method, which incorporates experimentally-based constitutive laws in a three-dimensional fluid-structure interaction framework. A transversely isotropic material constitutive model is used to characterize the mechanical behaviour of the MV tissue based on recent mechanical tests of healthy human mitral leaflets. Our results show good agreement, in terms of the flow rate and the closing and opening configurations, with measurements from in vivo magnetic resonance images. The stresses in the anterior leaflet are found to be higher than those in the posterior leaflet and are concentrated around the annulus trigons and the belly of the leaflet. The results also show that the chordae play an important role in providing a secondary orifice for the flow when the valve opens. Although there are some discrepancies to be overcome in future work, our simulations show that the developed computational model is promising in mimicking the in vivo MV dynamics and providing important information that are not obtainable by in vivo measurements.

Quasi-static image-based immersed boundary-finite element model of left ventricle under diastolic loading

Finite stress and strain analyses of the heart provide insight into the biomechanics of myocardial function and dysfunction. Herein, we describe progress toward dynamic patient-specific models of the left ventricle using an immersed boundary (IB) method with a finite element (FE) structural mechanics model. We use a structure-based hyperelastic strain-energy function to describe the passive mechanics of the ventricular myocardium, a realistic anatomical geometry reconstructed from clinical magnetic resonance images of a healthy human heart, and a rule-based fiber architecture. Numerical predictions of this IB/FE model are compared with results obtained by a commercial FE solver. We demonstrate that the IB/FE model yields results that are in good agreement with those of the conventional FE model under diastolic loading conditions, and the predictions of the LV model using either numerical method are shown to be consistent with previous computational and experimental data. These results are among the first to analyze the stress and strain predictions of IB models of ventricular mechanics, and they serve both to verify the IB/FE simulation framework and to validate the IB/FE model. Moreover, this work represents an important step toward using such models for fully dynamic fluid-structure interaction simulations of the heart. © 2014 The Authors. International Journal for Numerical Methods in Engineering published by John Wiley & Sons, Ltd.