A three-dimensional computational analysis of fluid-structure interaction in the aortic valve

Numerical analysis of the aortic valve has mainly been focused on the closing behaviour during the diastolic phase rather than the kinematic opening and closing behaviour during the systolic phase of the cardiac cycle. Moreover, the fluid-structure interaction in the aortic valve system is most frequently ignored in numerical modelling. The effect of this interaction on the valve's behaviour during systolic functioning is investigated. The large differences in material properties of fluid and structure and the finite motion of the leaflets complicate blood-valve interaction modelling. This has impeded numerical analyses of valves operating under physiological conditions. A numerical method, known as the Lagrange multiplier based fictitious domain method, is used to describe the large leaflet motion within the computational fluid domain. This method is applied to a three-dimensional finite element model of a stented aortic valve. The model provides both the mechanical behaviour of the valve and the blood flow through it. Results show that during systole the leaflets of the stented valve appear to be moving with the fluid in an essentially kinematical process governed by the fluid motion.     

Modeling active muscle contraction in mitral valve leaflets during systole: a first approach

The present study addresses the effect of muscle activation contributions to mitral valve leaflet response during systole. State-of-art passive hyperelastic material modeling is employed in combination with a simple active stress part. Fiber families are assumed in the leaflets: one defined by the collagen and one defined by muscle activation. The active part is either assumed to be orthogonal to the collagen fibers or both orthogonal to and parallel with the collagen fibers (i.e. an orthotropic muscle fiber model). Based on data published in the literature and information herein on morphology, the size of the leaflet parts that contain muscle fibers is estimated. These parts have both active and passive materials, the remaining parts consist of passive material only. Several solid finite element analyses with different maximum activation levels are run. The simulation results are compared to corresponding echocardiography at peak systole for a porcine model. The physiologically correct flat shape of the closed valve is approached as the activation levels increase. The non-physiological bulging of the leaflet into the left atrium when using passive material models is reduced significantly. These results contribute to improved understanding of the physiology of the native mitral valve, and add evidence to the hypothesis that the mitral valve leaflets not are just passive elements moving as a result of hemodynamic pressure gradients in the left part of the heart.     

Experimental measurement of dynamic fluid shear stress on the ventricular surface of the aortic valve leaflet

Aortic valve (AV) calcification is a highly prevalent disease with serious impact on mortality and morbidity. The exact causes and mechanisms of AV calcification are unclear, although previous studies suggest that mechanical forces play a role. It has been clinically demonstrated that calcification preferentially occurs on the aortic surface of the AV. This is hypothesized to be due to differences in the mechanical environments on the two sides of the valve. It is thus necessary to characterize fluid shear forces acting on both sides of the leaflet to test this hypothesis. The current study is one of two studies characterizing dynamic shear stress on both sides of the AV leaflets. In the current study, shear stresses on the ventricular surface of the AV leaflets were measured experimentally on two prosthetic AV models with transparent leaflets in an in vitro pulsatile flow loop using two-component Laser Doppler Velocimetry (LDV). Experimental measurements were utilized to validate a theoretical model of AV ventricular surface shear stress based on the Womersley profile in a straight tube, with corrections for the opening angle of the valve leaflets. This theoretical model was applied to in vivo data based on MRI-derived volumetric flow rates and valve dimension obtained from the literature. Experimental results showed that ventricular surface shear stress was dominated by the streamwise component. The systolic shear stress waveform resembled a half-sinusoid during systole and peaks at 64–71 dyn/cm², and reversed in direction at the end of systole for 15–25?ms, and reached a significant negative magnitude of 40–51 dyn/cm². Shear stresses from the theoretical model applied to in vivo data showed that shear stresses peaked at 77–92 dyn/cm² and reversed in direction for substantial period of time (108–110?ms) during late systole with peak negative shear stress of 35–38 dyn/cm².     

Stress related collagen ultrastructure in human aortic valves--implications for tissue engineering

Understanding the response of tissue structures to mechanical stress is crucial for optimization of mechanical conditioning protocols in the field of heart valve tissue engineering. In heart valve tissue, it is unclear to what extent mechanical loading affects the collagen fibril morphology. To determine if local stress affects the collagen fibril morphology, in terms of fibril diameter, its distribution, and the fibril density, this was investigated in adult native human aortic valve leaflets. Transmission electron microscopy images of collagen fibrils were analyzed at three locations: the commissures, the belly, and the fixed edge of the leaflets. Subsequently, the mechanical behavior of human aortic valves was used in a computational model to predict the stress distribution in the valve leaflet during the diastolic phase of the cardiac cycle. The local stresses at the three locations were related to the collagen fibril morphology. The fibril diameter and density varied significantly between the measured locations, and appeared inversely related. The average fibril diameter increased from the fixed edge, to the belly, and to the commissures of the leaflets, while fibril density decreased. Interestingly, these differences corresponded well with the level of stress at the locations. The presented data showed that large tissue stress is associated with greater average fibril diameter, lower fibril density, and wider fibril size distribution compared with low stress locations in the leaflets. The findings here provide insight in the effect of mechanical loading on the collagen ultrastructure, and are valuable to improve conditioning protocols for tissue engineering.     

Finite element analysis of the mitral apparatus: annulus shape effect and chordal force distribution

This study presents a three-dimensional finite element model of the mitral apparatus using a hyperelastic transversely isotropic material model for the leaflets. The objectives of this study are to illustrate the effects of the annulus shape on the chordal force distribution and on the mitral valve response during systole, to investigate the role of the anterior secondary (strut) chordae and to study the influence of thickness of the leaflets on the leaflets stresses. Hence, analyses are conducted with a moving and fixed saddle shaped annulus and with and without anterior secondary chordae. We found that the tension in the secondary chordae represents 31% of the load carried by the papillary muscles. When removing the anterior secondary chordae, the tension in the primary anterior chordae is almost doubled, the displacement of the anterior leaflet toward the left atrium is also increased. The moving annulus configuration with an increasing annulus saddle height does not give significant changes in the chordal force distribution and in the leaflet stress compared to the fixed annulus. The results also show that the maximum principle stresses in the anterior leaflet are carried by the collagen fibers. The stresses calculated in the leaflets are very sensitive to the thickness employed.     

Experimental technique of measuring dynamic fluid shear stress on the aortic surface of the aortic valve leaflet

Aortic valve (AV) calcification is a highly prevalent disease with serious impact on mortality and morbidity. The exact cause and mechanism of the progression of AV calcification is unknown, although mechanical forces have been known to play a role. It is thus important to characterize the mechanical environment of the AV. In the current study, we establish a methodology of measuring shear stresses experienced by the aortic surface of the AV leaflets using an in vitro valve model and adapting the laser Doppler velocimetry (LDV) technique. The valve model was constructed from a fresh porcine aortic valve, which was trimmed and sutured onto a plastic stented ring, and inserted into an idealized three-lobed sinus acrylic chamber. Valve leaflet location was measured by obtaining the location of highest back-scattered LDV laser light intensity. The technique of performing LDV measurements near to biological surfaces as well as the leaflet locating technique was first validated in two phantom flow systems: (1) steady flow within a straight tube with AV leaflet adhered to the wall, and (2) steady flow within the actual valve model. Dynamic shear stresses were then obtained by applying the techniques on the valve model in a physiologic pulsatile flow loop. Results show that aortic surface shear stresses are low during early systole (<5 dyn/cm2) but elevated to its peak during mid to late systole at about 18-20 dyn/cm2. Low magnitude shear stress (<5 dyn/cm2) was observed during early diastole and dissipated to zero over the diastolic duration. Systolic shear stress was observed to elevate only with the formation of sinus vortex flow. The presented technique can also be used on other in vitro valve models such as congenitally geometrically malformed valves, or to investigate effects of hemodynamics on valve shear stress. Shear stress data can be used for further experiments investigating effects of fluid shear stress on valve biology, for conditioning tissue engineered AV, and to validate numerical simulations.     

A three-dimensional mechanical analysis of a stentless fibre-reinforced aortic valve prosthesis

Failure of bioprosthetic and synthetic three-leaflet valves has been shown to occur as a consequence of high tensile and bending stresses, acting on the leaflets during opening and closing. Moreover, in the stented prostheses, whether synthetic or biological, the absence of contraction of the aortic base, due to the rigid stent, causes the leaflets to be subjected to an unphysiological degree of flexure, which is related to calcification. It is shown that the absence of the stent, which gives a flexible aortic base and leaflet attachment, and leaflet fibre-reinforcement result in reduced stresses in the weaker parts of the leaflets in their closed configuration. It is postulated that this leads to a decrease of tears and perforations, which may result in a improved long-term behaviour. The effect of a flexible leaflet attachment and aortic base of a synthetic valve is investigated with a finite element model. Different fibre-reinforced structures are analysed with respect to the stresses that are likely to contribute to the failure of fibre-reinforced prostheses and compared with the results obtained for a stented prosthesis. Results show that for the stentless models a reduction of stresses up to 75% is obtained with respect to stented models with the same type of reinforcement.     

Straightening of curved pattern of collagen fibers under load controls aortic valve shape

The network of collagen fibers in the aortic valve leaflet is believed to play an important role in the strength and durability of the valve. However, in addition to its stress-bearing role, such a fiber network has the potential to produce functionally important shape changes in the closed valve under pressure load. We measured the average pattern of the collagen network in porcine aortic valve leaflets after staining for collagen. We then used finite element simulation to explore how this collagen pattern influences the shape of the closed valve. We observed a curved or bent pattern, with collagen fibers angled downward from the commissures toward the center of the leaflet to form a pattern that is concave toward the leaflet free edge. Simulations showed that these curved fiber trajectories straighten under pressure load, leading to functionally important changes in closed valve shape. Relative to a pattern of straight collagen fibers running parallel to the leaflet free edge, the concave pattern of curved fibers produces a closed valve with a 40% increase in central leaflet coaptation height and with decreased leaflet billow, resulting in a more physiological closed valve shape. Furthermore, simulations show that these changes in loaded leaflet shape reflect changes in leaflet curvature due to modulation of in-plane membrane stress resulting from straightening of the curved fibers. This effect appears to play an important role in normal valve function and may have important implications for the design of prosthetic and tissue engineered replacement valves.     

Experimental measurement of dynamic fluid shear stress on the aortic surface of the aortic valve leaflet

Aortic valve (AV) calcification is a highly prevalent disease with serious impact on mortality and morbidity. Although exact causes and mechanisms of AV calcification are unclear, previous studies suggest that mechanical forces play a role. Since calcium deposits occur almost exclusively on the aortic surfaces of AV leaflets, it has been hypothesized that adverse patterns of fluid shear stress on the aortic surface of AV leaflets promote calcification. The current study characterizes AV leaflet aortic surface fluid shear stresses using Laser Doppler velocimetry and an in vitro pulsatile flow loop. The valve model used was a native porcine valve mounted on a suturing ring and preserved using 0.15% glutaraldehyde solution. This valve model was inserted in a mounting chamber with sinus geometries, which is made of clear acrylic to provide optical access for measurements. To understand the effects of hemodynamics on fluid shear stress, shear stress was measured across a range of conditions: varying stroke volumes at the same heart rate and varying heart rates at the same stroke volume. Systolic shear stress magnitude was found to be much higher than diastolic shear stress magnitude due to the stronger flow in the sinuses during systole, reaching up to 20 dyn/cm² at mid-systole. Upon increasing stroke volume, fluid shear stresses increased due to stronger sinus fluid motion. Upon increasing heart rate, fluid shear stresses decreased due to reduced systolic duration that restricted the formation of strong sinus flow. Significant changes in the shear stress waveform were observed at 90 beats/min, most likely due to altered leaflet dynamics at this higher heart rate. Overall, this study represents the most well-resolved shear stress measurements to date across a range of conditions on the aortic side of the AV. The data presented can be used for further investigation to understand AV biological response to shear stresses.     

Aortic valve leaflet mechanical properties facilitate diastolic valve function

This work was concerned with the numerical simulation of the behaviour of aortic valves whose material can be modelled as non-linear elastic anisotropic. Linear elastic models for the valve leaflets with parameters used in previous studies were compared with hyperelastic models, incorporating leaflet anisotropy with pronounced stiffness in the circumferential direction through a transverse isotropic model. The parameters for the hyperelastic models were obtained from fits to results of orthogonal uniaxial tensile tests on porcine aortic valve leaflets. The computational results indicated the significant impact of transverse isotropy and hyperelastic effects on leaflet mechanics; in particular, increased coaptation with peak values of stress and strain in the elastic limit. The alignment of maximum principal stresses in all models follows approximately the coarse collagen fibre distribution found in aortic valve leaflets. The non-linear elastic leaflets also demonstrated more evenly distributed stress and strain which appears relevant to long-term scaffold stability and mechanotransduction.     

A synthetic fiber-reinforced stentless heart valve

There is strong evidence that failure of bioprosthetic and synthetic valves occurs as a consequence of high tensile and bending stresses, acting on the leaflets during opening and closing. In stented prostheses, whether synthetic or biological, the absence of contraction of the aortic base causes the leaflets to be subjected to an unphysiological degree of flexure, which is also related to calcification. However, a stentless synthetic valve, which has a flexible aorta base, can be a good alternative for stented synthetic valves. Moreover, fiber-reinforcement is assumed to lead to a decrease of tears and perforation as a result of reduced stresses in the weaker parts of the leaflets in their closed configuration. The manufacturing method for a stentless, fiber-reinforced, synthetic valve is presented. Prototypes are tested in a pulse duplicator system. The results show that the mean systolic pressure difference is very low, while the high regurgitation (up to 26%) is probably caused by a too small coaptation area of the leaflets.     

Osteopontin controls endothelial cell migration in vitro and in excised human valvular tissue from patients with calcific aortic stenosis and controls

Calcific aortic stenosis (CAS) is a pathological condition of the aortic valve characterized by dystrophic calcification of the valve leaflets. Despite the high prevalence and mortality associated with CAS, little is known about its pathogenetic mechanisms. Characterized by progressive dystrophic calcification of the valve leaflets, the early stages of aortic valve degeneration are similar to the active inflammatory process of atherosclerosis including endothelial disruption, inflammatory cell infiltration, lipid deposition, neo-vascularization and calcification. In the vascular system, the endothelium is an important regulator of physiological and pathological conditions; however, the contribution of endothelial dysfunction to valvular degeneration at the cellular and molecular level has received little attention. Endothelial cell (EC) activation and neo-vascularization of the cusps characterizes all stages of aortic valvular degeneration from aortic sclerosis to aortic stenosis. Here we reported the role of osteopontin (OPN) in the regulation of EC activation in vitro and in excised tissue from CAS patients and controls. OPN is an important pro-angiogenic factor in several pathologies. High levels of OPN have been demonstrated in both tissue and plasma of patients with aortic valve sclerosis and stenosis. The characterization of valvular ECs as a cellular target for OPN will help us uncover the pathogenesis of aortic valve degeneration and stenosis, opening new perspectives for the prevention and therapy of this prevalent disease.     

Non-linear rotation-free shell finite-element models for aortic heart valves

Hyperelastic material models have been incorporated in the rotation-free, large deformation, shell finite element (FE) formulation of (Stolarski et al., 2013) and applied to dynamic simulations of aortic heart valve. Two models used in the past in analysis of such problem i.e. the Saint-Venant and May-Newmann–Yin (MNY) material models have been considered and compared. Uniaxial tests for those constitutive equations were performed to verify the formulation and implementation of the models. The issue of leaflets interactions during the closing of the heart valve at the end of systole is considered. The critical role of using non-linear anisotropic model for proper dynamic response of the heart valve especially during the closing phase is demonstrated quantitatively. This work contributes an efficient FE framework for simulating biological tissues and paves the way for high-fidelity flow structure interaction simulations of native and bioprosthetic aortic heart valves.     

Interlayer micromechanics of the aortic heart valve leaflet

While the mechanical behaviors of the fibrosa and ventricularis layers of the aortic valve (AV) leaflet are understood, little information exists on their mechanical interactions mediated by the GAG-rich central spongiosa layer. Parametric simulations of the interlayer interactions of the AV leaflets in flexure utilized a tri-layered finite element (FE) model of circumferentially oriented tissue sections to investigate inter-layer sliding hypothesized to occur. Simulation results indicated that the leaflet tissue functions as a tightly bonded structure when the spongiosa effective modulus was at least 25 % that of the fibrosa and ventricularis layers. Novel studies that directly measured transmural strain in flexure of AV leaflet tissue specimens validated these findings. Interestingly, a smooth transmural strain distribution indicated that the layers of the leaflet indeed act as a bonded unit, consistent with our previous observations (Stella and Sacks in J Biomech Eng 129:757–766, 2007) of a large number of transverse collagen fibers interconnecting the fibrosa and ventricularis layers. Additionally, when the tri-layered FE model was refined to match the transmural deformations, a layer-specific bimodular material model (resulting in four total moduli) accurately matched the transmural strain and moment-curvature relations simultaneously. Collectively, these results provide evidence, contrary to previous assumptions, that the valve layers function as a bonded structure in the low-strain flexure deformation mode. Most likely, this results directly from the transverse collagen fibers that bind the layers together to disable physical sliding and maintain layer residual stresses. Further, the spongiosa may function as a general dampening layer while the AV leaflets deforms as a homogenous structure despite its heterogeneous architecture.     

Validation of a fluid-structure interaction model of a heart valve using the dynamic mesh method in fluent

Simulations of coupled problems such as fluid-structure interaction (FSI) are becoming more and more important for engineering purposes. This is particularly true when modeling the aortic valve, where the FSI between the blood and the valve determines the valve movement and the valvular hemodynamics. Nevertheless only a few studies are focusing on the opening and closing behavior during the ejection phase (systole). In this paper, we present the validation of a FSI model using the dynamic mesh method of Fluent for the two-dimensional (2D) simulation of mechanical heart valves during the ejection phase of the cardiac cycle. The FSI model is successfully validated by comparing simulation results to experimental data obtained from in vitro studies using a CCD camera.     

The congenital bicuspid aortic valve can experience high-frequency unsteady shear stresses on its leaflet surface

The bicuspid aortic valve (BAV) is a common congenital malformation of the aortic valve (AV) affecting 1% to 2% of the population. The BAV is predisposed to early degenerative calcification of valve leaflets, and BAV patients constitute 50% of AV stenosis patients. Although evidence shows that genetic defects can play a role in calcification of the BAV leaflets, we hypothesize that drastic changes in the mechanical environment of the BAV elicit pathological responses from the valve and might be concurrently responsible for early calcification. An in vitro model of the BAV was constructed by surgically manipulating a native trileaflet porcine AV. The BAV valve model and a trileaflet AV (TAV) model were tested in an in vitro pulsatile flow loop mimicking physiological hemodynamics. Laser Doppler velocimetry was used to make measurements of fluid shear stresses on the leaflet of the valve models using previously established methodologies. Furthermore, particle image velocimetry was used to visualize the flow fields downstream of the valves and in the sinuses. In the BAV model, flow near the leaflets and fluid shear stresses on the leaflets were much more unsteady than for the TAV model, most likely due to the moderate stenosis in the BAV and the skewed forward flow jet that collided with the aorta wall. This additional unsteadiness occurred during mid- to late-systole and was composed of cycle-to-cycle magnitude variability as well as high-frequency fluctuations about the mean shear stress. It has been demonstrated that the BAV geometry can lead to unsteady shear stresses under physiological flow and pressure conditions. Such altered shear stresses could play a role in accelerated calcification in BAVs.     

Validation of a Fluid–Structure Interaction Model of a Heart Valve using the Dynamic Mesh Method in Fluent

Simulations of coupled problems such as fluid–structure interaction (FSI) are becoming more and more important for engineering purposes. This is particularly true when modeling the aortic valve, where the FSI between the blood and the valve determines the valve movement and the valvular hemodynamics. Nevertheless only a few studies are focusing on the opening and closing behavior during the ejection phase (systole). In this paper, we present the validation of a FSI model using the dynamic mesh method of Fluent for the two-dimensional (2D) simulation of mechanical heart valves during the ejection phase of the cardiac cycle. The FSI model is successfully validated by comparing simulation results to experimental data obtained from in vitro studies using a CCD camera.     

Patient-specific modeling of biomechanical interaction in transcatheter aortic valve deployment

The objective of this study was to develop a patient-specific computational model to quantify the biomechanical interaction between the transcatheter aortic valve (TAV) stent and the stenotic aortic valve during TAV intervention. Finite element models of a patient-specific stenotic aortic valve were reconstructed from multi-slice computed tomography (MSCT) scans, and TAV stent deployment into the aortic root was simulated. Three initial aortic root geometries of this patient were analyzed: (a) aortic root geometry directly reconstructed from MSCT scans, (b) aortic root geometry at the rapid right ventricle pacing phase, and (c) aortic root geometry with surrounding myocardial tissue. The simulation results demonstrated that stress, strain, and contact forces of the aortic root model directly reconstructed from MSCT scans were significantly lower than those of the model at the rapid ventricular pacing phase. Moreover, the presence of surrounding myocardium slightly increased the mechanical responses. Peak stresses and strains were observed around the calcified regions in the leaflets, suggesting the calcified leaflets helped secure the stent in position. In addition, these elevated stresses induced during TAV stent deployment indicated a possibility of tissue tearing and breakdown of calcium deposits, which might lead to an increased risk of stroke. The potential of paravalvular leak and occlusion of coronary ostia can be evaluated from simulated post-deployment aortic root geometries. The developed computational models could be a valuable tool for pre-operative planning of TAV intervention and facilitate next generation TAV device design.     

DIASTOLIC OPENING OF THE AORTIC VALVE (AV) IN A CASE OF AORTIC INSUFFICIENCY DUE TO AV FENESTRATION

A patient with severe aortic insufficiency due to fenestration of the non-coronary aortic valve leaflet is described. A preoperative echocardiogram demonstrated early closure of the mitral valve and early diastolic separation of the aortic valve leaflets. These findings disappeared after partial surgical correction and subsequent hemodynamic improvement. Premature opening of the aortic valve is common in severe aortic insufficiency.     

Improved prediction of the collagen fiber architecture in the aortic heart valve

Living tissues show an adaptive response to mechanical loading by changing their internal structure and morphology. Understanding this response is essential for successful tissue engineering of load-bearing structures, such as the aortic valve. In this study, mechanically induced remodeling of the collagen architecture in the aortic valve was investigated. It was hypothesized that, in uniaxially loaded regions, the fibers aligned with the tensile principal stretch direction. For biaxial loading conditions, on the other hand, it was assumed that the collagen fibers aligned with directions situated between the principal stretch directions. This hypothesis has already been applied successfully to study collagen remodeling in arteries. The predicted fiber architecture represented a branching network and resembled the macroscopically visible collagen bundles in the native leaflet. In addition, the complex biaxial mechanical behavior of the native valve could be simulated qualitatively with the predicted fiber directions. The results of the present model might be used to gain further insight into the response of tissue engineered constructs during mechanical conditioning.