Fluid-structure interaction simulation of artificial textile reinforced aortic heart valve: Validation with an in-vitro test

Prosthetic heart valves deployed in the left heart (aortic and mitral) are subjected to harsh hemodynamical conditions. Most of the tissue engineered heart valves have been developed for the low pressure pulmonary position because of the difficulties in fabricating a mechanically strong valve, able to withstand the systemic circulation. This necessitates the use of reinforcing scaffolds, resulting in a tissue-engineered textile reinforced tubular aortic heart valve. Therefore, to better design these implants, material behaviour of the composite, valve kinematics and its hemodynamical response need to be evaluated. Experimental assessment can be immensely time consuming and expensive, paving way for numerical studies. In this work, the material properties obtained using the previously proposed multi-scale numerical method for textile composites was evaluated for its accuracy. An in silico immersed boundary (IB) fluid structure interaction (FSI) simulation emulating the in vitro experiment was set-up to evaluate and compare the geometric orifice area and flow rate for one beat cycle. Results from the in silico FSI simulation were found to be in good coherence with the in vitro test during the systolic phase, while mean deviation of approximately 9% was observed during the diastolic phase of a beat cycle. Merits and demerits of the in silico IB-FSI method for the presented case study has been discussed with the advantages outweighing the drawbacks, indicating the potential towards an effective use of this framework in the development and analysis of heart valves.     

Replacement of the Aortic and Mitral Valves Using the Starr-Edwards Ball-Valve Prosthesis: A Report of 50 Cases

The aortic and mitral valves were replaced in 50 patients at the University of Alberta Hospital using the Starr-Edwards ball-valve prosthesis. The basis of the selection of 20 patients for isolated aortic valve replacement and 27 for mitral valve replacement using this type of prosthesis is presented, and the techniques of insertion of the aortic and mitral valve are described in detail. Of the 27 patients in whom the mitral valve was replaced by the Starr-Edwards prosthesis six died within 30 days of surgery and two after discharge from hospital at two and a half and four months, respectively. Left atrial thrombosis was the cause of death in four of these patients. In 20 patients in whom the aortic valve was replaced, four died in hospital and two died more than 30 days after returning home. Three of these six patients died from bleeding—the result of the use of anticoagulants. The difficulty in assessing whether or not anticoagulants are needed following replacement by a Starr-Edwards prosthesis is considered. It is felt, in our present state of knowledge, that anticoagulants should be used following mitral valve replacement but are probably not essential following replacement of the aortic valve. Two patients survived replacement of both aortic and mitral valves and have been followed up 18 months and seven months, respectively.     

Design of a system for the accelerated loading of heart valve prostheses

In this paper an equipment is described for the loading of heart valve prostheses under physiological pressure conditions at a frequency of 10 Hz. The system consists essentially of two reservoirs between which a housing is mounted for holding the valve prosthesis. The reservoirs are partly filled with liquid. The physiological pressure variation across the valve is obtained by pressure control within the two reservoirs. A phase difference between the pressures in the two reservoirs compensates for the mass-inertia effects which normally occur at these high frequencies. The system without a valve has been analysed on the basis of simplified relationships between pressure and flow. The predicted values for the phase difference between the flow and the pressure curves within the valve housing, have been verified experimentally for various values of phase difference and amplitude ratios of the pressure variations within the reservoirs. The agreement between theory and experiment is fair. For the system with a valve the experimentally observed patterns closely resemble the theoretically predicted ones. From experiments with a Bj?rk-Shiley ( 21ABP ) and a Hancock (242-A21) valve prosthesis it is concluded that the valves open and close completely and that the pressure and flow patterns around the valves mimic the essential features of the in-vivo signals.     

The late results of cardiac valve replacement using autologous fascia lata

Over an eight-month period beginning in November 1969, 53 patients received 63 fascia lata heart valves at the Toronto General Hospital. The late results of this form of valve substitution are reviewed. The fascia used to fashion the tricuspid valve underwent progressive thickening and contracture and this process caused failure of the prosthesis within months of insertion. The mitral fascial prosthesis failed in a similar manner although the process took longer. The aortic fascia lata valve, however, has not shown progressive thickening after 3½ years and it has so far retained its functional integrity. Indeed, we have been impressed by the excellent clinical results and absence of complications such as thromboembolism although anticoagulation has not been used. We therefore consider that fascia lata valves offer a suitable alternative to other forms of aortic valve substitution, but are unsatisfactory for tricuspid or mitral valve replacement.     

Disc erosion in Models 103 and 104 of Beall mitral valve prostheses

Three cases of severe disc variance and erosion of the Teflon-disc Beall mitral valve prosthesis (Models 103 and 104) are reported. In two patients, the Beall mitral valves were excised and replaced with two Bj?rk-Shiley mitral valves. The remaining patient did not survive, and at autopsy, the lens was found at the aortic bifurcation level. Because of this potentially lethal complication, careful follow-up of patients with Beall mitral valve prostheses (Models 103 and 104) is recommended.     

Patient-specific simulation of a stentless aortic valve implant: the impact of fibres on leaflet performance

In some cases of aortic valve leaflet disease, the implant of a stentless biological prosthesis represents an excellent option for aortic valve replacement (AVR). In particular, if compared with the implant of mechanical valves, it provides a more physiological haemodynamic performance and a reduced thrombogeneticity, avoiding the use of anticoagulants. The clinical outcomes of AVR are strongly dependent on an appropriate choice of both prosthesis size and replacement technique, which is, at present, strictly related to surgeon's experience and skill. This represents the motivation for patient-specific finite element analysis able to virtually reproduce stentless valve implantation. With the aim of performing reliable patient-specific simulations, we remark that, on the one hand, it is not well established in the literature whether bioprosthetic leaflet tissue is isotropic or anisotropic; on the other hand, it is of fundamental importance to incorporate an accurate material model to realistically predict post-operative performance. Within this framework, using a novel computational methodology to simulate stentless valve implantation, we test the impact of using different material models on both the stress pattern and post-operative coaptation parameters (i.e. coaptation area, length and height). As expected, the simulation results suggest that the material properties of the valve leaflets affect significantly the post-operative prosthesis performance.     

Numerical simulation of steady turbulent flow through trileaflet aortic heart valves—II. Results on five models

Turbulent flow simulations are run for five aortic trileaflet valve geometries, ranging from a valve leaflet orifice area of 1.1 cm² (Model A1—very stenotic) to 5.0 cm² (Model A5—natural valve). The simulated data compares well with experimental measurements made downstream of various aortic trileaflet valves by Woo (PhD Thesis, 1984). The location and approximate width and length of recirculation regions are correctly predicted. The less stenotic valve models reattach at the end of the aortic sinus region, 1.1 diameters downstream of the valve. The central jet exiting the less stenotic valve models is not significantly different from fully developed flow, and therefore recovers very quickly downstream of the reattachment point. The more stenotic valves disturb the flow to a greater degree, generating recirculation regions large enough to escape the sinuses and reattach further downstream. Peak turbulent shear stress values downstream of the aortic valve models which approximated prosthetic valves are 125 and 300 N m⁻², very near experimental observations of 150 to 350 N m⁻². The predicted Reynolds stress profiles also present the correct shape, a double peak profile, with the location of the peak occuring at the location of maximum velocity gradient, which occurs near the recirculation region. The pressure drop across model A2 (leaflet orifice area 1.6 cm²) is 20 mmHg at 1.6 diameters downstream. This compares well with values ranging from 19.5 to 26.2 mmHg for valves of similar orifice areas. The pressure drop decreases with decreasing valve stenosis, to a negligible value across the least stenotic valve model. Based on the good agreement between experimental measurements of velocity, shear stress and pressure drop, compared to the simulated data, the model has the potential to be a valuable tool in the analysis of heart valve designs.     

Telocytes in human heart valves

Valve interstitial cells (VICs) are responsible for maintaining the structural integrity and dynamic behaviour of the valve. Telocytes (TCs), a peculiar type of interstitial cells, have been recently identified by Popescu's group in epicardium, myocardium and endocardium (visit www.telocytes.com ). The presence of TCs has been identified in atria, ventricles and many other tissues and organ, but not yet in heart valves. We used transmission electron microscopy and immunofluorescence methods (double labelling for CD34 and c‐kit, or vimentin, or PDGF Receptor‐β) to provide evidence for the existence of TCs in human heart valves, including mitral valve, tricuspid valve and aortic valve. TCs are found in both apex and base of heart valves, with a similar density of 27–28 cells/mm² in mitral valve, tricuspid valve and aortic valve. Since TCs are known for the participation in regeneration or repair biological processes, it remains to be determined how TCs contributes to the valve attempts to re‐establish normal structure and function following injury, especially a complex junction was found between TCs and a putative stem (progenitor) cell.     

Osteogenesis inducers promote distinct biological responses in aortic and mitral valve interstitial cells

Both aortic and mitral valves calcify in pathological conditions; however, the prevalence of aortic valve calcification is high whereas mitral valve leaflet calcification is somewhat rare. Patterns of valvular calcification may differ due to valvular architecture, but little is known to that effect. In this study, we investigated the intrinsic osteogenic differentiation potential of aortic versus mitral valve interstitial cells provided minimal differentiation conditions. For the assessment of calcification at the cellular level, we used classic inducers of osteogenesis in stem cells: β-glycerophosphate (β-Gly), dexamethasone (Dex), and ascorbate (Asc). In addition to proteomic analyses, osteogenic markers and calcium precipitates were evaluated across treatments of aortic and mitral valve cells. The combination of β-Gly, Asc, and Dex induced aortic valve interstitial cells to synthesize extracellular matrix, overexpress osteoblastic markers, and deposit calcium. However, no strong evidence showed the calcification of mitral valve interstitial cells. Mitral cells mainly responded to Asc and Dex by cell activation. These findings provide a deeper understanding of the physiological properties of aortic and mitral valves and tendencies for calcific changes within each valve type, contributing to the development of future therapeutics for heart valve diseases.     

In-vitro fluid dynamic characteristics of the abiomed trileaflet heart valve prosthesis

The need for better and longer lasting trileaflet valves has led to the design and development of the Abiomed polymeric trileaflet valve prosthesis. In-vitro fluid dynamic studies on sizes 25 and 21 mm valves in the aortic position indicate an overall improvement in performance compared to the Carpentier-Edwards and Ionescu-Shiley tissue valves in current clinical use. The pressure drop studies yielded effective orifice areas of 1.99 and 1.54 cm2, and performance indices of 0.41 and 0.45 for the Nos. 25 and 21 valves, respectively. Leaflet photography studies indicated that the two valve sizes had maximum opening areas of 225 and 145 mm2, respectively, at a normal resting cardiac output. Steady and pulsatile flow velocity measurements with a laser-Doppler anemometer (LDA) system indicate that the flow field downstream of the Abiomed valve is jetlike and turbulent. Maximum mean square axial velocity fluctuations of 55 and 83 cm/s, and turbulent shear stresses of 220 and 450 N/m2 were measured in the immediate vicinity of the nos. 25 and 21 valves, respectively. The Abiomed valves studied had been originally configured for use in valved conduits, and it is therefore our opinion that further improvements can be made to the valve and stent design, which would enhance its fluid dynamic performance.     

In vitro pulsatile flow velocity and shear stress measurements in the vicinity of mechanical mitral heart valve prostheses

A three beam laser Doppler anemometer system was used to study the flow fields created by various types of mitral heart valve prostheses under physiological pulsatile flow conditions. The prosthetic valves studied were: Beall caged disc valve, Bjork-Shiley tilting disc valve, Medtronic-Hall tilting disc valve and St. Jude bileaflet valve. The results indicate that all four prosthetic valve designs studied create very disturbed flow fields with elevated turbulent shear stresses and regions of flow separation and/or stagnation. The observed elevated turbulent shear stresses could cause sublethal and/or lethal damage to red cells and platelets. The regions of flow separation and/or stagnation, could lead to thrombus formation and/or tissue overgrowth on the valve structure, as observed on clinically recovered prosthetic valves.     

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.     

Numerical models for the simulation of flexible artificial heart valves: part I--computational methods

A numerical model of the coupled motion of a flexing surface in a high Reynolds number flow is presented for the simulation of flexible polyurethane heart valves in the aortic position. This is achieved by matching a Lagrangian dynamic leaflet model with a panel method based flow solver. The two models are coupled via the time-dependent pressure field using the unsteady Bernoulli equation. Incorporation of sub-cycling in the dynamic model equations and fast pre conditioning techniques in the panel method solver yields efficient convergence and near real-time simulations of valve motion. The generality of dynamic model allows different material properties and/or geometries to be studied easily and interactively. This interactivity is realized by embedding the models within a design environment created using the software IRIS Explorer. Two flow domains are developed, an infinite domain and an internal domain using conformal mapping theory. In addition bending stress on the valve is computed using a simple stress model based on spline and circle equation techniques.     

Hemolytic anemia following insertion of Ionescu-Shiley mitral valve bioprosthesis.

Hemolytic anemia is a relatively common complication after the replacement of cardiac valves with mechanical prostheses; the prevalence rate varies from 38% to 85%, depending on the prosthesis implanted. However, cardiac valves fabricated from biologic material are associated with a reduced incidence of hemolytic anemia, and to the authors'' knowledge this report is the first to document hemolytic anemia in a patient who had the mitral valve replaced with an Ionescu-Shiley valve. The anemia was not associated with evidence of hemodynamically important mitral regurgitation and was ultimately controlled by iron and folate supplements.     

Earlier accumulation of calcium, phosphorus, and magnesium in the coronary artery in comparison with the ascending aorta, aortic valve, and mitral valve

To explore reasons for a high accumulation of Ca and P occurring in the coronary artery of Thai with aging, the authors investigated age-related changes of elements in the coronary artery, ascending aorta near the heart, and cardiac valves in single individuals, and the relationships in the elements between the coronary artery and either the ascending aorta or cardiac valves. After an ordinary dissection by medical students at Chiang Mai University was finished, the anterior descending arteries of the left coronary artery, ascending aortas, mitral valves, and aortic valves were resected from the subjects. The subjects consisted of 17 men and 9 women, ranging in age from 46 to 76 yr. The element content was analyzed by inductively coupled plasma-atomic emission spectrometry. The average content of Ca and P was the highest in the coronary artery and decreased in the order aortic valve, ascending aorta, and mitral valve. The Ca, P, and Mg content increased in the coronary artery in the fifties and in the ascending aorta, aortic valve, and mitral valve in the sixties. It should be noted that the accumulation of Ca, P, and Mg occurred earlier in the coronary artery than in the ascending aorta, aortic valve, and mitral valve. It was found that with respect to the Ca, P, Mg, and Na contents, the coronary artery correlated well with both the aortic valve and ascending aorta, especially with the aortic valve, but it did not correlate with the mitral valves. This finding suggests that the accumulation of Ca, P, Mg, and Na occurs in the coronary artery together with the aortic valve and ascending aorta, but not together with the mitral valve. Because regarding the accumulation of Ca, P, and Mg, the ascending aorta and aortic valve are preceded by the coronary artery, it is unlikely that the accumulation of Ca, P, and Mg spreads from the ascending aorta or aortic valve to the coronary artery.     

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.     

Numerical Models for the Simulation of Flexible Artificial Heart Valves: Part I - Computational Methods

Abstract

A numerical model of the coupled motion of a flexing surface in a high Reynolds number flow is presented for the simulation of flexible polyurethane heart valves in the aortic position. This is achieved by matching a Lagrangian dynamic leaflet model with a panel method based flow solver. The two models are coupled via the time-dependent pressure field using the unsteady Bernoulli equation.

Incorporation of sub-cycling in the dynamic model equations and fast pre conditioning techniques in the panel method solver yields efficient convergence and near real-time simulations of valve motion. The generality of dynamic model allows different material properties and/or geometries to be studied easily and interactively. This interactivity is realized by embedding the models within a design environment created using the software IRIS Explorer ^TM.

Two flow domains are developed, an infinite domain and an internal domain using conformal mapping theory. In addition bending stress on the valve is computed using a simple stress model based on spline and circle equation techniques.     

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.     

Experiments on dynamic behaviour of a Dacron aortic graft in a mock circulatory loop

Woven Dacron grafts are currently used for the surgical treatment of aortic aneurysm and acute dissection, two otherwise fatal pathologies when aortic wall rupture occurs. While Dacron is chosen for aortic grafts because of characteristics such as biocompatibility and durability, few data are available about the dynamic response of Dacron prosthetic devices and about their side effects on the cardiovascular system. In this study, a Dacron graft was subjected to physiological flow conditions in a specifically-developed mock circulatory loop. Experiments were conducted at different physiological pulsation-per-minute rates. Results show that, in comparison to an aortic segment of the same length, the prosthesis is extremely stiffer circumferentially, thus limiting the dynamical radial expansion responsible for the Windkessel effect in human arteries. The prosthesis is instead excessively compliant in the axial direction and develops preferentially bending oscillations. This very different dynamic behaviour with respect to the human aorta can alter cardiovascular pressure and flow dynamics resulting in long-term implant complications.