Tomographic PIV in a model of the left ventricle: 3D flow past biological and mechanical heart valves

Left ventricular flow is intrinsically complex, three-dimensional and unsteady. Its features are susceptible to cardiovascular pathology and treatment, in particular to surgical interventions involving the valves (mitral valve replacement). To improve our understanding of intraventricular fluid mechanics and the impact of various types of prosthetic valves thereon, we have developed a custom-designed versatile left ventricular phantom with anatomically realistic moving left ventricular membrane. A biological, a tilting disc and a bileaflet valve (in two different orientations) were mounted in the mitral position and tested under the same settings. To investigate 3D flow within the phantom, a four-view tomographic particle image velocimetry setup has been implemented. The results compare side-by-side the evolution of the 3D flow topology, vortical structures and kinetic energy in the left ventricle domain during the cardiac cycle. Except for the tilting disc valve, all tested prosthetic valves induced a crossed flow path, where the outflow crosses the inflow path, passing under the mitral valve. The biological valve shows a strong jet with a peak velocity about twice as high compared to all mechanical heart valves, which makes it easier to penetrate deeply into the cavity. Accordingly, the peak kinetic energy in the left ventricle in case of the biological valve is about four times higher than the mechanical heart valves. We conclude that the tomographic particle imaging velocimetry setup provides a useful ground truth measurement of flow features and allows a comparison of the effects of different valve types on left ventricular flow patterns.     

Comparison of the hemodynamic and thrombogenic performance of two bileaflet mechanical heart valves using a CFD/FSI model

The hemodynamic and the thrombogenic performance of two commercially available bileaflet mechanical heart valves (MHVs)--the ATS Open Pivot Valve (ATS) and the St. Jude Regent Valve (SJM), was compared using a state of the art computational fluid dynamics-fluid structure interaction (CFD-FSI) methodology. A transient simulation of the ATS and SJM valves was conducted in a three-dimensional model geometry of a straight conduit with sudden expansion distal the valves, including the valve housing and detailed hinge geometry. An aortic flow waveform (60 beats/min, cardiac output 4 l/min) was applied at the inlet. The FSI formulation utilized a fully implicit coupling procedure using a separate solver for the fluid problem (FLUENT) and for the structural problem. Valve leaflet excursion and pressure differences were calculated, as well as shear stress on the leaflets and accumulated shear stress on particles released during both forward and backward flow phases through the open and closed valve, respectively. In contrast to the SJM, the ATS valve opened to less than maximal opening angle. Nevertheless, maximal and mean pressure gradients and velocity patterns through the valve orifices were comparable. Platelet stress accumulation during forward flow indicated that no platelets experienced a stress accumulation higher than 35 dyne x s/cm2, the threshold for platelet activation (Hellums criterion). However, during the regurgitation flow phase, 0.81% of the platelets in the SJM valve experienced a stress accumulation higher than 35 dyne x s/cm2, compared with 0.63% for the ATS valve. The numerical results indicate that the designs of the ATS and SJM valves, which differ mostly in their hinge mechanism, lead to different potential for platelet activation, especially during the regurgitation phase. This numerical methodology can be used to assess the effects of design parameters on the flow induced thrombogenic potential of blood recirculating devices.     

Intraventricular flow patterns and stasis in the LVAD-assisted heart

Left ventricular assist device (LVAD) support disrupts the natural blood flow path through the heart, introducing flow patterns associated with thrombosis, especially in the presence of medical devices. The aim of this study was to quantitatively evaluate the flow patterns in the left ventricle (LV) of the LVAD-assisted heart, with a focus on alterations in vortex development and stasis. Particle image velocimetry of a LVAD-supported LV model was performed in a mock circulatory loop. In the Pre-LVAD flow condition, a vortex ring initiating from the LV base migrated toward the apex during diastole and remained in the LV by the end of ejection. During LVAD support, vortex formation was relatively unchanged although vortex circulation and kinetic energy increased with LVAD speed, particularly in systole. However, as pulsatility decreased and aortic valve opening ceased, a region of fluid stasis formed near the left ventricular outflow tract. These findings suggest that LVAD support does not substantially alter vortex dynamics unless cardiac function is minimal. The altered blood flow introduced by the LVAD results in stasis adjacent to the LV outflow tract, which increases the risk of thrombus formation in the heart.     

Effect of prosthetic mitral valve geometry and orientation on flow dynamics in a model human left ventricle

Pulsatile flow dynamics through bileaflet (St Jude and Duromedics), tilting disc (Bjork-Shiley and Omniscience), caged ball (Starr-Edwards), pericardial (Edwards) and porcine (Carpentier-Edwards) mitral valves in a model human left ventricle (LV) were studied. The model human ventricle, obtained from an in situ diastolic casting, was incorporated into a mock circulatory system. Measurements were made at various heart rates and flow rates. These included the transvalvular pressure drop and regurgitation in percent and cm3 beat-1. The effect of valve geometry and the orientation of the valve with respect to the valve annulus was analyzed using a flow visualization technique. Qualitative flow visualization study indicates certain preferred orientations for the tilting disc and bileaflet valve prostheses in order to obtain a smooth washout of flow in the LV chamber.     

The effect of vortex formation on left ventricular filling and mitral valve efficiency

A new mechanism for quantifying the filling energetics in the left ventricle (LV) and past mechanical heart valves (MHV) is identified and presented. This mechanism is attributed to vortex formation dynamics past MHV leaflets. Recent studies support the conjecture that the natural healthy left ventricle (LV) performs in an optimum, energy-preserving manner by redirecting the flow with high efficiency. Yet to date, no quantitative proof has been presented. The present work provides quantitative results and validation of a theory based on the dynamics of vortex ring formation, which is governed by a critical formation number (FN) that corresponds to the dimensionless time at which the vortex ring has reached its maximum circulation content, in support of this hypothesis. Herein, several parameters (vortex ring circulation, vortex ring energy, critical FN, hydrodynamic efficiencies, vortex ring propagation speed) have been quantified and presented as a means of bridging the physics of vortex formation in the LV. In fact, the diastolic hydrodynamic efficiencies were found to be 60, 41, and 29%, respectively, for the porcine, anti-anatomical, and anatomical valve configurations. This assessment provides quantitative proof of vortex formation, which is dependent of valve design and orientation, being an important flow characteristic and associated to LV energetics. Time resolved digital particle image velocimetry with kilohertz sampling rate was used to study the ejection of fluid into the LV and resolve the spatiotemporal evolution of the flow. The clinical significance of this study is quantifying vortex formation and the critical FN that can potentially serve as a parameter to quantify the LV filling process and the performance of heart valves.     

Device Thrombogenicity Emulator (DTE) − Design optimization methodology for cardiovascular devices: A study in two bileaflet MHV designs

Patients who receive prosthetic heart valve (PHV) implants require mandatory anticoagulation medication after implantation due to the thrombogenic potential of the valve. Optimization of PHV designs may facilitate reduction of flow-induced thrombogenicity and reduce or eliminate the need for post-implant anticoagulants. We present a methodology entitled Device Thrombogenicty Emulator (DTE) for optimizing the thrombo-resistance performance of PHV by combining numerical and experimental approaches. Two bileaflet mechanical heart valves (MHV) designs, St. Jude Medical (SJM) and ATS, were investigated by studying the effect of distinct flow phases on platelet activation. Transient turbulent and direct numerical simulations (DNS) were conducted, and stress loading histories experienced by the platelets were calculated along flow trajectories. The numerical simulations indicated distinct design dependent differences between the two valves. The stress loading waveforms extracted from the numerical simulations were programmed into a hemodynamic shearing device (HSD), emulating the flow conditions past the valves in distinct ‘hot-spot’ flow regions that are implicated in MHV thrombogenicity. The resultant platelet activity was measured with a modified prothrombinase assay, and was found to be significantly higher in the SJM valve, mostly during the regurgitation phase. The experimental results were in excellent agreement with the calculated platelet activation potential. This establishes the utility of the DTE methodology for serving as a test bed for evaluating design modifications for achieving better thrombogenic performance for such devices.     

Effect of turbulent models on left ventricle diastolic flow patterns simulation

Abstract

Vortex structures, as one of the most important features of cardiac flow, have a crucial impact on the left ventricle function and pathological conditions. These swirling flows are closely related to the presence of turbulence in left ventricle which is investigated in the current study. Using an extended model of the left heart, including a fluid-structure interaction (FSI) model of the mitral valve with a realistic geometry, the effect of using two numerical turbulent models, k-ε and Spalart-Allmaras (SA), on diastolic flow patterns is studied and compared with results from laminar flow model. As a result of the higher dissipation rate in turbulent models (k-ε and SA), vortices are larger and stronger in the laminar flow model. Comparing E/A ratio in the three models (Laminar, k-ε, and SA) with experimental data from healthy subjects, it is concluded that the results from k-ε model are more accurate.     

Heterogeneity of the myocardium. Function of the left and right ventricle under normal and pathological conditions

A number of important differences can be found between the left ventricle (LV) and right ventricle (RV) of the heart under physiological conditions. In anatomy, the most important is probably the architecture of the atrioventricular valve and its annulus. The LV has a mitral valve (with two cusps) and a firm annulus, while the RV has a tricuspid valve with a greater total area, but relatively small cuspid areas, and an elastic annulus. The difference in the blood supply is important. Owing to high intramural pressure, the coronary flow in the wall of the LV occurs only during the diastole; in the RV it is limited only in the presence of a significant increase in intracavitary pressure. The LV myocardium is functionally "accustomed" to short-term marked changes in the systolic load (in extreme static exercise the arterial pressure rises for a short time to three times the normal value), while the RV is adapted to changes in the diastolic load (marked filling changes associated with deep breathing, for instance). The difference in the response to a long-term volume load is difficult to evaluate: between a defect of the interatrial septum and aortic insufficiency there are too many differences. A long-term pressure load seems to be tolerated better by the right ventricle: patients with severe pulmonary stenosis and a pressure six times higher than the physiological value have lived 25 years and patients with isolated corrected L-transposition of the great arteries can reach 35 years without any signs of impaired RV function.(ABSTRACT TRUNCATED AT 250 WORDS)     

Defect in oxidative phosphorylation in LV papillary muscle mitochondria of patients undergoing mitral valve replacement

Mitochondria play a pivotal role in cellular metabolism, especially in energy production. Myocardial function depends on adenosine triphosphate (ATP) supplied by oxidation of several substrates. In the adult heart, this energy is obtained primarily from fatty acid oxidation through oxidative phosphorylation (OXPHOS). With this in view, we studied OXPHOS, Total-ATPase and cytochrome content in the mitochondria of the left ventricular (LV) papillary muscles in excised mitral valves of patients who underwent mitral valve replacement (MVR). The mitochondrial OXPHOS, cytochrome content and ATPase activity were studied in 70 patients (ranging from 22 to 40 years) operated on for mitral valve disease. Control study includes 25 normal mitral valves removed at necropsy from patients who died of extracardiac causes. In the presence of glutamate and succinate as substrates, the rate of mitochondrial oxygen consumption was significantly lower in LV papillary muscles of pathological mitral valves (P<0.001) by using with and without addition of ADP. The ADP/O ratio indices for glutamate and succinate were not significantly affected. Using glutamate as substrate, respiratory control index was significantly raised (P<0.05) as compared with control. A significant reduction of total cytochrome content and ATPase activity (P<0.001) was noted in LV papillary muscles of patients operated for mitral valve disease. Our results showed that OXPHOS, cytochromes 'a', 'b', 'c+c(1)' and ATP activity are significantly impaired in LV papillary muscles in patients with pathological mitral valve. Cardiac mitochondrial oxygen consumption is a very valuable tool to investigate the regulation of cardiac mitochondrial energy metabolism. There is increasing evidence that mitochondrial diseases, such as mitochondrial cardiomyopathy, valvular disease and some myopathies, can be responsive to treatment with metabolic intermediates such as coenzyme Q(10), thiamine, prednisone, and vitamin therapy.     

Evaluation of a transient,simultaneous, arbitrary Lagrange–Euler based multi-physics method for simulating the mitral heart valve

A transient multi-physics model of the mitral heart valve has been developed, which allows simultaneous calculation of fluid flow and structural deformation. A recently developed contact method has been applied to enable simulation of systole (the stage when blood pressure is elevated within the heart to pump blood to the body). The geometry was simplified to represent the mitral valve within the heart walls in two dimensions. Only the mitral valve undergoes deformation. A moving arbitrary Lagrange–Euler mesh is used to allow true fluid–structure interaction (FSI). The FSI model requires blood flow to induce valve closure by inducing strains in the region of 10–20%. Model predictions were found to be consistent with existing literature and will undergo further development.     

An Approach to the Stepwise Management of Severe Mitral Regurgitation with Optimal Cardiac Pacemaker Function

Right ventricular apical pacing may cause or worsen mitral regurgitation (MR). Potential mechanisms for this adverse sequelae include intraventricular dyssynchrony, altered papillary muscle function, pacing-induced cardiomyopathy with left ventricular dilation, and annular dilation. In contrast, biventricular (BiV) pacing may improve MR presumably by opposing the negative effects. Whether or not left ventricular lead location is important in treating mitral regurgitation in patients with pacemakers is unknown.We report a case of severe MR and left ventricular (LV) systolic failure in a patient with right ventricular pacing. Multiple potential etiologies for the worsening valve function were noted, and a stepwise iterative optimizing scheme that included basal lateral LV pacing improved mitral valve function and ameliorated heart failure symptoms.     

Doppler echocardiographic estimation of left ventricular end-diastolic pressure after MI in rats

The spectral Doppler mitral flow pattern, alone or combined with tissue Doppler mitral annulus velocity, can be used to predict left ventricular (LV) filling pressure in humans, whereas invasive hemodynamic measurements are still required in the rat. This study was undertaken to assess whether LV end-diastolic pressure (LVEDP) can be estimated using Doppler echocardiography in the rat after myocardial infarction (MI). Thirty-seven rats (23 rats with MI after left coronary artery ligation and 14 sham-operated rats) were evaluated 3 mo after surgery with echo-Doppler and invasive hemodynamic measurements. Pulse wave spectral Doppler at the mitral valve tip was used to measure the E wave, the E wave deceleration time (DT), and the A wave; spectral Doppler tissue imaging was used to measure the early diastolic lateral mitral annulus velocity (E(a)). We found weak correlations between LVEDP and the peak velocity of the early mitral inflow (E), E/peak velocity of the late mitral inflow, and DT, and strong correlations with E(a) and especially with E/E(a) [R(2) = 0.89, LVEDP (in mmHg) = 0.987E/E(a) - 4.229]. Longitudinal followup of a subgroup of rats with MI revealed a marked rise of E/E(a) between days 7 and 21 in rats with heart failure only. We conclude that Doppler echocardiography can be used for serial assessment of LV diastolic function in rats with MI.     

Left ventricular endocardial longitudinal and transverse changes during isovolumic contraction and relaxation: a challenge

Left ventricular (LV) longitudinal and transverse geometric changes during isovolumic contraction and relaxation are still controversial. This confusion is compounded by traditional definitions of these phases of the cardiac cycle. High-resolution sonomicrometry studies might clarify these issues. Crystals were implanted in six sheep at the LV apex, fibrous trigones, lateral and posterior mitral annulus, base of the aortic right coronary sinus, anterior and septal endocardial wall, papillary muscle tips, and edge of the anterior and posterior mitral leaflets. Changes in distances were time related to LV and aortic pressures and to mitral valve opening. At the beginning of isovolumic contraction, while the mitral valve was still open, the LV endocardial transverse diameter started to shorten while the endocardial longitudinal diameter increased. During isovolumic relaxation, while the mitral valve was closed, LV transverse diameter started to increase while the longitudinal diameter continued to decrease. These findings are inconsistent with the classic definitions of the phases of the cardiac cycle.     

Two-dimensional intraventricular flow pattern visualization using the image-based computational fluid dynamics

The image-based computational fluid dynamics (IB-CFD) technique, as the combination of medical images and the CFD method, is utilized in this research to analyze the left ventricle (LV) hemodynamics. The research primarily aims to propose a semi-automated technique utilizing some freely available and commercial software packages in order to simulate the LV hemodynamics using the IB-CFD technique. In this research, moreover, two different physiological time-resolved 2D models of a patient-specific LV with two different types of aortic and mitral valves, including the orifice-type valves and integrated with rigid leaflets, are adopted to visualize the process of developing intraventricular vortex formation and propagation. The blood flow pattern over the whole cardiac cycle of two models is also compared to investigate the effect of utilizing different valve types in the process of the intraventricular vortex formation. Numerical findings indicate that the model with integrated valves can predict more complex intraventricular flow that can match better the physiological flow pattern in comparison to the orifice-type model.     

Wave intensity analysis of left atrial mechanics and energetics in anesthetized dogs

The left atrium (LA) acts as a booster pump during late diastole, generating the Doppler transmitral A wave and contributing incrementally to left ventricular (LV) filling. However, after volume loading and in certain disease states, LA contraction fills the LV less effectively, and retrograde flow (i.e., the Doppler Ar wave) into the pulmonary veins increases. The purpose of this study was to provide an energetic analysis of LA contraction to clarify the mechanisms responsible for changes in forward and backward flow. Wave intensity analysis was performed at the mitral valve and a pulmonary vein orifice. As operative LV stiffness increased with progressive volume loading, the reflection coefficient (i.e., energy of reflected wave/energy of incident wave) also increased. This reflected wave decelerated the forward movement of blood through the mitral valve and was transmitted through the LA, accelerating retrograde blood flow in the pulmonary veins. Although total LA work increased with volume loading, the forward hydraulic work decreased and backward hydraulic work increased. Thus wave reflection due to increased LV stiffness accounts for the decrease in the A wave and the increase in the Ar wave measured by Doppler.     

Development of aortic and mitral valve continuity in the human embryonic heart

The anatomic relationship of the aortic and mitral valves is a useful landmark in assessing congenital heart malformations. The atrioventricular and semilunar valve regions originate in widely separated parts of the early embryonic heart tube, and the process by which the normal fibrous continuity between the aortic and mitral valves is acquired has not been clearly defined. The development of the aortic and mitral valve relationship was studied in normal human embryos in the Carnegie Embryological Collection, and specimens of Carnegie stages 13, 15, 17, 19, and 23, prepared as serial histologic sections cut in the sagittal plane, were selected for reconstruction. In stage 13, the atrioventricular valve area is separated from the semilunar valve area by the large bend between the atrioventricular and outflow-tract components of the single lumen heart tube created by the left interventricular sulcus. In stages 15 and 17, the aortic valve rotates into a position near the atrioventricular valves with development of four chambers and a double circulation. In stage 19, there is fusion of aortic and mitral endocardial cushion material along the endocardial surface of the interventricular flange, and this relationship is maintained in subsequent stages. Determination of three-dimensional Cartesian coordinates of the midpoints of valve positions shows that, while there is growth of intervalvular distances up to stage 17, the aortic to mitral distance is essentially unchanged thereafter. During the period studied, the left ventricle increases in length over threefold. The relative lack of growth in the saddle-shaped fold between the atrioventricular and outflow tract components of the heart, contrasting with the rapid growth of the outwardly convex components of most of the atrial and ventricular walls, may be attributed to the different mechanical properties of the two configurations. It is postulated that the pathogenesis of congenital heart malformations, which characteristically have failure of development of aortic and mitral valve continuity, may involve abnormalities of rotation of the aortic region or malpositioning of the fold in the heart tube.     

Dynamic modelling of prosthetic chorded mitral valves using the immersed boundary method

Current artificial heart valves either have limited lifespan or require the recipient to be on permanent anticoagulation therapy. In this paper, effort is made to assess a newly developed bileaflet valve prosthesis made of synthetic flexible leaflet materials, whose geometry and material properties are based on those of the native mitral valve, with a view to providing superior options for mitral valve replacement. Computational analysis is employed to evaluate the geometric and material design of the valve, by investigation of its mechanical behaviour and unsteady flow characteristics. The immersed boundary (IB) method is used for the dynamic modelling of the large deformation of the valve leaflets and the fluid-structure interactions. The IB simulation is first validated for the aortic prosthesis subjected to a hydrostatic loading. The predicted displacement fields by IB are compared with those obtained using ANSYS, as well as with experimental measurements. Good quantitative agreement is obtained. Moreover, known failure regions of aortic prostheses are identified. The dynamic behaviour of the valve designs is then simulated under four physiological pulsatile flows. Experimental pressure gradients for opening and closure of the valves are in good agreement with IB predictions for all flow rates for both aortic and mitral designs. Importantly, the simulations predicted improved physiological haemodynamics for the novel mitral design. Limitation of the current IB model is also discussed. We conclude that the IB model can be developed to be an extremely effective dynamic simulation tool to aid prosthesis design.     

A fully coupled fluid-structure interaction model of the secondary lymphatic valve

Abstract

The secondary lymphatic valve is a bi-leaflet structure frequent throughout collecting vessels that serves to prevent retrograde flow of lymph. Despite its vital function in lymph flow and apparent importance in disease development, the lymphatic valve and its associated fluid dynamics have been largely understudied. The goal of this work was to construct a physiologically relevant computational model of an idealized rat mesenteric lymphatic valve using fully coupled fluid-structure interactions to investigate the relationship between three-dimensional flow patterns and stress/deformation within the valve leaflets. The minimum valve resistance to flow, which has been shown to be an important parameter in effective lymphatic pumping, was computed as 268?g/mm⁴?s. Hysteretic behavior of the lymphatic valve was confirmed by comparing resistance values for a given transvalvular pressure drop during opening and closing. Furthermore, eddy structures were present within the sinus adjacent to the valve leaflets in what appear to be areas of vortical flow; the eddy structures were characterized by non-zero velocity values (up to ～4?mm/s) in response to an applied unsteady transvalvular pressure. These modeling capabilities present a useful platform for investigating the complex interplay between soft tissue motion and fluid dynamics of lymphatic valves and contribute to the breadth of knowledge regarding the importance of biomechanics in lymphatic system function.     

Numerical Modeling of Intraventricular Flow during Diastole after Implantation of BMHV

This work presents a numerical simulation of intraventricular flow after the implantation of a bileaflet mechanical heart valve at the mitral position. The left ventricle was simplified conceptually as a truncated prolate spheroid and its motion was prescribed based on that of a healthy subject. The rigid leaflet rotation was driven by the transmitral flow and hence the leaflet dynamics were solved using fluid-structure interaction approach. The simulation results showed that the bileaflet mechanical heart valve at the mitral position behaved similarly to that at the aortic position. Sudden area expansion near the aortic root initiated a clockwise anterior vortex, and the continuous injection of flow through the orifice resulted in further growth of the anterior vortex during diastole, which dominated the intraventricular flow. This flow feature is beneficial to preserving the flow momentum and redirecting the blood flow towards the aortic valve. To the best of our knowledge, this is the first attempt to numerically model intraventricular flow with the mechanical heart valve incorporated at the mitral position using a fluid-structure interaction approach. This study facilitates future patient-specific studies.